Thermal behavior of Bose-Einstein condensates of polar molecules

Basic Information

• Paper ID: 2502.08339
• Title: Thermal behavior of Bose-Einstein condensates of polar molecules
• Authors: J. Sánchez-Baena, G. Pascual, R. Bombín, F. Mazzanti, J. Boronat
• Classification: cond-mat.quant-gas
• Publication Date: October 17, 2025
• Paper Link: https://arxiv.org/abs/2502.08339

Abstract

This paper employs the finite-temperature extended Gross-Pitaevskii equation (TeGPE) to investigate the condensate of dipolar NaCs molecules under the conditions of the recent breakthrough experiment Bigagli et al., Nature 631, 289 (2024). The study reports the condensate fraction of the system and the density distribution following time-of-flight expansion at the coldest experimental conditions, finding excellent agreement with experimental measurements. Additionally, the peak density of the ground state is reported and compared with experimental estimates. The results based on the TeGPE formalism successfully describe Bose-Einstein condensation of polar molecules at finite temperature.

Research Background and Motivation

Problem Background

1. Historic Breakthrough: Recent progress in controlling and reducing reactive losses in ultracold dipolar molecules has ultimately enabled the first Bose-Einstein condensation (BEC) of dipolar molecules, representing a major breakthrough in the field.
2. Theoretical Challenges: Although the condensate remains in a dilute regime, it represents the first step toward realizing highly correlated dipole-dominated molecular systems. The experiment achieved only approximately 60% maximum condensate fraction, with significant thermal dissipation effects present. Existing zero-temperature extended Gross-Pitaevskii equations cannot describe this situation.
3. Method Validation Requirements: There is a need to validate the effectiveness of Bogoliubov theory and the finite-temperature extended Gross-Pitaevskii equation (TeGPE) for theoretical descriptions of these novel systems.

Research Significance

• Advantages of Dipolar Molecules: Dipolar molecules possess dipole strengths 10 to 1000 times larger than magnetic atoms (e.g., 164Dy), allowing substantial increases in scattering length and gas parameter without losing dipolar physics
• Temperature Effects: Temperature is a critical parameter in theoretical descriptions of molecular condensates, with thermal effects playing counterintuitive roles in dipolar systems
• Theoretical Benchmarking: Provides benchmark testing for theoretical tools for dipolar systems at finite temperature

Core Contributions

1. First Theoretical Description: Successfully describes the finite-temperature behavior of experimentally realized dipolar molecular BEC using TeGPE
2. Accurate Predictions: Calculations of condensate fraction achieve excellent agreement with experimental values
3. Density Distribution Analysis: Accurately predicts density distribution following 17 ms time-of-flight expansion
4. Peak Density Estimation: Provides theoretical estimates of ground-state peak density and comparison with experiments
5. Method Validation: Establishes the validity of Bogoliubov theory and TeGPE in dilute dipolar molecular systems

Methodology Details

Task Definition

Investigate the finite-temperature Bose-Einstein condensation behavior of NaCs dipolar molecules in harmonic traps, including:

• Variation of condensate fraction with temperature
• Theoretical predictions of density distributions
• Density evolution following time-of-flight expansion

Model Architecture

1. Temperature-Dependent Extended Gross-Pitaevskii Equation (TeGPE)

The condensate wave function ψ(r) satisfies:

$\mu\psi(r) = \left(-\frac{\hbar^2\nabla^2}{2m} + U(r) + \int dr'V(r-r')\left[|\psi(r')|^2 + \rho_d(r')\right] + H_{qu}(r) + H_{th}(r)\right)\psi(r)$

where:

• μ is the chemical potential
• U(r) is the harmonic trap potential
• V(r) includes contact and dipole-dipole interactions
• ρ_d(r) is the non-condensed molecular density
• H_qu(r) and H_th(r) are quantum and thermal fluctuation corrections, respectively

2. Interaction Potential

The total interaction potential is: $V(r) = \frac{4\pi\hbar^2 a}{m}\delta(r) + V_{dd}(r)$

