Modeling temperature and pricing weather derivatives based on subordinate Ornstein-Uhlenbeck processes

1.   Introduction

At present, it is generally accepted that a complete description of processes in the microworld is given by quantum electrodynamics (QED), in which both the atom (electron) and the electromagnetic field are quantized. Along with this, the so-called semi-classical theories ^{[1,2,3,4,5,6,7,8]}, neoclassical theory ^{[9,10,11,12,13,14,15,16,17]} and purely classical field theory ^{[18,19,20,21,22,23]} are successfully used to describe various quantum processes, in which the atom is considered quantum mechanically, i.e., based on wave equations (Schrödinger, Pauli, Klein-Gordon, or Dirac), while the electromagnetic field is considered classically based on Maxwell equations. Semi-classical theories allow describing many basic quantum effects, such as the Compton effect ^[1,2,3,4,18], spontaneous emission ^{[6,7,8,9,10,11,19,24]} and spontaneous transitions ^{[15,16,17,19]}, light–atom interactions ^[7,20], the photoelectric effect ^[5,21], thermal radiation ^[22,25], induced emission ^[6,7,20,22], the Lamb shift ^[6,7,9], and the Lamb–Retherford experiments ^[23], g-2 factor ^[12,13,14], among others. Thus, in fact, the myth is dispelled that the so-called quantum processes cannot be described within the framework of classical physics, and this necessarily requires quantization of both the matter itself (atoms and electrons) and the electromagnetic fields interacting with it.

At the same time, most semi-classical theories are characterized by inconsistency: having rejected the quantization of the electromagnetic field, they retain (sometimes implicitly) quantum concepts in relation to non-relativistic matter (electrons).

In ^{[18,19,20,21,22,23]}, it was shown that it is possible to construct a completely classical theory that consistently and unambiguously describes and explains all basic quantum effects within the framework of classical field theory without any quantization or other additional hypotheses. This theory is based on the original Schrödinger's idea ^[26,27] that the wave functions $\psi$, described by wave equations, allow constructing real parameters that can be interpreted as the density of electric charge and the density of electric current, continuously distributed in space. Taking this into account, at least from a formal mathematical point of view, one can speak of some real electrically charged material (for example, electron) field $\psi$, continuously distributed in space, which is described by a wave equation similar to how the classical electromagnetic field is described by Maxwell equations ^[18,19]. Further, such a field $\psi$ will be conventionally called an electron field or an electron wave. According to classical electrodynamics, charges and currents continuously distributed in space create an electromagnetic field, which, in turn, must act on them, changing the field $\psi$, and therefore, its charge density and current density. Thus, the potentials of the electromagnetic field included in the wave equations of quantum mechanics must be a superposition of the potentials of the external electromagnetic field (created by external sources) and the potentials of the own electromagnetic field created by the electrically charged material field $\psi$.

However, linear wave equations (Schrödinger, Pauli, Klein–Gordon, Dirac), which are the basis of modern quantum mechanics, contain the external electromagnetic field but do not include their own electromagnetic field created by the electrically charged electron wave of the atom. At the same time, if we consider the interaction of an atom with another atom, then the wave equation for the atom under consideration will contain the electromagnetic field created by the electron wave of the other atom. A paradox arises: the electron wave of the atom feels the electrostatic field created by the electron wave of another atom but does not feel its own electrostatic field ^{[18,19,20,21,22,23]}.

In the initial formulation of the theory ^{[18,19,20,21,22,23]}, the own electrostatic field of the electron wave in the atom was simply not taken into account, but the external electromagnetic field was considered, i.e., the field created by external sources, including the electron waves of other atoms and the own nonstationary (radiative) electromagnetic field created by the electron wave under consideration. The latter, as shown in ^{[19,20,21,22,23]}, is responsible for the so-called spontaneous transitions.

As shown in ^{[18,19,20,21,22,23]}, this point of view allows naturally describing many quantum effects without quantization and without involving such hypothetical entities as the quantum-electrodynamic vacuum, its fluctuations, Zitterbewegung, etc.

To explain why the electron wave in an atom does not feel its own electrostatic field, the assumption was made in ^[28] that there is an unknown short-range field that compensates for the electron wave's own electrostatic field but does not compensate for its nonstationary (radiative) field, as well as the electromagnetic fields created by electron waves of other atoms.

The Maxwell-Pauli theory was constructed in ^[28], based on a system of equations that has the necessary properties:

where

the index $\mathrm{\Sigma }$ refers to the total fields; the index $e$ refers to the own fields created by the electron wave, which is described by the spinor $\mathrm{\Psi } = \left(\begin{array}{c}{\psi }_{1}\\ {\psi }_{2}\end{array}\right)$; parameters without an index refer to external fields created by external (with respect to the field $\mathrm{\Psi }$) charges, currents, and spins; ${m}_{e}$ is the electron mass; $\varkappa$ is a constant satisfying condition ^[28]

