This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Long Light Storage Time in an Optical Fiber

Wui Seng Leong, Mingjie Xin, Chang Huang, Zilong Chen    Shau-Yu Lan sylan@ntu.edu.sg Division of Physics and Applied Physics, School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371, Singapore
Abstract

Light storage in an optical fiber is an attractive component in quantum optical delay line technologies. Although silica-core optical fibers are excellent in transmitting broadband optical signals, it is challenging to tailor their dispersive property to slow down a light pulse or store it in the silica-core for a long delay time. Coupling a dispersive and coherent medium with an optical fiber is promising in supporting long optical delay. Here, we load cold Rb atomic vapor into an optical trap inside a hollow-core photonic crystal fiber, and store the phase of the light in a long-lived spin-wave formed by atoms and retrieve it after a fully controllable delay time using electromagnetically-induced-transparency (EIT). We achieve over 50 ms of storage time and the result is equivalent to 8.7×\times10-5 dB μ\mus-1 of propagation loss in an optical fiber. Our demonstration could be used for buffering and regulating classical and quantum information flow between remote networks.

pacs:
Valid PACS appear here
preprint: APS/123-QED

Optical delay lines or optical buffers play important roles in long-distance quantum communication networks for storing, delaying, and, thus, exchanging information between different quantum nodes. The superb performance of optical fibers as transmission lines has driven the consideration of integrating these functionalities into the optical fibers themselves Th ; Geh . One of the challenges is maintaining the coherence of the medium in the fiber in which the information of light is encoded. The direct use of solid-core materials of telecommunication band fibers via stimulated Brillouin scattering and doped erbium ions to store light has been demonstrated to tens of nanoseconds Zhu ; Sag ; Jin ; Sag2 . Despite their tens of GHz bandwidth and telecom wavelength operation, the short-lived acoustic waves in the materials and spin coherence of doped erbium in the respective settings limit the performance of the light storage.

Alternatively, interfacing long-lived atomic spin states with optical fibers could provide a route to achieving long storage time. Loading room temperature atomic vapor into a hollow-core photonic crystal fiber Spr and cold atoms near the surface of an optical nanofiber Gou have achieved storage times of about tens of nanoseconds and few microseconds, respectively. They are limited by the transit time of atoms through the optical modes. Confining atoms in the evanescent field of an optical nanofiber Say ; Cor and guided field of a hollow-core fiber Bla ; Pet could eliminate decoherence from transit time but introduce additional decoherence from the confining fields. Both systems have shown storage time of about a few microseconds. To overcome the decoherence, we use a state insensitive optical potential to confine 85Rb atoms in the guiding mode of a 4-cm-long hollow-core photonic crystal fiber Deb ; Alh . The phase of the optical pulse is mapped onto the atomic spin wave formed by a pair of long-lived hyperfine ground states and retrieved with a controllable delay using EIT. We demonstrate light storage over 50 ms with a bandwidth of 1 MHz.

Refer to caption
Figure 1: Details of the experimental setup. NPBS: non-polarizing beam splitter. PBS: polarizing beam splitter. HWP: half-wave plate. QWP: quarter-wave plate. PM fiber: polarization-maintaining fiber. APD: avalanche photodiode. The QWPs are used to adjust the ellipticity of the dipole beams to induce the vector light shift. In order to combine the linearly polarized EIT beams and the elliptically polarized dipole beam, we use a 90%\% reflection and 10%\% transmission NPBS. A 50:50 NPBS is then used to separate the EIT beams from the counter-propagating dipole beam. We use a temperature-stabilized solid-state etalon to minimize leakage of the control beam to the APD. The depump light on the FF=3 to FF’=3 D2 line is used to prepare the atoms in the FF=2 state. The optical pump beam is resonant on the FF=2 to FF’=2 D1 line. The first control pulse is for converting the probe pulse into an atomic spin-wave and the second is for retrieving it.
Refer to caption
Figure 2: Input, slowed, and delayed light pulses through the hollow-core fiber. (a) Measured slowed light when the control beam is permanently on. The input pulse (black squares) is fitted by a Gaussian function with a FWHM of 304 ns, and the slowed pulse (blue circles) is fitted by a Gaussian function with a FWHM of 415 ns. The peak of the fitted Gaussian function of slowed light is delayed by 173 ns. (b) Numerical calculations of slowed pulses. The blue circles are calculations of the slowed pulse with non-uniform profile of the control beam, the probe beam, and the atom density in the fiber. The green squares are the calculations of the slowed light with uniform EIT beam profiles and atom density distribution. (c) The measured temporal profiles of the leaked and delayed light pulses. The control pulse is switched off when the probe pulse is propagating inside the atomic medium. The bandwidth of the delayed pulse is defined from its FWHM of 160 ns.

