Dynamical creation and detection of entangled many-body states in a chiral atom chain

TL;DR

This paper presents a method to detect and analyze entangled many-body states in chiral atom chains by monitoring emitted photons, enabling insights into complex quantum phases and dynamical behaviors relevant for current experimental setups.

Contribution

It introduces a photon monitoring technique to identify dark entangled states and probe dynamical phenomena in chiral light-matter systems, advancing experimental characterization methods.

Findings

Detection of dark entangled states via photon emission records

Observation of coexistence of dark entangled and mixed phases

Framework applicable to current experimental platforms

Abstract

Open quantum systems with chiral interactions can be realized by coupling atoms to guided radiation modes in waveguides or optical fibres. In their steady state these systems can feature intricate many-body phases such as entangled dark states, but their detection and characterization remains a challenge. Here we show how such collective phenomena can be uncovered through monitoring the record of photons emitted into the guided modes. This permits the identification of dark entangled states but furthermore offers novel capabilities for probing complex dynamical behavior, such as the coexistence of a dark entangled and a mixed phase. Our results are of direct relevance for current experiments, as they provide a framework for probing, characterizing and classifying dynamical features of chiral light-matter systems.

Equations

$$\begin{array}{rl}\dot{\rho }=& \mathcal{W}(\rho )=-i[H_{\text{laser}}+H_L+H_R,\rho ]\\ & +{\gamma }_L\mathcal{D}(J_L)\rho +{\gamma }_R\mathcal{D}(J_R)\rho +\Gamma \sum _{j=1}^N\mathcal{D}({\sigma }_j^{-})\rho .\end{array}$$

$$\begin{array}{rl}& H_L=-i\frac{{\gamma }_L}{2}\sum _{j<l}(e^{-ik|x_j-x_l|}{\sigma }_j^{+}{\sigma }_l^{-}+\text{h.c})\\ & H_R=-i\frac{{\gamma }_R}{2}\sum _{j>l}(e^{-ik|x_j-x_l|}{\sigma }_j^{+}{\sigma }_l^{-}+\text{h.c}),\end{array}$$

$$\begin{array}{r}\dot{\rho }=-i\left[H,\rho \right]+\gamma \mathcal{D}(J)\rho +\Gamma \sum _{j=1}^N\mathcal{D}({\sigma }_j^{-})\rho ,\end{array}$$

$$H=\Omega \left(J+J^{†}\right)-\sum _j{\delta }_j{\sigma }_j^{+}{\sigma }_j^{-}-\frac{i}{2}\Delta \gamma \sum _{j>l}({\sigma }_j^{+}{\sigma }_l^{-}-{\sigma }_l^{+}{\sigma }_j^{-}),$$

$${\mathcal{W}}_{s_L,s_R}(\bullet )=\mathcal{W}(\bullet )-J_L\bullet J_L^{†}(1-e^{-s_L})-J_R\bullet J_R^{†}(1-e^{-s_R}).$$

$$Z_t(s)=\sum _K\text{Tr}\left[{\rho }^{(K)}\right]e^{-sK}=\text{Tr}\left[\sum _K{\rho }^{(K)}{e}^{-sK}\right]=\text{Tr}\left[{\rho }_s\right].$$

Full text

††preprint: APS/123-QED

Introduction. The ability to interface quantum emitters with optical systems opens novel routes for investigating non-equilibrium phenomena in open condensed matter physics Lee and provides, potentially, a platform to perform quantum information processing Weim ; Stann ; Lvov ; Asen . In recent years, the open quantum dynamics of chiral systems, where the emission of photons into a waveguide presents a broken left-right symmetry, has been the object of intense investigation Dutta ; Scheu ; Lod ; Verm ; Arno2 . This propagation-direction-dependent light-matter interaction has been observed in a variety of systems, for instance atoms coupled to the evanescent field of a waveguide Say , and quantum dots in photonic nano-structures Coles .

From the quantum many-body physics perspective, it has been demonstrated that chiral coupling induces steady state entanglement between emitters. In particular, when a chain of spins is coupled to a waveguide, under precise conditions, peculiar phases of matter can form, among them pure entangled many-body dark states that are protected from decoherence Pich . However, it still remains an experimental challenge to realize the specific conditions required for the creation of these entangled dark states. Furthermore, it is generally a difficult task to detect and characterize such states once created. This poses obstacles for exploring and using of the rich variety of quantum states and many-body phases for practical purposes, e.g. in optical quantum computing protocols Met .

