arXiv:0808.2557v1 [quant-ph] 19 Aug 2008 ulse nie . . . online: Published . . . Received: 3 2 ihe .Hartmann J. Michael Cavity Coupled Arrays in Phenomena Many-Body Quantum 1 .Apiain 24 ...... 12 ...... 21 . . . 3 ...... 3 . 1 ...... [email protected] . . 16 . . e-mail: ...... Applications . . . . . Realisations . . for 9. . . Platforms . Arrays . Promising Cavity . . Coupled . . in . 11 . 8. Models . 4 . Spin . Effective ...... 7. . . Diagrams . . Hamiltonian Phase . Jaynes-Cummings . . . the . on 6. . . Based Arrays Cavity Models . . Coupled . in . . Model . 5. Bose-Hubbard . . The . . . 4. Arrays . Cavity . Coupled . . . 3. QED Cavity . . 2. Introduction 1. Contents words: Key 2018 25, October Reviews, Photonics & Laser PACS: OS h lcetLbrtr,Ipra olg odn Pr London, College Imperial Laboratory, Blackett The QOLS, nttt o ahmtclSine,Ipra olg Lon College Imperial Sciences, Mathematical for Institute ehiceUiesta ¨nhn hskDprmn I,Ja Department Universit¨at M¨unchen,Technische Physik fidvda ie n h cesblt finhomogeneous els. of accessibility the addressabil and full sites individual the of including simula physics, as many-body properties arraysquantum interesting Such Bose- offer cavities models. the coupled Heisenberg as of anisotropic such the systems, and such Hubbard in models many-body pa of matic realisation the ph on different proposals review several we in architectures, on structures the and atomic-like in micro-cavities to progress interacting coupling of experimental arrays recent of many- on fabrication quantum Based create systems. to body electrodynamics th quantum concerning cavity studies of theoretical recent rev this survey In we manner. controlled article more a in experi them of accessing possibility tally the and physics, only condensed-matter encountered in appearance originally the phenomena, by many-body characterised quantum physics, quantu in atomic regimes and physical pave new optics have of exploration years the past for the way in the gained systems optical and atomic Abstract 42.50.Dv,42.50.Pq,03.67.Bg,75.10.Jm h nraiglvlo xeietlcnrlover control experimental of level increasing The unu aybd oes oaios aiyQD htnb photon QED, cavity polaritons, models, many-body quantum 1 , 2 , 3 , ∗ ennoGSL Brand G.S.L. Fernando , osof tors ao ˜ radig- ysical e rnkSr,878Grhn,Germany Garching, 85748 Str., Franck mes mod- use e men- o,S72G ntdKingdom United 2PG, SW7 don, their iew 2 ity of m , d 3 neCnotRa,S72W ntdKingdom United 2BW, SW7 Road, Consort ince n atnB Plenio B. Martin and n ihlgt h hsc fsc ytm a usually can and effects systems collective 24 ignoring such by As of other understood quanta. well physics each very light the with . be and light, strongly . atoms interact with . few naturally and . a not of . do scale atoms . typ- regime, the . quantum on the . in ically matter . and light the . with of concerned interaction . are physics atomic . and optics . Quantum . . . . Introduction 1...... Summary 10. light: Crystal ylsrpmigo utbeaos(e)pae ntecaviti the in placed (red) atoms by suitable place of indu pumping in be laser can frozen by repulsion be Such might (bottom). repulsion (top) mutual cavities their of array an through lockade htn ylo)ta riaiymv freely move ordinarily that (yellow) Photons 2 , 3 oyih iewl epoie ytepublisher the by provided be will line Copyright oyih iewl epoie ytepublisher the by provided be will line Copyright ced es. 1 2 Michael J. Hartmann et. al.: Quantum Many-Body Phenomena in Coupled Cavity Arrays

treating the interactions perturbatively. This is a much sim- can be understood by looking at a few quantum levels. pler and cleaner situation than the one encountered in the In this context, arrays of interacting Josephson junctions condensed matter context. There, strong interactions among have been used to reproduce properties of bosonic parti- the basic constituents, such as nuclei and electrons, lead to cles. For example, the Mott insulator phase [14] and the the emergence of completely new physics when one con- superconductor-to-Mott insulator transition [15] have been siders a mesoscopic or macroscopic number of interacting observed in Josephson junctions arrays. particles. Thus, even though the fundamental interactions Finally, theoretical proposals for using trapped ions for between the constituent particles are usually known, it is realising many-body models have been put forwarded [16, challenging to fully describe the properties of such sys- 17,18,19] and a first benchmark experiment in this direc- tems. tion has been reported recently [20]. The higher level of isolation of quantum optical sys- This article is devoted to review theoretical proposals tems have permitted, over the past decades, the develop- for the use of the atom-light interaction in coupled micro- ment of experimental techniques for their manipulation cavity arrays to create strongly correlated many-bodymod- with an unprecedented level of control. Highlights include els. This research was initiated in [21] and shortly after in the inventionof trapping and cooling techniques[1], which [22,23] (see also [24]). Since then a substantial number of have found applications e.g. in spectroscopy and precision publications have extended these ideas in a variety of di- measurements and which culminated in the realisation of a rections. These include photonic Mott insulators [25,26], Bose-Einstein condensate (BEC) of neutral atoms [2]; the spin models [27,28,29,30,31,32], dynamical effects [33, manipulationof the state of individual atomic and photonic 34,35,36,37], multi component models [38], phase dia- systems at the quantum level [3,4] and its use in the study grams [39,40,41,42,43], alternative physical realisations of the foundations of quantum theory [4,5]; and more re- [44,45,46,47,48], experimentalsignatures [49,50] and quan- cently the realisation of primitives for the processing of tum informationapplications [51,52,53,54,55,56]. In fact, quantum information and for the implementation of quan- historically, some quantum information proposals [51] were tum computation in several different quantum atomic and presented before the theory for effective many-body mod- optical set-ups [6]. els [21,22,23]. This progress raised the possibility of engineering strong Recent experimentalprogress in the fabrication of micro- interactions in these systems and of realising such strongly cavity arrays and the realisation of the quantum regime correlated models in large scales. This is of interest from a in the interaction of atomic-like structures and quantised fundamentalpoint of view, as in this way we can form new electromagnetic modes inside those [57,58,59,60,61,62, phases of matter for atoms, photons, or even molecules 63,64] opened up the possibility of using them as quan- that do not exist outside the laboratory. Moreover, due to tum simulators of many-bodyphysics. The first motivation the high level of experimental control over them, it is pos- for this study is the desire to make photons or combined sible to study many-body physics in a much cleaner man- photonic-atomicexcitations enter new quantum states, which ner than it is possible employing condensed-matter sys- do not naturally appear in nature. The second is that this tems, in which the strong interaction with the environment set-up offer advantages over other proposals for the reali- and the very short space- and time-scales involved inhibit sation of strongly interacting many-body models in quan- a precise control over their static and dynamic properties. tum optical systems. Due to the small separations between The use of one quantum mechanical system to simulate neighbouring sites, it is experimentally very challenging to the physics of a different one is an idea that dates back access individual sites in Josephson junctions arrays and to Feynman [7] and has attracted considerable interest re- optical lattices. Neighbouring sites in coupled cavity ar- cently in the field of quantum information science [6]. rays, on the other hand, are usually separated by dozens A very successful example in this direction is the study of micrometres and can therefore be accurately accessed of cold trapped atoms in optical lattices [8]. By loading a by optical frequencies. This allows e.g. the measurement BEC of neutral atoms in a periodic optical lattice, formed of local properties as well as the exploration of inhomoge- by the off-resonant dipole interaction of atoms with over- neous systems. lapping standing waves of counter propagating lasers, the As the field of coupled cavity arrays is still quite young otherwise weak interactions of the atoms can be increased and many preprints are currently being reviewed, we have to such an extend that a strongly correlated many-body restricted ourselves here to only present details of work model is formed [9]. This idea led to several important that has been published in peer reviewed journals before experiments, such as the realisation of the Bose-Hubbard completion of the present text. This has been done in order model [10] and of a Tonks-Girardeau gas [11] and is cur- to avoid any interferences with the policies of journals that rently an active area of theoretical and experimental re- consider the work in question for publication. search [8,12]. The organisation of the paper is as follows: In section Also of note is the early study of strongly correlated 2 we discuss the use of cavities to enhance the interac- many-body models in Josephson junctions arrays [13]. Al- tion of light and matter and in section 3 we develop the though a Josephson junction is a mesoscopic solid state basic description for an array of interacting cavities. The device, it behaves in many respects just like the type of following section 4 reviews the theory for engineering ef- systems studied in quantum optics, as the relevant physics fective Bose-Hubbard Hamiltonians with polaritons and

Copyright line will be provided by the publisher PSfrag replacements γ

Laser & Photon. Rev. 1, No. 1 (2008) 3 g Ω ∆ rate γ and leaking out of photonsof the cavity mode at rate δ κ, see figure 1. The cavity decay rate κ is connected to the ω quality factor of the cavity Q by Q = ω /(2κ). Although ω C g other processes, such as thermal motion of the atom or 1 κ 2 dephasing due to background electromagnetic fields, also 3 contribute to the losses of the system, their effect are usu- 4 γ ally much smaller and can be disregarded in many appli- cations. The strong coupling regime of cavity QED is reached Figure 1 Cavity QED: A two level atom interacts via a dipole when the cooperativity factor, given by the the vacuum coupling g with the photons in the cavity. Excitations a lost via Rabi frequency to the square over the product of the spon- spontaneous emission at a rate γ and cavity decay at a rate κ. taneous emission and cavity decay rates,

g2 ξ := , (3) photons. In section 5 we discuss a second class of effec- 2κγ tive Hamiltonians that exhibits a phase diagram similar to is much larger than one: ξ 1. In this regime the coher- the Bose-Hubbard model: Jaynes-Cummings based many ≫ body models. The phase diagrams of both classes are then ent part of the evolution dominates over the decoherence reviewed in section 6. Effective spin lattice Hamiltonians processes and quantum dynamics of the joint atom cav- can also be generated in coupled cavity arrays and we re- ity mode system can be observed. Such a regime has first view these approaches in section 7. In section 8 we then been achieved in seminal experiments for microwave [4] discuss several experimental platforms which appear to be and later optical frequencies [65], both using a single atom suitable for a future implementationof the theory approach inside a Fabry-P´erot cavity formed between two minia- presented in this review. We finish with a summary in sec- ture spherical mirrors. In these works, the strong coupling tion 10. regime has been employed e.g. to create conditional quan- tum dynamics and entanglement between the photonic and atomic degrees of freedom [4,66] and to perform quantum non-demolition measurements [4]. Recently a strong cou- 2. Cavity QED pling regime with very high cooperativity factors has also been achieved with artificial atoms formed by supercon- Usually photons interact only over a very short time, and ducting structures that interact with the electromagnetic hence very weakly, with atoms. This interaction can be modes of stripline resonator [62]. substantially increased if the atoms and photonsare trapped Most importantly for our purposes the strong coupling inside a resonator. In the simplest situation, we can con- regime led to a realisation of the photon blockade effect sider the coupling of a single light mode with a single two [67], where, due to the strong interaction of the cavity level system inside a cavity. Using the dipole and the ro- mode with atoms, a single photon can modify the reso- tating wave approximations their dynamics is governed by nance frequency of the cavity mode in such a way that a the Jaynes-Cummings Hamiltonian, second photoncan not enter the cavity before the first leaks out. H = ~ω a†a + ~ω e e + ~g(a† g e + a e g ), (1) C 0| ih | | ih | | ih | In this review we are interested in the situation where several cavities operating in the strong coupling regime where ωC ,ω0 are the frequencies of the resonant mode of the cavity and of the atomic transition, respectively, g is are coupled to each other. It turns out that the Fabry-P´erot Jaynes-Cummings coupling between the cavity mode and architecture is not very suitable for the coupling of sepa- rate cavities. Nonetheless, several new cavity QED set-ups the two level system, a† is the creationoperatorof a photon in the resonant cavity mode, and g , e are the groundand have recently emerged in which (i) the strong coupling excited states of the two level system.| i | Thei Rabi frequency regime has been achieved and (ii) the construction of ar- g is given by rays of coupled cavities has been realised, or at least seem reasonable to be so. These new systems, which we briefly ω review in section 8, are the strongest candidates for realis- g = df(r ) C , (2) 0 2ǫ V ing the proposals that we review in this paper. r 0 mode

where Vmode is the mode volume of the cavity, f(r0) is the mode-function at the position r0 of the two level system, 3. Coupled Cavity Arrays d is the electric dipole transition moment, and ǫ0 is the free-space permittivity. In this section we present the theoretical description for In addition to the coherent interaction, there are two the coupling of cavities in an array. Photon hopping can main loss processes that affect the dynamics of the sys- occur between neighbouring cavities due to the overlap of tem: spontaneous emission from level e to level g at a the spacial profile of the cavity modes, see figure 2 . In | i | i

Copyright line will be provided by the publisher 4 Michael J. Hartmann et. al.: Quantum Many-Body Phenomena in Coupled Cavity Arrays

4. The Bose-Hubbard Model in Coupled Cavity Arrays

