motivation

Physical Motivation

Why a general framework for continuous ambient energy conversion in open systems is needed, and what the P(t) master formula sets out to describe.

Last updated:

The Schubart Master Formula, written P(t), is a general mathematical framework for describing the continuous conversion of ambient field energy into usable electrical power in open physical systems. Before introducing the formalism itself, it is worth stating plainly what problem the framework is meant to organise. A great many energy-harvesting concepts — thermoelectric, photovoltaic, piezoelectric, radio-frequency, and more speculative field-coupling schemes — share a common structure: an ambient flux impinges on a converter, interacts with it through some coupling channel, and a fraction of the transferred energy is rectified into a directed, usable current. The purpose of P(t) is to express that shared structure once, in a technology-neutral way, so that specific devices can be treated as special cases rather than as unrelated phenomena.

A unifying description, not a device

It is essential to be clear from the outset about what P(t) is and is not. P(t) is a bookkeeping framework for the power balance of an open converter. It is not a claim that any particular device delivers a particular output, and it does not by itself assert that any specific coupling channel is strong enough to be practically useful. Those are empirical questions that must be settled by measurement for each proposed application. What the framework does provide is a consistent set of variables and a consistent way of combining them, so that competing proposals can be compared on the same terms and so that the assumptions behind any claimed output are made explicit.

In this sense the master formula plays a role analogous to the way a transport integral or a rate equation organises a field: it tells you which quantities matter and how they enter, while leaving the numerical values of those quantities to be supplied by theory or experiment. A reader who wants the compact symbolic statement can consult the equation directly; the meaning of each symbol is catalogued in the notation reference and expanded in variable definitions.

The value of such a unifying object is not merely aesthetic. When a field accumulates many superficially different proposals — each with its own notation, its own implicit assumptions, and its own claimed figure of merit — it becomes very difficult to tell whether two claims are compatible, whether one is a special case of another, or whether a claimed output is even dimensionally possible. A single master formula forces every proposal onto a common axis. To place a device within the framework, one must state its flux distribution, its coupling function, its symmetry-breaking mechanism, its transport channel, and its rectification efficiency; the moment those are stated, the proposal becomes checkable. Proposals that cannot be cast in these terms are not thereby refuted, but the burden shifts to explaining which ingredient the standard chain is missing.

Relation to established transport and rate formalisms

The structure of P(t) echoes a lineage of well-established physical expressions. A Boltzmann transport integral accumulates a current by weighting a distribution function against a velocity and a scattering operator over all momentum states. A radiative-transfer expression accumulates a radiant power by integrating a specific intensity against an absorption coefficient over energy and solid angle. Fermi’s golden rule accumulates a transition rate by weighting a matrix element against a density of final states. In each case a rate or a power is built as an integral of a distributed source against one or more coupling factors. The Schubart Master Formula sits squarely in this tradition: it is a rate-of-energy-delivery integral in which the “source” is the ambient flux and the “coupling” is a product of interaction, directedness, transport, and rectification factors. Recognising this lineage is reassuring, because it means the framework is not an exotic construction but a particular instance of a pattern that has organised transport physics for a century.

The scope of “ambient field energy”

A word is needed about how broadly the phrase ambient field energy should be read. The framework does not privilege any particular field. It applies wherever a distributed flux crosses a boundary and can, in principle, be coupled and rectified: electromagnetic radiation across the spectrum, acoustic and vibrational fields, thermal gradients, and — as one frequently discussed and still speculative example — putative couplings to weakly interacting particle fluxes such as those invoked in neutrinovoltaic proposals. Listing these together is not an endorsement of any of them; it is a statement about the generality of the formalism. The mathematics is indifferent to which field supplies the flux, and that indifference is the point: a single structure covers established and speculative resources alike, with the crucial difference that the parameters entering the formula are measured in the former case and merely assumed in the latter. Throughout this reference the neutrinovoltaic case is used only as an illustration of a proposed application, never as a demonstrated one.

The ambient flux as a resource

The motivating observation is that physical environments are permeated by radiation and field fluxes across an enormous range of energies and directions. Any converter placed in such an environment is continuously bathed in incident flux, denoted Φ(E,θ,t)\Phi(E,\theta,t), which depends on energy EE, incidence angle θ\theta, and time tt. Most of this flux passes through ordinary matter without depositing usable, directed momentum, because the relevant coupling is weak and because randomly directed interactions average to zero net force. The central question the framework poses is therefore not whether flux is present — it always is — but how much of it can, in principle, be coupled, rectified, and delivered as power.

Framing the problem this way immediately isolates the three obstacles any real converter must overcome. First, the coupling strength between the incident field and the converter material, captured by the scattering cross-section σ(E)\sigma(E), is typically small. Second, the interactions must break directional symmetry so that momentum transfer Δp\Delta p is net-directed rather than isotropic. Third, the transferred momentum must be transported and rectified efficiently, which the framework encodes through the phonon velocity vphv_{\text{ph}} and the conversion efficiency ηconv\eta_{\text{conv}}. Each of these is a genuine physical bottleneck, and the master formula makes each one an explicit, separable factor.

Isolating the obstacles as separate factors is more than a notational convenience. It means that a weakness in any single stage can be diagnosed without re-deriving the whole model, and that improvements in different stages can be reasoned about independently to first order. If a proposed converter fails, the framework tells you where it fails: the flux may be adequate but the coupling negligible; the coupling may be strong but the symmetry unbroken; the directed momentum may be generated but lost before it can be rectified. This diagnostic power — the ability to attribute a shortfall to a specific factor — is one of the main practical reasons for writing the balance in factorised, integral form rather than as a single lumped coefficient that hides where the physics actually lives.

Why start from an integral

