Medical, Pharma, Engineering, Science, Technology and Business

^{1}Wigner Research Center of the Hungarian Academy of Sciences Konkoly-Thege út 29-33, 1121 Budapest, Hungary

^{2}ELI-HU Nonprofit Kft, Dugonics Tér 13, H-6720 Szeged, Hungary

^{3}University of Pécs, Institute of Physics, Ifjúság útja 6 H-7624 Pécs, Hungary

^{4}Sapientia University, Department of Bioengineering, Libertatii sq. 1, 530104 Miercurea Ciuc, Romania

- *Corresponding Author:
- Barna IF

Wigner Research Center of the

Hungarian Academy of Sciences Konkoly

Thege út 29-33, 1121 Budapest, Hungary+36-1-392-2222/3504

Tel:+36-1-392-2598

Fax:

**E-mail:**[email protected]

**Received Date:** June 16, 2017; **Accepted Date:** July 05, 2017; **Published Date:** July 10, 2017

**Citation: **Barna IF, Pocsai MA, Mátyás L (2017) Analytic Solutions of the Madelung
Equation. J Generalized Lie Theory Appl 11: 271. doi: 10.4172/1736-4337.1000271

**Copyright:** © 2017 Barna IF, et al. This is an open-access article distributed under
the terms of the Creative Commons Attribution License, which permits unrestricted
use, distribution, and reproduction in any medium, provided the original author and
source are credited.

**Visit for more related articles at** Journal of Generalized Lie Theory and Applications

We present analytic self-similar solutions for the one, two and three dimensional Madelung hydrodynamical equation for a free particle. There is a direct connection between the zeros of the Madelung fluid density and the magnitude of the quantum potential.

Quantum mechanics; Partial differential equations; Madelung equation; Planck’s constant

Finding classical physical basements of quantum mechanics is a great challenge since the advent of the theory. Madelung was one among the firsts who gave one explanation, this was the hydrodynamical foundation of the Schrödinger equation [1,2]. His exponential transformation simply indicates that one can model quantum statistics hydrodynamically. Later, it became clear that the Madelung Anzatz is just the complex Cole-Hopf transformation [3,4] which is sometimes used to linearize non-linear partial differential equations(PDEs).

The transformed equation has an attractive feature that the Planck’s constant appears only once, as the coefficient of the quantum potential or pressure. Thus, the fluid dynamicist can gather experience of its effects by translating some of the elementary situations of the quantum theory into their corresponding fluid mechanical statements and vice versa.

The quantum potential also appears in the de Broglie-Bohm pilot wave theory [5,6] (in other context) which is a non-mainstream attempt to interpret quantum mechanics as a deterministic non-local theory. In the case of the Euler equation goes over the Hamilton-Jacobi equation.

As an interesting peculiarity Wallstrom showed with mathematical means that the initial-value problem of the Madelung equation is not well-defined and additional conditions are needed [7].

Nowadays, hydrodynamical description of quantum mechanical systems is a popular technical tool in numerical simulations. Review articles on quantum trajectories can be found in a booklet edited by Huges in 2011 [8].

From general concepts as the second law of thermodynamics a weakly non-local extension of ideal fluid dynamics can be derived which leads to the Schrödinger-Madelung equation as well [9].

In our following study we investigate the Madelung equation with the self-similar Ansatz and present analytic solutions with discussion.

This way of investigation is a powerful method to study the global properties of the solutions of various non-linear PDEs [10]. Self-similar Ansatz describes the intermediate asymptotic of a problem: it is hold when the precise initial conditions are no longer important, but before the system has reached its final steady state. This is much much simpler than the full solutions and so easier to understand and study in different regions of parameter space. A final reason for studying them is that they are solutions of a system of ordinary differential equations(ODEs) and hence do not suffer the extra inherent numerical problems of the full PDEs. In some cases self-similar solutions helps to understand diffusion-like properties or the existence of compact supports of the solution.

In the last years we successfully applied the multi-dimensional generalization of the self-similar Ansatz to numerous viscous fluid equations [11,12] ending up with a book chapter in ref. [13].

To our knowledge there are no direct analytic solutions available for the Madelung equation. Baumann and Nonnenmacher [14] exhaustively investigated the Madelung equation with Lie transformations and presented numerous ODEs, however non exact and explicit solutions are presented in a transparent way. Additional numerous studies exist where the non-linear Schrödinger equation is investigated with the Madelung Ansatz ending up with solitary wave solutions, [15] however that is not the field of our present interest.

Following the original paper of Madelung [2] the time-dependent Schrödinger equation reads:

(1)

where Ψ,U,m,h are the wave function, potential, mass and Planck’s constant, respectively. Taking the following Ansatz where *ρ*(*x*,*t*) and *S*(*x*,*t*) are time and space dependent functions. Substituting this trial function into eqn. (1) going through the derivations the real and the complex part gives us the following continuity and Euler equations with the form of:

(2)

