Superradiance from lattice-confined atoms inside hollow core fibre

Abstract

Unravelling superradiance, also known as superfluorescence, relies on an ensemble of phase-matched dipole oscillators and the suppression of inhomogeneous broadening. Here we report a superradiance platform that combines an optical lattice free from the ac Stark shift and a hollow-core photonic crystal fibre, enabling an extended atom-light interaction over 2 mm free from the Doppler effect. This system allows control of the atom spatial distribution and spectral homogeneity whilst efficiently coupling the radiation field to an optical fibre. The experimentally-observed and theoretically-corroborated temporal, spectral and spatial dynamic behaviours of the superradiance, e.g., superradiance ringing and density-dependent frequency shift, demonstrate a unique interplay between the trapped atoms and the fibre-guided field with multiple transverse modes. Our theory indicates that the resulting temporal evolution of the guided light shows a minimal beam radius of 3.1 µm which is three times smaller than that of the lowest-loss fibre mode.

Introduction

Superradiance (SR) has been the subject of active research since the pioneering work of Dicke in 1954¹. Such a collective effect arises from an ensemble of two-level emitters spontaneously locking their dipole oscillations in phase², giving rise to an enhanced decay rate γ_SR many times faster than the spontaneous emission rate γ₀ of individual emitters and a strong peak intensity of radiation in quadratic proportion to the excited-state population. The other SR signature manifests as a collective frequency shift and broadening, which are underpinned by fascinating physical phenomena such as the energy-level shift due to virtual photon exchange (collective Lamb shift)³ or van der Waals dephasing². In addition to its undisputed importance in fundamental science, SR also has substantial applications in various fields. For example, lasing based on SR has been proved to be a potentially outstanding frequency reference because of its strong immunity to environmental perturbations^4,5. In quantum information, SR is utilized in realizing quantum memories^6,7 and single photon sources⁸. Also, SR has been recently used for cooling mechanisms by means of cooperative emission^9,10.

Despite the observations of SR in various gaseous and solid-state samples^{11,12,13,14,15,16,17,18,19}, capturing the full picture of SR properties within a single experimental platform still remains an ongoing target. The chief challenge in unravelling the SR dynamics is to find the best trade-off between maximising the number of phase-locked emitters coupled to the common radiation field and preserving the inter-emitter phase correlation. The enhancement of SR emission rate scales as the former while the latter is a prerequisite for the SR generation. The conventional SR schemes adopt an ensemble of emitters densely packed within a sub-radiation-wavelength volume, i.e., Dicke regime^11,12,13,14, to achieve the SR. Recently, a large number of atoms confined inside optical cavities^{4,5,6,7,8,9,10} or coupled to waveguides¹⁹ have been a successful alternative for SR generation. This is based on the fact that the fraction of photon energy emitted into the SR mode Nη/(1 + Nη) approaches unity for the atomic number N ≫ 1 even when the single-atom cooperativity parameter η = γ_SR/γ₀, which is a geometric parameter characterizing the emissive coupling of an atom to the cavity mode, is much less than unity⁸.

Research at the interdisciplinary frontier between the fields of cold atoms and photonic materials reports promising results in controlling ultracold atoms with photonic waveguides²⁰, which holds the potential of offering a strong cooperative enhancement while keeping the low atomic density. Indeed, ultracold atoms are especially attractive materials since they possess the distinct feature of being precisely controllable by optical lattice²¹. In particular, when ultracold atoms are confined in an optical lattice tuned to the magic wavelength, the optical excitation of atomic transition is not only free from the Doppler shift but also free from the lattice-field-induced ac Stark shift²², offering an ideal platform for investigating SR based on well-isolated atomic systems. In addition, low-loss photonic waveguides allow a strongly enhanced one-dimensional (1D) atom-light interaction, as they allow an interaction length much longer than the Rayleigh range $z_{\mathrm{R}} = \pi w_0^2/\lambda _0$ at the wavelength λ₀ while still keeping a small light-beam waist of radius w₀. Successful combination of ultracold atoms and photonic waveguide has been exemplified by the works with nanofibers^23,24, photonic waveguides¹⁹, and hollow-core fibers²⁵. The long-range atom-atom interaction is accomplished via guided photons of the so-called alligator photonic crystal waveguide¹⁹, and the resulting single-atom cooperativity reaches η = 0.34. Ultracold atoms are successfully loaded inside a hollow-core photonic crystal fibre (HCPCF) and trapped in a 1D magic-wavelength optical lattice²⁵, revealing 10-kHz-wide linewidth of the transition without being affected by wall-atom collisions.

In the present work, we report on the observation of SR from ultracold atoms confined in the magic optical lattice inside the HCPCF²⁶. We couple an HCPCF-transmitted SR to a single-mode fibre to measure the signal by employing a highly sensitive heterodyne technique, which allows investigating temporal, spectral and spatial dynamic behaviours of SR with a single-shot measurement. We observe a ringing and accelerated decay that is up to 55 times faster than the spontaneous emission of single atoms. By analysing the inphase/quadrature components of the signal, we measure a density-dependent frequency shift. Finally, we discuss the multimode dynamics of SR in the HCPCF, where a superposition of specific guided modes contributes to the photon-mediated SR, complemented with theoretical predictions. The dense cigar-shaped atomic cloud restricts the guided SR field mainly within a diameter three times smaller than that of the fundamental fibre mode. For an in-fibre atom number of N = 2.5 × 10⁵ with the single-atom cooperativity η = 4.1 × 10^−5 available for the HCPCF, the N-atom cooperativity parameter⁸ Nη can be as high as 10, leading to a near-unity coupling efficiency to the guided mode of the hollow-core fibre.

Results

Experimental scheme

