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

New Journal of Physics

Classical field simulation of vortex lattice melting in a two-dimensional fast rotating Bose gas

Sálvio Jacob Bereta, Lucas Madeira, M. A. Caracanhas, Hélène Perrin, Romain Dubessy

Published Jun 1, 2026• 0 citations
Publisher

Abstract

Abstract We present a classical field simulation study of the thermal melting of a two-dimensional (2D) vortex lattice in a rotating Bose gas, focusing on the role of finite-size effects on the melting temperature. This work constitutes a numerical continuation of the recent experimental investigation reported in (Sharma 2024 Phys. Rev. Lett. 133 143401), which addressed the thermal melting of a vortex lattice in a quasi-2D Bose gas. Using the stochastic projected Gross–Pitaevskii equation in a harmonic plus quartic trap, we simulate the finite-temperature equilibrium state and extract vortex configurations from density snapshots. Clear signatures of the two-step Kosterlitz–Thouless–Halperin–Nelson–Young melting scenario are identified. Our simulations enable a detailed characterization of the crystalline, hexatic, and liquid phases through correlation functions quantifying the translational and orientational order and through defect statistics. Finite-size effects are shown to play a crucial role at lower rotation frequencies, affecting the proliferation of lattice defects.

DOI
10.1088/1367-2630/ae7ab2

Related papers

Classical field simulation of vortex lattice melting in a two-dimensional fast rotating Bose gas

Sálvio Jacob Bereta, Lucas Madeira, M. A. Caracanhas +2 more

We present a classical field simulation study of the thermal melting of a two-dimensional vortex lattice in a rotating Bose gas, focusing on the role of finite-size effects on the melting temperature. This work constitutes a numerical continuation of the recent experimental investigation reported in [Physical Review Letters 133, 143401 (2024)], which addressed the thermal melting of a vortex lattice in a quasi-two-dimensional Bose gas. Using the stochastic projected Gross-Pitaevskii equation in a harmonic plus quartic trap, we simulate the finite-temperature equilibrium state and extract vortex configurations from density snapshots. Clear signatures of the two-step Kosterlitz--Thouless--Halperin--Nelson--Young melting scenario are identified. Our simulations enable a detailed characterization of the crystalline, hexatic, and liquid phases through correlation functions quantifying the translational and orientational order and through defect statistics. Finite-size effects are shown to play a crucial role at lower rotation frequencies, affecting the proliferation of lattice defects.

arXiv:2603.01911•Mar 2, 2026•0 citations

Thermal Melting of a Vortex Lattice in a Quasi-Two-Dimensional Bose Gas

Rishabh Sharma, David Rey, Laurent Longchambon +3 more

We report the observation of the melting of a vortex lattice in a fast rotating quasi-two dimensional Bose gas, under the influence of thermal fluctuations. We image the vortex lattice after a time-of-flight expansion, for increasing rotation frequency at constant atom number and temperature. We detect the vortex positions and study the order of the lattice using the pair correlation function and the orientational correlation function. We evidence the melting transition by a change in the decay of orientational correlations, associated with a proliferation of dislocations. Our findings are consistent with the hexatic to liquid transition in the Kosterlitz, Thouless, Halperin, Nelson, and Young scenario for two-dimensional melting.

arXiv:2404.05460•Physical Review Letters•Sep 30, 2024•13 citations

Vortex dynamics on a shell-trapped Bose-Einstein condensate with dipole-dipole interactions and thermal melting of a two-dimensional vortex lattice

Sálvio Jacob Bereta

There are many efforts to understand better the dipolar long-range interaction, the shell potential confinement, and the vortices in Bose-Einstein condensates (BECs). Furthermore, combining these three areas suggests that it can give rise to new phenomena. In this thesis, we studied the dynamics of vortices in a shell-trapped dipolar BEC. The long-range dipolar interaction induces an inhomogeneous density distribution on the shell, affecting the vortex velocity field. We analyzed the dynamics of one and two pairs of vortices symmetrically located in the shell to separate the vortex-vortex interaction contribution from that of the density gradient due to the inhomogeneous background. In addition to the numerical studies, a Lagrangian formalism allowed us to derive semi-analytical expressions which are in agreement with the results of the Gross-Pitaevskii equation (GPE) simulations, giving us better insights into the influence of the core in the vortex dynamics outside the validity of the point-vortex model. In the second part of this thesis, we describe an actual experimental implementation of the shell BEC. Using the Stochastic Projected GPE simulations, we determine the configuration of a finite temperature vortex lattice produced in this potential over the action of gravity. More specifically, we studied the transition from a perfectly hexagonal (triangular) lattice to the hexatic and liquid phases by detecting the vortex positions to calculate the angular and pair correlation functions of this system for different values of temperatures and trap rotation frequencies. Our results for the vortex lattice melting are consistent with the hexatic to liquid transition in the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) scenario and can be extended to predict the melting temperature of other finite-size systems.

•Oct 30, 2024•0 citations

Bose-Einstein condensation from a rotating thermal cloud: Vortex nucleation and lattice formation

Ashton S. Bradley, Crispin W. Gardiner, Matthew J. Davis

