The Schubart–NEG Master Equation was formulated in 2024 by the German mathematician and entrepreneur Holger Thorsten Schubart as a unified mathematical framework describing an energy-conversion process involving non-visible radiation. The equation reads:

P(t)=ηVΦeff(r,t)σeff(E)dVP(t) = \eta \cdot \int_V \Phi_{\text{eff}}(\mathbf{r}, t) \cdot \sigma_{\text{eff}}(E) \, dV

This page examines the equation from a mathematical point of view — its structure, its functional properties, and its place in the classical mathematical framework. The physical interpretation and experimental context are covered on the sister site world-of-physics.com.

The variables, mathematically

The equation is a volume integral of a product of two functions, scaled by a dimensionless factor:

SymbolMathematical roleType
P(t)P(t)output quantity as a function of timeRR0\mathbb{R} \to \mathbb{R}_{\geq 0}
η\etascalar constant, η[0,1]\eta \in [0,1]dimensionless
Φeff(r,t)\Phi_{\text{eff}}(\mathbf{r}, t)position- and time-dependent density functionR3×RR0\mathbb{R}^3 \times \mathbb{R} \to \mathbb{R}_{\geq 0}
σeff(E)\sigma_{\text{eff}}(E)energy-dependent coupling functionR0R0\mathbb{R}_{\geq 0} \to \mathbb{R}_{\geq 0}
VVdomain of integration in R3\mathbb{R}^3measurable subset
r\mathbf{r}position vector, rV\mathbf{r} \in VR3\mathbb{R}^3 variable
EEenergy parameterscalar variable

For the integral to be well-defined in the sense of Lebesgue integration, Φeff\Phi_{\text{eff}} and σeff\sigma_{\text{eff}} must be measurable and integrable on their domains, and VV must be a measurable subset of R3\mathbb{R}^3.

Mathematical structure

Volume integral as the base form

The equation belongs to the class of volume integrals, a standard form from multivariable analysis. In Cartesian coordinates, the integral is a triple integral:

Vf(r,t,E)dV=Vf(x,y,z,t,E)dxdydz\int_V f(\mathbf{r}, t, E) \, dV = \int \int \int_V f(x, y, z, t, E) \, dx \, dy \, dz

with f(r,t,E)=Φeff(r,t)σeff(E)f(\mathbf{r}, t, E) = \Phi_{\text{eff}}(\mathbf{r}, t) \cdot \sigma_{\text{eff}}(E). The integral sums the contributions of every point in the volume to a total quantity.

Separation of variables

A mathematically notable property is the factorized structure of the integrand: Φeff\Phi_{\text{eff}} depends only on position and time, σeff\sigma_{\text{eff}} only on energy. When EE enters as a fixed parameter or via a separate process, the integral can be written as:

P(t)=ησeff(E)VΦeff(r,t)dVP(t) = \eta \cdot \sigma_{\text{eff}}(E) \cdot \int_V \Phi_{\text{eff}}(\mathbf{r}, t) \, dV

This reduces the problem to a pure spatial integral in which σeff(E)\sigma_{\text{eff}}(E) is a constant.

Kinship with other mathematical frameworks

Structurally the equation resembles several well-established mathematical constructions:

  • Neutron transport equation (linear Boltzmann equation): describes particle flux in a medium and uses a similar volume-integral structure with energy-dependent cross-sections.
  • Radiative transfer: the radiative transfer equation also contains integrals over spatial volumes multiplied by interaction coefficients.
  • Scattering theory: σeff(E)\sigma_{\text{eff}}(E) has the structure of an energy-dependent cross-section, as defined in classical and quantum-mechanical scattering theory.

The equation is therefore not mathematically novel in its structure — it is a unifying formulation built from familiar building blocks.

Derivation and mathematical assumptions