Figure 1a shows the schematic of our experiment. A 32-mm-long hypocycloid-core-shaped Kagome-lattice HCPCF²⁶ is horizontally placed in a vacuum chamber. The inner diameter of the hollow-core fibre reaches 40 μm. The optical fibre guides light via the inhibited-coupling mechanism²⁷ and supports multiple transverse modes whose transverse profiles are very close to those of a dielectric capillary²⁸. The Gaussian-like fundamental LP₀₁ core mode has a 1/e² waist radius w₀ = 11.8 μm (1/e radius is 8.34 μm). ⁸⁸Sr atoms are laser cooled and trapped in a two-stage magneto-optical trap (MOT)²⁹ at z ≈ −1 mm (the origin of z-coordinates is set at the entrance end of the HCPCF) and loaded into an optical lattice which passes through the HCPCF (see Methods). The counter-propagating linearly-polarized lattice lasers match the LP₀₁ mode and operate at the wavelength λ_L = 813 nm. A bias magnetic field ${\mathbf{B}}_0 = \left( {0.14\,{\mathrm{mT}}} \right){\hat{\mathbf{e}}}_x$ is applied to define the quantization axis. We focus on the ¹S₀ − ³P₁(m = 0) intercombination transition of ⁸⁸Sr at the frequency ω₀ = 2π × 435 THz with the natural linewidth γ₀ = 2π × 7.5 kHz and the dipole moment D = 0.087 ea₀ (elementary charge e and Bohr radius a₀). The lattice laser polarization is tuned to cancel out the ac Stark shift difference between ¹S₀ and ³P₁(m = 0) states caused by the optical lattice^21,25.

Fig. 1

Experimental setup. a Schematic of experimental setup. Ultracold ⁸⁸Sr atoms are loaded from a two-stage magneto-optical trap (MOT) into an optical lattice and then transported inside the hollow-core photonic crystal fibre (HCPCF) at z = 2 mm by the moving lattice. The magnetic field B₀ = (0.14mT) ${\hat{\mathbf{e}}}_x$ is applied to define the quantization axis of the system. The superradiance (SR) emission strongly couples to the guided modes of the fibre while the spontaneously emitted photons hardly match the fibre modes. The HCPCF guided SR is coupled to a single-mode fibre (SMF) that is mainly mode matched to the LP₀₁-mode component in SR. A part of the SR output is detected by a photomultiplier tube (PMT) to measure the total intensity of SR including the higher order modes. b Fluorescence images of atoms inside the HCPCF by exciting the ¹S₀ − ¹P₁ transition at 461 nm, which are taken by an electron multiplying charge-coupled device (EMCCD) camera through the fibre wall. The FWHM of the atomic distribution is l_z = 0.87 mm for the unexpanded cloud and l_z = 2.1 mm for the expanded one. c A pump laser at 689 nm with a linewidth less than 1 kHz is split into two paths. One is coupled to a SMF and used as a local-oscillator E_LO cos ω_LOt, while the other is used as the pump light with the amplitude E_p at the frequency ω_p = ω₀ + Δ_p shifted by an acousto-optic modulator (AOM). The SR field from the HCPCF is coupled to the SMF and balanced-heterodyned with the local field. The resultant signal $V_{{\mathrm{RF}}}\left( t \right) \propto {\mathrm{Re}}\left\{ {\left[ {\int} {{\int} {E_{{\mathrm{LO}}}\left( {\mathbf{r}} \right)} } {E_{{\mathrm{SR}}}^\prime \left( {{\mathbf{r}},t} \right)d^2{\mathbf{r}}} \right]e^{ - i{\mathrm{\Omega }}_0t}} \right\}$ at around Ω₀ = 2π × 50 MHz is band-pass-filtered (BPF) and recorded by a digital storage oscilloscope (DSO). d Frequency relations. In respect to the ¹S₀ − ³P₁(m = 0) transition frequency ω₀, the SR, pump-light, and local-oscillator frequencies are given by ω_SR = ω₀ + Δ_SR, ω_p = ω₀ + Δ_p, and ω_LO = ω₀ − Ω₀, respectively

The loaded atoms are further transported to the position z ≈ 2 mm inside the hollow-core fibre by a moving lattice with its depth of 300 μK²⁵. By turning off one of the two lattice lasers, the atomic gas, which is radially guided by a travelling-wave lattice laser, spreads axially to a full width at half maximum (FWHM) l_z. Then, the lattice laser is switched on again to recapture the atoms. We consider two cloud widths, i.e. unexpanded (l_z = 0.87 mm) and expanded (l_z = 2.1 mm) clouds (see Fig. 1b), corresponding to the number of lattice sites N_L = l_z/(λ_L/2) = 2.1 × 10³ and 5.2 × 10³, respectively (see Methods). Each lattice site forms a pancake-shaped harmonic potential, where the atoms are distributed within a transverse 1/e radius of r_a = 1.7 μm in the x–y plane and a longitudinal 1/e width of l_a = 54 nm along the fibre axis. Varying the number of loaded atoms N tunes the on-site density of atoms $\rho = N_{\mathrm{a}}/\left( {\pi r_{\mathrm{a}}^2l_{\mathrm{a}}} \right)$ with the average number of atoms in each lattice site N_a = N/N_L. We note that the multiple photon scattering caused by the periodic spatial distribution of atoms³⁰ hardly affects the atom-light interaction as the lattice period λ_L/2 does not fulfil the Bragg condition for the transition wavelength λ₀.

Superradiance ringing

