Extreme single-excitation subradiance from two-band Bloch oscillations in atomic arrays

Luojia Wang; Da-Wei Wang; Luqi Yuan; Yaping Yang; Xianfeng Chen

doi:10.1364/PRJ.506450

Abstract

Atomic arrays provide an important quantum optical platform with photon-mediated dipole–dipole interactions that can be engineered to realize key applications in quantum information processing. A major obstacle for such applications is the fast decay of the excited states. By controlling two-band Bloch oscillations of single excitation in an atomic array under an external magnetic field, here we show that exotic subradiance can be realized and maintained with orders of magnitude longer than the spontaneous decay time in atomic arrays with the finite size. The key finding is to show a way for preventing the wavepacket of excited states scattering into the dissipative zone inside the free space light cone, which therefore leads to the excitation staying at a subradiant state for an extremely long decay time. We show that such operation can be achieved by introducing a spatially linear potential from the external magnetic field in the atomic arrays and then manipulating interconnected two-band Bloch oscillations along opposite directions. Our results also point out the possibility of controllable switching between superradiant and subradiant states, which leads to potential applications in quantum storage.

1. INTRODUCTION

Light–matter interaction in subwavelength scale is of broad interest in quantum information processing and quantum metrology [1,2]. Such interactions in subwavelength atomic arrays [3 –28] can induce long-range nonlinear dipole–dipole interactions by mediating photons in radiation modes [3,5,29,30] or via long-range van der Waals interactions between Rydberg states [31]. Further controls of atomic arrays by magnetic fields bring remarkable physical phenomena, including super- and subradiant states [32], photon storage and retrieval [33 –35], subradiance-protected quantum state transport [36], and many others [37 –45]. These quantum information processing methods utilize subradiant quantum states with inhibited spontaneous decay and exhibit useful applications in quantum storage [46]. Previous research shows that single-excitation subradiant states in one-dimensional (1D) arrays of $N$ atoms typically have decay rates scaling as $N^{- 3}$ with a lifetime of 7 orders of magnitude greater than that of a single atom for $N \sim 200$ [6,11,15]. Nevertheless, it generally desires a quantum state lasting long enough for the purpose of quantum storage [46]. A challenging problem is that subradiant states can diffuse into the superradiant subspace through the interaction or boundary effect, which limits their applications in quantum information processing.

Here we show that a judicious coherent control of finite-size atomic arrays can prevent single-excitation subradiant states from entering the superradiant subspace and thus realize subradiance with a lifetime of orders of magnitude longer than that of a single atom with the help of an external magnetic field. Specifically, we consider a well-established model with atomic arrays under the magnetic field, which lifts the degeneracy of excited states in each atom through the Zeeman shift (labeled as $| + ⟩$ and $| - ⟩$ ) [37,38,41]. Such systems with an infinite size can support subradiant collective excited states outside of the dissipative zone in the momentum space [6,11,15]. However, for the same system with a finite size, such excitation inevitably spreads to the boundaries and the wavepacket gets scattered into the dissipative zone. To overcome such obstacles, we instead apply a spatially linear magnetic field, which brings linear potentials with opposite slopes on $| + ⟩$ and $| - ⟩$ states [see Fig. 1(b)]. Similar 1D lattices under linear potentials support Bloch oscillations [47 –58]. Fundamentally different from the previous work, Bloch oscillations in atomic arrays under a linear magnetic field support two-band Bloch oscillations with inter-band interactions such that the wavepacket of excited states oscillates on two bands along opposite directions alternatively. This unique picture can force the wavepacket spatially localized within a few atoms while preventing the excitation from entering into the dissipative zone. Hence it supports subradient oscillations surviving extremely long time with finite atoms. We also show the way of switching between superradiant and subradiant oscillations by tuning the external magnetic field, which provides the opportunity for reading out information from subradiant quantum states [59 –62].