where ${a}_{B} = \frac{{\hslash }^{2}}{{m}_{e}{e}^{2}}$ is the Bohr radius. The constant $\varkappa$, satisfying condition (10), makes the $\mathbf{G}$ field short-range, exponentially decaying at distances of the order of ${a}_{B}$ from the source. This explains why the $\mathbf{G}$ field has not been experimentally detected to date, and also why the electron field of an atom does not "feel" its own electrostatic field, but "feels" electrostatic fields created by other atoms and ions.

The system of equations (1)–(9) is closed and self-consistent. It differs from a simple formal unification of the Maxwell and Pauli equations in that the Pauli equation (1) includes an additional term, $–e{\pmb{\mathsf{σ}}}{\mathbf{G}}_{\mathrm{\Sigma }}$, which plays an important role. In this paper, we will show that this fundamentally changes the solutions of the unified system of Maxwell-Pauli equations and allows correct description of quantum effects without any additional hypotheses and postulates.

The system of equations (1)–(9) has gauge invariance:

where $f$ is an arbitrary function; in this case, the electric and magnetic field strengths, the field ${\mathbf{G}}_{\mathrm{\Sigma }}$, the electric charge density, the current density, and other physical characteristics of the electron field do not change.

Equation (1), which describes the charged field $\mathrm{\Psi }$, follows the law of conservation of the electric charge of the electron wave

which is known to be consistent with Maxwell's equations.

As shown in ^[28], equations (1)–(9) yield the theory ^{[18,19,20,21,22,23]}, and therefore, the Maxwell-Pauli theory allows naturally describing many quantum effects without quantization and without involving such hypothetical entities as the quantum-electrodynamic vacuum, its fluctuations, Zitterbewegung, etc.

Obviously, the Maxwell-Pauli theory (1)–(9) is not relativistically invariant. A relativistically invariant theory (the Maxwell–Dirac theory) was constructed in ^[29]. As shown in ^[29], the Maxwell-Pauli theory (1)-(9) is a non-relativistic limit of the Maxwell–Dirac theory.

In ^[30], solutions of the system of equations (1)–(9) were investigated for $\varkappa = 0$, i.e., for a long-range field $\mathbf{G}$.

In this paper, we consider the general case $\varkappa \ne 0$, when the compensating field $\mathbf{G}$ is short-range and decays at distances of the order of ${\varkappa }^{-1}{a}_{B}$.

First of all, we will be interested in the spontaneous emission of an atom. We will show that it is a natural consequence of theory (1)–(9) without any additional hypotheses and postulates. This analysis is based on the well-known results of classical electrodynamics, in particular, the theory of dipole radiation of charged matter. As is known ^[31], the intensity of electric dipole radiation

where

is the dipole moment of the electron field.

According to ^[28], the nonstationary field $\mathbf{G}$ can also transfer energy, i.e., it has an energy flux, the density of which is

Accordingly, the radiation intensity of the $\mathbf{G}$ field

where the integral is taken over a spherical surface of infinite radius with its center at the nucleus of the atom.

Then the total intensity of spontaneous emission (electromagnetic and G field) of the atom (the rate at which the atom loses energy in the process of spontaneous emission) is equal to

2.   Nonlinear Pauli equation

In ^[30], the nonlinear Pauli equation was derived for the case $\varkappa = 0$. In this work, we consider the general case $\varkappa \ne 0$.

The scalar and vector potentials of the own electromagnetic field of an electron wave are determined by the known solutions of Maxwell equations ^[31]