All atoms are initialized in the ¹S₀ ground state. We apply an ${\hat{\mathbf{e}}}_x$ - polarized light pulse E_p with a duration τ_p = 500 ns to pump the atoms. The pump-laser beam is matched to the LP₀₁ fibre core mode, and the in-fibre amplitude E_p is set to E_p = ħπ/Dτ_p with ħ being the reduced Planck constant. The pump-laser frequency ω_p is stabilized to an optical cavity made of an ultralow-expansion (ULE) glass with a drift rate ~1 kHzh^−1 and detuned Δ_p = ω_p − ω₀ from the atomic resonance ω₀. The observed SR field at the fibre output end can be written as ${\mathbf{E}}_{{\mathrm{SR}}}\left( {{\mathbf{r}},t} \right) = \left( {{\hat{\mathbf{e}}}_x/2} \right)E_{{\mathrm{SR}}}\left( {{\mathbf{r}},t} \right)e^{ - i\omega _0t} + {\mathrm{c}}.{\mathrm{c}}$. with the radial coordinate ${\mathbf{r}} = x{\hat{\mathbf{e}}}_x + y{\hat{\mathbf{e}}}_y$. We observe that the SR field propagates in the same direction as the pump light, which is consistent with the directed spontaneous emission of the timed Dicke state^3,31. The radiation power $P_{{\mathrm{SR}}}\left( t \right) \propto {\int} {{\int} {|{\mathbf{E}}_{{\mathrm{SR}}}({\mathbf{r}},t)|^2} } d^2{\mathbf{r}}$ is typically on the order of nW, requiring the low noise detection. To address this issue, we utilize the balanced-heterodyne technique³² (see Fig. 1c). The SR output is coupled to a single-mode fibre (SMF) and mixed with a local oscillator ${\mathbf{E}}_{{\mathrm{LO}}}\left( {{\mathbf{r}},t} \right) = {\hat{\mathbf{e}}}_xE_{{\mathrm{LO}}}\left( {\mathbf{r}} \right){\mathrm{cos}}\omega _{{\mathrm{LO}}}t$ by a 50:50 fibre coupler (see Methods). The power and the frequency of local oscillator are on the order of 1 mW and ω_LO = ω₀ − Ω₀ with Ω₀ = 2π × 50 MHz respectively. The relation among ω₀, ω_p, and ω_LO is summarized in Fig. 1d. Furthermore, the SMF acts as a spatial mode filter to mainly extract the LP₀₁ component from the SR output (see Methods) as given by ${\mathbf{E}}_{{\mathrm{SR}}}^\prime \left( {{\mathbf{r}},t} \right) = \left( {{\hat{\mathbf{e}}}_x/2} \right)E_{{\mathrm{SR}}}^\prime \left( {{\mathbf{r}},t} \right)e^{ - i\omega _0t} + {\mathrm{c}}.{\mathrm{c}}.$, where a prime, such as E^', denotes the LP₀₁ component.

The heterodyne signal V_RF(t) is band-pass filtered with a bandwidth of 4 MHz to improve the signal-to-noise ratio (SNR) of the measurement (see Methods), and then is recorded by a digital storage oscilloscope (DSO). Upper panel of Fig. 2a shows a typical trace of the single-shot heterodyne signal V_RF(t), where a yellow-shaded region corresponds to a 500-ns-long pump pulse. The envelope exhibits a ringing behaviour, where the first pulse-shaped envelope with a 0.7-μs width corresponds to the pumping process while the following bursts are caused by the reabsorption and reemission of radiation between different parts of atomic cloud^33,34. The temporal widths of the bursts after the pump pulse are much shorter than the spontaneous emission lifetime 1/γ₀ = 21 μs of single atom, manifesting the occurrence of SR. In particular, the first-SR burst with a temporal $1/\sqrt e$ width of $\gamma _{{\mathrm{bw}}}^{ - 1} = 0.38$ μs gives a photon-emission enhancement factor of γ_bw/γ₀ = 55. Lower panel of Fig. 2a shows a temporal variation of SR power ∝|V_RF(t)|² in the logarithmic scale. Thanks to the heterodyne detection, we achieve the measurement dynamic range over five orders of magnitude. Such sensitive detection scheme will allow accessing orders of magnitude weaker subradiant process, which usually falls behind the SR process^35,36, by extending the measurement time sufficiently longer than $\gamma _0^{ - 1}$. Figure 2b shows the normalized burst width γ_bw/γ₀ as a function of the number of atoms N, indicating the characteristic feature of the SR or cooperative emission that is proportional to N. A linear fit to the data determines the single-atom cooperativity $\eta _{{\mathrm{bw}}} = \frac{{\gamma _{{\mathrm{bw}}}}}{{N\gamma _0}} = 4.9 \times 10^{ - 4}.$

Fig. 2

Superradiant ringing. a Ringing behaviour of the filtered heterodyne signal V_RF(t) (black) after applying the resonant (Δ_p = 0) pump pulse with the duration τ_p = 0.5 µs (indicated by a yellow-shaded region) and the amplitude E_p = ℏπ/(Dτ_p). N ≈ 9.4 × 10⁴ atoms are distributed over l_z = 0.85 mm along the fibre axis. The first 0.7-μs-long peak corresponds to the pumping process, which is followed by four superradiance (SR) bursts within 10 μs. The inset shows the zoom up to reveal the carrier oscillation at about Ω₀. The time-dependent SR power proportional to $V_{{\mathrm{RF}}}^2\left( t \right)$ is given by the red curve. The SR decay rate γ_bw of the first-SR burst is evaluated by the reciprocal of the temporal width at $1/\sqrt e$ of the peak amplitude. For N ≈ 9.4 × 10⁴, γ_bw reaches 55γ₀. b Experimentally measured burst width γ_bw as a function of the number of atoms N. The curve fitting gives the linear relation γ_bw/γ₀ = (4.9 × 10^−4) × N, manifesting the SR feature. The error bars, which are given for three representative data, show the standard deviations