We develop a stochastic Gross-Pitaveskii theory suitable for the study of Bose-Einstein condensation in a rotating dilute Bose gas. The theory is used to model the dynamical and equilibrium properties of a rapidly rotating Bose gas quenched through the critical point for condensation, as in the experiment of Haljan et al. [Phys. Rev. Lett. 87, 210403 (2001)]. In contrast to stirring a vortex-free condensate, where topological constraints require that vortices enter from the edge of the condensate, we find that phase defects in the initial noncondensed cloud are trapped en masse in the emerging condensate. Bose-stimulated condensate growth proceeds into a disordered vortex configuration. At sufficiently low temperature the vortices then order into a regular Abrikosov lattice in thermal equilibrium with the rotating cloud. We calculate the effect of thermal fluctuations on vortex ordering in the final gas at different temperatures, and find that the BEC transition is accompanied by lattice melting associated with diminishing long-range correlations between vortices across the system.

arXiv:0712.3436•Physical Review A•Mar 14, 2008•108 citations

Vortex pairing in two-dimensional Bose gases

Christopher J. Foster, P. B. Blakie, Matthew J. Davis

Recent experiments on ultracold Bose gases in two dimensions have provided evidence for the existence of the Berezinskii-Kosterlitz-Thouless (BKT) phase via analysis of the interference between two independent systems. In this work we study the two-dimensional quantum degenerate Bose gas at finite temperature using the projected Gross-Pitaevskii equation classical field method. Although this describes the highly occupied modes of the gas below a momentum cutoff, we have developed a method to incorporate the higher momentum states in our model. We concentrate on finite-sized homogeneous systems in order to simplify the analysis of the vortex pairing. We determine the dependence of the condensate fraction on temperature and compare this to the calculated superfluid fraction. By measuring the first order correlation function we determine the boundary of the Bose-Einstein condensate and BKT phases, and find it is consistent with the superfluid fraction decreasing to zero. We reveal the characteristic unbinding of vortex pairs above the BKT transition via a coarse-graining procedure. Finally, we model the procedure used in experiments to infer system correlations [Hadzibabic et al., Nature 441, 1118 (2006)], and quantify its level of agreement with directly calculated in situ correlation functions.

arXiv:0912.1675•Physical Review A•Feb 24, 2010•64 citations

Quantized vortices in a Bose-Einstein condensate: thermal activation and dynamic nucleation

Sabine Stock

This thesis manuscript presents experimental results on Bose-Einstein condensates in optical lattices and in the regime of fast rotation. In an optical lattice an array of inde-pendent 2-dimensional Bose-Einstein condensates is produced from an initially 3-dimen-sional cloud. The number of 2D BECs in this array can be controlled by superimposing a magnetic field gradient along the axis of the lattice and then selectively evaporating the sites. When producing a small number of sites and letting them interfere, dislocations in the interference pattern were observed, which must be due to nontrivial phase defects, as for example free vortices. Such vortices and vortex-antivortex pairs play an important role in the superfluid to normal phase transition in a 2-dimensional gas, the Kosterlitz-Thouless transition. With 3-dimensional condensates the limit of fast rotation is studied in a hybrid potential contain-ing a small quartic term. In this trap rotation frequencies larger than the harmonic trapping fre-quency can be reached. These experiments are a first step to a possible observation of strongly correlated states close to the regime of the fractional quantum Hall effect.