In this work, we show that the counting statistics of the photons emitted into the waveguide provides a means to detect and characterize the many-body state of a chiral atom chain. Our approach exploits that the guided photons do not only induce interactions between the emitters, but also carry information about their quantum state. In an experiment, measuring the output field of the waveguide thus allows to infer the occurrence of dark steady states. More generally, this perspective allows to probe intricate dynamical phenomena such as transitions and the coexistence between dynamical phases in the steady state. In the latter case, entangled states may occur as fluctuations in an intermittent dynamics and are heralded through characteristic features of the time-resolved photon count signal. To establish this correspondence on a formal level we utilize the scaled cumulant generating function formalism together with the large deviation principle touchette2009 ; Garr . This does not only allow to develop an understanding of the dynamical non-equilibrium behavior of a chirally coupled atom chain, but also permits to systematically assess the effect of inevitable imperfections, such as the emission of photons into unguided modes.

Model. The general setup that we consider is depicted in Fig. 1. Here, $N$ spin-$1/2$ particles form a chain parallel to a one-dimensional waveguide, placed at positions $x_{j}$, with $j=1,\dots,N$. Note that in the following we explicitly refer to two-level atoms with internal states $\ket{g}$ and $\ket{e}$, but the same description applies, e.g. to quantum dots in photonic crystals. The atoms are trapped at a distance $d$ from the center of the waveguide and separated by a distance $a$ from each other. Each atom is driven by a coherent driving field with frequency $\omega_{p}$, Rabi frequency $\Omega$, and a (possibly spatially dependent) detuning $\delta_{j}=\omega_{p}-\omega_{j}$, with $\omega_{j}$ being the transition frequency between the ground and the excited state. The Hamiltonian that governs the dynamics of the laser excitation (under the rotating wave approximation) is $H_{\textrm{laser}}=\sum_{j=1}^{N}\left[-\delta_{j}\sigma^{+}_{j}\sigma^{-}_{j}+\Omega(\sigma^{+}_{j}+\sigma^{-}_{j})\right]$, where $\sigma^{-}_{j}=\ket{g}_{j}\!\bra{e}$.

The atoms are coupled both to the left and right propagating guided modes of the waveguide as well as to the unguided modes of the electromagnetic field outside the waveguide. Two photon detectors are placed at the extremes of the waveguide which allows in principle to resolve the photons emitted into both directions. Under the Born-Markov approximation, the open quantum dynamics of the atoms can be described by a Lindblad master equation:

[TABLE]

Here, $H_{L}$ and $H_{R}$ describe induced coherent dipole-dipole interactions due to the coupling between the atom chain and the left and right propagating modes of the waveguide, respectively. They read

[TABLE]

with $k=\omega_{p}/c$. The incoherent collective spin decay into the guided modes is described by the dissipators $\gamma_{L,R}\mathcal{D}(J_{L,R})$, where $J_{L,R}=\sum_{j}e^{\pm ikx_{j}}\sigma^{-}_{j}$ are jump operators describing the emission of a photon into the guided modes, and $\mathcal{D}(\chi)\rho\equiv\chi\rho\chi^{\dagger}-\frac{1}{2}\left\{\chi^{\dagger}\chi,\rho\right\}$. The last term of Eq. (1) represents the coupling of the atoms to the modes of the free (unguided) electromagnetic field with rate $\Gamma$ 111Note, that this expression is only valid when the separation between the atoms is sufficiently large, $ka\gg 1$. When $ka<1$, the emission into the unguided modes becomes collective instead of local, and the exchange of virtual photons gives rise to coherent interactions between the atoms Lehm ..

In the remainder of the paper we assume, for concreteness, that the spacing between the atoms, $a$, is commensurate with the wavelength of the laser field: $\lambda=\frac{2\pi}{k}$, such that $kx_{j}=2\pi n$ with $n=1,2,\dots$ for all $j=1,\dots,N$. Under these conditions, the master equation (1) simplifies considerably and takes the form

[TABLE]

with $\gamma=\gamma_{R}+\gamma_{L}$ and

[TABLE]

where $J=\sum_{j}\sigma_{j}^{-}$ and $\Delta\gamma=\gamma_{R}-\gamma_{L}$. For now (we will return to this at a later stage) we also neglect decay into the unguided modes, i.e., we consider $\Gamma=0$.

We are here particularly interested in the case of chiral interactions, i.e. when the couplings with the left and the right propagating modes are not symmetric and hence $\Delta\gamma\neq 0$. This is realized by an appropriate choice of the laser polarization and transition dipole moment of each atom: as shown in Ref. LeKi , the dipole moments are required to have a real and imaginary part (elliptically polarized light) in order for the coupling to the guided modes to acquire a chiral character. Note that, in the non-chiral case, i.e. when $\Delta\gamma=0$, Eq. (3) with $\Gamma=0$ describes a fully symmetric Dicke model Rod .

