Recursive formulation of Madelung continuity equation leads to propagation equation D. M. Heim1, 2, a) 1)Fraunhofer Institute for Industrial Mathematics, 67663 Kaiserslautern, Germany 2)Key Lab of Information of Yunnan Province, Yunnan University, Kunming, Yunnan, 650091 China

(Dated: 27 December 2018)

We apply a recursive approach to the continuity equation of the Madelung fluid resulting in a propagation equation for particle probability densities. This propaga- tion equation can be used to propagate particle distributions in the presence of a Madelung pressure field. We show that the derived propagation equation goes over into the guidance equation of the de Broglie-Bohm theory in the limit of well located single particles. As an example, we propagate particles that enter the lower slit of a double-slit experiment, while the Madelung fluid enters both slits. arXiv:1611.07518v4 [quant-ph] 26 Dec 2018

a)Electronic mail: [email protected]

1 I. INTRODUCTION

The concept of a hydrodynamical model of was first proposed in 1926 by Madelung1, followed by the proposal of Korn in 19272. This concept revealed to be very fruitful in a number of applications and is now established in several branches of quantum mechanics3, such as Bose-Einstein-condensation4,5, condensed matter physics6,7 and quantum cosmology8,9. It has been also applied as a useful tool to solve linear and non-linear partial differential equations10–13. We extend the Madelung fluid description by applying a recursive approach to its continuity equation which leads to a propagation equation. The Madelung equations arise when the quantum mechanical in polar decomposition √ ψ ≡ P eiS (1) is inserted into the Schr¨odingerequation

∂  2  i ψ = − ~ ∇2 + V ψ (2) ~∂t 2m

and the imaginary part

∂ P = − ~ ∇(P ∇S) (3) ∂t m

and the real part √ 2 " 2 # ∂ ~ ∇ P 2 ~ S = √ − (∇S) − V (4) ∂t 2m P

are separated. Here, we used the probability density P = P (r, t) to find a particle of mass m at the position r and time t, and the phase S = S(r, t). By defining a velocity field

v = v(r, t) ≡ ~ ∇S, (5) m

the imaginary part (3) can be written in the form of a continuity equation

∂ P = −∇(v P ) (6) ∂t

2 while the gradient of the real part (4) ∂ 1 v + v∇v = − ∇ (V + Q) (7) ∂t m is similar to the Euler equation of classical hydrodynamics14. Here, we defined the quantum pressure field √ 2 ∇2 P Q ≡ − ~ √ . (8) 2m P Consequently, we obtain a set of equations which describe the motion of an inviscid fluid and can be used to predict the outcome of non-relativistic quantum experiments3. We extend the Madelung fluid description by re-writing its continuity equation (6) as a propagation equation

P = Gb[v] P0, (9) which only involves an initial probability density P0 ≡ P (r, t0) and the velocity field (5) as argument of the propagator Gb, which we derive in this publication. The propagation equation (9) is useful to describe the motion of particles in the presence of a Madelung pressure field. We show that this equation interestingly goes over into the guidance equation of the de Broglie-Bohm theory15 for single particles. As an example, we apply the propagation equation in the case of a double-slit experiment, which illustrates the connections of our formalism with quantum mechanics and the de Broglie-Bohm theory.

II. PROPAGATION EQUATION

In order to derive the propagation equation (9), we integrate the continuity equation (6), which leads us with the abbreviations vn ≡ v(r, tn) and Pn ≡ P (r, tn) to

t Z P = P0 − dt1 ∇(v1 P1). (10)

t0 Similar to deriving a formal solution of the Schr¨odingerequation16, we reinsert equa- tion (10) multiple times into the right-hand-side of itself which leads to the series

t t t Z Z Z 1 P = P0 − dt1 ∇ (v1 P0) + dt1 ∇v1 dt2 ∇ (v2 P0) − ... (11)

t0 t0 t0

3 and finally to the propagation equation (9), that is

P = Gb[v] P0, (12) where we defined the propagator   Z t  ˆ Gb[v] ≡ T exp − dt1 ∇v1 . (13) t0 The Tˆ is the well-known time-ordering operator, which is defined by  h i Aˆ(t) Bˆ(t0) for t ≥ t0 Tˆ Aˆ(t) Bˆ(t0) ≡ . (14) Bˆ(t0) Aˆ(t) for t0 ≤ t