A cold 85Rb atomic ensemble is prepared by a magneto-optical trap (MOT) \sim5 millimeters above a 4-cm-long open-end hollow-core photonic crystal fiber, as shown in Fig. 1. After sub-Doppler cooling, the atomic ensemble is loaded into an optical trap guided by the fiber with a temperature of 10 μ\muK. The wavelength of the optical dipole beam is 821 nm with 250 mW power, and the measured optical depth (DD) is about 40 by probing all the Zeeman states of the 85Rb D2 line FF=3 to FF’=4 transition. When the atoms are in the fiber, we optically pump them into the FF=2, mm=0 state by a π\pi-polarized light and a repump light coupled through the fiber. The quantization axis is defined by a 2 G magnetic field along the fiber axis and the magnetic field is also used to lift the degeneracy of the Zeeman sublevels.

Our EIT scheme is formed by a three-level Λ\Lambda configuration Fle operating on the hyperfine clock states of 85Rb for long coherence times. The probe beam is resonant on ||1\rangle=|F|F=2, mm=0\rangle to ||3\rangle=|F|F’=3, mm=±\pm1\rangle and the control beam is resonant on ||2\rangle=|F|F=3 mm=0\rangle to ||3\rangle. A 795 nm diode laser is tuned to the control beam frequency, and a small portion of the power is split to set up the probe beam by an electro-optical modulator around 3 GHz followed by an etalon to filter out the undesired sidebands. Both beams are perpendicularly linearly polarized, and co-propagate inside the fiber. After exiting the fiber, the probe beam is separated from the control beam by polarization and etalon filtering and detected with an avalanche photodiode. The control and probe beams power inside the fiber are 300 nW and 40 nW, respectively, and the frequencies and the power of both beams are controlled by two independent and phase coherent acousto-optical modulators (AOMs). The 80 MHz radio frequency for the probe AOM is modulated by a Gaussian pulse from a waveform synthesizer. The 80 MHz radio frequency for the control AOM is modulated by a square pulse.

When the atoms are under free fall inside the fiber, we measure the group delay of the probe pulse, as shown in Fig. 2(a). The intensity of the input probe pulse can be approximated with a Gaussian temporal profile of IpI_{\textrm{p}}(tt)=Ip0I_{\textrm{p0}}exp(-4(ln2)t2t^{2}/Tp2T_{\textrm{p}}^{2}), where Ip0I_{\textrm{p0}} is the maximum intensity, TpT_{\textrm{p}}=304 ns is the FWHM duration, and tt is the temporal coordinate. A Gaussian function is fitted to the transmitted probe pulse to determine the group delay TdT_{\textrm{d}}=174 ns with 24%\% efficiency of transmission.

Refer to caption
Figure 3: Storage efficiency as a function of time. (a) The grey circles are the measurements of the storage efficiency without light shift cancellation. The function AA(1+(TT/TcT_{\textrm{c}})2)-3/2 is fitted to the data with TcT_{\textrm{c}}=823(10) μ\mus, where AA is the efficiency at TT=0. The decay time is limited by the differential ac Stark shift of the dipole potential on the atomic spin coherence. For comparison, the curve in purple represents the 0.03 dB μ\mus-1=0.15 dB km-1 loss of 1550 nm light propagating in a solid core fiber. The error bars are the standard errors of 40 experimental runs. (b) Measurements of the efficiency with light shift cancellation (orange circles). An exponential decay function is fitted to the data with a 1/ee time of 20(1) ms. The curve in purple is the same as in (a). The error bars are the standard errors of 40 experimental runs. (c) Measurements of normalized storage efficiency for a free-falling atomic spin-wave with light-shift cancellation. The control and probe pulses are sent when atoms are in motion inside the fiber.