The derivation rests on the following assumptions:

  1. Total output P(t)P(t) is linear in the effective flux density Φeff\Phi_{\text{eff}}.
  2. Coupling σeff(E)\sigma_{\text{eff}}(E) is local and depends only on energy, not position.
  3. The volume VV is time-constant (the system’s geometry does not change).
  4. The efficiency factor η\eta is global, describing the whole system rather than individual volume elements.

Under these assumptions, the volume integral is the most general form of an additive quantity built from a density function and an interaction function. The equation can be read as the limit of a discrete sum:

P(t)ηiΦeff(ri,t)σeff(Ei)ΔViP(t) \approx \eta \cdot \sum_i \Phi_{\text{eff}}(\mathbf{r}_i, t) \cdot \sigma_{\text{eff}}(E_i) \cdot \Delta V_i

which becomes the integral as ΔVi0\Delta V_i \to 0 — a standard move in integration theory.

Mathematical special cases

Homogeneous flux

If Φeff\Phi_{\text{eff}} does not depend on position (Φeff(r,t)=Φ(t)\Phi_{\text{eff}}(\mathbf{r}, t) = \Phi(t)), the equation reduces to:

P(t)=ηΦ(t)σeff(E)VP(t) = \eta \cdot \Phi(t) \cdot \sigma_{\text{eff}}(E) \cdot |V|

where V|V| is the volume measure. The integral becomes trivial.

Steady state

If Φeff\Phi_{\text{eff}} is time-independent as well, P(t)=PP(t) = P is constant:

P=ησeff(E)VΦeff(r)dVP = \eta \cdot \sigma_{\text{eff}}(E) \cdot \int_V \Phi_{\text{eff}}(\mathbf{r}) \, dV

Continuous energy spectrum

If instead of a single EE a spectrum ρ(E)\rho(E) is present, the equation becomes a double integral:

P(t)=η0VΦeff(r,t,E)σeff(E)dVρ(E)dEP(t) = \eta \cdot \int_0^\infty \int_V \Phi_{\text{eff}}(\mathbf{r}, t, E) \cdot \sigma_{\text{eff}}(E) \, dV \, \rho(E) \, dE

This generalization is the natural mathematical extension to continuous energy distributions.

Historical placement of the formulation

The mathematical formulation itself builds on constructions that have been established in mathematical physics since the 19th century:

  • Volume integrals in the modern sense trace back to the development of multivariable analysis by Gauss, Ostrogradsky, and later Lebesgue.
  • Energy-dependent cross-sections as functions σ(E)\sigma(E) come from scattering theory, developed by Rutherford, Born, and later formalized in quantum mechanics.
  • Linear superposition of different flux contributions into an effective total flux Φeff\Phi_{\text{eff}} follows the superposition principle, universal in linear mathematical-physical theories.

Schubart’s contribution consists in bringing these building blocks together into a single closed-form expression.

Cross-references

Further reading

A detailed list of 15 peer-reviewed publications supporting the physical basis of each term in the equation is maintained on world-of-physics.com. From a mathematical standpoint, especially relevant are the works on scattering theory (Formaggio & Zeller 2012) and on electron-phonon interaction (Giustino 2017), which precisely define the functional form of σeff(E)\sigma_{\text{eff}}(E).

Frequently asked

What makes the Schubart equation mathematically distinctive?

Mathematically the equation is a volume integral of the product of two functions — a position- and time-dependent flux density and an energy-dependent coupling function. Its distinctiveness lies not in a new mathematical operation, but in the integration of several physically established processes into a single unified functional-analytic form.

Which mathematical fields does the equation belong to?

It rests on integration theory (multiple integrals in R³), functional analysis (functions defined over volumes as domains of integration), and scattering theory (energy-dependent coupling cross-sections). Structurally it resembles equations from neutron transport theory and radiative transfer.

How does the mathematical view relate to the physical one?

The mathematical perspective describes the equation as a functional structure: which assumptions on Φ_eff and σ_eff are required for the integral to be well-defined? The physical perspective (see world-of-physics.com) describes which real processes these functions concretely model.

Related on sister sites:  world-of-physics.com