The dipole-dipole interaction: $V_{dd}(r) = \frac{d_{eff}^2}{4\pi\epsilon_0}\frac{1-3\cos^2\theta}{r^3}$

3. Quantum Fluctuation Correction (Lee-Huang-Yang Term)

$H_{qu}(r) = \frac{32}{3\sqrt{\pi}}g\sqrt{a^3Q_5(a_{dd}/a)}|\psi(r)|^3$

4. Thermal Fluctuation Correction

$H_{th}(r) = \int\frac{dk}{(2\pi)^3}\frac{V(k)}{e^{\beta\varepsilon_k}-1}\left[|u_k|^2 + |v_k|^2 - 2|u_k||v_k|\right]$

Bogoliubov Amplitude Calculation

Using local density approximation (LDA) and improved Bogoliubov-de Gennes equations:

$E_k u_k = \left(\frac{\hbar^2k^2}{2m} + \tilde{F}(\rho_0)\right)u_k + \rho_0 u_k V(k) - \rho_0 v_k V(k)$

$-E_k v_k = \left(\frac{\hbar^2k^2}{2m} + \tilde{F}(\rho_0)\right)v_k + \rho_0 v_k V(k) - \rho_0 u_k V(k)$

where $\tilde{F}(\rho_0) = \text{Max}\{F(\rho_0), 0\}$ accounts for curvature effects of the condensate at trap boundaries.

Technical Innovations

1. Improved Thermal Correction Calculation: Resolves divergence issues of traditional LDA methods at trap boundaries by considering condensate curvature effects
2. Density Functional Approach: Extends thermal density results from uniform systems to non-uniform trap systems
3. Time Evolution Simulation: Uses split-operator methods to simulate time-of-flight expansion

Experimental Setup

System Parameters

• Molecular Species: NaCs dipolar molecules
• Dipole Length: a_dd = 1250a_0
• s-wave Scattering Length: a = 1500a_0
• Dipole Parameter: ε_dd = a_dd/a = 0.83 < 1 (contact interaction dominated)
• Trap Frequencies: (ω_x, ω_y, ω_z) = 2π(23, 49, 58) Hz
• Molecular Number: N = 289-884
• Temperature Range: T = 6-12.6 nK

Comparison Methods

• Experimentally measured condensate fraction (obtained through bimodal distribution fitting)
• Experimental time-of-flight expansion density distribution
• Experimentally estimated peak density

Implementation Details

• Imaginary time evolution to solve TeGPE ground state
• Split-operator method for real-time evolution
• Iterative solution of self-consistent condensate wave function and thermal density

Experimental Results

Main Results

1. Condensate Fraction Comparison

Theoretically calculated condensate fractions are in complete agreement with experimental measurements within error margins:

2. Density Distribution Prediction

At the coldest experimental point (N=289, T=6nK), the theoretical prediction of density distribution following 17 ms time-of-flight expansion achieves excellent agreement with experimental measurements.

3. Peak Density Analysis

• TeGPE-predicted total density peak is relatively close to experimental estimates
• Considering that experimental estimates neglect thermal cloud effects, the pure condensate peak density is in complete agreement with experimental results

Ablation Studies

1. Bimodal Fitting Verification

Using bimodal distribution fitting on theoretical density distributions:

• Fitted condensate fraction is slightly lower than direct TeGPE calculation
• Fitting quantitatively captures the width and center position of condensate and thermal cloud

2. Dipole Orientation Effects

Comparison of results for z-axis and x-axis dipole orientations:

• Density distributions show only minor differences
• Condensate fraction variation approximately 0.01, indicating insensitivity to dipole orientation

Key Findings

1. Theory Validity: TeGPE maintains high precision in contact-interaction-dominated dipolar systems
2. Universality: The system is in the universal regime where physics is primarily controlled by scattering length a
3. Importance of Thermal Effects: Theory remains accurate even when condensate fraction is as low as 20%
4. Critical Temperature: The critical temperature for the N=289 molecular system is approximately 9.8 nK

Related Work

Dipolar Quantum Gas Research

• Magnetic Atom Experiments: Primarily using Dy and Er atoms, but dipolar effects limited to low gas parameter regions
• Supersolidity: Observed in 164Dy atoms driven by dipolar interactions
• Droplet Formation: Self-bound droplets driven by dipolar interactions