Frequency shift

The signal V_RF(t) carries also frequency information. Indeed, the amplitude of SR light output from the SMF may be rewritten as $E_{{\mathrm{SR}}}^\prime \left( {{\mathbf{r}},t} \right) = A_{{\mathrm{SR}}}^\prime \left( {{\mathbf{r}},t} \right)e^{ - i\Delta _{{\mathrm{SR}}}t}$, where the complex function $A_{{\mathrm{SR}}}^\prime \left( {{\mathbf{r}},t} \right)$ is related to the envelope of V_RF(t) and $\Delta _{{\mathrm{SR}}} = \omega _{{\mathrm{SR}}} - \omega _0$ corresponds to the shift of the centre frequency ω_SR of the SR light from the atomic-transition frequency ω₀ (Fig. 1d). The frequency shift Δ_SR can be extracted by I/Q demodulating the carrier signal of V_RF(t) recorded by DSO (see Fig. 3a) into the in-phase V_I(t) and quadrature-phase V_Q(t) components, i.e. V_RF(t) = V_I(t) cos Ω₀t + V_Q(t) sin Ω₀t. Figure 3b shows an example of an amplitude $V\left( t \right) = \sqrt {V_{\mathrm{I}}\left( t \right)^2 + V_{\mathrm{Q}}\left( t \right)^2}$ and a phase θ(t) = ArcTan[V_Q(t)/V_I(t)] obtained from the I/Q signals, where we repeat the measurement of V_RF(t) ten times and average the data. The linearly varying phase θ(t) in time indicates the frequency shift of the emitted light Δ_SR = dθ/dt, which is nearly constant over 15 μs except the discontinuities at the nodes of the amplitudes that cause π-phase jump. The same I/Q signals are parametrically mapped onto the I/Q plane as shown in Fig. 3c, where the purple line corresponds to the trajectory from t = 1.48–8.38 μs. By taking a certain rotating frame Δ_SR = dθ/dt, the trajectory can be projected to a line as indicated by the blue colour (see Methods). We thus determine Δ_SR ≈ −2π × 93(10) kHz for the unexpanded atomic cloud with N = 9.4 × 10⁴(ρ = 8.8 × 10¹³â€‰cm^−3). A slight departure of the phase θ(t) from dashed lines in Fig. 3b, which corresponds to finite width of the blue trajectory in Fig. 3c, may indicate a frequency chirp in the emitted signal. Further investigation on the chirping effect will be given elsewhere. In the following, we analyse Δ_SR in the latter method as illustrated in Fig. 3c, assuming Δ_SR to be constant.

Fig. 3

Frequency shift of superradiant emission. a In-phase V_I(t) and quadrature-phase V_Q(t) components are obtained by demodulating the V_RF(t) signal stored in digital storage oscilloscope. b Example of an amplitude $V\left( t \right) = \sqrt {V_{\mathrm{I}}\left( t \right)^2 \, + \, V_{\mathrm{Q}}\left( t \right)^2}$ (black line) and a phase θ(t) = ArcTan[V_Q(t)/V_I(t)] (red line) measured for N ≈ 9.4 × 10⁴. A linear phase shift in time indicates the constant frequency shift of the emitted light. Dashed lines with dθ/dt = 2π × 93kHz are shown for guides to eyes. c Mapping V_I(t) and V_Q(t) into an I/Q plane. In order to get rid of the influence of the pump field and to ensure enough signal-to-noise ratio, V_I(t) and V_Q(t) within 1.48 μs ≤ t ≤ 8.38 μs are adopted to derive the frequency shift to be 93(10)kHz. d Dependence of the frequency shift Δ_SR = ω_SR − ω₀ on the average on-site atomic density ρ under the resonant-pump (Δ_p = 0) condition. Black (red) symbols correspond to unexpanded (expanded) atomic clouds. The error bars are given by the frequency where the value of the evaluation function is doubled. e. Δ_SR as a function of the pump-pulse detuning Δ_p measured for different numbers of atoms N. We apply unexpanded atomic clouds with l_z = 0.85 mm

We investigate the conditions that may affect the frequency shift Δ_SR of the SR field. Figure 3d shows the dependence of Δ_SR on the atomic density ρ for the resonant pumping Δ_p = 0. It is seen that Δ_SR is negative and |Δ_SR| grows up linearly with ρ. A linear fit finds the coefficient for the frequency shift to be Δ_SR/ρ ~ −1 × 10^−9 Hz cm³, which is consistent with the values measured in the absorption spectrum^25,37, suggesting the emission and absorption spectra are subject to the same density shift. We also measure Δ_SR as a function of the pump-field detuning Δ_p with the same amplitude and duration of the pump pulse as that of the π-pulse with Δ_p = 0. As shown in Fig. 3e, for a given density ρ, Δ_SR stays nearly constant over a wide range of Δ_p, which indicates that the shift Δ_SR arises entirely from the density shift and is independent of the pumping frequency ω_p. For N = 2 × 10⁵ atoms confined in the unexpanded lattice, the mean site-occupancy is about 10² and the interatomic separation is 300 nm, leading to a collective Lamb shift plus a Lorentz-Lorenz shift³⁸ of about γ₀. The frequency shift observed in experiment well exceeds any of them. Indeed, the resonant dipole-dipole interactions (RDDIs) between the atoms in the same lattice site primarily contribute to the observed density-dependent frequency shift. The linear dependence between Δ_SR and ρ is well reproduced by our numerical simulation (see below and Supplementary Note 5).

Efficiency of superradiance

