Superradiance from lattice-confined atoms inside hollow core fibre

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. Superradiance, a collective light emission produced by an ensemble of excited atoms has been extensively investigated. The authors study the collective emission from ultracold strontium atoms trapped inside a hollow core fibre, which may offer a miniaturized platform for exploring superradiance.

S uperradiance (SR) has been the subject of active research since the pioneering work of Dicke in 1954 1 . Such a collective effect arises from an ensemble of two-level emitters spontaneously locking their dipole oscillations in phase 2 , giving rise to an enhanced decay rate γ SR many times faster than the spontaneous emission rate γ 0 of individual emitters and a strong peak intensity of radiation in quadratic proportion to the excitedstate 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) 3 or van der Waals dephasing 2 . 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 8 . 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 solidstate 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 phaselocked 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 subradiation-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 19 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 /γ 0 , which is a geometric parameter characterizing the emissive coupling of an atom to the cavity mode, is much less than unity 8 .
Research at the interdisciplinary frontier between the fields of cold atoms and photonic materials reports promising results in controlling ultracold atoms with photonic waveguides 20 , 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 21 . 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 22 , 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 R ¼ πw 2 0 =λ 0 at the wavelength λ 0 while still keeping a small light-beam waist of radius w 0 . Successful combination of ultracold atoms and photonic waveguide has been exemplified by the works with nanofibers 23,24 , photonic waveguides 19 , and hollow-core fibers 25 . The long-range atom-atom interaction is accomplished via guided photons of the so-called alligator photonic crystal waveguide 19 , and the resulting single-atom cooperativity reaches η = 0.34. Ultracold atoms are successfully loaded inside a hollowcore photonic crystal fibre (HCPCF) and trapped in a 1D magicwavelength optical lattice 25 , 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 26 . We couple an HCPCF-transmitted SR to a singlemode 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 photonmediated 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 5 with the single-atom cooperativity η = 4.1 × 10 −5 available for the HCPCF, the N-atom cooperativity parameter 8 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 Kagomelattice HCPCF 26 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 27 and supports multiple transverse modes whose transverse profiles are very close to those of a dielectric capillary 28 . The Gaussianlike fundamental LP 01 core mode has a 1/e 2 waist radius w 0 = 11.8 μm (1/e radius is 8.34 μm). 88 Sr atoms are laser cooled and trapped in a two-stage magneto-optical trap (MOT) 29 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 01 mode and operate at the wavelength λ L = 813 nm. A bias magnetic field B 0 ¼ 0:14 mT ð Þ e x is applied to define the quantization axis. We focus on the 1 S 0 − 3 P 1 (m = 0) intercombination transition of 88 Sr at the frequency ω 0 = 2π × 435 THz with the natural linewidth γ 0 = 2π × 7.5 kHz and the dipole moment D = 0.087 ea 0 (elementary charge e and Bohr radius a 0 ). The lattice laser polarization is tuned to cancel out the ac Stark shift difference between 1 S 0 and 3 P 1 (m = 0) states caused by the optical lattice 21,25 .
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 25 . 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 3 and 5.2 × 10 3 , respectively (see Methods). Each lattice site forms a pancakeshaped 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 ρ ¼ N a = πr 2 a l a À Á 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 30 hardly affects the atom-light interaction as the lattice period λ L /2 does not fulfil the Bragg condition for the transition wavelength λ 0 .
Superradiance ringing. All atoms are initialized in the 1 S 0 ground state. We apply anê 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 01 fibre core mode, and the in-fibre amplitude E p is set to E p = ħπ/Dτ p with ħ being the reduced Planck constant. The pumplaser 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 − ω 0 from the atomic resonance ω 0 . The observed SR field at the fibre output end can be written as E SR r; t ð Þ ¼ê x =2 ð ÞE SR r; t ð Þe Àiω 0 t þ c:c. with the radial coordinate r ¼ xê x þ yê 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 SR t ð Þ / R R jE SR ðr; tÞj 2 d 2 r is typically on the order of nW, requiring the low noise detection. To address this issue, we utilize the balanced-heterodyne technique 32 (see Fig. 1c). The SR output is coupled to a single-mode fibre (SMF) and mixed with a local oscillator E LO r; t ð Þ ¼ê x E LO r ð Þcosω 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 = ω 0 − Ω 0 with Ω 0 = 2π × 50 MHz respectively. The relation among ω 0 , ω p , and ω LO is summarized in Fig. 1d. Furthermore, the SMF acts as a spatial mode filter to mainly extract the LP 01 component from the SR output (see Methods) as given by E 0 SR r; t ð Þ ¼ê x =2 ð ÞE 0 SR r; t ð Þe Àiω 0 t þ c:c:, where a prime, such as E ' , denotes the LP 01 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 yellowshaded 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/γ 0 = 21 μs of single atom, manifesting the occurrence of SR. In particular, the first-SR burst with a temporal 1= ffiffi e p width of γ À1 bw ¼ 0:38 μs gives a photon-emission enhancement factor of γ bw /γ 0 = 55. Lower panel of Fig. 2a shows a temporal variation of SR power ∝|V RF (t)| 2 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 γ À1 0 . Figure 2b shows the normalized burst width γ bw /γ 0 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 Hollow-core photonic crystal fiber (HCPCF) a Schematic of experimental setup. Ultracold 88 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 = (0.14mT)ê 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 01 -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 1 S 0 − 1 P 1 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 ω LO t, while the other is used as the pump light with the amplitude E p at the frequency ω p = ω 0 + Δ 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 MHz is band-pass-filtered (BPF) and recorded by a digital storage oscilloscope (DSO). d Frequency relations. In respect to the 1 S 0 − 3 P 1 (m = 0) transition frequency ω 0 , the SR, pump-light, and local-oscillator frequencies are given by ω SR = ω 0 + Δ SR , ω p = ω 0 + Δ p , and ω LO = ω 0 − Ω 0 , respectively 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 0 where the complex function A 0 SR r; t ð Þ is related to the envelope of V RF (t) and Δ SR ¼ ω SR À ω 0 corresponds to the shift of the centre frequency ω SR of the SR light from the atomic-transition frequency ω 0 (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 Figure 3b shows an example of an amplitude V t ð Þ ¼ 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 4 (ρ = 8.8 × 10 13 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.
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 3 , 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 5 atoms confined in the unexpanded lattice, the mean site-occupancy is about 10 2 and the interatomic separation is 300 nm, leading to a collective Lamb shift plus a Lorentz-Lorenz shift 38 of about γ 0 . 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 multitransverse-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 01 , LP 02 , and LP 03 modes, where the LP 01 component dominates the weight. This light beam is detected by an avalanche photodiode (APD) so as to measure the coupling efficiency κ 0 , which is mostly determined by the LP 01 mode, i.e. the efficiency of transferring the pumping power to a specific mode (see 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 4 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 Ω 0 . The time-dependent SR power proportional to V 2 RF t ð Þ 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= ffiffi ffi e p of the peak amplitude. For N ≈ 9.4 × 10 4 , γ bw reaches 55γ 0 . b Experimentally measured burst width γ bw as a function of the number of atoms N. The curve fitting gives the linear relation γ bw /γ 0 = (4.9 × 10 −4 ) × N, manifesting the SR feature. The error bars, which are given for three representative data, show the standard deviations 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 4 . 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 2 γ 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γ 0 and hardly matches the HCPCF modes. Thus, the efficiency κ can be formulated as κ ¼ χ ηN 1þηN , where η = γ SR /γ 0 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 η Ω = Similarly, we obtain the coupling efficiency κ 0 of the superposition state of the fibre LP 01 , LP 02 , and LP 03 modes, as shown in Fig. 4b by open square symbols. It is found that κ 0 is maximized around N m ≈ 10 5 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 01 -and high-order-mode SR fields, i.e. γ SR ¼ γ 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 01 -mode beam radius w 0 , the higherorder fibre modes (such as LP 02 and LP 03 ) with their central-peak radii smaller than w 0 get enhanced prior to LP 01 and hence γ f ð Þ SR < γ SR . For a larger N, less power is transferred to LP 01 and γ f ð Þ SR becomes smaller, resulting in a reduced efficiency κ 0 . Our numerical simulation qualitatively confirms the similar dependence of κ and κ 0 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 39 , 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) 2 and derive the infibre light field directly. The theoretical model takes into account the effects of the inhomogeneous atomic distribution, the multitransverse-fibre-mode propagation, and the resonant dipoledipole 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 42 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 higherorder 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 01 , LP 02 , and LP 03 ) are maximized at the central point of the fibre core while the others in Group II (i.e. LP a;b 11 , LP a;b 21 , and LP a;b 31 ) do not have central peaks (see Supplementary Note 4), in the simulation. The Gaussian-like fundamental mode LP 01 has a waist radius w 0 = 11.8 μm, at which the intensity drops to 1/e 2 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Ṽ cal t ð Þ / A 0 SR r; t ð Þ 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 bỹ V cal t ð Þ, 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Ṽ cal t ð Þ 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 V cal t ð Þ 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 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 > τ 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 κ 0 of the SR coupling to the fundamental LP 01 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 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 04 ) with a central-peak radius smaller than that of LP 03 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.
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 01 -like mode intensity profile at t = 0 μs. The doughnut-  , which is corresponding to the experimental measurement. In comparison, the envelope of the raw signalṼ cal t ð Þ (red) before entering the band-pass filter is also presented. Here, we have used the coupling efficiencies of the fibre LP 01 , LP 02 , and LP 03 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 01 , LP 02 , and LP 03 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 01 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 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 0 . This profile results from a linear combination of the LP 01 , LP 02 , and LP 03 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 singlemode 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 /γ 0 ≈ 5.4 × 10 −4 N, 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 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 (κ 0 ) 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 κ 0 . Nevertheless, the numerical and experimental results still exhibit similar behaviours.
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 43 Fig. 6 Numerical simulation for the resonantly-pumped unexpanded atomic cloud. a Superradiance (SR) spectrum S(ω) with the atom number N = 9.4 × 10 4 . The spectrum peak is shifted to the red side of ω 0 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 3 , 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 κ 0 of the SR coupling to the fundamental LP 01 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 real photon exchanges for the atom-atom interaction 44 . Confining two atomic clouds at two ends of a long fibre allows the fibre-based remote macro-entanglement and quantum information transmission 45 . Especially, substituting the long-lived 1 S 0 − 3 P 0 clock transition for the intercombination transition enables the miniaturization of the physics package of optical lattice clocks 25 . This compact and transportable setup may facilitate the high-level comparison of optical clocks and time keeping 46 .

