Magnetic-field-induced deflection of nonlocal light bullets in a Rydberg atomic gas

Xiujia Dong; Yao Ding; Zhengyang Bai; Guoxiang Huang

doi:10.3788/COL202220.041902

Abstract

Light bullets (LBs) are localized nonlinear waves propagating in high spatial dimensions. Finding stable LBs and realizing their control are desirable due to the interesting physics and potential applications. Here, we show that nonlocal LBs generated in a cold Rydberg atomic gas via the balance among the dispersion, diffraction, and giant nonlocal Kerr nonlinearity contributed by long-range Rydberg-Rydberg interaction can be actively manipulated by using a weak gradient magnetic field. Nonlocal LBs are generated by a balance among dispersion, diffraction, and large nonlocal Kerr nonlinearities contributed by long-range Rydberg-Rydberg interactions. Here, we find that active manipulation can be achieved by weak gradient magnetic fields in cold Rydberg atomic gases. Especially, the LBs may undergo significant Stern–Gerlach deflections, and their motion trajectories can be controlled by adjusting the magnetic-field gradient. The results reported here are helpful not only for understanding unique properties of LBs in nonlocal optical media but also for finding ways for precision measurements of magnetic fields.

Keywords

nonlinear optics optical soliton Rydberg atom

1. Introduction

In recent years, much attention has been paid to the investigation on electromagnetically induced transparency (EIT) in cold Rydberg atomic gases^[1–26]. This is rooted in the fact that atomic Rydberg states have long coherent lifetimes and strong long-range interaction (called Rydberg-Rydberg interaction) between remote atoms. The Rydberg-Rydberg interaction makes the atomic gases be nonlocal optical media, and they can be effectively mapped to strong photon-photon interactions via EIT. As a result, strong nonlinearities at very low light intensity can be realized^[5,6,14], which opens up a new avenue to study nonlinear and quantum optics and realize novel photon devices, such as single-photon switches^[27–31], optical transistors^[4,32], photon memories^[10,11], and single-photon sources^[33].

Light bullets (LBs)^[34] are solitary nonlinear waves localized in $m$ spatial dimensions and one time dimension [ $(m + 1) D$ ; $m = 1, 2, 3$ ]. In recent years, the study of LBs has attracted intensive theoretical and experimental interests^[35] because of their rich nonlinear physics and technological applications^[36,37]. However, the generation of stable high-dimensional LBs is a topic not solved for a long time. It has been shown recently that stable $(3 + 1) D$ nonlocal LBs can be realized in Rydberg atomic gases; such LBs have extremely low generation power and ultra slow propagation velocity^[17]. Different from the non-interaction system, the central element is the co-existence of giant nonlocal and local optical Kerr nonlinearities. The former features a fast (sub-microsecond) response^[38], which is complemented by the latter, whose response is relatively slow (in the order of microseconds). In conjunction with tunable dispersion and diffraction, this allows us to precisely control dynamics of LBs.

In this work, we propose a scheme to realize the active control of the nonlocal LBs in a Rydberg atomic gas. We show that the $(3 + 1) D$ LBs generated in such a system via EIT can be manipulated by using a gradient magnetic field. In particular, the LBs can undergo significant Stern–Gerlach deflections even when the magnetic-field gradient is weak, and their motion trajectories can be adjusted through the changing of the magnetic-field gradient. Our work contributes to the efforts for understanding the unique properties and realizing the active controls of high-dimensional LBs and also for finding new techniques for precision measurements of magnetic fields.

Sign up for Chinese Optics Letters TOC. Get the latest issue of Chinese Optics Letters delivered right to you！Sign up now

2. Model

The system under study is a cold three-level atomic gas working with a Rydberg-EIT scheme, shown in Fig. 1(a). Here, the levels $| 1 〉$ , $| 2 〉$ , and $| 3 〉$ are ground, intermediate, and Rydberg states, respectively; a weak, pulsed (with time duration $τ_{0}$ ) probe laser field (center angular frequency $ω_{p}$ , center wavenumber $k_{p} = ω_{p} / c$ , half-Rabi frequency $Ω_{p}$ ) couples to the transition between $| 1 〉$ and $| 2 〉$ ; a strong, continuous-wave control laser field (angular frequency $ω_{c}$ , wavenumber $k_{c} = ω_{c} / c$ , half-Rabi frequency $Ω_{c}$ ) couples to the transition between $| 2 〉$ and $| 3 〉$ . The total electric field of the system reads $E = E_{p} + E_{c}$ , with $E_{α} = e_{α} E_{α} \exp [i (k_{α} \cdot r - ω_{α} t)] + h.c.$ . Here, with $r = (x, y, z)$ , $e_{α}$ and $E_{α}$ are the unit polarization vector and field amplitude for the $α$ field ( $α = p, c$ ), respectively. To suppress Doppler effect, the probe and the control fields are assumed to counter-propagate along the $z$ direction. Figure 1(b) shows the geometry of the system.

