Stochastic Electrodynamics: Lessons from Regularizing the Harmonic Oscillator

: In this paper, the harmonic oscillator problem in Stochastic Electrodynamics is revisited. Using the exact shape of the Lorentz damping term prevents run-away effects. After introducing a cut-off in the stochastic power spectrum and regularizing the stochastic force, all relevant integrals are dominated by resonance effects only and results are derived that stem from those in the quantum ground state. For an orbit with speciﬁc position and momentum at an initial time, the average energy and the average rate of energy change are evaluated, which stem with each other. Resonance effects are highlighted along the way. An outlook on the hydrogen ground state problem is provided.


Introduction
Quantum mechanics (QM) is a statistical theory, which does no less and no more than provide the Born probabilities for the outcomes of experiments.Though many interpretations have been put forward and various ontologies have been sought within the quantum theory, a deep analysis of the dynamics of a quantum measurement in a rich enough, but solvable, model has strengthened the case that QM provides no more than statistics.It is even capable of being connected to individual measurements, thus addressing the celebrated "measurement problem" [1,2].However, in this philosophy, QM does not provide a framework to describe individual systems in detail.
In practice, in each single experiment, something occurs inside the apparatus which triggers the specific outcome for its pointer.It is our thesis that we do not have a theory for such events.The theory can not be QM, which provides only the statistics of outcomes; it is not classical electrodynamics, nor any other known theory.Nature produces individual events "every day, the whole day", and we are still lacking a theory for these.The challenge to find this theory is as large as it was to find special relativity and quantum mechanics in the early 1900s.
The theory of Stochastic Electrodynamics (SED) is a bold and attractive attempt to describe nature at a deeper level than by quantum theory, aiming to capture individual events.As such, it is an attempt to provide a sub-quantum mechanics.SED is a theory for classical particles embedded in a classical electromagnetic vacuum, with the intensity of its modes being compatible with the quantum zero-point spectrum.In particular, each vacuum mode is a plane wave with energy 1  2 hω, where Planck's constant sets the energy scale of the fluctuation spectrum and appears as a new constant of nature.The SED theory has been well-developed; the reader is referred to the books [3,4].Before addressing deeper questions, one, first, has to demonstrate that SED explains various properties of QM at the statistical level.This has been shown, on general grounds, within a certain flow of arguments [3,4].Hence, specific cases are required to support the derivation and for gaining deeper insights.Progress has been made for harmonic oscillators [5][6][7][8][9][10][11], where a proper ground state emerges.Several phenomena related to oscillators have been explained as well, such as van der Waals forces, the Casimir effect, and the Unruh effect.Additionally, the leading logarithm of the Lamb shift between the hydrogen 1s and 2p states, has emerged in the harmonic oscillator and the hydrogen problems; however, the finite part is not in agreement for the harmonic oscillator problem, and it has not been worked out for the hydrogen problem.
In a recent publication, we considered the hydrogen ground state in SED [12].To lay the basis for a possible reformulation of that problem, we revisit, here, the harmonic oscillator problem, paying special attention to the prevention of peculiarities related to damping and ultraviolet divergences.We outline the problem in Section 2 and consider the average dynamical properties in Section 3. In Section 4, we consider the average progression from a general initial condition.We conclude with a summary.

The Basis of SED
The harmonic oscillator problem for an electron in SED is set by the stochastic differential equation where (in Gaussian units) α ≈ 1/137 is the fine-structure constant.We study this problem on both space and time scales, inspired by the Kepler problem for the hydrogen atom [12].Hence, we consider a typical frequency ω 1 and typical deviation r 1 = √ h/mω 1 such that mω 2 1 r 2 1 = hω 1 .This leads us to define the adimensional ("ad") time, distance, and adimensional charge β, Noting that the dimension of E satisfies Defining ω 0 = ω/ω 1 and dropping the index "ad", we derive the adimensional equation of motion Though β and ω 0 are dimensionless, the combination ωτ c → ω 0 β 2 makes β 2 appear to have the dimension of time.
In the hydrogen problem, with a nucleus of charge Ze one would have ω 2 0 r → r/r 3 .Furthermore, hω 1 = Ze 2 /r 1 = Z 2 α 2 mc 2 implies that ω 1 = Z 2 α 2 mc 2 /h and r 1 = h/Zαmc, and We consider a time window, during which the orbit remains, for many revolutions, in the plane of the temporary orbit; which is possible due to the small value of β.

Statistics of the Stochastic Electric Field
The stochastic electric field E contains polarization indices.It is handy not to be bothered by this technical complication.In our present units, where E ad (t ad ) = (3/2)(c 3 /h) E(t)/ω 2  1 , we will keep the explicit factor √ 3/2 after E(t) is expressed in t ad .Hence, in adimensional units, its spectrum can equally well be presented by the stochastic Gaussian field with the properties [13] Its correlator correctly represents the statistics of the EM field: Notice the factor of ( √ 3/2 ) 2 .This notation and outcome can be verified from the frequency-discretization derived in [12].It is expected that the cutoff τ c can be taken to zero at the end.In the hydrogen problem, the small value of τ c = α 2 Z 2 corresponds, in physical units, to the short Compton time h/mc 2 .
With E = −∂ t A, the correlation function of the vector potential A follows from At the end of Section 4, we will consider the primitive C of A, and its correlation function where ω co ∼ α 3 log 1/α is an adimensional low-frequency cut-off, anticipated below (51), and γ E is Euler's constant.