Methods
Loading ultracold atoms into HCPCF. Ultracold 88 Sr atoms are prepared via the two-stage MOT 29 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 01 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 a0 .5 and n r~3 3, 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 01 mode of the HCPCF, which is about 50%. The Gaussian-like ground mode ϕ 0 (r) of the SMF approximates as ϕ 0 r ð Þ ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi 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 RF t ð Þ / Re R R E LO r ð ÞE 0 SR r; t ð Þd 2 r Â Ã e ÀiΩ 0 t È É , 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 ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2KP SR = πeΔf ð Þ p 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 Ω 0 /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 t ð Þe iΔ SR tþθ 0 , where q denotes the distance from the point V I (t) + iV Q (t) to the origin of the coordinate and θ 0 is the initial phase, because the phase components are demodulated with the frequency of ω 0 (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 t ð Þe i Δ SR Àω 0 ð Þ tþ θ 0 Àθ 0 ð Þ . By assigning ω 0 ¼ Δ SR and θ′ = θ 0 , 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 Δ SR % À2π 93ð10Þ kHz. Here the uncertainty is estimated by the frequency where the value of the evaluation function is doubled.

Data availability
All data and computer code supporting the findings of this study are available from the corresponding author on reasonable request. Received: 7 June 2019; Accepted: 8 October 2019;