Figure 1.(a) Schematic of a 1D V-type atomic array under a magnetic field $B (y)$ . (b) Schematic of the corresponding lattice where the $n$ th atom has non-degenerate excited states $| \pm_{n} ⟩$ with frequency shifted by $\pm μ B_{n} / ℏ$ . Blue lines indicate a linear trend of $B_{n}$ as in Eq. (2). (c) Band structures with different constant magnetic fields $B_{c}$ with values $B_{c} = 0, 4 ℏ γ_{0} / μ, 8 ℏ γ_{0} / μ, 12 ℏ γ_{0} / μ$ , respectively. Black arrows indicate values of $B_{c}$ for corresponding bands. Decay rates of modes are color coded. Probabilities of eigenstates on (d) band II and (e) band I projected on $| + ⟩$ states for different $B_{c}$ . Here, $a = 0.1 λ$ .

2. MODEL

The developments of state-of-the-art technologies in various experimental platforms such as trapped neutral atoms at the subwavelength scale [19 –28,63 –68] pave the way for exploring exotic phenomena in atomic arrays. We study such arrays of atoms that are equally spaced at a distance $a$ along the $y$ -axis [ $y_{n} = n a$ for the $n$ th atom in Fig. 1(a)]. Each atom has a ground state $| g ⟩$ and degenerate excited states. We apply an off-plane magnetic field $B (y)$ in the $z$ -direction to lift the degeneracy of two excited states $| \pm ⟩$ . Under the single-excitation limit, the dynamics of the wave packet of the excited states is described by the non-Hermitian Hamiltonian after integrating out photonic modes [37,38,41], $H_{eff} = ℏ \sum_{n = 1}^{N} \sum_{α = \pm} (ω_{A} - i \frac{γ_{0}}{2}) | α_{n} ⟩ ⟨ α_{n} | + \sum_{n = 1}^{N} μ B_{n} (| +_{n} ⟩ ⟨ +_{n} | - | -_{n} ⟩ ⟨ -_{n} |) + \frac{3 π ℏ γ_{0} c}{ω_{A}} \sum_{n \neq m} \sum_{α, β = \pm} G_{α β} (y_{n} - y_{m}) | α_{n} ⟩ ⟨ β_{m} |,$ (1)where $ω_{A}$ is the atomic transition frequency, $γ_{0}$ is the atomic decay rate in the free space, $μ B_{n} / ℏ$ is the Zeeman shift for the $n$ th atom with the magnetic moment $μ$ , and $G_{α β} (y_{n} - y_{m})$ is the free-space dyadic Green’s function describing the electric field at $y_{n}$ emitted by the atom located at $y_{m}$ . With the relation $| \pm ⟩ = \mp (| x ⟩ \pm i | y ⟩) / \sqrt{2}$ , one obtains $G_{+ +} = G_{- -} = - e^{i k_{0} r} \cdot (k_{0}^{2} r^{2} - i k_{0} r + 1) / (8 π k_{0}^{2} r^{3})$ and $G_{+ -} = G_{- +} = e^{i k_{0} r} (k_{0}^{2} r^{2} + 3 i k_{0} r - 3) / (8 π k_{0}^{2} r^{3})$ , where $r = | y_{n} - y_{m} |$ and $k_{0} = ω_{A} / c = 2 π / λ$ with $λ$ being the wavelength and $c$ being the vacuum speed of light [38,41]. Equation (1) shows a tight-binding lattice model where excitations on sites at two arms ( $| \pm ⟩$ ) have on-site potentials $V_{\pm, n} \equiv \pm μ B_{n}$ and are connected by complex and long-range photon-mediated hoppings. In particular, we consider the magnetic field $B_{n} = n B_{0},$ (2)where $B_{0}$ is a constant. Equation (2) gives $V_{\pm, n} = \pm n μ B_{0}$ , which leads to the effective constant force $F_{\pm} = \mp μ B_{0} / a$ in opposite directions on the two excited states. Two effective electric fields in opposite directions on two separated arms in the lattice bring Bloch oscillations on two arms exhibiting symmetric patterns. However, once the two arms are connected [by terms including $G_{\pm \mp}$ in Eq. (1)], the dynamics of the Bloch-oscillation wavepackets on the two arms is influenced by each other, which can be analyzed in the field-free Bloch band picture.

Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you！Sign up now

