A part of the results have been published in:

1. Suppression of the quantum-mechanical collapse by repulsive interactions in bosonic gasesHidetsugu SakaguchiDepartment of Applied Science for Electronics and Materials, Interdisciplinary Graduate School of Engineering Sciences, Kyushu University, Fukuoka, JapanBoris A. MalomedDepartment of Physical Electronics, Faculty of Engineering, Tel Aviv University, Tel Aviv, Israel

2. A part of the results have been published in:

3. (1) Introduction and the modelThe attractive potential ~1/r2plays the well-known critical role in the standard quantum mechanics, leading to the quantum collapse, alias “fall onto the center” in 3D and 2D geometries. The scaled form of the Schrödinger equation with this potential (the harmonictrapping potential, Ω²r²/2,is included too) is

4. The manifestation of the criticality: The 3D Schrödinger equation with potential -1/r2 has no ground state for U0 > (U0)cr ≡ ¼, and does have it at U0 < ¼. The 2D Schrödinger equation has no ground state for anyU0 > 0. This phenomenon of the quantum collapse is also called “fall onto the center” in the book Quantum Mechanics by Landau and Lifshitz.

5. The physical interpretation of both the 3D and 2DSchrödinger equations with this potential: the interaction between a particle carrying a permanent electricdipole momentd, and an electric chargeQ placed at the origin. If the orientation of the moment is locked to the local electric field (E), the effective interaction potential is precisely -(1/2)U0/r2, withU0 = 2|Q|d :

6. In particular, if Q = 1 is the electron’s charge, and the mass of the particle is taken as ~ 10 proton masses, the critical value of the attraction strength U0 = ¼, corresponds, in physical units, to a verysmall value of the dipole moment, d ~ 10-5 Debye, hence the case of U0 > ¼ is essential indeed.

7. In the 2D geometry (which corresponds to a pancake-shaped BEC configuration), another interpretation is possible too, in addition to the above-mentioned one: the interaction of a particle carrying a magnetic dipole moment with a current filament (e.g., an electron beam) which transversely intersects the 2D layer, assuming that the orientation of the magnetic moment is locked to the local magnetic field induced by the current.

8. A solution to the quantum anomaly problem in this quantum-mechanical setting (i.e., the nonexistence of the ground state) was proposed outside of the framework of quantum mechanics – in terms of the linear quantum-field theory. Essentially, that solution postulates that the ground state is created by the field-theory renormalization procedure. The so created ground state features an arbitrary spatial scaleimposed by the renormalization: K. S. Gupta and S. G. Rajeev, Phys. Rev. D 48, 5940 (1993); H. E. Camblong, L. N. Epele, H. Fanchiotti, and C. A. Garcia Canal, Phys. Rev. Lett. 85, 1590 (2000).

9. Our objective is to propose a different solution of the quantum-collapse problem: taking into account, in the mean-field approximation, collisions between dipolar particles attracted to the center by the potential -1/r2, i.e., replacing the Schrödinger equation by the corresponding Gross-Pitaevskii equation. This approach will create a ground state, with the spatial scaleuniquely determined by physical parameters of the system.

10. The corresponding Gross-Pitaevskii equation, with the contact repulsive nonlinearity: The realization of ultracold quantum gases formed by molecules carrying electric dipolemoments may be possible, using (for instance) Li-Cs bosonic molecules: J. Deiglmayr, A. Grochola, M. Repp, K. Mörtlbauer, C. Glück, J. Lange, O. Dulieu, R. Wester, and M. Weidemüller, Phys.Rev. Lett. 101, 133004 (2008).

11. Taking into account the dipole-dipole interactions in the Gross-Pitaevskii equation:

12. For m ~ 10 proton masses and as ~ 10 nm, the renormalization of the scattering length is essential for d ~ 1 Debye.

13. (2) Presentation of results (i) The standard description of the quantum collapse will be recapitulated. (ii)The central point: In the 3D case, it is demonstrated that the repulsive cubic nonlinearity suppresses the quantum collapse, i.e., it creates the ground state which is missing in the linear Schrödinger equation. (iii) A quantum phase transition is found at U0 = 1 ≡ 4 (U0)cr . (iv) If the potential -(1/2)U0/r2 is not strong enough to cause the quantum collapse in 3D, the action of the cubic nonlinearity in the combination with the harmonic trapping potential gives rise to a tri-stability of bound states.

14. (v)It will be shown that the cubic nonlinearity is insufficiently strong to suppress the quantum collapse in the 2D setting. However, the quintic nonlinearitydoes it (this results is a more formal one). (vi) It will be demonstrated (in a brief form) that the nonlinearity also suppresses the quantum collapse and creates the missing ground state in the case when a strong external uniform polarizing field reduces the symmetry of the 3D setting from spherical to cylindrical.

15. (i) Recapitulation of the description of the quantum collapse in the linear Schrödinger equation Solutions to the 3D and 2DlinearSchrödinger equations (including the harmonic trap), with angular momentumℓ, are sought for as, respectively,

