# Self-Similarity Analysis of the Nonlinear Schrödinger Equation in the Madelung Form

Authors: Imre F. Barna; Mihály A. Pocsai; L. Mátyás
Journal: Advances in Mathematical Physics
 (2018)

---

## Abstract

In the present study a particular case of Gross-Pitaevskii or nonlinear Schrödinger equation is rewritten to a form similar to a hydrodynamic Euler equation using the Madelung transformation. The obtained system of differential equations is highly nonlinear. Regarding the solutions, a larger coefficient of the nonlinear term yields stronger deviation of the solution from the linear case.


---

## Body

## 1. Introduction

In 1926, at the advent of the quantum mechanics Madelung [1, 2] gave a hydrodynamical foundation of the Schrödinger equation and raised interpretation questions which are open till today. Later, people realized that the Madelung Anzatz is just the complex Cole-Hopf transformation [3, 4] which is sometimes applied to linearize nonlinear partial-differential equations (PDEs). The Madelung form of the Schrödinger equation has a remarkable property that Planck’s constant (representing quantum mechanics) appears only once, namely, as the coefficient of the quantum potential. This quantity also appears in the nonmainstream attempt of quantum mechanics (the de Broglie-Bohm pilot wave theory [5] (in other context)) which interprets quantum mechanics as a deterministic nonlocal theory. The Madelung transformation may create a common language between scientist of fluid dynamics and quantum mechanics and can be fruitful to give an additional interpretation of quantum mechanics. In the last decades, hydrodynamical description of quantum mechanical systems become a popular technical tool in numerical simulations. Detailed reviews in this field can be found in a booklet [6] or in Wyatt’s monograph [7].The application of the original Madelung (or to say Cole-Hopf) transformation grew out the linear quantum theory and were applied to more sophisticated field theoretical problems. Recently, Chavanis and Matos [8] applied the Madelung Ansatz to the Klein-Gordon-Maxwell-Einstein equation in curved space-time possibly describing self-gravitating Bose-Einstein condensates, coupled to an electromagnetic field. This study clearly presents the wide applicability of the original Madelung Ansatz to refined quantum field theories with the possible hydrodynamic interpretation picture.As an analytical tool of investigation we use the self-similar Ansatz of [9]. It is a powerful method to study the global properties of the dispersive solutions of various nonlinear partial-differential equations (PDEs) In other words, the intermediate asymptotic of the problem is described. Therefore this solution is valid when the precise initial conditions are no longer important; however, the final steady state is not yet reached. For real physical systems this region attracts the most attention. In contrast to the full solution it is much simpler and easier to understand for the entire parameter space. The second reason for studying them is that they are solutions of a system of ordinary differential equations (ODEs) (not a PDE) and thus do not suffer the full inherent numerical problems of the full PDEs. Most of the cases this tool helps to understand some of the global properties of the solution, e.g., shock waves or the existence of compact supports. First, we applied this Ansatz to the modified Cattaneo-Venotte heat-conduction problem [10] and obtained a compact solution with finite propagation speed. Later, we successfully generalized the Ansatz to multispace dimensions and studied numerous viscous fluid equations [11, 12] ending up with a book chapter of [13]. As main results we are now able to interpret the birth of Benard cells in heated fluids [14].In our former study [15] we investigated the Schrödinger equation in the Madelung form with the self-similar Ansatz and presented analytic solutions with discussion. Unlike to classical fluids the obtained density of the Madelung fluid has infinite number of zeroes which is a fingerprint of the quantum mechanical origin. Now we analyze the Madelung form of nonlinear Schrödinger (NLS) equation with the same Ansatz and compare the obtained results to the previous linear one. As far as we know, there are no direct clear-cut analytic solutions existing for the Madelung form of the NLS equation by now. It is important to mention that celebrated standard numerical algorithms exist for solving the NLS [16, 17]. We hope that our present study helps to attract attention to this question. Parallel to the well-known quantum mechanical picture there is a looming hydrodynamical world hidden behind the Schrödinger equation from the very beginning which is poorly understood.


## 2. Theory and Results