We first study the band structure of the infinite atomic arrays with constant magnetic field $B_{c}$ [37,38]. The single-excitation Bloch modes are given by $| Ψ_{k_{y}} ⟩ = \sum_{n} (C_{+, k_{y}} | +_{n} ⟩ + C_{-, k_{y}} | -_{n} ⟩) e^{i k_{y} n a} = {(\begin{array}{l} C_{+, k_{y}} & C_{-, k_{y}} \end{array})}^{T},$ (3)where $k_{y}$ is the Bloch wave vector. The Hamiltonian in the momentum space on two arms has the form $\frac{H_{k}}{ℏ} = (\begin{matrix} ω_{A} - i \frac{γ_{0}}{2} + χ_{k_{y}}^{+ +} & χ_{k_{y}}^{+ -} \\ χ_{k_{y}}^{- +} & ω_{A} - i \frac{γ_{0}}{2} + χ_{k_{y}}^{- -} \end{matrix}) + (\begin{matrix} i \frac{μ B_{0}}{ℏ a} \frac{d}{d k_{y}} & 0 \\ 0 & - i \frac{μ B_{0}}{ℏ a} \frac{d}{d k_{y}} \end{matrix}),$ (4)where $χ_{k_{y}}^{α β} = \frac{3 π γ_{0} c}{ω_{A}} \sum_{n \neq 0} e^{i k_{y n a}} G_{α β} (n a)$ and diagonalization of the first term plus the shifts caused by the constant magnetic field gives the two-band structure. Our results show that the oscillatory mode for an extended Gaussian distribution of the single-excitation wavepacket approximately evolves in a single band. The second term in Eq. (4) would transform the distribution on one band into two bands and gives the acceleration force in one band and the tunneling strength to the other band, which surpasses the former one if the two arms have comparable probabilities. We thus bring in a phenomenological force averaging on two arms, $⟨ F ⟩ \equiv \frac{μ B_{0}}{a} \sum_{n} ⟨ Ψ | (| +_{n} ⟩ ⟨ +_{n} | - | -_{n} ⟩ ⟨ -_{n} |) | Ψ ⟩ = \frac{μ B_{0}}{a} \sum_{n} (P_{+} - P_{-})$ , to estimate the acceleration force.

We choose $a = 0.1 λ$ , which brings the first Brillouin zone $k_{y} \in [- 5 k_{0},5 k_{0}]$ . We plot band structures in $k_{y}$ -space for different $B_{c}$ in Fig. 1(c). For each $B_{c}$ , there are two bands, and it exhibits large collective decay, i.e., the dissipative zone in the free-space light cone ( $| k_{y} | < k_{0}$ ). Outside of the dissipative zone, subradiant states are supported. For $B_{c} = 0$ , there are two degenerate points near $k_{y} = \pm 2.3 k_{0}$ , while two bands are separated in the entire $k_{y}$ -space for other $B_{c}$ . We denote the lower (upper) band as band I (II) and define $P_{\pm}^{I (II)}$ as the probability intensity of an eigenstate on band I (II) projected onto the arm of $| \pm ⟩$ states. $P_{+}^{II}$ and $P_{+}^{I}$ are shown in Figs. 1(d) and 1(e), respectively. One notes that $P_{+}^{I (II)} + P_{-}^{I (II)} = 1$ . When there is no magnetic field, $P_{\pm}^{I} = P_{\pm}^{II} = 0.5$ in the entire $k_{y}$ -space, indicating that probabilities for exciting $| \pm ⟩$ states are equivalent due to the degeneracy of $| \pm ⟩$ . However, for $m \neq 0$ , it shows asymmetric distributions of $P_{\pm}^{I (II)}$ . In particular, for $B_{c} > 0$ on band I (II), we can see that $P_{+}^{I} < 0.5$ ( $P_{+}^{II} > 0.5$ ), meaning that less (more) density of the state is located on the arm of $| + ⟩$ states. Distributions of $P_{\pm}^{I (II)}$ are highly symmetric because the influence of $| \pm ⟩$ is symmetric under the positive/negative magnetic field. Moreover, we find that, for $B_{c} > 0$ at $k_{y} = \pm 2.3 k_{0}$ , it shows $P_{+}^{I} = 0$ ( $P_{+}^{II} = 1$ ) on band I (II). We note that band structures in Fig. 1 are calculated by simply summing the Green’s function in real space over all lattice sites, which artificially induces the inaccurate eigenvalues at $k_{y}$ near $\pm k_{0}$ in the free-space light cone due to slow convergence. However, this inaccuracy in band structures does not affect the later analysis throughout this paper. Also, if desired, it can be overcome by summations of Green’s function in momentum space with an appropriate regularizing method [37,38].