Photon detection and cumulant generating function. The statistics of the photon count at the two detectors (see Fig. 1) carries information about the state of the atom chain as well as its dynamical behavior. The counting statistics is fully characterized by the cumulants (e.g. the mean value, the fluctuations, etc.), which are encoded in the scaled cumulant generating function (SCGF). We now briefly summarize how to obtain the SCGF in practice. We start with the density matrix $\rho^{(K_{L},K_{R})}(t)$, which represents the state of the atomic system at time $t$ conditioned to the detection of exactly $K_{L}$ and $K_{R}$ photons on the left and right detector, respectively. The corresponding photon detection probability is thus given by $P_{t}(K_{L},K_{R})=\textrm{Tr}\left[\rho^{(K_{L},K_{R})}(t)\right]$ and its generating function is calculated via the two-dimensional Laplace transform $Z_{t}(s_{L},s_{R})=\sum_{K_{L},K_{R}=0}^{\infty}P_{t}(K_{L},K_{R})e^{-K_{L}s_{L}}e^{-K_{R}s_{R}}$. In the limit of long times, i.e. times much longer than any microscopic time scale of the dynamics, it assumes a so-called large deviation form: $Z_{t}(s_{L},s_{R})\simeq e^{t\theta(s_{L},s_{R})}$. The SCGF is here given by $\theta(s_{L},s_{R})$, which also happens to be the largest real eigenvalue of the transformed master operator Garr

[TABLE]

This permits a direct calculation of the SCGF, whose derivatives with respect to $s_{L}$ and $s_{R}$, taken at $s_{L}=s_{R}=0$, yield the moments of the photon number distribution. For example, the average number of photons emitted into the left propagating mode is given by $k_{L}=-\partial_{s_{L}}\theta(s_{L},s_{R})\!\mid_{s_{L}=s_{R}=0}$, while the corresponding fluctuations are given by $\partial^{2}_{s_{L}}\theta(s_{L},s_{R})\!\mid_{s_{L}=s_{R}=0}$. Hence, a sudden change of the structure of the SCGF in the vicinity of $s_{L}=s_{R}=0$ is concomitant with large cumulants, e.g. strong fluctuations, which is a manifestation of non-trivial dynamical behavior Bark ; IgMad .

Intermittent entangled dark states. To be more concrete we consider a system of $N=6$ atoms in the chiral regime ($\Delta\gamma\neq 0$), described by the master equation (3) with $\Gamma=0$. All atoms are driven by a laser field with spatially uniform detunings $\delta_{j}=\delta$. To simplify things further, we count the total emission $K\!=\!K_{L}\!+\!K_{R}$ into both detectors, without resolving the left and right photons separately [$s_{L}\!=\!s_{R}\!=\!s$, $\theta(s_{L},s_{R})\rightarrow\theta(s)$]. The resulting SCGF is displayed in Fig. 2(a) as a function of the laser Rabi frequency $\Omega$. For sufficiently large $\Omega$, $\theta(s)$ shows a sudden change around $s=0$. This signals a sharp increase of the fluctuations in the number of photons arriving at the detectors per unit of time.

We show below that this is a consequence of the fact that the stationary state for large $\Omega$ is actually one where two states (or phases) with drastically different dynamical properties coexist. To analyze the nature of this phenomenon, let us come back to the generating function of the photon count distribution function:

[TABLE]

According to this formula, $Z_{t}(s)$ is given as the trace of a biased (unnormalized) density matrix $\rho_{s}$: For $s<0$, it is biased towards conditional states, $\rho^{(K)}$, with large numbers of photon count $K$, while for $s>0$ the states with low $K$ are the ones dominantly contributing to $\rho_{s}$. Note that $\rho_{s}$ can in practice be found as the eigenstate of the dynamical generator given by Eq. (4) corresponding to the largest eigenvalue $\theta(s)$.

The fact that the SCGF changes quickly in the vicinity of $s=0$, as seen in Fig. 2(a), therefore signals that the biased state $\rho_{s}$ changes drastically when $s$ is varied across zero. In the case at hand the steady state turns out to be characterized by the coexistence of a dark state with zero photon count $K$ (inactive phase) and a bright (active) phase with large $K$. This is seen in Fig. 2(b), where we show the (biased) mean number of emitted photons per unit time, $k(s)$ 222Note here that, for $s<0$, we find another crossover in $k(s)$: This is a finite-size effect as we confirmed by numerical analysis..