Figure 1.(a) Excitation scheme of the Rydberg EIT. |1〉, |2〉, and |3〉 are, respectively, the ground, intermediate, and Rydberg states; Ω_p (Ω_c) is the half-Rabi frequency of the probe (control) laser field; Γ₁₂ (∼MHz) and Γ₂₃ (∼kHz) are, respectively, decay rates from |2〉 to |1〉 and |3〉 to |2〉; Δ₂ = ω_p − (ω₂ − ω₁) and Δ₃ = ω_p + ω_c − (ω_c − ω₁) are, respectively, the one- and two-photon detunings. ℏV(r − r^′) is the vdW interaction between the two atoms in Rydberg states, respectively, located at r and r^′. (b) Geometry of the system. The probe and control fields counter-propagate in the Rydberg atomic gas. (c) Normalized χ_p,2⁽³⁾ (i.e., the coefficient of nonlocal Kerr nonlinearity) as a function of coordinate x, with the solid black (dashed red) line representing its real part Re(χ_p,2⁽³⁾) [imaginary part Im(χ_p,2⁽³⁾)] for coordinate y = 0 (see text for more details). Evolution of a (d1) nonlocal LB and (d2) vortex in the system.

For realizing the active control on the LBs, a weak gradient magnetic field is assumed to act on the atomic gas, with the form $B (x, y) = \hat{z} B (r_{⊥}) = \hat{z} (B_{1} x + B_{2} y),$ (1)where $r_{⊥} = (x, y,0)$ , $\hat{z} = (0, 0, 1)$ is the unit vector in the $z$ direction, $B_{1}$ and $B_{2}$ characterize the gradients of the magnetic field in the $x - y$ plane. Due to the presence of the magnetic field, each level of the atoms is split into a series of Zeeman sub-levels with energy $Δ E_{α}^{Zeeman} = μ_{B} g_{F}^{α} m_{F}^{α} B$ , where $μ_{B}$ , $g_{F}^{α}$ , and $m_{F}^{α}$ are the Bohr magneton, gyromagnetic factor, and magnetic quantum number of level $| α 〉$ , respectively. As a result, the one- and two-photon detunings $Δ_{2}$ and $Δ_{3}$ are changed into $Δ_{2} (r_{⊥}) = (ω_{p} - ω_{c} - ω_{21}) + μ_{21} B (r_{⊥})$ and $Δ_{3} (r_{⊥}) = (ω_{p} - ω_{31}) + μ_{31} B (r_{⊥})$ , with $μ_{α β} = μ_{B} (g_{F}^{α} m_{F}^{α} - g_{F}^{β} m_{F}^{β}) / ℏ$ .