We next introduce the simulation method for studying the dynamics of the single excitation wavepacket in the atomic arrays. We use the Schrödinger equation $d | Ψ (t) ⟩ / d t = - i H_{eff} | Ψ (t) ⟩ / ℏ$ and the excitation wavepacket of the atomic array $| Ψ (t) ⟩ = \sum_{n} [C_{+, n} (t) | +_{n} ⟩ + C_{-, n} (t) | -_{n} ⟩] e^{- i ω_{A} t}$ , where ${| C_{\pm, n} |}^{2}$ gives the excitation probability of the $| \pm ⟩$ state in the $n$ th atom. Two hundred one atoms ( $n = - 100, \dots,0, \dots,100$ ) are considered in the arrays. $B_{0} = 0.2 ℏ γ_{0} / μ$ is used in simulations. We assume that the wavepacket of the system is initially prepared at a superposition state on two arms following $C_{\pm, n} (0) = ψ_{\pm} \exp [i k_{c} a n - {(n - n_{c})}^{2} / 200]$ , where $n_{c}$ is the spatial center of the initial excitation, $k_{c}$ is the initial momentum, and $ψ_{\pm}$ gives the initial ratio of excitation amplitudes on two arms of the $| \pm ⟩$ states. Such an initial state in the weak excitation limit can be pumped with phase-controlled schemes [32,33,59,69] or by applying a spatial modulation of the atomic detuning [18].

3. RESULTS

For finite atomic arrays, the key ingredient for maintaining a long-standing subradiant state is to prevent the wavepacket from entering into the dissipative zone either by oscillations or scattering due to the boundary effects. One might notice that the well-established Bloch oscillation may lead the wavepacket oscillation within a finite spatial region and propagating unidirectionally in the $k_{y}$ -space [50,53,70], so by frequently altering the direction of the constant force, it is possible to accomplish the task for avoiding entering the dissipative zone. Nevertheless, here we show that one can realize such a task in atomic arrays under a spatially linear but temporally constant magnetic field [Eq. (2)], which leads to effective constant forces in opposite directions on the two excited states. Once the two arms are connected via photon-mediated dipole–dipole interactions, the effective electric force that acts on the single-excitation wavepacket depends on the wavepacket distributions between the two arms. Our results show that an extended Gaussian distribution of the single-excitation wavepacket approximately evolves in a single band driven by the phenomenological force $⟨ F ⟩$ , which can be estimated by the eigenstate distributions of the moving wavepacket in a band structure under a constant magnetic field felt by the wavepacket center. The ability of altering the effective force direction during Bloch oscillations provides the mechanism to build subradiant states that avoid entering the dissipative zone inside the free-space light cone and are spatially localized far away from the boundary of the finite array.