Equation (12) describes the propagation of an initial probability density P0 to the final density P by using the velocity field v.

III. PROPAGATION OF SINGLE PARTICLE

In this section, we show that the propagation equation (12) goes over into the guidance equation of the de Broglie-Bohm theory in the limit of well located single particles. The de Broglie-Bohm theory was originally proposed by de Broglie in 192717 and refined by Bohm in 195215,18. This theory supplements quantum mechanics by classical particle trajectories. This is done by regarding the real part (4) of the Schr¨odingerequation as Hamilton-Jacobi equation ∂ (∇ S)2 S + ~ + V + Q = 0, (15) ∂t ~ 2m where Q is called the and ~S the Hamilton-Jacobi function. The trajec- tory ξ(t) of a single particle is then defined in analogy to the Hamilton-Jacobi formalism by solving the first-order differential equation dξ(t) 1 = ∇ S, (16) dt m ~ which is known as the guidance equation19. To reproduce this equation from our formalism, we replace the probability density P in the propagation equation (12) by the Dirac delta function

P → δ[r − ξ(t)], (17)

4 where ξ(t) represents the trajectory of a single particle. After multiplying the resulting equation by r and integrating over all space we obtain Z ξ(t) = dV r Gb[v] δ[r − ξ(t0)]. (18)

To rewrite the probability density Gb[v] δ[r − ξ(t0)] we perform the steps from equation (12) to equation (10) in reverse which leads us to

t Z Z ξ(t) = ξ(t0) − dt1 dV r ∇ {v1 δ[r − ξ(t1)]} . (19)

t0

Here, we perform a partial integration to obtain

t Z ξ(t) = ξ(t0) + dt1 v[ξ(t1), t1]. (20)

t0

By inserting the velocity field (5) into (20) and taking the derivative with respect to t we arrive at the guidance equation (16) of the de Broglie-Bohm theory. We compare the application of the propagation equation (12) with (20) in the next section.

IV. PROPAGATION OF PARTICLE DISTRIBUTION

In the following, we propagate particle probability densities in a double-slit experiment which we sketched in Fig.1. We compare the quantum mechanical probability density with trajectories arising from the Broglie-Bohm theory and with a probability density resulting from a thought experiment, where we consider a Madelung fluid that enters both slits while the particles pass only the lower slit.

A. Quantum mechanics

The wave function

ψ (y, t) + ψ (y, t) ψ(y, t) = u √ l , (21) 2

which describes the possibility of finding particles behind the slits in Fig.1, consists of the

sum of the wave functions ψu(y, t) passing the upper and ψl(y, t) passing the lower slit.

5 FIG. 1. Double-slit experiment, where particles can pass two slits denoted by u and l with a distance 2Y . Each slit has a width σ and distance d from the screen where the particles impinge.

We represent these wave functions by Gaussians

1  (y ∓ Y )2  ψ (y, t) = √ exp − (22) u,l 4 p 2 2π σ + i~t/(2mσ) 4σ + 2i~t/m which solve the Schr¨odinger equation (2) for V = 0. Here, 2Y is the slit separation and σ describes how much the probability density expands in a single slit. In Fig.2(a), we used the wave function (21) to depict the quantum mechanical probability density |ψ(y, t)|2 for parameters similar to the ones of the experiment performed by J¨onsson et al.20. The slit separation is taken to be 2Y = 1 µm and σ = 0.1 µm. After a short propagation time t, the probability density displays the well known interference fringes.

B. De Broglie-Bohm

To calculate de Broglie-Bohm trajectories, we select 30 positions yj(0) for the time t0 = 0 at the location of the lower slit in Fig.1. We propagate each of these positions by using equation (20) in the form

t Z yj(t) = yj(0) + dt1 v[y(t1), t1]. (23) 0

Here, we used the wave function (21) in the velocity field (5) which leads us to

∂ v(y, t) = ~ arg[ψ(y, t)]. (24) m ∂y

The resulting trajectories are shown in Fig.2(b) for the same parameters as we used in Fig.2(a).