We numerically calculate the group delay and the efficiency of the transmitted light by integrating over the inhomogeneous profiles of the atomic cloud in the dipole trap, the probe beam, and the control beam inside the fiber from the center of the fiber to the wall of the fiber WW=31.5 μ\mum Hsi . The intensity of the probe pulse after transmitting through the atomic ensemble is obtained as

0W2πr×(Ωpβ×exp(γ21DΓΩc2+4γ31γ21)×exp(2(ln2)(tTdβTp)2))2dr,\begin{split}\int_{0}^{W}2\pi r\times(\frac{\Omega_{\textrm{p}}}{\beta}\times\exp(\frac{-\gamma_{\textrm{21}}D\Gamma}{\Omega_{\textrm{c}}^{2}+4\gamma_{31}\gamma_{21}})\\ \times\exp(-2(\ln 2)(\frac{t-T_{\textrm{d}}}{\beta T_{\textrm{p}}})^{2}))^{2}dr,\end{split} (1)

where the factor

β=1+32(ln2)DΓ(γ31Ωc2+2γ312γ21+4γ213)Tp2(Ωc2+4γ31γ21)3,\displaystyle\beta=\sqrt{1+\frac{32(\ln 2)D\Gamma(\gamma_{31}\Omega_{\textrm{c}}^{2}+2\gamma_{31}^{2}\gamma_{21}+4\gamma_{21}^{3})}{T_{\textrm{p}}^{2}(\Omega_{\textrm{c}}^{2}+4\gamma_{31}\gamma_{21})^{3}}}, (2)

and the group delay Td=DΓ(Ωc2+4γ212)/(Ωc2+4γ31γ21)2T_{\textrm{d}}=D\Gamma(\Omega_{c}^{2}+4\gamma_{21}^{2})/(\Omega_{\textrm{c}}^{2}+4\gamma_{31}\gamma_{21})^{2}. The normalized Rabi frequency of the probe pulse is assumed to have a Gaussian distribution Ωp\Omega_{\textrm{p}}=exp\exp⁡(-r2r^{2}/R2R^{2}), where rr is the radial coordinate of the fibre core and RR=22 μ\mum is the 1/e2e^{2} mode field radius. The ground state coherence γ21\gamma_{21}, which is mainly due to the inhomogeneous differential ac Stark shift of the control and dipole beams, is fixed at 2π×\pi\times(2×\times105 Hz) based on the measurement of EIT linewidth and ground and excited states coherence γ31\gamma_{31} is fixed at half of the 85Rb D1 line PP state decay rate Γ\Gamma= 2π×\pi\times(5.75×\times106 Hz). The optical depth DD then follows the Gaussian distribution D=L×n×exp(2r2/R2)×σ×exp(2r2/R2)D=L\times n\times\exp⁡(-2r^{2}/R^{2})\times\sigma\times\exp⁡(-2r^{2}/R^{2}), where the length of the ensemble LL along the fiber axis is approximated with 5 mm Xin , the number density nn is approximated with 1.3×\times1011 cm-3 from the maximum measured DD=40, and σ\sigma=0.65×\times10-9 cm2 is the scattering cross-section of the D1 line FF = 3, mm = 0 to FF’ = 3, mm = ±\pm1 transition at the center of the dipole trap.

When the cloud size in the radial direction is assumed to have the same distribution as the light in the fiber, the simulated delayed pulse is within 10%\% of the experimental delayed pulse in terms of delay time and efficiency. The uniform profile calculation uses the peak intensity of control and probe beams, and the peak optical depth. The ground state coherence γ21\gamma_{21} is chosen as 2π×\pi\times(103 Hz) from the coherence time measurement of the microwave Ramsey interferometer before light shift cancellation of the optical lattice described in the following section.