To see the long-standing subradiant state, we excite the wavepacket centered at $n = 0$ , i.e., $n_{c} = 0$ , with $ψ_{+} = ψ_{-} = 0.168$ and $k_{c} = 1.5 k_{0}$ , which satisfies $\sum_{n} [{| C_{+, n} (0) |}^{2} + {| C_{-, n} (0) |}^{2}] = 1$ . The evolutions in the simulation are plotted for the time up to $10 T_{B}$ , where $T_{B} = 2 π ℏ / μ B_{0} = 10 π / γ_{0}$ . Figures 2(a1) and 2(b1) depict the dynamics of ${| C_{+, n} (t) |}^{2}$ and ${| C_{-, n} (t) |}^{2}$ , respectively, and one sees a sustained Bloch oscillation pattern with the shape of the wavepacket deforming gradually. The excitation of the $| + ⟩$ ( $| - ⟩$ ) states in the atomic arrays is mainly located in the finite region with $y < 0$ ( $y > 0$ ) so there is no boundary effects due to finite atoms.

Figure 2.Bloch oscillations for Gaussian excitations initially centered at $n_{c} = 0$ with $k_{c} = 1.5 k_{0}$ on band I (left) and $k_{c} = 4 k_{0}$ on band II (right), shown by temporal evolution for excitation probabilities of (a1), (a2) ${| C_{+, n} |}^{2}$ ; (b1), (b2) ${| C_{-, n} |}^{2}$ ; (c1), (c2) $P^{I} (k_{y})$ ; (d1), (d2) $P^{II} (k_{y})$ ; (e1), (e2) $P_{+}$ , $P_{-}$ , and $P_{t}$ . Here, $a = 0.1 λ$ and $μ B_{0} / ℏ = 0.2 γ_{0}$ .

We then perform fast Fourier transform (FFT) on the simulation results $C_{\pm, n} (t)$ and calculate the projection of the excited wavepacket in $k_{y}$ -space onto the two bands [excitation probabilities $P^{I} (k_{y}, t)$ and $P^{II} (k_{y}, t)$ , respectively]. Note that the Bloch band picture is valid only with constant magnetic field. Here, the linear magnetic field and resulting constant force produce perturbation, so we can calculate the projection by taking into account the uniform potentials at the value of $V_{\pm, n} = \pm \bar{n} μ B_{0}$ , where $\bar{n} \equiv ⟨ Ψ | \sum_{n} (| +_{n} ⟩ n ⟨ +_{n} | + | -_{n} ⟩ n ⟨ -_{n} |) | Ψ ⟩$ is the mean position of the wavepacket. This analysis in the Bloch band picture with uniform on-site potentials reveals dynamical features different from a conventional Bloch oscillation. Figures 2(c1) and 2(d1) show the evolution of excitation probabilities of band I and band II and $P^{I}$ and $P^{II}$ at each $k_{y}$ . The choice of $ψ_{\pm}$ makes the initial excitation on band I, and throughout the evolution, $P^{II}$ is about 1 order of magnitude smaller than $P^{I}$ . Remarkably, one sees that $P^{I}$ does not evolve unidirectionally on $k_{y}$ with time, which gives the fundamental difference from the conventional Bloch oscillations. The evolution direction of $P^{I}$ reverses every time (referred to as the reverse time) before it enters the free-space light cone in the $k_{y}$ -axis [dash lines in Figs. 2(c1) and 2(d1)], which leads to the important consequence of avoiding a large collective decay of the excited wavepacket. We further show the total excitation probabilities of the $| + ⟩$ state and the $| - ⟩$ state ( $P_{+}$ , $P_{-}$ ) in the real space and the total excitation probability $P_{t} = P_{+} + P_{-}$ in Fig. 2(e1). We find that the times when $P_{+} = P_{-}$ are exactly the same as the reverse times because the phenomenological force $⟨ F ⟩$ becomes zero and subsequently changes its direction. In other words, when the center of the wavepacket arrives at $n = 0$ , $P_{+} \approx P_{-}$ and hence $⟨ F ⟩$ changes its direction, resulting in the reverse of $P^{I} (k_{y})$ in the momentum space. Moreover, $P_{t}$ exhibits no decay in this dissipative atomic system, which is actually lasting $\sim 10^{8} γ_{0}^{- 1}$ given by the simulation with a much longer evolution time. Therefore, as long as the excited wavepacket in the atomic arrays does not enter the dissipative zone inside the free-space light cone, the life time of the excitation gives the subradiant feature, which is much longer than the free-space decay time. Such extreme subradiant oscillations can be optimized by choosing initial excitations. In particular, we find that $P_{t}$ can last more than $\sim 10^{13} γ_{0}^{- 1}$ for an initial state centered at $k_{c} = 3 k_{0}$ on band I. Based on the eigenvalue analysis of the effective Hamiltonian for an array under a linear magnetic field, one can determine the decay rate of the most subradiant eigenstate, which is consistent with the property of the optimized subradiant Bloch oscillations (see Appendix A for details).