6 FIG. 2. Application of different formalisms in the double-slit setup shown in Fig.1. In (a) we depicted the quantum mechanical probability density |ψ(y, t)|2 of the wave function (21), (b) shows de Broglie-Bohm trajectories calculated from equation (23) and the probability density in (c) is calculated by using the propagation equation (26). From (c) one can see that this propagation equation can describe continuous probability densities as in quantum mechanics (a), but also offers the possibility to propagate only parts of the initial probability density as it can be done in the de Broglie-Bohm theory (b).

As already reported by Bohm and Hiley21, the trajectories in such a double-slit experiment are arranged along the valleys of the interference pattern. Differently from their publication, we selected only trajectories emerging from the lower slit in order to stress the similarity between our formalism and the de Broglie-Bohm theory which we discuss in more detail in the next section. However, to reproduce experimental results, the distribution of initial positions should match the density |ψ(y, 0)|2 which covers both slits. C. Madelung fluid

For the calculation of the probability density resulting from a Madelung pressure field

(8), we use the wave function ψl(y, t) defined in (22) to describe an initial probability density

2 ρ0 ≡ |ψl(y, 0)| (25) located at the lower slit in Fig.1. With the help of the propagation equation (12), we

7 determine its time evolution

ρ = Gb[v] ρ0, (26)

where we use the the wave function (21) in the velocity field (24) to describe the Madelung fluid. We depicted the resulting probability density ρ in Fig.2(c) for the same parameters as we used in Figs.2(a) and2(b). With increasing time, the probability density ρ still covers only values in the lower part behind the double-slit, but interference can be observed. Similar patterns were created by Zheng et al.22 with the help of variant simulations. In Fig.2(c) it looks as if the upper maximum is not connected to the others, but there is just a low probability to find the particles. In a de Broglie-Bohm interpretation of quantum mechanics this would mean that the particles move fast at these positions.

To reproduce results of physical double-slit experiments, the initial probability density ρ0 must match the squared modulus |ψ(y, 0)|2 of the quantum mechanical wave function which is used in the velocity field (24). This is not fulfilled by the initial probability density (25) of our thought experiment. Nevertheless, this thought experiment is very interesting from a mathematical point of view. It shows that our formalism combines characteristics of quantum mechanics with char- acteristics of the de Broglie-Bohm theory. This can be seen from Fig.2, where the result of our formalism (c) describes continuous probability densities as in quantum mechanics (a), but also offers the possibility to propagate parts of the initial probability density as in the de Broglie-Bohm theory (b).

V. SUMMARY AND CONCLUSION

In summary, we have rewritten the Madelung continuity equation as a propagation equa- tion which can be used to propagate particle probability densities in a Madelung fluid. We have shown that the derived propagation equation goes over into the guidance equation of the de Broglie-Bohm theory, where only single particles are propagated. As an example, we illustrated similarities of our formalism with quantum mechanics and with the de Broglie- Bohm theory by propagating particle probability densities in a double-slit experiment. A possible application of our formalism is to simplify the calculation of multiple trajecto- ries in the de Broglie-Bohm theory. Here, instead of calculating each trajectory individually,

8 the derived propagation equation can be used to propagate a distribution of multiple tra- jectories at once. Since our approach extends the use of the Madelung equations, it can also be a helpful tool for its fields of application, e.g. the description of Bose-Einstein condensates, the physics of condensed matter and quantum cosmology. In addition, our recursive approach to the continuity equation can be helpful not only for probabilities, but also for systems that describe mass, energy, electric charge or heat flows.

VI. ACKNOWLEDGMENT

I thank Prof. J. Zheng for initiating my interest in formulations of quantum mechanics and for many stimulating discussions on this topic. Furthermore, I thank Prof. W. P. Schleich, S. V. Aksenov, J. Fischbach, C. K. Tempel and R. Heese for critical and very fruitful discussions and O. Heim for revising the manuscript. Financial support by the National Natural Science Foundation of China NSFC (No. 61362014) and the Overseas Higher-level Scholar Project of Yunnan Province, China, is gratefully acknowledged.

REFERENCES

