Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2601.18764

A general variational approach for equilibrium phase boundaries of trapped spin-1 Bose-Einstein condensates

Sahil Satapathy, Projjwal Kanti Kanjilal, A. Bhattacharyay

Published Jan 26, 2026• 0 citations
arXivHTMLPDF

Abstract

We develop a simple and general variational method to estimate the solutions of the Gross-Pitaevskii equations and obtain the corresponding density profiles for all spin states of a trapped spin-1 Bose-Einstein condensate. We further employ this approach to obtain the complete phase diagram of the system under quasi-one-dimensional harmonic confinement, with ferromagnetic or antiferromagnetic spin interactions. We identify a suitable scaling that collapses all phase diagrams for different system sizes (i.e., total particle number) into a universal (system size-independent) phase diagram. The complete phase diagram for a confined system shows some significant qualitative differences compared to that of a condensate with homogeneous density distribution. The phase diagrams reported here could help identify the important parameter regimes in which phase transitions in the confined system, in general, occur. This knowledge of the region of phase boundaries can enable a reliable investigation of the instabilities near the boundaries that drive phase transitions.

DOI
10.48550/arxiv.2601.18764

Related papers

A general variational approach for equilibrium phase boundaries of trapped spin-1 Bose–Einstein condensates

Sahil Satapathy, Projjwal Kanti Kanjilal, A. Bhattacharyay

Abstract We develop a simple and general variational method to estimate the solutions of the Gross–Pitaevskii equations and obtain the corresponding density profiles for all spin states of a trapped spin-1 Bose–Einstein condensate. We further employ this approach to obtain the complete phase diagram of the system under quasi-one-dimensional harmonic confinement, with ferromagnetic or antiferromagnetic spin interactions. We identify a suitable scaling that collapses all phase diagrams for different system sizes (i.e., total particle number) into a universal (system size-independent) phase diagram. The complete phase diagram for a confined system shows some significant qualitative differences compared to that of a condensate with homogeneous density distribution. The phase diagrams reported here could help identify the important parameter regimes in which phase transitions in the confined system, in general, occur. This knowledge of the region of phase boundaries can enable a reliable investigation of the instabilities near the boundaries that drive phase transitions.

Journal of Physics B Atomic Molecular and Optical Physics•Jul 9, 2026•0 citations

Multicomponent states for trapped spin-1 Bose-Einstein condensates in the presence of a magnetic field

Projjwal Kanti Kanjilal, A. Bhattacharyay

In the presence of a magnetic field, multicomponent ground states appear in trapped spin-1 Bose-Einstein condensates for both ferromagnetic and antiferromagnetic types of spin-spin interaction. We aim to produce an accurate analytical description of the multicomponent states which is of fundamental importance. Despite being in the so-called regime of the Thomas-Fermi approximation (condensates with large particle number), the scenario of multicomponent states is problematic under this approximation due to the large variation in densities of the subcomponents. We generalize the variational method that we have introduced [P. K. Kanjilal and A. Bhattacharyay, Eur. Phys. J. Plus 137, 547 (2022)] to overcome the limitations of the Thomas-Fermi approximation. We demonstrate that the variational method is crucial in identifying multicomponent ground states. A comparison of the results of the variational method, which is multimodal by construction, with that of the single-mode approximation is also presented in this paper to demonstrate a marked improvement in accuracy over the single-mode approximation. We have also looked into the phase transition between phase-matched and polar states in a trapped condensate using the variational method and have identified substantial change in the phase boundary. The correspondence of the phase diagram of the trapped case with the homogeneous one identifies other limitations of the Thomas-Fermi approximation as opposed to the more accurate variational method.

arXiv:2301.06856•Physical Review A•Nov 27, 2023•2 citations

Two-species Bose–Einstein condensate in an optical lattice: analytical approximate formulae

R. Cipolatti, Leonardo Villegas-Lelovsky, Ming-Chiang Chung +1 more

Abstract Employing a general variational method and perturbation theory, we derived explicit solutions for the description of one-dimensional two species Bose–Einstein condensates confined by a harmonic trap potential in an optical lattice. We consider the system of two coupled Gross–Pitaevskii equations (GPE) and derive explicit expressions for the chemical potentials and wavefunctions in terms of the atom–atom interaction parameters and laser intensity. We have compared our results with the numerical solutions of the GPE and performed a quantitative analysis for the both considered methods. We underline the importance of the obtained explicit solutions to characterize the density profile or degree of miscibility of the two components.