We further consider the efficiency of SR that weighs the energy transfer from the pump field to the SR light. In the spontaneous emission, the emitted photons are randomly oriented in free space and thereby hardly coupled to the fibre modes. In contrast, the SR shows a well-defined direction depending upon the sample’s geometry. The multi-transverse-mode propagation in the HCPCF may support an in-fibre beam which matches the atomic cloud aligned along the fibre axis. As a result, the light power output from the fibre provides the attenuated pumping pulse by the absorption of atoms and the SR field emitted from the atoms. We investigate the efficiency of SR by measuring the total radiation power P_SR(t) for the applied pump power P_p(t).

The measurement procedure is illustrated in Fig. 1a. The light output from the HCPCF is split into two paths by a beam splitter. In one arm, a photomultiplier tube (PMT) measures the total light intensity to determine the coupling efficiency κ of the SR to multi-spatial modes while in the other arm, the beam enters the SMF. The light output from the SMF is in a superposition state of the fibre LP₀₁, LP₀₂, and LP₀₃ modes, where the LP₀₁ component dominates the weight. This light beam is detected by an avalanche photodiode (APD) so as to measure the coupling efficiency κ₀, which is mostly determined by the LP₀₁ mode, i.e. the efficiency of transferring the pumping power to a specific mode (see Methods). The photon-counting mode is applied for both measurements, where the neutral-density (ND) filters are inserted in front of the detectors to attenuate the signal down to the photon-counting level.

Figure 4a shows a typical data obtained by the PMT after averaging 25,000 measurements with N = 9.4 × 10⁴. The green curve displays the result in the absence of the atoms, i.e. the power of the pump-light P_p(t), while the red line denotes the measurement with the presence of lattice-confined atoms, i.e. the power of the light passing through the atoms and superradiantly emitted by the atoms P_SR(t). The blue line corresponds to the difference ΔP(t) = P_SR(t) − P_p(t), where the positive (negative) area is proportional to the total emitted P_em (absorbed P_ab) power. The efficiency κ is then given by the ratio κ = |P_em/P_ab|. Figure 4b shows the dependence of κ on the atomic number N for the unexpended (black solid squares) and expanded (red solid circles) atomic clouds. It is seen that in both cases κ goes up monotonically and approaches unity asymptotically as N is increased. The absorbed energy in the pumping process can be converted into either the SR radiation or the spontaneously emitted radiation. The former, whose power scales as N²Î³_SR with the average SR emission rate γ_SR of single atom, is collectively coupled to the HCPCF with an efficiency χ. In contrast, the latter is proportional to Nγ₀ and hardly matches the HCPCF modes. Thus, the efficiency κ can be formulated as $\kappa = \chi \frac{{\eta N}}{{1 + \eta N}}$, where η = γ_SR/γ₀ denotes the single-atom cooperativity. The curve fitting to filled symbols in Fig. 4b gives χ ≈ 1 and η ≈ 4.1 × 10^−5, which reasonably agrees with the geometrically estimated value η_Ω = NA²/4 ≈ 10^−4 with the numerical aperture NA = λ₀/πw₀ ≈ 0.02 for coupling to the guided fibre mode.

Fig. 4

Efficiency of superradiance. a Time evolution of radiation measured by a photomultiplier tube (PMT), where the pump pulse resonantly (Δ_p = 0) excites the unexpanded atomic cloud with N ≈ 9.4 × 10⁴ and l_z = 0.85 mm. The red (green) curve corresponds to the PMT signal with (without) atoms inside the fibre. The non-overlapping areas within the pumping regimes t ≤ τ_p and $t \, > \, \tau _{\mathrm{p}}$ denote the energies absorbed and emitted by the atoms, respectively. The difference between two curves is given by the blue line. b Total efficiency κ of the superradiance (SR) coupling to multiple transverse modes (shown by filled symbols) and efficiency κ₀ of the SR coupling to the fundamental LP₀₁ mode (empty symbols) as a function of the atomic number N for unexpanded (l_z = 0.85 mm, shown by black) and expanded (l_z = 2.1 mm, shown by red) atomic clouds. Two solid lines give the single-atom cooperativity η = 4.1 × 10^−5 (filled black), η = 3.2 × 10^−5 (filled red) with χ ≈ 1. The error bars representing the standard deviation are given for four data points

Similarly, we obtain the coupling efficiency κ₀ of the superposition state of the fibre LP₀₁, LP₀₂, and LP₀₃ modes, as shown in Fig. 4b by open square symbols. It is found that κ₀ is maximized around N_m ≈ 10⁵ and then decreases as N is further increased, suggesting that fibre higher-order modes take up more SR energy for N > N_m. The SR rate γ_SR may be separated into two parts, which respectively contribute to the fundamental LP₀₁- and high-order-mode SR fields, i.e. $\gamma _{{\mathrm{SR}}} = \gamma _{{\mathrm{SR}}}^{\left( \mathrm{f} \right)} + \gamma _{{\mathrm{SR}}}^{(\mathrm{h})}$, and $\kappa _0 = \left( {\gamma _{{\mathrm{SR}}}^{\left( {\mathrm{f}} \right)}/\gamma _{{\mathrm{SR}}}} \right)\kappa$. Since the atoms emit the guided light within the cross-sectional area of the atomic cloud, whose radius r_a is much smaller than the LP₀₁-mode beam radius w₀, the higher-order fibre modes (such as LP₀₂ and LP₀₃) with their central-peak radii smaller than w₀ get enhanced prior to LP₀₁ and hence $\gamma _{{\mathrm{SR}}}^{\left( {\mathrm{f}} \right)}\, < \, \gamma _{{\mathrm{SR}}}$. For a larger N, less power is transferred to LP₀₁ and $\gamma _{{\mathrm{SR}}}^{\left( {\mathrm{f}} \right)}$ becomes smaller, resulting in a reduced efficiency κ₀. Our numerical simulation qualitatively confirms the similar dependence of κ and κ₀ on N (see below and Supplementary Note 5). We find that η is about ten times smaller than η^(bw) that is derived based on the first-SR burst. This may be attributed to the fact that the first-SR burst does not behave exponentially.