Figure 2(b) compares the group delay of the probe pulse in the fiber with the free space scenario of uniform control beam, probe beam, and atoms density profiles. The inhomogeneous profiles of the control and the atom density in the fiber result in more than a factor of 2 reduction of the group delay and the efficiency. Figure 2(c) shows the temporal profile of the probe pulse after turning off the control beam when it is propagating inside the medium and turning on the control beam again after a storage time TT. The retrieved pulse has a FWHM duration of 160 ns, which corresponds to about 1 MHz bandwidth. We are able to store and retrieve 10%\% of the input probe pulse in the fiber.

For the delayed light measurements, we use a moving optical lattice to guide the atoms into the hollow-core fiber and stop them inside the fiber with a stationary optical lattice. This is to avoid loss and decoherence of the atoms due to their motion in the fiber. The moving optical lattice beams are formed by a pair of counter-propagating fields (100 mW each) in the fiber with a velocity vv=Δf×λ/2\Delta f\times\lambda/2=2.05 cm s-1, where Δf\Delta f=50 kHz is the frequency detuning of the two lattice fields, and λ\lambda= 821 nm is the lattice wavelength. We transport about 10%\% of the atoms into the fiber compared to the loading with just a single optical dipole beam. The dominant loss is due to heating of the atoms during the transporting process which results in escape of atoms trapped in the higher vibrational energy levels of the moving optical lattice. Figure 3(a) shows the storage efficiency of the light as a function of time. We plot the efficiency of 1550 nm light propagating in a solid-core fiber loop for comparison. Our result exceeds the performance of 1550 nm light after 600 μ\mus of delay time.

Refer to caption
Figure 4: Efficiency versus delay time with different dynamic decoupling (DD) pulse separation times. The black squares represent an average of 40 experimental runs of light delay efficiency and the error bars indicate the standard errors. The red circles are the fit of the simulation to the experimental data. The lines are used to guide the eyes. The efficiencies are normalized to the first points for each TddT_{\textrm{dd}}. The ratio of the Gaussian waists of the ensemble to the control beam in radial direction α\alpha and the frequency difference of the microwave and the EIT two-photon detuning Δ\Delta are set as free parameters to fit the data. Due to the radial distribution of the EIT beams and atoms, the Rabi frequencies of the control and the probe pulses are set with a Gaussian distribution in the radial direction in the simulation, where their maximum Rabi frequencies excite a π\pi/2 pulse with 250 ns duration. The duration of the microwave π\pi pulses is 37 μ\mus and the Rabi frequency is assumed uniform across all the atoms. The inset shows the measurement of contrast versus time using microwave Ramsey interferometer with DD sequence. The curve is an exponential fit to the experimental data.

We further improve the performance by minimizing the differential ac Stark shift of |1|1\rangle and |2|2\rangle states which is the primary decoherence source of the atomic spin coherence Xin . This is carried out by introducing an elliptical-polarized dipole beam and an external magnetic field Der to create a vector light shift to cancel the scalar light shift. Figure 3(b) shows the storage efficiency as a function of the delay time. We apply an external magnetic field of 2.08 G and optimize the efficiency by changing the ellipticity of lattice beams through quarter-wave plates. An exponential decay function is fitted to the data with the 1/ee decay time of 20 ms, which is limited by the inhomogeneous stray magnetic field of the experimental apparatus Xin . Figure 3(c) shows experimental results of the normalized efficiency of light storage in a free-falling atomic spin-wave for comparison. The storage time is shorter than for the atoms trapped in a stationary optical lattice. We attribute this to the inhomogeneity of the magnetic field and dipole potential along the trajectory of atoms that perturb the spin coherence and the relative motion of atoms, respectively.

