Quantitative theory-experiment reconciliation for hybrid soliton-vortex-ring structures in the double light-roadblock Bose-Einstein condensate experiment

In a landmark 2005 experiment, Ginsberg, Brand, and Hau created pairs of density defects in a magnetically trapped sodium Bose-Einstein condensate (BEC) using 'slow-light' pulses halted at a double light roadblock, and observed compound excitations in which vortex rings are embedded in solitonic shells — hybrid soliton-vortex-ring structures. Their Gross-Pitaevskii simulations in imposed cylindrical symmetry reproduced the principal features but exhibited quantitative discrepancies with the data, and the detailed analysis of stationary spherical low-density shell states in traps was deferred to a companion publication. Later literature has resolved the deferred stationary-state question: the companion analysis of Carr and Clark (Phys. Rev. A 74, 043613 (2006)) and the Bogoliubov-de Gennes stability and decay analysis of Wang, Kevrekidis, Carretero-Gonzalez, and Frantzeskakis (Phys. Rev. A 93, 023630 (2016)) established the existence, stability, and vortex-mediated decay channels of spherical shell solitons in trapped condensates. What remains open is the experiment-specific core: a symmetry-unrestricted quantitative account of the double-roadblock experiment that explains the simulation-experiment discrepancies (accelerated evolution, reduced soliton curvature, attenuated collective-mode breathing, late-time shearing), establishes the role of laser-beam torque and broken cylindrical symmetry, and maps the experimental parameter window (density, trap frequencies, defect separation and depth) in which hybrid structures form, survive, and decay.

Background

A Bose-Einstein condensate (BEC) is a dilute gas of bosonic atoms cooled below the temperature at which a macroscopic fraction of atoms occupies a single quantum state; at mean-field level its order parameter $\psi(\mathbf{r},t)$ obeys the Gross-Pitaevskii equation (GPE), a nonlinear Schrödinger equation in which repulsive contact interactions produce a local sound speed and a healing length $\xi$ that sets the minimum length scale over which the density can change. A dark (or grey) soliton is a non-dispersive density depression on a planar wavefront whose amplitude fixes its velocity and the phase jump across it; in three dimensions such planar solitons are subject to the transverse 'snake' (Kadomtsev-Petviashvili) instability, which bends the wavefront and breaks it into vortex rings (VRs) — closed filaments of quantized circulation with depleted cores around which the phase winds by $2\pi$. BEC experiments in 2001 established both the creation of solitons from laser-imprinted or shock-like density defects and their snake-instability decay into VRs (Anderson et al., Phys. Rev. Lett. 86, 2926; Dutton et al., Science 293, 663). The Dutton-Hau experiment introduced the 'light roadblock': electromagnetically induced 'slow light' pulses, compressed by orders of magnitude in spatial extent, are brought to a halt inside the condensate where a thin wire masks the coupling laser beam, ejecting atoms in narrow density defects from which solitons are shed.