Discussion

To understand the subtle mechanism behind the above experimental observations, we need to simulate the collective atom-light interplay under the experimental detection conditions. A similar in-fibre SR model has been studied in a recent theoretical work³⁹, where the physical system is simplified and only the collective decay of excited atoms, rather than the photon emission, is focused on. Here, we model the practical physical system by employing Maxwell-Bloch equations (MBEs)² and derive the in-fibre light field directly. The theoretical model takes into account the effects of the inhomogeneous atomic distribution, the multi-transverse-fibre-mode propagation, and the resonant dipole-dipole interactions^40,41 among atoms in the same lattice site (see Supplementary Note 1–4). The fibre modes are calculated via the full-vector finite element method⁴² and the 4th order Runge-Kutta technique is applied to solve MBEs. Choosing fibre modes is crucial for the numerical simulation. Since the SR light is emitted within the cross section of the atomic cloud, the higher-order fibre modes, whose intensity distributions have a central peak with a radius larger than or similar to the radial radius r_a of the atomic cloud, should be considered. We consider up to nine transverse fibre modes, among which the intensity patterns of the modes in Group I (i.e. LP₀₁, LP₀₂, and LP₀₃) are maximized at the central point of the fibre core while the others in Group II (i.e. ${\mathrm{LP}}_{{\mathrm{11}}}^{{\mathrm{a}},{\mathrm{b}}}$, ${\mathrm{LP}}_{{\mathrm{21}}}^{{\mathrm{a}},{\mathrm{b}}}$, and ${\mathrm{LP}}_{{\mathrm{31}}}^{{\mathrm{a}},{\mathrm{b}}}$) do not have central peaks (see Supplementary Note 4), in the simulation. The Gaussian-like fundamental mode LP₀₁ has a waist radius w₀ = 11.8 μm, at which the intensity drops to 1/e² of the maximum value. The maximum loss (caused by the limited light confinement strength of the cladding design) experienced by the selected modes over the atom-light interaction region is estimated to be lower than 10^−2 dB and is thus negligible. Besides these selected modes, other transverse modes may also propagate inside the actual fibre. However, including more fibre modes requires large amount of computer memory and computational time and may also mismatch the numerical-simulation and experimental results (see below).

Figure 5a displays the simulated SR temporal trace in the same conditions as Fig. 2a. The heterodyne signal V_cal(t) is presented in the black colour while the red line shows the envelope of the signal $\tilde V_{{\mathrm{cal}}}\left( t \right) \propto A_{{\mathrm{SR}}}^\prime \left( {{\mathbf{r}},t} \right)$ without the band-pass filtering. It is seen that the ringing behaviour is well replicated. The intervals between two adjacent nodes are consistent with the experimental measurement, proving the validity of simulation. As illustrated by $\tilde V_{{\mathrm{cal}}}\left( t \right)$, the first node occurs within the pumping process, meaning the SR starts before the pump pulse end. It is also found that the band-pass filter affects the measurement in two aspects: the rapid rising edges of the first two pulses in $\tilde V_{{\mathrm{cal}}}\left( t \right)$ are stretched, which explains the fact that the envelope peak corresponding to the pumping pulse in V_cal(t) is lower than that of the first-SR burst in V_cal(t); and the variation of V_cal(t) falls behind $\tilde V_{{\mathrm{cal}}}\left( t \right)$ by a retarded time of about 0.2 μs. It is worth noting that the multimode propagation is crucial to this ringing phenomenon. The superposition of multiple transverse modes can sustain an in-fibre beam radius similar to the transverse 1/e radius r_a of the atomic cloud, magnifying the collective atom-light interplay (see below and Supplementary Movie 1). We should also note that when extra transverse modes (for example, LP₀₄) with a central-peak radius smaller than that of LP₀₃ are involved in Group I in the numerical simulation, the separations between inter-adjacent envelop nodes of V_cal(t) are strongly shortened, which is inconsistent with the experimental measurement and manifests the negligible effect of these high-order modes in the atom-light interface.

Fig. 5

Numerical simulation of the ringing behaviour. a Numerical result of heterodyne signal V_cal(t) (black) for N ≈ 9.4 × 10⁴ assuming the resonant (Δ_p = 0) pump pulse with the duration τ_p = 0.5 µs and the amplitude E_p = ℏπ/(Dτ_p), which is corresponding to the experimental measurement. In comparison, the envelope of the raw signal $\tilde V_{{\mathrm{cal}}}\left( t \right)$ (red) before entering the band-pass filter is also presented. Here, we have used the coupling efficiencies of the fibre LP₀₁, LP₀₂, and LP₀₃ modes to SMF c_LP01 = 0.60, c_LP02 = 0.33, and c_LP03 = 0.05. b Time evolution of fibre-mode weights. The weights of components LP₀₁, LP₀₂, and LP₀₃ in the output SR field E_SR(r, t) vary in time. Examples of SR pattern at several selected times are presented. The pattern deviates from the Gaussian-like distribution at each node of the SR ringing shown in Fig. 5a. Away from ringing nodes, the transverse beam radius of the output SR is about three times smaller than that of the fundamental LP₀₁ fibre mode. c. Numerical first-SR-burst width γ_bw vs. the number of atoms N. The solid line shows the curve fitting. All system parameters are same to Fig. 5a