One way to prolong the delay time is to execute the dynamic decoupling (DD) method Sagi ; Dud , a series of NN microwave population inverting π\pi pulses separated by cycling time TddT_{\textrm{dd}} is applied to the atoms during the storage period to rapidly rephase the spin coherence, where NN is an even integer. We study the DD method under our fiber delay line condition experimentally and theoretically. Figure 4 shows the delay efficiency versus time measurements with TddT_{\textrm{dd}}=1, 2.5, and 5 ms. We observe the oscillation of the efficiency over time due to the frequency difference Δ\Delta between the microwave frequency and the two-photon detuning between the control and the probe pulses. We confirm the oscillation by simulating the rotation of the Bloch vectors of the atoms under DD.

In the simulation, we set the residual inhomogeneous broadening of the dipole beam δ\delta=2π×\pi\times(50 Hz) from the data obtained in Fig. 3(b) and introduce α\alpha as the ratio of Gaussian waists of the atomic ensemble to the control beam in radial direction and set it as a free parameter together with Δ\Delta. We fit the simulation to the data and find the parameters α\alpha and Δ\Delta for each TddT_{\textrm{dd}}. The decay of the efficiency is due to the spatial inhomogeneity of the control beam propagating in the fiber, which is determined by α\alpha. Further details can be found in the Appendix.

In the Bloch sphere picture, the inhomogeneous control pulse prepares Bloch vectors on the Bloch sphere with a distribution for the DD sequence. To verify it, we use the same DD sequence on a microwave Ramsey interferometer, replacing the EIT probe and control pulses with microwave pulses which are assumed to be homogeneous across the whole atomic ensemble. We find that the 1/ee decay time τ\tau of the contrast of the pure microwave sequence is two times longer than the storage time under the same DD; see the inset of Fig. 4(b).

Light storage in the hollow-core fiber shows three orders of magnitude slower decay rate than the propagation loss of light in the solid-core fiber. Moreover, the delay time could be varied precisely by merely turning on and off the control pulse as an optical switch without modifying any hardware as in the solid-core fiber loop delay line. The low initial storage efficiency (3%\sim 3\%) in our experiment compared to the free space case is due to the inhomogeneous profile of the control beam and atom density inside the fiber. We have achieved the time-bandwidth product (TBP) to about 5×\times104, a figure of merit for quantifying the performance of the optical storage capacity. The solid-core fiber systems Zhu ; Sag ; Jin ; Sag2 have achieved TBP to about 800 with a few percent of efficiency and atomic systems Spr ; Gou ; Say ; Cor ; Bla ; Pet have achieved TBP to about 10 with about 30%\% efficiency. Our efficiency can be improved by increasing the optical depth Bla ; Hsi or using other light storage methods such as Autler-Townes splitting to ease the demand of high optical depth Sag3 . The near-infrared wavelength of our light can be converted into the telecom band in the same medium using a higher level of atomic transitions Cha ; Rad and be extended to single photons for quantum network applications Kim .

We thank Alex Kuzmich and Yi-Hsin Chen for reading the manuscript. This work is supported by the Singapore National Research Foundation under Grant No. NRFF2013-12, Nanyang Technological University under start-up grants, and Singapore Ministry of Education under Grants No. Tier 1 RG107/17.

*

Appendix A Details of simulations of the light storage efficiency under dynamical decoupling

A.1 Assumptions

