Advances in Mathematical Physics

Volume 2018, Article ID 7087295, 5 pages

https://doi.org/10.1155/2018/7087295

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

^{1}Wigner Research Center for Physics of the Hungarian Academy of Sciences, Konkoly-Thege Miklós ú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, Libertătii Sq. 1, 530104 Miercurea Ciuc, Romania

Correspondence should be addressed to Imre F. Barna; uh.atm.rengiw@ermi.anrab

Received 16 July 2018; Accepted 12 September 2018; Published 1 October 2018

Academic Editor: Emmanuel Lorin

Copyright © 2018 Imre F. Barna et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### 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.

#### 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]:Here, represents the coupling constant characterized by the s-wave scattering length a. The particle mass and the chemical potential are 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 where the density of the condensate is and the superfluid velocity is . (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: 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 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 of* second quantum pressure* due to the nonlinear origin of the GP equation. Not, that these are most general vectorial equation for the velocity field 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 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]: where , and 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 . 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:where prime means derivation with respect to the new variable . The relations among the similarity exponents are the following: can be taken arbitrary, and . At this point we have to emphasize, that a particular form of the external potential can be added to (1) without leading to contradiction among the similarity exponents. The form in 2D is the following: where 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 ofwhere all the similarity exponents have the fixed value of . 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:

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:where and 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 this gives , . Now the sum of the second and third equation of (5) can be integrated leading to For the sake of simplicity we set to zero. Note, that the variables and are still coupled in this single equation. However, (8) is quadratic in and the solutions can be expressed with the well-know formula of Even this includes both functions of and . 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: with the additional restraint that 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 which is the same as for the linear case.

Figure 1 shows the analytical solution of (7) and numerically solution of Eq. (10) for various nonlinear parameters with the same initial conditions of and . The 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 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.