As we have pointed out above, the multi-transverse-fibre-mode propagation plays an important role in the atom-light coupling. Figure 5b depicts the time-dependent weights of several fibre modes and SR patterns at selected times. The more detailed time evolution of SR pattern can be found in Supplementary Movie 1. The results reveal a unique feature whereby the sequence of the pump absorption and the SR bursts is imprinted in the time evolution of the intensity profile of the fibre-guided modes. The initial state of the pump coupling into the fibre is illustrated by the LP₀₁-like mode intensity profile at t = 0 μs. The doughnut-shaped pattern during the pumping pulse denotes the strong absorption within the cross section of atomic cloud. The SR events are represented by a Gaussian-like pattern with a beam radius smaller than w₀. This profile results from a linear combination of the LP₀₁, LP₀₂, and LP₀₃ modes. Hence, the SR in multimode HCPCF can sustain a strong atom-light interplay, which differs from the situation of the atoms coupled to a single-mode fibre or a cavity^4,5. In addition, around each node E_SR(r, t) ≈ 0 of SR, the transverse profile becomes different from the Gaussian-like distribution, indicating the strong influence from the non-central-peak fibre modes in Group II. Only when away from the nodes, the modes in Group I get enhanced prior to the modes in Group II. Moreover, from the numerical simulation one may derive the SR decay rate γ_bw corresponding to the first-SR burst. Our theoretical result predicts a linear dependence of γ_bw on the number of atoms N (see Fig. 5c), γ_bw/γ₀ ≈ 5.4 × 10^−4N, consistent with the experimental measurement in Fig. 2b.

Our numerical model also reproduces other experimental results as illustrated in Fig. 6. Performing the Fourier transform on the numerically-simulated amplitude E_SR(r, t), one obtains the SR spectrum $S\left( \omega \right) \propto \left| {\int_{\tau _{\mathrm{p}}}^\infty {\left[ {{\int} {{\int} {E_{{\mathrm{SR}}}\left( {{\mathbf{r}},t} \right)d^2{\mathbf{r}}} } } \right]} e^{i\left( {\omega - \omega _0} \right)t}dt} \right|^2$ (see Fig. 6a), from which the frequency shift Δ_SR is read out by deducting the envelope-oscillation frequency introduced by the ringing behaviour. The dependence of Δ_SR on the average on-site atomic density ρ is plotted in Fig. 6b, which agrees with the experimental results. We also theoretically reproduce the counterparts of Fig. 4a, b in Fig. 6c, d. In comparison, the numerical results of the efficiencies of SR coupling to all fibre modes (κ) and to the fundamental mode (κ₀) are both lower than the experimental results. This is mainly because the calculated transverse modes are not exactly the same as those existing in the actual fibre. A slight mismatch can lead to an apparent difference. Furthermore, the insufficient number of the fibre transverse modes in Group II joining in the atom-light interaction may give rise to an erroneous increment of the power absorption within the pumping process, resulting in the reduced κ and κ₀. Nevertheless, the numerical and experimental results still exhibit similar behaviours.

Fig. 6

Numerical simulation for the resonantly-pumped unexpanded atomic cloud. a Superradiance (SR) spectrum S(ω) with the atom number N = 9.4 × 10⁴. The spectrum peak is shifted to the red side of ω₀ by 2π × 152 kHz. Deducting the envelope-oscillation frequency of the SR field, which is 2π × 54 kHz, gives Δ_SR = −2π × 98 kHz. b Frequency shift Δ_SR vs. the average on-site atomic density ρ. The filled black symbols correspond to the numerical results whose curve fitting, i.e. Δ_SR/ρ = 2π × 1.07 × 10^−12 kHz cm³, is given by the solid line. c Numerically-derived time-dependent powers of pump beam P_p(t) (green), SR field P_SR(t) (red) and their difference ΔP(t) (blue). The blue line denotes the energies emitted (absorbed) by the atoms when the ΔP(t) is positive (negative). d Total efficiency κ of the SR coupling to multiple transverse modes (filled symbols) and efficiency κ₀ of the SR coupling to the fundamental LP₀₁ mode (empty symbols) as a function of the atomic number N. The curve fitting of κ leads to η = 3.6 × 10^−5 and the efficiency χ = 0.87

In summary, we have demonstrated the multimode dynamics of collective spontaneous emission based on an engineered in-fibre atom-light-interface architecture. This platform may be applicable to explore more exotic quantum phenomena such as photonic band gaps in 1D ordered atomic structure, quantum jumps of many-body systems with long-range interaction, the open-system Dicke-model phase transition⁴³, and the competing process between virtual and real photon exchanges for the atom-atom interaction⁴⁴. Confining two atomic clouds at two ends of a long fibre allows the fibre-based remote macro-entanglement and quantum information transmission⁴⁵. Especially, substituting the long-lived ¹S₀−³P₀ clock transition for the intercombination transition enables the miniaturization of the physics package of optical lattice clocks²⁵. This compact and transportable setup may facilitate the high-level comparison of optical clocks and time keeping⁴⁶.

Methods

Loading ultracold atoms into HCPCF

Ultracold ⁸⁸Sr atoms are prepared via the two-stage MOT²⁹ and then captured at the antinodes of the optical lattice that is formed by a pair of linearly-polarized counter-propagating λ_L = 813 nm lasers. Both lattice beams match the fundamental LP₀₁ fibre mode with a power of 200 mW. The lattice antinodes are moved along the axial direction (i.e. the z-axis) at a velocity of v(t) = δ_L(t)λ_L/2 by introducing time-dependent detuning |δ_L(t)| ≤ 90 kHz between two lattice lasers. The atomic cloud is transported by a distance of 3.25 mm within a moving time of about 160 ms at a maximum velocity of 32.5 mms^−1 The temperature of atoms is measured to be ~4.0 μK by the Doppler width of the spectrum, which suggests the averaged vibrational occupations of n_a ~ 0.5 and n_r ~ 33, for the axial and radial motions, respectively.