We assume the following for the dynamical decoupling (DD) simulations:

  1. 1.

    The effective Rabi frequency ΩR(r)\Omega_{R}(r) of the EIT control and probe pulse for an atom at radial position rr is Gaussian distributed according the fundamental mode field profile of the hollow-core fiber as ΩR(r)=ΩRmaxe2r2/R2\Omega_{R}(r)=\Omega^{max}_{R}e^{-2r^{2}/R^{2}}, where R=22μR=22~{}\mum is the 1/e21/e^{2} mode field radius and ΩRmax=π/(500ns)\Omega^{max}_{R}=\pi/(500~{}\mathrm{ns}).

  2. 2.

    The EIT fields are two-photon resonant on the unperturbed atomic transition.

  3. 3.

    The Rabi frequency ΩM\Omega_{M} of the microwaves is uniform across the ensemble with a π\pi-pulse time of 37 μ\mus.

  4. 4.

    The microwave detuning from the unperturbed transition frequency is Δ\Delta. Δ\Delta is used as one of the fit parameters to the data.

  5. 5.

    The inhomogeneous broadening δ(r)\delta(r) caused by the residual optical lattice differential ac Stark shift is also Gaussian distributed as δ(r)=δmaxe2r2/R2\delta(r)=\delta_{max}e^{-2r^{2}/R^{2}}, where δmax=2π×(50Hz)\delta_{max}=2\pi\times(50~{}\mathrm{Hz}) from the data in Fig. 3(b) of the main paper.

  6. 6.

    Initially, at t=0t=0, the phase of the microwaves is the same as the EIT fields, i.e., over time tt, the microwave phase relative to an atom at position rr is (Δ+δ(r))t(\Delta+\delta(r))t.

  7. 7.

    The atoms are Gaussian distributed according to ρ(r)=exp(r2/2Ratom2)/Ratom2\rho(r)=\mathrm{exp}\left(-r^{2}/2R^{2}_{atom}\right)/R_{atom}^{2} where RatomR_{atom} is the rms radius of the atomic ensemble. The probability distribution is normalized: 0ρ(r)r𝑑r=1\int_{0}^{\infty}\rho(r)rdr=1. In the simulation, the atom positions are fixed for the duration of the dynamical decoupling sequence. The ratio α\alpha of the ensemble size to the mode field radius, i.e., αRatom/R\alpha\equiv R_{atom}/R, is used as the other fit parameter to the data.

A.2 Simulation on a single atom

The goal is to determine the collective state of the atoms after the EIT pulse and microwave dynamical decoupling pulses with the atoms evolving phase in the optical lattice.

Beginning with an atom at position rr in the ground state, represented by the Bloch vector 𝐮0=(0,0,1)\mathbf{u}_{0}=(0,0,-1), we model the EIT pulse as a two-photon Raman π/2\pi/2-pulse acting on the Bloch vector. The rotation matrix representing the π/2\pi/2 pulse takes into account the inhomogeneous Rabi frequency ΩR(r)\Omega_{R}(r) and detuning δ(r)\delta(r). In the Bloch sphere picture, the inhomogeneous EIT pulse prepares Bloch vectors on the Bloch sphere with a distribution for the DD sequence.

After the EIT pulse, the atom evolves phase δ(r)Tdd\delta(r)T_{\mathrm{dd}} in the optical lattice for DD cycling time TddT_{\mathrm{dd}}. In the Bloch sphere picture, the Bloch vector rotates about the zz-axis by angle δ(r)Tdd\delta(r)T_{\mathrm{dd}}. A microwave π\pi-pulse is then applied, with Rabi frequency ΩM\Omega_{M}, detuning Δ+δ(r)\Delta+\delta(r), and phase (Δ+δ(r))t(\Delta+\delta(r))t, where tt is the time at which the π\pi-pulse is applied. Again, the radial dependence encodes the inhomogeneity of the ensemble with respect to the microwaves. The DD sequence repeats the elementary sequence [free-evolution TddT_{\mathrm{dd}}π\pi-pulse] for an even number of times. The Bloch vector of the atom 𝐉=(Jx,Jy,Jz)\mathbf{J}=(J_{x},J_{y},J_{z}) is numerically determined after each repeating unit of the DD sequence.

A.3 Weighting over ensemble

The ensemble has a radial distribution across the fiber mode profile leading to a distribution of Rabi frequency ΩR(r)\Omega_{R}(r) and detuning δ(r)\delta(r). As there is a one-to-one mapping between the radial position rr and the normalized fiber mode intensity I(r)=e2r2/R2I(r)=e^{-2r^{2}/R^{2}}, the single atom simulation 𝐉\mathbf{J} is actually a function 𝐉(I)\mathbf{J}(I) of the normalized mode intensity II.

To simulate the effect of this inhomogeneous distribution, we weight the simulated single atom Bloch vectors 𝐉(I)\mathbf{J}(I) by the intensity distribution ρ(I)\rho(I) that the ensemble samples as in the following