1E. Madelung, “Eine anschauliche Deutung der Gleichung von Schr¨odinger,” Naturwis- senschaften 14, 1004 (1926). 2A. Korn, “Schr¨odingersWellenmechanik und meine mechanischen Theorien,” Z. Phys. 44, 745 (1927). 3R. E. Wyatt, Quantum Dynamics with Trajectories: Introduction to Quantum Hydrody- namics, Interdisciplinary Applied Mathematics, Vol. 28 (Springer, New York, 2005). 4L. Pitaevskii and S. Stringari, Bose-Einstein Condensation, International Series of Mono- graphs on Physics, Vol. 116 (Clarendon Press, Gloucestershire, 2003). 5M. A. Khamehchi, K. Hossain, M. E. Mossman, Y. Zhang, T. Busch, M. M. Forbes, and P. Engels, “Negative-mass hydrodynamics in a -orbit–coupled Bose-Einstein conden- sate,” Phys. Rev. Lett. 118, 155301 (2017). 6R. Fedele, D. Anderson, and M. Lisak, “How the coherent instabilities of an intense high- energy charged-particle beam in the presence of nonlocal effects can be explained within

9 the context of the Madelung fluid description,” Eur. Phys. J. B 49, 275 (2006). 7L. Wei and Y.-N. Wang, “Quantum ion-acoustic waves in single-walled carbon nanotubes studied with a quantum hydrodynamic model,” Phys. Rev. B 75, 193407 (2007). 8S¸. A. B¨uy¨uka¸sıkand O. K. Pashaev, “Madelung representation of damped parametric quantum oscillator and exactly solvable Schr¨odinger-Burgersequations,” J. Math. Phys. 51, 122108 (2010). 9I. Carusotto and R. Balbinot, “: Acoustic Hawking radiation,” Nat. Phys. 12, 897 (2016). 10G. Baumann and T. F. Nonnenmacher, “Lie transformations, similarity reduction, and solutions for the nonlinear Madelung fluid equations with external potential,” J. Math. Phys. 28, 1250 (1987). 11A. Visinescu, D. Grecu, R. Fedele, and S. De Nicola, “Madelung fluid description of the generalized derivative nonlinear Schr¨odingerequation: Special solutions and their stability,” Theor. Math. Phys. 160, 1066 (2009). 12C. Rogers and B. Malomed, “On Madelung systems in nonlinear optics: A reciprocal invariance,” J. Math. Phys. 59, 051506 (2018). 13J. K¨oppe, M. Patzold, W. Grecksch, and W. Paul, “Quantum Hamilton equations of motion for bound states of one-dimensional quantum systems,” J. Math. Phys 59, 062102 (2018). 14L. D. Landau and E. M. Lifshitz, Hydrodynamik, Lehrbuch der theoretischen Physik, Vol. 6 (Akademie Verlag, Berlin, 1991). 15D. Bohm, “A suggested interpretation of the quantum theory in terms of ”hidden” vari- ables. I,” Phys. Rev. 85, 166 (1952). 16A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinski, Methods of in Statistical Physics, Selected Russian Publications in the Mathematical Sciences, Vol. 2 (Englewood Cliffs, NJ: Prentice-Hall, 1964). 17Electrons et Photons: Rapports et Discussions du Cinquieme Conseil de Physique tenu a Bruxelles du 24 au 29 October 1927 sous les auspices de l’Institut International Physique Solvay (Gauthier-Villars, Paris, 1928). 18D. Bohm, “A suggested interpretation of the quantum theory in terms of ”hidden” vari- ables. II,” Phys. Rev. 85, 180–193 (1952).

10 19P. R. Holland, The quantum theory of motion: An account of the de Broglie-Bohm causal interpretation of quantum mechanics (Cambridge University Press, 1995). 20C. J¨onsson,“Elektroneninterferenzen an mehreren k¨unstlich hergestellten Feinspalten,”Z. Phys. 161, 454 (1961). 21D. J. Bohm and B. J. Hiley, “On the intuitive understanding of nonlocality as implied by quantum theory,” Found. Phys. 5, 93–109 (1975). 22J. Zheng and C. Zheng, “Variant simulation system using quaternion structures,” J. Mod. Opt. 59, 484 (2012).

11