Stochastic quantum mechanics

Last updated

Stochastic mechanics is a framework for describing the dynamics of particles that are subjected to an intrinsic random processes as well as various external forces. The framework provides a derivation of the diffusion equations associated to these stochastic particles. It is best known for its derivation of the Schrödinger equation as the Kolmogorov equation for a certain type of conservative (or unitary) diffusion, [1] [2] and for this purpose it is also referred to as stochastic quantum mechanics.

Contents

The derivation can be based on the extremization of an action in combination with a quantization prescription. [2] This quantization prescription can be compared to canonical quantization and the path integral formulation, and is often referred to as Nelson’s stochastic quantization or stochasticization. [3] As the theory allows for a derivation of the Schrödinger equation, it has given rise to the stochastic interpretation of quantum mechanics. This interpretation has served as the main motivation for developing the theory of stochastic mechanics. [1]

The first relatively coherent stochastic theory of quantum mechanics was put forward by Hungarian physicist Imre Fényes. [4] [5] [6] [7] Louis de Broglie [8] felt compelled to incorporate a stochastic process underlying quantum mechanics to make particles switch from one pilot wave to another. [7] The theory of stochastic mechanics is ascribed to Edward Nelson, who independently discovered a derivation of the Schrödinger equation within this framework. [1] [2] This theory was also developed by Davidson, Guerra, Ruggiero, Pavon and others. [7]

Stochastic interpretation of quantum mechanics

The stochastic interpretation interprets the paths in the path integral formulation of quantum mechanics as the sample paths of a stochastic process. [9] It posits that quantum particles are localized on one of these paths, but observers cannot predict with certainty where the particle is localized. The only way to locate the particle is by performing a measurement. An observer can only predict probabilities for the outcomes of such a measurement based on their earlier measurements and their knowledge about the forces that are acting on the particle.

This interpretation is well-known from the context of statistical mechanics, [9] and Brownian motion in particular. Hence, according to the stochastic interpretation, quantum mechanics should be interpreted in a way similar to Brownian motion. [1] However, in the case of Brownian motion, the existence of a probability measure (called the Wiener measure [10] ) that defines the statistical path integral is well established, and this measure can be generated by a stochastic process called the Wiener process. [11] On the other hand, proving the existence of a probability measure that defines the quantum mechanical path integral faces difficulties, [12] [13] and it is not guaranteed that such a probability measure can be generated by a stochastic process. Stochastic mechanics is the framework concerned with the construction of such stochastic processes that generate a probability measure for quantum mechanics.

For a Brownian motion, it is known that the statistical fluctuations of a Brownian particle are often induced by the interaction of the particle with a large number of microscopic particles. [14] [15] [16] [17] In this case, the description of a Brownian motion in terms of the Wiener process is only used as an approximation, which neglects the dynamics of the individual particles in the background. Instead it describes the influence of these background particles by their statistical behavior.