17. (ii) The suppression of the quantum collapse by the repulsive cubic nonlinearity in 3D (no harmonic trap, Ω = 0).

18. Using this approximation, one can calculate the norm of the solution, N, and thus derive an approximate μ(N) characteristic for the ground state:

19. The family of the ground-state solutions was then found in a numerical form, as a solution of the radial equation for the 3D wave function. The solutions, bothanalytical and numerical, exist equally well for U0< ¼ and U0> ¼, i.e., in cases when the linear Schrödinger equationdoes and does not have the ground state. Numerical simulations confirm that all the ground states are stable.

20. Examples of the numerically foundground states of the 3D nonlinear equation without the external trap, Ω = 0(continuous lines), and its comparison to the corresponding simplest analytical prediction (dashed lines):Radial profiles μ(N) curves μ(N) curvesfor U0 = 0.8 > ¼ for U0 = 0.8 > ¼ for U0 = 0.1 < ¼(here, N = 6.26)

22. (iii) The quantum phase transition at U0 = 1 ≡ 4 (U0)cr

24. (iv) The effect of the additional harmonic trapping potential (Ω² > 0) in the 3D case The harmonic trap (Ω² > 0) distorts the branch of the ground states which was created by the nonlinearity, and additionally supports extensions of the twobound modes that were found in the exact form as solutions to the linear Schrödinger equation with Ω² > 0. These additional modes exist precisely at U0< ¼ (where the linear equationdoes not give rise to the quantum collapse). The two additional modes were shown (by dint of direct simulations) to be stable. Thus, the nonlinear equation gives rise to the tri-stability at U0< ¼.

25. The illustration of the tri-stability at Ω² = 0.1 andU0 = 0.2 < ¼ . At N→ ∞, the asymptotic form of all the three branches of μ(N)is given by the Thomas-Fermi approximation:

26. (v) The suppression of the quantum collapse in the 2D setting by the quintic repulsive term In the 2D situation, the cubic repulsive term formally gives rise to the wave function, which replaces the collapsing one, with the asymptotic form ~ 1/r at r → 0 (the same as in 3D) . However, unlike the 3D case, the norm of such a wave function diverges (logarithmically) at r → 0 in 2D, i.e., the cubic term is not strong enough to create the missing physically acceptable 2D ground state.

27. The suppression of the 2D quantum collapse is provided by a quintic term, which may account for (elastic) triple collisions in the quantum gas: F. K. Abdullaev, A. Gammal, L. Tomio, and T. Frederico, Phys. Rev. A 63, 043604 (2001).

30. Examples of the radial profiles of the 2Dground state created by the quintic nonlinear term, and respective μ(N) curves (solid and dashed curves depict numerical findings and the analytical approximation, respectively):The radial profiles The μ(N) curves for U0= +0.05 for U0= +0.05 and U0 = -0.18

31. (vi) The anisotropic extension of the 3D setting: If the local orientation of the electric dipoles is frozen into a strong external uniform field, the symmetry of the potential of the interaction of the dipoles with the central electric charge is reduced from spherical to cylindrical, U(r ,θ) = -(U0/2)r-²cosθ, where θis the angular coordinate in 3D. The cylindrically-symmetric 3Dground-state wave function is created by the nonlinearity in this case too: H. Sakaguchi and B.A. Malomed, Phys. Rev. A 84, 033616 (2011).

32. (3) CONCLUSIONS (i) The problem of the “fall onto the center/quantum collapse/quantum anomaly”, induced by the attractive potential –U0/r ² in 3D and 2D settings, can be resolved by the self-repulsive cubic (in 3D) or quintic (in 2D) terms, which are added to the respective linear Schrödinger equation, transforming it into the Gross-Pitaevskii equation. The nonlinear term creates the stable ground state, which was absent in the linear equation. The spatial scale of the newly created ground state is uniquely determined by physical parameters of the system. These 3D and 2D models may be realized as BECof molecules or atoms composed of dipolar molecules or atoms.

33. (ii) The quantum phase transition, signaled by the non-analyticity in the dependence of the wave-function structure at r → 0 onU0 , occurs atU0 ≡ 4 (U0)cr . (iii) In the case when the linear equationdoes not yet give rise to the collapse, and the harmonic trapping potential is added, the corresponding nonlinear equation gives rise to the tri-stability of the bound states in the 3D case.

34. (iv) The ground state is created as well by the self-repulsive nonlinearity in the 3D model with the symmetry reduced from spherical to cylindrical by the strong external uniform field, which polarizes the dipole moments of the particles. (v) Very recently, we have extended the analysis for the two-component system, including miscible and immiscible cases (γ < 1 and γ > 1, respectively):

35. (4) OULTLOOK (i) In the 3D case, it is interesting to extend the analysis of the nonlinear equation for the vortical modes carrying the angular momentum. (ii) A challenging issue is to extend the analysis to the quantum gas of fermionic particles. 非常感谢您的关注!