Theoretical Method Development

• Zero-Temperature eGPE: Fails in strongly correlated regions, requiring quantum Monte Carlo methods
• Finite-Temperature Theory: Relatively limited applications in dipolar systems
• Bogoliubov Theory: Successful applications in weakly interacting regions

Molecular BEC Progress

• Reactive Loss Control: Development of techniques such as microwave shielding
• First Realization: Breakthrough experiment by Bigagli et al.
• Theoretical Description: First systematic theoretical analysis in this work

Conclusions and Discussion

Main Conclusions

1. Theory Validation: TeGPE successfully describes the finite-temperature behavior of the first dipolar molecular BEC
2. Accurate Predictions: Achieves excellent agreement with experiments in condensate fraction, density distribution, and peak density
3. Method Reliability: Establishes the validity of Bogoliubov theory in dilute dipolar molecular systems
4. Physical Insight: Demonstrates that the system is in the universal regime, primarily controlled by scattering length

Limitations

1. Applicable Range: Valid only in dilute regions (ρa³ < 7.2×10⁻⁴)
2. Strongly Correlated Limit: Theory fails in dipole-dominated strongly correlated regions
3. Interaction Details: Uses simplified pseudopotential approximation
4. Time Evolution: Assumes thermal cloud remains unchanged during time-of-flight expansion

Future Directions

1. Strong Dipole Regions: Study molecular condensates in dipole-interaction-dominated regimes
2. Droplet Formation: Expect critical molecular droplet numbers to be 10 times smaller than atoms
3. Supersolidity: Explore whether phase coherence can be maintained between molecular clusters
4. Thermal Effects: Investigate whether heating can induce supersolidity in molecular systems
5. Quantum Monte Carlo: Use exact methods to further improve theoretical descriptions

In-Depth Evaluation

Strengths

1. Timeliness: Follows experimental breakthroughs closely, providing the first systematic theoretical analysis
2. Accuracy: Impressive agreement between theoretical predictions and experimental results
3. Completeness: Covers multiple aspects including condensate fraction, density distribution, and time-of-flight evolution
4. Methodological Innovation: Improves traditional LDA methods, resolving divergence issues
5. Physical Insight: Reveals universal properties of the system and applicable theoretical ranges

Weaknesses

1. Theoretical Limitations: Acknowledges theory failure in strongly correlated regions, requiring more precise methods
2. Simplified Assumptions: Uses pseudopotential approximation and static thermal cloud assumption
3. Parameter Dependence: Some experimental parameters (e.g., dipole orientation) not fully determined
4. Limited Scope: Applicable only to weakly dipolar interaction regions

Impact

1. Academic Value: Establishes theoretical foundation for emerging molecular BEC field
2. Practical Value: Provides theoretical guidance for subsequent experimental design
3. Methodological Contribution: Improved TeGPE methods applicable to other dipolar systems
4. Field Advancement: Lays foundation for exploring strongly correlated molecular physics

Applicable Scenarios

1. Dilute Molecular BEC: Applicable to contact-interaction-dominated dipolar molecular systems
2. Finite-Temperature Analysis: Particularly suitable for systems with significant thermal effects
3. Experimental Comparison: Provides theoretical benchmarks for molecular BEC experiments
4. Parameter Optimization: Guides selection and optimization of experimental parameters

References

1. Bigagli et al., "Observation of bose–einstein condensation of dipolar molecules," Nature 631, 289–293 (2024)
2. Sánchez-Baena et al., "Heating a dipolar quantum fluid into a solid," Nature Communications 14, 1868 (2023)
3. Langen et al., "Dipolar droplets of strongly interacting molecules," Phys. Rev. Lett. 134, 053001 (2025)
4. Karman et al., "Double microwave shielding," arXiv:2501.08095 (2025)

This paper represents important progress in molecular BEC theoretical research, successfully applying mature finite-temperature theoretical methods to this emerging field and establishing a solid theoretical foundation for future exploration of more complex molecular quantum many-body phenomena.