- Open Access
Generalized Gross-Pitaevskii Equation for 2D Bosons with Attractive Interactions
Phys. Rev. Lett. 136, 243402 – Published 18 June, 2026
DOI: https://doi.org/10.1103/2831-t3jk
Abstract
We introduce a generalized Gross-Pitaevskii equation that provides a nonlinear framework for studying 2D attractive Bose systems. Its defining feature is the logarithmic density dependence of the coupling constant, which breaks the scale invariance inherent in the standard mean-field equations. This framework allows straightforward calculations of the system properties arising from the quantum anomaly. As a first illustration, we study universal bound states in free space, commonly referred to as quantum droplets. Then, we analyze breathing modes and quench dynamics in trapped systems, paving the way for a systematic exploration of nonequilibrium phenomena in 2D attractive Bose systems. Finally, we predict the existence of universal excited states, including vortex configurations, which may be more accessible to experimental investigation than the ground state. Our results provide a robust theoretical foundation for studying both static and dynamical properties of finite systems, and offer guidance for the design of future experiments.
Physics Subject Headings (PhySH)
Article Text
Supplemental Material
References (72)
- R. Y. Chiao, E. Garmire, and C. H. Townes, Phys. Rev. Lett. 13, 479 (1964).
- V. E. Zakharov, Sov. Phys. JETP 35, 908 (1972).
- V. E. Zakharov and E. A. Kuznetsov, Zh. Eksp. Teor. Fiz. 91, 1310 (1986).
- L. Pitaevskii, Phys. Lett. A 221, 14 (1996).
- L. P. Pitaevskii and A. Rosch, Phys. Rev. A 55, R853 (1997).
- R. Jackiw, Delta Function Potentials in Two-Dimensional and Three-Dimensional Quantum Mechanics, edited by A. Ali and P. Hoodbhoy (World Scientific, Singapore, 1991), pp. 35–53.
- M. Olshanii, H. Perrin, and V. Lorent, Phys. Rev. Lett. 105, 095302 (2010).
- H.-W. Hammer and D. T. Son, Phys. Rev. Lett. 93, 250408 (2004).
- S. G. Rajeev, arXiv:hep-th/9905120.
- S. G. Rajeev, arXiv:hep-th/9902025.
- Z.-H. Luo, W. Pang, B. Liu, Y.-Y. Li, and B. A. Malomed, Front. Phys. 16, 32201 (2020).
- L. Chomaz, I. Ferrier-Barbut, F. Ferlaino, B. Laburthe-Tolra, B. L. Lev, and T. Pfau, Rep. Prog. Phys. 86, 026401 (2022).
- S. Mistakidis, A. Volosniev, R. Barfknecht, T. Fogarty, T. Busch, A. Foerster, P. Schmelcher, and N. Zinner, Phys. Rep. 1042, 1 (2023).
- D. S. Petrov, SciPost Phys. 19, 037 (2025).
- A. Tononi, G. E. Astrakharchik, and D. S. Petrov, AVS Quantum Sci. 6, 023201 (2024).
- F. Brauneis, A. G. Volosniev, and H.-W. Hammer, Phys. Rev. A 111, 063308 (2025).
- K. Iida and G. Baym, Phys. Rev. D 66, 014015 (2002).
- M. Eto, M. Nitta, and N. Yamamoto, Phys. Rev. Lett. 104, 161601 (2010).
- D. J. Gross and F. Wilczek, Phys. Rev. Lett. 30, 1343 (1973).
- H. D. Politzer, Phys. Rev. Lett. 30, 1346 (1973).
- J. C. Collins and M. J. Perry, Phys. Rev. Lett. 34, 1353 (1975).
- A. Kurkela, P. Romatschke, and A. Vuorinen, Phys. Rev. D 81, 105021 (2010).
- M. Schick, Phys. Rev. A 3, 1067 (1971).
- S. I. Shevchenko, Sov. J. Low Temp. Phys. 18, 223 (1992).
- E. H. Lieb, R. Seiringer, and J. Yngvason, Commun. Math. Phys. 224, 17 (2001).
- E. B. Kolomeisky, T. J. Newman, J. P. Straley, and X. Qi, Phys. Rev. Lett. 85, 1146 (2000).
- J. Levinsen and M. M. Parish, Strongly interacting two-dimensional Fermi gases, in Annual Review of Cold Atoms and Molecules (World Scientific, Singapore, 2015), pp. 1–75.
- See Supplemental Material at http://link.aps.org/supplemental/10.1103/2831-t3jk for (S1) additional information about the flow-equations method, (S2) description of the numerical method used to solve the generalized GPE, (S3) additional information about variational data presented Fig. 2(b), (S4) derivative of the energy presented in Fig. 2(b), and (S5) detailed derivation of Eqs. (13) and (14). The Supplemental Material contains additional Refs. [29–43].
- A. G. Volosniev and H.-W. Hammer, Phys. Rev. A 96, 031601 (2017).
- F. Brauneis, T. G. Backert, S. I. Mistakidis, M. Lemeshko, H.-W. Hammer, and A. G. Volosniev, New J. Phys. 24, 063036 (2022).
- F. Brauneis, A. Ghazaryan, H.-W. Hammer, and A. G. Volosniev, Commun. Phys. 6, 224 (2023).
- F. Brauneis, H.-W. Hammer, S. M. Reimann, and A. G. Volosniev, Phys. Rev. A 111, 013303 (2025).
- W. Bao and Q. Du, SIAM J. Sci. Comput. 25, 1674 (2004).
- X. Antoine and R. Duboscq, Comput. Phys. Commun. 185, 2969 (2014).
- Y. Saad, Iterative Methods for Sparse Linear Systems, 2nd ed., Other Titles in Applied Mathematics (SIAM, Philadelphia, 2003).
- Krylovkit.jl, https://jutho.github.io/KrylovKit.jl (2022).
- W. Bao and H. Wang, J. Comput. Phys. 217, 612 (2006).
- G. M. Koutentakis, S. I. Mistakidis, and P. Schmelcher, Atoms 10, 3 (2022).
- O. Bohigas, A. Lane, and J. Martorell, Phys. Rep. 51, 267 (1979).
- F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 71, 463 (1999).
- L. Pitaevskii and S. Stringari, Bose-Einstein Condensation and Superfluidity (Oxford University Press, New York, 2016).
- H. E. Camblong, A. Chakraborty, P. Lopez-Duque, and C. R. Ordóñez, J. Math. Phys. (N.Y.) 64, 092302 (2023).
- J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed. (Cambridge University Press, Cambridge, England, 2020).
is a parameter that we shall use to quantify the interaction strength. For simplicity, we omit the parentheses in the numerator and write it simply as in what follows.
- B. Bazak and D. S. Petrov, New J. Phys. 20, 023045 (2018).
It may be tempting to disregard the term due to its typically minor energy contribution, just as it is typically done for repulsive systems [47]. However, in finite systems and time-dependent scenarios, it is essential for preserving energy conservation.
- A. Y. Cherny and A. A. Shanenko, Phys. Rev. E 64, 027105 (2001).
- B. Bakkali-Hassani and J. Dalibard, in Proceedings of the International School of Physics “Enrico Fermi”, Course 211—Quantum Mixtures with Ultra-Cold Atoms, edited by R. Grimm, M. Inguscio, and S. Stringari (IOS Press, Amsterdam, 2024).
In practice, calculations are performed in the presence of a very weak external potential, which enhances the stability of our numerical approach. See Supplemental Material for further information about the numerical method [28].
- I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008).
- S. Kehrein, The Flow Equation Approach to Many-Particle Systems (Springer, Berlin Heidelberg, 2006).
- K. Tsukiyama, S. K. Bogner, and A. Schwenk, Phys. Rev. Lett. 106, 222502 (2011).
- A. G. Volosniev and H.-W. Hammer, New J. Phys. 19, 113051 (2017).
- F. Brauneis, H.-W. Hammer, M. Lemeshko, and A. Volosniev, SciPost Phys. 11, 008 (2021).
This weak dependence stands in contrast to the behavior observed in a bubble trap [15], which features a transition from a uniform density to a localized droplet. In that case, the transition involves spontaneous symmetry breaking, potentially leading to dynamics that differ significantly from those presented here. The existence of a vanishing energy Bogoliubov mode in Eq. (3) for a bubble trap confirms this assertion. (For a homogeneous density, the Bogoliubov analysis can be performed as in Ref. [15]).
- Y. Guo and R. Seiringer, Lett. Math. Phys. 104, 141 (2013).
For an experimental validation of this result, the quasi-2D nature of realistic cold-atom setups must be taken into account. Effects arising from the quasi-2D geometry also break the SO(2,1) symmetry, leading to additional corrections to the breathing-mode frequency [7, 58]. These effects are, however, extrinsic and can be disentangled from the intrinsic symmetry breaking discussed here by varying the trapping frequency along the tightly confined direction.
- Y. Hu and Z. Liang, Phys. Rev. Lett. 107, 110401 (2011).
- M. D. Lee and S. A. Morgan, J. Phys. B 35, 3009 (2002).
- M. S. Hasan, T. Fogarty, J. Li, A. Ruschhaupt, and T. Busch, Phys. Rev. Res. 6, 023114 (2024).
- A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalattore, Rev. Mod. Phys. 83, 863 (2011).
Notably, it is already known that the 2D three-body (four-body) system supports a single excited state with zero angular momentum and energy [63, 64] ( [65]). According to Eq. (3) one could expect that the large- limit for such excited states reads , where corresponds to the excited state of the generalized GPE with . It is also plausible that these excited states are beyond the reach of our mean-field-like approach. Classifying such states in 2D attractive systems presents a compelling direction for future research.
- L. W. Bruch and J. A. Tjon, Phys. Rev. A 19, 425 (1979).
- E. Nielsen, D. V. Fedorov, and A. S. Jensen, Few-Body Syst. 27, 15 (1999).
- L. Platter, H.-W. Hammer, and U.-G. Meissner, Few-Body Syst. 35, 169 (2004).
- H. A. Haus, Appl. Phys. Lett. 8, 128 (1966).
- V. Kruglov, Y. Logvin, and V. Volkov, J. Mod. Opt. 39, 2277 (1992).
- M. Suchorowski, F. Brauneis, H.-W. Hammer, M. Tomza, and A. Volosniev, 10.5281/zenodo.20275167 (2026).
- M. Suchorowski, Julia Generalized GPE Solver for 2D Attractive Bose Gas, 10.5281/zenodo.17593677 (2025).
Note that in many sources, including Ref. [14], is denoted as , unlike the notation used in the main text of this Letter.
- V. Efimov, Phys. Lett. 33B, 563 (1970).
- E. Braaten and H. W. Hammer, Phys. Rep. 428, 259 (2006).