where the following substitution has to be made . The *ρ* is the density of the investigated fluid and v is the velocity field. Madelung also showed that this is a rotation-free flow. The transformed equations has an attractive feature that the Planck’s constant appears only once, at the coefficient of the quantum potential or pressure, which is the first term of the right hand side of the second equation. Note, that these are most general vectorial equation for the velocity field *v* which means that one, two or three dimensional motions can be investigated as well. In the following we will consider the two dimensional flow motion *v*=(*u*,*v*) in Cartesian coordinates without any external field *U*=0. The functional form of the three and one dimensional solutions will be mentioned briefly as well.

We are looking for the solution of eqn. (2) with self-similar Ansatz which is well-known from ref. [10]

(3)

where *f*, *g* and *h* are the shape functions of the density and the velocity field, respectively. The similarity exponents *α*,*β*,*δ*,*ε* are of primary physical importance since *α*,*δ*,*ε* represents the damping of the magnitude of the shape function while *β* represents the spreading. More about the general properties of the Ansatz can be found in our former papers [11,12]. Except some pathological cases all positive similarity exponents mean physically relevant dispersive solutions with decaying features at *x*,*y*,*t*→∞. Substituting the Ansatz (3) into (2) and going through some algebra manipulation the next ODE system can be expressed for the shape functions.

(4)

The first continuity equation can be integrated giving us the mass as a conserved quantity and the parallel solution for the velocity fields *η*=2(*g*+*h*)+ *c*_{0} where *c*_{0} is the usual integration constant, which we set to zero. (A non-zero *c*_{0} remains an additive constant in the final ODE (5) as well.) It is interesting, and unusual (in our practice) that even the Euler equation can be integrated once giving us the conservation of momenta. For classical fluids this is not the case. After some additional algebraic steps a decoupled ODE can be derived for the shape function of the density.

(5)

All the similarity exponents have the fixed value of which is usual for regular heat conduction, diffusion or for Navier-Stokes equations [13]. Note, that the two remaining free parameters are the mass of the particle *m* and which is the Planck’s constant divided by 2π. For a better transparency we fix . This is consistent with experience of regular quantum mechanics that quantum features are relevant at small particle masses.

At this point it is worth to mention, that the obtained ODE for the density shape function is very similar to eqn. (5) for different space dimensions, the only difference is a constant in the last term. For one, two or three dimensions the denominator has a factor of 1,2 or 3, respectively.

An additional space dependent potential *U* (like a dipole, or harmonic oscillator interaction) in the original Schrödinger equation would generate an extra fifth term in eqn. (5) like, *f*(*η*),*η*^{2}. Unfortunately, no other analytic closed form solutions can be found for such terms (**Figures 1** and **2**).

The solution of (5) can be expressed with the help of the Bessel functions of the first and second kind [16] and has the following form of:

(6)

where *c*_{1} and *c*_{2} are the usual integration constants. The correctness of this solutions can be easily verified via back substitution into the original ODE.

To imagine the complexity of these solutions **Figures 1** and **2** present *f*(*η*) for various *m*, *c*_{1} and *c*_{2} values. It has a strong decay with a stronger and stronger oscillation at large arguments. The function is positive for all values of the argument, (which is physical for a fluid density), but such oscillatory profiles are completely unknown in regular fluid mechanics [13]. The most interesting feature is the infinite number of zero values which cannot be interpreted physically for a classical real fluid.

The presented form of the shape function cannot be simplified further, only *Y _{v}*s can be expresses with the help of

(7)

Both *J _{v}* and

(8)

Now,

(9)

**Figure 3** shows the projection of the real part wave function to the *x*,*t* sub-space. At small times the oscillations are clear to see, however at larger times the strong damping is evident.

For arbitrary quantum systems, the wave function can be evaluated according to the Schrödinger equation, however we never know directly how large is the quantum contribution to the classical one. Now, it is possible for a free particle to get this contribution. (The Schrödinger equation gives the Gaussian wave function for a freely propagating particle). With the Madelung Ansatz we got the classical fluid dynamical analogue of the motion with the physical parameters *ρ*(*x*,*y*,*t*),*v*(*x*,*y*.*t*) which can be calculated analytically via the self-similar Ansatz thereafter original wave function Ψ(*x*,*y*,*t*) of the quantum problem can be evaluated as well. The magnitude of the quantum potential *Q* directly informs us where quantum effects are relevant. This can be evaluated from the classical density of the Madelung eqn. (2) via:

(10)

**Figure 4** shows the shape function of quantum potential *Q*() comparing to the shape function of the density *f*(). Note, that where the density has zeros the quantum potential is singular. Such singular potentials might appear in quantum mechanics, however the corresponding wave function should compensate the effect. This question is analyzed in the book of Holland [17] for various other quantum systems.

