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

Vortex transport and screening in a closed two-dimensional Bose field: friction, transverse force and diffusion of quantized vortices, and the polarisability of the thermal vortex-pair gas, measured in an energy-conserving classical field with no reservoir and no fit parameter (version 2, revised after peer review)

Xavier Callens

Published Oct 7, 2026• 0 citations
Publisher

Abstract

A quantized vortex in a thermal two-dimensional Bose gas is dragged by the normal fluid (friction α), pushed sideways by it (transverse coefficient α′) and kicked by it (diffusion η). We measure the three coefficients in the projected Gross-Pitaevskii equation — microcanonical, no damping term, no noise put in by hand — with sub-grid vortex tracking and estimators that are identities of the dissipative point-vortex model, machine-checked in Lean 4. Every hypothesis, gate and decision rule was committed before the corresponding run; the instrument failed its T = 0 control twice before it passed, an earlier gate verdict built on the flawed instrument is withdrawn, and both are reported. Results. At T/TBKT = 0.14, 0.27, 0.43: α = 0.0062(4), 0.014(2), 0.02(1), size- and estimator-independent, α ≈ 0.24 ρn/ρ. α′ is zero within 1% once the T = 0 compressible pair speed is subtracted; the Iordanskii value ρn/ρ is excluded at 18σ and the registered Iordanskii-like prediction is refuted. The Einstein relation is inconclusive by the registered rule (one admitted temperature at ηK/α = 2.1 ± 0.7; informally 1.4-3x, not 100x). Single pairs stop shrinking while the phonon band takes the impulse they shed (a self-generated counterflow with a published precedent in the truncated GPE); the decisive box-size control is not run. The single non-equilibrium pair that corrupts the superfluid-stiffness estimator in finite boxes decays by 4-27x given Lz-scaled time (six of six L = 192 trajectories admitted, one unfinished at t = 14 500). The pre-registered mode-by-mode dielectric relation fails in the thermal pair gas (coherence 0.06-0.28), while the vortex-charge structure factor equals the Kosterlitz-Thouless polarisability Σd2 from the measured pair sizes to 5-21%. The Josephson offset ηK − 1 = +0.3 survives at L = 192 in equilibrium, independent of fit window and torus geometry. Formal companions (Lean 4 / Mathlib, standard axioms): energy dissipation of any skew-plus-gradient vortex flow, the d2 law and lifetime bound, the centre-velocity identity, the wind stall, the zero-flux/Einstein equivalence, the pair-polarisation form and response ceiling, the charge-density/matching bound, the Kosterlitz-flow invariant, the bounded increments of quasi-periodic signals. Pre-registrations, amendments, gate records, data and code are in the software repository (concept DOI 10.5281/zenodo.22855581). No physics is claimed as new; the numbers, the refutations and the method are.

DOI
10.5281/zenodo.23202454

Related papers

The friction on a quantized vortex in a truncated Gross-Pitaevskii field follows the temperature, not the normal density: a pre-registered test of alpha proportional to rho_n against the cutoff, and the correction of a published reading

Xavier Callens

The mutual friction alpha on a quantized vortex in a two-dimensional Bose field is commonly taken to scale with the normal fraction rho_n/rho. In an earlier paper (doi:10.5281/zenodo.23225293) we measured alpha(T) in the projected Gross-Pitaevskii equation at one cutoff and found alpha ~ 0.24 rho_n/rho; we read the coefficient as a mode-averaged transport cross-section of the vortex (1.4 healing lengths) and registered its cutoff dependence as the open question, with the Born phonon-scattering law (coefficient proportional to the cutoff) as the rival. Here the same energy-conserving instrument, with its T = 0 known answer reproduced to 0.3 % at each cutoff, is run at two further cutoffs (k_c xi = 2 pi/3 and 2 pi). The registered prediction of a cutoff-independent coefficient fails and so does its rival: the coefficient is 0.31 +- 0.03, 0.232 +- 0.014, 0.102 +- 0.011 at k_c xi = 2.1, 3.1, 6.3. What does not move is alpha itself: alpha/T = 0.054 +- 0.003 over six arms and three cutoffs (chi2 = 1.3 for 5 degrees of freedom; 0.060 +- 0.002 with the regression estimator, because the energy estimator, the primary one, is biased low by diffusion -- a flaw of the registered synthetic gate exposed by a second, Rust implementation of the estimators). The proportionality to rho_n was a degeneracy of a single cutoff, where rho_n is proportional to T; the 'geometric core size' reading is withdrawn. The temperature law was found after the data were read; it is registered as a hypothesis, with its prediction for the next arm fixed before that arm was read (reported in the paper). Pre-registrations, controls, scripts, analysis file and Lean companion (FrictionKinetic.lean) are in the software record.