In a weakly interacting Bose system, the macroscopic wave functionΨ appears as the order parameter of the Bose-Einstein condensation, obeying the Gross-Pitaevskii (GP) equation [18–21]:(1)iħ∂Ψ∂t=-ħ2∇22m+nΨ2-μΨ.Here, n=4πħ2a/m represents the coupling constant characterized by the s-wave scattering length a. The particle mass and the chemical potential are m and μ, respectively. Due to the second (nonlinear) term of the right hand side of this equation, it is definitely different from the linear Schrödinger equation [15]; therefore the following analysis with the presented results are independent as well.The role and the properties of the GP system has been investigated from various viewpoints. Without completeness we mention some of them. The question of collective excitations of the trapped Bose gases was analyzed by [22]. Later, it is considered in the study of self-organization of the Bose-Einstein (BE) condensates [23].Certain generalized versions of (1) have been also applied to the study of BE condensate of polaritons [24–32].Considering the Madelung Ansatz for the order parameter with the amplitude and phaseΨ=ρeiϕ where the density of the condensate is ρ and the superfluid velocity is v=(ħ/m)∇ϕ. (The expression “superfluid velocity” is taken from Tsubota et al. [33] and means, in this content, that the fluid is inviscid.)The obtained hydrodynamical equations, the continuity and the Euler, are as follows:(2)∂ρ∂t+∇·ρv=0,∂v∂t+v·∇v=ħ22m2∇▵ρρ-∇mρnρ22,where we consider free motion (no extra potential energy term U(r) is used) and the chemical potential μ was set to zero as well. It is easy to show that an additional chemical potential term μ gives no contribution to the final hydrodynamical equations. The role of the external potential energy term U(r) will be mentioned later.The first term of the right hand side of the second equation is the quantum pressure term. The second term makes the difference to the original Madelung equations; this is a kind ofsecond quantum pressure due to the nonlinear origin of the GP equation. Not, 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. The functional form of the three- and one-dimensional equation system will be mentioned briefly, too.We search the solution of (2) in the form of the self-similar Ansatz [9]: (3)ρx,y,t=t-αfx+ytβ≔t-αfη,ux,y,t=t-δgη,vx,y,t=t-ϵhη,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 α,δ,ϵ represent the damping of the magnitude of the shape function. The exponent β represents the spreading of the solution. Additional technical details about the general features of the Ansatz can be found in our former papers [12, 14]. Except some extraordinary cases only all positive similarity exponents correspond to physically relevant dispersive solutions which mean decaying features at x,y,t→∞. Usually negative exponents describe exploding solutions. Substituting the Ansatz (3) into (2) and executing additional algebraic manipulation the following ODE system can written down for the shape functions:(4)-αf-α2f′η+f′g+fg′+f′h+fh′=0,-32α-1g-α2g′η+gg′+hg′=ħ22m2f′32f3-f′f″f2+f″′2f-nmf′,-32α-1h-α2h′η+gh′+hh′=ħ22m2f′32f3-f′f″f2+f″′2f-nmf′,where prime means derivation with respect to the new variable η. The relations among the similarity exponents are the following: α can be taken arbitrary, β=α/2 and ϵ=δ=3β-1=3α/2-1. At this point we have to emphasize, that a particular form of the external potential U(x,y,t) can be added to (1) without leading to contradiction among the similarity exponents. The form in 2D is the following: U(x,y,t)=1/t2β·v(η) where v(η) can be any kind of function with existing first derivative.(Of course, there is no evidence that the obtained modified (4) could be solved with analytical means.)Note, that the ODE system of the linear Schrödinger equation [15] has a much simpler form of(5)-12f-12f′η+f′g+fg′+f′h+fh′=0,-12g-12g′η+gg′+hg′=ħ22m2f′32f3-f′f″f2+f″′2f,-12h-12h′η+gh′+hh′=ħ22m2f′32f3-f′f″f2+f″′2f,where all the similarity exponents have the fixed value of +1/2. Our experience tells that this is usual for regular heat conduction, diffusion or for noncompressible Navier-Stokes equations [13]. For the linear case after some additional trivial algebraic steps a decoupled ODE can be derived for the shape function of the density: (6)2f″f-f′2+m2η2f22ħ2=0.For one, two, or three dimensions the multiplicative factor ofħ in the last term is different. It is 1, 2, or 3, respectively.With the help of Bessel functions [34] the solution of (6) can be expressed:(7)fη=πη64c1J1/42mη28-c2Y1/42mη282,where c1 and c2 are the usual integration constants. All the additional properties of this solution can be found in our former paper [15]. (There is a misprint in the final [15] at Eq. 7, the parenthesis is missing in front of the exponent 2. This is the correct form.)The most general equation of (4) cannot be fully integrated by analytical means. The first ODE of such systems [13] explicitly means mass conservation. Furthermore in the most cases it can be integrated analytically. However, the first equation of (4) is not a total derivative for any α. This is highly unusual among our investigated systems till now. As a second case let us consider (3/2α-1)=α/2 this gives α=1, β=δ=ϵ=1/2. Now the sum of the second and third equation of (5) can be integrated leading to (8)g22+gh-η2+h22+hη2+2nfm-ħ22m2-f′24f2+f″2f+c1=0.For the sake of simplicity we set c1 to zero. Note, that the variables g,f and h are still coupled in this single equation. However, (8) is quadratic in g(η) and the solutions can be expressed with the well-know formula of (9)g1,2=η2-h±η24-4nfm+ħ22m2-f′24f2+f″2f.Even this includes both functions of h and f. To avoid this problem let us fix the discriminant to zero, which means that the velocity shape function will be single valued. The remaining ODE reads as follows: (10)2ff″-f′2+f2η2m22ħ2-8nmfħ=0with the additional restraint that 4f2≠0 meaning that the physical density of a fluid must be positive. Note, that without the last term we get back (6) with the solution of (7). We tried numerous variable transformations unfortunately, we could not find any closed solution for the full equation (10). (Note that if the coefficient of the second term would be twice as the first term there would be a possible solution for the full equation including the third order term.) It is also clear from (9) that the velocity shape functions have the trivial form of g+h=η/2=x+y/2t which is the same as for the linear case.Figure1 shows the analytical solution of (7) and numerically solution of Eq. (10) for various nonlinear parameters n with the same initial conditions of f(0)=1 and f′(0)=0. The ħ=m=1 units are used. All solutions have a strong damping with stronger and stronger oscillations at large arguments. The functions are 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. Remembering the restraint 4f2≠0 of the calculation the solution function falls into infinite number of finite intervals. We think that this is a clear fingerprint that the obtained Euler equations have a quantum mechanical origin.Figure 1
Various shape functions of the density. The blue curve is the analytic solution of equation (7) with n=0 nonlinearity parameter. Black, red, and green functions are numerical results of equation (10) for n=0.078,0.11 and 0.13 nonlinearity parameters. (All curves are for the same initial conditions of f(0)=1 and f′(0)=0.) Note that larger parameters mean larger deviation from the analytic result.The major difference to the analytic solution is the following the∫0∞fρdρ can be evaluated (now of course just numerically up to an arbitrary ρmax) giving us a convergent integral. This clearly shows that the wave function of the corresponding quantum mechanical system (see below) is the element of the L2 space. We investigated numerous initial conditions all giving the same result of a finite integral. Of course, this is just an empirical statement and not a rigorous mathematical proof.Theρ function could be interpreted as the fluid mechanical analogue of the real part of the wave function of the free GP quantum mechanical particle. For the linear case the square root of (7) is a kind of far analogue of the original free Gaussian wave packet. Unfortunately, we cannot find any direct transformation between these two functions.To obtain the complete original wave function, the imaginary part has to be evaluated as well. It is simple fromη=x+y/t1/2=2(g+h) that (11)S=mħ∫r0r1vdr=mħx+y24t.Now(12)Ψx,y,t=ρx,y,teiSx,y,t=t-1·fx+yt1/2eim/ħx+y2/4t,where f(η) is now a numerical function. Figure 2 shows the projection of the real part of the wave function to the x,t subspace. At small times there are clear oscillations, but at larger times the strong damping is evident. If the time is fixed (t=t0 so to say) there are infinite number of spatial oscillations which can be interpreted as a quantum mechanical heritage of the system.Figure 2
The projection of the real part of the wave function (12) along the x-t plane.


## 3. Summary and Outlook

We gave a short historical review of the interpretation and development of the Madelung equation. As second point we introduced (the not so well known) self-similar Ansatz which is a robust tool to investigate the local and global properties of nonlinear PDEs. Finally, the Madelung form of the free NLS equation was investigated and the results were compared to the ones obtained from the linear Schrödinger equation. We think that the highly oscillatory fluid density (with the infinite number of zeroes) are a clear fingerprint of the quantum mechanical origin of the system. We also believe that the noticeable quadratic argument of the solutions might get an interpretation in the future. However, the presented approach is still evolving and works are in progress. We plan to apply our method to the Madelung equations of more general nonlinear Schrödinger equations like [35].



---

*Source: 7087295-2018-10-01.xml*