Under electric-dipole and rotating-wave approximations, the Hamiltonian of the atomic gas including the Rydberg-Rydberg interaction is given by $\hat{H} = N_{a} \int d^{3} r \hat{H} (r, t)$ , with $\hat{H} (r, t)$ the Hamiltonian density, given by $\hat{H} = \sum_{α = 1}^{3} ℏ ω_{α} {\hat{S}}_{α α} (r, t) - ℏ [Ω_{p} {\hat{S}}_{12} (r, t) + Ω_{c} {\hat{S}}_{23} (r, t) + h.c.] + N_{a} \int d^{3} r^{'} {\hat{S}}_{33} (r^{'}, t) ℏ V (r^{'} - r) {\hat{S}}_{33} (r, t) .$ (2)

Here, $N_{a}$ is atomic density, $Ω_{p} \equiv (e_{p} \cdot p_{21}) E_{p} / ℏ$ and $Ω_{c} \equiv (e_{c} \cdot p_{32}) E_{c} / ℏ$ are, respectively, the half-Rabi frequencies of the probe and control fields (with $p_{α β}$ the electric-dipole matrix element associated with the transition from $| β 〉$ to $| α 〉$ ), ${\hat{S}}_{α β} \equiv | β 〉 〈 α | e^{i [(k_{β} - k_{α}) \cdot r - (ω_{β} - ω_{α} + Δ_{β} - Δ_{α}) t]}$ are atomic transition operators ( $α, β = 1, 2, 3$ ), and $ℏ V (r - r^{'}) \equiv - ℏ C_{6} / {| r - r^{'} |}^{6}$ is the van der Waals (vdW) interaction potential (with $C_{6}$ the dispersion coefficient) between the Rydberg atoms located, respectively, at the positions $r$ and $r^{'}$ ^[17].

The dynamics of the atoms is controlled by the Heisenberg equation of motion for the atomic operators ${\hat{S}}_{α β}$ , i.e., $i ℏ \partial {\hat{S}}_{α β} / \partial t = [\hat{H}, {\hat{S}}_{α β}]$ . Taking expectation values on both sides of this equation, we obtain the optical Bloch equation involving one- and two-body reduced density matrices (DMs), which can be cast into the form $\frac{\partial ρ}{\partial t} = - \frac{i}{ℏ} [\hat{H}, ρ] - Γ [ρ] .$ (3)

Here, $ρ (r, t)$ is the reduced DM in the single-particle basis ${| 1 〉, 2 〉,3 〉}$ , with the matrix elements defined by $ρ_{α β} (r, t) \equiv 〈 {\hat{S}}_{α β} (r, t) 〉$ ; $Γ$ is a $3 \times 3$ relaxation matrix describing the spontaneous emission and dephasing. Due to the existence of the Rydberg-Rydberg interaction, two-body reduced DM [i.e., $ρ_{twobody}$ with DM elements $ρ_{α β, μ ν} (r^{'}, r, t)$ ] is involved in this equation, denoting the contribution from the Rydberg-Rydberg interaction. The explicit expression of Eq. (3) is presented in Supplementary Material^[39].

The dynamics of the probe field is described by the Maxwell equation, which, under the slowly varying envelope approximation, reads $i (\frac{\partial}{\partial z} + \frac{1}{c} \frac{\partial}{\partial t}) Ω_{p} + \frac{c}{2 ω_{p}} \nabla_{⊥}^{2} Ω_{p} + κ_{12} ρ_{21} = 0,$ (4)where $\nabla_{⊥}^{2} = \partial^{2} / \partial x^{2} + \partial^{2} / \partial y^{2}$ , $κ_{12} = N_{a} ω_{p} | p_{12} |^{2} / (2 ε_{0} c ℏ)$ , with $c$ as the light speed in vacuum. For a relatively weak probe field, the population in atomic levels changes not much when the probe field is applied to the system, and hence a perturbation expansion beyond mean-field approximation for many-atom correlations can be employed to solve the Bloch equation, Eq. (3)^[6,14,17]. The expression of the nonlinear optical susceptibility of the exact probe field to the third order of the perturbation expansion is given by $χ_{p} ≃ χ_{p}^{(1)} + χ_{p,1}^{(3)} | E_{p} |^{2} + \int d r_{⊥}^{'} χ_{p, 2}^{(3)} (r_{⊥} - r_{⊥}^{'}) {| E_{p} (r_{⊥}^{'}) |}^{2},$ (5)where $χ_{p}^{(1)} = N_{a} | p_{12} |^{2} a_{21}^{(1)} / (ε_{0} ℏ)$ , $χ_{p,1}^{(3)} = N_{a} | p_{12} |^{4} a_{21, 1}^{(3)} / (ε_{0} ℏ^{3})$ , and $χ_{p, 2}^{(3)} = \int d z N_{a}^{2} | p_{12} |^{4} a_{21, 2}^{(3)} / (ε_{0} ℏ^{3})$ ^[14,17]. In the above expression, the third term is the nonlocal Kerr nonlinear susceptibility, contributed from the Rydberg-Rydberg interaction, while the second term is the local Kerr nonlinear susceptibility, contributed from the non-zero two-photon detuning (i.e., $Δ_{3} \neq 0$ ).

To be concrete in the following calculations, we choose strontium atoms ( ${}^{88}{Sr}$ ), although our approach is valid for other Rydberg atomic gases. The energy levels shown in Fig. 1(a) are selected as $| 1 〉 = | 5 s^{2} {{}^{1}S}_{0} 〉$ , $| 2 〉 = | 5 s 5 p^{1} P_{1} 〉$ , $| 3 〉 = | 5 s n s^{1} S_{0} 〉$ , with $n$ the principal quantum number^[39]. The spontaneous emission rates of the atoms are given by $Γ_{12} = 2 π \times 16 MHz$ , $Γ_{23} = 2 π \times 16.7 kHz$ , so one has $γ_{21} = Γ_{12} / 2$ , $γ_{31} = Γ_{23} / 2$ , and $γ_{32} = (Γ_{12} + Γ_{23}) / 2$ . For this choice, the vdW interaction in ${{}^{1}S}_{0}$ states is isotropically attractive ( $C_{6} > 0$ ), which is important to realize self-focusing nonlocal Kerr nonlinearity.

The result shown in Fig. 1(c) is $χ_{p, 2}^{(3)}$ as a function of $x$ for $y = 0$ . The real part [ $Re (χ_{p, 2}^{(3)})$ ] and imaginary part [ $Im (χ_{p, 2}^{(3)})$ ] are plotted on the condition that $Δ_{2} ≫ Γ_{12}$ (i.e., the system works in the dispersive nonlinearity regime) by the solid black line and the dashed red line, respectively. We see that $Im (χ_{p, 2}^{(3)})$ is much smaller than $Re (χ_{p, 2}^{(3)})$ , which means that the optical absorption of the probe field is negligible, resulting from the EIT effect and the condition of large one-photon detuning; moreover, $Re (χ_{p, 2}^{(3)})$ is an attractive potential well, and there is a saturation near $x = 0$ , which is due to the Rydberg blockade effect (with blockade radius $\sim 7 μm$ ) that suppresses the excitation of atoms to the Rydberg state and hence causes the nonlinear kernel $χ_{p, 2}^{(3)}$ to saturate to a finite value. By virtue of the strong Rydberg-Rydberg interaction, the nonlocal optical nonlinearities can reach $\int d r_{⊥}^{'} χ_{p, 2}^{(3)} (r_{⊥}^{'}) \sim 10^{- 8} m^{2} V^{- 2}$ , which are many orders of magnitude larger than that of conventional EIT systems^[6,14].

3. (3 + 1)D Nonlinear Envelope Equation

Our main aim is to implement an active control of LBs in the system. To make the related physical mechanism transparent, we first derive the equation describing the nonlinear evolution of the probe-field envelope. For a modulated plane-wave of the probe field, we assume $Ω_{p} \sim F \exp [i (K z - ω t)]$ ^[40]. The equation for the envelope $F$ in the presence of the magnetic field can be derived by means of the multiple-scale perturbation method, similar to that carried out in Ref. [41]. We obtain $i \frac{\partial u}{\partial s} = - (\frac{\partial^{2}}{\partial ξ^{2}} + \frac{\partial^{2}}{\partial η^{2}}) u - g_{d} \frac{\partial^{2} u}{\partial σ^{2}} - V_{m} (ξ, η) u + [W_{1} {| u |}^{2} + \int d {\tilde{r}}_{⊥} W_{2} ({\tilde{r}}_{⊥} - {\tilde{r}}_{⊥}^{'}) | u ({\tilde{r}}_{⊥}^{'}, s {) |}^{2}] u + i d_{0} u .$ (6)

Here, $s = z / (2 L_{diff})$ ; ${\tilde{r}}_{⊥} = (ξ, η) = (x, y) / R_{0}$ ; $σ \equiv (t - z / V_{g}) / τ_{0}$ [ $V_{g} \equiv {(\partial K / \partial ω)}^{- 1}$ is group velocity, with $K \equiv K (ω)$ the linear dispersion relation]; $u = (F / U_{0}) \exp (- α_{0} z)$ , with $α_{0} \equiv Im (K)$ a decay constant; $g_{d} \equiv - L_{diff} K_{2} / τ_{0}^{2}$ , $W_{1} \equiv - b k_{p}^{2} R_{0}^{2} U_{0}^{2} χ_{p, 1}^{(3)}$ , $W_{2} \equiv - b k_{p}^{2} R_{0}^{4} U_{0}^{2} χ_{p, 2}^{(3)}$ , and $d_{0} \equiv - 2 L_{diff} / L_{A}$ are dimensionless coefficients of dispersion, local Kerr nonlinearity, nonlocal Kerr nonlinearity, and absorption, respectively. In these coefficients, $b = ℏ / {| p_{12} |}^{2}$ , $K_{2} \equiv \partial^{2} K / \partial ω^{2}$ describes group-velocity dispersion, $L_{diff} \equiv k_{p} R_{0}^{2}$ and $L_{A} \equiv 1 / α_{0}$ are, respectively, the typical diffraction and absorption lengths, and $U_{0}$ and $R_{0}$ are, respectively, the typical half-Rabi frequency and transverse size of the probe field. Since we are interested only in the dispersive nonlinearity regime of the system, where the $L_{A}$ is much larger than the other typical lengths, and hence $d_{0}$ is very small, the imaginary parts of the coefficients in Eq. (6) are negligible.

In Eq. (6), $V_{m} \equiv - k_{p}^{2} R_{0}^{2} χ_{p}^{(1)}$ is a dimensionless linear potential contributed by the gradient magnetic field. It has the form $V_{m} (ξ, η) = V_{1} ξ + V_{2} η$ , with $V_{1} = \frac{κ_{12} R_{0} L_{diff} [(ω + Δ_{3}) Δ_{3} μ_{21} + {| Ω_{c} |}^{2} μ_{31}]}{{[(ω + Δ_{2}) (ω + Δ_{3}) - {| Ω_{c} |}^{2}]}^{2}} B_{1},$ (7a) $V_{2} = \frac{κ_{12} R_{0} L_{diff} [(ω + Δ_{3}) Δ_{3} μ_{21} + {| Ω_{c} |}^{2} μ_{31}]}{{[(ω + Δ_{2}) (ω + Δ_{3}) - {| Ω_{c} |}^{2}]}^{2}} B_{2} .$ (7b)

We then consider the formation of LBs when the gradient magnetic field is absent (i.e., $V_{m} = 0$ ). In this case, stable $(3 + 1) D$ LBs and vortices can form, with the result by a numerical simulation shown in Fig. 1(d). From the figure, we see that the $(3 + 1) D$ LB (upper part) and vortex (lower part) relax to self-cleaned forms quite close to the unperturbed ones^[42], and their shapes undergo no apparent change during propagation. The physical parameters used in the simulation are chosen as $Δ_{2} = - 15 Γ_{12}$ , $Δ_{3} = - 0.02 Γ_{12}$ , $R_{0} = 10 μm$ , $τ_{0} = 9 \times 10^{- 7} s$ , $N_{a} = 3 \times 10^{10} {cm}^{- 3}$ , $U_{0} = 0.3 Γ_{12}$ , and $C_{6} ≃ 2 π \times 81.6 MHz {μm}^{6}$ (for the principal quantum number $n = 60$ ). With these parameters, we obtain $L_{diff} = 1.36 mm$ , $L_{A} = 907 mm$ , $R_{b} = 6 μm$ , $g_{d} = 0.134$ , and $d_{0} = - 0.03$ . Such an LB can form in a very short distance and generate extremely low light power ( $\sim 1 nW$ ), which is due to the giant Kerr nonlinearity (contributed by the Rydberg-Rydberg interaction) and the ultraslow propagation velocity of the probe pulse ( $\sim 2.3 \times 10^{- 6} c$ , contributed from the EIT effect).

4. Manipulation of LBs

We now turn to investigate what will happen for a nonlocal LB when an external gradient magnetic field is present. As a first step, we consider Eq. (6) in the absence of the Kerr nonlinearity (i.e., $W_{1} = W_{2} = 0$ ). Using the transformation $u = u^{'} \exp [i (V_{1} ξ^{'} + V_{2} η^{'} + V_{1}^{2} s^{2} / 3 + V_{2}^{2} s^{2} / 3) s]$ , with $ξ^{'} = ξ - V_{1} s^{2} / 2$ and $η^{'} = η - V_{2} s^{2} / 2$ , Eq. (6) is converted into the form $i \partial u^{'} / \partial s = - (1 / 2) (\partial^{2} / \partial ξ^{' 2} + \partial^{2} / \partial η^{' 2} + g_{d} \partial^{2} / \partial τ^{2}) u^{'}$ . It is easy to obtain the expression of the central position of the probe pulse in the $ξ - η$ plane, which is given by $(ξ, η) = (V_{1} s^{2} / 2, V_{2} s^{2} / 2)$ . Returning to the original variables, the central position of the pulse reads $(x, y) = κ_{12} \frac{[(ω + Δ_{3}) Δ_{3} μ_{21} + {| Ω_{c} |}^{2} μ_{31}]}{{[(ω + Δ_{2}) (ω + Δ_{3}) - {| Ω_{c} |}^{2}]}^{2}} \frac{R_{0}^{2}}{L_{diff}} z^{2} (B_{1}, B_{2}) .$ (8)

We see that, due to the presence of the magnetic field, the motion of the linear wave is changed, and its trajectory in the $x - y$ plane has a deflection with a quadratic dependence on the propagation coordinate $z$ ; moreover, the trajectory can be controlled by tuning the gradient of the magnetic field, i.e., by manipulating the parameters $B_{1}$ and $B_{2}$ .

In the presence of the Kerr nonlinearities, it is hard to get an exact expression for the motion trajectory of the probe pulse. In this situation, however, one can obtain the trajectory deflection by resorting to a numerical simulation for solving Eq. (6). Figure 2(a) shows the result of the 3D motion trajectory of an LB as a function of $x / R_{0}$ , $y / R_{0}$ , and $z / (2 L_{diff})$ , in the presence of the gradient magnetic field with $(B_{1}, B_{2}) = (3.2,0) mG {cm}^{- 1}$ . The corresponding trajectory in the $x - z$ plane is illustrated in Fig. 2(b). We see clearly that the LB experiences a deflection due to the existence of the magnetic field. Shown in Fig. 2(c) is the result of the 3D motion trajectory of the LB for an increased magnetic-field gradient in the $x$ direction by taking $(B_{1}, B_{2}) = (6.4,0) mG {cm}^{- 1}$ . One sees that the trajectory of the LB is changed significantly due to the increase of the magnetic field.

Figure 2.Stern–Gerlach deflections of nonlocal LBs. (a) 3D motion trajectory of an LB as a function of x/R₀, y/R₀, and z/(2L_diff) in the presence of the gradient magnetic field (B₁,B₂) = (3.2, 0) mG cm⁻¹; (c) 3D motion trajectory of the LB for (B₁, B₂) = (6.4, 0) mG cm⁻¹. (b) and (d) are trajectories of the LB in the x–z plane, corresponding, respectively, to panels (a) and (c).

In addition, richer motion trajectories of the LB can be obtained by using different magnetic fields. To prove this, we consider a time-varying gradient magnetic field of the form $B (x, t) = \hat{z} B (x, t) = \hat{z} B_{0} \cos (ω_{0} t) x,$ (9)where $ω_{0}$ characterizes the motion period of the magnetic field in time. Figure 3(a) shows the motion trajectory of the LB under the action of such a magnetic field. We see that the trajectory of the LB follows the variation of the magnetic field. Illustrated in Fig. 3(b) is the corresponding sinusoidal trajectory of the LB in the $x - z$ plane. Obviously, one can use various magnetic fields to manipulate the motion of LBs; conversely, the trajectory deflections of the LBs may be exploited to measure external magnetic fields.

Figure 3.Motion trajectory of the LB in the presence of a time-varying gradient magnetic field. (a) Trajectory of the LB as a function of x/R₀, y/R₀, and z/(2L_diff) when the time-varying gradient magnetic field of Eq. (9) is present. (b) The corresponding sinusoidal trajectory of the LB in the x–z plane.

5. Conclusion

We have shown that nonlocal LBs created in a cold Rydberg atomic gas can be actively manipulated by using a weak gradient magnetic field. In particular, the LBs can experience significant Stern–Gerlach deflections when a weak external magnetic field is applied, and their motion paths may be controlled through the adjustment of the magnetic-field gradient. The results reported here are useful not only for understanding novel properties of the LBs in nonlocal optical media but also for finding new ways for precision measurements of magnetic fields.

References

[1] M. Saffman, T. G. Walker, K. Mϕlmer. Quantum information with Rydberg atoms. Rev. Mod. Phys., 82, 2313(2010).

[2] J. D. Pritchard, K. J. Weatherill, C. S. Adams. Nonlinear optics using cold Rydberg atoms. Annual Review of Cold Atoms and Molecules, 1, 301(2012).

[3] O. Firstenberg, C. S. Adams, S. Hofferberth. Nonlinear quantum optics mediated by Rydberg interactions. J. Phys. B, 49, 152003(2016).

[4] C. Murray, T. Pohl. Quantum and nonlinear optics in strongly interacting atomic ensembles. Advances in Atomic, Molecular, and Optical Physics, 65, 321(2016).

[5] J. D. Pritchard, D. Maxwell, A. Gauguet, K. J. Weatherill, M. P. A. Jones, C. S. Adams. Cooperative atom-light interaction in a blockaded Rydberg ensemble. Phys. Rev. Lett., 105, 193603(2010).

[6] S. Sevinçli, N. Henkel, C. Ates, T. Pohl. Nonlocal nonlinear optics in cold Rydberg gases. Phys. Rev. Lett., 107, 153001(2011).

[7] D. Yan, C. Cui, Y. Liu, L. Song, J. Wu. Normal and abnormal nonlinear electromagnetically induced transparency due to dipole blockade of Rydberg excitation. Phys. Rev. A, 87, 023827(2013).

[8] A. V. Gorshkov, J. Otterbach, M. Fleischhauer, T. Pohl, M. D. Lukin. Photon-photon interactions via Rydberg blockade. Phys. Rev. Lett., 107, 133602(2011).

[9] O. Firstenberg, T. Peyronel, Q. Liang, A. V. Gorshkov, M. D. Lukin, V. Vuletić. Attractive photons in a quantum nonlinear medium. Nature, 502, 71(2013).

[10] D. Maxwell, D. J. Szwer, D. Paredes-Barato, H. Busche, J. D. Pritchard, A. Gauguet, K. J. Weatherill, M. P. A. Jones, C. S. Adams. Storage and control of optical photons using Rydberg polaritons. Phys. Rev. Lett., 110, 103001(2013).

[11] L. Li, A. Kuzmich. Quantum memory with strong and controllable Rydberg-level interactions. Nat. Commun., 7, 13618(2016).

[12] D. Tiarks, S. Schmidt, G. Rempe, S. Dürr. Optical π phase shift created with a single-photon pulse. Sci. Adv., 2, e1600036(2016).

[13] H. Busche, P. Huillery, S. W. Ball, T. Ilieva, M. P. A. Jones, C. S. Adams. Contactless nonlinear optics mediated by long-range Rydberg interactions. Nat. Phys., 13, 655(2016).

[14] Z. Bai, G. Huang. Enhanced third-order and fifth-order Kerr nonlinearities in a cold atomic system via Rydberg-Rydberg interaction. Opt. Express, 24, 4442(2016).

[15] H. Wu, Y. Li, Z. Yang, S. Zheng. Quantum signature for laser-driven correlated excitation of Rydberg atoms. Phys. Rev. A, 95, 013842(2017).

[16] Q. Y. Liang, A. V. Venkatramani, S. H. Cantu, T. L. Nicholson, M. J. Gullans, A. V. Gorshkov, J. D. Thompson, C. Chin, M. D. Lukin, V. Vuletić. Observation of three-photon bound states in a quantum nonlinear medium. Science, 359, 783(2018).

[17] Z. Bai, W. Li, G. Huang. Stable single light bullets and vortices and their active control in cold Rydberg gases. Optica, 6, 309(2019).

[18] S. Bai, X. Han, J. Bai, Y. Jiao, J. Zhao, S. Jia, G. Raithel. Cesium nDJ+6S1/2 Rydberg molecules and their permanent electric dipole moments. Phys. Rev. Res, 2, 033525(2020).

[19] D. Ding, H. Busche, B. Shi, G. Guo, C. S. Adams. Phase diagram and self-organizing dynamics in a thermal ensemble of strongly interacting Rydberg atoms. Phys. Rev. X, 10, 021023(2020).

[20] C. Fan, D. Rossini, H. Zhang, J. Wu, M. Artoni, G. C. La Rocca. Discrete time crystal in a finite chain of Rydberg atoms without disorder. Phys. Rev. A, 101, 013417(2020).

[21] M. Zhou, J. Liu, P. Sun, Z. An, J. Li, X. Bao, J. Pan. Experimental creation of single Rydberg excitations via adiabatic passage. Phys. Rev. A, 102, 013706(2020).

[22] K. Liao, H. Tu, S. Yang, C. Chen, X. Liu, J. Liang, X. Zhang, H. Yan, S. Zhu. Microwave electrometry via electromagnetically induced absorption in cold Rydberg atoms. Phys. Rev. A, 101, 053432(2020).

[23] M. Jing, Y. Hu, J. Ma, H. Zhang, L. Zhang, L. Xiao, S. Jia. Quantum superhet based on microwave-dressed Rydberg atoms. Nat. Phys., 16, 911(2020).

[24] Y. Liu, Y. Sun, Z. Fu, P. Xu, X. Wang, X. He, J. Wang, M. Zhan. Infidelity induced by ground-Rydberg decoherence of the control qubit in a two-qubit Rydberg-blockade gate. Phys. Rev. Appl., 15, 054020(2021).

[25] J. Sinclair, D. Angulo, N. Lupu-Gladstein, K. Bonsma-Fisher, A. M. Steinberg. Observation of a large, resonant, cross-Kerr nonlinearity in a cold Rydberg gas. Phys. Rev. Res., 1, 033193(2019).

[26] C. Chen, F. Yang, X. Wu, C. Shen, M. K. Tey, L. You. Two-color optical nonlinearity in an ultracold Rydberg atom gas mixture. Phys. Rev. A, 103, 053303(2021).

[27] M. Moreno-Cardoner, D. Goncalves, D. E. Chang. Quantum nonlinear optics based on two-Rydberg atom arrays(2021).

[28] S. Baur, D. Tiarks, G. Rempe, S. Dürr. Single-photon switch based on Rydberg blockade. Phys. Rev. Lett., 112, 073901(2014).

[29] W. Li, I. Lesanovsky. Coherence in a cold-atom photon switch. Phys. Rev. A, 92, 043828(2015).

[30] H. Gorniaczyk, C. Tresp, P. Bienias, A. Paris-Mandoki, W. Li, I. Mirgorodskiy, H. P. Büchler, I. Lesanovsky, S. Hofferberth. Enhancement of Rydberg-mediated single-photon nonlinearities by electrically tuned Fórster resonances. Nat. Commun., 7, 12480(2016).

[31] C. R. Murray, T. Pohl. Coherent photon manipulation in interacting atomic ensembles. Phys. Rev. X, 7, 031007(2017).

[32] D. Tiarks, S. Baur, K. Schneider, S. Dürr, G. Rempe. Single-photon transistor using a Förster resonance. Phys. Rev. Lett., 113, 053602(2014).

[33] F. Ripka, H. Kübler, R. Löw, T. Pfau. A room-temperature single-photon source based on strongly interacting Rydberg atoms. Science, 362, 446(2018).

[34] Y. Silberberg. Collapse of optical pulses. Opt. Lett., 15, 1282(1990).

[35] Y. S. Kivshar, G. P. Agrawal. Optical Solitons: From Fibers to Photonic Crystals(2006).

[36] M. Belić, N. Petrović, W.-P. Zhong, R.-H. Xie, G. Chen. Analytical light bullet solutions to the generalized (3 + 1)-dimensional nonlinear Schrödinger equation. Phys. Rev. Lett., 101, 123904(2008).

[37] Y. V. Kartashov, B. A. Malomed, L. Torner. Solitons in nonlinear lattices. Rev. Mod. Phys., 83, 247(2011).

[38] Q. Zhang, Z. Bai, G. Huang. Fast-responding property of electromagnetically induced transparency in Rydberg atoms. Phys. Rev. A, 97, 043821(2018).

[39] S. Mauger, J. Millen, M. P. A. Jones. Spectroscopy of strontium Rydberg states using electromagnetically induced transparency. J. Phys. B, 40, F319(2007).

[40] The frequency and wave number of the probe field are given by ωp+ω and kp+K(ω), respectively. Thus, ω=0 corresponds to the center frequency of the probe field.

[41] Z. Chen, G. Huang. Trapping of weak signal pulses by soliton and trajectory control in a coherent atomic gas. Phys. Rev. A, 89, 033817(2014).

[42] In the related numerical simulation, a 10% random disturbance has been added into the initial condition.