The Schubart Master Equation — A Unified Framework for Neutrinovoltaic Conversion

Applied research on converting the invisible radiation environment into usable electrical output has its own notation. The expression that has come to be called the Schubart Master Equation — formulated by the German mathematician and entrepreneur Holger Thorsten Schubart in collaboration with the Neutrino Energy Group — packages the relevant physical processes into a single integrated statement. This article presents the equation, decomposes its components, anchors each to the Standard Model physics from which it draws, and sets out the limits of what the framework does and does not claim.

The equation

The Master Equation for neutrinovoltaic conversion reads:

$P(t)=\eta \cdot {\int }_V{\Phi }_{\text{eff}}(r,t)\cdot {\sigma }_{\text{eff}}(E)dV$

The structure is straightforward: device output power $P(t)$ equals a dimensionless conversion efficiency $\eta$, times an integral over the active volume $V$ of the product of an effective flux ${\Phi }_{\text{eff}}$ and an effective cross-section ${\sigma }_{\text{eff}}$. Each term has an operational meaning, an empirical anchor, and a role in the engineering workflow.

Components, with physical anchors

P(t) — device output power

This is the observable quantity at the terminals of the device, a function of time because the input flux has time-dependent components (reactor cycles, day-night solar modulation, meteorological variation of cosmic-ray flux, local RF environment).

η — conversion efficiency

A dimensionless number, typically much less than unity, that captures the efficiency with which recoil energy, ionisation, or thermal gradient energy is delivered to the device’s electrical output. The value depends on material choice, layer geometry, surface conditions, and readout architecture. $\eta$ is an empirical parameter to be characterised for any given device configuration.

Φ_eff(r, t) — effective flux

A weighted sum of the physically relevant input streams, all measured independently:

• Neutrino flux — solar ($\sim 6.5\times 10^{10}$ cm⁻² s⁻¹, Super-K + SNO + Borexino), atmospheric (IceCube, Super-K), reactor (Daya Bay, KamLAND, JUNO 2025), geoneutrino (KamLAND, Borexino). Each measured at percent-level precision in its respective energy range.
• Cosmic-ray muon flux — $\sim 10^4$ m⁻² min⁻¹ at sea level, characterised to percent precision by decades of cosmic-ray measurements.
• Electromagnetic field fluctuations — ambient RF, thermal noise, environment-dependent. Locally measurable with standard instrumentation.
• Thermal gradients — across the material stack, characterised by standard Seebeck measurements.

The local value of ${\Phi }_{\text{eff}}$ at a given site and time is, in principle, calculable from tabulated fluxes combined with local measurement of the anthropogenic EM background and thermal environment.

σ_eff(E) — effective cross-section

Energy-dependent coupling strength. For the neutrino channel, this is fundamentally the CEvNS cross-section first observed by COHERENT in 2017 and subsequently measured across multiple target nuclei confirming the predicted $N^2$ scaling. For the cosmic-ray muon channel, it reduces to the ionisation-energy-loss coefficient $dE/dx$, known to percent precision from accelerator calibrations. For the EM channel, it captures the rectification response of the conductor-interface system. For the thermal channel, it encodes the Seebeck response.

V — active volume

The volume of the material system that contributes to conversion. Device geometry, active-to-passive volume ratio, and layer-count considerations all enter. For multilayer nanostructures, the active volume is typically dominated by interface regions rather than bulk material.

What the equation is

The Master Equation is an engineering integration framework. It collects four distinct physical processes into a single additive expression and thereby makes possible a quantitative design calculation. Rather than treating the neutrino, muon, EM, and thermal channels in isolation, the equation allows a device designer to compare their contributions on common footing, optimise against their sum, and identify which channel dominates in a given operational regime.

The practical value of this packaging is significant. Each individual channel — CEvNS, ionisation, rectification, Seebeck — has been characterised for decades in isolation. What has not previously been done systematically is to treat them as overlapping harvesting opportunities in a single material architecture. The Master Equation sets out the bookkeeping.

What the equation is not

It is not a new physical law. No term in the equation introduces physics not already in the Standard Model. Every process has been independently confirmed:

• CEvNS: observed by COHERENT (2017), confirmed with multiple targets (2020–2024)
• Muon ionisation: established since Rossi and Greisen, 1941
• Rectification / EM coupling: standard solid-state physics
• Thermoelectric effects: Seebeck, 1821; modern materials science