𝐉01𝐉(I)ρ(I)𝑑I.\langle\mathbf{J}\rangle\equiv\int^{1}_{0}\mathbf{J}(I)\rho(I)dI. (3)

This defines the collective Bloch vector 𝐉\langle\mathbf{J}\rangle that we use to calculate the retrieval efficiency. Using the Gaussian atom distribution provided in the assumptions, it can be shown that the equivalent intensity distribution is given by

ρ(I)=14α2I14α21.\rho(I)=\frac{1}{4\alpha^{2}}I^{\frac{1}{4\alpha^{2}}-1}. (4)

For the numerical simulations, we approximate the integral using the trapezium rule.

A.4 Retrieval efficiency

The retrieval efficiency is determined by the collective coherence of the ensemble, which corresponds to the component of the collective Bloch vector in the equatorial plane J2Jx2+Jy2J_{\perp}^{2}\equiv\langle J_{x}\rangle^{2}+\langle J_{y}\rangle^{2}, where Jx\langle J_{x}\rangle is the xx-component of the collective Bloch vector 𝐉\langle\mathbf{J}\rangle; similarly for Jy\langle J_{y}\rangle. The normalized efficiency is then (J/J,0)2(J_{\perp}/J_{\perp,0})^{2} where J,0J_{\perp,0} is the collective coherence of the ensemble right after the EIT storage process.

A.5 Fitting simulation to data

To fit the simulation to the data, we use the method of least squares. Specifically, we define the sum of squares deviation between the simulation and experimental data points as

SS=i(simulation(Δ,α)idatai)2errori2\mathrm{SS}=\sum_{i}\frac{(\mathrm{simulation(\Delta,\alpha)}_{i}-\mathrm{data}_{i})^{2}}{\mathrm{error}_{i}^{2}} (5)

where errori is the standard error of the ii-th experimental data point. We then numerically minimize SS\mathrm{SS} over the free parameters Δ\Delta and α\alpha to obtain fits to the data sets.