The Harmonic Oscillator in SED and Its Steady State
The harmonic oscillator problem in SED has been studied by many leaders in the field; see, for example, [5][6][7][8][9][10][11], it is also discussed, at length, in the book [3].This is typically done by taking, in frequency integrals, the contributions from resonances and not bothering much about high-or low-frequency peculiarities.It is our purpose to clarify where regularizations are needed and which form they should have, in order to derive these physically relevant results in a proper fashion.
Let the stochastic field and the particle position have the frequency representation The equation of motion in time and frequency read, respectively, The damping term, if taken from Equation (4), would be D ω ≈ iβ 2 ω 3 r ω .Fortunately, it has been derived, from first principles, in Equation (3.110) of [3].In our notation, this exact result reads where we assumed a similar exponential cut-off exp(−|ω |τ c ) as in the stochastic spectrum.This structure is a truncated convolution, so that it still leads to a product in Fourier space.With (ṙ) ω = −iωr ω , one gets with the function D ω , not to be confused with D ω or the integration measure Dω of Equation ( 6), given by where PV denotes the principal value.Transforming to the time domain, this yields where θ is the Heaviside step function.Indeed, for u < 0, the contour can be closed in the upper half of the complex ω plane, where iωD ω is analytic, as is evident from the middle expression in the first line of ( 14).Hence, D(u) = 0 for u < 0. Likewise, for u > 0, the contour can be closed in the lower half-plane, yielding D(u) = D(u).This explains the causality relation D(t − s) = 0 for s > t.
After replacing D → D in Equation ( 12) for D(t), we may extend the s and u integrals from −∞ to ∞.
For small ωτ c , the last expression in ( 14) yields The δ m term corresponds to a mass renormalization, due to the presence of the electromagnetic field modes.In the units of the hydrogen problem, δm e /m e = δ m = 4α/3π = 0.0031, independent of Z (as it should).The term iβ 2 ω corresponds to the Lorentz damping term D ≈ β 2 ... r in (11).This approximation is known to have run-away solutions r ∼ exp t/β 2 , artifacts that are absent in our exact treatment.For large |ω|τ c , one finds that the mass renormalization drops out and, instead, we have that D ω → −(4β 2 /πτ c )/(ωτ c ) 2 becomes negligible.
The solution of the inhomogeneous equation of motion (11) now reads The poles of G ω determine the complex eigenfrequencies.They follow from Contrary to the approximation (16), which formally allows D ω ∼ −1 at a large imaginary ω = i|ω| with |ω| ∼ (1 + δ m )/β 2 + β 2 ω 2 0 , the exact function D ω is small for all real and complex ω, so there appear no such spurious eigenvalues and run-away solutions ∼ exp(t/β 2 ) in the exact treatment.As expected on physical grounds, there are only solutions near ±ω 0 .For small β and δ m , they read with the approximations giving the leading terms in β 2 .

The Steady State
We now get, from ( 6) and ( 17), the steady value where the first expression can be verified from the frequency-discretization of [13].The integral is finite for τ c → 0 and dominated by the narrow resonance region around ω 0 , with the result There is a similar resonance around ω ≈ ω 0 , yielding ṙ2 (t) = 3 2 ω 0 and 2 ω 0 , which is the ground state energy of the 3d quantum oscillator.However, the large ω limit is only suppressed by the exponential.To evaluate its contribution to the leading order, ω 2 0 can be set to zero, while D ω ∼ −(β 2 /τ c )/(ωτ c ) 2 can be neglected.The remaining integral in (26) is trivial, and brings For ω 0 = 0, this result follows immediately from the free-particle solution (ṙ) ω = βA ω /(1 + D ω ), leading to ṙ(t) ≈ βA(t) and ṙ2 (t) ≈ 3β 2 C AA (0).The large term β 2 C AA (0) ∼ 1/αZ 2 comes from large frequencies ∼ 1/τ c , which are cut off-but not enough-by the factor exp(−|ω|τ c ).
To further suppress the large ω contributions, we propose to subtract the free propagator G 0 ω , The expression r ω = G ω E ω can be written as r ω = G ω E ω , with the renormalized stochastic field Now, the extra factor (ω 0 /ω) 4 |1 + D ω | 2 assures enough suppression of the large ω contributions, so that the resonance at ω 0 is dominant.This yields to the leading order In the new expression for the position fluctuations, the correction φ inherits the logarithmic divergency at ω = 0 from the free particle.This, likely, is addressed by accounting for soft photon emission, leading to a small β 2 log 1/β ∼ α 3 log 1/α "Lamb" correction.The average energy now agrees, to leading order, with that of the ground state of the QM oscillator,