Ginsberg, Brand, and Hau (Phys. Rev. Lett. 94, 040403 (2005)) created two symmetrically placed defects with a double light roadblock in a cloud of $3.6\times10^{6}$ sodium atoms, so that the emitted solitons and vortex rings collide; the collision region early on displays a 'checkerboard' pattern of alternating solitonic density depressions. They observed a new class of excitations: hybrid structures in which VRs are embedded in roughly hemispherical solitonic shells, closed ellipsoidal solitonic shells incorporating four VRs of alternating circulation arranged 'like lines of constant latitude on a globe', and a single central low-density shell, with VR propulsion, attraction, and annihilation governing the dynamics. Generalized-GPE simulations (defect creation, in-trap evolution, and post-release expansion in three dimensions with imposed cylindrical symmetry) reproduced the principal observed features and identified the shells as soliton-VR compounds, but with acknowledged quantitative discrepancies: evolution slightly faster than the data beyond $\Delta t_{\rm trap}\approx1.3$ ms, reduced soliton curvature, and less dramatic collective-mode breathing (oscillation of the cloud's overall shape on its low-lying collective modes), attributed tentatively to the different initial defect separation used in the theory and to the cylindrical-symmetry approximation. The observed break-up of the closed shell at $\Delta t_{\rm trap}=6.6$ ms was noted as evidence of sensitivity to torque from the laser beams, an effect outside the axisymmetric simulations. The Letter also posed the question whether low-density shells can exist as stationary states, noting that the GPE admits spherical stationary solutions consisting of a near-uniform inner sphere surrounded by low-density shells resembling the spherical Bessel standing waves $\sin(kr)/r$ of the linear Schrödinger equation, that trapped versions with finite atom number were expected ('details will be given elsewhere'), and that the observed hybrid shells 'might play a role in the decay of the stationary solutions'.

That stationary-state question has since been answered in the literature. The deferred analysis appeared as the companion work the Letter itself cites: Carr and Clark, 'Vortices and ring solitons in Bose-Einstein condensates' (Phys. Rev. A 74, 043613 (2006)), which constructs the stationary radially symmetric states — ring solitons in two dimensions and spherical shells in three, with nonlinear-Bessel asymptotics — in free space and confined geometries, and performs a Bogoliubov-de Gennes (BdG) linear stability analysis — BdG theory linearizes the Gross-Pitaevskii equation about a solution to obtain its collective excitation spectrum and mode growth rates, and is the standard tool for deciding linear stability — showing that ring solitons in a harmonic trap are formally unstable but can have instability times far exceeding experimental timescales. Herring, Carr, Carretero-Gonzalez, Kevrekidis, and Frantzeskakis (Phys. Rev. A 77, 023625 (2008)) systematically continued the two-dimensional radially symmetric trapped states and their stability; Theocharis et al. (Phys. Rev. Lett. 90, 120403 (2003)) had earlier shown that two-dimensional ring dark solitons decay via the snake instability into vortex necklaces, a scenario observed experimentally in 2023 by Tamura, Chen, and Hung (Phys. Rev. X 13, 031029) in a box-trapped superfluid. Komineas and Brand (Phys. Rev. Lett. 95, 110401 (2005)) showed that inelastic soliton-vortex-ring collisions in cylindrical traps produce spherical shell-like structures related to the observed ones. Most directly, Wang, Kevrekidis, Carretero-Gonzalez, and Frantzeskakis (Phys. Rev. A 93, 023630 (2016)) analyzed spherical shell dark solitons in three-dimensional isotropically trapped condensates: their Bogoliubov-de Gennes spectra, their stability sufficiently close to the linear limit, and their instabilities leading to 'vortex line and vortex ring cages' — i.e., the vortex-mediated decay channels the 2005 Letter had conjectured. Mateo and Brand (New J. Phys. 17, 125013 (2015)) mapped the stability and dispersion of the three-dimensional solitary-wave families (solitonic vortices, vortex rings, intersecting-line 'Chladni solitons') in cylindrical traps, and the review of Frantzeskakis (J. Phys. A 43, 213001 (2010)) summarizes this state of the field.

By contrast, the experiment-specific items acknowledged in the 2005 Letter have no dedicated published treatment. Later experimental progress on soliton-vortex compounds used different platforms — density-engineered oscillating soliton/vortex rings (Shomroni et al., Nat. Phys. 5, 193 (2009)), vortex-ring phase imprinting (Bai et al., Results in Physics 22, 103828 (2021)), and box-quench ring dark solitons (Tamura et al., 2023) — and later simulations of planar-defect decay without imposed symmetry (Scherpelz et al., Phys. Rev. A 91, 033621 (2015)) and of the spontaneous symmetry-breaking instability of interacting vortex rings (Ikuta, Sugano, and Saito, Phys. Rev. A 99, 043610 (2019)) addressed generic rather than roadblock-specific dynamics. The quantitative reconciliation of the original experiment — including the laser-torque symmetry breaking and the parameter map of hybrid-structure formation — therefore remains an unfinished objective.

Problem Statement

Achieve a complete quantitative theoretical account of the hybrid soliton-vortex-ring structures observed in the double light-roadblock experiment of Ginsberg, Brand, and Hau. The stationary-state part of the original question — the existence, form, stability, and vortex-mediated decay channels of spherical low-density shell solutions in traps — has been resolved in later literature (Carr and Clark 2006; Wang et al. 2016), and the surviving core consists of three linked deliverables, all referring to the experiment as performed (magnetically trapped sodium condensate, trap frequencies $\omega_z=2\pi\times21$ Hz and $\omega_x=\omega_y=3\omega_z$, approximately $3.6\times10^{6}$ atoms, defects created by slow-light pulses halted at a wire-masked double light roadblock, in-trap evolution followed by approximately 20 ms of free expansion (the experimental value is 19.9 ms) before slice imaging, in which only a thin central slab, 20–30 micrometers thick, of the expanded cloud is imaged):

1. Symmetry-unrestricted simulation-experiment reconciliation. Construct a three-dimensional Gross-Pitaevskii-level description that does not impose cylindrical symmetry, including a model of the defect-creation step that represents the two-component slow-light coupling, the wire-masked coupling beam, and the probe-beam intensity profiles, and determine whether it quantitatively reproduces the observed evolution. In particular, explain the acknowledged discrepancies between the axisymmetric simulations and the data — evolution slightly faster than observed beyond $\Delta t_{\rm trap}\approx1.3$ ms, reduced soliton curvature, and less pronounced collective-mode breathing — and identify which physical effects (initial defect separation and shape, imperfect experimental symmetry, laser-beam torque) cause each discrepancy.

2. Symmetry breaking and torque. Establish the role of broken cylindrical symmetry and of external torque from the laser beams in producing the experimentally observed shearing and break-up of the closed low-density shell at $\Delta t_{\rm trap}=6.6$ ms: does the shearing require the laser-induced torque, or does it arise from an intrinsic instability of the soliton-vortex-ring compound?

3. Parameter-space map. Delineate, as a function of the experiment's control parameters (condensate density, trap frequencies, initial defect separation, and defect depth), which may be varied independently and should be explored over ranges that bracket their experimentally realized values, the regimes in which hybrid soliton-vortex-ring structures form, their lifetimes, and their decay pathways — replacing the Letter's qualitative remark that its parameters 'cover a regime in which VRs are produced but parent solitons are not destroyed' with a systematic and explained map.

All three deliverables are jointly required. The observed structures (checkerboard soliton patterns, closed shells with embedded vortex rings of alternating circulation, the shearing at late times, and the collective-mode breathing) should be quantitatively accounted for by the theory — the soliton positions and curvature, the shell geometry and the times at which the structures form and break up, and the collective-mode breathing amplitude reproduced to within the resolution and uncertainty of the reported experimental data — the sources of the simulation-data discrepancies identified, and the parameter regimes of hybrid-structure formation, lifetime, and decay delineated with their mechanisms explained.

References

1. N. S. Ginsberg, J. Brand, and L. V. Hau, 'Observation of Hybrid Soliton Vortex-Ring Structures in Bose-Einstein Condensates', Phys. Rev. Lett. 94, 040403 (2005), DOI 10.1103/PhysRevLett.94.040403, https://doi.org/10.1103/PhysRevLett.94.040403 (preprint arXiv:cond-mat/0408464, https://arxiv.org/abs/cond-mat/0408464)
2. B. P. Anderson, P. C. Haljan, C. A. Regal, D. L. Feder, L. A. Collins, C. W. Clark, and E. A. Cornell, 'Watching Dark Solitons Decay into Vortex Rings in a Bose-Einstein Condensate', Phys. Rev. Lett. 86, 2926 (2001), DOI 10.1103/PhysRevLett.86.2926, https://doi.org/10.1103/PhysRevLett.86.2926
3. Z. Dutton, M. Budde, C. Slowe, and L. V. Hau, 'Observation of Quantum Shock Waves Created with Ultra-Compressed Slow Light Pulses in a Bose-Einstein Condensate', Science 293, 663 (2001), DOI 10.1126/science.1062527, https://doi.org/10.1126/science.1062527
4. L. D. Carr and C. W. Clark, 'Vortices and ring solitons in Bose-Einstein condensates', Phys. Rev. A 74, 043613 (2006), DOI 10.1103/PhysRevA.74.043613, https://doi.org/10.1103/PhysRevA.74.043613 (preprint arXiv:cond-mat/0408460, https://arxiv.org/abs/cond-mat/0408460)
5. S. Komineas and J. Brand, 'Collisions of Solitons and Vortex Rings in Cylindrical Bose-Einstein Condensates', Phys. Rev. Lett. 95, 110401 (2005), DOI 10.1103/PhysRevLett.95.110401, https://doi.org/10.1103/PhysRevLett.95.110401 (preprint arXiv:cond-mat/0504072, https://arxiv.org/abs/cond-mat/0504072)
6. G. Theocharis, D. J. Frantzeskakis, P. G. Kevrekidis, B. A. Malomed, and Yu. S. Kivshar, 'Ring Dark Solitons and Vortex Necklaces in Bose-Einstein Condensates', Phys. Rev. Lett. 90, 120403 (2003), DOI 10.1103/PhysRevLett.90.120403, https://doi.org/10.1103/PhysRevLett.90.120403
7. D. J. Frantzeskakis, 'Dark solitons in atomic Bose-Einstein condensates: from theory to experiments', J. Phys. A: Math. Theor. 43, 213001 (2010), DOI 10.1088/1751-8113/43/21/213001, https://doi.org/10.1088/1751-8113/43/21/213001 (preprint arXiv:1004.4071, https://arxiv.org/abs/1004.4071)
8. G. Herring, L. D. Carr, R. Carretero-Gonzalez, P. G. Kevrekidis, and D. J. Frantzeskakis, 'Radially symmetric nonlinear states of harmonically trapped Bose-Einstein condensates', Phys. Rev. A 77, 023625 (2008), DOI 10.1103/PhysRevA.77.023625, https://doi.org/10.1103/PhysRevA.77.023625 (preprint arXiv:0709.2193, https://arxiv.org/abs/0709.2193)
9. A. Munoz Mateo and J. Brand, 'Stability and dispersion relations of three-dimensional solitary waves in trapped Bose-Einstein condensates', New J. Phys. 17, 125013 (2015), DOI 10.1088/1367-2630/17/12/125013, https://doi.org/10.1088/1367-2630/17/12/125013 (preprint arXiv:1510.01465, https://arxiv.org/abs/1510.01465)
10. W. Wang, P. G. Kevrekidis, R. Carretero-Gonzalez, and D. J. Frantzeskakis, 'Dark spherical shell solitons in three-dimensional Bose-Einstein condensates: Existence, stability and dynamics', Phys. Rev. A 93, 023630 (2016), DOI 10.1103/PhysRevA.93.023630, https://doi.org/10.1103/PhysRevA.93.023630 (preprint arXiv:1601.02176, https://arxiv.org/abs/1601.02176)
11. I. Shomroni, E. Lahoud, S. Levy, and J. Steinhauer, 'Evidence for an oscillating soliton/vortex ring by density engineering of a Bose-Einstein condensate', Nat. Phys. 5, 193–197 (2009), DOI 10.1038/nphys1177, https://doi.org/10.1038/nphys1177
12. P. Scherpelz, K. Padavic, A. Murray, A. Glatz, I. S. Aranson, and K. Levin, 'Generic equilibration dynamics of planar defects in trapped atomic superfluids', Phys. Rev. A 91, 033621 (2015), DOI 10.1103/PhysRevA.91.033621, https://doi.org/10.1103/PhysRevA.91.033621 (preprint arXiv:1410.0067, https://arxiv.org/abs/1410.0067)
13. M. Ikuta, Y. Sugano, and H. Saito, 'Symmetry-breaking instability of leapfrogging vortex rings in a Bose-Einstein condensate', Phys. Rev. A 99, 043610 (2019), DOI 10.1103/PhysRevA.99.043610, https://doi.org/10.1103/PhysRevA.99.043610
14. W.-K. Bai, J.-C. Xing, T. Yang, W.-L. Yang, and W.-M. Liu, 'Nonlinear dynamics of a Bose-Einstein condensate excited by a vortex ring phase imprinting', Results in Physics 22, 103828 (2021), DOI 10.1016/j.rinp.2021.103828, https://doi.org/10.1016/j.rinp.2021.103828
15. H. Tamura, C.-A. Chen, and C.-L. Hung, 'Observation of Self-Patterned Defect Formation in Atomic Superfluids - from Ring Dark Solitons to Vortex Dipole Necklaces', Phys. Rev. X 13, 031029 (2023), DOI 10.1103/PhysRevX.13.031029, https://doi.org/10.1103/PhysRevX.13.031029 (preprint arXiv:2211.08575, https://arxiv.org/abs/2211.08575)