As one may notice, the initial choice of the excited wavepacket determines how the wavepacket evolves in $k_{y}$ -space and hence how the wavepacket decays. For example, if we choose $ψ_{+} = 0.168$ , $ψ_{-} = - 0.168$ , and $k_{c} = 4 k_{0}$ to excite band I or $ψ_{+} = 0.168$ , $ψ_{-} = - 0.168$ , and $k_{c} = 2 k_{0}$ to excite band II, we find the total probability of the wavepacket undergoes no decay even though the evolution patterns are different (see Fig. 5 in Appendix A). However, in the case of $ψ_{+} = ψ_{-} = 0.168$ and $k_{c} = 4 k_{0}$ , with band II being excited, the dramatically different dynamics is as shown in Figs. 2(a2)–2(e2). One notes that the excitation undergoes a rapid decay after the partial oscillation during the initial time duration $\sim 11 γ_{0}^{- 1}$ . The fundamental difference here is that the excitation probability $P^{II}$ moves into the free-space light cone along the $k_{y}$ -axis before $⟨ F ⟩$ changes its direction. This result leads to a very fast decay within $\sim 8 γ_{0}^{- 1}$ for the total probability $P_{t}$ dropping to $\sim 0.02$ . In real space, the excitation on the $| - ⟩$ arm decreases fast, while the excitation on the $| + ⟩$ arm increases initially with an oscillation towards the positive $y$ -axis but decays rapidly afterwards.

The two-band Bloch oscillation in the atomic arrays under the linear magnetic field brings unique opportunity for manipulating the decay of the excited wavepacket. For instance, we can control the collective decay exhibiting quantization-like decay with stable plateaus, as shown in Fig. 3(c1). To achieve it, we choose $ψ_{+} = ψ_{-} = 0.168$ and $k_{c} = k_{0}$ , so the wavepacket on band I centered at $k_{0}$ is initially excited. The simulation results are summarized in Figs. 3(a1)–3(c1). The wavepacket experiences decay at $t \sim 0$ because the initial wavepacket is partially inside the free-space light cone. However, the wavepacket initially moves toward the positive direction in the $k_{y}$ -axis [see $P^{I} (k_{y}) + P^{II} (k_{y})$ in Fig. 3(b1)], so $P_{t}$ ends in a stable plateau at $t \sim 2 γ_{0}^{- 1}$ with a value $\sim 0.67$ . At $t \sim 27 γ_{0}^{- 1}$ , the excited wavepacket is reached at $k_{y} = - k_{0}$ , and a reversal happens due to the change of $⟨ F ⟩$ , which results in a rapid drop of $P_{t}$ to $\sim 0.48$ . Periodic oscillations are exhibited, together with the oscillations of $P_{n} = {| C_{+, n} |}^{2} + {| C_{-, n} |}^{2}$ , in the real space. The total excitation probability $P_{t}$ therefore shows the overall decay tendency with the periodic plateaus [see Fig. 3(c1)]. Initial excitations inside the free-space light cone can exhibit similar oscillations but quicker decay. Particularly, for an initial wavepacket at $k_{c} = - 0.5 k_{0}$ on band II, $P_{t}$ decays to $\sim 0.02$ within $t \sim 0.85 γ_{0}^{- 1}$ , or a collective superradiant decay rate of $\sim 4.7 γ_{0}$ . Moreover, excitations with $n_{c} \neq 0$ can bring different Block oscillation phenomena because the central momenta at reverse times are changed during the evolution of the wavepacket (see Appendix A for details).