arXiv:1508.04341•Journal of Physics A Mathematical and Theoretical•Feb 22, 2016•5 citations

Scaling equations for Bose-Einstein condensate dynamics across all interaction regimes

D. C. Marinica, C. Puertas González, T. Estrampes +5 more

We derive a unified set of scaling equations for Bose-Einstein condensates in time-dependent harmonic traps, connecting the weakly interacting Gaussian regime to the strongly interacting Thomas-Fermi regime. The exact Gross-Pitaevskii ground-state density is taken as a fixed profile in rescaled coordinates, with its dynamics described by three scaling factors corresponding to compression or expansion along the three spatial directions. The resulting equations contain no adjustable parameter, since their coefficients are determined once from the initial state. They recover analytically the Gaussian variational and Thomas-Fermi scaling equations in their respective limits and satisfy an axis-resolved virial theorem. The approach also yields the global phase and the low-lying collective-mode frequencies of the condensate. We benchmark the model against three-dimensional Gross-Pitaevskii simulations and against measured expansion energies. The model remains accurate across all interaction regimes, in strongly anisotropic traps, and along time-dependent sequences including a relatively fast quench, up to the point where the underlying frozen-profile hypothesis, shared by all three scaling approaches, breaks down.

cond-mat.quant-gasCross-listedarXiv:2609.13497•Sep 15, 2026•0 citations

Phase diagram of a rotating Bose-Einstein condensate with anharmonic confinement

Andrew David Jackson, G. M. Kavoulakis, Emil Lundh

We examine the phase diagram of an effectively repulsive Bose-Einstein condensate of atoms that rotates in a quadratic-plus-quartic potential. With use of a variational method we identify the three possible phases of the system as a function of the rotational frequency of the trap and of the coupling constant. The derived phase diagram is shown to be universal and partly exact in the limit of weak interactions and small anharmonicity. The variational results are found to be consistent with numerical solutions of the Gross-Pitaevskii equation.

arXiv:cond-mat/0401292•Physical Review A•May 18, 2004•57 citations

Stochastic Gross-Pitaevskii theory for a spin-1 Bose gas: Application to superfluidity in two dimensions

Andrew P. C. Underwood, P. B. Blakie

This paper develops and implements the stochastic projected Gross-Pitaevskii equation for spin-1 Bose gases, addressing key considerations for numerical simulations. As an application of the theory we explore equilibrium phases in a two-dimensional spin-1 gas, where quasi-long-range order emerges via a Berezinskii-Kosterlitz-Thouless transition. Our analysis includes the definition of superfluid densities for both mass and spin degrees of freedom in a manner suitable for implementation within a stochastic projected Gross-Pitaevskii equation simulation. We present a finite-temperature phase diagram for the ferromagnetic spin-1 Bose gas and identify three distinct superfluid phases: two exhibiting conventional Berezinskii-Kosterlitz-Thouless-like behavior and an additional phase that simultaneously supports independent mass and spin superflows. As temperature increases, the stability region of this additional phase shrinks. We provide a characterization of the phase transitions through consideration of the spin-component densities and the unbinding of multiple types of vortices. This work provides a foundation for further studies of nonequilibrium and finite-temperature phenomena in spinor Bose gases.

