Symmetry and Self-Bound Droplets in Dipolar Molecular Gases

Basic Information

• Paper ID: 2510.04634
• Title: Symmetry and Self-Bound Droplets in Dipolar Molecular Gases
• Author: D. Baillie (Department of Physics, University of Otago, New Zealand)
• Classification: cond-mat.quant-gas (Condensed Matter Physics - Quantum Gases)
• Publication Date: October 13, 2025
• Paper Link: https://arxiv.org/abs/2510.04634

Abstract

Recent experiments with degenerate molecular gases dressed by elliptically polarized microwave fields have enabled new control of dipolar interactions via engineered anisotropy. We reveal a symmetry structure of the dipolar interaction that generates degeneracies among the interaction parameters, enabling a classification of spatial symmetries and equilibrium shapes of the gases. Exploiting these symmetries, we analyze solutions including beyond-meanfield quantum fluctuations, and develop a complementary variational theory. We map out the phase diagram of self-bound droplets and characterize their widths, energies, and densities.

Research Background and Motivation

Problem Background

1. Development of Ultracold Molecular Gases: Ultracold diatomic alkaline-earth molecular gases have long attracted attention due to their strong long-range anisotropic dipolar interactions. Early fermionic studies were hindered by chemical reaction losses until the realization of Bose molecular condensates in 2024.
2. Discovery of Self-Bound Droplets: Self-bound dipolar droplets observed in magnetic atoms represent stable structures in free space without external confinement, stabilized by quantum fluctuations. Recently, self-bound droplets have also been reported in molecular gases.
3. New Progress in Microwave Field Control: Experiments with elliptically polarized microwave fields have enabled control of molecular gas interactions beyond the cylindrical symmetry of atomic dipolar interactions, achieving higher anisotropy.

Research Motivation

1. Theoretical Gap: Existing theories of molecular Bose gases are primarily limited to cylindrically symmetric interactions, with very few works considering elliptically engineered anisotropy.
2. Insufficient Understanding of Symmetry: Lack of systematic understanding of the symmetry structure of microwave-modified dipolar interactions.
3. Missing Phase Diagram: Complete phase diagram of self-bound droplets and their properties are needed.

Core Contributions

1. Revelation of Symmetry Structure: Discovery of the D₃ symmetry group structure of microwave-modified dipolar interactions, with each interaction parameter combination belonging to a six-fold degenerate state.
2. Establishment of Classification Framework: Classification of spatial symmetries and equilibrium shapes of gases based on symmetry, determining conditions for preservation or breaking of cylindrical symmetry.
3. Development of Theoretical Methods: Application of symmetry to the extended Gross-Pitaevskii equation (eGPE), incorporating beyond-meanfield quantum fluctuation effects.
4. Construction of Variational Theory: Development of variational solutions capturing key anisotropic effects, providing analytical expressions for interaction energies.
5. Complete Phase Diagram Mapping: Phase diagram of self-bound droplets in both thermodynamic limit and finite particle number regimes, characterizing their widths, energies, and densities.

Detailed Methodology

Task Definition

Investigation of formation conditions, symmetry properties, and physical characteristics of self-bound droplets in dipolar molecular gases modified by elliptically polarized microwave fields, establishing a complete theoretical framework.

Interaction Model

The effective dipole-dipole potential of microwave-shielded molecules is: $U(r) = \frac{3g_s}{4\pi r^3}[\epsilon_0(1-3\cos^2\theta) + \sqrt{3}\epsilon_2\sin^2\theta\cos 2\varphi]$

where:

• $g_s = 4\pi\hbar^2 a_s/m$ is the s-wave interaction strength
• $\epsilon_0$ is the relative dipolar interaction strength
• $\epsilon_2 \propto \sin 2\xi$ is related to microwave ellipticity
• $\theta$ is the angle between $\mathbf{r}$ and the z-axis, $\varphi$ is the azimuthal angle

Symmetry Analysis

The interaction is rewritten as: $U(r) = \frac{3g_s}{4\pi r^5}[(\epsilon_0+\sqrt{3}\epsilon_2)x^2 + (\epsilon_0-\sqrt{3}\epsilon_2)y^2 - 2\epsilon_0 z^2]$

This reveals key symmetries:

1. Six-fold Degeneracy: Each state belongs to a six-fold multiplet related through coordinate permutations
2. D₃ Group Structure: Parameter transformation $\epsilon' = T\epsilon$, where T is an element of the dihedral group D₃
3. Symmetry Axes: $\epsilon_2=0$, $\epsilon_0=\epsilon_2/\sqrt{3}$, $\epsilon_0=-\epsilon_2/\sqrt{3}$

Extended Gross-Pitaevskii Equation

$\left[-\frac{\hbar^2\nabla^2}{2m} + V(\mathbf{x}) + \Phi(\mathbf{x}) + \gamma_{QF}|\psi|^3\right]\psi = \mu\psi$

where the quantum fluctuation coefficient is: $\gamma_{QF} = \frac{32g_s a_s^{3/2}Q_5(\epsilon_0,\epsilon_2)}{3\sqrt{\pi}}$

$Q_5(\epsilon_0,\epsilon_2) = \int\frac{d\Omega_k}{4\pi}[1+\bar{U}(k)]^{5/2}$

Variational Method

Employing a Gaussian variational wavefunction: $\psi_v(\mathbf{x}) = \sqrt{\frac{N}{\pi^{3/2}l_x l_y l_z}} \exp\left(-\frac{1}{2}\left(\frac{x^2}{l_x^2}+\frac{y^2}{l_y^2}+\frac{z^2}{l_z^2}\right)\right)$

The interaction energy is expressed as a combination of anisotropic functions $f(l_x/l_z, l_y/l_z)$ and $f_2(l_x/l_z, l_y/l_z)$.

Experimental Setup

Theoretical Calculation Parameters

• Particle Number: $N = 1000, 5000$
• Interaction Parameter Range: $-3 \leq \epsilon_0 \leq 4$, $0 \leq \epsilon_2 \leq 3$
• Length Unit: s-wave scattering length $a_s$
• Energy Unit: $E_0 = \hbar^2/(ma_s^2)$

Numerical Methods

1. Grid Method: Non-uniform grids adapted to density distribution
2. Momentum Space Transform: Zero-padded grids for spherical truncation
3. Numerical Integration: Numerical integration for quantum fluctuation coefficients

Comparative Analysis

Focus on three typical cases:

• Prolate symmetry: $(\epsilon_0 > 0, \epsilon_2 = 0)$
• Oblate symmetry: $(\epsilon_0 < 0, \epsilon_2 = 0)$
• Asymmetric: $\epsilon_0 = 0$

Experimental Results

Phase Diagram Structure

1. Thermodynamic Limit Boundary: Self-bound droplet existence condition is $\epsilon_0 + \sqrt{3}|\epsilon_2| > 1$ or $\epsilon_0 < -1/2$
2. Finite Particle Number Effects: Phase boundaries approach the thermodynamic limit as N increases. eGPE-computed boundaries are closer to the thermodynamic limit than variational Gaussian predictions.
3. Significant Energy Differences: Energy in the $(\epsilon_0 < 0, \epsilon_2 = 0)$ and $\epsilon_0 = 0$ cases is notably lower than in the $(\epsilon_0 > 0, \epsilon_2 = 0)$ case.

Large N Limit Analytical Results

Energy expression: $\frac{E}{NE_0} = -\frac{50\pi^2(\epsilon_{lim}-1)^3}{3[64Q_5(\epsilon_0,\epsilon_2)]^2}$

Peak density: $n_{peak}a_s^3 = \frac{25\pi(\epsilon_{lim}-1)^2}{[64Q_5(\epsilon_0,\epsilon_2)]^2}$

where $\epsilon_{lim} = \max_k(-\bar{U}(k))$.

Droplet Property Characterization

1. Shape Features:
   • Prolate: $\epsilon_0 = 3, \epsilon_2 = 0$ shows prolate contours and liquid-like plateau density
   • Oblate: $\epsilon_0 = -1.5, \epsilon_2 = 0$ exhibits oblate contours
   • Asymmetric: $\epsilon_0 = 0$ shows different contours along three axes
2. Width Evolution: From N=1000 to N=5000, long widths increase significantly, while short widths increase marginally or sometimes decrease
3. Density and Energy: Peak density rises sharply with increasing $|\epsilon_0|$, but decreases for $|\epsilon_0| > 2$ due to quantum fluctuations