Figure 3.Bloch oscillations for Gaussian excitations initially centered at $n_{c} = 0$ and $k_{c} = k_{0}$ on band I with a static magnetic field (left) and a controllable magnetic field (right), whose zero point is shifted to $y = 30 a$ over a time period from 0 to $20 γ_{0}^{- 1}$ and back to $y = 0$ over a time period from $180 γ_{0}^{- 1}$ to $200 γ_{0}^{- 1}$ , shown by temporal evolution of (a1), (a2) $P_{n}$ ; (b1), (b2) $P^{I} + P^{II}$ ; (c1), (c2) $P_{+}$ , $P_{-}$ , and $P_{t}$ . Other parameters are the same as those in Fig. 2.

Taking advantage of the relation between reverse times and the atomic position in the linear magnetic field, one can control the momentum of the wavepacket at the reverse time. The magnetic field acts as an effective field on the single particle and breaks the spatial translation symmetry, which leads to the particle movement in the momentum space. Change of the magnetic field would alter the effective force and the momentum space evolution of the wavepacket. We give an example in Figs. 3(a2)–3(c2), where the spatial distribution of the linear magnetic field relative to the atomic position is tuned to control the evolution of the wavepacket that is initially excited at $k_{0}$ on band I [the same as that in Figs. 3(a1)–3(c1)]. In the numerical simulation, we tune the zero point of the magnetic field from $y = 0$ to $y = 30 a$ over a time period from 0 to $20 γ_{0}^{- 1}$ . As a result, after $P_{t}$ drops to the first plateau, the wavepacket oscillates around the new center $y = 30 a$ [Fig. 3(a2)] and reverses at $\sim \pm 1.4 k_{0}$ in the momentum space to maintain subradiance [Fig. 3(b2)]. The excitation can also be driven to the dissipative zone by moving the zero point of the magnetic field back to $y = 0$ over a time period from $180 γ_{0}^{- 1}$ to $200 γ_{0}^{- 1}$ and decay radiatively, shown by descending plateaus in the total excitation probability.

4. DISCUSSION AND CONCLUSION

The proposed atomic arrays can be realized in various experimental systems. Dependent on different experimental platforms, atomic arrays can be composed with atom ensembles in the gas phase, such as neutral particles, including Rydberg atoms trapped in optical lattices [19,21 –23,63,64] and optical tweezers [20,24 –28,65 –68], or atom-like particles in the solid phase, such as color defects in diamond nanophotonic devices [71,72], and excitons in atomically thin semiconductors [73]. Considering the possible imperfection of atomic positions and the inhomogeneous broadening in the atomic array, numerical results show that positional disorder suppresses the positive influence of the magnetic field on subradiance (see Appendix B for details). This limitation is expected to be overcome by the development of experimental control techniques. On the other hand, the effect of decay suppression can be found in arrays with different atomic numbers and interatomic distances to possibly decrease experimental challenges. The performance of subradiance relies on the magnitude of the linear magnetic field. The parameters used above present better initial state selections and characteristic performance of Bloch oscillations, which may further be optimized for extreme subradiance. Moreover, it has been shown recently that dispersionless or trapped subradiant photon states can be created in one-dimensional emitter chains [74]. Different from Ref. [74], our study shows that, in the model with the external magnetic field it is possible to decrease the subradiant decay rate to smaller than that in the atomic arrays without the magnetic field and also points out the possibility of quantum-state switching from the external-field control.