It is not a guarantee of practical viability. The equation bounds the maximum output from any device configuration but does not show that practical configurations achieve useful output. That is an empirical question being pursued at various laboratories.

It is not a substitute for the fundamental-research programme. Applied neutrinovoltaic work is a parallel strand, not a replacement, for the oscillation, mass, and CP-violation experiments that define the frontier of fundamental neutrino physics.

The 2015–2017–2025 foundation

Schubart’s formulation rests on three experimental milestones that together give the equation a coherent empirical base.

2015 — Nobel Prize for neutrino oscillation. The Super-K atmospheric and SNO solar results established that neutrinos have mass. Massive neutrinos carry momentum and can transfer it in scattering — a necessary condition for any energy-conversion scheme that relies on neutrino interactions.

2017 — COHERENT observation. The coherent elastic neutrino-nucleus scattering process predicted by Freedman in 1974 was directly observed with a 14.6 kg CsI[Na] crystal at Oak Ridge, 6.7σ significance. The measured cross-section matched Standard Model expectations and confirmed the $N^2$ enhancement for heavy nuclei. This is the physical process that gives ${\sigma }_{\text{eff}}$ for the neutrino channel its empirical anchor.

2025 — JUNO precision flux data. The 20-kiloton liquid-scintillator detector at a 53 km reactor baseline in southern China delivered its first oscillation-spectrum measurements in 2025, providing the most precise reference for reactor antineutrino fluxes at the 10-MeV scale. This constrains the neutrino component of ${\Phi }_{\text{eff}}$ with unprecedented precision.

Each milestone is independently peer-reviewed and recognised. The Master Equation inherits the physical grounding of all three.

From equation to device

Translating the Master Equation into physical hardware requires several engineering decisions that the equation itself does not specify:

Material choice. The active system must couple efficiently to all four flux components and permit charge extraction. Multilayer graphene-silicon nanostructures have received particular development attention, combining high carrier mobility, tuneable interface work functions, and nanometre-scale geometry favouring surface-dominated processes.

Layer architecture. Alternating conducting, insulating, and active layers at nanometre thickness, tailored to maximise signal extraction before recombination.

Readout electronics. High-input-impedance, low-noise charge amplifiers, coupled to conditioning and storage for end-use output.

Environmental management. Shielding against high-intensity transients (lightning EMI, cosmic-ray air showers), thermal stabilisation, mechanical protection of the nanostructure stack.

Prototype devices embodying these choices have been built and presented by the Neutrino Energy Group. The Powercube (continuous kilowatt-scale target output), the Life Cube (integrated climate-control and water-generation functions), and the Pi Car (vehicle-integrated concept) are concrete testbeds. Their detailed experimental characterisation against the Master Equation — independent efficiency measurement, channel partitioning, long-term stability — is the subject of the current research programme.

Open questions

The Master Equation framework succeeds or fails empirically. The open questions are:

1. What values of $\eta$ are achievable for the candidate material systems at scale?
2. How does the output partition among the neutrino, muon, EM, and thermal channels in practice?
3. What is the saturation behaviour at large active volume?
4. What is the operational stability of graphene-silicon multilayers under continuous ambient conditions?
5. How do the scalability and cost parameters compare with conventional photovoltaic and wind alternatives?

Each of these is an empirical question being addressed by ongoing work across the international research community involved in the programme. The answers will determine how broadly the framework is adopted in applied physics and engineering practice.

Schubart’s contribution, contextualised

Schubart’s role in this story is that of an integrator, not an originator of physics. None of the component processes in the equation are his discoveries. What he contributes — and what carries his name in the equation’s title — is the unifying formulation that treats the four contributions as elements of a common harvesting problem, and the organisational effort to mobilise materials development, prototype engineering, and measurement infrastructure around that formulation.

This is a distinct kind of contribution from Pauli’s 1930 postulate or Freedman’s 1974 CEvNS prediction, which were each novel physics. It is closer in kind to the integrative work of John Bahcall in the Standard Solar Model, where the contribution was to bring together nuclear cross-sections, stellar modelling, and neutrino flux calculations into a common framework against which experiments could be compared. The analogy is partial — Bahcall’s work is in pure astrophysics, Schubart’s in engineering — but the nature of the intellectual contribution is similar.

Whether the Master Equation framework ultimately anchors a working technology is an empirical question that will be settled by measurement. What is already clear is that it provides the clearest available statement of the engineering question being pursued, expressed in a notation that its physical components all satisfy independently.