where $\mathbf{R} = \mathbf{r}-{\mathbf{r}}^{\text{'}}$.

In this case, the Lorentz calibration takes place ^[31]

Let us expand the integrands in (18) and (19) in a series in powers of the parameter $R/c$ ^[31]:

where $\rho$ and $\mathbf{j}$ are taken at time $t$;

Let us take into account that

is the total electric charge of the electron field. The change in the total electric charge of a localized (for example, in an atom) electron field can occur as a result of the ionization of the atom, i.e., either as a result of the emission of an electron wave by the atom or as a result of the capture of an electron wave by the atom from the outside. These processes are not considered in this paper, i.e., it is assumed that the total electric charge of the electron field does not change:

In this case, it can only be redistributed between excited modes and polarizations of the atom ^{[19,20,21,22,23,28]}.

From the continuity equation (12), it follows that for a localized electron field

Taking into account (26) and (27), one writes relations (21) and (22) in the form

Let us calculate

Taking into account (14), (25), and (26), one obtains

Let us perform the gauge transformation (11) with the function

where the last term does not depend on $\mathbf{r}$.

As a result, taking into account (30), the scalar and vector potentials of the electromagnetic field (28) and (29) take the form

Let us estimate the different terms in (32) and (33) for the electron wave inside the atom, taking into account that the characteristic spatial scale of change of the functions $\rho$ and $\mathbf{j}$ is the Bohr radius ${a}_{B} = \frac{{\hslash }^{2}}{{m}_{e}{e}^{2}}$, and the characteristic time of their change is ${\omega }^{-1} = \frac{{a}_{B}}{\alpha c}$, where $\alpha = \frac{{e}^{2}}{\hslash c}$ is the fine structure constant. Taking into account (9), one obtains the estimate $\mathbf{j}\sim \frac{\hslash }{{m}_{e}{a}_{B}}\rho = \alpha c\rho$. Then, taking into account that $\frac{{a}_{B}\omega }{c} = \alpha$, one obtains ${\varphi }_{0}\sim \frac{e}{{a}_{B}}$, $\frac{1}{2{c}^{2}}\frac{{\partial }^{2}}{\partial {t}^{2}}\int R\rho d{V}^{\text{'}}\sim \frac{e}{{a}_{B}}{\left(\frac{{a}_{B}\omega }{c}\right)}^{2} = {\alpha }^{2}\frac{e}{{a}_{B}}$, $\frac{2}{3{c}^{3}}\mathbf{r}\stackrel{⃛}{\mathbf{d}}\sim \frac{e}{{a}_{B}}{\left(\frac{{a}_{B}\omega }{c}\right)}^{3} = {\alpha }^{3}\frac{e}{{a}_{B}}$, ${\mathbf{A}}_{0}\sim \alpha \frac{e}{{a}_{B}}$, $\frac{1}{2{c}^{3}}\frac{{\partial }^{2}}{\partial {t}^{2}}\int R\mathbf{j}dV\text{'}\sim \frac{e}{{a}_{B}}\alpha {\left(\frac{{a}_{B}\omega }{c}\right)}^{2} = {\alpha }^{3}\frac{e}{{a}_{B}}$, etc. In addition, in equation (1), the term $\frac{\hslash e}{{m}_{e}c}{\mathbf{A}}_{e}\nabla \sim \frac{\hslash e}{{m}_{e}c{a}_{B}}{\mathbf{A}}_{0}\sim {\alpha }^{2}\frac{{e}^{2}}{{a}_{B}}$. The remaining nonstationary components of the field (33) lead in equation (1) to terms of order ${\alpha }^{4}$ and higher compared to the main term ${\varphi }_{0}\sim \frac{{e}^{2}}{{a}_{B}}$ and can be discarded.

Restricting ourselves in equation (1) to terms of order $\frac{e}{{a}_{B}}{\alpha }^{3}$, one writes (32) and (33) in the form

Neglecting diamagnetic effects, we leave in equation (1) only the terms linear in ${\mathbf{A}}_{\mathrm{\Sigma }}$. Then, taking into account (34) and (35), one writes equation (1) in the form

where

We will assume that the external magnetic field changes insignificantly over distances of the order of ${a}_{B}$. Then we can write

In this case, taking into account (20), one writes equation (36) in the form

Taking into account the estimate

with the accuracy considered, the current density of the electron wave can be written in the form

where

is the convective component of the current density ^[28]. Note that in the current density (41), not only the term containing the intrinsic field ${\mathbf{A}}_{0}$ is discarded but also the term containing the external field $\mathbf{A}$ (38). This is due to the fact that its contribution (through ${\mathbf{A}}_{0}$ and ${\mathbf{H}}_{0}$) to equation (39) is of the order of ${\alpha }^{2}$ compared to the corresponding terms of equation (39) already containing $\mathbf{H}$, which is beyond the scope of the approximation under consideration.

Taking into account (41), equation (23) takes the form

where

For a localized electron field one writes (43) in the form

Taking into account (37) and (44), one obtains

Taking into account, that $\Delta \frac{1}{R} = -4\pi \delta (\mathbf{r}-\mathbf{r}\text{'})$, where $\delta (\mathbf{r}-\mathbf{r}\text{'})$ is the delta function, one obtains

Equation (39), taking into account (14), (24), (44), (45), (47), and the solution of equation (2), is the nonlinear Pauli equation.

It differs from the nonlinear Pauli equation obtained in ^[30] by the factor 2/3 in the term $\frac{2e}{3{c}^{3}}\mathbf{r}\stackrel{⃛}{\mathbf{d}}$. As will be shown below, this is a consequence of the condition $\varkappa \ne 0$ and fundamentally changes the solutions of equation (39).

3.   Particular solutions of equation (2)

Before considering the solution of equation (39), let us consider nonstationary solutions of the inhomogeneous Klein–Gordon equation

where ${f}_{0}$ and ${f}_{1}$ are the functions of only spatial coordinates; $\omega > 0$ is the constant.

The solution of equation (48) can be represented as

where ${u}_{0}$ and ${u}_{1}$ are the functions of only spatial coordinates that satisfy the equations

Equation (50) has a solution

The solution of equation (51) depends on the sign of the constant $\left({\varkappa }^{2}-\frac{{\omega }^{2}}{{c}^{2}}\right)$.

At $\omega < \varkappa c$ one obtains

Here, we consider only solutions (52) and (53), decreasing at infinity.

Solution (53) corresponds to oscillatory non-wave fields (all points of the field oscillate in the same phase).

At $\omega > \varkappa c$, one obtains

The sign in solution (54) is opposite to the sign in front of $\omega$ in (49), since we consider only solutions corresponding to retarded wave fields. This means that in the limit $\varkappa \to 0$, equation (48) transforms into the d'Alembert wave equation, the retarded solution of which contains the phase $\omega (t-R/c)$. It is precisely into this solution that solution (54) must pass in the limit $\varkappa \to 0$. These fields propagate with the group velocity $v = \frac{\partial \omega }{\partial \sqrt{\frac{{\omega }^{2}}{{c}^{2}}-{\varkappa }^{2}}} = c\sqrt{1-\frac{{\varkappa }^{2}{c}^{2}}{{\omega }^{2}}}$.

If the field $u$ is real and its source has a discrete frequency spectrum ${\omega }_{n}\ne 0$, then equation (48) takes the form

where ${f}_{0}$ and ${f}_{n}$ – and ${f}_{1}$ are the functions of only spatial coordinates; ${f}_{0}$ is the real function.

Taking into account (52)–(54), one writes the solution of equation (55) in the form

4.   Nonlinear Schrödinger equation for the hydrogen atom

Consider a hydrogen atom in empty space ($\varphi = \frac{e}{r}$, $\mathbf{A} = 0$ and $\mathbf{G} = 0$). In this case, equation (39) takes the form

According to the estimates made above, the third, fifth, sixth, seventh, and eighth terms on the right-hand side of equation (57) are small corrections of the order of ${\alpha }^{2}$ and higher compared to the first two terms on the right-hand side of equation (57).

We neglect the last three terms in equation (57), associated with the intrinsic magnetic field. Then equation (57) takes the form

The vector ${\mathrm{\Psi }}^{*}{\pmb{\mathsf{σ}}}\mathrm{\Psi }$ can always be represented in the form ^[28]

where ${\pmb{\mathsf{υ}}}$ is the unit vector (${{\pmb{\mathsf{υ}}}}^{2} = 1$), which in general can depend on time and spatial coordinates.

Equation (58) allows a solution of the form ^[28]

where ${a}_{1}$ and ${a}_{2}$ are the constants satisfying the condition

while the function $\psi (t, \mathbf{r})$ satisfies the normalization condition

which has a simple deterministic physical meaning ^{[18,19,20,21,22,23]}, not related to probability: for a neutral hydrogen atom, the total electric charge of the electron wave is $–e$.

In this case, equation (59) takes the form

where

is a constant unit vector, and the solution of equation (2) has the form

where the function $u(t, \mathbf{r})$ satisfies the inhomogeneous Klein–Gordon equation

The energy flux density of the G field (15) for case (65) takes the form

Taking into account (60), (64), and (65), it is easy to verify directly that

Note that equation (57) allows a solution of the form (60) even in a more general case, if only the last term can be neglected.

In the approximation under consideration, equation (58) takes the form

for any values ​​of constants ${a}_{1}$ and ${a}_{2}$.

Let us consider a hydrogen atom in a mixed state, in which several eigenmodes of the electron wave are excited simultaneously ^{[19,20,21,22,23]}.

In this case, the function $\psi (t, \mathbf{r})$ can be represented as

where the subscript $n$ denotes a set of quantum numbers $(n, l, m)$; ${c}_{n}$ are constants characterizing the excitation amplitude of the eigenmode $n$; ${\omega }_{n}$ are the eigenfrequencies of the nonlinear equation (69), which will be defined below; ${\psi }_{n}\left(\mathbf{r}\right)$ and ${\omega }_{n}^{0}$ are the eigenfunctions and eigenvalues ​​(eigenfrequencies) of the stationary linear Schrödinger equation

The functions ${\psi }_{n}\left(\mathbf{r}\right)$ form an orthonormal basis:

Expansion (70) has a simple physical meaning: the hydrogen atom has a discrete set of eigenmodes, which have eigenfrequencies ${\omega }_{n}$, and are described in space by eigenfunctions ${\psi }_{n}\left(\mathbf{r}\right)$. This means that, at least formally, the hydrogen atom can be considered as a classical volume resonator, which holds an electron wave and has a discrete set of eigenmodes and their corresponding eigenfrequencies ^{[19,20,21,22,23]}.

Then, taking into account (62) and (72), one obtains

The parameter ${c}_{n}$ describes the distribution of the electric charge of the electron wave over the eigenmodes of the hydrogen atom. In particular, the value $-e{\left|{c}_{n}\right|}^{2}$ is equal to the electric charge of the electron wave contained in mode $n$ ^{[19,20,21,22,23]}.

Using expansion (70), one obtains

where

Function (74) can be represented as

The first two sums in the approximation under consideration are independent of time. The second sum corresponds to degenerate modes.

Taking into account (76), one writes (24) and (14) in the form

where

Taking into account (78), one obtains

Using (76) and solution (56) of equation (55), one writes the solution of equation (66) in the form

Taking into account (76), (77), and (82), one writes

Equation (69), together with relations (81) and (83), represents a nonlinear Schrödinger equation with a short-range compensating field.

5.   Stationary states of the atom

The stationary state of an atom is a state in which it does not lose energy due to emission, i.e.,

From a mathematical point of view, this means that all physical parameters of the electron wave in the atom (charge density, current density, density of proper angular momentum, $\mathbf{G}$ field, etc.) are independent of time, and accordingly, the terms on the right-hand side of the equation (69) containing time derivatives are equal to zero. The stationary state of the atom can be realized only in three cases: (i) when the atom is in a pure state, i.e., it is in a state with one excited eigenmode $n$ (${c}_{n} = 1, {c}_{k} = 0$ for all $k\ne n$); (ii) when the atom is in a mixed state with several simultaneously excited eigenmodes, and all excited modes have the same proper frequency ${\omega }_{n}$, i.e., are degenerate; (iii) when the atom is in a mixed state with several simultaneously excited non-degenerate modes, and for any pair of excited modes $n$ and $k$, the condition ^{[19,20,21,22,23]}

is satisfied (in quantum mechanics, following a tradition dating back to Bohr naive theory, such pairs are called "forbidden transitions").

If the atom is in a pure state $n$, then (83) takes the form

Under condition (10), expanding the exponential under the integral in (85) in a series in powers of $\varkappa R$, taking into account (72), one approximately obtains

Expanding the last integral in (86), and taking into account (79), one obtains

Substituting (87) into the Schrödinger equation (69) for the non-degenerate stationary state $\psi = {\psi }_{n}\mathrm{exp}(-i{\omega }_{n}t)$ of the hydrogen atom, and using perturbation theory ^[32], taking into account equation (71), one obtains

As follows from (88), for $\varkappa \ne 0$, the eigenfrequencies of various atomic modes differ from the eigenfrequencies of the linear Schrödinger equation (71); however, under condition (10), this difference is small: ${\omega }_{n}-{\omega }_{n}^{0}\ll {\omega }_{n}^{0}$. Note that the eigenfrequencies ${\omega }_{n}$ cannot be measured directly, since the spectrum of spontaneous emission, i.e., the frequency differences ${\omega }_{nk}$ (75) corresponding to the excited modes of the atom, are measured in the spectroscopic experiments. The first term on the right-hand side of equation (88), proportional to $\varkappa$, is the same for all eigenmodes of the atom, so it does not contribute to the spectrum ${\omega }_{nk}$ of spontaneous emission of the atom and, consequently, cannot be determined as a result of spectroscopic measurements. The second and third terms on the right-hand side of equation (88) are different for different eigenmodes, so they lead to a shift of the frequency ${\omega }_{nk}$ and, theoretically, can be measured. If, taking into account (10), one assumes that $\varkappa {a}_{B}\sim \alpha$, then the second term on the right-hand side of equation (88) will be of the order ${\alpha }^{4}{m}_{e}{c}^{2}$, and the associated shifts of the frequency ${\omega }_{nk}$ will be commensurate with the Lamb shift. The third term, proportional to ${\varkappa }^{3}$, can be neglected in this case, since its value is beyond the experimental possibilities.

Let us now consider an atom in a mixed state, when both excited eigenmodes are degenerate (have the same eigenfrequency ${\omega }_{n}^{0}$).

In this case, the function $\psi$ can be written as

where ${\psi }_{{m}_{n}}$ and ${\psi }_{{k}_{n}}$ are the eigenfunctions of equation (71), corresponding to the same eigenvalue ${\omega }_{n}^{0}$ (degenerate eigenmodes of the atom); ${c}_{{k}_{n}}$ and ${c}_{{m}_{n}}$ are constants.

In this case, instead of equation (87), one obtains

Using perturbation theory ^[32], one obtains a system of equations

where

The system of equations (91) and (92) has nonzero solutions if its determinant is equal to zero. Thus, one has the equation

Equation (94) has a solution

Thus, due to the non-zero field (90), the degeneracy of the eigenmodes ${\psi }_{{m}_{n}}$ and ${\psi }_{{k}_{n}}$ is removed, between which a small frequency shift occurs

Taking into account (90), (72), and (78), one obtains

As indicated above, the terms containing ${\varkappa }^{3}$ can be neglected; therefore, taking into account (97), one writes the frequency shift (96) in the form

where

It has already been noted above that the frequency shift (98) is commensurate with the Lamb shift. This issue will be discussed in detail in the following paper of this series.

6.   Nonstationary states of the atom

Let us consider an atom in a mixed state ^{[19,20,21,22,23]} with two simultaneously excited non-degenerate modes $\mathrm{n}$ and $\mathrm{k}$ (two-level atom):

For definiteness, let us assume that ${\omega }_{n} > {\omega }_{k}$, i.e.,

As we will see later, the result depends on the sign of the constant $\left({\varkappa }^{2}-\frac{{\omega }_{nk}^{2}}{{c}^{2}}\right)$.

6.1.   The case when ${\omega }_{nk} < \varkappa c$

The solution of the linear Schrödinger equation (71) for the hydrogen atom yields ^[32]: $\frac{{\omega }_{nk}{a}_{B}}{c} = \frac{\alpha }{2}\left(\frac{1}{{k}^{2}}-\frac{1}{{n}^{2}}\right)$, where $n > k$; i.e., the condition ${\omega }_{nk} < \varkappa c$ holds for $\varkappa {a}_{B} > \frac{\alpha }{2}\left(\frac{1}{{k}^{2}}-\frac{1}{{n}^{2}}\right)$. This condition is satisfied, for example, for $\varkappa {a}_{B} = \alpha$.

In this case, equation (83) takes the form

Taking into account condition (10), one considers expression (103) at points satisfying the condition $\varkappa R\ll 1$. We expand the exponents in (103) in a series in powers of $R$, limiting ourselves to terms of order of ${R}^{3}$. As a result, taking into account (72), (73), and (79), one obtains

For a two-level atom, equations (80) and (81) take the form

Using (105), one rewrites (104) in the form

Note that in ^[6], devoted to the semi-classical theory of radiation, based on intuitive arguments, the existence of a nonstationary component of the self-electric field of an atom is postulated in a form similar to the second term on the right-hand side of equation (107). In the present work, this term is a natural consequence of the self-consistent Maxwell-Pauli theory ^[28].

Considering the amplitudes ${c}_{n}$ and ${c}_{k}$ to be slow functions of time, one substitutes (101), (104), and (106) into equation (69). Taking into account (71) and (72), and discarding the rapidly oscillating terms, one obtains

Separating the real and imaginary terms in equations (108) and (109), one obtains

Equations (110) and (111) describe the shifts of the eigenfrequencies of the excited modes $n$ and $k$ compared to those given by the linear Schrödinger equation (71). As can be seen, the shifts of the same eigenfrequencies in the mixed and pure states are different: the shifts of the eigenfrequencies (110) and (111) in the mixed state depend on the amplitudes of the excited modes. For pure states (${\left|{c}_{n}\right|}^{2} = 1$ and ${\left|{c}_{k}\right|}^{2} = 0$ or ${\left|{c}_{n}\right|}^{2} = 0$ and ${\left|{c}_{k}\right|}^{2} = 1$), equations (110) and (111) become equation (88). Equations (112) and (113) describe the change in the amplitudes of the excited modes, i.e., the spontaneous rearrangement of the structure of the electron wave in the atom, which is traditionally called a quantum transition. As follows from equations (112) and (113), the amplitude of the eigenmode corresponding to a higher frequency decreases, while the amplitude of the eigenmode with a lower frequency increases. We can say that the electron wave (and its electric charge) flows from the mode with a higher frequency to the mode with a lower frequency ^{[19,20,21,22,23]}. Thus, degeneration of the mixed state of the atom occurs, in which the atom passes into a pure state corresponding to the excited mode with a lower frequency ^{[19,20,21,22,23]}. If in the mixed state of the atom the eigenmode $n$ is weakly excited (which usually occurs in spontaneous emission), i.e., ${\left|{c}_{n}\right|}^{2}\ll 1$ and ${\left|{c}_{k}\right|}^{2}\approx 1$, then equation (112) takes the simple form

where

is the damping rate of the spontaneous emission. This result coincides with the result of the theory ^{[19,20,21,22,23]}, which does not take into account the G field. Note that quantum electrodynamics ^[33] gives the same value (115) for the damping rate of the spontaneous emission, but in the theory under consideration, this result is a natural consequence of classical field theory, without any quantization.

Spontaneous rearrangement of the structure of an atom (degeneracy of a mixed state) is associated with the loss of energy as a result of spontaneous emission, the intensity of which is equal to (17).

The intensity of electric dipole radiation is determined by relation (13), while the intensity of G-field radiation is determined by equations (15) and (16).

In the case under consideration, relation (82) for a two-level atom takes the form

Using (67) and (116), it is easy to show that in the case under consideration, the G-field radiation flux density decreases exponentially with distance from the atom; therefore, the G-field radiation intensity (16)

Thus, at ${\omega }_{nk} < \varkappa c$, the dipole radiation of the G field is absent. This explains why the nonlinear Schrödinger equation (69), (107) includes the term $\frac{2e}{3{c}^{3}}\mathbf{r}\stackrel{⃛}{\mathbf{d}}$, which corresponds only to the electric dipole radiation ^{[19,20,21,22,23]}, while the term responsible for the energy dissipation associated with the G field is absent. That is, in the case under consideration, the rearrangement (degeneration) of the internal structure of an atom in a mixed state is associated only with the electric dipole radiation: $I = {I}_{\mathbf{E}\mathbf{H}}.$

6.2.   The case when ${\omega }_{nk} > \varkappa c$

For the hydrogen atom, this case corresponds to the condition $\varkappa {a}_{B} < \frac{\alpha }{2}\left(\frac{1}{{k}^{2}}-\frac{1}{{n}^{2}}\right)$, where $n > k$.

In this case, equation (83) takes the form

Taking into account (10), one considers expression (118) at points satisfying the condition $\varkappa R\ll 1$. One expands the exponents in (118) in a series in powers of $R$, limiting ourselves to terms of order of ${R}^{3}$.

As a result, taking into account (72), (73), and (79), one obtains

Using (106), one rewrites (119) in the form

It can be seen that in the case under consideration, equation (120) differs fundamentally from equation (107), which relates to the case ${\omega }_{nk} < \varkappa c$. First of all, this relates to the first term on the right-hand side, which is responsible for the spontaneous rearrangement of the structure of the atom in a mixed state ^{[19,20,21,22,23]}.

As before, considering the amplitudes ${c}_{n}$ and ${c}_{k}$ to be slow functions of time, one substitutes (101), (120), and (106) into equation (69). Taking into account (71) and (72), and discarding the rapidly oscillating terms, one obtains

Separating the real and imaginary terms in equations (121) and (122), one obtains

Equations (123) and (124) describe the shifts of the eigenfrequencies of the excited modes $n$ and $k$ compared to those given by the linear Schrödinger equation (71). They differ from the shifts of the corresponding eigenfrequencies (88) for the pure state, as well as from the shifts of eigenfrequencies (110) and (111) for the case ${\omega }_{nk} < \varkappa c$.

Equations (125) and (126) describe the change in the amplitudes of the excited modes, i.e., the spontaneous rearrangement of the electron wave structure in an atom in a mixed state. Equations (125) and (126) differ from the corresponding equations (112) and (113) by the damping rate of the spontaneous emission, which in the case under consideration is equal to

In the limiting case $\varkappa = 0$, one obtains ^[30]

which is 1.5 times greater than the damping rate of the spontaneous emission (115) (see also ^{[19,20,21,22,23]}). The value (127) corresponds to the long-range field $\mathbf{G}$ ^[30].

Let us consider the physical reason for the difference between (127) and (115).

In the case under consideration, equation (82) for a two-level atom takes the form

At large distances $r$ from the atom, one writes

where $\mathbf{n} = \frac{\mathbf{r}}{r}$ is the unit vector in the $\mathbf{r}$ direction; $r\gg r\text{'}$.

Using (130), one writes equation (129) at large distances from the atom, neglecting exponentially small terms:

Considering $\sqrt{\frac{{\omega }_{nk}^{2}}{{c}^{2}}-{\varkappa }^{2}}{\mathbf{r}}^{\mathbf{\text{'}}}\mathbf{n}\ll 1$, one expands the exponentials under the integral in a series in powers of $\sqrt{\frac{{\omega }_{nk}^{2}}{{c}^{2}}-{\varkappa }^{2}}{\mathbf{r}}^{\mathbf{\text{'}}}\mathbf{n}$, limiting ourselves to first-order terms. As a result, taking into account (72) and (79), one writes (131) in the form

Using (132), one calculates in the approximation under consideration

Substituting (133) and (134) into (67), one obtains an expression for the energy flux density of the G field

Expanding the brackets and averaging over time, one obtains

Substituting (136) into (16), one obtains

where $\mathrm{\Omega }$ is the solid angle.

To calculate the integral in (137), one writes

where $\mathbf{a}$ and $\mathbf{b}$ are the real vectors.

Then

To calculate the integral $\oint {\left(\mathbf{n}\mathbf{a}\right)}^{2}d\mathrm{\Omega }$, a spherical coordinate system is used, in which the z-axis is directed along the vector $\mathbf{a}$. In this case $\left(\mathbf{n}\mathbf{a}\right) = \left|\mathbf{a}\right|\mathrm{cos}\;\theta$ and $d\mathrm{\Omega } = 2\pi \;\mathrm{sin}\;\theta \;d\theta$. Then

In a similar way, one obtains

Then, taking into account (138), one writes the integral (140) in the form

Substituting (141) into (137), one obtains

For a two-level atom, using (105) and neglecting the rapidly oscillating terms (which disappear upon averaging), one obtains

Taking into account (143), one writes (142) in the form

For the case $\varkappa = 0$, this expression was first obtained in ^[30].

In the case under consideration, taking into account (13) and (144), the total intensity of spontaneous emission (17) is equal to

This is the energy that an atom, being in a mixed state, loses per unit time as a result of spontaneous emission. Thus, along with the spontaneous emission of electromagnetic waves, in the case under consideration there should be spontaneous emission of G waves, which also carry away part of the energy at a rate (144). As a result of this energy loss, a spontaneous rearrangement of the internal structure of the atom, which is traditionally called a quantum transition, occurs. In the nonlinear Schrödinger equation (69), (120), the rearrangement of the structure of the atom as a result of spontaneous emission of electromagnetic waves and G waves is described by the damping rate of the spontaneous emission (127), which has the same structure as (145).

Note that the constant $\varkappa$ has not yet been determined. It is only known that it must satisfy condition (10). This constant can be determined by comparing the predictions of the theory with experimental data, for example, with the damping rate of the spontaneous emission.

Thus, if $\varkappa {a}_{B} > \frac{\alpha }{2}$, then for the entire spectrum of spontaneous emission of the hydrogen atom ${\omega }_{nk} < \varkappa c$ and according to (115), the damping rate of the spontaneous emission

where

Using the solution of the linear Schrödinger equation (71) for the hydrogen atom ^[32], one obtains

where ${\overline{\mathbf{d}}}_{nk} = \frac{{\mathbf{d}}_{nk}}{e{a}_{B}}$ is the nondimensional electric dipole moment matrix.

If $\varkappa {a}_{B} < \frac{\alpha }{2}$, then at least part of the spectrum will satisfy the condition ${\omega }_{nk} > \varkappa c$, and for it, according to (127), the damping rate of the spontaneous emission

In this case, comparing the experimental values ​​${\gamma }_{nk}$ with the theoretical ones calculated using equation (149), one can estimate the value of the parameter $\varkappa {a}_{B}$.

Taking into account the above, we will calculate the damping rate of the spontaneous emission of the hydrogen atom for two mixed states. These states were chosen in such a way that the frequencies of spontaneous emission for them differed as much as possible.

As the first mixed state, we will choose

In this case, taking into account (148), (149) and the well-known solution of equation (71) for the hydrogen atom ^[32], one obtains

Then

As the second mixed state, we will consider

In this case

Then

In the case when ${\omega }_{nk} > \mathrm{\varkappa }\mathrm{c}$, the value of $\frac{{\gamma }_{nk}}{{\omega }_{nk}}$ may differ from $\frac{{\gamma }_{nk}^{0}}{{\omega }_{nk}}$ by up to one and a half times, which can be recorded in the experiment. It may turn out that for the entire spectrum of the hydrogen atom, ${\omega }_{nk} < \mathrm{\varkappa }\mathrm{c}$. In this case, to register this effect, it is necessary to analyze the damping rate of the spontaneous emission for other atoms emitting at higher frequencies than the hydrogen atom.

7.   Conclusions

We have shown that spontaneous emission and so-called spontaneous transitions, in contrast to linear orthodox quantum mechanics and quantum electrodynamics, in the self-consistent Maxwell-Pauli theory ^[28,29], are a direct and natural consequence of the theory, and their explanation and description do not require additional hypotheses, such as, for example, the existence of a quantum-mechanical vacuum and its fluctuations or Zitterbewegung.

Indeed, if the atom is in a mixed state, then, as follows from Eqs. (112) and (113), as well as from Eqs. (125) and (126), the amplitude of the eigenmode corresponding to a higher frequency decreases, while the amplitude of the eigenmode with a lower frequency increases. In this case, in full accordance with classical electrodynamics, emission of electromagnetic energy occurs. Thus, a mixed state of an atom in which several eigenmodes are excited simultaneously is unstable: the electron wave (and its electric charge) spontaneously flows from the excited mode with a higher eigenfrequency to the excited mode with a lower eigenfrequency. The exceptions are mixed states with ${\omega }_{nk} = 0$ (i.e., with degenerate excited modes) and mixed states with ${\left|{\mathbf{d}}_{nk}\right|}^{2} = 0$ (the so-called "selection rule"). Note that the eigenfrequencies (eigenvalues) of the linear Schrödinger equation (71) are negative. Therefore, the eigenfrequency of the ground mode of the hydrogen atom is lower than the frequency of any other eigenmode of the hydrogen atom. Thus, if a hydrogen atom is in a mixed state in which one of the excited eigenmodes is the ground mode, then such a mixed state of the atom is unstable and the atom will spontaneously transit to the ground state. This means that the ground state of the atom is absolutely stable.

An analysis of the complete system of equations (1)–(9) showed that the compensating short-range $\mathbf{G}$ field plays a significant role in this process, and the constant $\varkappa$ included in the inhomogeneous Klein–Gordon equation (2) determines the dynamics of spontaneous transitions. Thus, if the frequency spectrum of spontaneous radiation ${\omega }_{nk}$ is limited from above, ${\omega }_{nk}\le {\omega }_{max}$ and, in addition, ${\omega }_{max} < \varkappa c$, then the field $\mathbf{G}$ of such an atom will not only quickly (exponentially) decay at distances greater than ${a}_{B}$ but will also not create radiation that could be detected at macroscopic distances from the atom. In this case, macroscopic devices are not capable of directly detecting the field $\mathbf{G}$, and its existence can only be judged by indirect data.

In the case when ${\omega }_{nk} > \varkappa c$, despite the fact that the stationary compensating G field is short-range and decreases exponentially already at distances of the order of ${a}_{B}$ from the center of the atom, G waves can exist that carry away energy and angular momentum ^[28]. Unlike the previous case, this means that in this case, although the stationary G field still cannot be recorded by instruments at macroscopic distances from the source, the G waves emitted by the atom can be detected at macroscopic distances from their source.

If the frequency spectrum of spontaneous emission ${\omega }_{nk}$ is limited from above, ${\omega }_{nk}\le {\omega }_{max}$ and, in addition, ${\omega }_{max} > \varkappa c$, then the atom can emit G waves at frequencies ${\omega }_{nk}$ lying in the range $[\varkappa c, {\omega }_{max}]$ but cannot emit G waves at frequencies ${\omega }_{nk} < \varkappa c$.

Thus, although we cannot detect the stationary G field with macroscopic instruments due to its short-range, the results obtained give us hope that under certain conditions, we will be able to detect the wave G field (G waves) with macroscopic instruments.

To do this, it is necessary to study the emission of G waves and their interaction with matter (atoms), which can be done by searching for the corresponding solutions of equations (1)–(9) or their relativistic analogue ^[29].

This issue will be discussed in subsequent articles.

In conclusion, we note that the theory under consideration is a purely classical field theory. At the same time, we should mention other purely field models for quantum phenomena; see, e.g., ^[34] and references herein.

Use of AI tools declaration

The author declares he has not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

This work was done on the theme of the State Task No. 124012500440-9.

Conflict of interest

The author declares that there is no conflict of interest.