Although realising Hamiltonian (5) might be interesting for quantum information propagation [71,72], it is not the Figure 2 An array of coupled cavities. Photon hopping occurs type of model that one would like to quantum simulate. In- due to the overlap (shaded green) of the light modes (green lines) deed, the model is harmonic and can be very easily solved of adjacent cavities. in terms of Bloch waves, which allows for a simple under- standing of all its basic properties. The situation changes dramatically if we add an on-site interaction term. The in- order to model such a process, we follow references [21, terplay of tunnelling and interaction leads to interesting many-body physics. A paradigmatic model for such a sit- 68,69] and describe the array of cavities by a periodic di- uation is the Bose-Hubbard Hamiltonian, which reads electric constant, ǫ(r) = ǫ(r + Rn), where r is a given n three dimensional vector, R is a constant and labels tu- H = µ b† b J (b† b ′ + h.c.) BH k k − k k pels of integers. In the Coulomb gauge the electromag- k netic field is represented by a vector potential A satisfying X X (ǫ(r) A)=0. We can expand A in Wannier functions, + U bk† bk(bk† bk 1), (7) ∇ · R n − ωR, each localised in a single cavity at location = R . Xk We describe this single cavity by the dielectric function where b† creates a boson at site k, J is the hopping rate, U ǫR(r). Then, from Maxwell’s equations one can derive the k the on-site interaction strength, and µ the chemical poten- eigenvalue equation tial. This model, which was first discussed by Fisher et al r 2 in Ref. [73,74] as the bosonic counterpart of the Hubbard ǫR( ) ωC wR ( wR)=0 , (4) model [75], contains two different phases at zero tempera- c2 − ∇× ∇× ture that are separated by a quantumphase transition where 2 where the eigenvalue ωC is the square the resonance fre- the ration U/J crosses a critical value. Let us briefly de- quency of the cavity, which is independent of R due to the scribe the two phases qualitatively. periodicity. Assuming that the Wannier functions decay When J U, the tunnelling term dominates and the sufficiently fast outside the cavity, only Wannier modes of lowest energy≫ state of the system is a condensate of de- nearest neighbour cavities have a non-vanishing overlap. localised bosons: the system is in a superfluid phase. If Introducing the creation and annihilation operators of we start to increase a repulsive on-site interaction, adding the Wannier modes, aR† and aR, the Hamiltonian of the morethan one particleto onesite requires energythat grows field can be written as quadratically with the particle number. It will thus be en- ergetically favourable for the bosons to evenly distribute 1 = ω aR† aR + across all lattice sites and avoid moving around in order H C 2 R   to minimise the repulsion. Above a specific value of U/J, X the system ceases to be a superfluid and becomes a Mott +2ωCα aR† aR′ + h.c. . (5) insulator, with a well defined number of localised parti- R R′ < X, >   cles in each site1. The superfluid-to-Mott insulator phase Here is the sum of all pairs ofcavities whichare transition is of a quantum character and is driven by quan- nearest neighbours. Since α 1, we neglected counter tum fluctuations, as it occurs even at zero temperature. An rotatingP terms which contain≪ products of two creation or order parameter for such a transition is the variance of the 2 two annihilation operators of Wannier modes in deriving numberof bosonsin a given site, (bk† bk bk† bk ) .Inthe (5). The coupling parameter α is given by [69], superfluid phase, the number of bosonsh −h in eachi sitei fluctu- ates and thus the variance has a non-zerovalue. In the Mott 3 ⋆ α = d r (ǫR(r) ǫ(r)) wRwR′ ; R R′ = R, insulator phase, in turn, the variance is close to zero, as the − | − | Z bosons tend to be localised in each site. (6) As we saw in section 3, the coupling between neigh- and has been obtainednumerically for specific models[70]. bouring cavities naturally leads to the tunnelling term for In deriving Hamiltonian (5) we have assumed that all photons in arrays of coupled cavities. Hence, for a reali- the cavities have the same resonant frequency and that the sation of a Bose-Hubbard model a way to effectively gen- overlap of the Wannier modes α is the same for all cavity- erate an interaction term is needed. Indeed, up to harmonic cavity interactions. In practise of course there will always terms which can easily be compensated for, the Bose-Hubbard be some disorder in the array and the resonance frequen- interaction term has the form of a Kerr self non-linearity cies, ωC (R), and tunnelling rates, 2ωC(R)α((R) will dif- fer from cavity to cavity. The disorder in the array can even 1 In the case of incommensurate filling factors, i.e. if the num- give rise to interesting effects such as the emergence of ber of particles is not a multiple of the number of lattice sites, the glassy phases, see 6. system will always be in a superfluid phase.

Copyright line will be provided by the publisher Laser & Photon. Rev. 1, No. 1 (2008) 5 PSfrag replacements

4 ∆ 3 δ ωC g24

Ω g13 ωC ε 2

1 Figure 3 An array of cavities as described by our model. Photon hopping occurs due to the overlap (shaded green) of the light Figure 4 The level structure and the possible transitions of one modes (green lines) of adjacent cavities. Atoms in each cavity atom, ωC is the frequency of the cavity mode, Ω is the Rabi fre- (brown), which are driven by external lasers (blue) give rise to an quency of the driving by the laser, g13 and g24 are the parameters on site potential. Reprinted with permission from [21]. of the respective dipole couplings and δ, ∆ and ε are detunings. Reprinted with permission from [21].

2 2 (b† ) b . Although such an interaction naturally appears in k k the atoms in the cavity reads, some media as a result of a non-zero third order electric susceptibility, its effect is negligible on the level of indi- N I j j j vidual quanta, which is a reason for the difficulty of real- H = εσ22 + δσ33 + (∆ + ε)σ44 (8) ising nonlinear optics for individual quantum systems. Us- j=1 X   ing the enhanced light-matter interaction in cavity QED, it N is possible to engineer much stronger nonlinearities. Intu- j j j + Ω σ23 + g13 σ13 a† + g24 σ24 a† + h.c. , itively the strong interactions of the light mode with atoms j=1 inside the cavity, under particular circumstances, mediates X   strong nonlinear interaction among the photons of the cav- where σj = k l projects level l of atom j to level k kl | j ih j | ity mode. The strength of the nonlinearity can be increased of the same atom, ωC is the frequency of the cavity mode, even further if instead of consideringphotonsas the bosonic Ω is the Rabi frequency of the driving by the laser and g13 particles of the model, one considerspolaritons, joint photonic- and g24 are the parameters of the dipole coupling of the atomic excitations. In the next section we review the pro- cavity mode to the respective atomic transitions. posal for realising the Bose-Hubbard Hamiltonian using In a cavity array, N atoms in each cavity couple to the such polaritons. cavity mode via the interaction HI and photons tunnel be- tween neighbouring cavities as described in equation (5). Hence the full Hamiltonian describing this system reads

I 1 4.1. Polaritonic Bose-Hubbard model H = HR + ωC aR† aR + R R 2 X X   The polaritonic Bose-Hubbard model is hosted in an array +2ωCα aR† aR′ + h.c. (9) of cavities that are filled with atoms of a particular 4 level R R′ < , >   structure, which are driven with an external laser, see fig- X ures 3 and 4. Therebythe laser drivesthe atoms in thesame and can emulate a Bose-Hubbard model for polaritons as manner as in Electromagnetically Induced Transparency we will see in the following. (EIT) [76]: The transitions between levels 2 and 3 are cou- Assuming that all atoms interact in the same way with pled to the laser field and the transitions between levels the cavity mode, the description can be restricted to Dicke 2-4 and 1-3 couple via dipole moments to the cavity reso- type dressed states, in which the atomic excitations are de- nance mode. Levels 1 and 2 are assumed to be metastable localised among all the atoms. In the case where g24 = and their spontaneous emission rates are thus negligible. 0 and ε = 0 , level 4 of the atoms decouples from the dressed-state excitation manifolds [79]. The Hamiltonian It has been shown by Imam˘oglu and co-workers, that (8) can then be expressed in terms of the following cre- this atom cavity system can exhibit a very large nonlin- ation (and annihilation) operators: earity [77,78,79,80], and a similar nonlinearity has been observed experimentally [67]. 1 p† = gS† Ωa† , (10) Considering the level structure of figure 4, in a rotating 0 B 12 − frame with respect to 1   H0 = ωC a†a + 2 + 2 A δ N j j j p† = ΩS12† + ga† ± S13† , j=1 ωC σ22 + ωCσ33 +2ωCσ44 , the Hamiltonian
of ± A(A δ) ± 2 s ±   P   Copyright line will be provided by the publisher 6 Michael J. Hartmann et. al.: Quantum Many-Body Phenomena in Coupled Cavity Arrays

2 2 2 2 where g = √Ng13, B = g + Ω , A = √4B + δ , all terms that rotate at frequencies 2µ0 (µ+ µ ) or 1 N j 1 N j − ± − S† = σ and S† = σ . µ0 µ+ or µ0 µ can be neglected, which eliminates all 12 √N j=1 21 p13 √N j=1 31 ± ± − interactions that would couple p† and S† to the remaining The operators p†, p† and p† describe polaritons, quasi 0 14 P 0 + P polariton species. particles formed by combinations− of atom and photon ex- For g gΩ/B2 ∆ , the coupling to level 4 can be citations. In the subspace spanned by symmetric Dicke 24 treated perturbatively.| | ≪ This | | results in an energy shift of 2U states and in the limit of large atom numbers, N 1, ≫ with they satisfy bosonic commutation relations, g2 Ng2 Ω2 U = 24 13 (15) − ∆ 2 2 2 (Ng13 + Ω ) [p ,p ]=0 and p ,p† = δ for j, l =0, +, . (11) j l j l jl − h i and in an occupation probability of the state of one S14† excitation of 2U/∆, which determines an effective decay p0†, p+† and p† thus describe independent bosonic particles. − − † In terms of these polaritons, the Hamiltonian (8) for g24 = rate for the polariton p0 via spontaneous emission from 0 and ε =0 reads, level 4. Note that U > 0 for ∆ < 0 and vice versa. In a similar way, the two photon detuning ε leads to an energy I 2 2 H = µ p†p +µ p† p +µ p† p , (12) shift of εg B for the polariton p†, which plays the role g24=0,ε=0 0 0 0 + + + − 0 − − − of a chemical potential in the effective Hamiltonian.   where the frequencies are given by µ = 0, µ = (δ Hence, provided (14) holds, the Hamiltonian for the 0 + − A)/2 and µ = (δ + A)/2. dark state polariton p0† can be written as − The polaritons p0 only contain atomic contributions in 2 2 the two metastable states 1 and 2 but not in level 3, which 2 g H = U p† (p ) + ε p†p , (16) shows spontaneous emission. Hence, these polaritons do eff 0 0 B2 0 0 not populate radiating atomic levels and are thus called   in the rotating frame. dark state polaritons [76]. Let us now look at an array of interacting cavities, each As we will discuss now, the dynamics of the dark state under the conditions discussed above. The first term of the polaritons p0 decouples from the dynamics of the remain- Hamiltonian (5) has already been incorporated in the po- ing species p+ and p for a suitable parameter regime. − lariton analysis for one individual cavity. The second term, In this regime the frequency µ0 is sufficiently separated which describes photon tunnelling, transforms into the po- from µ+ and µ , such that neither the interaction with the − lariton picture via (10). To distinguish between the dark atomic level 4 nor the photon tunnelling between cavities state polaritons in different cavities, we introduce the no- can induce any mixing between the polariton species p0 tation † to label the polariton † in the cavity at position and p . pR p0 ± R. The photon hopping translates to polariton hoppings according to, 4.1.1. Polariton-polariton interactions and polariton Ω2 aR† aR′ pR† pR′ (17) tunnelling ≈ B2 + ”terms for other polariton species” .

The polaritons p0† interact with each other via a self Kerr Contributions of different polaritons decouple due to the nonlinearity. To derive this interaction, one first needs to separation of their frequencies µ0, µ+ and µ . As a con- write the full Hamiltonian HI , (8), in the polariton basis, − † N j N j sequence the Hamiltonian for the polaritons pR takes on expressing the operators j=1 σ22 and a† j=1 σ24 in the form (7), with terms of p†, p† and p† . The coupling of the polaritons to 0 + P P 2 − 2ωCΩ the level 4 of the atoms via the dipole moment g24 reads, (18) J = 2 2 α , Ng13 + Ω N j gΩ 2 U as given by equation (15) and µ = ǫg2/B2. g σ a + h.c. g S† p + h.c. , 24  42  ≈− 24 B2 14 0 j=1 X     (13) 1 N j 4.1.2. Polariton lifetime: spontaneous emission and cavity where S† = σ . In deriving (13), a rotating 14 √N j=1 41 decay wave approximation was used: In a frame rotating with re- P spect to (12), the polaritonic creation operators rotate with We now discuss the loss sources in the system. Levels 1 the frequencies µ0, µ+ and µ . Furthermore, the operator and 2 are metastable and therefore have negligible decay − S14† rotates at the frequency 2µ0, (c.f. (8)). Hence, pro- rates on the relevant time scales. The decay mechanism for vided that the dark state polariton then originates mainly from pho- tons that leak out of the cavity and from the small, but non- g24 , ε , ∆ µ+ , µ (14) zero, population of level 4, which leads to spontaneous | | | | | | ≪ | | | −|

Copyright line will be provided by the publisher Laser & Photon. Rev. 1, No. 1 (2008) 7

emission. The resulting effective decay rate for the dark A polariton model which generates an effective Lieb- state polariton reads Liniger Hamiltonian and is capable of reproducingthe Tonks Girardeau regime was proposed in [45]. The scheme em- Ω2 g2 g2Ω2 ploys atoms of the same level structure as in figure 4 that Γ = κ + Θ(n 2) 24 γ , (19) 0 B2 − ∆2B4 4 couple to light modes of a tapered optical fibre. Further- more quantum phase transitions of Dicke-model in cavity where κ is the cavity decay rate, γ4 the spontaneous emis- QED have been analysed [82]. sion rate from level 4, n is the average number of pho- tons, and Θ the Heaviside step function. Assuming that g13 = g24, the maximal achievable rate U/Γ0 can be read- 4.1.4. Two component Bose-Hubbard model ily seen to be g / 4κΘ(n 2)γ . Hence, for large co- 13 − 4 operativity factors, ξ 1, ratios of U/Γ0 1 can be p ≫ ≫ Before we conclude this section, let us briefly review the achieved, which ensures that the dynamics of the Bose- proposal of Ref. [38] for creating a two component polari- Hubbard model can be observed under realistic conditions. tonic Bose-Hubbard model. This model, which describes the dynamics of two independent bosonic particles in each site, reads 4.1.3. The phase transition (j) ′ Heff = µj nk Jj,l jk† lk + h.c. − ′ Of great interest is of course whether the quantum phase k,jX=b,c ;j,l=b,c   transition from a Mott insulator to a superfluid state could + U n(j) n(j) 1 + U n(b)n(c) , (20) be observed with the present polariton approach. In a sys- j k k − b,c k k tem with on average one polariton per cavity, the Mott in- k,jX=b,c   Xk sulator state is characterised by the fact that the local state where bk† (ck† ) create polaritons of the type b(c) in the cav- in each cavity is a single polariton Fock state with van- (b) (c) ishing fluctuations of the polariton number. The superfluid ity at site k, nk = bk† bk and nk = ck† ck. µb and µc are phase in contrast shows polariton number fluctuations. the polariton energies, Ub, Uc and Ub,c their on-site inter- Figure 5 shows numerical simulations of the full dy- actions and Jb,b, Jc,c and Jb,c their tunnelling rates. This namics of an array of three cavities as described by equa- model exhibits interesting quantum many-body phenom- tion (9), including spontaneous emission and cavity decay, ena such as spin density separation [83], spin order in the and studies the number of polaritons in one cavity/site, Mott regime [84] and phase separation [85]. The Hamiltonian (20) can be realised in the same set- n = p†p for site l, and the number fluctuations, F = l l l l up as above, by considering the atom-photon interactions h2 i 2 (p†p ) p†p . A comparison between the full model h l l i − h l li in a dispersive regime with δ Ω,g. Then, the polariton (9) and the effective model (7) is done by considering the operators (10) read, ≫ differences in the occupation numbers, δn = [n ] l l cavities − [nl]BH, and number fluctuations, δFl = [Fl]cavities [Fl]BH, 1 at each time step [21]. Initially there is exactly− one po- p0† = B gS12† Ωa† µ0 =0 − 2 lariton in each cavity and parameters for toroidal micro- 1   B B p† ΩS12† + ga† S13† µ = , 9 B δ − δ cavities from [81], g24 = g13 =2.5 10 s 1, γ4 = γ3 = − ≈ − − 2 7 1 5× 1 −  1  B 1.6 10 s− and κ = 0.4 10 s− , are assumed. As p+† S13† + ΩS12† + ga† µ+ = δ + × × ≈ δ δ the system is driven from a Mott insulator state to a super-   1 fluid state by ramping up the driving laser Ω, the particle to leading order in δ− . Note that now, not only the po- number fluctuations increase significantly. The numerics lariton p0, but also p† do not experience loss from sponta- in figure 5 show good agreement of the full dynamics as − neous emission to leading order in δ 1. One can therefore described by equation (9) and the dynamics of the corre- − define two dark state polaritons species sponding Bose-Hubbard model (7) and thus confirm the possibility of observing the Mott insulator-to-superfluid 1 1 transition in such a system. The oscillatory behaviour is b† = gS† Ωa† ; c† = ΩS† + ga† . B 12 − B 12 related to the fact that the initial state of the system is, due    (21) to the nonzero tunnelling rate J, not its ground state. This From an analysis similar to the one outlined above for the initial state has been chosen since it’s preparation in an ex- one component case, it can be shown that the dynamics of periment is expected to be easier compared to other states these two polaritons species is given by the two component and figure 5 thus shows the dynamics that is expected to Bose-Hubbard model (20), where the parameters are given 2 2 2 2 2 be observed in an experiment. g Ω gΩ g24g Ω by, Jbb = α B2 , Jcc = α B2 , Jbc = α B2 , Ub = B4∆ , The phase diagram for this model and its dependence 2 2 2 2 2 2 2− g24g Ω g24(g Ω ) on the atom number N have been analysed in [39], where Uc = 4 2 and Ubc = 4 − 2 . For fur- − B (∆+2B /δ) − B (∆+B /δ) also a glassy phase has been predicted for a system with ther details of the derivation the reader is referred to Ref. disorder, see section 6 for details. [38].

Copyright line will be provided by the publisher PSfrag replacements

8 Michael J. Hartmann et. al.: Quantum Many-Body Phenomena in Coupled Cavity Arrays

2 a10 b 0.7 106 Ω 0.6

0.5 105 0.4 101

MHz 0.3

0.2 U 0 1 J . 100 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 t in 10−6 s t in 10−6 s c0.02 d 1 0.015

0.01 0.8 0.005 n1 0 0.6

−0.005 F1 0.4 −0.01

−0.015 0.2 −0.02

−0.025 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 −6 PSfrag replacements PSfrag replacements−6 t in 10 s t in 10 s Ω g13 Ω Figure 5 The Mott insulator to superfluid transition for 3 polaritonsg13in 3 cavities compared to 3 particles in a 3 site Bose-Hubbard 1 model. a: Log-plot of U and J and linear plot of the time dependent Ω (inset). b: Number fluctuations for10 polaritons in cavity 1, F1, for 2 a single quantum jump trajectory. c: Difference between number fluctuations for polaritons in a cavity and10 number fluctuations in a pure 3 Bose-Hubbard model, δF1, for a single quantum jump trajectory. d: Expectation value and fluctuations10 for the number of polaritons 9 −3 7 −1 in one cavity/site according to the effective model (7) with damping (19). g24 = g13 = 2.5 10 s10 1, γ4 = γ3 = 1.6 10 s , 5 −1 10 −1 7 −1 × −−2 × κ = 0.4 10 s , N = 1000, ∆ = 2.0 10 s and 2ωC α = 1.1 10 s . The Rabi10 frequency of the driving laser is × 10 −1 × 12 −1 × −1 increased from initially Ω = 7.9 10 s to finally Ω = 1.1 10 s . Deviations from the pure10 Bose-Hubbard model are about × × 0 2%. Reprinted with permission from [21]. 10 1 10 1 10 2 10 3 1 1 An examplehow the parametersof the effectiveHamil- 1010 10 −4 0 −4 0 tonian (20) vary as a function of the intensity of the driv- 10 10 10 10 −3−1 −1 ing laser Ω is given in figure 6, wherethe parametersof the 1010 10 −2−2 −2 atom cavity system are chosen to be g = g , N = 1000, 1010 10 24 13 −3 −1 −3 , √ and . Fig- 1010 10 ∆ = g13/20 δ = 2000 Ng13 α = g13/10 0 Ω − 10 Ω ure 6 shows the interactions U , U and U , the tunnelling 1 2 3 g13 b c bc 1 10 10 10 1 2 3 g13 10 10 10 10 rates Jbb, Jcc and Jbc and µc µb as a functionof Ω/g13. For g Ω one has U | −U |, U and J J ≈ | bc| ≪ | b| | c| bb ≈ cc ≈ Jbc. Whenever µc µb < Jbc , b† polaritons get con- Figure 6 Left: The polariton interactions Ub (dashed line), Uc | − | | | (grey line) and (solid line) in units of as a function verted into c† polaritons and vice versa via the tunnelling Ubc g13 of Ω/g13. Right: The tunnelling rates Jbb (dashed line), Jcc Jbc. With the present choice of α and δ, this happens for | | | | (grey line) and Jbc (solid line) together with µc µb (dot- 0.16g<Ω< 1.6g. | | | − | ted line) in units of g13 as a function of Ω/g13. The parame- ters of the system are g24 = g13, N = 1000, ∆ = g13/20, − δ = 2000√Ng13 and α = g13/10. Reprinted with permission 4.2. Frozen light: photonic Bose-Hubbard from [38]. models

In the last section we have discussed that, as soon as a than unity, ξ 1 is achieved, it is possible to engineer po- strong coupling regime with a cooperativity factor larger laritonic Bose-Hubbard≫ models in coupled cavity arrays.

Copyright line will be provided by the publisher Laser & Photon. Rev. 1, No. 1 (2008) 9

The prospect that even higher cooperativity factors (of the a b 0.998 0.06 order of 102 103 [81]) can be realised in some cav- − 0.05 ity QED setups, opens up the possibility of engineering 0.9975 0.04 a Bose-Hubbard model also for the cavity mode photons. 0.997 0.03 Such a setting is actually conceptually simpler than the po- 0.02 laritonic and has already been outlined before. Given that 0.9965 the photonic tunnelling term is a consequence of the inter- 0.01 0.996 0 action of neighbouring cavities, to realise a photonic Bose- 0 0.5 1.0 1.5 2.0 0 0.5 1.0 1.5 2.0 −6 −6 Hubbard model one needs to find a way to efficientlyPSfrag replacements cre- c t in 10 s d t in 10 s 0.02 1 ate large photon Kerr nonlinearities using the interaction U 0.01 of the atoms with the cavity mode. In this section we re- J 0.8 0 view two proposals to this aim. 0.6 −0.01 To create a pure photonic Bose-Hubbard model would 0.4 −0.02 be very interesting as in this way new regimes for light, 0.2 which do not occur naturally, could be engineered. For ex- −0.03 −0.04 0 ample, in the photonic Mott insulator exactly one photon 0 0.5 1.0 1.5 2.0 0 0.5 1.0 1.5 2.0 −6 −6 exists in each cavity, provided the whole structure contains n1 t in 10 s t in 10 s on average one photon per cavity. Moreover, the photons are localised in the cavity they are in and are not able to Figure 7 Mott insulator state for three photons in three cavi- hop between different cavities. In such a situation pho- ties: a: photon number n1, b: photon number fluctuation F1, c: tons behave as strongly correlated particles that are each differences between full and effective Bose-Hubbard model, δn1 ”frozen” to their lattice site, a system that would corre- (solid line) and δF1 (dashed line) and d: n1 and F1 as given by spond to a crystal formed by light. a solution of a master equation for the three site Bose-Hubbard The first proposal for generating strong Kerr photon model (7) including damping Γ . Reprinted with permission from nonlinearities, which we discuss here, is the EIT based set- [25]. ting analysed in the last section. As shown in [25], in the regime characterised by g Ω and g g ∆Ω , the ≪ 24 ≪ | | dark sate polariton p0 becomes a photon and the nonlin- Figure 8 The level struc- earity, c.f. equation (15), is given by ture for the generation of the effective Hamiltonian (26). g24g g 2 2 U = g24(a†) a . (22) (g, ∆ ) and (Λ, ∆ ) are − ∆Ω Ω 1 2 | | the Rabi frequencies and de- The effective decay rate for the photons subject to the non- tunings of the cavity-atom linearity is given by and laser-atom interactions. Reprinted with permission 2 2 2 2 Γ = κ + Θ(n 2)g g ∆− Ω− γ , (23) from [26]. − 24 4 where cavity decay is the main source of loss since g24g ∆Ω . ≪ | | state. The cavity mode couples dispersively only to the 0 Taking, for example, the case of toroidal micro-cavities, − for which the achievable parameters are predicted to be 2 transition and the levels 0 and 1 are coupled via a far- 9 1 7 1 detuned Raman transition. Finally, the detuning associated g24 = 2.5 10 s− , γ4 = 1.6 10 s− , and κ = 0.4 5 1 × × × with the lasers and with the cavity mode are assumed to 10 s− [81], and setting the parameters such that g/Ω = be sufficiently different from each other. The origin of the g24g/ ∆Ω = 0.1, leads to U/Γ = 625, which would be sufficient| to| observe a Mott insulator of photons, as shown nonlinearity can be understood entirely by considering the in figure 7. Stark shifts due to the dispersive interactions. A second proposal to realise large Kerr nonlinearities The full Hamiltonian of the system, in the interaction is due to Ref. [26]. There a new configuration, involving a picture, reads three levels lambda structure and only one coupling to the i∆1t √ i∆2t cavity mode (see figure 8), was shown to produce nonlin- H = ge− aS20 + 2Λe− S2+ + h.c. (24) earities as large as the EIT based setting discussed before. where S := 2 0 and S := 2 + . As A possible advantage of this approach is that, as the only 20 k | kih k| 2+ k | kih k| shown in figure 8, (g, ∆1) and (Λ, ∆2) are the Rabi fre- three levels are used, it might be easier to implement it in quencies and detuningsP of the cavity-atomP and laser-atom some of the new cavity QED systems, such as quantum interactions, respectively. In the dispersive regime, charac- dots in photonic crystal cavities or Cooper pair boxes in terised by √Ng 1 and √NΛ 1, where furthermore stripline cavities. ∆1 ≪ ∆2 ≪ The relevant atomic level structure, depicted in figure √Ng, √NΛ ∆2 ∆1 , the processes driven by the 8, is a Λ system with two metastable states and an excited cavity-atom and≪ laser-atom | − | interactions can be treated in-

Copyright line will be provided by the publisher 10 Michael J. Hartmann et. al.: Quantum Many-Body Phenomena in Coupled Cavity Arrays

dependently and the excited state will hardly be populated. excitations in level 2 is equal to the initial number of dark Hence an effective Hamiltonian for the two metastable states state polaritons, while the other polariton species will not and the cavity mode can be found by adiabatically elimi- be transferred to the atomic level 2 (c.f. (10)). In this state nating the excited state. Dropping out terms proportional cavity loss is eliminated and the atomic excitations are to the identity, it is given by very long lived. They can then be measured with high ef- ficiency and precision using state-selective resonance flu- g2 Θ orescence as proposed in [86] and [87]. H1 = a†aS00 + (S10 + S01), (25) ∆1 2

2 with Θ := 2Λ /∆2. In a second interaction picture with respect to H0 = 2 2 g int g Θ 4.3.2. Two-component model a†a it reads H1 = a†a (S+ + S +)+ S3, 2∆1 2∆1 − − 2 where S+ := +k k , S + = (S+ )†, and, as − k | ih− | − − before, S := N + + . In the 3 Pk=1 k k k k The number statistics for both polariton species for the basis the system can| thusih be| viewed − |− ih− as an| ensemble{|±i} of P model (20), b† and c†, in one cavity can again be mea- two level atoms driven by a laser with a photon-number- sured using state selective resonance fluorescence. In the dependent Rabi frequency.Considering the dispersive regime 2 two-component case the STIRAP can however not be ap- of this system, √Ng , the atoms prepared in the 2∆1 Θ 1 plied as in the single component case because the energies ≪ |−i2 state will experience a Stark shift proportional to (a†a) , µb and µc are similar and the passage would thus need to which gives rise to a strong Kerr nonlinearity. The effec- be extremely slow to be adiabatic. 4 g tive Hamiltonian is be given by Heff = 2 a†aa†aS3. For two components, one can do the measurements as 4∆1Θ Therefore, if all atoms are prepared in the state, one follows. First the external driving laser Ω is switched off. obtains an essentially absorption free Kerr nonlinearity|−i given Then the roles of atomic levels 1 and 2 are interchanged in by each atom via a Raman transition by applying a π/2-pulse. 2 To this end the transitions 1 3 and 2 3 are driven √Ng √Ng ↔ ↔ Hkerr = ga†aa†a. (26) with two lasers (both have the same Rabi frequency Λ) in 2∆1Θ 2∆1 2 two-photon resonance for a time T = πδΛ/ Λ (δΛ is the detuning from atomic level 3). The configuration| | is shown in figure 9a. This pulse results in the mapping 1j 2j 4.3. Measurements for all atoms j. | i ↔ | i Next another laser, Θ, that drives the transition 1 4 In this section we discuss ways how information about ↔ the state of the polaritonic Bose-Hubbard models can be is switched on, see figure 9b. Together with the coupling obtained in experiments. The possibilities to access the g24, this configuration can be described in terms of three system in measurements are somewhat complementary to polaritons, q0†, q+† and q† , in an analogous way to p0†, p+† − those for cold atoms in optical lattices [8]. Whereas time of and p† , where now the roles of the atomic levels 1 and flight images give access to global quantities such as phase 2 and− the levels 3 and 4 are interchanged. Hence, if one coherence in optical lattices, coupled cavity arrays allow chooses Θ = Ω the π/2-pulse maps the b† onto the dark to access the local particle statistics in measurements. Un- state polaritons of the new configuration, q†, whereas for like optical lattices they would thus allow to obtain direct 0 Θ = Ω it maps the c onto q†. The driving laser is then experimental evidence for the exact integer particle num- † 0 adiabatically− switched off, Θ 0, and the corresponding ber in each lattice site in a Mott Insulator regime. We first → discuss the measurement scheme for the one-component STIRAP process maps the q0† completely onto atomic exci- modeland then turn to its extension to forthe two-component tations of level 1. This process can be fast since the detun- model. ing ∆ is significantly smaller than δ and hence the energies of all polariton species q0†, q+† and q† are well separated. Another π/2-pulse finally maps the− excitations of level 1 4.3.1. One-component model onto excitations of level 2, which can be measured by state selective resonance fluorescence in the same way as for the For the model (7), the number of polaritons in one individ- one-component model. ual cavity can be measured via state selective resonance The whole sequence of π/2-pulse, STIRAP process fluorescence. To that end, the polaritons are transformed and another π/2-pulse can be done much faster than the into purely atomic excitations by adiabatically switching timescale set by the dynamics of the Hamiltonian (20) and 2 d off the drivinglaser. The process is adiabatic if gB− Ω b† or c† can be mapped onto atomic excitations in a time dt ≪ µ+ , µ [76], which means it can be fast enough to pre- in which they are not able to move between sites. The pro- |vent| polaritons| −| from hopping between cavities during the cedure thus allows to measure the instantaneous local par- switching. Hence, in each cavity, the final number of atomic ticle statistics of each species separately.

Copyright line will be provided by the publisher PSfrag replacements PSfrag replacements b a g24 LaserΘ & Photon. Rev. 1, No. 1 (2008) 11

Λ ∆ Λ − − b 4 (E2 2E1 )/g a 3 ∆ 0.6 − ε δΛ δΛ 3 0.5 ωC ε 0.4 ωC Λ g24 0.3 Λ Θ 2 PSfrag replacements 0.2 2 0.1 4 1 ∆/g 1 1 2 3 4

Figure 9 a: Configuration of the π/2-pulse. Two driving lasers Figure 10 Dependence of the effective repulsion E− 2E− on 2 − 1 in two-photon transition with identical Rabi frequencies Λ couple the detuning ∆. The effective repulsion vanishes for large ∆ to the atomic transitions 1 3 and 2 3. b: Configuration for the STIRAP process. A↔ driving laser↔ couples to the 1 4 transition with Rabi frequency Θ. The cavity mode couples↔ to transitions 2 4 and 1 3, where the coupling to 1 3 is which plays the role of an effective on-site repulsion. This ineffective and↔ not shown.↔ Reprinted with permission from↔ [38]. repulsion can be tuned via the detuning ∆ as figure 10 shows. The full effective many-body model in a cavity array that is based on the effective interaction (29) reads, 5. Models Based on the Jaynes-Cummings JC Hamiltonian H = HR J aR† aR′ + h.c. , (30) − R Many-body dynamics on a lattice that can either be domi- X X   nated by the tunnelling between adjacent lattice sites or by where R labels the site of a cavity. Each cavity contains an on-site interaction can also be generated in an array of one atom interacting via the Jaynes-Cummings interaction coupled cavities with Jaynes-Cummings type coupling be- with the cavity mode and photons tunnel between neigh- tween cavity photons and two level atoms in each cavity. bouring cavities at a rate J. As shown in [22] and [23], the anharmonicityof the energy Figure 11 gives evidence for the Mott insulator and su- levels of the dressed states of a Jaynes-Cummings Hamil- perfluid regimes of Hamiltonian (30). The probabilities to tonian acts as an effective on-site interaction. To see this let find one or two excitations in one cavity are shown for an us recall the Jaynes-Cummings Hamiltonian for one atom, array where initially one cavity had no and all others one excitation. For detuning , double occupations never JC ∆ =0 H = ω a†a + ω e e + g(a† g e + a e g ), (27) C 0| ih | | ih | | ih | occur, indicating a Mott insulator behaviour. For detunings ∆ g, on the other hand, occupation numbers larger than where ωC ,ω0 are the frequencies of the resonant mode of ≫ the cavity and of the atomic transition, respectively, g is one occur as to be expected for a superfluid regime. In the Jaynes-Cummings coupling between the cavity mode and deep Mott insulator regime, the model (30) can also repro- duce effective spin models [?], see section 7.1. the two level system, a† is the creationoperatorof a photon in the resonant cavity mode, and g , e are the groundand Effective many-body physics based on the Hamilto- excited states of the two level system.| i | i nian (30) can not only be observed for a single two level system in each cavity, as assumed in (27), but also for Since the energy of a photon ωC and the atomic tran- setups with several two level systems per cavity. Such a sition energy ω0 are much greater than the atom photon coupling g, the number of excitations is conserved for the model can describe photonic crystal micro-cavities doped Hamiltonian (27). Hence it can be diagonalised for each with substitutional donor or acceptor impurities. This ap- manifold with a fixed number of excitations n separately. proach to implementing effective many-bodymodels, which can have suitable parameters, has been proposed in [44]. The energy eigenvalues for n excitations read E0 =0 and The phase transitions of a model with several two level 2 ∆ 2 ∆ systems in each cavity have also been studied in [39,47]. E± = nω + ng + (28) n C 2 ± 4 To allow for the observation of coherenteffective many- r body dynamics, the effective on-site repulsion (29) has to for n 1, where ∆ = ω ω . 0 C be much stronger than the rates of spontaneous emission Let≥ us now consider− an array of cavities where the and cavity decay. Since the dressed states are superposi- atom photoninteraction in each cavity is described by equa- tions of a photonic and an atomic component, they are tion (27) and assume that we have prepared the system in vulnerable to both decay mechanisms. Spontaneous emis- a state with one excitation of energy E− in each cavity. 1 sion can be suppressed by increasing the detuning as Since the lowest energy for two excitations in one cavity is ∆ this causes the occupation of the excited atomic level to E−, moving one additional excitation to a cavity requires 2 decrease. An increase of however also causes the re- an extra energy of ∆ pulsion to decrease as shown in figure 10. The effective 2 2 many-body dynamics is thus only observable if each cav- 2 ∆ 2 ∆ ∆ E2− 2E1− =2 g + 2g + , (29) ity operates in the strong coupling regime, see section 2. − r 4 − r 4 − 2

Copyright line will be provided by the publisher 12 Michael J. Hartmann et. al.: Quantum Many-Body Phenomena in Coupled Cavity Arrays

Evolution for initial state 11...1011...1 1 1.4

3 sites, no decay 0.9 4 sites, no decay 1.2 5 sites, no decay N=1, δ=0 6 sites, no decay 3 sites,κ,γ=0.1A 4 sites,κ,γ=0.1A 0.8 1 5 sites,κ,γ=0.1A 6 sites,κ,γ=0.1A 0.7 7 sites, no decay 0.8

0.6 0.6 Order parameter

0.5 0.4 Mott Superfluid

0.4 0.2

0.3 N=1,δ=100g 0 −2.5 −2.25 −2 −1.75 −1.5 −1.25 −1 −0.75 −0.5 −0.25 0 0.25 0.5 Detuning between atom and cavity (logarithmic scale) 0.2

0.1Probabilty of exciting to the N−th polaritonic level in middle site N=2,δ=100g N=2,δ=0 0 Figure 12 The order parameter Fl as a function of the detun- 0 2 4 6 8 10 12 14 16 Time in 1/g 5 x 10 ing between the hopping photon and the doped two-level sys- tem in logarithmic units of the matter-light coupling g. Simu- Figure 11 Probabilities for finding one (thick solid line) or two lations include results for 37 sites, with and without dissipation due to spontaneous emission and cavity leakage. The graphs are (thick dashed line) excitations in the middle cavity which was obtained by diagonalisation of the Hamiltonian (30). initially taken to be empty (in the resonant regime). In this case, as there were initially fewer excitations than cavities, oscillations occur for the single excitation. The double one is never occupied due to the blockade effect-Mott insulator phase (thick dashed For the characterisation of the various phases, possible line). For the detuned case however, hopping of more than one signatures in the visibility of interference fringes of pho- excitation is allowed (superfluid regime) which is evident as the tons that could be emitted by the cavities have been dis- higher excitation manifolds are being populated now (dashed line cussed [39,50]. The visibility of the interference pattern and thin solid line). Reprinted with permission from [22]. is given in terms of the photonV number distribution in mo- 1 L 2πi(j l)/L mentum space, S(k) = L j,l=1 e − aj†al ,(L is the number of sites in a one dimensional array)h byi = P V As for the polaritonic Bose-Hubbard model (section (max(S) min(S)) / (max(S)+ min(S)). In the super- fluid regime,− approaches unity whereas it is very small 4.1), the transition between Mott insulator and superfluid V phase can be characterised by the fluctuations of the num- for a Mott insulator. ber of excitations in one single cavity. Here, the number of excitations in cavity l is measured by the operator = Nl 6. Phase Diagrams al†al + el el and their fluctuations Fl are given by the | ih | 2 2 ground state variance of l, Fl = l . Figure As already indicated above one of the key questions for N hNl i − hN i 12 shows the number fluctuations Fl, which play the role effective many-body systems in coupled cavity arrays is of an order parameter, as a function of the atom cavity de- which quantum phase transitions these can show. The phase tuning ∆. Whereas the fluctuations are very small in the diagram of both, polaritonic Bose-Hubbard models and Mott insulator regime (small detuning), they become sig- Jaynes-Cummingsinteraction based models, has been anal- nificantly larger in the superfluid phase (large detuning). ysed in various settings. Approaches based on different techniques have been carried out for one- and two-dimensional Various approaches have shown that this model indeed models. exhibits a quantum phase transition that closely resem- The phase boundary between a Mott insulator and a bles the Mott insulator to superfluid transition of the Bose- superfluid phase can be determined in a grand canonical Hubbard model [22,23,39,44,40,42], see also section 6. approach. To this end, a chemical potential µ is introduced On the other hand deviations of the phase transitions in the and the grand canonical Hamiltonian reads, model (30) from the Mott insulator to superfluid transition HGC = H µ , (31) in a Bose-Hubbard model have also been found. − N In this analysis exact diagonalisations of small mod- where the operator counts the total number of excita- els of very few cavities have been performed [41]. Here, tions and is given byN insulator states have been found, where the number of ex- citations in one site does not fluctuate and the effective = aR† aR + eR eR (32) N | ih | particles are either atomic (for ∆/g < 1 and J < ∆ ) R − | | | | X   or polaritonic (∆/g . 1 and J . ∆ ). Superfluid states for Jaynes-Cummings interaction based models (30) and on the other hand show local| excitation| | | number fluctua- tions and can either contain polaritonic (∆/g < 1 and = a† aR + 2R 2R + 3R 3R +2 4R 4R − R J ∆ ) or photonic (∆ < 0, J > ∆ & g or ∆ > N R | ih | | ih | | ih | |0, |J ≈& | g|) effective particles. | | | | X  (33) | |

Copyright line will be provided by the publisher Laser & Photon. Rev. 1, No. 1 (2008) 13

-0.5 for EIT based models (9), which lead to polaritonic Bose- ρ -0.25 N = 1 = 3 N = 2 ρ= 3 -0.75 Hubbard models. Here, an excitation of the atomic level β / ρ= 2 β / 4 is equivalent to two polariton excitations since an atom ρ= 2 ω) -0.5 ω) - -1 needs to absorb two cavity photons to make a transition - µ ( µ

-0.75 ρ ( from level 1 to level 4, whereas the transition from level 1 = 1 -1.25 ρ= 1 to 2 or 3 only requires one cavity photon. The occupations -1 of level 4 thus carry a double weight. -1.5 0 0.08 0.16 β 0.24 0 0.08 0.16 β 0.24 We note that the chemical potential introduced in t/ t/ µ 4 4 equation (31) is conceptually different from the chemical L = 12 L = 12 potential in equation (16). The approach in section 4.1 3 L = 24 3 L = 24 ε β L = 36 β L = 36 / L = 48 / L = 48 considers a closed system without particle exchange with κ L = 64 κ 2 L = 80 2 the environment and ε is an energy shift that affects the L = 96 polaritons. The grand canonical approachin this section, in 1 1 contrast, considers a situation in which particle exchange 0 0 with the surrounding is permitted and is used because of 0 0.08 0.16 β 0.24 0 0.08 0.16 β 0.24 its convenience for determining the phase diagrams. t/ t/ The total number of excitations is conserved for both models, [H, ]=0. The boundary between Mott insula- Figure 13 Upper panels: Phase diagram for the Hamiltonian tor and superfluidN phases is determined by the value of µ model (30), with N = 1, 2 atoms inside each cavity at ǫ = ω. for which adding or removing a particle does not require Lower panels: System compressibility κ for the first lobe (i.e., energy. ρ = 1), for different system sizes L, with N = 1 (left) and N = 2 (right). Reprinted with permission from [39].

6.1. Phase diagrams for one dimensional -0.25 -0.75 models ρ= 2 Ω

/ N = 1 N = 2

-0.5 Ω

/ -1 µ ρ The phase diagrams for the one dimensional versions of -0.75 µ = 2 -1.25 both, polaritonic Bose-Hubbard and Jaynes-Cummings in- -1 ρ= 1 -1.5 ρ= 1 teraction based models, have been obtained by DMRG cal- -1.25 culations (see [88] for a review on DMRG). These results -1.75 -1.5 were presented in [39], where in addition the existence of 0 0.1 0.2 0.3Ω 0.4 0 0.1 0.2 0.3Ω 0.4 a polaritonic glassy phase in setups with disorder has been t/ t/ 3 3 predicted. L = 12 L = 8 Ω Ω

/ L = 24 / L = 16 The authors consider a linear array of L cavities with L = 36 L = 24 κ κ various numbers of atoms, N, in each cavity. Figure 13 2 L = 48 2 L = 32 shows the phase diagram and compressibility for the model (30). Here ω denotes the frequency of the cavity-photons 1 1 (ω ωC), ǫ the energy of the excited atomic levels (ǫ 0 0 ≡ ≡ 0 0.1 0.2 0.3 0.4 0 0.1 0.2 0.3 0.4 ω0), β the atom-photon coupling (β g) and ρ the number t/Ω t/Ω of excitations per cavity. The compressibility≡ κ is defined by κ = ∂ρ/∂µ. The Mott insulator is characterised by a vanishing compressibility whereas κ is finite for a super- Figure 14 Upper panels: Phase diagram of model (9). The de- fluid. Figure 14 shows the phase diagram and compress- tuning parameters are set to zero and g/Ω = 1. Lower panels: ibility for the model (9). System compressibility κ for the first lobe, for different system For a Bose-Hubbard Hamiltonian, the lobe width at sizes L, with N = 1 (left) and N = 2 (right). Reprinted with permission from [39]. t = 0 is given by the on-site repulsion U and the criti- cal value t∗ of the hopping strength at which the various lobes shrink into a point is proportionalto U. In contrast to the Bose-Hubbard model, in both models, (30) and (9), for fixed N, the critical values t∗ are not proportional to the a correlated hopping of the polaritons, where the hopping lobe width at t = 0. Furthermore, the ratio between the depends on the occupation of the cavity. For larger occu- upper and the lower slopes of the lobes at small hopping is pations, the polaritons have a stronger tendency to tunnel greater than the one predicted for the Bose-Hubbard model between cavities and consequently the superfluid regime is (2npolariton for the lower and 2(npolariton + 1) for the up- already reached for lower values of the hopping term pref- − per slope [89]); this discrepancy disappears on increasing actor t. Deviations of the critical hopping t∗ from the crit- the number of atoms inside the cavity. In terms of an ef- ical hopping in a Bose-Hubbard model vanish for N & 10 fective Bose-Hubbard model, this may be interpreted as atoms per cavity [39].

Copyright line will be provided by the publisher 14 Michael J. Hartmann et. al.: Quantum Many-Body Phenomena in Coupled Cavity Arrays

For the detection of the Mott insulator-superfluidquan- tum phase transition, the polariton and photon fluctuations in a block of linear dimension M (in units of the lattice constant of the array) have been studied [49]. Whereas po- lariton fluctuations are independent of the block size M in the Mott insulator regime, they grow logarithmically with M in the superfluid phase. The prefactor of the logarithm is thereby related to the compressibility of the system. For the photon fluctuations, the scaling of the leading terms is always linear in M and non-critical while the critical be- haviour is encoded in the subleading scaling [49].

6.2. Phase diagrams for two dimensional models

The phase diagrams for the two dimensional version of the Jaynes-Cummings based model (30) have been calculated by various approaches, among those a mean field approach [23] and a variational cluster approach [40].

6.2.1. Mean-field phase diagram for the 2d Jaynes-Cummings based model Figure 15 Slices showing the superfluid order parameter, Ψ, as a function of the photon hopping frequency, here denoted by κ (κ J), and the chemical potential, µ. For detuning, ∆ = 0 (a), A mean-field phase diagram for the two dimensional ver- ∆ =≡ 2β (β denotes the Jaynes-Cummings coupling, β g)(b) sion of the model (30) on a triangular lattice has been ≡ and ∆ = 2β (c) (Note that the figure defines ∆ = ωC ω0 obtained in [23]. Introducing the superfluid order param- and denotes− the cavity resonance by ω, contrary to the text).− Mott eter Ψ = aR and using the decoupling approximation h i insulator lobes are indicated by the regions of Ψ = 0, the low- aR† aR′ = aR† aR′ + aR† aR′ + aR† aR′ , one obtains est of which are shown. Dominating the left-hand edge (where the mean-fieldh i Hamiltonian,h i h ih i photonic repulsion dominates over hopping) is the Mott insu- lator phase (denoted MI), and the superfluid phase is found on MF JC 2 the right-hand edge (denoted SF). The white contour in a corre- H = HR zJΨ aR† + aR + zJ Ψ R − | | sponds to the region where Ψ becomes non-zero, delineating the X h   quantum phase transition. Reprinted with permission from [23]. µ aR† aR + eR eR , (34) − | ih |  i where z counts the number of nearest neighbours and is here set to z =3 since a triangular lattice is assumed [23]. Here again, a chemical potential µ has been introduced in (34) and the phase boundary between superfluid and Mott insulator is given by the value of µ, for which adding or removing one excitation to or from the system does not cost energy. The mean-field phase diagram for the two dimensional model (34) on a triangular lattice is shown in figure 15. Here, Ψ = 0 corresponds to stable Mott lobes, where the number of excitations in each cavity increases with µ. The regions on the right correspond to Ψ = 0, where the sys- tem will be found in a superfluid phase,6 i.e. within the Figure 16 Plateaus with constant density, ρ, indicating regions mean field approximation, the ground state corresponds to with a definite number state excitations, as a function of the a coherent state of photons in each cavity. The number of chemical potential µ and photon hopping frequency κ for detun- photons in each Mott lobe can be confirmed by consider- ing ∆ = 0. Regions with varying ρ have coherent states as the ing the average number of excitations per site in the grand ground state configuration. Reprinted with permission from [23]. canonical ensemble ρ = ∂Eg/∂µ, where Eg is the ground state energy. Figure 16 shows ρ as a function of the hop- ping and chemical potential. The plateaus indicate regions

Copyright line will be provided by the publisher Laser & Photon. Rev. 1, No. 1 (2008) 15

∆ = 0 ∆ = 0 T/g = 0.00 (L = 6) 0.6 MI (n = 2) (a) 1D, T = 0, (a) 1D, c 0.5 T/g = 0.01 (L = 2) L = 2 c 0.5 c T/g = 0.01 (L = 4) 0.4 c Lc = 4 T/g = 0.02 (Lc = 2) 0.4 L = 6 g SF c g T/g = 0.03 (L = 2) L = 8 0.3 c

µ / c µ / 0.3 III MI (n = 1) 0.2 0.2 II I 0.1 0.1 SF 0 0 0 0.04 0.08 0.12 0.16 0.2 0 0.04 0.08 0.12 0.16 0.2 L = 1 x 1 T/g = 0 (b) 2D, T = 0, ∆ = 0 c (b) 2D, ∆ = 0 0.5 Lc = 2 x 1 0.5 T/g = 0.005 T/g = 0.01 Lc = 2 x 2 (L = 1 x 1) SF c T/g = 0.02 0.4 L = 8 0.4 c T/g = 0.03 g g 0.3 III 0.3 µ / µ / II MI (n = 1) 0.2 I 0.2

0.1 0.1 SF 0 0 0 0.01 0.02 0.03 0.04 0.05 0.06 0 0.01 0.02 0.03 0.04 0.05 0.06 t / g (c) T = 0, ∆ = 0 0.5 1D, Lc = 8 (z = 2) 2D, L = 6 (z = 3) 0.4 c 2D, L = 8 (z = 4) Figure 18 Mott-like regions for different for the Hamiltonian g c T 0.3

µ / (30). Mott lobes shrink with increasing temperature T . Reprinted SF 0.2 MI (n = 1) with permission from [40]. 0.1 SF 0 0 0.04 0.08 0.12 0.16 0.2 t / g Phase diagrams for the model (30) in two dimensions Figure 17 Ground-state phase diagram for the Hamiltonian (30): have also been obtained by the stochastic series expansion (a) 1D, (b) 2D and (c) different coordination numbers z. Lines are quantum Monte Carlo method [43], where it has been pre- guides to the eye. Reprinted with permission from [40]. dicted that the universality class of coupled cavity arrays may be different from the Bose-Hubbard model. of constant density and incompressibility which are char- acteristic for a Mott insulator phase. 6.3. Bose-glass phases A mean-field phase diagram for settings with more than A glassy phase for polaritons in a one dimensional cavity one two level system in each cavity has been calculated in array with disorder has been predicted in [39]. This phase [44], where 8 two level systems par cavity have been con- is characterised by a finite compressibility, gapless exci- sidered. Such a model describes photonic crystal micro- cavities doped with substitutional donor or acceptor impu- tation spectrum and zero superfluid density and has been predicted to exist in a Bose-Hubbard model for relative rities. fluctuations of the interaction, δUeff / Ueff , uniformly dis- tributed in an interval of length 2ǫ =h 0.5i in a system of L = 200 sites at filling ρ =1.01 for 0.078 . t / U . 6.2.2. Phase diagrams for 2d models from a variational eff h eff i cluster approach 0.133. Figure 19 shows relative fluctuations of the interac- tion δUeff / Ueff as a function of atom number fluctua- The phase diagram of model (30) for the one and two di- tions over anh intervali δN in the cavities. For the model mensional case and the stability of its Mott insulator states (30), a glassy phase should be visible for , e.g. N = 100 with respect to non-zero temperatures have been analysed particles per cavity and δN 19 (such to haveh ǫ i=0.25) 3 ≃ 2 in [40]. The authors employ a variational cluster approach in the range 7.51 10− . t/β . 1.401 10− . that is expected to become exact with increasing cluster The observability× of quantum phase× transitions from size Lc. Figure 17 shows the phase diagram in one dimen- a superfluid to a Bose-glass to a Mott insulator that is sion, two dimensions and for different coordination num- driven by an decrease of the photon tunnelling rate J has bers z. µ is the chemical potential, g the Jaynes-Cummings also been predicted for a two dimensional model in [44]. coupling and t the tunnellingrate of photonsbetween neigh- Here, photonic crystal micro-cavities doped with substi- bouring cavities. Figure 18 shows the unit filling Mott lobes tutional donor or acceptor impurities have been consid- for both, the one and two dimensional model (30) at finite ered, where fluctuations in the number of impurities (∆N), temperatures T . the impurity-photon coupling (∆gk) and the detunings be-

Copyright line will be provided by the publisher 16 Michael J. Hartmann et. al.: Quantum Many-Body Phenomena in Coupled Cavity Arrays

2 Outside the isovalue surface, only superfluid and Bose glass 1.4 = 1 > = 10 phases exist. eff = 100 U 1.5 1.2 = 1000 > ) / < i ( eff

eff 1 U 7. Effective Spin Models in Coupled Cavity 0.8 1 0 25 50 75 100 /

U In the previous sections we have seen that under appro- δ 0.5 priate conditions interesting bosonic many-body Hamilto- nians can be created and probed in coupled cavity arrays. 0 -1 0 1 2 3 4 There the atoms were used for detection and most impor- 10 10 10 10 10 10 δN tantly to generate an interaction between photons in the same cavity, either actively, as in the polaritonic case, or in an indirect manner, as for the creation of large Kerr non- Figure 19 Relative standard deviation of the effective on-site 4 linearities. In this section we focus on a complementary interaction strength Ueff averaged over L = 10 cavities, as a regime, in which the tunnelling of photons is used to me- function of the atom number fluctuations δN; from left to right diate interactions for atoms in neighbouring cavities. We N = 1, 10, 100, 1000. Dashed line indicates an effective review proposals to realise effective spin lattice Hamilto- variationh i of the interaction strength equal to the standard devia- nians with individual atoms in micro-cavities that are cou- tion of a random interaction uniformly distributed in the interval pled to each other [27,29,31]. Ueff (i)/ Ueff [ ǫ, +ǫ] with ǫ = 0.25. In the inset an exam- ple of theh variationi ∈ − on the on-site effective interaction is shown; Interacting spin systems have a key role in quantum N = 100, δN = 20. Reprinted with permission from [39]. statistical mechanics, condensed-matter physics, and more h i recently in quantum information science. On the one hand they have a simple microscopic description, which allows the use of techniques from statistical mechanics in their analysis. On the other, they have a very rich structure [90] and can hence be employed to model and learn about phe- nomena of real condensed-matter systems. In quantum in- formation science, spin chains can be harnessed to prop- agate and manipulate quantum information, as well as to implement quantum computation. We first discuss how spin models that conserve the total number of excitations can be generated as limiting cases of Bose-Hubbard type mod- els. Then we review an approach where the full anisotropic Heisenberg Hamiltonian can be engineered for atoms that interact via the exchange of virtual photons. Finally we mention schemes that investigate magnetic properties for photons of different polarisations. Figure 20 The isovalue surface that separates the existence and nonexistence of an equilibrium polaritonic quantum phase transi- 7.1. Spin models that conserve the total number tion from the Bose glass to Mott insulator states. Reprinted with permission from [44]. of excitations A certain class of effective spin Hamiltonians can be gen- erated as a limiting case of the effective many-bodyHamil- tonians discussed in sections 4 and 5. Whenever the engi- neered on-site interactions are repulsive and strong enough tween impurity and photon resonances (∆ωph) are expected in experimental realisations. The information whether an to suppress the occupation of states with more than one equilibrium polaritonic quantum phase transition from a excitation almost completely, the states with zero and one excitation can be identified with and respectively Bose glass to Mott insulator states exists is obtained by | ↑i | ↓i plotting the isosurface in which the critical polariton tun- [22] and the dynamics can effectively be described by the nelling rate at the Bose glass Mott insulator phase bound- Hamiltonian z x x y y ary is set to be 10 times larger than the polariton loss rate, H = Bσ + J(σ σ ′ + σ σ ′ ), (35) 6 j j j j j assuming a cavity Q-factor 10 , an impurity number N = j 3 and an impurity photon detuning ∆ = 12g. Inside the X X x/y/z isovalue surface, quantum phase transitions from a super- where the second sum runs over nearest neighbours. σj fluid to a Bose glass to Mott insulator states can be ob- are the Pauli matrices at site j, J a coupling coefficient and served by continuously decreasing the tunnelling rate J. B a magnetic field.

Copyright line will be provided by the publisher Laser & Photon. Rev. 1, No. 1 (2008) 17

e The Hamiltonian (35) conserves the numberPSfrag of spins replacements δb in the state and can thus be solved efficiently. Much ∆a richer spin Hamiltonians| ↑i can be generated when putting atoms of a Λ level structure into the cavities. δa ∆b gb

7.2. The anisotropic Heisenberg model Ω ωb ωa b Ω a b A particular interesting spin Hamiltonian is the anisotropic ga ω Heisenberg or XYZ model, given by ab

z a H = Bσj (36) j X y y Figure 21 Level structure, driving lasers and relevant cou- x x z z x x y y + Jxσj σj′ + Jyσj σj′ + Jzσj σj′ . plings to the cavity mode to generate effective σ σ - and σ σ - ′ couplings for one atom. The cavity mode couples with strengths ga and gb to transitions a e and b e respectively. | i ↔ | i | i ↔ | i This model contains a rich phase diagram and can be used One laser with frequency ωa couples to transition a e | i ↔ | i to study quantummagnetism as a limiting case ofthe fermionic with Rabi frequency Ωa and another laser with frequency ωb to Hubbard model, which is believed to contains the main b e with Ωb. The dominant 2-photon processes are indi- cated| i ↔ in | fainti grey arrows. Reprinted with permission from [27]. features of high Tc superconductors [90]. From a quan- tum information perspective, such a model can be used to create cluster states, highly entangled multipartite quan- tum states which are universal for quantum computation by single-site measurements [91]. 7.2.1. XX and YY interactions The Hamiltonian (36) can be generated in an array of coupled cavities where each cavity contains one atom of a XX and YY interactions can be generated in an array of Λ level structure, c.f. figure 21 and 22. The two spin polar- coupled cavities with one 3-level atom in each cavity (see figure 21). The two long lived levels, a and b , represent isations and are represented by the two longlived | i | i atomic levels| ↑i of the| ↓iΛ level structures. Together with exter- the two spin states. The cavity mode couples to the tran- sitions a e and b e , where e is the excited nal lasers, the cavity mode that couples to the atom inside | i ↔ | i | i ↔ | i | i each cavity can induce Raman transitions between these state of the atom. Furthermore, two driving lasers couple to the transitions a e respectively b e . For the two longlived levels. Due to a large enough detuning be- | i ↔ | i | i ↔ | i tween laser and cavity mode, these transitions can only sake of simplicity we consider here a one-dimensional ar- create virtual photons in the cavity mode which mediate ray, but the generalisation to higher dimensions is straight an interaction with another atom in a neighbouring cav- forward. The Hamiltonian of the atoms reads ity. With appropriately chosen detunings, both the excited N atomic levels and photon states have vanishing occupation H = ω e e + ω b b , (37) and can be eliminated from the description. As a result, the A e| j ih j | ab| j ih j | j=1 dynamics is confined to only two states per atom, the long- X lived levels, and can be described by a spin 1/2 Hamilto- where the index j counts the cavities, ωe is the energy of nian. The small occupation of photon states and excited the excited level and ωab the energy of level b . The en- atomic levels also strongly suppresses spontaneous emis- ergy of level a is set to zero. The Hamiltonian| i that de- sion and cavity decay. As discussed later, this proposal can scribes the photons| i in the cavity modes is given by equa- be realised as soon as the strong coupling regime of cavity tion. (5). For simplicity of notation we set JC =2ωCα and QED is achieved. for convenience we assume periodic boundary conditions, We note that similar ideas have been recently proposed where HC can be diagonalised via the Fourier transform to realise spin systems of larger spins [31] and with inter- N actions that depend on the direction of the respective lattice 1 2πl N N edge [29]. a = eikj a , k = , l k √ j N − 2 ≤ ≤ 2 N j=1 In the next subsections we first discuss how effective X σxσx, σyσy and σzσz interactions together with an ef- fective magnetic field Bσz can be engineered and then for N odd to give HC = k ωkak† ak with ωk = ωC + explain how these can be combined to generate the full 2JC cos(k). Finally the interaction between the atoms and anisotropic Heisenberg model (36). the photons as well as theP driving by the lasers are de-

Copyright line will be provided by the publisher PSfrag replacements 18 Michael J. Hartmann et. al.: Quantum Many-Body Phenomena in Coupled Cavity Arrays

scribed by e ∆˜b − N ∆˜a δb Ωa iωat ∆b H = e− + g a e a (38) AC 2 a j | j ih j | j=1 δa X   ∆a

Ωb iωbt + e− + gbaj ej bj + h.c. . 2 | ih | Λb    gb ω ν Ωb Here ga and gb are the couplings of the respective transi- tions to the cavity mode, Ωa is the Rabi frequency of one laser with frequency ω and Ω the Rabi frequency of a Ω g a b ω a a b second laser with frequency ωb [92]. The complete Hamil- Λ a ωab tonian is then given by H = HA + HC + HAC . By switching to an interaction picture with respect to ωb H = H + H δ N b b , where δ = ω a 0 A C − 1 j=1 | j ih j | 1 ab − (ω ω )/2, the excited atom levels e and the photons a b j Figure 22 Level structure, driving lasers and relevant couplings can− be adiabatically eliminatedP [93].| Toi obtain effective to the cavity mode to generate effective σzσz-couplings for one spin-spin interactions, terms up to 2nd order need to be atom. The cavity mode couples with strengths ga and gb to tran- considered while dropping fast oscillating terms. For this sitions a e and b e respectively. Two lasers with fre- approach the detunings ∆a ωe ωa, ∆b ωe ωb quencies| i ↔and | i couple| i ↔ with | i Rabi frequencies respectively k ≡ k− ≡ − − ω ν Ωa (ωab δ1), δ ωe ωk and δ ωe ωk (ωab − a ≡ − b ≡ − − − Λa to transition a e and Ωb respectively Λb to b e . δ1) have to be large compared to the couplings Ωa, Ωb,ga The dominant 2-photon| i ↔ | processesi are indicated in faint| i grey ↔ |ar-i and gb. Furthermore, the parameters must be such that the rows. Reprinted with permission from [27]. dominant Raman transitions between levels a and b are those that involve one laser photon and one cavity photon each (c.f. figure 21). To avoid excitations of real photons via these transitions one requires ∆ δk , ∆ δk a − b b − a ≫ Ωagb Ωbga For Ωa = (∆aga/∆bgb)Ωb with Ωa and Ωb real, , . ± x x y y 2∆a 2∆b the interaction is either purely σ σ (+) or purely σ σ Hence whenever the atom emits or absorbs a virtual ( ) and the Hamiltonian (40) becomes the Ising model in photon into or from the cavity mode, it does a transition − a transverse field, whereas the isotropic XY model (Jx = from level a to b or vice versa. If one atom emits a vir- | i | i Jy, i.e. J2 =0) is obtained for either Ωa 0 or Ωb 0. tual photon in such a process that is absorbed by a neigh- The effective magnetic field B in turn can,→ independently→ bouring atom, which then also does a transition between of Jx and Jy, be tunedto assume any value between B a to b , an effective spin-spin interaction has happened. | | ≫ | i | i Jx , Jy and B Jx , Jy by varying δ1. Hence the Dropping irrelevant constants, the resulting effective Hamil- system| | | can| be| driven| ≪ through | | | quantum| phase transition by tonian reads varying the Rabi frequencies and detunings of the driving N lasers. Now we proceed to discuss effective ZZ interac- z + tions. Hxy = Bσj +(J1σj σj−+1 +J2σj−σj−+1 +h.c.), (39) j=1 X where σz = b b a a , σ+ = b a and B, J j | j ih j | − | jih j | j | jih j | 1 and J2 are functions of the various Rabi frequencies and 2 ⋆ detunings . If J2 = J2, this Hamiltonian reduces to the 7.2.2. ZZ interactions XY model,

N z x x y y Hxy = Bσj + Jxσj σj+1 + Jyσj σj+1 , (40) j=1 z z X Effective σ σ interactionscan be generated with the same atomic level configuration as XX and YY terms, but now with J = (J + J )/2 and J = (J J )/2. x 1 2 y 1 − 2 only one laser with frequency ω mediates atom-atom cou- 2 δ1 pling via virtual photons. A second laser with frequency ν To second order they are given by B = 2 2 2 2 − 1 |Ωb| |Ωb| |Ωa| 2 2 is used to tune the effective magnetic field via a Stark shift. 2 2 ∆b 4∆ 4(∆ −∆ ) γbgb γ1ga 4∆b b a b The atoms together with their couplings to cavity mode h 4 “ − − − 2 −2 2 2 2 ga γ2 |Ωa| gb |Ωb| ga , 2 2 and and lasers are shown in figure 22. Again, we consider the +γ1 ∆ (a b) J1 = 4 ∆ + ∆ b − ↔ a b ” ⋆ i “ ” one-dimensional case as an example. The generalisation γ2 ΩaΩbgagb 1 1 J2 = , where γa,b = , γ1 = 2 ∆a∆b N k ωa,b−ωk to higher dimensions is straightforward. The Hamiltonians 1 1 1 Pexp(ik) and γ2 = . H of the atoms and H of the cavity modes thus have the N k (ωa+ωb)/2−ωk N k (ωa+ωb)/2−ωk A C P P Copyright line will be provided by the publisher PSfrag replacements PSfrag replacements

b a

5 5

15 15

25 25

Laser & Photon. Rev. 1, No. 1 (2008) 19 0

0.2 a 1 b 1 same form as above, whereas HAC now reads ) ) 1 1 0.8 ↓ ↓ ( N ( 0.8 Ωa iωt Λa iνt 0.6 , p , p ) HAC = e− + e− + gaaj ej aj ) 1 2 2 1 0.4