•Jan 13, 2006•4 citations
Abstract We present a classical field simulation study of the thermal melting of a two-dimensional (2D) vortex lattice in a rotating Bose gas, focusing on the role of finite-size effects on the melting temperature. This work constitutes a numerical continuation of the recent experimental investigation reported in (Sharma 2024 Phys. Rev. Lett. 133 143401), which addressed the thermal melting of a vortex lattice in a quasi-2D Bose gas. Using the stochastic projected Gross–Pitaevskii equation in a harmonic plus quartic trap, we simulate the finite-temperature equilibrium state and extract vortex configurations from density snapshots. Clear signatures of the two-step Kosterlitz–Thouless–Halperin–Nelson–Young melting scenario are identified. Our simulations enable a detailed characterization of the crystalline, hexatic, and liquid phases through correlation functions quantifying the translational and orientational order and through defect statistics. Finite-size effects are shown to play a crucial role at lower rotation frequencies, affecting the proliferation of lattice defects.
We present a classical field simulation study of the thermal melting of a two-dimensional vortex lattice in a rotating Bose gas, focusing on the role of finite-size effects on the melting temperature. This work constitutes a numerical continuation of the recent experimental investigation reported in [Physical Review Letters 133, 143401 (2024)], which addressed the thermal melting of a vortex lattice in a quasi-two-dimensional Bose gas. Using the stochastic projected Gross-Pitaevskii equation in a harmonic plus quartic trap, we simulate the finite-temperature equilibrium state and extract vortex configurations from density snapshots. Clear signatures of the two-step Kosterlitz--Thouless--Halperin--Nelson--Young melting scenario are identified. Our simulations enable a detailed characterization of the crystalline, hexatic, and liquid phases through correlation functions quantifying the translational and orientational order and through defect statistics. Finite-size effects are shown to play a crucial role at lower rotation frequencies, affecting the proliferation of lattice defects.
We report the observation of the melting of a vortex lattice in a fast rotating quasi-two dimensional Bose gas, under the influence of thermal fluctuations. We image the vortex lattice after a time-of-flight expansion, for increasing rotation frequency at constant atom number and temperature. We detect the vortex positions and study the order of the lattice using the pair correlation function and the orientational correlation function. We evidence the melting transition by a change in the decay of orientational correlations, associated with a proliferation of dislocations. Our findings are consistent with the hexatic to liquid transition in the Kosterlitz, Thouless, Halperin, Nelson, and Young scenario for two-dimensional melting.
There are many efforts to understand better the dipolar long-range interaction, the shell potential confinement, and the vortices in Bose-Einstein condensates (BECs). Furthermore, combining these three areas suggests that it can give rise to new phenomena. In this thesis, we studied the dynamics of vortices in a shell-trapped dipolar BEC. The long-range dipolar interaction induces an inhomogeneous density distribution on the shell, affecting the vortex velocity field. We analyzed the dynamics of one and two pairs of vortices symmetrically located in the shell to separate the vortex-vortex interaction contribution from that of the density gradient due to the inhomogeneous background. In addition to the numerical studies, a Lagrangian formalism allowed us to derive semi-analytical expressions which are in agreement with the results of the Gross-Pitaevskii equation (GPE) simulations, giving us better insights into the influence of the core in the vortex dynamics outside the validity of the point-vortex model. In the second part of this thesis, we describe an actual experimental implementation of the shell BEC. Using the Stochastic Projected GPE simulations, we determine the configuration of a finite temperature vortex lattice produced in this potential over the action of gravity. More specifically, we studied the transition from a perfectly hexagonal (triangular) lattice to the hexatic and liquid phases by detecting the vortex positions to calculate the angular and pair correlation functions of this system for different values of temperatures and trap rotation frequencies. Our results for the vortex lattice melting are consistent with the hexatic to liquid transition in the Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) scenario and can be extended to predict the melting temperature of other finite-size systems.
We develop a stochastic Gross-Pitaveskii theory suitable for the study of Bose-Einstein condensation in a rotating dilute Bose gas. The theory is used to model the dynamical and equilibrium properties of a rapidly rotating Bose gas quenched through the critical point for condensation, as in the experiment of Haljan et al. [Phys. Rev. Lett. 87, 210403 (2001)]. In contrast to stirring a vortex-free condensate, where topological constraints require that vortices enter from the edge of the condensate, we find that phase defects in the initial noncondensed cloud are trapped en masse in the emerging condensate. Bose-stimulated condensate growth proceeds into a disordered vortex configuration. At sufficiently low temperature the vortices then order into a regular Abrikosov lattice in thermal equilibrium with the rotating cloud. We calculate the effect of thermal fluctuations on vortex ordering in the final gas at different temperatures, and find that the BEC transition is accompanied by lattice melting associated with diminishing long-range correlations between vortices across the system.
Recent experiments on ultracold Bose gases in two dimensions have provided evidence for the existence of the Berezinskii-Kosterlitz-Thouless (BKT) phase via analysis of the interference between two independent systems. In this work we study the two-dimensional quantum degenerate Bose gas at finite temperature using the projected Gross-Pitaevskii equation classical field method. Although this describes the highly occupied modes of the gas below a momentum cutoff, we have developed a method to incorporate the higher momentum states in our model. We concentrate on finite-sized homogeneous systems in order to simplify the analysis of the vortex pairing. We determine the dependence of the condensate fraction on temperature and compare this to the calculated superfluid fraction. By measuring the first order correlation function we determine the boundary of the Bose-Einstein condensate and BKT phases, and find it is consistent with the superfluid fraction decreasing to zero. We reveal the characteristic unbinding of vortex pairs above the BKT transition via a coarse-graining procedure. Finally, we model the procedure used in experiments to infer system correlations [Hadzibabic et al., Nature 441, 1118 (2006)], and quantify its level of agreement with directly calculated in situ correlation functions.
This thesis manuscript presents experimental results on Bose-Einstein condensates in optical lattices and in the regime of fast rotation. In an optical lattice an array of inde-pendent 2-dimensional Bose-Einstein condensates is produced from an initially 3-dimen-sional cloud. The number of 2D BECs in this array can be controlled by superimposing a magnetic field gradient along the axis of the lattice and then selectively evaporating the sites. When producing a small number of sites and letting them interfere, dislocations in the interference pattern were observed, which must be due to nontrivial phase defects, as for example free vortices. Such vortices and vortex-antivortex pairs play an important role in the superfluid to normal phase transition in a 2-dimensional gas, the Kosterlitz-Thouless transition. With 3-dimensional condensates the limit of fast rotation is studied in a hybrid potential contain-ing a small quartic term. In this trap rotation frequencies larger than the harmonic trapping fre-quency can be reached. These experiments are a first step to a possible observation of strongly correlated states close to the regime of the fractional quantum Hall effect.