•Oct 9, 2026•0 citations

C-field description of the quasi-equilibrium Bose gas: Persistent currents and kelvon-induced vortex decay

S. J. Rooney

Recent experiments at the University of Arizona have demonstrated superb control and imaging of flattened superfluids revealing a rich landscape of vortex dynamics. Vortex dipoles have been created and their dynamics observed [Neely et al. Phys. Rev. Lett. 104, 160410 (2010)], and highly charged persistent currents have been formed through a clever combination of stirring and dissipation [T. W. Neely, PhD thesis, (2010)]. Many ultra-cold gas systems, including these, are in a regime that we may characterize as quasi-equilibrium: despite some external forcing, a large fraction of the gas remains in thermal equilibrium. In this thesis we use the stochastic projected Gross-Pitaevskii equation (SPGPE) to make quantitative predictions of the quasi-equilibrium Bose gas. \n\nWe quantitatively compare the SPGPE with experiment by modeling the persistent current formation experiment performed at the University of Arizona [T. W. Neely, PhD thesis, (2010)]. We determine all SPGPE parameters for the toroidal system prior to simulation, enabling quantitative modeling of the experiment with no fitting parameters, giving a true test of the SPGPE theory. We find the SPGPE gives quantitative agreement with experiment, accurately predicting the size of the persistent current as well as the decay time of the vortices. We also observe the crucial role that thermal fluctuations have on enabling this agreement, showing the comprehensive SPGPE treatment is necessary to make quantitatively accurate calculations of quasi-equilibrium systems. This is the first quantitative agreement with experimental observations of vortex dynamics obtained with a theory of dissipation from first principles using no fitting parameters.\n\nWe then systematically study the Kelvin mode excitations on a vortex line in a three dimensional Bose-Einstein condensate using the SPGPE. We give a quantitative measure of the magnitude of vortex bending caused by the activation of Kelvin modes, and find that vortex bending can be suppressed by tightening the confinement along the direction of the vortex line. This leads to a strong suppression of the vortex decay rate as the system enters a regime of two-dimensional vortex dynamics, characterized by a critical oblateness.\n\nIn our final application of the SPGPE we simulate the decay of a vortex dipole. We model a recent experiment of Neely et al. [Neely et al. Phys. Rev. Lett. 104, 160410 (2010)], finding the SPGPE predicts a dipole lifetime consistent with experimental observations. We also show the experiment lies in the two-dimensional regime of vortex dynamics, validating our critical oblateness calculation.

•Jan 1, 2012•0 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

Vortex motion in trapped Bose--Einstein condensates

Enikő J. M. Madarassy, Carlo Barenghi

We have performed numerical simulations of various vortex configurations in a trapped BoseEinstein condensate by solving the two-dimensional Gross-Pitaevskii equation in the presence of a simple model of interaction between the condensate and the finite temperature thermal cloud that surrounds it. In that interaction the non-condensed thermal cloud acts as a source of dissipation with a damping effect of excitations. In the case ora single vortex and a vortex - anti vortex pair, we have found that the path of the vortices depends on the initial position, the initial separation distance if the case of two vortices and dissipation. This motion is periodic and it was found that sound waves are created by vortex motion; the intensity was stronger when the initial vortex separation distance was smaller. We have calculated the sound energy as the difference between the kinetic energy and the vortex energy. With no dissipation the vortices followed the same path with a slight oscillation due to the sound waves. We found that the smaller the initial vortex separation distance do is, the larger the sound production. The period, frequency, translation speed, sound energy and vortex energy were measured for different initial separation distances do and for different dissipation parameters 'Y. In the case of motion of one vortex, the connection between the dissipation 'Y and the friction coefficients, a and a' was studied as well. To create a simple turbulent state, we put eight pairs of vortex - anti vortex at random positions in the condensate with initial separation distance do =1.8 between them. 'Ve have studied the decay rate of the total energy, kinetic energy, quantum energy, trap energy and the z - component of the angular momentum together with the increase rate of the internal energy. Finally, we finished our investigation by putting randomly vortex - anti vortex pairs and studied the decrease of the number of vortices with time t. We found that the decrease is exponential.

•Aug 1, 2008•0 citations

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

Mutual friction and diffusion of two-dimensional quantum vortices

Zain Mehdi, J. J. Hope, Stuart S. Szigeti +1 more

