Long Light Storage Time in an Optical Fiber
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.710-5 dB s-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 hereOptical 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.


A cold 85Rb atomic ensemble is prepared by a magneto-optical trap (MOT) 5 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 K. The wavelength of the optical dipole beam is 821 nm with 250 mW power, and the measured optical depth () is about 40 by probing all the Zeeman states of the 85Rb D2 line =3 to ’=4 transition. When the atoms are in the fiber, we optically pump them into the =2, =0 state by a -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 configuration Fle operating on the hyperfine clock states of 85Rb for long coherence times. The probe beam is resonant on 1==2, =0 to 3=’=3, =1 and the control beam is resonant on 2==3 =0 to 3. 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 ()=exp(-4(ln2)/), where is the maximum intensity, =304 ns is the FWHM duration, and is the temporal coordinate. A Gaussian function is fitted to the transmitted probe pulse to determine the group delay =174 ns with 24 efficiency of transmission.

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 =31.5 m Hsi . The intensity of the probe pulse after transmitting through the atomic ensemble is obtained as
(1) |
where the factor
(2) |
and the group delay . The normalized Rabi frequency of the probe pulse is assumed to have a Gaussian distribution =(-/), where is the radial coordinate of the fibre core and =22 m is the 1/ mode field radius. The ground state coherence , which is mainly due to the inhomogeneous differential ac Stark shift of the control and dipole beams, is fixed at 2(2105 Hz) based on the measurement of EIT linewidth and ground and excited states coherence is fixed at half of the 85Rb D1 line state decay rate = 2(5.75106 Hz). The optical depth then follows the Gaussian distribution , where the length of the ensemble along the fiber axis is approximated with 5 mm Xin , the number density is approximated with 1.31011 cm-3 from the maximum measured =40, and =0.6510-9 cm2 is the scattering cross-section of the D1 line = 3, = 0 to ’ = 3, = 1 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 is chosen as 2(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 . 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 ==2.05 cm s-1, where =50 kHz is the frequency detuning of the two lattice fields, and = 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 s of delay time.

We further improve the performance by minimizing the differential ac Stark shift of and 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/ 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 microwave population inverting pulses separated by cycling time is applied to the atoms during the storage period to rapidly rephase the spin coherence, where 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 =1, 2.5, and 5 ms. We observe the oscillation of the efficiency over time due to the frequency difference 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 =2(50 Hz) from the data obtained in Fig. 3(b) and introduce 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 . We fit the simulation to the data and find the parameters and for each . The decay of the efficiency is due to the spatial inhomogeneity of the control beam propagating in the fiber, which is determined by . 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/ decay time 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 () 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 5104, 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 AutlerTownes 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.
The effective Rabi frequency of the EIT control and probe pulse for an atom at radial position is Gaussian distributed according the fundamental mode field profile of the hollow-core fiber as , where m is the mode field radius and .
-
2.
The EIT fields are two-photon resonant on the unperturbed atomic transition.
-
3.
The Rabi frequency of the microwaves is uniform across the ensemble with a -pulse time of 37 s.
-
4.
The microwave detuning from the unperturbed transition frequency is . is used as one of the fit parameters to the data.
-
5.
The inhomogeneous broadening caused by the residual optical lattice differential ac Stark shift is also Gaussian distributed as , where from the data in Fig. 3(b) of the main paper.
-
6.
Initially, at , the phase of the microwaves is the same as the EIT fields, i.e., over time , the microwave phase relative to an atom at position is .
-
7.
The atoms are Gaussian distributed according to where is the rms radius of the atomic ensemble. The probability distribution is normalized: . In the simulation, the atom positions are fixed for the duration of the dynamical decoupling sequence. The ratio of the ensemble size to the mode field radius, i.e., , 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 in the ground state, represented by the Bloch vector , we model the EIT pulse as a two-photon Raman -pulse acting on the Bloch vector. The rotation matrix representing the pulse takes into account the inhomogeneous Rabi frequency and detuning . 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 in the optical lattice for DD cycling time . In the Bloch sphere picture, the Bloch vector rotates about the -axis by angle . A microwave -pulse is then applied, with Rabi frequency , detuning , and phase , where is the time at which the -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 – -pulse] for an even number of times. The Bloch vector of the atom 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 and detuning . As there is a one-to-one mapping between the radial position and the normalized fiber mode intensity , the single atom simulation is actually a function of the normalized mode intensity .
To simulate the effect of this inhomogeneous distribution, we weight the simulated single atom Bloch vectors by the intensity distribution that the ensemble samples as in the following
(3) |
This defines the collective Bloch vector 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
(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 , where is the -component of the collective Bloch vector ; similarly for . The normalized efficiency is then where 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
(5) |
where errori is the standard error of the -th experimental data point. We then numerically minimize over the free parameters and 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).