After reviewing the historical development and interpretation of the Madelung equation we introduced the self-similar Ansatz which is a not-so-well-known but powerful tool to investigate nonlinear PDEs. The free particle Madelung equation was investigated in two dimensions with this method, (the one and three dimensional solutions were mentioned as well.) We found analytic solution for the fluid density, velocity field and the original wave function. All can be eliminated with the help of the Bessel functions. The classical fluid density has interesting properties, oscillates and has infinite number of zeros which is quite unusual and has not yet been seen in such analysis.

- Madelung E (1926) Eine anschauliche Deutung der Gleichung von Schrödinger. Naturwissenschaften 14: 1004.
- Madelung E (1927) Zur Frage der Elektronenaffinität des Wasserstoffs. Z Phys 40: 331-332.
- Hopf E (1950) Commun Pure Appl Math 3: 201.
- Cole JD (1951) On a quasi-linear parabolic equation occurring in aerodynamics. Quart Appl Math 9: 225-336.
- de Broglie L (1927) L'univers à cinq dimensions et la mécanique ondulatoire. Journ de Phys et la Rad 8: 65-73.
- Bohm D (1952) Phys Rev 15: 166.
- Wallstrom TC (1994) Phys Lett A 184: 229.
- Hughes KH, Parlant G (2011) Quantum Trajectories (CCP6, Daresbury Laboratory, 2011).
- Ván P, Fülöp T (2006) Weakly non-local fluid mechanics: the Schrödinger equation. Proc R Soc A, pp: 462-541.
- Sedov LI (1993) Similarity and Dimensional Methods in Mechanics CRC Press.
- Barna IF, Mátyás L (2014) Analytic solutions for the three-dimensional compressible Navier–Stokes equation. Fluid Dyn Res 46: 055508.
- Barna IF, Mátyás L (2015) Chaos Solitons and Fractals 78: 249.
- Campos D (2016) Handbook on Navier-Stokes Equations, Theory and Applied Analysis, Nova Publishers, New York.
- Baumann G, Nonnenmacher TF (1987) Lie transformations, similarity reduction, and solutions for the nonlinear Madelung fluid equations with external potential. Journ Math Phys 28: 1250.
- Grecu D, Alexandru T, Grecu AT, De Nicola S (2008) Solitary Waves in a Madelung Fluid Description of Derivative NLS Equations. Journal of Nonlinear Mathematical Physics 15: 209-219.
- Olver FWJ, Lozier DW, Boisvert RF, Clark CW (2010 NIST Handbook of Mathematical Functions Hardback and CD-ROM. Cambridge University Press.
- Holland PR (1995) The quantum theory of motion. Cambridge University Press.

Select your language of interest to view the total content in your interested language

- Adomian Decomposition Method
- Algebra
- Algebraic Geometry
- Analytical Geometry
- Applied Mathematics
- Axioms
- Balance Law
- Behaviometrics
- Big Data Analytics
- Binary and Non-normal Continuous Data
- Binomial Regression
- Biometrics
- Biostatistics methods
- Clinical Trail
- Combinatorics
- Complex Analysis
- Computational Model
- Convection Diffusion Equations
- Cross-Covariance and Cross-Correlation
- Deformations Theory
- Differential Equations
- Differential Transform Method
- Fourier Analysis
- Fuzzy Boundary Value
- Fuzzy Environments
- Fuzzy Quasi-Metric Space
- Genetic Linkage
- Geometry
- Hamilton Mechanics
- Harmonic Analysis
- Homological Algebra
- Homotopical Algebra
- Hypothesis Testing
- Integrated Analysis
- Integration
- Large-scale Survey Data
- Latin Squares
- Lie Algebra
- Lie Superalgebra
- Lie Theory
- Lie Triple Systems
- Loop Algebra
- Matrix
- Microarray Studies
- Mixed Initial-boundary Value
- Molecular Modelling
- Multivariate-Normal Model
- Noether's theorem
- Non rigid Image Registration
- Nonlinear Differential Equations
- Number Theory
- Numerical Solutions
- Operad Theory
- Physical Mathematics
- Quantum Group
- Quantum Mechanics
- Quantum electrodynamics
- Quasi-Group
- Quasilinear Hyperbolic Systems
- Regressions
- Relativity
- Representation theory
- Riemannian Geometry
- Robust Method
- Semi Analytical-Solution
- Sensitivity Analysis
- Smooth Complexities
- Soft biometrics
- Spatial Gaussian Markov Random Fields
- Statistical Methods
- Super Algebras
- Symmetric Spaces
- Theoretical Physics
- Theory of Mathematical Modeling
- Three Dimensional Steady State
- Topologies
- Topology
- mirror symmetry
- vector bundle

- 7th International Conference on Biostatistics and Bioinformatics

September 26-27, 2018 Chicago, USA - Conference on Biostatistics and Informatics

December 05-06-2018 Dubai, UAE - Mathematics Congress - From Applied to Derivatives

December 5-6, 2018 Dubai, UAE

- Total views:
**660** - [From(publication date):

August-2017 - Jul 17, 2018] - Breakdown by view type
- HTML page views :
**606** - PDF downloads :
**54**

Peer Reviewed Journals

International Conferences 2018-19