The atom-light interaction length can be extended via expanding the atomic cloud inside the fibre by shutting off either of lattice lasers to remove the periodic lattice potential. Then, the atomic cloud experiences a Gaussian-like optical potential in the x − y plane and an axial expansion with the velocity satisfying the Maxwell-Boltzmann distribution. After a short duration, the lattice laser, which has been turned off, is switched on again. The atomic cloud with an expanded axial width l_z is recaptured by the lattice potential. During the expanding process, the atom loss is negligible. Figure 1b shows the distribution of the atomic cloud with and without 80-ms-long expansion. We then adiabatically reduced the lattice power to 20 mW, which is one tenth of the initial power, to moderate the density shift and residual light shift for observing SR.

Mode coupling between HCPCF and SMF

A portion of SR light output from the HCPCF enters the SMF (see Fig. 1a). The optical path between HCPCF and SMF is aligned to maximize the coupling efficiency to the fundamental LP₀₁ mode of the HCPCF, which is about 50%. The Gaussian-like ground mode ϕ₀(r) of the SMF approximates as $\phi _0\left( {\mathbf{r}} \right) = \sqrt {2/\left( {\pi \tilde w_0^2} \right)} e^{ - \left| {\mathbf{r}} \right|^2/\tilde w_0^2}$ with the beam radius of $\tilde w_0 = 0.4w_0$. The coupling efficiency between LP₀₁ and ϕ₀(r) is calculated to be $c_{{\mathrm{LP}}_{{\mathrm{01}}}} = \left| {{\int} {{\int} {\psi _0^ \ast \left( {\mathbf{r}} \right)\phi _0\left( {\mathbf{r}} \right)d^2{\mathbf{r}}} } } \right|^2 = 0.60$, consistent with the experimental measurement (~50%). Here ψ₀(r) corresponds to the transverse distribution of LP₀₁. Similarly, the coupling efficiencies of selected high-order modes of HCPCF to SMF are given by $c_{{\mathrm{LP}}_{{\mathrm{02}}}} = 0.33$, $c_{{\mathrm{LP}}_{{\mathrm{03}}}} = 0.05$, and $c_{{\mathrm{LP}}_{{\mathrm{11}}}^{{\mathrm{a}},{\mathrm{b}}}} = c_{{\mathrm{LP}}_{{\mathrm{21}}}^{{\mathrm{a}},{\mathrm{b}}}} = c_{{\mathrm{LP}}_{{\mathrm{31}}}^{{\mathrm{a}},{\mathrm{b}}}} = 0.00$. Thus, the fundamental mode LP₀₁ contributes dominantly in energy at the output end of SMF.

Balanced-heterodyne detection

As depicted in Fig. 1c, two Si PIN photodiodes are used to convert the optical-frequency interference down to a radio-frequency (RF) signal $V_{{\mathrm{RF}}}\left( t \right) \propto {\mathrm{Re}}\left\{ {\left[{\int} {{\int} {E_{{\mathrm{LO}}}\left( {\mathbf{r}} \right)} } {E_{{\mathrm{SR}}}^\prime \left( {{\mathbf{r}},t} \right)d^2{\mathbf{r}}} \right]e^{ - i{\mathrm{\Omega }}_0t}} \right\}$, where Re{…} gives the real part. Setting the power of the local optical oscillator E_LO at P_LO ≈ 1 mW, whose shot noise well exceeds the circuit noise, the measurement SNR is approximately given by $\sqrt {2KP_{{\mathrm{SR}}}/\left( {\pi e{\mathrm{\Delta }}f} \right)}$ where K = 0.42 AW^−1 is the photodiodes’ responsivity. We further adopt a band-pass filter (see Fig. 1c) with a FWHM bandwidth of Δf = 4 MHz and a centre frequency of Ω₀/2π = 50 MHz to improve the SNR for observing the time-varying radiation amplitude E_SR(r, t).

Deriving radiation frequency shift Δ_SR

The 10-times-averaged V_RF(t) is demodulated into the in-phase V_I(t) and quadrature-phase V_Q(t) components. We then map the point V_I(t) + iV_Q(t) within the period from 1.48 μs to 8.38 μs into a complex coordinate system (see the purple line in Fig. 3c). Assuming the frequency shift Δ_SR is time independent, the plot can be formulated as $V\left( t \right)e^{i\Delta _{{\mathrm{SR}}}t + \theta _0}$, where $V\left( t \right) = \sqrt {V_{\mathrm{I}}^2\left( t \right) + V_{\mathrm{Q}}^2\left( t \right)}$ denotes the distance from the point V_I(t) + iV_Q(t) to the origin of the coordinate and θ₀ is the initial phase, because the phase components are demodulated with the frequency of ω₀ (see Figs. 1d and 3a). If we rotate the frame around the origin of the coordinate at a frequency of ω′ and with a phase offset of θ', the plot may be re-expressed as $V\left( t \right)e^{i\left( {{\mathrm{\Delta }}_{{\mathrm{SR}}} - \omega^\prime } \right)t + \left( {\theta _0 - \theta^\prime } \right)}$. By assigning $\omega^\prime = \Delta _{{\mathrm{SR}}}$ and θ′ = θ₀, the point moves only on the real axis. To find Δ_SR from two phase components, we prepare an evaluation function that integrates the square of the imaginary part with respect to time and adopt ω′ that minimizes the evaluation function related to Δ_SR. The blue line in Fig. 3c shows the plot in the rotating frame that suppresses the imaginary part and Δ_SR is determined to be $\Delta _{{\mathrm{SR}}} \approx - 2\pi \times 93(10)\,{\mathrm{kHz}}$. Here the uncertainty is estimated by the frequency where the value of the evaluation function is doubled.