Analyzing the biased state $\rho_{s}$ for positive and negative values of $s$ provides direct access to the properties of the two phases that coexist in the steady state. The inactive phase ($\rho_{s}$ for $s>0$) is characterized by the absence of photon emission into the waveguide. It is expected to be a so-called dark state Hu , which is pure and in which ”photons are trapped” in the atomic system due to quantum interference such that the steady state is virtually uncoupled from the modes of the waveguide. Indeed, calculating the purity of the biased state, $\mathrm{Tr}\left(\rho^{2}_{s}\right)$ [see Fig. 2(c)], shows that for large $\Omega$ and $s>0$ a pure state is assumed. Conversely, a mixed state is found for $s<0$. Interestingly, the dark state is entangled in a very peculiar way: In Fig. 2(d) we show the concurrence Woot ${\cal C}\left(\rho_{s}^{(AB)}\right)$ of the two-atom reduced (and biased) density matrix $\rho_{s}^{(AB)}$ for the $AB$-pairs shown above the panel. Note, that other pairs are not entangled, i.e. this inactive phase is one of highly entangled dimers.

Fig. 2 nicely illustrates the strength of the applied analysis method: large fluctuations, i.e. a steep change of the photon count rate $k(s)$ at $s=0$, are indeed the consequence of two coexisting phases in the steady state. Experimentally, this manifests itself in a highly intermittent photon count, as shown in Fig. 2(e) (middle panel). However, the theoretical investigation of the biased states $\rho_{s}$ reveals much more. It shows that with increasing $\Omega$ the steady state of the chiral atom chain is at first unimodal and enters a mixed entangled phase from an unentangled pure (product) state. Increasing $\Omega$ further, the steady state acquires an at first highly suppressed (only appearing at large negative $s$) component of a mixed phase. The weight of this phase increases until it coexists with the entangled phase, which has in turn become a pure dark state.

We now turn briefly to a case which was recently investigated in Ref. Pich . Here the atoms are excited with a laser field whose detuning is alternating, $\delta_{j}=(-1)^{j}\delta$. This may be challenging to achieve experimentally but has the advantage that the steady state can be calculated analytically: it is a dark state formed by dimers, which takes the form $\rho_{\mathrm{dimer}}=\ket{\psi}\bra{\psi}$, where $\ket{\psi}=\bigotimes_{j=1}^{N/2}\ket{D}_{2j-1,2j}$ and $\ket{D}\propto\ket{gg}+\beta\left(\ket{ge}-\ket{eg}\right)$ with $\beta=-2\Omega/(2\delta+i\Delta\gamma)$. The corresponding SCGF, the biased photon emission rate as well as the purity of the biased state are shown in Fig. 3. The plots are similar to the situation depicted in Fig. 2, i.e. there is an inactive (dark state) phase and an active (mixed) phase that is closing in as $\Omega$ is increased. However, unlike in the previous case, both phases cannot be brought into coexistence, i.e. the mixed phase always remains strongly suppressed. For an initial state in which all atoms are in their ground state $\left|g\right>$ it is only visible in the transient behavior of the photo current, as shown in Fig. 3(d).

Emission into unguided modes. In an experiment, photons are inevitably emitted also into unguided modes Gob at a finite rate $\Gamma$. We investigate this effect in Fig. 4 for the two detuning patterns discussed in Fig. 2 [$\delta_{j}=\delta$] and Fig. 3 [$\delta_{j}=(-1)^{j}\delta$], respectively. We find that an increasing value of $\Gamma$ generally smoothes the sharp features in the (biased) emission rate $k(s)$ and that the mixedness of the steady state increases. While the behavior of the case with homogeneous detuning [Figs. 4(a,b)] is only qualitatively affected, the situation changes quite drastically in the alternating case [Figs. 4(c,d)]. Here we see that a small $\Gamma$ removes the exact dark steady state and brings it into coexistence with a mixed phase. The result is an intermittent photon emission dynamics which is no longer transient, as is shown in the inset of Fig. 4(d). I.e. just like in Fig. 2 the pure dimer dark state is probed dynamically through fluctuations while the stationary state becomes mixed.

Finally we remark, that the coupling strength between the atoms (or generically the quantum emitters) and the guided modes can be controlled by employing a cavity, which is coupled itself to the guided modes. In this way, as showed in Tiec , it is possible to engineer a coupling to the unguided modes of one order of magnitude smaller than the one to the guided ones.

Conclusion. The counting statistics of photons emitted into guided modes yields direct insights into the steady state dynamics of chiral atom chains. This connection does not only allow to infer the emergence of dark state phases, but — more importantly — provides also in-situ access to dynamical properties and fluctuations. We demonstrated this by focussing on a regime of coexistence between an entangled dark state and a mixed bright phase. Beyond its appeal from an experimental perspective, our approach also shows the conceptual strength of the SCGF method and the usefulness of the correspondingly biased quantum state for unravelling complex dynamical behavior in open quantum systems. In the future it would be interesting to explore other photon detection protocols, e.g. homodyne detection, and the possibility to directly create the biased quantum state via post-selection Gamm ; Tan or a ”deformed dynamics” Fed . The latter is particularly relevant when considering the emission into unguided modes, where dark entangled phases merely emerge as a fluctuation in the steady state dynamics.