| ih | a a ( j=1   ( 0.6 p X p (41) 0.2 0.4 0 Ωb iωt Λb iνt 0 10 20 30 40 50 60 0 10 20 30 40 50 60 + e− + e− + gbaj ej bj + h.c. . −6 −6 2 2 | ih | t in 10 s t in 10 s    Figure 23 The occupation probability p(a1) of state a1 (solid Here, Ωa and Ωb are the Rabi frequencies of the driving | i line) and the probability p( 1) of spin 1 to point down (dashed laser with frequency ω on transitions a e and b ↓ 6 | i → | i | i → line) for the parameters ωe = 10 GHz, ωab = 30GHz, ∆a = e , whereas Λa and Λb are the Rabi frequencies of the ˜ driving| i laser with frequency ν on transitions a e 30GHz, ∆b = 60GHz, ωC = ωe ∆b + 2GHz, ∆a = 15GHz, GHz, − GHz, GHz, and b e . | i → | i Ωa = Ωb = 2 Λa = Λb = 0.71 ga = gb = 1 J = 0.2GHz and δ = 0.0165GHz (plot a) respectively δ = In| i an → interaction | i picture with respect to H = H + C 1 1 0 A 0.0168GHz (plot b). Both,− the occupation of the excited atomic HC , the excited atom levels ej and the photons can be − † states ej ej and the photon number a a are always smaller adiabatically eliminated [93].| Again,i the detunings ∆ h| ih |i h i a ≡ than 0.03. Reprinted with permission from [27]. ωe ω, ∆b ωe ω ωab, ∆˜a ωe ν, ∆˜b ω −ν ω , δ≡k ω − ω −and δk ω ≡ ω −ω have to≡ e − − ab a ≡ e − k b ≡ e − k − ab be large compared to the couplings Ωa, Ωb,Λa,Λb,ga and gb, whereas now Raman transitions between levels a and 7.2.3. The complete effective model b should be suppressed. Hence parameters must be such that the dominant2-photonprocesses are those that involve Making use of the Suzuki-Trotter formula, the two Hamil- one laser photon and one cavity photon each but where the tonians (40) and (42) can now be combined to one effec- atom does no transition between levels a and b (c.f. fig- tive Hamiltonian. To this end, the lasers that generate the ure 22). To avoid excitations of real photons in these pro- Hamiltonian (40) are turned on for a short time interval k k Ωaga Ωbgb cesses, one requires ∆a δ , ∆b δ , . dt ( Hxy dt 1) followed by another time interval dt − a − b ≫ 2∆a 2∆b || || · ≪ Whenever two atoms exchange a virtual photon in this ( Hzz dt 1) with the lasers that generate the Hamil- || || · ≪ scheme, both experience a Stark shift that depends on the tonian (42) turned on. This sequence is repeated until the state of the partner atom. This conditional Stark shift plays total time range to be simulated is covered. the role of an effective σzσz-interaction. Dropping irrele- The effective Hamiltonian simulated by this procedure vant constants, the effective Hamiltonian reads is Hspin = Hxy + Hzz, which is precisely the Heisen- berg anisotropic model of equation. (36) with Btot = B + N ˜ ˜ z z z B. The time interval dt should thereby be chosen such Hzz = Bσj + Jzσj σj+1 . (42) 1 1 1 1 1 1 that Ω− ,g− dt1, dt2 Jx− ,Jy− ,Jz− ,B− and j=1 1 ≪ ≪ X   B˜− , so that the Trotter sequence concatenates the effec- where B˜ and Jz are determined by the detunings and Rabi tive Hamiltonians HXY and HZZ . The procedure can be frequencies of the employed lasers3. Here again, the in- generalised to higher order Trotter formulae or by turning on the sets of lasers for time intervals of different length. teraction Jz and the field B˜ can be tuned independently, The validity of all above approximations is shown in either by varying Ωa and Ωb for Jz or by varying Λa and 2 2 figure 23, where numerical simulations of the dynamics Λb for B˜. In particular, Λa and Λb can for all values of | | | | generated by the full Hamiltonian H are compared it to Ω and Ω be chosen such that either J B˜ or J B˜. a b z z the dynamics generated by the effective model (40). The Both, XX and YY interactions and ZZ≪ interactions≫ can present example considers two atoms in two cavities, ini- be combined to engineer the full Hamiltonian (36). This tially in the state 1 ( a + b ) a , and calculates the can be done via the Suzuki-Trotter formula. √2 | 1i | 1i ⊗| 2i 2 occupation probability p(a1) of the state a1 which cor- 3 ˜ 1 |Λb| | i To second order they are give by B = 2 ˜2 responds to the probability of spin 1 to point down, p( 1). 16∆b −2 h ↓ 2 2 |Ω | Figure 23a shows and for an effective Hamil- ˜ |Λa| |Λb| j 2 p(a1) p( 1) 4∆b ˜ ˜ ˜ j=a,b ˜ + 4˜γjbgj + ↓ − ∆a−∆b − ∆b − ∆j −∆b tonian (36) with Btot = 0.135MHz, J = 0.065MHz, “ 2 2 2 “ 2 ”” x |Ωb| |Ωa| |ΩbP| |Λj | 2 2 J =0.007MHz and J =0.004MHz and hence Btot > 16∆ 4∆b ∆ −∆ ∆ j=a,b ˜ + 4γjbgj y z b − a b − b − ∆j −∆b | | “ 4 P “ ⋆ ⋆ ”2 Jx , whereas figure 23b shows p(a1) and p( 1) for an ef- 2 gb Ωb gb Ωaga +4γbb (a b) and Jz = γ2 | | ↓ ∆b − ↔ 4∆b − 4∆a fective Hamiltonian (36) with the same Jx, Jy and Jz but ” 1 i 1 ˛ 1 exp(ik˛) with γ1 = , γ2 = ˛ ˛ , Btot = 0.025MHz and hence Btot < Jx [71]. N k ω−ωk ˛ N k ω−ωk˛ Discrepancies− between numerical| | results| | for the full 1 P 1 γab 1 P 1 γaa = γbb = N k ω−ω , = N k ω±ω −ω , and the effective model are due to higher order terms for k γba ff ab k P P the parameters , ˜, , and , which lead to relative γ˜ab 1 1 1 1 B B Jx Jy Jz = N k ν±ω −ω , γ˜aa =γ ˜bb = N k ν−ω . γ˜ba ff ab k k corrections of up to 10% in the considered cases. Despite P P

Copyright line will be provided by the publisher PSfrag replacements PSfrag replacements a b

Ps

EvN/ ln 2

20 Michael J. Hartmann et. al.: Quantum Many-Body Phenomena in Coupled Cavity Arrays

0 0.2 a b this lack of accuracy of the second order approximations, 0.4 1 1 0.6 the effective model is indeed a spin-1/2 Hamiltonian as 0.8 0.999 occupations of excited atomic and photon states are negli- 0.8

ln 2 0.6 0.998 s / gible. P

0.996vN 0.4 0.997 0.997E 0.2 0.996 0.998 0 7.2.4. Experimental requirements 0.999 0 10 20 30 40 50 60 0 10 20 30 40 50 60 −6 −6 t in 10 s t in 10 s For an experimental implementation, the parameters of the ˜ Figure 24 a: The von Neumann entropy EvN of the reduced den- effective Hamiltonian, Jx, Jy, Jz, B and B have to be sity matrix of 1 effective spin in multiples of ln 2 and b: the pu- much larger than rates for decay mechanisms via the pho- rity of the reduced state of the effective spin model for 3 cavities tons or the excited states e . With the definitions Ω = | ji where Jz = 0.021MHz. The plots assume that no spontaneous max(Ωa, Ωb), g = max(ga,gb), ∆ = min(∆a, ∆b), the emission took place. For a spontaneous emission rate of 0.1MHz occupation of the excited levels ej and the photon num- (g = 1GHz), the probability for a decay event in the total time | i 2 ber np can be estimated to be ej ej Ω/2∆ and range is 1.5%. Hence, cluster state generation fails with proba- 2 h| ih |i ≈ | | np (Ωg/2∆)γ1 , whereas the couplings Jx, Jy and Jz bility 0.005 n for n 1/0.005 = 200 cavities, irrespectively ≈ | | 2 × ≪ are approximately (Ωg/2∆) γ2. of the lattice dimension. Reprinted with permission from [27]. Spontaneous emission| from| the levels e at a rate Γ | j i E and cavity decay of photons at a rate ΓC thus lead to ef- fective decay rates z 1 x 1 where Sj = (a† aj a† aj ), Sj = (a† aj +a† aj ), 2 j ↑ j ↓ 2 j ↓ j ↑ 2 2 y 1 ↑ − ↓ ↑ ↓ Γ1 = Ω/2∆ ΓE and Γ2 = (Ωg/2∆)γ1 ΓC . † † and the indices and label | | | | Sj = 2i (aj aj aj aj ) j j′ the cavities.↑ ↓ − ↓ ↑ The effective dynamics of Hamiltonian (36) can thus be observed if

Γ (Ωg/2∆) 2γ and Γ (Ωg/2∆) 2γ , 7.4. Cluster state generation 1 ≪ | | 2 2 ≪ | | 2 2 2 which implies ΓE JC g /δ and ΓC JC (JC < Effective spin models in arrays of coupled cavities can be δ/2), where, δ = (ω≪+ ω )/2 ω for the≪ XX and YY used for the generation of cluster states [91], which form, | a b − C | interactions and δ = ω ωC for the ZZ interactions and together with the local addressability, a platform for one- | − | 1 2 the approximations γ1 δ− and γ2 JC δ− have way quantum computation [94,95]. The experimentalveri- been used. Since photons| | ≈ should be more| | ≈likely to tunnel fication that the desired entanglement has indeed been cre- to the next cavity than decay into free space, ΓC JC ated is thereby a non-trivial task in its own right [96]. 2 2 ≪ should hold in most cases. For ΓE JC g /δ , to hold, One way to generate cluster states is via the effective ≪ cavities with a high ratio g/ΓE are favourable. Since δ > Hamiltonian (42) [27]. A similar rout has been presented 2JC, the two requirements together imply that the cavities in [30]. To this end, all atoms are initialised in the states should have a high cooperativity factor ξ. ( aj + bj )/√2, which canbe donevia a STIRAP process [76],| i and| theni evolved under the Hamiltonian (42) for t = π/4Jz. 7.3. Magnetic dynamics of polarised light Figure 24 shows the von Neumann entropy of the re- duced density matrix of one effective spin EvN and the In contrast to identifying the internal states of the atoms purity of the reduced density matrix of the effective spin with spin polarisations, magnetic dynamics can also be chain Ps for a full three cavity model. Since EvN log22 ≈ studied by considering cavities which support two distinct for t 50µs while the state of the effective spin model ≈ 2 modes of light with different polarisations, and . In [28] remains highly pure (Ps = tr[ρ ] > 0.995) the degree of the dynamics of the photon polarisations in↑ such a↓ scheme entanglement will be very close to maximal. Thus the lev- was analysed. The authors consider cavities that are doped els aj and bj have been driven into a state which is, up | i | i with atoms of a V level structure, where two distinct inter- to local unitary rotations, very close to a three-qubitcluster nal transitions only couple to photons of different polari- state. sations. The interactions of the photons with the atoms are A second way to generate cluster states is to imple- chosen such that double occupancies are suppressed by the ment controlled phase gates between the qubits in differ- blockade effects discussed in section 5 and the dynamics ent cavities [52]. First all qubits are initialised in the state of the photons with creation operators a† and a† for the + = ( 1 + 0 )/√2, where 1 and 0 refer to states with | i | i | i | i | i two polarisations obeys the effective Hamiltonian↑ ↓ [28] one or zero excitationsin a cavity. The phase gates are then implemented by isolating chains of three qubits via a Stark x x y y z z shift that turns surrounding qubits off resonance. The iso- H = J(Sj Sj′ + Sj Sj′ )+ J ′Sj Sj′ , (43) ′ lated three-qubit chains evolve for a time t = π/(2√2J)

Copyright line will be provided by the publisher Laser & Photon. Rev. 1, No. 1 (2008) 21

Figure 25 A 2d cavity array. Con- secutive application of gates A, B, C and D, depicted in parts (a), (b), (c) and (d), each isolate chains of 3 qubits, realising controlled- phases and SWAPs (swapping quan- tum/classical gates) between the qubits at either end of the chain (the dashed lines indicate where the controlled-phase gates have been ap- plied). Reprinted with permission from [52].

(J is the photon tunnelling rate) and the central qubit is with small volume micro-cavities have emerged in recent subsequently measured. Local unitaries on the remaining years. These are routinely produced in large arrays and by two qubits that depend on the outcome of this measure- construction can allow for efficient photon tunnelling be- ment complete the controlled phase gate. Figure 25 shows tween neighbouring cavities. In addition, the strong cou- the entire sequence of gates. pling regime with cooperativity factors much larger than Further ways to create entangled states include mea- unity has been achieved in several of these devices. These surement based schemes. One version uses two atoms in two features together make such micro-cavity arrays very two cavities that interact via the exchange of photons [55]. promising candidates for the generation of the effective Both atoms can emit photons of two polarisations (left cir- many-body models discussed in this review. In the follow- cular and right-circular), depending on which excited state ing we discuss some platforms in a bit more detail. they are in. Hence the two atoms are prepared such that they each emit a single photon, one with left-circular and one with right-circular polarisation. After a time t = π/4J 8.1. Photonic band gap cavities (J being the photon tunnelling rate), when photons have tunnelled between the cavities with probability 1/2, the One promising device consists of quantum dots coupled to two atoms are employed to perform a non-demolitionmea- photonic band gap defect nanocavities. Photonic crystals surement to determine the number of photons in each cav- are structures with periodic dielectric properties which af- ity. Provided each cavity contains one photon the state is fect the motion of photons, in a similar manner as the peri- collapsed to a two photon polarisation entangled state thus odicity of a semiconductor affects the motion of electrons providing a heralded source of entangled photons. [97]. Such structures can exhibit band gaps in frequency Apart from quantum information motivated applica- space in which no photon of particular frequencies can tions coupled cavity arrays with local addressability can propagate in the material (for an introductory review see also be used to study novel physical setups. Among those [98]). are interface between regions in different condensed mat- A nanocavity can be created in a photonic-band-gap ter phases, see section 9. material by producing a localised defect in the structure of the crystal, in such a way that light of a particular fre- quency cannot propagate outside the defect area. Large ar- 8. Promising Platforms for Realisations rays of such nanocavities have been produced [99]. Photon hopping between neighbouring cavities has been observed In this section, we present some cavity QED setups which in the microwave and optical domains [100,101], where its have favourable features for the realisation of the effec- description in the form as given in equation (5) was found tive many-body models described in this review. These to be accurate. Quantum dots, semiconductors whose exci- platforms are thus promising candidates for experimental tations are confined in all three directions, behave in many demonstrations of the theory presented here. The selection respects like atoms and can be effectively addressed by we present here obviously reflects our personal view. It is lasers [102]. Quantum dots can be created inside photonic of course possible that devices which we do not explicitly crystal nanocavities [103] and made to interact with the mention here may in the end provide the first realisation of cavity mode to form a standard cavity QED system. The polaritonic many-body systems as discussed in previous strong coupling regime with cooperativity factors up to sections. ξ 10 has already been achieved in such systems [61]. A major part of all cavity QED experiments have been Cooperativity≈ factors of ξ 106 are predicted to be re- performed with Fabry-P´erot cavities formed by rather mas- alisable [81]. Due to the extremely∼ small volume of the sive mirrors. With such mirrors, coupling two cavities such nanocavity, the Jaynes-Cummings coupling coefficient can that photons can tunnel between them seems quite diffi- be extremely large [60,104], while spontaneous emission cult. On the other hand several new cavity QED structures rates are very low. Figure 26 shows a quantum dot in a

Copyright line will be provided by the publisher 22 Michael J. Hartmann et. al.: Quantum Many-Body Phenomena in Coupled Cavity Arrays

Figure 26 A quantum dot in a photonic band gap cav- Figure 27 SEM picture of a ity. a, AFM topography of a toroidal microresonator with photonic crystal nanocavity an intrinsic cavity Q of 108. aligned to a hill of material Reprinted with permission on the surface arising from from [57]. a quantum dot buried 63 nm below. b, Electric field in- tensity of the photonic crys- tal cavity mode showing that atoms are placed close to the cavity surface and interact the location of the buried with the evanescent field. This has been seen in an exper- QD, indicated by the teal iment where atoms fell by the cavity surface [58] and a dot, overlaps the field max- photon blockade effect generated by one atom has been imum. Reprinted with per- mission from [61]. demonstrated [108]. Theoretical predictions [81] show that very large cooperativity factors of ξ 106 are expected to be achieved in such systems. ∼ For realising effective Hamiltonians such as (7), (30) photonic band gap cavity. For the applications discussed or (36), some aspects still remain challenging for micro in this review current cavity decay rates of photonic band disc and micro toroidal cavities. One aspect is that atoms gap cavities remain a limiting factor. This is mainly due to need to be trappedfor a sufficiently long time close enough the fact that most structures employ two dimensional pho- to the cavity surface to maintain a sufficient strong cou- tonic crystals as placing a quantum dot into a three dimen- pling regime. A second aspect is that all cavities of an ar- sional structure is difficult. The difficulties to control the ray need to be tuned into resonance with each other with properties of individual quantum dots and to place them at enough accuracy. This can be achieved by varying their desired locations could even pose a more significant prob- temperature and thus shrinking or expanding them, but is lem. Current experiments pick a quantum dot from an en- still a technical challenge. Finally there can be modifica- semble and then fabricate a photoniccrystal tailored to this tions of the photon hopping as compared to the form as- dot around it. Such an approach can obviously not be used sumed in equation (5). A challenge here would be to pre- for producing a whole array of sufficient quality. vent that significant photon population is localised in the One mightthink abouttrappingatoms in photonicband fibre instead of in the cavities themselves. Nonetheless, gap cavities instead of letting them interact with quantum even if the fibre contains appreciable excitations, interest- dots. Since the frequencies of currently available photonic ing many-body physics is expected to happen as the fibre band gaps do not fit to atomic transition frequencies, such would merely play the role of an additional lattice site with a strategy would require the development of new photonic different local terms than the other sites. band gap materials. Figure 28 shows an array of 10 micro disc cavities that are all coupled by one taperedoptical fibre and mountedon an atom chip.The purposeofthe atomchip is to trapatoms 8.2. Silicon structures: micro discs and micro close to the cavity surfaces and, in that way, maintain a strong coupling regime for substantial times. toroids

A second class of micro-cavities that are routinely pro- duced in large arrays and which have very high Q-factors 8.3. Fibre based cavities are silicon structures of either a disc [59] or a toroidal [57] shape, see figure 27. These cavities are closed silicon discs Strong coupling experiments with many atoms in one cav- or toroids that typically are carried by a thin pillar in their ity and at the same time large cooperativityfactors for each centre. They are routinely produced in large arrays [105]. individual atom have been performed with fibre Fabry- The light is trapped in whispering gallery modes that are P´erot cavities [64]. These cavities are formed by two tips localised close to the outer surface of the structure and of optical fibres that are facing each other and are mounted have a small mode volume [106] . Very efficient and highly on an atom chip to trap an manipulate the atoms that inter- tunable coupling of photons into and from those cavities is act with the cavity modes. Current chips contain two such possible via tapered optical fibres that are placed close to cavities as shown in figure 29 and cooperativity factors as their surface [107]. Here the evanescent field of both cav- high as ξ = 145 have been achieved [64]. ity and fibre overlap and photons can hop between both. Photons are coupled into and from the cavity directly The photon hopping rate can thereby be controlled very via the optical fibres. Hence as for fibre coupled micro accurately via the separation between cavity and fibre. discs or toroids this may lead to modifications of the form These cavities can interact with atoms in a strong cou- of the photon hopping as compared to equation (5). On the pling regime with cooperativity factors ξ 50 [58] if the other hand, the atoms are here trapped on the atom chip ∼

Copyright line will be provided by the publisher Laser & Photon. Rev. 1, No. 1 (2008) 23

Figure 30 Left: Single on-chip Fabry-P´erot cavity formed by an etched micro-mirror and a tip of an optical fibre. Right: Two di- mensional array of cavity micro-mirrors. Reprinted with permis- sion from [110]. Figure 28 (a) Top-view optical image of a fibre taper aligned with and array of ten micro-disks (the remaining four micro-disks are out of the field of view of the image). (b) SEM image of the array of ten micro-disks. (c) Fibre taper transmission vs wave- regular arrays of spherical mirrors can be created with con- length when the taper is aligned with the micro-disk array in (a). The letters match specific micro-disks in (a) with the correspond- trollable sizes and spacings. These cavities which have a ing resonances in (c). Reprinted with permission from [59]. very small mode volume are illustrated in figure 30. The extremelysmall mode volume leads to a large single-photon Rabi frequency, and light is conveniently coupled into and out of the cavities via the fibres themselves. These cavi- ties to date achieve a finesse F 5000, corresponding to Purcell factors (or cooperativities)∼ ξ 40. ∼ As with the other cavities that involve optical fibres modifications of the photon hopping as compared to equa- tion (5) are expected here as well. Current strong coupling experiments have also been only performed with falling atoms [63].

8.5. Superconducting stripline resonators Figure 29 a, Layout of an atom chip that contains two fibre Fab- ryP´erot cavities. b, Close-up view of the two fibre FabryP´erot (FFP) optical cavities that are mounted on the chip. Cavity modes A promising experimental platform that operates in the are drawn to scale in red. The BEC is produced in a magnetic trap microwave regime are cavity QED system formed by a and positioned in the FFP1 mode (’B’). c, Geometry of the FFP1 Cooper-pair box coupled to a superconducting transmis- cavity. d, Overlay of three CCD time-of-flight (TOF) absorption sion line resonator. A Cooper-pair box is formed by two images, showing the anisotropic expansion of a BEC having in- superconducting islands separated by a Josephson junc- teracted for 50 ms with the cavity field. The optical fibres are tion. Although it is a macroscopic device, it behaves in outlined for clarity. Reprinted with permission from [64]. many ways like an atom as the relevantstructure can be de- scribed by a two level quantum system. A superconducting transmission line cavity, on the other hand, is a quasi-one- and can thus be held in a location of strong interaction with dimensional coplanar waveguide resonator, formed on a the cavity modes for a long time. chip with a resonance frequency in the microwave range. The strong-couplingregime in such a system has been achieved [62] and recently two Cooper-pairboxes have been strongly coupled to the cavity mode [111,112]. Very large cooper- 8.4. On-chip Fabry-Perot´ cavities ativity factors have been achieved in this set-up [62,112] and the nonlinearity of the Jaynes-Cummings model has Micro-fabricated arrays of Fabry-P´erot type cavities with been measured [113]. Although coupling between differ- large cooperativity factors have been presented in [109]. ent cavities has not been realised, up to ten Cooper-pair Each plano-concave resonator consists of an etched spher- boxes interacting in the strong coupling regime with the ical mirror and a coated single-mode optical fibre, allow- resonant cavity mode has been predicted to be achievable ing open access to the intra-cavity field mode. Here, large using the architecture of Ref. [112].

Copyright line will be provided by the publisher 24 Michael J. Hartmann et. al.: Quantum Many-Body Phenomena in Coupled Cavity Arrays

2

1.5

j 1 n

0.5

0 3000 2000 40 30 1000 20 10 0 t 0 site

Figure 31 a, A superconducting niobium coplanar waveguide Figure 32 The particle density across the chain as function of resonator of 24 mm length. It is coupled by a capacitor at each time. Sites 1 to 20 start in a Mott insulator regime and the rest in end of the resonator (see b) to an input and output feed line. b, a superfluid regime. , , , ˜ , , The capacitive coupling to the input and output lines. c, False U = 1.0 J = 0.1 µ = 0 U = 0.2 µ˜ = 0 ˜ and . Reprinted with permission from [33]. colour electron micrograph of a Cooper pair box (blue) fabri- J = 1.0 JI = 0.1 cated onto the silicon substrate (green) into the gap between the centre conductor (top) and the ground plane (bottom) of a res- onator (beige) using electron beam lithography and double angle evaporation of aluminium. Reprinted with permission from [62]. are shown as a function of time. The migration takes place whenever the difference between the chemical potentials of both regions is less than the maximal energy of any eigenmode of the superfluid [33]. 9. Applications

Structures of coupled arrays of micro-cavities have sev- 10. Summary eral interesting potential applications. Particularly useful in this respect is the individual addressability of single In this review, we have tried to give an overview of the cur- cavities with optical lasers. The feasibility of single site rent status of approaches to engineer effective many body addressing is based on the size of the cavities which can dynamics in arrays of cavities that are coupled via photon easily be larger than an optical wavelength. Furthermore tunnelling. The schemes that have been developed so far the inter-cavity coupling could be realised by optical fibres can be divided into two main categories, effective Bose- which would move the cavities even further apart. Hubbard type models and effective spin models. Among the potential applications are the implementa- In Bose-Hubbard type models, quantum phase transi- tion of tools for quantum information processing. These tions between superfluid phases of polaritons or photons include the generation of entangled states such as cluster (coherent light) and Mott insulators can be studied. Here, states [27,52,54,55] (see section 7.4), the controlled prop- the experimentally challenging and hence physically new agation of polaritonic excitations in cavity arrays [53,35, and interesting part is to create a Mott insulator in which 36,37,50], quantum gates [51] and the generation of single polaritons or photons are by their effective mutual repul- photons in parallel [56]. sion forced to stay in one lattice site (cavity) and not able Perhaps most importantly, the local controllability of- to move. This situation would in correspond to a crystal fered by coupled cavity arrays can also give rise to novel formed by light or ”frozen light”. physical effects and we will mention here one example of A key issue on the experimental side is to generate these. Local manipulations of an effective many body sys- a Kerr nonlinearity or effective interaction between po- tem could be applied to engineer a boundary between a laritons/photons that is significantly larger than the decay Mott insulator and a superfluid region of a Bose-Hubbard mechanisms present in the cavity array. model at unit filling. Initially both regions are decoupled Existing approaches for effective spin models mainly and cooled to their respective ground states. This situation employ two metastable levels of an atom trapped in each can for example be prepared in a cavity array that is oper- cavity to represent the two spin polarisations and . Ra- ated in the way as discussed in section 4 and thus hosts a man transitions mediated by external driving↑ lasers↓ that polaritonic Bose-Hubbard Hamiltonian (7). After switch- give rise to interactions between neighbouring cavities via ing on a small tunnelling rate between both regions, all the exchange of virtual photons are then used to generate particles of the Mott region migrate to the superfluid area. effective spin-spin interactions. This behaviour is shown for a one dimensional chain in As for Bose Hubbard type models a main experimen- figure 32, where the particle densities for all lattice sites tal challenge is to make atom-photon interactions signif-

Copyright line will be provided by the publisher Laser & Photon. Rev. 1, No. 1 (2008) 25

icantly stronger than the decay mechanisms in the cavity [8] I. Bloch, J. Dalibard and W. Zwerger, Many-Body Physics array. with Ultracold Gases, Rev. Mod. Phys. 80, 885 (2008) Effective many body system in a coupled cavity arrays [9] D. Jaksch, C. Bruder, J.I. Cirac, C.W. Gardiner and P. Zoller, are of great interest from variouspoints of view. First of all Cold Bosonic Atoms in Optical Lattices, Phys. Rev. Lett. 81, they open up a way to study condensedmatter physicswith 3108 (1998). light and hence to engineer new quantum states of light. [10] M. Greiner, O. Mandel, T. Esslinger, T.W. H¨ansch and I. Furthermore they offer a great versatility for applications Bloch, Quantum phase transition from a superfluid to a Mott insulator in a gas of ultracold atoms, Nature 415, 39 (2002). as quantum simulators due to the possibilities to manip- [11] B. Paredes, A. Widera, V. Murg, O. Mandel, S. Flling, I. ulate them via external driving lasers. Most importantly Cirac, G.V. Shlyapnikov, T.W. H¨ansch and I. Bloch, Tonks- they would, by construction, allow for control, manipu- Girardeau gas of ultracold atoms in an optical lattice, Nature lations and measurements of individual lattice sites. The 429, 277 (2004). situation is somehow complementary to cold atoms in op- [12] M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, A. tical lattices, where very large lattices can be created and Sen, U. Sen, Ultracold atomic gases in optical lattices: mim- experiments provide good measurement access to global icking condensed matter physics and beyond, Advances in quantities such as the phase coherence of the superfluid. Physics 56, 243 (2007). A polaritonic many body system in a cavity array on [13] Y. Makhlin, G. Sch¨on, and A. Shnirman, Quantum-state en- the other hand will for the foreseeable future only com- gineering with Josephson-junction devices, Rev. Mod. Phys. prise much smaller lattices, but provide good access for 73, 357 (2001). the measurement of local quantities. In contrast to optical [14] H.S. van der Zant, F.C. Fritschy, W.J. Elion, L.J. Geerligs lattices it would thus allow to directly verify the integer oc- and J.E. Mooij, Field-induced superconductor-to-insulator cupation number of each lattice site in the Mott insulator transitions in Josephson-junction arrays, Phys. Rev. Lett. 69, regime. 2971 (1992). [15] A. van Oudenaarden and J.E. Mooij, One-Dimensional Mott The local addressability of coupled cavity arrays gives Insulator Formed by Quantum Vortices in Josephson Junc- rise to a number of novel application. These contain quan- tion Arrays, Phys. Rev. Lett. 76, 4947 (1996). tum information processing application but also novel phys- [16] F. Mintert and C. Wunderlich, Ion-Trap Quantum Logic ical situations that give rise to features such as particle mi- Using Long-Wavelength Radiation, Phys. Rev. Lett. 87, gration. 257904 (2001). [17] D. Porras and J.I. Cirac, Effective Quantum Spin Systems with Trapped Ions, Phys. Rev. Lett. 92, 207901 (2004). Acknowledgements This work is part of the QIP-IRC supported [18] D. Porras and J.I. Cirac, Bose-Einstein Condensation and by EPSRC (GR/S82176/0), the Integrated Project Qubit Applica- Strong-Correlation Behavior of Phonons in Ion Traps, Phys. tions (QAP) supported by the IST directorate as Contract Num- Rev. Lett. 93, 263602 (2004). ber 015848’, the EU STREP project HIP and was supported by [19] A. Retzker, R. Thompson, D. Segal and M.B. Plenio, Double the EPSRC grant EP/E058256/1, the Alexander von Humboldt wells, scalar fields and quantum phase transitions in ions Foundation, the Conselho Nacional de Desenvolvimento Cient´ıfico traps, arXiv:0801.0623 e Tecnol´ogico (CNPq), the Royal Society and the DFG Emmy [20] A. Friedenauer, H. Schmitz, J.T. Gl¨uckert, D. Porras and T. Noether grant HA 5593/1-1. Sch¨atz, Simulating the Quantum Magnet, arXiv:0802.4072 [21] M.J. Hartmann, F.G.S.L. Brand˜ao and M.B. Plenio, Strongly interacting polaritons in coupled arrays of cavities, Nature References Phys. 2, 849 (2006). [22] D.G. Angelakis, M.F. Santos and S. Bose, Photon-blockade- induced Mott transitions and XY spin models in coupled [1] H.J. Metcalf and P. van der Straten, Laser Cooling and Trap- cavity arrays, Phys. Rev. A 76, R031805 (2007). ping. Springer, New York, (1999). [23] A.D. Greentree, C. Tahan, J.H. Cole and L.C.L. Hollen- [2] A.J. Leggett, Bose-Einstein condensation in the alkali gases: berg, Quantum phase transitions of light, Nature Phys. 2, Some fundamental concepts, Rev. Mod. Phys. 73, 307 856 (2006). (2001). [24] F. Illuminati, Quantum optics: Light does matter, Nature [3] D. Leibfried, R. Blatt, C. Monroe, D. Wineland, Quantum Phys. 2, 803 (2006). dynamics of single trapped ions, Rev. Mod. Phys., 75, 281 [25] M.J Hartmann and M.B Plenio, Strong Photon Nonlinear- (2003). ities and Photonic Mott Insulators, Phys. Rev. Lett. 99, [4] J.M. Raimond, M. Brune and S. Haroche, Manipulating 103601 (2007). quantum entanglement with atoms and photons in a cavity, [26] F.G.S.L. Brand˜ao, M.J. Hartmann and M.B. Plenio, Light- Rev. Mod. Phys. 73, 565 (2001). shift-induced photonic nonlinearities, New J. Phys. 10, [5] A. Zeilinger, Experiment and the foundations of quantum 043010 (2008). physics , Rev. Mod. Phys. 71, S288 (1999). [27] M.J. Hartmann, F.G.S.L. Brand˜ao and M.B. Plenio, Effec- [6] F. De Martini and C. Monroe (editors), Experimental Quan- tive Spin Systems in Coupled Microcavities, Phys. Rev. Lett. tum Computation and Information, IOS Press, Amsterdam, 99, 160501 (2007). (2002). [28] A.-C. Ji, X.C. Xie and W.M. Liu, Quantum Magnetic Dy- [7] R.P. Feynman, Simulating physics with computers, Inter. J. namics of Polarized Light in Arrays of Microcavities, Phys. Theor. Phys. 21, 467 (1982). Rev. Lett. 99, 183602 (2007).

Copyright line will be provided by the publisher 26 Michael J. Hartmann et. al.: Quantum Many-Body Phenomena in Coupled Cavity Arrays

[29] A. Kay and D.G. Angelakis, Reproducing spin lat- [49] D. Rossini, R. Fazio and G. Santoro, Photon and polariton tice models in strongly coupled atom-cavity systems, fluctuations in arrays of QED-cavities, Europhys. Lett. in arXiv:0802.0488 press, arXiv:0806.0942. [30] P. Li, Q. Gong and G. Guo, Effective generation of Ising [50] M.X. Huo, Y. Li, Z. Song and C.P. Sun, Atomic entan- interaction and cluster states in coupled micro-cavities, glement versus visibility of photon interference for quan- arXiv:0712.3864. tum criticality of a hybrid system, Phys. Rev. A 77, 022103 [31] J. Cho, D.G. Angelakis and S. Bose, Simulation of high-spin (2008). Heisenberg chains in coupled cavities, arXiv:0802.3365 [51] D.G. Angelakis, M.F. Santos, V. Yannopapas and A. Ekert, [32] J. Cho, D.G. Angelakis and S. Bose, Fractional Quantum A proposal for the implementation of quantum gates with Hall State in Coupled Cavities, arXiv:0807.1802 photonic-crystal coupled cavity waveguides, Phys. Lett. A. [33] M.J. Hartmann and M.B. Plenio, Migration of Bosonic Par- 362, 377 (2007). ticles across a Mott Insulator to a Superfluid Phase Interface, [52] D.G. Angelakis and A. Kay, Weaving light-matter qubits Phys. Rev. Lett. 100, 070602 (2008). into a one way quantum computer, New J. Phys. 10, 023012 [34] C.D. Ogden, E.K. Irish and M.S. Kim, Dynamics in a (2008). coupled-cavity array, arXiv:0804.2882. [53] S. Bose, D.G. Angelakis and D. Burgarth, Transfer of a Po- [35] L. Zhou, J. Lu and C.P. Sun, Coherent control of photon laritonic Qubit through a Coupled Cavity Array, Journ. of transmission: Slowing light in a coupled resonator waveg- Mod. Opt. 54, 2307 (2007). uide doped with Λ atoms, Phys. Rev. A 76, 012313 (2007). [36] F.M. Hu, L. Zhou, T. Shi and C.P. Sun, Coupled cavity QED [54] D.G. Angelakis, Mancini and S. Bose, Steady state for coherent control of photon transmission: Green-function entanglement between hybrid light-matter qubits, approach for hybrid systems with two-level doping Phys. arXiv.org:0711.1830 Rev. A 76, 013819 (2007). [55] J. Cho, D.G. Angelakis and S. Bose, Heralded generation [37] L. Zhou, Y.B. Gao, Z. Song and C.P. Sun, Coherent out- of two-photon polarization entanglement with coupled cav- put of photons from coupled superconducting transmission ities, Phys. Rev. A 78, 022323 (2008). line resonators controlled by charge qubits, Phys. Rev. A 77, [56] N. Na and Y. Yamamoto, Generation of indistinguishable 013831 (2008). single photons and polarization-entangled photon-pairs via [38] M.J. Hartmann, F.G.S.L. Brand˜ao and M.B. Plenio, A po- polaritonic superfluid to Mott insulator quantum phase tran- laritonic two-component BoseHubbard model, New J. Phys. sition, arXiv:0804.1829 10, 033011 (2008). [57] D.K. Armani, T.J. Kippenberg, S.M. Spillane and K.J. Va- [39] D. Rossini and R. Fazio, Mott-Insulating and Glassy Phases hala, Ultra-high-Q toroid micro-cavity on a chip, Nature of Polaritons in 1D Arrays of Coupled Cavities, Phys. Rev. 421, 925 (2003). Lett. 99, 186401 (2007). [58] T. Aoki, B. Dayan, E. Wilcut, W.P.Bowen, A.S. Parkins, [40] M. Aichhorn, M. Hohenadler, C. Tahan and P. B. Little- H.J. Kimble, T.J. Kippenberg and K.J. Vahala, Observation wood, Quantum Fluctuations, Temperature, and Detuning of strong coupling between one atom and a monolithic mi- Effects in Solid-Light Systems, Phys. Rev. Lett. 100, 216401 croresonator, Nature 443, 825 (2006). (2008). [59] P.E. Barclay, K. Srinivasan, O. Painter, B. Lev and H. [41] E.K. Irish, C.D. Ogden, and M.S. Kim, Polaritonic charac- Mabuchi, Integration of fiber-coupled high-Q SiNx micro- teristics of insulator and superfluid states in a coupled-cavity disks with atom chips Appl. Phys. Lett. 89, 131108 (2006). array, Phys. Rev. A 77, 033801 (2008). [60] Y. Akahane, T. Asano, B.-S.Song and S. Noda High-Q pho- [42] M.I. Makin, J.H. Cole, C. Tahan, L.C.L. Hollenberg and tonic nanocavity in a two-dimensional photonic crystal, Na- A.D. Greentree, Quantum phase transitions in photonic cav- ture 425, 944 (2003). ities with two-level systems, Phys. Rev. A 77, 053819 [61] K. Hennessy, A. Badolato, M. Winger, D. Gerace, M. (2008). Atat¨ure, S. Gulde, S. F¨alt, E.L. Hu and A. Imamo˘glu, Quan- [43] J. Zhao, A.W. Sandvik and K. Ueda, Insulator to superfluid tum nature of a strongly coupled single quantum dot – cavity transition in coupled photonic cavities in two dimensions, system, Nature 445, 896 (2007). arXiv:0806.3603 [44] N. Na, S. Utsunomiya, L. Tian and Y. Yamamoto, Strongly [62] A. Wallraff, D.I. Schuster, A. Blais, L. Frunzio, R.- correlated polaritons in a two-dimensional array of photonic S. Huang, J. Majer, S. Kumar, S.M. Girvin and R.J. crystal micro-cavities, Phys. Rev. A 77, 031803(R) (2008). Schoelkopf, Strong coupling of a single photon to a su- [45] D.E. Chang, V. Gritsev, G. Morigi, V. Vuletic, M.D. Lukin perconducting qubit using circuit quantum electrodynamics, and E.A. Demler, Crystallization of strongly interacting Nature 431, 162 (2004). photons in a nonlinear optical fiber arXiv:0712.1817 [63] M. Trupke, J. Goldwin, B. Darqui´e, G. Dutier, S. Eriksson, [46] M. Paternostro, G.S. Agarwal and M.S. Kim, Solitons in in- J. Ashmore and E. A. Hinds, Atom Detection and Photon teracting Dicke models of coupled cavities with two-level Production in a Scalable, Open, Optical Microcavity, Phys. systems, arXiv:0707.0846 Rev. Lett. 99, 063601 (2007). [47] S.-C. Lei and R.-K. Lee, Quantum phase transitions of light [64] Y. Colombe, T. Steinmetz, G. Dubois, F. Linke, D. Hunger in the Dicke-Bose-Hubbard model, Phys. Rev. A 77, 033827 and J. Reichel, Strong atomfield coupling for BoseEinstein (2008). condensates in an optical cavity on a chip, Nature 450, 272 [48] Y.-F. Xiao, J. Gao, X.-B. Zou, J.F. McMillan, X. Yang, Y.- (2007). L. Chen, Z.-F. Han, G.-C. Guo and C.W. Wong, Coupled [65] R.J. Thompson, G. Rempe and H.J. Kimble, Observation of quantum electrodynamics in photonic crystal nanocavities, normal-mode splitting for an atom in an optical cavity, Phys. arXiv.org:0707.2632 Rev. Lett. 68, 1132 (1992).

Copyright line will be provided by the publisher Laser & Photon. Rev. 1, No. 1 (2008) 27

[66] Q.A. Turchette, C.J. Hood, W. Lange, H. Mabuchi and [84] L.-M. Duan, E. Demler and M. D. Lukin, Controlling Spin H.J. Kimble, Measurement of Conditional Phase Shifts for Exchange Interactions of Ultracold Atoms in Optical Lat- Quantum Logic, Phys. Rev. Lett. 75, 4710 (1995). tices, Phys. Rev. Lett. 91, 090402 (2003). [67] K.M. Birnbaum, A. Boca, R. Miller, A.D. Boozer, T.E. [85] M.A. Cazalilla and A.F. Ho, Instabilities in Binary Mixtures Northup and H.J. Kimble, Photon blockade in an optical of One-Dimensional Quantum Degenerate Gases, Phys. cavity with one trapped atom, Nature 436, 87 (2005). Rev. Lett. 91, 150403 (2003). [68] N. Stefanou and M. Modinos, Impurity bands in photonic [86] A. Imamo˘glu, High Efficiency Photon Counting Using insulators, Phys. Rev. B 57, 12127 (1998). Stored Light, Phys. Rev. Lett. 89, 163602 (2002). [69] A. Yariv, Y. Xu, R.K. Lee and A. Scherer, Coupled-resonator [87] D.F.V.James and P.G. Kwiat, Atomic-Vapor-Based High Ef- optical waveguide: a proposal and analysis, Optics Lett. 24, ficiency Optical Detectors with Photon Number Resolution, 711 (1999). Phys. Rev. Lett. 89, 183601 (2002). [70] M. Marzani and D. Vanderbilt, Maximally localized gener- [88] U. Schollw¨ock, The density-matrix renormalization group, alized Wannier functions for composite energy bands, Phys. Rev. Mod. Phys. 77, 259 (2005). Rev. B 56, 12847 (1997). [89] J.K. Freericks and H. Monien, Strong-coupling expansions [71] M.J. Hartmann, M. Reuter, and M.B. Plenio, Excitation and for the pure and disordered Bose-Hubbard model, Phys. entanglement transfer versus spectral gap, New J. Phys. 8, Rev. B 53, 2691 (1996). 94 (2006). [90] T. Moriya and K. Ueda, Antiferromagnetic spin fluctuation [72] M.B. Plenio, J. Hartley and J. Eisert, Dynamics and manip- and superconductivity, Rep. Prog. Phys. 66, 1299 (2003). ulation of entanglement in coupled harmonic systems with [91] H.J. Briegel and R. Raussendorf, Persistent Entanglement many degrees of freedom, New J. Phys. 6, 36 (2004). in Arrays of Interacting Particles, Phys. Rev. Lett. 86, 910 [73] M.P.A. Fisher, P.B. Weichman, G. Grinstein and D.S. Fisher, (2001). Boson localization and the superfluid-insulator transition, [92] A.S. Sørensen and K. Mølmer, Entangling atoms in bad cav- Phys. Rev. B 40, 546 (1989). ities, Phys. Rev. A 66, 022314 (2002). [93] D.F.V. James, Quantum Computation with Hot and Cold [74] M.P.A. Fisher and G. Grinstein and S. M. Girvin, Presence Ions: An Assessment of Proposed Schemes, Fortschr. Phys. of quantum diffusion in two dimensions: Universal resis- 48, 823 (2000). tance at the superconductor-insulator transition, Phys. Rev. [94] R. Raussendorf and H.J. Briegel, A One-Way Quantum Lett. 64, 587 (1990). Computer, Phys. Rev. Lett. 86, 5188 (2001). [75] J. Hubbard, Electron Correlations in Narrow Energy Bands, [95] D. Gross and J. Eisert, Novel Schemes for Measurement- Proc. Roy. Soc. A 276 238 (1963). Based Quantum Computation, Phys. Rev. Lett. 98, 220503 [76] M. Fleischhauer, A. Imamoglu and J.P. Marangos, Electro- (2006). magnetically induced transparency: Optics in coherent me- [96] K.M.R. Audenaert and M.B. Plenio, When are correlations dia, Rev. Mod. Phys. 77, 633 (2005). quantum? – verification and quantification of entanglement [77] H. Schmidt and A. Imamoglu, Giant Kerr nonlinearities ob- by simple measurements, New J. Phys. 8, 266 (2006). tained by electromagnetically induced transparency, Opt. [97] S.G. Johnson and J.D. Joannopoulos, Photonic Crystals: Lett. 21, 1936 (1996). The Road from Theory to Practice (Kluwer, 2002). [78] A. Imamoglu, H. Schmidt, G. Woods and M. Deutsch, [98] D.G. Angelakis, P.L. Knight and E. Paspalakis, Photonic Strongly Interacting Photons in a Nonlinear Cavity, Phys. crystals and inhibition of spontaneous emission: an intro- Rev. Lett. 79, 1467 (1997); P. Grangier, D.F. Walls and K.M. duction, Cont. Phys. 45, 303 (2004). Gheri, Comment: Strongly Interacting Photons in a Nonlin- [99] H. Altug and J. Vuckoic, Two-dimensional coupled photonic ear Cavity, Phys. Rev. Lett. 81, 2833 (1998); A. Imamoglu, crystal resonator arrays, Appl. Phys. Lett. 84, 161 (2004). H. Schmidt, G. Woods and M. Deutsch, Erratum: Strongly [100] M. Bayindir, B. Temelkuran, and E. Ozbay, Propagation of Interacting Photons in a Nonlinear Cavity, Phys. Rev. Lett. photons by hopping: A waveguiding mechanism through lo- 81, 2836 (1998). calized coupled cavities in three-dimensional photonic crys- [79] M.J. Werner and A. Imamoglu, Photon-photon interactions tals, Phys. Rev. B 61, R11855 (2000). in cavity electromagnetically induced transparency, Phys. [101] M. Bayindir, B. Temelkuran and E. Ozbay, Tight-Binding Rev. A 61, 011801(R) (1999). Description of the Coupled Defect Modes in Three- [80] K.M. Gehri, W. Alge and P. Grangier, Quantum analysis of Dimensional Photonic Crystals, Phys. Rev. Lett. 84, 2140 the photonic blockade mechanism, Phys. Rev. A 60, R2673 (2000). (1999). [102] F.H. Julien and A. Alexandrou, Controlling Artificial [81] S.M. Spillane, T.J. Kippenberg, K.J. Vahala, K.W. Goh, Atoms, Science 20, 1429 (1998). E. Wilcut and H.J. Kimble, Ultrahigh-Q toroidal microres- [103] A. Badolato, K. Hennessy, M. Atat¨ure, J. Dreiser, E. Hu, onators for cavity quantum electrodynamics, Phys. Rev. A P.M. Petroff, Atac Imamo˘glu, Deterministic Coupling of 71, 013817 (2005). Single Quantum Dots to Single Nanocavity Modes, Sience [82] F. Dimer, B. Estienne, A. S. Parkins and H. J. Carmichael, 308, 1158 (2005). Proposed realisation of the Dicke-model quantum phase [104] B.-S. Song, S. Noda, T. Asano and Y. Akahane, Ultra-high- transition in an optical cavity QED system, Phys. Rev. A Q photonic double-heterostructure nanocavity, Nature Ma- 75, 013804 (2007). terials 4, 207 (2005). [83] A. Recati, P. O. Fedichev, W. Zwerger and P. Zoller, Spin- [105] A.M. Armani, A. Srinivasan and K.J. Vahala, Soft Litho- Charge Separation in Ultracold Quantum Gases, Phys. Rev. graphic Fabrication of High Q Polymer Microcavity Arrays, Lett. 90, 020401 (2003). Nano Letters 7, 6 (2007).

Copyright line will be provided by the publisher 28 Michael J. Hartmann et. al.: Quantum Many-Body Phenomena in Coupled Cavity Arrays

[106] T.J. Kippenberg, S.M. Spillane, B. Min and K.J. Vahala, Demonstration of ultra-high-Q small mode volume toroid micro-cavities on a chip, Appl. Phys. Lett. 85, 6113 (2004). [107] S.M. Spillane, T.J. Kippenberg, O.J. Painter and K.J. Va- hala, Ideality in a Fiber-Taper-Coupled Microresonator Sys- tem for Application to Cavity Quantum Electrodynamics, Phys. Rev. Lett. 91, 043902 (2003). [108] B. Dayan, A.S. Parkins, T. Aoki, E.P. Ostby, K.J. Vahala and H.J. Kimble, A Photon Turnstile Dynamically Regu- lated by One Atom, Science 319, 1062 (2008). [109] M. Trupke, E.A. Hinds, S. Eriksson, E.A. Curtis, Z. Mok- tadir, E. Kukharenka and M. Kraft, Microfabricated high- finesse optical cavity with open access and small volume, Appl. Phys. Lett. 87, 211106 (2005). [110] http://www3.imperial.ac.uk/ccm/research/coldatoms/cavity [111] M.A. Sillanp¨a¨a, J.I. Park and R.W. Simmonds, Coherent quantum state storage and transfer between two phase qubits via a resonant cavity, Nature 449, 438 (2007). [112] J. Majer, J.M. Chow, J.M. Gambetta, J. Koch, B.R. John- son, J.A. Schreier, L. Frunzio, D.I. Schuster, A.A. Houck, A. Wallraff, A. Blais, M.H. Devoret, S.M. Girvin and R.J. Schoelkopf, Coupling superconducting qubits via a cavity bus Nature 449, 443 (2007). [113] J.M. Fink, M. G¨oppl, M. Baur, R. Bianchetti, P.J. Leek, A. Blais and A. Wallraff, Climbing the JaynesCummings lad- der and observing its √n nonlinearity in a cavity QED sys- tem, Nature 454, 315 (2008).

Copyright line will be provided by the publisher