Key Findings

1. Constant Aspect Ratio: For the $\epsilon_0 = 0$ case, in the $l_x \ll l_z$ limit, $\lambda = l_y/l_z \approx 1.964$ is constant
2. Quantum Fluctuation Dominance: Near phase boundaries, quantum fluctuation terms dominate energy and density behavior
3. Symmetry Verification: Numerical results fully verify theoretically predicted symmetry relations

Related Work

Theoretical Development of Dipolar Gases

Existing theories focus on:

1. Atomic Magnetic Gases: Extensive studies of cylindrically symmetric dipolar interactions
2. Molecular Condensates: Limited theoretical work on cylindrical symmetry
3. Supersolid Phases: Recent studies on elliptical anisotropy

Experimental Progress Comparison

1. Atomic Systems: Lifetime limitations of rapidly rotating magnetic fields
2. Molecular Systems: Long lifetimes and strong interactions achieved through microwave shielding
3. Anisotropy Control: Advantages of molecular systems beyond cylindrical symmetry

Conclusions and Discussion

Main Conclusions

1. Symmetry Framework Established: Successfully revealed the D₃ symmetry structure of microwave-modified dipolar interactions, providing an organizing principle for understanding the interaction landscape
2. Complete Phase Diagram Mapped: Comprehensive characterization of self-bound droplet existence regions in both thermodynamic limit and finite particle number regimes
3. Analytical Theory Developed: Established widely applicable variational theory providing analytical expressions for interaction energies

Limitations

1. Local Density Approximation: The local density approximation for quantum fluctuations may have limited accuracy under strong anisotropy
2. Imaginary Part Treatment: Simple discarding of imaginary parts of quantum fluctuations requires more refined theory
3. Finite Temperature Effects: Finite temperature effects on droplet stability are not considered

Future Directions

1. Excitation Spectrum Research: Analysis of collective excitation spectra of self-bound molecular droplets
2. Experimental Verification: Direct verification of symmetry predictions through parameter space scanning of $(\epsilon_0, \epsilon_2)$
3. Dynamical Evolution: Investigation of dynamical processes in droplet formation and evolution

In-Depth Evaluation

Strengths

1. Strong Theoretical Innovation: First systematic revelation of symmetry structure in microwave-modified dipolar interactions, providing a new theoretical framework for the field
2. Complete Methodology: Combination of analytical theory, numerical computation, and variational methods provides a comprehensive theoretical toolkit
3. Deep Physical Insight: Symmetry analysis greatly simplifies the complex interaction landscape, making the physical picture clearer
4. High Practical Value: Provides clear theoretical guidance for experimental design, particularly symmetry predictions for parameter scanning

Weaknesses

1. Approximation Treatment: Local density approximation and imaginary part treatment of quantum fluctuations may affect accuracy in strong interaction regions
2. Missing Experimental Comparison: As purely theoretical work, lacks direct comparison with latest experimental results
3. Missing Dynamics: Focuses only on ground state properties, does not address dynamical stability and excited states of droplets

Impact

1. Field Contribution: Provides important symmetry framework for ultracold molecular gas theory, expected to become an important reference in the field
2. Practical Value: Directly guides ongoing molecular droplet experiments
3. Reproducibility: Clear theoretical derivations and detailed numerical method descriptions ensure good reproducibility

Applicable Scenarios

1. Ultracold Molecular Experiments: Directly applicable to current microwave-shielded molecular gas experiments
2. Theoretical Model Research: Symmetry framework can be generalized to other systems with similar anisotropy
3. Numerical Computation Optimization: Symmetry relations can significantly reduce computational load and improve efficiency

References

The paper cites important literature in the field, including:

• Pioneering work on ultracold molecular gases 1-12
• Theoretical and experimental development of dipolar droplets 13-25
• Latest experiments on molecular self-bound droplets 26
• Theoretical foundations of microwave shielding techniques 5-7, 40-41

This paper is of significant importance in theoretical physics, providing powerful theoretical tools for understanding and predicting novel quantum phases in ultracold molecular gases. Its symmetry analysis methods and complete theoretical framework are expected to have profound impacts on subsequent research in this field.