Average Progression of a Specific Orbit
Following specific orbits in time reveals the structure of the dynamics.Due to the stochastic force, this can only be done numerically.The idea to look at the average progression of a collection of orbits, starting at some initial time at a general initial position and speed, was put forward, by de la Peña [14] and Puthoff [15], for the hydrogen problem.Our numerical results [13,16] motivated us to revisit this average progression of a set of orbits [12].To understand better what happens in this approach, we derive it, here, for the harmonic oscillator.
Suppose that the orbit has reached, at time t = 0, a position r 0 and speed v 0 , and has energy The situation is described by adding a proper homogeneous solution to the previous inhomogeneous one.The solution with r(0) = r 0 and ṙ(0) = v 0 reads, for t > 0, where we have already inserted As required, it holds that r ω (0) = v ω (0) = 0. Now, the expectation value of r 2 already involves an integral of the form Dω ω 2 |G ω | 2 , which, as we saw in Equation ( 26), is neatly dominated by the resonances at ±ω 0 .Now consider the "state" as the ensemble of orbits with these initial conditions, progressing under all possible realizations of the stochastic field, with the appropriate weights.The average energy of this state is calculated, as above, by To leading order in β, the classical orbit shrinks due to the radiation, The frequency integrals are dominated by the resonances at ±ω 0 .We need the rule for functions a and b which are smooth around ω 0 .The first terms of r ω and v ω in Equation (30) yield, in the integral (31), the steady state result again, while the absolute-squares of their second terms are of the type a in (33).The cross-terms, finally, have one factor e −γt , but pick up another one because they are of the b-type in (33).Putting this all together, the total energy reads, if we neglect the γ and δ m corrections (except in the exponent), as which shows equilibration towards the average energy (28).The rate of energy change, exhibits that orbits with E 0 > 3 2 ω 0 have the tendency to lose energy, and orbits with E 0 > 3 2 ω 0 to gain energy, again demonstrating the stability of the ground state.The adimensional characteristic timescale is 1/(2γ) = 1/(2β 2 ω 2 0 ), which reads 1/(2ω 2 τ) in the physical units of Equation (1).

Average Change of Energy of a Specific Orbit at Moderate Times
For the hydrogen problem, we have found numerical and analytical support for the thesis that SED does not produce a proper stationary state, but instead leads to self-ionization.For analyzing such cases, the best one can do is to evaluate the average rate of energy change of sets of specific orbits.Let us, therefore, derive this here, in detail, for the harmonic oscillator, with the goal of confirming Equation (35).
The change of energy is, by definition, In the term Ėfield , the integrand diverges before regularization as ω 3 exp(−|ω|τ c ) for large |ω|.
To regularize this, we need to insert two factors ω 2 0 G 0 ω , corresponding to evaluating −β E • ṙ with, again, both its factors E ω replaced by E ω .For ṙ, this was already done-see Equation (30).

The Energy Gain Term
Taking the contributions from ω < 0 and ω > 0 together, we evaluate The propagator G has poles in the lower half plane, G ≈ 1/[(ω For t > 0, we can close the contour in the upper quarter plane, and integrate along the imaginary ω = iz axis Ėfield = 6γ where we skipped the γ ∼ β 2 corrections in v iz and in G iz ≈ −1/(z 2 + ω 2 0 ).Without needing the cut-off term, which has become cos zτ c here, we can express this as The sin ω 0 t term has a logarithmic singularity at y = 0. Averaging over one period, from t to t + P with period P = 2π/ω 0 , this singularity drops out, yielding to leading order in γ The integral decays at 1/ω 0 t, so that decays, in an algebraic and damped oscillatory fashion, to its ensemble average.

The Energy Loss Term
For the energy loss by radiation, the contribution from the classical orbit with damping Averaged over one period, this brings, for β 1 and γ ω 0 , the radiation damping effect which depends on r 0 and v 0 only through E 0 .The fluctuation contribution to this takes the form With D ω = β 2 ( π 2τ c + iω) ∼ γ, this expression is formally of order γ 2 ; however, the resonances connected to contributions from many previous orbits, bring a factor, when summed, of the order ω 0 /γ.The integral is dominated by the imaginary part of D ω .Indeed, this integral yields Ė (fl) rad = 3γω 0 (1 − e −2γt ). (46)

The Energy Balance for a Specific Orbit
We see that the rate of energy gain from the stochastic field quickly goes to its steady value 3γω 0 , while the fluctuation contribution to the energy radiation starts out as zero, and grows as 1 − e −2γt .
Combining (42) with ( 44) and ( 46), the average rate of total energy change of the state starting at t = 0 with energy E 0 is to leading order in the small parameter γ ∼ This is in accordance with (35), derived from the time-dependent energy.

Renormalized Force in the Temporal Description
In order to avoid high-frequency peculiarities, we had to renormalize the noise-see Equation (25).The renormalized equation of motion,