Because the ambient flux is distributed over energy and angle, the total power cannot be obtained from a single representative value; it must be accumulated over the full spectrum and over all incidence directions. This is why the framework is written as an integral over EE and θ\theta rather than as an algebraic product. The integral form also makes the model honest about spectral dependence: a coupling that is strong at one energy but negligible elsewhere contributes only where σ(E)\sigma(E) and Φ\Phi overlap. The consequences of the open, non-equilibrium setting are developed in open thermodynamic systems and non-equilibrium energy conversion, and the multi-channel character of realistic flux is treated in multi-channel energy flux.

To make the reasoning concrete, consider the differential contribution to the power from a single spectral–angular cell. Writing the integrand as I(E,θ,t)\mathcal{I}(E,\theta,t), the master formula is

P(t)=∬I(E,θ,t)  dE dθ,I=Φ(E,θ,t) σ(E) Δp vph ηconv.P(t) = \iint \mathcal{I}(E,\theta,t)\;dE\,d\theta, \qquad \mathcal{I} = \Phi(E,\theta,t)\,\sigma(E)\,\Delta p\,v_{\text{ph}}\,\eta_{\text{conv}} .

Each factor answers one physical question, and the questions are logically ordered: Is there flux here? (Φ\Phi), does the material couple to it? (σ\sigma), is the coupling directional? (Δp\Delta p), is the momentum transported? (vphv_{\text{ph}}), and is it rectified? (ηconv\eta_{\text{conv}}). A “no” at any stage sends the local integrand to zero, and only cells that answer “yes” to all five contribute to the accumulated total. This is the sense in which the integral is not a mathematical flourish but a faithful ledger of where power can and cannot originate.

An illustrative order-of-magnitude sketch

It helps to see, in rough terms, why the framework does not automatically promise large numbers. The following figures are illustrative placeholders chosen only to show how the factors combine; they are not measured device parameters and imply no claim about any real converter. Suppose a converter presents an effective coupling area of order σ∼10−28 m2\sigma \sim 10^{-28}\,\mathrm{m^2} across a band where the ambient number flux is of order Φ∼1010 m−2 s−1\Phi \sim 10^{10}\,\mathrm{m^{-2}\,s^{-1}} per relevant channel, that each rectified interaction contributes a directed momentum of order Δp∼10−27 kg m s−1\Delta p \sim 10^{-27}\,\mathrm{kg\,m\,s^{-1}} carried at a lattice velocity vph∼103 m s−1v_{\text{ph}} \sim 10^{3}\,\mathrm{m\,s^{-1}}, and that the conversion efficiency is a modest ηconv∼10−2\eta_{\text{conv}} \sim 10^{-2}. Multiplying these order-of-magnitude figures and treating the spectral–angular integral as contributing a bandwidth-and-solid-angle factor of order unity gives an interaction rate per unit area of order Φ σ∼10−18 s−1\Phi\,\sigma \sim 10^{-18}\,\mathrm{s^{-1}} — an extraordinarily small number. The point of the exercise is precisely this: the framework transparently exposes when a proposed channel is many orders of magnitude too weak to matter, rather than hiding that fact inside a single fitted constant. Whether any real channel escapes this smallness is an empirical question the formalism refuses to prejudge; it only makes the arithmetic honest.

The dimensional bookkeeping behind such estimates is set out in dimensional analysis, and readers can vary the factors interactively in the playground.

What the framework deliberately does not claim

Because the master formula is technology-neutral, it is important to state clearly the claims it does not make, so that its generality is not mistaken for advocacy. It does not assert that any particular ambient field carries enough rectifiable energy to power a device; that depends on the numerical value of Φ\Phi, which the environment supplies. It does not assert that any material achieves a large coupling σ\sigma at the relevant energies; that is a materials question. It does not assert that any known mechanism generates a substantial directed momentum Δp\Delta p for a given field; that is precisely the hard, open question for speculative applications. And it does not assert a value for the efficiency ηconv\eta_{\text{conv}}. What the formula does claim is only that if these quantities take certain values, then the delivered power is their integrated product — a conditional, falsifiable statement. This conditional character is a strength: it means the framework can host both established converters, whose parameters are measured, and proposed ones, whose parameters are hypotheses awaiting test, without conflating the epistemic status of the two. A reader evaluating any specific claim should therefore always ask which of the five factors are measured and which are assumed.

What the remaining chapters build

The remainder of this reference builds the formula piece by piece, defines every quantity with its units, delimits the system boundary across which the balance is taken, and describes how the framework can be explored numerically. Throughout, speculative applications — neutrinovoltaics being one frequently discussed example — are labelled as proposed rather than established, so that the mathematical content stands independently of any single physical realisation. The thermodynamic footing appears in open thermodynamic systems; the symmetry-breaking mechanics in non-equilibrium energy conversion; the formal integral in mathematical architecture; and the full derivation in full derivation.

Further reading

  • S. R. de Groot and P. Mazur, Non-Equilibrium Thermodynamics (North-Holland, 1962; Dover reprint 1984). The classic account of flux–force relations in open, driven systems.
  • I. Prigogine, Introduction to the Thermodynamics of Irreversible Processes, 3rd ed. (Wiley-Interscience, 1967). Foundational treatment of entropy production and stationary non-equilibrium states.
  • L. D. Landau and E. M. Lifshitz, Statistical Physics, Part 1, 3rd ed. (Pergamon, 1980). Standard reference for the statistical foundations invoked throughout this framework.
  • N. W. Ashcroft and N. D. Mermin, Solid State Physics (Holt, Rinehart and Winston, 1976). Background for the lattice-transport factors (vphv_{\text{ph}}) that appear in the integrand.
↑