References

  • (1) L. Thévenaz, Slow and fast light in optical fibres. Nat. Photonics 2, 474 (2008).
  • (2) G. M. Gehring, R. W. Boyd, A. L. Gaeta, D. J. Gauthier, and A. E. Willner, Fiber-based slow-light technologies. J. Lightwave Technol. 26, 3752 (2008).
  • (3) Z. Zhu, D. J. Gauthier, and R. W. Boyd, Stored light in an optical fiber via stimulated Brillouin scattering. Science 318, 1748 (2007).
  • (4) E. Saglamyurek, J. Jin, V. B. Verma, M. D. Shaw, F. Marsili, S. W. Nam, D. Oblak, and W. Tittel, Quantum storage of entangled telecom wavelength photons in an erbium-doped optical fibre. Nat. Photonics 9, 83 (2015).
  • (5) J. Jin, E. Saglamyurek, M. G. Puigibert, V. Verma, F. Marsili, S. W. Nam, D. Oblak, and W. Tittel, Telecom-wavelength atomic quantum memory in optical fiber for heralded polarization qubits. Phys. Rev. Lett. 115, 140501 (2015).
  • (6) E. Saglamyurek, M. G. Puigibert, Q. Zhou, L. Giner, F. Marsili, V. B. Verma, S. W. Nam, L. Oesterling, D. Nippa, D. Oblak, and W. Tittel, A multiplexed light-matter interface for fibre-based quantum networks. Nat. Commun. 7, 11202 (2016).
  • (7) M. R. Sprague, P. S. Michelberger, T. E. F. Champion, D. G. England, J. Nunn, X.-M. Jin, W. S. Kolthammer, A. Abdolvand, P. St. J. Russell, and I. A. Walmsley, Broadband single-photon-level memory in a hollow-core photonic crystal fiber. Nat. Photonics 8, 287 (2014).
  • (8) B. Gouraud, D. Maxein, A. Nicolas, O. Morin, and J. Laurat, Demonstration of a memory for tightly guided light in an optical nanofiber. Phys. Rev. Lett. 114, 180503 (2015).
  • (9) C. Sayrin, C. Clausen, B. Albrecht, P. Schneeweiss, A. Rauschenbeutel, Storage of fiber-guided light in a nanofiber trapped ensemble of cold atoms. Optica 2, 353 (2015).
  • (10) N. V. Corzo, J. Raskop, A. Chandra, A. S. Sheremet, B. Gouraud, and J. Laurat, Waveguide-coupled single collective excitation of atomic arrays. Nature 566, 359 (2019).
  • (11) F. Blatt, L. S. Simeonov, T. Halfmann, T. Peters, Stationary light pulses and narrowband light storage in a laser-cooled ensemble loaded into a hollow-core fiber. Phys. Rev. A 94, 043833 (2016).
  • (12) T. Peters, T.-P. Wang, A. Neumann, L. S. Simeonov, and T. Halfmann, Stopped and stationary light at the single-photon level inside a hollow-core fiber. Opt. Express 28, 5340 (2020).
  • (13) B. Debord, M. Alharbi, T. Bradley, C. Fourcade-Dutin, Y. Y. Wang, L. Vincetti, F. Gérôme, and, F, Benabid, Hypocycloid-shaped hollow-core photonic crystal fiber Part I: Arc curvature effect on confinement loss. Opt. Express 21, 28597 (2013).
  • (14) M. Alharbi, T. Bradley, B. Debord, C. Fourcade-Dutin, D. Ghosh, L. Vincetti, F. Gérôme, and, F, Benabid, Hypocycloid-shaped hollow-core photonic crystal fiber Part II: Cladding effect on confinement and bend loss. Opt. Express 21, 28609 (2013).
  • (15) M. Fleischhauer, A. Imamoglu, and J. P. Marangos, Electromagnetically induced transparency: optics in coherent media. Rev. Mod. Phys. 77, 633 (2005).
  • (16) Y.-F. Hsiao, P.-J. Tsai, H.-S. Chen, S.-X. Lin, C.-C. Hung, C.-H. Lee, Y.-H. Chen, Y.-F. Chen, I. A. Yu, and Y.-C. Chen, Highly efficient coherent optical memory based on electromagnetically induced transparency. Phys. Rev. Lett. 120, 183602 (2018).
  • (17) M. Xin, W. S. Leong, Z. Chen, and S.-Y. Lan, Transporting long-lived quantum spin coherence in a photonic crystal fiber. Phys. Rev. Lett. 122, 163901 (2019).
  • (18) A. Derevianko, Theory of magic optical traps for Zeeman insensitive clock transitions in alkali-metal atoms. Phys. Rev. A 81, 051606(R) (2010).
  • (19) Y. Sagi, I. Almog, and N. Davidson, Process tomography of dynamical decoupling in a dense cold atomic ensemble. Phys. Rev. Lett. 105, 053201 (2010).
  • (20) Y. O. Dudin, L. Li, and A. Kuzmich, Light storage on the time scale of a minute. Phys. Rev. A 87, 031801(R) (2013).
  • (21) E. Saglamyurek, T. Hrushevskyi, A. Rastogi, K. Heshami, and L. J. LeBlanc, Coherent storage and manipulation of broadband photons via dynamically controlled Autler–Townes splitting. Nat. Photonics 12, 774 (2018).
  • (22) T. Chanelière, D. N. Matsukevich, S. D. Jenkins, T. A. B., Kennedy, M. S. Chapman, and A. Kuzmich, Quantum telecommunication based on atomic cascade transitions. Phys. Rev. Lett. 96, 093604 (2006).
  • (23) A. G. Radnaev, Y. O. Dudin, R. Zhao, H. H. Jen, S. D., Jenkins, A. Kuzmich, and T. A. B. Kennedy, A quantum memory with telecom wavelength conversion. Nat. Physics 6, 894 (2010).
  • (24) H. J. Kimble, The quantum internet. Nature 453, 1023–1030 (2008).