PHYSICAL REVIEW A 95, 033818 (2017)
Atom-light interactions in quasi-one-dimensional nanostructures: A Green’s-function perspective
A. Asenjo-Garcia,1,2,* J. D. Hood,1,2 D. E. Chang,3 and H. J. Kimble1,2 1Norman Bridge Laboratory of Physics MC12-33, California Institute of Technology, Pasadena, California 91125, USA 2Institute for Quantum Information and Matter, California Institute of Technology, Pasadena, California 91125, USA 3ICFO-Institut de Ciencies Fotoniques, Barcelona Institute of Science and Technology, 08860 Castelldefels (Barcelona), Spain (Received 2 November 2016; published 17 March 2017)
Based on a formalism that describes atom-light interactions in terms of the classical electromagnetic Green’s function, we study the optical response of atoms and other quantum emitters coupled to one-dimensional photonic structures, such as cavities, waveguides, and photonic crystals. We demonstrate a clear mapping between the transmission spectra and the local Green’s function, identifying signatures of dispersive and dissipative interactions between atoms. We also demonstrate the applicability of our analysis to problems involving three-level atoms, such as electromagnetically induced transparency. Finally we examine recent experiments, and anticipate future observations of atom-atom interactions in photonic band gaps.
DOI: 10.1103/PhysRevA.95.033818
I. INTRODUCTION have been interfaced with photonic crystal waveguides. Most of these experiments have been performed in conditions where As already noticed by Purcell in the first half of the past the resonance frequency of the emitter lies outside the band century, the decay rate of an atom can be either diminished gap. However, very recently, the first experiments of atoms [38] or enhanced by tailoring its dielectric environment [1–3]. and superconducting qubits [39] interacting with evanescent Likewise, by placing more than one atom in the vicinity of modes in the band gap of photonic crystal waveguides have photonic nanostructures, one can curtail or accelerate their been reported. collective decay. In addition to modifying the radiative decay, Within this context, it has become a necessity to understand nanophotonic structures can be employed to spatially and the rich spectral signatures of atom-like emitters interacting spectrally engineer atom-light interactions, thus obtaining through the guided modes of quasi-one-dimensional nanopho- fundamentally different atom dynamics from those observed tonic structures within a unified framework that extends in free space [4]. beyond those of cavity [40] or waveguide QED [41]. In In the past decade, atoms and other quantum emitters this work, we employ a formalism based on the classical have been interfaced with the electromagnetic fields of a electromagnetic Green’s function [42–46] to characterize the plethora of quasi-1D nanostructured reservoirs, ranging from response of atoms that interact by emitting and absorbing high-quality optical [5–10] and microwave [11,12] cavities photons through the guided mode of the nanostructure. Since to dielectric [13–20], metallic [21–24], and superconducting the fields in the vicinity of the structure might have complex [25,26] waveguides. Photonic crystal waveguides, periodic spatial and polarization patterns, the full Green’s function dielectric structures that display a band gap where light is only known analytically for a handful of systems (such propagation is forbidden [27,28], have been proposed as as planar multilayer stacks [47], infinite nanofibers [48–50], promising candidates to study long- and tunable-range co- and a few more [44]) and beyond that one has to resort herent interactions between quantum emitters [29–32]. Due to numerical solvers of Maxwell’s equations. However, in to the different character of the guided modes at various quasi-1D nanostructures, one can isolate the most relevant frequencies within the band structure of the photonic crystal, guided mode and build a simple prescription for the 1D Green’s the interaction of the quantum emitters with the nanostructure function that accounts for the behavior of this mode, greatly can be remarkably distinct depending on the emitter resonance simplifying the problem. frequency. Far away from the band gap, where light propagates, In the first part of the article, we summarize the procedure the guided modes resemble those of a conventional waveguide. to obtain an effective atom-atom Hamiltonian, in which the Close to the band gap, but still in the propagating region, guided-mode fields are effectively eliminated and the atom the fields are similar to those of a quasi-1D cavity, whereas interactions are written in terms of Green’s functions [42–44]. inside the band gap the fields become evanescent, decaying We then apply this formalism to a collection of atoms in exponentially. different quasi-one-dimensional dielectric environments, and All these regimes have been recently explored in the analyze the atomic transmission and reflection spectra in terms laboratory, where atoms [33–35] and quantum dots [4,36,37] of the eigenvalues of the matrix consisting of the Green’s functions between every pair of atoms. We show that, in the linear (low-saturation) regime, asymmetry in the transmis- *[email protected] sion spectra and frequency shifts are signatures of coherent atom-light interactions, whereas symmetric line shapes reveal Published by the American Physical Society under the terms of the dissipation. We also tackle the problem of three-level atoms Creative Commons Attribution 4.0 International license. Further coupled to quasi-1D nanostructures under electromagnetically distribution of this work must maintain attribution to the author(s) induced transparency (EIT) conditions, and derive analytical and the published article’s title, journal citation, and DOI. expressions for the polaritonic band structure of such systems
2469-9926/2017/95(3)/033818(16) 033818-1 Published by the American Physical Society ASENJO-GARCIA, HOOD, CHANG, AND KIMBLE PHYSICAL REVIEW A 95, 033818 (2017) in terms of the eigenvalues of the Green’s-function matrix. does in principle cover the interesting case of chiral atom- Finally, based on the rapid technical advances in fabrication of light interactions [54,55], where the time-reversal symmetry both optical and microwave structures, we project observable breaking is due to the atomic states. signatures that can be made in the next generation of In analogy to its classical counterpart, the electric field experiments of atoms and superconducting qubits interacting operator at frequency ω can be written in terms of bosonic in the band gap of photonic crystal waveguides. annihilation (creation) operators fˆ (fˆ†)as[42] 2 h ¯ 0 II. ATOM-LIGHT INTERACTIONS IN TERMS Eˆ (r,ω) = iμ0 ω dr Im{ (r ,ω)} G(r,r ω) OF GREEN’S FUNCTIONS π · ˆ + = ˆ + + ˆ − Much effort has gone into developing a quantum formalism f(r ,ω) H.c. E (r,ω) E (r,ω), (2) to describe atoms coupled to radiation. A conventional tech- where Eˆ +(−)(r,ω) is the positive (negative) frequency compo- nique is to express the field in terms of a set of eigenmodes nent of the field operator, H.c. stands for Hermitian conjugate, of the system, with corresponding creation and annihilation † and the total field operator reads Eˆ (r) = dωEˆ (r,ω). Within operators a and a [41]. This canonical quantization technique ˆ is well suited for approximately closed systems such as high-Q this quantization framework, f(r,ω) is associated with the cavities and homogeneous structures such as waveguides, both degrees of freedom of local material polarization noise, which { } of which have simple eigenmode decompositions. However, accompanies the material dissipation Im (r,ω) as required the application of this quantization scheme to more involved by the fluctuation-dissipation theorem [46]. This expression nanostructures is not straight forward. Further, the formalism guarantees the fulfillment of the canonical field commutation is not suited for dispersive and absorbing media as the relations, even in the presence of material loss. The appearance commutation relations for the field operators are not conserved of the Green’s function reveals that the quantumness of the [51]. system is encoded in either the correlations of the noise ˆ Instead, here we describe atom-light interactions using a operators f or in any other quantum sources (such as atoms), quantization scheme based on the classical electromagnetic but the field propagation obeys the wave equation and as such Green’s function, valid for any medium characterized by the spatial profile of the photons is determined by the classical a linear and isotropic dielectric function (r,ω), closely propagator. following the work of Welsch and colleagues [42–44,46]. In We now want to investigate the evolution of N identical two- the following, we employ this formalism to derive an atom- level atoms of resonance frequency ωA that interact through a atom Hamiltonian in which the field is effectively eliminated, guided-mode probe field of frequency ωp. Within the Born- yielding an expression that only depends on atomic operators. Markov approximation, we trace out the photonic degrees Moreover, once the dynamics of the atoms is solved, the of freedom, obtaining an effective atom-atom Hamiltonian electric field at every point along the quasi-one-dimensional [43,56,57]. This approximation is valid when the atomic cor- structure can be recovered through an expression that relates relations decay much slower than the photon bath correlations, the field to the atomic operators. or, in other words, when the Green’s function is characterized Classically, the field E(r,ω) at a point r due to a source by a broad spectrum, which can be considered to be flat over the current j(r ,ω)atr is obtained by means of the propagator atomic linewidth. Then, the atomic density matrixρ ˆA evolves ρ˙ =−i/h H ,ρ + L ρ of the electromagnetic field, the dyadic Green’s function (or according to ˆA ( ¯)[ ˆA] [ˆA][40]. Within the rotating wave approximation, and in the frame rotating with the Green’s tensor), as E(r,ω) = iμ0ω dr G(r,r ,ω) · j(r ,ω). probe field frequency, the Hamiltonian and Lindblad operators In particular, for a dipole source p located at r0, the current read is j(r,ω) =−iωp δ(r − r0), and the field reads E(r,ω) = 2 μ0ω G(r,r0,ω) · p. The tensorial structure of the Green’s function accounts for the vectorial nature of the electromag- N N H =− i − ij i j netic field, as a dipole directed along the xˆ direction can create h¯ A σˆee h¯ J σˆegσˆge a field polarized not only along xˆ, but also along yˆ and zˆ. i=1 i,j=1 Throughout this paper, the Green’s tensor will be also denoted N − i ∗ + i as the Green’s function. − d · Eˆ (ri )ˆσ + d · Eˆ (ri )ˆσ , (3a) p ge p eg The Green’s function G(r,r ,ω) is the fundamental solution i=1 of the electromagnetic wave equation and obeys [52] N ij L = i j − i j − i j 2 [ˆρA] 2ˆσgeρˆAσˆeg σˆegσˆgeρˆA ρˆAσˆegσˆge , ω 2 ∇ × ∇ × G(r,r ,ω) − (r,ω) G(r,r ,ω) = δ(r − r )1, i,j=1 c2 (1) (3b) ˆ = − where (r,ω) is the medium relative permittivity. For a scalar where Ep is the guided-mode probe field, and A ωp permittivity, Lorentz reciprocity holds and, then, GT(r,r ,ω) = ωA is the detuning between the guided-mode probe field and G(r ,r,ω), where T stands for transpose (and operates on the the atom. The dipole moment operator is expressed in terms = ∗ j + j polarization indexes). This formalism ignores the possibility of the dipole matrix elements as pˆ j d σˆeg d σˆge, where j of systems made nonreciprocal by the material response σˆeg =|eg| is the atomic coherence operator between the = | | [53] (e.g., nonsymmetric permittivity tensors), although it ground and excited states of atom j, and d g pˆ j e is the
033818-2 ATOM-LIGHT INTERACTIONS IN QUASI-ONE- . . . PHYSICAL REVIEW A 95, 033818 (2017) dipole matrix element associated with that transition. The spin- This expression can be understood as a generalized input- exchange and decay rates are output equation, where the total guided-mode field is the ˆ + sum of the probe, i.e., free, field Ep (r) and the field ij = 2 ∗ · · J μ0ωp/h¯ d Re G(ri ,rj ,ωp) d, (4a) rescattered by the atoms. The quantum nature of these equations has been treated before when deriving a gener- ij = 2 ∗ · · 2μ0 ωp/h¯ d Im G(ri ,rj ,ωp) d. (4b) alized input-output formalism for unstructured waveguides Note that the dispersive and dissipative atom-atom cou- [58,59]. plings are given in terms of the total Green’s function of the medium. For a given dielectric geometry, G(ri ,rj ,ωp) can be calculated either numerically or analytically to obtain III. TRANSMISSION AND REFLECTION quantitative predictions for the spin exchange and decay matrix IN QUASI-1D SYSTEMS elements. On the other hand, given some basic assumptions A. Atomic coherences in the low-saturation regime about the regime of interest, one can construct a simple We now explore the behavior of the atoms under a effective model for G(ri ,rj ,ωp). This enables one to broadly coherent, continuous-wave probe field. In the single-excitation capture a number of physical systems and gain general insight, manifold and low-saturation (linear) regime (σˆee=0), the and we take this approach here. atoms behave as classical dipoles. Then, the Heisenberg In particular, we assume that there is a single 1D equations for the expectation value of the atomic coher- guided band to which the atoms predominantly couple, ences (σˆ =σ ) are linear on the atomic operators, and and explicitly separate its contribution from the Green’s eg eg read function, G(ri ,rj ,ωp) = G1D(ri ,rj ,ωp) + G (ri ,rj ,ωp). The second term in the sum contains the atom-atom interactions N mediated by all other free-space and guided modes. We further σ˙ i = i + i σ i + i + i g σ j , (7) ge A 2 ge i ij ge assume that only collective interactions via the explicitly j=1 separated 1D channel are important, while G (ri ,rj ,ωp) only = ∗ · + provides independent single-atom decay and energy shifts. where i d Ep (ri )/h¯ is the guided-mode Rabi frequency ij = ij + ij = ij + Then, we can write J J1D J δij and 1D δij , (with Ep =Eˆ p), and where δij is the Kronecker delta. In particular, in free 2 ∗ space, is simply 0 = (2μ0 ω /h¯) d · Im G0(ri ,ri ,ωp) · ij ij 2 ∗ p gij = J + i /2 = μ0ω h¯ d · G1D(ri ,rj ,ωp) · d = 3| |2 3 1D 1D p d ωp d /3πh ¯ 0c , where G0 is the vacuum’s Green’s function [i.e., the solution to Eq. (1) when (r,ω) = 1]. depends only on the Green’s function of the guided mode. Depending on the geometry and dielectric response of the For long times, the coherences will damp out to a steady state nanostructure, and on the atom position, can be larger (˙σ i = 0). The solution for the atomic coherences is then ge or smaller than 0. J accounts for frequency shifts due to other guided and nonguided modes, and is in general spatially σ =−M −1 with M = ( + i /2)1 + g. (8) dependent. For an atom placed in free space, J is the Lamb ge A shift, which renormalizes the atomic resonance frequency, In the above equation, σ = (σ 1 ,...,σN ) and = and is considered to be included in the definition of the free ge ge ge ( ,..., ) are vectors of N components, and M is an space ω . The presence of a nanostructure shifts the atomic 1 N A N × N matrix that includes the dipole-projected matrix g resonance frequency. We will for simplicity consider this shift of elements g . Significantly, the matrix is not Hermitian, identical for every atom and assume its value to be zero. ij as there is radiation loss. However, due to reciprocity, the Once the dynamics of the atomic coherences are solved for, Green’s-function matrix is complex symmetric [GT(r,r ,ω) = one can reconstruct the field at any point in space. Generalizing G(r ,r,ω)], and g inherits this property if the dipole matrix Eq. (6.16) of Ref. [43] for more than a single atom, the elements are real, which will be a condition enforced from now evolution of the bosonic field operator is given by = on. Complex symmetric matrices can be diagonalized, gvξ 2 λ v with ξ = 1 ...N, where λ and v are the eigenvalues ˙ ω 1 ξ ξ ξ ξ fˆ (r,ω) =−iωfˆ(r,ω) + Im{ (r,ω)} and eigenvectors of g, respectively. Since the first term of M c2 πh ¯ 0 is proportional to the identity, M and g share the same set of N eigenvectors. × ∗ · j G (r,rj ,ω) d σˆge, (5) The eigenmodes represent the spatial profile of the collec- j=1 tive atomic excitation, i.e., the dipole amplitude and phase at each atom. However, as the matrix g is non-Hermitian, where the atoms act as sources for the bosonic fields. We can the eigenmodes are not orthonormal in the regular sense, formally integrate this expression and plug it into the equation but instead follow different orthogonality and completeness for the field [Eq. (2)]. After some algebra, and performing N T · = ⊗ T = Markov’s approximation, we arrive at the final expression for prescriptions, namely vξ vξ δξ,ξ and ξ=1 vξ vξ 1, the field operator, which is simply where T indicates transpose instead of the customary con- jugate transpose [60]. After inserting the completeness re- N lation into Eq. (8), we find that the expected value of ˆ + = ˆ + + 2 · j E (r) Ep (r) μ0ωp G(r,rj ,ωp) d σˆge. (6) the atomic coherences in the steady state in terms of j=1 the eigenvalues and eigenvectors of the quasi-1D Green’s
033818-3 ASENJO-GARCIA, HOOD, CHANG, AND KIMBLE PHYSICAL REVIEW A 95, 033818 (2017) function is are T · vξ N σ =− v , (9) t(A) 1 ge ( + J ) + i( + )/2 ξ = 1 − ξ∈mode A ξ,1D ξ,1D t ( ) g (r ,r ) 0 A ββ right left ξ=1 where Jξ,1D = Re λξ and ξ,1D = 2Imλξ are the frequency T T g (rright) · vξ v · gαβ (rleft) shifts and decay rates corresponding to mode ξ, and the sum × αβ ξ , (11a) ( + J ) + i( + )/2 is performed over mode number from 1 to N. The scalar A ξ,1D ξ,1D T · = N N product in the numerator vξ j=1 vξ,j j describes the 1 coupling between the probe field and a particular collective r( ) = r ( ) − A 0 A g (r ,r ) atomic mode. Both the frequency shifts and decay rates, as ββ left left ξ=1 well as the spatial profile eigenstates of g, are frequency T T g (rleft) · vξ v · gαβ (rleft) dependent. × αβ ξ , + + + (11b) The dynamics of the atoms can be understood in terms of (A Jξ,1D) i( ξ,1D)/2 the eigenmodes of g, where the real and imaginary parts of the where t0(A) and r0(A) are the transmission and reflection eigenvalues correspond to cooperative frequency shifts and de- coefficients for the 1D photonic structure when no atoms are { } cay rates of the collective atomic modes ξ . As the modes are present. non-normal, the observables cannot be expressed as the sum over all different mode contributions but, instead, any mea- surable quantity will show signatures of interference between C. Simplified expression for the transmission different modes. Although it could be considered a mathemat- For linearly polarized, transverse guided modes with an ical detail, the fact that the modes of a system are non-normal approximately uniform transverse field distribution, one can has deep physical consequences. For instance, non-normal further simplify the expression for the transmission coefficient. dynamics is responsible for phenomena as different as the We find a product equation that only depends on the eigenval- Petermann excess-noise factor observed in lasers [61–63]or ues of the Green’s-function matrix g, and not on their spatial the transient growth of the shaking of a building after an structure (i.e., the eigenfunctions). Following Appendix B,we earthquake [64]. obtain N + t(A) = A i /2 B. Transmission and reflection coefficients t ( ) ( + J ) + i( + )/2 0 A ξ=1 A ξ,1D ξ,1D Having previously calculated the linear response of an ensemble of atoms to an input field, we now relate the response N ≡ to observable outputs, i.e., the reflected and transmitted fields. tξ (A). (12) One can calculate the total field from Eq. (6), by substituting ξ=1 in the solution of Eq. (9) for the atomic coherences σge. The total transmission coefficient can thus be written = For a dipole moment directed along α (i.e., d dαˆ ), the as the product of the transmission coefficients of each of β-polarization component of the field reads the collective atomic modes. Noticeably, when looking at N T T + the transmission spectrum of atoms that interact through g (r) · vξ v · E E+(r) = E+ (r) − αβ ξ p,α , the guided mode of a quasi-1D nanostructure, there is a β p,β ( + J ) + i( + )/2 ξ=1 A ξ,1D ξ,1D redundancy between the eigenfunctions and eigenvalues, and one is able to obtain an expression that does not depend on the (10) former (i.e., all the relevant information about the geometry is + contained in the collective frequency shifts and decay rates). where the j component of the electric field vector Ep,α, jj + + In particular, for a single atom located at x with J ≡ J which reads E = αˆ · E (r ), no longer represents dif- j 1D 1D p,α,j p j jj ≡ ferent polarization components, but the dipole-projected and 1D 1D, the eigenvalues are directly proportional to field evaluated at the atoms’ positions rj .Thej compo- the local Green’s function, and nent of vector g (r)isg (r) = g (r,r ) = (μ ω2d2/h¯)ˆα · + αβ αβ,j αβ j 0 p t(A) = A i /2 G r r · ˆ . (13) 1D( , j ,ωp) β, where j runs over the atom number. In t0(A) (A + J1D) + i( + 1D)/2 particular, the scalar product gT (r) · v = N g (r)v αβ ξ i=j βα,j ξ,j =||2 represents how much the mode ξ contributes to the field The transmittance T t can be recast into a Fano-like line emitted by the atoms. shape [65]as We now assume that the atomic chain and the main T (q + χ)2 2 1 axis of the nanostructure are oriented along xˆ.Inor- = + , (14) T 1 + χ 2 + 1 + χ 2 der to connect the above expression to the transmission 0 1D + = + + =− + and reflection coefficients, we evaluate the field Eβ (r) where χ 2(A J1D)/( 1D ) and q 2J1D/( 1D at the positions r = rright and r = rleft, which are con- ) is the so-called asymmetry parameter. For sidered to be immediately outside the atomic chain, and 1D, the second term is negligible and the normal- only differ in the x component. Following Appendix A, ized transmittance is a pure Fano resonance, with q = the normalized transmission and reflection coefficients −Re{G1D(rj ,rj ,ωp)}/Im{G1D(rj ,rj ,ωp)}. Fano resonances
033818-4 ATOM-LIGHT INTERACTIONS IN QUASI-ONE- . . . PHYSICAL REVIEW A 95, 033818 (2017)
(a) (b) e 5 J1D=-2 1D g 1D 4
0 3
T/T 2 J =- J =0 1D 1D 1 1D 0 -10 -5010 5
A / ’
FIG. 1. (a) Sketch of an atom interacting with the guided mode of a structured 1D nanostructure. The single-atom decay rate is 1D,and 2 the decay into nonguided modes is characterized by . (b) Normalized transmission spectra (|t/t0| ) for a single atom for different values of the ratio between the real and imaginary parts of the guided-mode Green’s function, following Eq. (13). The decay rate into the guided modes is taken to be 1D = for all cases. arise whenever there is interference between two different of length L and effective area A, the JC Hamiltonian and its transport channels. For instance, in a cavity far from resonance, corresponding Lindblad operator read there is interference arising from all the possible optical paths that contribute to the transmission signal due to reflec- N N H =− † − i + † i + i tions at the mirrors, whereas in an unstructured waveguide h¯ caˆ aˆ h¯ A σˆee h¯ i aˆ σˆge σˆegaˆ there is no such interference and thus the line shape is i=1 i=1 Lorentzian. +hη¯ (aˆ + aˆ †), (15a) For a single atom, there is a clear mapping between the N spectrum line shape and the local 1D Green’s function. For a L [ˆρ] = 2ˆσ i ρˆσˆ j − σˆ i σˆ j ρˆ − ρˆσˆ i σˆ j nanostructure with a purely imaginary self Green’s function 2 ge eg eg ge eg ge i,j=1 G(xi ,xi ) (such as a waveguide or a cavity at resonance), κ the spectrum is Lorentzian, and centered around the atomic + c (2aˆρˆaˆ † − aˆ †aˆρˆ − ρˆaˆ †aˆ), (15b) frequency. However, if the real part is finite, one would observe 2 a frequency shift of the spectrum, which becomes asymmetric. where aˆ is the cavity-field annihilation operator,ρ ˆ is the Figure 1(b) shows how the normalized transmission spectrum density matrix for the atoms and the cavity field, η is a for a single atom becomes more and more asymmetric for frequency that represents the amplitude of the classical driving higher ratios J1D/ 1D. Also, there is an appreciable blueshift field, c = ωp − ωc is the detuning between the driving of the spectral features. (probe) and the cavity fields, and κc is the cavity-field decay. = = We would like to remark that the Markov approximation The√ atom-cavity coupling is i cos(kcxi ), where has thus far been employed in our analysis, as every Green’s d ωc/(¯h 0LA) is modulated by a function that depends on the function is considered to be a complex constant over frequency atoms’ positions and the cavity wave vector kc. The Heisenberg ranges larger than the linewidth of the atoms. However, the equations of motion for the field and atomic operators are expressions for the transmission and reflection coefficients are valid also in the non-Markovian regime. We analyze this issue κ N ˙ = − c − i − in more detail in Appendix C. aˆ ic aˆ i i σˆge iη, (16a) 2 = i 1 i i i i IV. APPLICATION TO SEVERAL ONE-DIMENSIONAL σˆ˙ = iA − σˆ + i i σˆ − σˆ a.ˆ (16b) ge 2 ge ee gg PHOTONIC STRUCTURES When κ and < min{ ,κ }, the cavity field can be In this section, we analyze the transmission spectra of c c c adiabatically eliminated, and the field operator reexpressed in atoms placed along common quasi-1D nanostructures, such terms of the atomic ones, i.e., as cavities, waveguides, and photonic crystals. N 1 i aˆ˙ = 0 → aˆ = η + i σˆ . + i κc ge A. Standing-wave cavities c 2 i=1 To begin with, we want to illustrate the connection between Introducing this expression back into the equation for the the Green’s-function formalism and the well-known Jaynes- atomic operator, one can deduce a master equation for the Cummings (JC) model [40,66]. For N atoms in a driven cavity atomic density matrixρ ˆA. The new Hamiltonian and Lindblad
033818-5 ASENJO-GARCIA, HOOD, CHANG, AND KIMBLE PHYSICAL REVIEW A 95, 033818 (2017) operators read just as those of Eqs. (3a) and (3b), but for a field. Depending on the detuning between the probe field and classical driving field, and with spin exchange and decay rates the cavity resonance, g(ωp) can be purely imaginary, yielding into the cavity mode given by [67] dissipative atom-atom interactions, or can have both real and imaginary parts, resulting in both dissipative and dispersive 2 J ij =− c cos(k x ) cos(k x ), (17a) couplings. 1D 2 + 2 c i c j c κc 4 The matrix g is separable (has rank one) as it can be written as the tensor product of just one vector by 2κ ij = c cos(k x ) cos(k x ). (17b) itself. The matrix has one eigenstate describing a super- 1D 2 + 2 c i c j c κc 4 position of atomic coherences that couples to the cav- = N ii = It can thus be seen that the Markovian approximation to ity (a “bright mode”), with eigenvalue λB i=1 g max + max N 2 arrive at these equations is equivalent to the absence of (J1D i 1D /2) i=1 cos (kcxi ). This atomic collective ex- strong-coupling effects within the JC model. citation follows spatially the mode profile of the cavity, i.e., i ∝ − The last step for connecting this simple model with our σge cos(kcxi ). The matrix g has also N 1 decoupled formalism is to calculate the Green’s function of a cavity (“dark”) modes of eigenvalue 0. Because these dark modes ij ij have a zero decay rate into the cavity mode, it is also and confirm that J1D and 1D are precisely those obtained within the JC framework. The Green’s function of a quasi-1D impossible to excite them employing the cavity field. The cavity formed by partially transmitting mirrors of reflection optical response is thus entirely controlled by the bright mode, coefficient r (chosen to be real) is [68] and the transmission is simply 2 + ic | − | t(A) A i /2 ikp xi xj = . (23) G1D(xi ,xj ,ωp) [e N N − 2 2ikpL t ( ) + ii + + ii 2vgωpA(1 r e ) 0 A A i=1 J1D i i=1 1D /2 + + − + + reikp(L xi xj ) + reikp[L (xi xj )] Remarkably, this expression is valid no matter the separation −| − | between the atoms or whether they form an ordered or +r2eikp(2L xi xj )], (18) disordered chain. The transmission spectrum corresponds to where vg is the group velocity. For high-Q standing-wave that of a “superatom,” where the decay rates and the frequency cavities, i.e., with r 1, and choosing vg = c, the Green’s shifts are enhanced (N-fold if all the diagonal components of function can be approximated as g are equal) compared to those of a single atom. This result 2ic 1 replicates the well-known expressions for conventional cavity G (x ,x ,ω ) 1D i j p 2 2ik L QED. ωpA 1 − r e p × cos(k x ) cos(k x ). (19) p i p j B. Unstructured waveguides The cavity is resonant at a frequency ωc with corresponding Another paradigm that has been investigated frequently is = wave vector kc, chosen to be such that kcL 2πm, with that of “waveguide QED” [41]. The simple model of such a m being an integer. Close to resonance, one can write system consists of a single guided mode with translational = + − 2 2ikpL kp kc δk, and assume that δkL 1. Then 1 r e invariance, and where the dispersion relation is well approx- 1 − r2 − 2ir2δkL, and the Green’s function is simply imated as linear around the atomic resonance frequency. In a 2 1D translationally invariant system, a source simply emits a − c cos(kcxi ) cos(kcxj ) G1D(xi ,xj ,ωp) , (20) plane wave whose phase at the detection point is proportional ωpLA c + iκc/2 to the distance of separation. Therefore, the elements of the 2 where κc = (1 − r )c/L is the cavity linewidth. Therefore, Green’s-function matrix g depend on the distance between the the atoms’ spin-exchange and decay rates are given by atoms, and read 2 2 1D | − | ij μ0ωpd c = ikp xi xj J = Re G (x ,x ,ω ) =− , gij i e . (24) 1D 1D i j p i j 2 + 2 2 h¯ c κc 4 Remarkably, the self Green’s function in a waveguide is purely (21a) imaginary. The coherent interactions between atom i and 2μ ω2d2 ij = 0 p = κc atom j are dictated by the Hamiltonian [given by Eq. (3a)], 1D Im G1D(xi ,xj ,ωp) i j , 2 + 2 and are proportional to Re{gij }=−( 1D/2) sin kp|xi − xj |, h¯ c κc 4 whereas the dissipation is given by the Lindblad operator (21b) [given by Eq. (3b)], which is proportional to Im{gij }= which is precisely what is obtained within the ( 1D/2) cos kp(xi − xj )[24,56]. It is thus clear that by care- Jaynes-Cummings model. Let us now look at the transmission fully tuning the distance between the emitters, one can engineer spectrum of N atoms in a cavity. As we have just demonstrated, fully dissipative interactions. If the atoms form a regular coefficients of the dipole-projected Green’s-function matrix g chain and are spaced by a distance d such that kpd = nπ, read where n is an integer number, the matrix g has only one nonzero eigenvalue λ = iN /2 associated with the bright g = g(ω ) cos(k x ) cos(k x ), (22) B 1D ij p c i c j atomic mode. This situation is analogous to the case of atoms = max + max max max where g(ωp) J1D i 1D /2, where J1D and 1D are the interacting in an on-resonance cavity. Therefore, there will not spin-exchange and decay rates at the antinode of the cavity be any collective frequency shift, and the line shape will be a
033818-6 ATOM-LIGHT INTERACTIONS IN QUASI-ONE- . . . PHYSICAL REVIEW A 95, 033818 (2017)
(a) 2
o
o o o
o o o o o o o o
o o o o o o o o
o o o o o o o o
o o o o
o o o o
o o o o
o
o
o o
o o
o o
o o
o o
o o o o
o o o o
o
o o o
o o
o o
o
o o o
o o o o
o o o o
o
o
o o
o o o o
o
o
o
o
o o o o
o o o o
o
o o o
o
o o o
o o
o
o o o o o o o
o o
o o o o
o o o o
o
o o
o
o o o o
o o o o o o o o o o o o o o o o
o o o o
o o o o
o o o o
o o o o o o o o
o o o o o o o o
o o o o
o
o
o o
o o
o o
o o o o
o o o o o o o o
o o
o o
o o
o o o
o o o o o o o
o o o
o
o
o o o o o o o
o
o
o
o
o o o o
o
o o o
o o o o
o o o o
o o o
1D o o o o o o o
o o
o o o
o
o o o o
o o o o o o o o
o
o o
o
o o o o o o o o o o o o o o o o
o o o o
o o o o
o o o o o o o o
o o o o
o o o o
o
o
o
o
o o o o o o o o o o o
o o o o o
o o o o
o o o o
o
o
o o
o
o o o o o o o o o o
o
o
o
o o
o o o o
o o
o
o
o o o o o o o o
o o o o
o
o o
o
o o o o
o o o o o o o o
o o o o
o
o o o o o o o
o
o
o o
o o o o o o o o
o o o o
o o o o
o o o o o
o o o
o
o o
o o o o o
o o o o o o o o
o o o o
o o o o o o o o
o
o o o
o o o o
o o o o
o o o
o o o o
o o o
o o o o o o o o o o o o o o o o o o o
o o o o
o
o o o
o
o o
o
o o o o o o o o
o o o o
o
o o
o o o o o
o o o o o o o o o o o o
o o o o
o o o o
o o o o
o o o o o
o o
o
o
o o
o
o o o o o o o o
o o o o
o o o o o
o o o
o o o o
o o o o o o o o
o
o o o o o o o
o o o o
o o o o
o o o
o
o o o o o o o
o
o o o o o o o o
o o o o
o o o o
o o o o o o o o
o
o o
o
o o o o
o o o o
o
o o
o o o o o
o o o o o o o o
o o o o
o o o o
o o o
o
o o o o o
o
o
o o o o o
o o
o o
o
o o o
o o o o
o o o
/ o
o
o
o o o o o o
o
o o
0 o
o o o o
o
o o o
o o o o
o
o
o o
o o o o
o o o o
o o o o
o o o o
o
o o o
o o o o o o o o o o o o
o o o o
o o o o o o o o
o o
o o o o o o o o o o
o o o o
o
o
o o o o o o
o o o
o o o o o
o o o o
o o o o
o o o o
o o o o o o
o o
o o o o o o o o o o o o
o o o o o o o o o o o o
o o o o o o o o o o o o o o o o
o o o o
o o o o
o o o
o
o
o o o o o o o o o o o
o o o 1D o o o o o o o o
o
o o o o o o o o o o o o
o
o o o
o o o
o
o o o o
o o o o
o o o o
o o o o o o o J o
o o o o o o o o
o o o o o o o o o o o o o o o o o o o o -2 0 0.5 1 1.5 2 (b) 6 oo ooo ooo ooo oo oo oo oo oo oo oo oo oo o o o o o o o o o o o o o o o o 4 o o o o o o o o 1D o o o o o o o o ooo ooo ooo ooo ooo ooo ooo ooo o oo oo ooooo oo oo o o oo oo ooooo oo oo o o oo oo ooooo oo oo o ooo oo ooooo oo oo o oo oooo oooo oo oo oooo oooo oo oo oooo oooo oo oo oooo oooo oo / oo oo oo oo oo oo oo oo oo oo oo oo oo oo oo oo 2 o o o oo oo o o o o o o oo oo o oo o o o oo oo o o o o o o oo oo o o o 1D o o o o o o o o o o o o o o o o o o o o o o o o o o o oooo o o o o o o oooo o o o o o o oooo o o o o o o oooo o o o o oo ooo oo o o oo ooo oo o o oo ooo oo o o oo ooo oo o o oo o oo oo o oo o o oo o oo oo o oo o o ooo oo oo o oo o o oo o oo oo o oo o o oo ooooo ooooo oo o o oo ooooo ooooo oo o o oo ooooo ooooo oo o o oo ooooo ooooo oo o 0 ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 010.5 1.5 2
d / p
FIG. 2. (a) Frequency shifts and (b) decay rates of the collective modes of a regular chain of 5 atoms placed along a waveguide normalized to the single-atom decay rate into the guided mode 1D, as a function of the distance d between the atoms in units of the probe wavelength.