arXiv:2502.18792•Physical Review A•Jul 17, 2025•3 citations
We develop a simple and general variational method to estimate the solutions of the Gross-Pitaevskii equations and obtain the corresponding density profiles for all spin states of a trapped spin-1 Bose-Einstein condensate. We further employ this approach to obtain the complete phase diagram of the system under quasi-one-dimensional harmonic confinement, with ferromagnetic or antiferromagnetic spin interactions. We identify a suitable scaling that collapses all phase diagrams for different system sizes (i.e., total particle number) into a universal (system size-independent) phase diagram. The complete phase diagram for a confined system shows some significant qualitative differences compared to that of a condensate with homogeneous density distribution. The phase diagrams reported here could help identify the important parameter regimes in which phase transitions in the confined system, in general, occur. This knowledge of the region of phase boundaries can enable a reliable investigation of the instabilities near the boundaries that drive phase transitions.
Abstract We develop a simple and general variational method to estimate the solutions of the Gross–Pitaevskii equations and obtain the corresponding density profiles for all spin states of a trapped spin-1 Bose–Einstein condensate. We further employ this approach to obtain the complete phase diagram of the system under quasi-one-dimensional harmonic confinement, with ferromagnetic or antiferromagnetic spin interactions. We identify a suitable scaling that collapses all phase diagrams for different system sizes (i.e., total particle number) into a universal (system size-independent) phase diagram. The complete phase diagram for a confined system shows some significant qualitative differences compared to that of a condensate with homogeneous density distribution. The phase diagrams reported here could help identify the important parameter regimes in which phase transitions in the confined system, in general, occur. This knowledge of the region of phase boundaries can enable a reliable investigation of the instabilities near the boundaries that drive phase transitions.
In the presence of a magnetic field, multicomponent ground states appear in trapped spin-1 Bose-Einstein condensates for both ferromagnetic and antiferromagnetic types of spin-spin interaction. We aim to produce an accurate analytical description of the multicomponent states which is of fundamental importance. Despite being in the so-called regime of the Thomas-Fermi approximation (condensates with large particle number), the scenario of multicomponent states is problematic under this approximation due to the large variation in densities of the subcomponents. We generalize the variational method that we have introduced [P. K. Kanjilal and A. Bhattacharyay, Eur. Phys. J. Plus 137, 547 (2022)] to overcome the limitations of the Thomas-Fermi approximation. We demonstrate that the variational method is crucial in identifying multicomponent ground states. A comparison of the results of the variational method, which is multimodal by construction, with that of the single-mode approximation is also presented in this paper to demonstrate a marked improvement in accuracy over the single-mode approximation. We have also looked into the phase transition between phase-matched and polar states in a trapped condensate using the variational method and have identified substantial change in the phase boundary. The correspondence of the phase diagram of the trapped case with the homogeneous one identifies other limitations of the Thomas-Fermi approximation as opposed to the more accurate variational method.
Abstract Employing a general variational method and perturbation theory, we derived explicit solutions for the description of one-dimensional two species Bose–Einstein condensates confined by a harmonic trap potential in an optical lattice. We consider the system of two coupled Gross–Pitaevskii equations (GPE) and derive explicit expressions for the chemical potentials and wavefunctions in terms of the atom–atom interaction parameters and laser intensity. We have compared our results with the numerical solutions of the GPE and performed a quantitative analysis for the both considered methods. We underline the importance of the obtained explicit solutions to characterize the density profile or degree of miscibility of the two components.
We derive a unified set of scaling equations for Bose-Einstein condensates in time-dependent harmonic traps, connecting the weakly interacting Gaussian regime to the strongly interacting Thomas-Fermi regime. The exact Gross-Pitaevskii ground-state density is taken as a fixed profile in rescaled coordinates, with its dynamics described by three scaling factors corresponding to compression or expansion along the three spatial directions. The resulting equations contain no adjustable parameter, since their coefficients are determined once from the initial state. They recover analytically the Gaussian variational and Thomas-Fermi scaling equations in their respective limits and satisfy an axis-resolved virial theorem. The approach also yields the global phase and the low-lying collective-mode frequencies of the condensate. We benchmark the model against three-dimensional Gross-Pitaevskii simulations and against measured expansion energies. The model remains accurate across all interaction regimes, in strongly anisotropic traps, and along time-dependent sequences including a relatively fast quench, up to the point where the underlying frozen-profile hypothesis, shared by all three scaling approaches, breaks down.
We examine the phase diagram of an effectively repulsive Bose-Einstein condensate of atoms that rotates in a quadratic-plus-quartic potential. With use of a variational method we identify the three possible phases of the system as a function of the rotational frequency of the trap and of the coupling constant. The derived phase diagram is shown to be universal and partly exact in the limit of weak interactions and small anharmonicity. The variational results are found to be consistent with numerical solutions of the Gross-Pitaevskii equation.
This paper develops and implements the stochastic projected Gross-Pitaevskii equation for spin-1 Bose gases, addressing key considerations for numerical simulations. As an application of the theory we explore equilibrium phases in a two-dimensional spin-1 gas, where quasi-long-range order emerges via a Berezinskii-Kosterlitz-Thouless transition. Our analysis includes the definition of superfluid densities for both mass and spin degrees of freedom in a manner suitable for implementation within a stochastic projected Gross-Pitaevskii equation simulation. We present a finite-temperature phase diagram for the ferromagnetic spin-1 Bose gas and identify three distinct superfluid phases: two exhibiting conventional Berezinskii-Kosterlitz-Thouless-like behavior and an additional phase that simultaneously supports independent mass and spin superflows. As temperature increases, the stability region of this additional phase shrinks. We provide a characterization of the phase transitions through consideration of the spin-component densities and the unbinding of multiple types of vortices. This work provides a foundation for further studies of nonequilibrium and finite-temperature phenomena in spinor Bose gases.