In summary, we have explored single-excitation dynamics in 1D atomic arrays under a linear magnetic field. Resulting two-band Bloch oscillations are studied. We found that the evolution of the wavepacket in the momentum space depends on the phenomenological force $⟨ F ⟩$ , which is determined by local properties of the band structure and thus influenced by the position of the wavepacket in the real space. Compared to the conventional Bloch oscillations, the evolution direction of the momentum reverses as long as the wavepacket passes the position where $⟨ F ⟩$ changes its sign. Taking advantage of this unique property, we show the capability of generating exotic subradiant oscillations that lasts orders of magnitude greater than the atomic spontaneous decay time. Our study takes advantage of the external control from the magnetic field and therefore points toward fundamental opportunities in switching quantum states between superradiant and subradiant states and realizing an extreme subradiant quantum state that could be useful in important applications for quantum storage.

Acknowledgment

Acknowledgment. L.Y. thanks the sponsorship from the Yangyang Development Fund.

APPENDIX A: BLOCH OSCILLATIONS FOR DIFFERENT INITIAL STATES

a = 0.1 λ

Figure 4.Bloch oscillations for a Gaussian excitation initially centered at $n_{c} = 0$ with $k_{c} = 1.5 k_{0}$ (left) and $k_{c} = 3 k_{0}$ (right) on band I. Temporal evolution of (a1), (a2) $P_{n} = {| C_{+, n} |}^{2} + {| C_{-, n} |}^{2}$ and (b1), (b2) $P^{I} (k_{y}) + P^{II} (k_{y})$ in an array with atomic number $N = 201$ within a time period of 300 $γ_{0}^{- 1}$ after $t = 10^{8} γ_{0}^{- 1}$ (left) and after $t = 10^{13} γ_{0}^{- 1}$ (right). (c1), (c2) The total probability $P_{t}$ versus order of magnitude of time for $N = 51$ (blue solid line), $N = 101$ (green dash-dotted line), $N = 151$ (red dash line), and $N = 201$ (black dotted line). $a = 0.1 λ$ and $μ B_{0} / ℏ = 0.2 γ_{0}$ .

Figure 5.Bloch oscillations started from $n_{c} = 0$ for Gaussian excitations initially centered at $k_{c} = 4 k_{0}$ on band I (left), $k_{c} = 2 k_{0}$ on band II (middle), and $k_{c} = - 0.5 k_{0}$ on band II (right), shown by temporal evolutions of (a1)–(a3) $P_{n}$ , (b1)–(b3) $P^{I} (k_{y}) + P^{II} (k_{y})$ , (c1), (c2) total excitation probabilities $P_{+}$ (red solid line), $P_{-}$ (black dash line), and $P_{t}$ (green dash-dotted line). An exponential function is used to fit the decay of the total excitation probability in (c3) (blue dotted line). Here, $N = 201$ , $a = 0.1 λ$ , and $μ B_{0} / ℏ = 0.2 γ_{0}$ .

Figure 6.Bloch oscillations started from $n_{c} = 50$ for a Gaussian excitation initially centered at $k_{c} = 1.5 k_{0}$ on band I, shown by the temporal evolution of (a) $P_{n}$ , (b) $P^{I} + P^{II}$ , (c) $P_{+}$ (red solid line), $P_{-}$ (black dash line), and $P_{t}$ (green dash-dotted line). Here $N = 201$ , $a = 0.1 λ$ , and $μ B_{0} / ℏ = 0.2 γ_{0}$ .

APPENDIX B: DECAY RATES OF SINGLE-EXCITATION EIGENSTATES

B_{0}

Figure 7.Numerically calculated decay rates in descending order of single-excitation eigenstates in an atomic array under a linear magnetic field with varying $μ B_{0} / ℏ γ_{0}$ indicated by different colors for (a) $a = 0.1 λ$ , $N = 201$ ; (b) $a = 0.1 λ$ , $N = 101$ ; and (c) $a = 0.2 λ$ , $N = 201$ .

Figure 8.Decay rates of the most subradiant single-excitation modes with (a) increasing spatial disorder characterized by the deviation width $δ a$ and (b) increasing Doppler broadening $Δ_{D}$ for $B_{0} = 0$ (blue dash line) and $B_{0} = 0.2 ℏ γ_{0} / μ$ (red solid line). $a = 0.1 λ$ and $N = 201$ .