Dissipation of quantum vortex motion is fundamental to superfluid dynamics and quantum turbulence, yet there is currently a large gap between theory and experiments with ultracold atoms. Here we present a microscopic open quantum systems theory of thermally damped vortex motion in oblate atomic superfluids that includes previously neglected energy-damping interactions between superfluid and thermal atoms. This mechanism couples strongly to vortex core motion and causes dissipation of vortex energy due to mutual friction, as well as Brownian motion of vortices due to thermal fluctuations. We derive an analytic expression for the dimensionless mutual friction coefficient that gives excellent quantitative agreement with experimentally measured values, without any fitted parameters. Our work closes an existing two orders of magnitude gap between dissipation theory and experiments, previously bridged by fitted parameters, and provides a microscopic origin for the mutual friction and diffusion of quantized vortices in two-dimensional atomic superfluids.

arXiv:2205.04065•Physical Review Research•Mar 16, 2023•10 citations
A quantized vortex in a thermal two-dimensional Bose gas is dragged by the normal fluid (friction α), pushed sideways by it (transverse coefficient α′) and kicked by it (diffusion η). We measure the three coefficients in the projected Gross-Pitaevskii equation — microcanonical, no damping term, no noise put in by hand — with sub-grid vortex tracking and estimators that are identities of the dissipative point-vortex model, machine-checked in Lean 4. Every hypothesis, gate and decision rule was committed before the corresponding run; the instrument failed its T = 0 control twice before it passed, an earlier gate verdict built on the flawed instrument is withdrawn, and both are reported. Results. At T/TBKT = 0.14, 0.27, 0.43: α = 0.0062(4), 0.014(2), 0.02(1), size- and estimator-independent, α ≈ 0.24 ρn/ρ. α′ is zero within 1% once the T = 0 compressible pair speed is subtracted; the Iordanskii value ρn/ρ is excluded at 18σ and the registered Iordanskii-like prediction is refuted. The Einstein relation is inconclusive by the registered rule (one admitted temperature at ηK/α = 2.1 ± 0.7; informally 1.4-3x, not 100x). Single pairs stop shrinking while the phonon band takes the impulse they shed (a self-generated counterflow with a published precedent in the truncated GPE); the decisive box-size control is not run. The single non-equilibrium pair that corrupts the superfluid-stiffness estimator in finite boxes decays by 4-27x given Lz-scaled time (six of six L = 192 trajectories admitted, one unfinished at t = 14 500). The pre-registered mode-by-mode dielectric relation fails in the thermal pair gas (coherence 0.06-0.28), while the vortex-charge structure factor equals the Kosterlitz-Thouless polarisability Σd2 from the measured pair sizes to 5-21%. The Josephson offset ηK − 1 = +0.3 survives at L = 192 in equilibrium, independent of fit window and torus geometry. Formal companions (Lean 4 / Mathlib, standard axioms): energy dissipation of any skew-plus-gradient vortex flow, the d2 law and lifetime bound, the centre-velocity identity, the wind stall, the zero-flux/Einstein equivalence, the pair-polarisation form and response ceiling, the charge-density/matching bound, the Kosterlitz-flow invariant, the bounded increments of quasi-periodic signals. Pre-registrations, amendments, gate records, data and code are in the software repository (concept DOI 10.5281/zenodo.22855581). No physics is claimed as new; the numbers, the refutations and the method are.
The mutual friction alpha on a quantized vortex in a two-dimensional Bose field is commonly taken to scale with the normal fraction rho_n/rho. In an earlier paper (doi:10.5281/zenodo.23225293) we measured alpha(T) in the projected Gross-Pitaevskii equation at one cutoff and found alpha ~ 0.24 rho_n/rho; we read the coefficient as a mode-averaged transport cross-section of the vortex (1.4 healing lengths) and registered its cutoff dependence as the open question, with the Born phonon-scattering law (coefficient proportional to the cutoff) as the rival. Here the same energy-conserving instrument, with its T = 0 known answer reproduced to 0.3 % at each cutoff, is run at two further cutoffs (k_c xi = 2 pi/3 and 2 pi). The registered prediction of a cutoff-independent coefficient fails and so does its rival: the coefficient is 0.31 +- 0.03, 0.232 +- 0.014, 0.102 +- 0.011 at k_c xi = 2.1, 3.1, 6.3. What does not move is alpha itself: alpha/T = 0.054 +- 0.003 over six arms and three cutoffs (chi2 = 1.3 for 5 degrees of freedom; 0.060 +- 0.002 with the regression estimator, because the energy estimator, the primary one, is biased low by diffusion -- a flaw of the registered synthetic gate exposed by a second, Rust implementation of the estimators). The proportionality to rho_n was a degeneracy of a single cutoff, where rho_n is proportional to T; the 'geometric core size' reading is withdrawn. The temperature law was found after the data were read; it is registered as a hypothesis, with its prediction for the next arm fixed before that arm was read (reported in the paper). Pre-registrations, controls, scripts, analysis file and Lean companion (FrictionKinetic.lean) are in the software record.
Recent experiments at the University of Arizona have demonstrated superb control and imaging of flattened superfluids revealing a rich landscape of vortex dynamics. Vortex dipoles have been created and their dynamics observed [Neely et al. Phys. Rev. Lett. 104, 160410 (2010)], and highly charged persistent currents have been formed through a clever combination of stirring and dissipation [T. W. Neely, PhD thesis, (2010)]. Many ultra-cold gas systems, including these, are in a regime that we may characterize as quasi-equilibrium: despite some external forcing, a large fraction of the gas remains in thermal equilibrium. In this thesis we use the stochastic projected Gross-Pitaevskii equation (SPGPE) to make quantitative predictions of the quasi-equilibrium Bose gas. \n\nWe quantitatively compare the SPGPE with experiment by modeling the persistent current formation experiment performed at the University of Arizona [T. W. Neely, PhD thesis, (2010)]. We determine all SPGPE parameters for the toroidal system prior to simulation, enabling quantitative modeling of the experiment with no fitting parameters, giving a true test of the SPGPE theory. We find the SPGPE gives quantitative agreement with experiment, accurately predicting the size of the persistent current as well as the decay time of the vortices. We also observe the crucial role that thermal fluctuations have on enabling this agreement, showing the comprehensive SPGPE treatment is necessary to make quantitatively accurate calculations of quasi-equilibrium systems. This is the first quantitative agreement with experimental observations of vortex dynamics obtained with a theory of dissipation from first principles using no fitting parameters.\n\nWe then systematically study the Kelvin mode excitations on a vortex line in a three dimensional Bose-Einstein condensate using the SPGPE. We give a quantitative measure of the magnitude of vortex bending caused by the activation of Kelvin modes, and find that vortex bending can be suppressed by tightening the confinement along the direction of the vortex line. This leads to a strong suppression of the vortex decay rate as the system enters a regime of two-dimensional vortex dynamics, characterized by a critical oblateness.\n\nIn our final application of the SPGPE we simulate the decay of a vortex dipole. We model a recent experiment of Neely et al. [Neely et al. Phys. Rev. Lett. 104, 160410 (2010)], finding the SPGPE predicts a dipole lifetime consistent with experimental observations. We also show the experiment lies in the two-dimensional regime of vortex dynamics, validating our critical oblateness calculation.
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.
We have performed numerical simulations of various vortex configurations in a trapped BoseEinstein condensate by solving the two-dimensional Gross-Pitaevskii equation in the presence of a simple model of interaction between the condensate and the finite temperature thermal cloud that surrounds it. In that interaction the non-condensed thermal cloud acts as a source of dissipation with a damping effect of excitations. In the case ora single vortex and a vortex - anti vortex pair, we have found that the path of the vortices depends on the initial position, the initial separation distance if the case of two vortices and dissipation. This motion is periodic and it was found that sound waves are created by vortex motion; the intensity was stronger when the initial vortex separation distance was smaller. We have calculated the sound energy as the difference between the kinetic energy and the vortex energy. With no dissipation the vortices followed the same path with a slight oscillation due to the sound waves. We found that the smaller the initial vortex separation distance do is, the larger the sound production. The period, frequency, translation speed, sound energy and vortex energy were measured for different initial separation distances do and for different dissipation parameters 'Y. In the case of motion of one vortex, the connection between the dissipation 'Y and the friction coefficients, a and a' was studied as well. To create a simple turbulent state, we put eight pairs of vortex - anti vortex at random positions in the condensate with initial separation distance do =1.8 between them. 'Ve have studied the decay rate of the total energy, kinetic energy, quantum energy, trap energy and the z - component of the angular momentum together with the increase rate of the internal energy. Finally, we finished our investigation by putting randomly vortex - anti vortex pairs and studied the decrease of the number of vortices with time t. We found that the decrease is exponential.
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.
Dissipation of quantum vortex motion is fundamental to superfluid dynamics and quantum turbulence, yet there is currently a large gap between theory and experiments with ultracold atoms. Here we present a microscopic open quantum systems theory of thermally damped vortex motion in oblate atomic superfluids that includes previously neglected energy-damping interactions between superfluid and thermal atoms. This mechanism couples strongly to vortex core motion and causes dissipation of vortex energy due to mutual friction, as well as Brownian motion of vortices due to thermal fluctuations. We derive an analytic expression for the dimensionless mutual friction coefficient that gives excellent quantitative agreement with experimentally measured values, without any fitted parameters. Our work closes an existing two orders of magnitude gap between dissipation theory and experiments, previously bridged by fitted parameters, and provides a microscopic origin for the mutual friction and diffusion of quantized vortices in two-dimensional atomic superfluids.