The stochastic interpretation of quantum mechanics is agnostic about the origin of the quantum fluctuations of a quantum particle. It introduces the quantum fluctuations as the result of new stochastic law of nature called the background hypothesis. [2] This hypothesis can be interpreted as a strict implementation of the statement that `God plays dice’, but it leaves open the possibility that this dice game is replaced by a hidden variable theory, as in the theory of Brownian motion. [18]

The remainder of this article deals with the definition of such a process and the derivation of the diffusion equations associated to this process. This is done in a general setting with Brownian motion and Quantum mechanics as special limits, where one obtains respectively the heat equation and the Schrödinger equation. The derivation heavily relies on tools from Lagrangian mechanics and stochastic calculus.

Stochastic quantization

The postulates of Stochastic Mechanics can be summarized in a stochastic quantization condition that was formulated by Nelson. [2] For a non-relativistic theory on this condition states: [19]

Using the decomposition , and the fact that has finite variation, one finds that the quadratic variation of and is given by

Hence, by Lévy's characterization of Brownian motion, and describe two correlated Wiener processes with a drift described by the finite variation process , a diffusion constant scaling with , and a correlation depending on the angle . The processes are maximally correlated in the quantum limit, associated to and corresponding to , whereas they are uncorrelated in the Brownian limit, associated to and corresponding to ,

The term stochastic quantization to describe this quantization procedure was introduced in the 1970s. [21] Nowadays, stochastic quantization more commonly refers to a framework developed by Parisi and Wu in 1981. Consequently, the quantization procedure developed in stochastic mechanics is sometimes also referred to as Nelson's stochastic quantization or stochasticization. [3]

Velocity of the process

The stochastic process is almost surely nowhere differentiable, such that the velocity along the process is not well-defined. However, there exist velocity fields, defined using conditional expectations. These are given by

and can be associated to the Itô integral along the process . Since the process is not differentiable, these velocities are, in general, not equal to each other. The physical interpretation of this fact is as follows: at any time the particle is subjected to a random force that instantaneously changes its velocity from to . As the two velocity fields are not equal, there is no unique notion of velocity for the process . In fact, any velocity given by

with represents a valid choice for the velocity of the process . This is particularly true for the special case denoted by , which can be associated to the Stratonovich integral along .

Since has a non-vanishing quadratic variation, one can additionally define second order velocity fields given by

The time-reversibility postulate imposes a relation on these two fields such that . Moreover, using the structure relation by which the quadratic variation is fixed, one finds that . It follows that in the Stratonovich formulation the second order part of the velocity vanishes, i.e. .

The real and imaginary part of the velocities are detnoted by

Using the existence of these velocity fields, one can formally define the velocity processes by the Itô integral . Similarly, one can formally define a process by the Stratonovich integral and a second order velocity process by the Stieltjes integral . Using the structure relation, one then finds that the second order velocity process is given by . However, the processes and are not well-defined: the first moments exist and are given by , but the quadratic moments diverge, i.e. . The physical interpretation of this divergence is that in the position representation the position is known precisely, but the velocity has an infinite uncertainty.

Stochastic action

The stochastic quantization condition states that the stochastic trajectory must extremize a stochastic action , but does not specify the stochastic Lagrangian . This Lagrangian can be obtained from a classical Lagrangian using a standard procedure. Here, we consider a classical Lagrangian of the form

Here, are coordinates in phase space (the tangent bundle), is the Kronecker delta describing the metric on , denotes the mass of the particle, the charge under the vector potential , and is a scalar potential. Moreover, Einstein's summation convention is assumed.

An important property of this Lagrangian is the principle of gauge invariance. This can be made explicit by defining a new action through the addition of a total derivative term to the original action, such that

where and . Thus, since the dynamics should not be affected by the addition of a total derivative to the action, the action is gauge invariant under the above redefinition of the potentials for an arbitrary differentiable function .

In order to construct a stochastic Lagrangian corresponding to this classical Lagrangian, one must look for a minimal extension of the above Lagrangian that respects this gauge invariance. [22] In the Stratonovich formulation of the theory, this can be done straightforwardly, since the differential operator in the Stratonovich formulation is given by

Therefore, the Stratonovich Lagrangian can be obtained by replacing the classical velocity by the complex velocity , such that

In the Itô formulation, things are more complicated, as the total derivative is given by Itô's lemma:

Due to the presence of the second order derivative term, the gauge invariance is broken. However, this can be restored by adding a derivative of the vector potential to the Lagrangian. Hence, the stochastic Lagrangian is given by a Lagrangian of the form

The stochastic action can be defined using the Stratonovich Lagrangian, which is equal to the action defined by the Itô Lagrangian up to a divergent term: [2]

The divergent term can be calculated and is given by

where are winding numbers that count the winding of the path around the pole at . [23]

As the divergent term is constant, it does not contribute to the equations of motion. For this reason, this term has been discarded in early works on stochastic mechanics. [2] However, when this term is discarded, stochastic mechanics cannot account for the appearance of discrete spectra in quantum mechanics. This issue is known as Wallstrom's criticism, [24] [25] and can be resolved by properly taking into account the divergent term.

There also exists a Hamiltonian formulation of stochastic mechanics. [22] [26] It starts from the definition of canonical momenta:

The Hamiltonian in the Stratonovich formulation can then be obtained by the first order Legendre transform:

In the Itô formulation, on the other hand, the Hamiltonian is obtained through a second order Legendre transform: [27]

Euler-Lagrange equations

The stochastic action can be extremized, which leads to a stochastic version of the Euler-Lagrange equations. In the Stratonovich formulation, these are given by

For the Lagrangian, discussed in previous section, this leads to the following second order stochastic differential equation in the sense of Stratonovich:

where, the field strength is given by . This equation serves as a stochastic version of Newton's second law.

In the Itô formulation, the stochastic Euler-Lagrange equations are given by

This leads to a second order stochastic differential equation in the sense of Itô, given by a stochastic version of Newton's second law in the form

Hamilton-Jacobi equations

The equations of motion can also be obtained in a stochastic generalization of the Hamilton-Jacobi formulation of classical mechanics. In this case, one starts by defining Hamilton's principal function. For the Lagrangian , this function is defined as

where it is assumed that the process obeys the stochastic Euler-Lagrange equations. Similarly, for the Lagrangian , Hamilton's principal function is defined as

where it is assumed that the process obeys the stochastic Euler-Lagrange equations. Due to the divergent part of the action, these principal functions are subjected to the equivalence relation

By varying the principal functions with respect to the point one finds the Hamilton-Jacobi equations. These are given by

Note that these look the same as in the classical case. However, the Hamiltonian, in the second Hamilton-Jacobi equation is now obtained using a second order Legendre transform. Moreover, due to the divergent part of the action, there is a third Hamilton-Jacobi equation, which takes the form of the non-trivial integral constraint

For the given Lagrangian the first two Hamilton-Jacobi equations yield

These two equations can be combined, yielding

Using that , this equation, subjected to the integral condition and the initial condition or terminal condition , can be solved for . The solution can then be plugged into the Itô equation

which can be solved for the process . Thus, when an initial condition (for the future directed equation labeled with ) or terminal condition (for the past directed equation labeled with ) is specified, one finds a unique stochastic process that describes the trajectory of the particle.

Diffusion Equation

The key result of stochastic mechanics is that it derives the Schrödinger equation from the postulated stochastic process. In this derivation, the Hamilton-Jacobi equations

are combined, such that one obtains the equation

Subsequently, one defines the wave function

Since Hamilton's principal functions are multivalued, one finds that the wave functions are subjected to the equivalence relations

Furthermore, the wave functions are subjected to the complex diffusion equations

Thus, for any for any process that solves the postulates of stochastic mechanics, one can construct a wave function that obeys these diffusion equations. Due to the equivalence relations on Hamilton's principal function, the opposite statement is also true: for any solution of these complex diffusion equations, one can construct a stochastic process that is a solution of the postulates of stochastic mechanics. A similar result has been established by the Feynman-Kac theorem.

Finally, one can construct a probability density

which describes transition probabilities for the process . More precisely, describes the probability of being in the state given that the system ends up in the state . Therefore, the diffusion equation for can be interpreted as the Kolmogorov backward equation of the process . Similarly, describes the probability of being in the state given that the system ends up in the state , when it is evolved backward in time. Therefore, the diffusion equation for can be interpreted as the Kolmogorov backward equation of the process when it is evolved towards the past. By inverting the time direction, one finds that describes the probability of being in the state given that the system starts in the state , when it is evolved forward in time. Thus, the diffusion equation for can also be interpreted as the Kolmogorov Forward equation of the process when it is evolved towards the future.

Mathematical aspects

Limiting cases

The theory contains various special limits:

In the Brownian limit with initial condition or terminal condition (which implies ), the processes and are decoupled, such that the dynamics of the auxiliary process can be discarded, and is described by a real Wiener process. In all other cases with , the processes are coupled to each other, such that the auxiliary process must be taken into account in deriving the dynamics of .

Time-reversal symmetry

The theory is symmetric under the time reversal operation .

In the Brownian limits, the theory is maximally dissipative, whereas the quantum limits are unitary, such that

Canonical commutation relations

The diffusion equation can be rewritten as

where is a Hamiltonian operator. This allows to introduce position and momentum operators as

such that the Hamiltonian has its familiar shape

These operators obey the canonical commutation relation

See also

Notes

    1. 1 2 3 4 E. Nelson 1966.
    2. 1 2 3 4 5 6 7 E. Nelson 1985.
    3. 1 2 E. Nelson 2014.
    4. I. Fenyes 1946.
    5. I. Fenyes 1948, pp. 336–346.
    6. M.P. Davidson 1979, p. 1.
    7. 1 2 3 L. de la Peña & A.M. Cetto 1996, p. 36.
    8. L. de Broglie 1967.
    9. 1 2 F. Guerra 1981.
    10. Wiener 1923.
    11. M. Kac 1949.
    12. R.H. Cameron 1960, pp. 126–140.
    13. Yu.L. Daletskii 1962, pp. 1–107.
    14. A. Einstein 1905.
    15. J. Perrin 1911.
    16. J.L. Doob 1942.
    17. E. Nelson 1967.
    18. F. Kuipers, 2023 & p. 61.
    19. F. Kuipers 2023, p. 9.
    20. F. Kuipers 2023, p. 23.
    21. K. Yasue 1979.
    22. 1 2 J.C. Zambrini 1985.
    23. F. Kuipers 2023, p. 34.
    24. T.C. Wallstrom 1989.
    25. T.C. Wallstrom 1994.
    26. M. Pavon 1995.
    27. Q. Huang & J.C. Zambrini 2023.

    Related Research Articles

    Continuum mechanics is a branch of mechanics that deals with the deformation of and transmission of forces through materials modeled as a continuous medium rather than as discrete particles. The French mathematician Augustin-Louis Cauchy was the first to formulate such models in the 19th century.

    <span class="mw-page-title-main">Heat equation</span> Partial differential equation describing the evolution of temperature in a region

    In mathematics and physics, the heat equation is a certain partial differential equation. Solutions of the heat equation are sometimes known as caloric functions. The theory of the heat equation was first developed by Joseph Fourier in 1822 for the purpose of modeling how a quantity such as heat diffuses through a given region.

    <span class="mw-page-title-main">Fokker–Planck equation</span> Partial differential equation

    In statistical mechanics and information theory, the Fokker–Planck equation is a partial differential equation that describes the time evolution of the probability density function of the velocity of a particle under the influence of drag forces and random forces, as in Brownian motion. The equation can be generalized to other observables as well. The Fokker-Planck equation has multiple applications in information theory, graph theory, data science, finance, economics etc.

    <span class="mw-page-title-main">Poisson bracket</span> Operation in Hamiltonian mechanics

    In mathematics and classical mechanics, the Poisson bracket is an important binary operation in Hamiltonian mechanics, playing a central role in Hamilton's equations of motion, which govern the time evolution of a Hamiltonian dynamical system. The Poisson bracket also distinguishes a certain class of coordinate transformations, called canonical transformations, which map canonical coordinate systems into canonical coordinate systems. A "canonical coordinate system" consists of canonical position and momentum variables that satisfy canonical Poisson bracket relations. The set of possible canonical transformations is always very rich. For instance, it is often possible to choose the Hamiltonian itself as one of the new canonical momentum coordinates.

    In mathematics, the Hodge star operator or Hodge star is a linear map defined on the exterior algebra of a finite-dimensional oriented vector space endowed with a nondegenerate symmetric bilinear form. Applying the operator to an element of the algebra produces the Hodge dual of the element. This map was introduced by W. V. D. Hodge.

    <span class="mw-page-title-main">Path integral formulation</span> Formulation of quantum mechanics

    The path integral formulation is a description in quantum mechanics that generalizes the stationary action principle of classical mechanics. It replaces the classical notion of a single, unique classical trajectory for a system with a sum, or functional integral, over an infinity of quantum-mechanically possible trajectories to compute a quantum amplitude.

    In Hamiltonian mechanics, a canonical transformation is a change of canonical coordinates (q, p) → that preserves the form of Hamilton's equations. This is sometimes known as form invariance. Although Hamilton's equations are preserved, it need not preserve the explicit form of the Hamiltonian itself. Canonical transformations are useful in their own right, and also form the basis for the Hamilton–Jacobi equations and Liouville's theorem.

    Creation operators and annihilation operators are mathematical operators that have widespread applications in quantum mechanics, notably in the study of quantum harmonic oscillators and many-particle systems. An annihilation operator lowers the number of particles in a given state by one. A creation operator increases the number of particles in a given state by one, and it is the adjoint of the annihilation operator. In many subfields of physics and chemistry, the use of these operators instead of wavefunctions is known as second quantization. They were introduced by Paul Dirac.

    In quantum mechanics, the canonical commutation relation is the fundamental relation between canonical conjugate quantities. For example,

    In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics.

    In theoretical physics, supersymmetric quantum mechanics is an area of research where supersymmetry are applied to the simpler setting of plain quantum mechanics, rather than quantum field theory. Supersymmetric quantum mechanics has found applications outside of high-energy physics, such as providing new methods to solve quantum mechanical problems, providing useful extensions to the WKB approximation, and statistical mechanics.

    In differential geometry, the four-gradient is the four-vector analogue of the gradient from vector calculus.

    In quantum mechanics, the probability current is a mathematical quantity describing the flow of probability. Specifically, if one thinks of probability as a heterogeneous fluid, then the probability current is the rate of flow of this fluid. It is a real vector that changes with space and time. Probability currents are analogous to mass currents in hydrodynamics and electric currents in electromagnetism. As in those fields, the probability current is related to the probability density function via a continuity equation. The probability current is invariant under gauge transformation.

    In atomic, molecular, and optical physics and quantum chemistry, the molecular Hamiltonian is the Hamiltonian operator representing the energy of the electrons and nuclei in a molecule. This operator and the associated Schrödinger equation play a central role in computational chemistry and physics for computing properties of molecules and aggregates of molecules, such as thermal conductivity, specific heat, electrical conductivity, optical, and magnetic properties, and reactivity.

    In classical mechanics, holonomic constraints are relations between the position variables that can be expressed in the following form:

    In classical mechanics, Appell's equation of motion is an alternative general formulation of classical mechanics described by Josiah Willard Gibbs in 1879 and Paul Émile Appell in 1900.

    In mathematics — specifically, in stochastic analysis — the infinitesimal generator of a Feller process is a Fourier multiplier operator that encodes a great deal of information about the process.

    In pure and applied mathematics, quantum mechanics and computer graphics, a tensor operator generalizes the notion of operators which are scalars and vectors. A special class of these are spherical tensor operators which apply the notion of the spherical basis and spherical harmonics. The spherical basis closely relates to the description of angular momentum in quantum mechanics and spherical harmonic functions. The coordinate-free generalization of a tensor operator is known as a representation operator.

    The Landau kinetic equation is a transport equation of weakly coupled charged particles performing Coulomb collisions in a plasma.

    Projection filters are a set of algorithms based on stochastic analysis and information geometry, or the differential geometric approach to statistics, used to find approximate solutions for filtering problems for nonlinear state-space systems. The filtering problem consists of estimating the unobserved signal of a random dynamical system from partial noisy observations of the signal. The objective is computing the probability distribution of the signal conditional on the history of the noise-perturbed observations. This distribution allows for calculations of all statistics of the signal given the history of observations. If this distribution has a density, the density satisfies specific stochastic partial differential equations (SPDEs) called Kushner-Stratonovich equation, or Zakai equation. It is known that the nonlinear filter density evolves in an infinite dimensional function space.

    References

    Papers

    Books