1. Introduction
A recent paper [
1] discussed the formal similarities and differences between the following diffusion (or heat) equation
which describes the spatial and temporal change of a
real distribution function
due to a random process—where
D is a
real diffusion coefficient and
is the Laplace-operator—and the following time-dependent Schrödinger equation (particularly for the free motion, i.e., potential,
)
which describes the same evolution of a
complex state function
whose absolute value squared, i.e.,
, is a quantum mechanical probability distribution and
(with
and the reduced Planck constant
) now being a purely
imaginary coefficient. Despite the formal similarities, both equations describe very different processes.
The solutions
of Equation (
2) are also called wave functions because they behave comparably to solutions of the following wave equation:
for a wave
that shows properties such as diffraction and dispersion, where
n is the diffraction index and
c is the speed of light in a vacuum. In contrast with Equation (
1), this equation is invariant under time reversal as
t appears quadratic on the lhs of (
3).
So, although mathematically the Schrödinger equation has a structure similar to that of (
1)—i.e., a parabolic differential equation—its solutions behave similarly to those of (
3): a hyperbolic differential equation with a second derivative in time. The difference between Equations (
1) and (
2) is that
is in general a
complex function and the imaginary unit
appears explicitly in Equation (
2). (These fundamental differences and their mathematical and physical consequences were already pointed out explicitly by Schrödinger himself in [
2]—see also [
3] for an English translation and comments.) However, this is an important—if not the most important—difference between classical and quantum physics, as stressed by C. N. Yang in their lecture on the occasion of Schrödinger’s 100th anniversary (see [
4]), since “complex numbers become a conceptual element of the very foundations of physics”. This is reflected by the fact that
occurs explicitly in the fundamental equations of quantum mechanics not only in Schrödinger’s Equation (
2), but also in Heisenberg’s commutation relation
in their matrix mechanics (with
position and
momentum operators).
There are attempts in the literature to turn a complex space into a real space with twice as many dimensions, arguing that complex numbers are just a convenient mathematical tool, but not essential for the quantum mechanical description of nature. However, this is not true, as real and imaginary parts, or amplitude and phase of the complex wave function, are not independent of each other, but are uniquely coupled. This can be seen clearly in the hydrodynamic formulation (see
Section 2) where the continuity Equation (
6) for the density,
, and thus for the amplitude,
, contains the phase of the wave function via
; additionally, this is shown in Equation (
8) where the phase,
S, contains the amplitude via
. This is not taken into account by doubling the number of degrees of freedom. The viewpoint that complex quantities are essential for quantum mechanics has also recently been confirmed experimentally (see, e.g., [
5,
6]).
Even though Heisenberg and Schrödinger arrived from different starting points, they share the same guideline: the relation between classical and quantum mechanics should be mediated via the mean value which is the principle of correspondence [
7] and mathematically expressed by the so-called Ehrenfest theorem [
8]. As Heisenberg pointed out, the difference lies in the treatment of kinematic quantities, such as position and velocity, while Heisenberg associates them with time-dependent complex matrices; for Schrödinger, the time-dependence was contained in the complex wave function. However, in the beginning, Schrödinger did not consider the imaginary unit in their first communication [
9] where, in the definition of their wave function via the action
S according to
, the constant
K is chosen to be
and the wave function itself is assumed to be real. It was after some correspondence with Lorentz (see, e.g., [
10]), he had to admit that their wave function is complex. So, in [
11]—published a few months later—he uses a complex form for the wave function in Equation (
2), and below Equation (
4) he explicitly uses …“komplexe Wellenfunktion
”.
The importance of
in quantum mechanics and Schrödinger’s struggle with complex quantities in wave mechanics has also been discussed recently in the literature [
12,
13,
14]. It is not only the wave function, but also the operators—essential to Schrödinger’s wave mechanical form of quantum mechanics—which can be complex or imaginary, such as the momentum operator in position representation,
(with
Nabla operator).
What about the time reversal invariance of the Schrödinger Equation (
2)? As pointed out, even recently in several articles on this subject [
15,
16,
17,
18,
19], a simple replacement of
t by
would change the form of (
2) as in the case of the diffusion equation. Time reversal in classical physics also means that velocity or momentum change sign (magnetic fields are not included in our discussion), but the time,
t, does not occur explicitly in the quantum mechanical operator,
; however, a change of sign can be achieved by taking the conjugate complex of this operator. This is sometimes used as an explanation for the so-called Wigner time-reversal invariance. (A more precise formulation requires the time inversion operator to be anti-linear. Wigner actually used more appropriately the terminology “reversal of motion” instead of “time reversal” [
20].) In other words, in quantum mechanics, it is not only
t that has to be changed to
, but the conjugate complex also has to be taken; meaning
and
. In this case, the form of Equation (
2) remains invariant under this transformation.
However, considering the latter condition for the time reversal, it is actually only the phase of the complex wave function—when written in polar form as —that is affected by this transformation, not the amplitude, , providing that t occurs in in an even power. Does this mean that the time-dependent Schrödinger equation contains an element of broken time-reversal invariance as the (real) diffusion equation does? If so, how can this formally be seen in the equations?
To answer these questions, we consider the Gaussian wave packet solution of Equation (
2) (in one dimension), where the probability density,
, also Gaussian, can be written in the following form:
where
is the mean value of position and
describes the position uncertainty and, in this case, the width of the Gaussian function. It is well known that a Gaussian function is completely determined by its maximum—here
(position mean value)—and its width—here
—where, in our case, both parameters can be time-dependent. In the case of the free motion,
, the width is spreading quadratically in time according to
Let us consider an initial wave packet with its maximum at that follows the classical trajectory, moving with the (constant) velocity and having the initial width . After a time interval, , the maximum has reached the following position: . Now, applying time reversal, the maximum fulfills that which Lohschmidt described in their argument against Boltzmann’s -function, i.e., the system will go (in configuration space) though the same states as before, only in reverse order. So, after another time interval, , the maximum of the Gaussian will reach the following initial position: .
However, what happens to the width? Starting the time-reversed evolution with the width, , the Gaussian wave packet will continue spreading according to and will arrive after a second period, , with an even broader width, at the initial position—i.e., the wave packet will not shrink to its initial width. (It should be pointed out that this continuation of the spreading is not connected with any perturbance due to the interaction with an environment that might be assumed to take place in order to reverse the motion. The time reversal is purely a “Gedankenexperiment”, where no interaction with any environment is taken into account, only the Schrödinger equation for an isolated system is considered.) !
Obviously, the amplitude of the wave function, and thus the density, , behaves differently from the phase of the wave function that is affected by the time reversal, similar to velocity in classical physics. Can this different behavior also somehow be expressed in the corresponding equations of motion for amplitude and phase—bearing in mind that Heisenberg’s kinematic quantities, in particular velocity, are to be considered complex?
For this purpose, we apply Madelung’s hydrodynamic formulation of quantum mechanics [
21] that allows us to separate the time-dependent Schrödinger equation into two (coupled) equations: one continuity equation for the amplitude and a modified Hamilton–Jacobi equation for the phase. This corresponds to an Euler-type equation for a velocity field. However, the continuity equation does not contain any irreversible diffusion current. In the following, this problem is solved by introducing complex kinematic quantities, inspired by Heisenberg’s original idea. This naturally causes the appearance of a diffusive current in the continuity equation with an already known diffusion coefficient. The reason that the diffusive current is absent in the Madelung approach is that it cancels with the imaginary part of the now-complex velocity field.
The paper is organized as follows: In
Section 2, a short overview of the current hydrodynamic formulation of quantum mechanics is presented and the inconsistencies related to the diffusion process are specified. To solve this problem, in
Section 3.1, a complex notation of kinematic quantities is introduced and applied in
Section 3.2 to obtain the missing diffusion term in the continuity equation. Using the Gaussian wave packet solutions of the time-dependent Schrödinger equation as example, in
Section 3.3, the relation between the real and imaginary parts of the introduced complex velocity field and phase and amplitude of the wave function is demonstrated.
Section 3.4 shows that the introduction of the complex velocity field does not change the form of the Euler equation for the phase of the wave packet, in agreement with the correct evolution of the wave packet maximum. Finally, in
Section 4, the results are summarized and our conclusions are drawn.
2. Conventional Quantum Hydrodynamics
Shortly after Schrödinger published their papers on wave mechanics, Madelung showed [
21] that the
complex Schrödinger equation is equivalent to a set of
two real hydrodynamic equations.
A continuity equation for the density,
, has the following form:
where
denotes a velocity field related with the action
S in the phase of the wave function via
The second hydrodynamic equation is obtained from the following evolution of the phase:
that is a modified Hamilton–Jacobi Equation. (Note that the real Equations (
6) and (
8) are invariant under time reversal, if
t is replaced by
and
S, by
, leading to
using (
7).) Taking its gradient yields
showing formal similarity with the Euler equation
with
being the force density per mass and
P the thermodynamic pressure. Comparison suggests the connection between the so-called “quantum potential”
and the pressure
P via
As Takabayasi points out in [
22], there should be a stress tensor rather than a pressure term, what makes sense in particular for
(such as for a spreading wave packet in quantum mechanics), as in this case these two quantities are not identical. However, in our quantum mechanical case, a proper explanation of
in terms of the complex quantities—introduced in the next Section—will be given.
Madelung’s formulation of wave mechanics was also later independently used by David Bohm in their deterministically inspired version of quantum mechanics [
23,
24] where he claimed the existence of real paths of (quantum) particles that can be obtained by integration of Equation (
7). As we have shown recently [
25], this deterministic viewpoint is incorrect and has to be replaced by a probabilistic one that differs from the usual probabilistic viewpoint taken in quantum mechanics. Nevertheless, Bohmian mechanics can still be helpful in the treatment of quantum systems, e.g., tunneling problems, particularly when performing numerical simulations [
26,
27].
Looking at Madelung’s fluid as a classical one, it can be associated with an irrotational, inviscid and compressible one with no obvious quantum effects at first sight. One may argue that these effects are embedded in the term related to the “quantum potential”,
. However, apart from this rather misleading terminology (
actually originates from the kinetic energy operator, as will be shown below, and is not a potential at all), the corresponding term has still a more classical meaning, only indicating that the fluid is compressible, since
depends on the density. As mentioned above, there is the well-known diffusive effect of the wave packet spreading; however, there is no diffusive current in the continuity Equation (
6), only a convective one.
On the other hand, from a fundamental point of view, it should be expected that the (kinetic) quantities in a hydrodynamic formulation of quantum mechanics would be of a complex nature, unlike the classical counterpart. As Heisenberg [
28] wrote in their pioneering paper of 1925:
“ The Einstein–Bohr frequency condition (which is valid in all cases) already presents such a complete departure from classical mechanics, or rather (using the viewpoint of wave theory) from the kinematics underlying this mechanics, that even for the simplest quantum theoretical problems the velocity of classical mechanics simply cannot be maintained.”
This was a milestone for future quantization methods—i.e., mappings onto classical-like equations: only kinematic quantities should be of a different nature than their classical counterparts. From this viewpoint, the mapping of a quantum system onto hydrodynamical equations should be of the following form:
where the flow,
, and the density current,
, should be of a
complex nature, since their origin is kinematic, whilst
, not being of a kinematic quantity, must remain a positive real quantity.
It is important to emphasize that we are not claiming that the hydrodynamical Equations (
6) and (
9) are inaccurate, but rather that they are only one perspective of the quantization principle. Therefore, the explanation of the diffusive effects of wave packets is unanswered within the current hydrodynamic formulation due to the absence of a diffusive current in the continuity equation. We will now show how this can be taken into account in a complex formulation.
4. Conclusions
In classical mechanics, time reversal means the following replacements: , or . According to Loschmidt’s arguments, if a system reaches a certain position in configuration space after a certain period of time, this system will return to its initial position after the same period of time if time reversal is applied according to the abovementioned rules.
Concerning quantum mechanics, the situation is different as time does not appear explicitly in the operator for momentum or velocity, rather it appears as the imaginary unit
. This is related to a major difference between classical and quantum mechanics. According to Wigner [
20], one should not talk about “time reversal” in quantum mechanics but about “motion reversal”. The latter, however, is not simply the replacement of
t by
, but also involves the application of complex conjugation. This, in turn, affects not just the operators that can contain the imaginary unit, i, but also the complex wave function in the Schrödinger equation. Writing this wave function in polar form,
—consistent with Madelung—shows that complex conjugation only affects the phase of the wave function, not its amplitude,
.
Looking at the time evolution of the wave packet solution of the time-dependent Schrödinger equation, specifically for free motion, , it is obvious that motion reversal in the abovementioned way does not revert to the initial state after a second period of time evolution. In agreement with the classical trajectory, the maximum of the wave packet reaches its initial position after motion reversal. However, the width does not shrink to its initial one but becomes much broader due to continuing dispersion—such as a diffusion effect—also in the time-reversed segment of the evolution. This irreversible effect should somehow also be manifested in the equation of motion for density, .
According to Madelung [
21], the complex Schrödinger equation can be rewritten as two (coupled) real equations; one for the phase, the other for the amplitude of the wave function. Using the classical definition of time reversal as mentioned above, both equations seem invariant under this operation. However, the typical quantum mechanical aspect of complex quantities is eliminated in this formulation.
Following Heisenberg, in the quantum theory kinematic quantities, such as velocity, should be replaced by complex ones. Proceeding accordingly in Madelung’s (continuity) equation for the density, replacing the real velocity,
, with the complex,
—as defined in (
14)—this equation can only be fulfilled if a diffusion term (with imaginary diffusion coefficient) is added to compensate for the contribution from the imaginary part,
, of the complex velocity, leading to a Smoluchowski-type equation. This additional diffusion term expresses, at least formally, the spreading behavior of the wave packet.
In Madelung’s second equation for the phase, the replacement of by complex does not change the form of the corresponding Euler equation; this demonstrates that the Wigner “time reversal” also affects the phase of the wave function and therefore leads to a consistent treatment. This is supported by the fact that the maximum of the wave packet—related to the mean value of position following the classical trajectory—returns to its initial position after time inversion and complex conjugation. Complex conjugation of the phase changes the sign of the classical velocity (or momentum) that appears in the phase, thus having the same effect as time reversal in the classical case.
In conclusion, the irreversible behavior of the quantum mechanical wave packet observed under motion reversal formally becomes obvious only if the complex kinematic quantities are applied.
In this sense, Equations (
21) and (
37) constitute the correct mapping from a quantum system onto a hydrodynamic formulation showing the effects involved more clearly. That does not mean that Madelung’s Equations (
6) and (
9) are incorrect but rather that they display only one aspect of the quantization principle.
Though the mean value of the imaginary part of the complex velocity vanishes, it does not that mean
has no further influence on the quantum system. This becomes obvious when considering the uncertainty product of position and momentum. As has been show in (
32), using (
29), the minimum uncertainty
depends entirely on
and its derivative and can be expressed via the absolute value of the (imaginary) diffusion coefficient as
.
There are certainly other currents related to quantum mechanics that could be considered from the viewpoint of our approach, but this shall not be the subject of this work.
In this paper, only non-relativistic quantum mechanics and the corresponding currents are considered. Interesting questions might be the following: could the asymmetry in the spreading behavior be related to the absence of Lorentz symmetry? Could a relativistic formulation could resolve this problem? This would require the consideration of Dirac currents and, maybe, a transition from complex numbers to quaternions; this is quite beyond the scope of this work, but might be an interesting target for future investigations.