arXiv:0803.3963v2 [cond-mat.other] 20 May 2008 ee tmwt states with atom level scribing promising. band-gap most seem periodic 5] dot [1, materials center/quantum (PBG) color and condensa [12] 11], [1, arrays cavities b “transmon” linked or supercondu boxes stripline microcavities of Cooper-pair arrays particular, small 10], In [2, fibers optical [1]. implem many possible are The states. tations initial create to Tempo- pumping n coherent using and studied be accessible, distances. may are dynamics long quantum functions equilibrium over correlation spatial transported excellent and as be ral serve can system they potential the “flying” that excellent of have part and probes, as experimental equilibrium from [4]. in far vices by and variation, construction local asymmetry tuning, for in dividual potential device the and sites and con- individual gration, advantages of some addressing have atom direct but cold cerning [9], of (CAOL) reminiscent lattices exhibit—are optical they photons of state The new explored. brings complete. fully system be not to the is yet of analogy Bose-Hubbard BH matter- studied mixed the widely Yet, the (BHM). to model similar s to be signa- in and (SF) reproducible particular be systems superfluid a should the transition, that (MI) physics, insulator clear Mott many-body is quantum It atoms of or ture context. electrons many-body like a much behave in to component polaritons) a cavity (as interactio photons of work -photon allow further cavities strong between Engineered by hopping followed 8]. and 7, 3], 6, 2, 5, [4, [1, this accomplished how shown be recently propos- might have Three arrays quantum cavity-QED on devices. sim- for based quantum quantum als testbed advanced of a other potential and as the ulators for both and exciting, dynamics is many-body photons of tem h rseto elzn ual,srnl orltdsy correlated strongly tunable, a realizing of prospect The eew ou ntesmls oi-ih oe 1,de- [1], model solid-light simplest the on focus we Here MI intriguing the for systems—so-named “Solid-light” L 2 hoyo odne atr aeds aoaoy Univers Laboratory, Cavendish Matter, Condensed of Theory pia ircvte ahcnann igetwo- single a containing each microcavities optical unu lcutos eprtr n euigEfcsin Effects Detuning and Temperature Fluctuations, Quantum xeietlyacsil eprtrsfrraitccav realistic for temperatures accessible experimentally ekyitrcigpoos h otsaewt est one density with parti state bosonic Mott The the of photons. nature a interacting by the weakly separated changing bands for hole allows and detuning particle by featur dominated non-mean-field are important exactly phase exhibit fluctuations dimensions quantum two account and into taking approach, ASnmes 13.c 34.q 82.h 42.50.Ct 78.20.Bh, 73.43.Nq, 71.36.+c, numbers: PACS h uefli oMt nuao rniini aiypolari cavity in transition insulator Mott to superfluid The aksAichhorn, Markus |↓i 1 nttt o hoeia hsc n srpyis Unive Astrophysics, and Physics Theoretical for Institute , |↑i 3 ozAlnHmlo n. 81N ara r,Alntn VA Arlington, Dr., Fairfax N. 3811 Inc., Hamilton Allen Booz eaae yenergy by separated 1 atnHohenadler, Martin ǫ The . cting uch on- de- en- te- so ns s- te y - 2 hre Tahan, Charles t couplings ity t fCmrde abig B H,Uie Kingdom United 0HE, CB3 Cambridge Cambridge, of ity Here aitna ed ( reads Hamiltonian rtr)gvsrs ofraino oaios(obndat number whose (combined excitations) polaritons photon of formation to rise gives erators) h euig ahcvt sdsrbdb h well-known the by described coupling is Hamiltonian cavity (JC) Each Jaynes-Cummings detuning. the scnevdadculst h hmclpotential chemical the to couples and conserved is en h oaio density polariton amplitude the with define hopping photon nearest-neighbor consider rmioae lses(with clusters isolated from nryadset and energy δ for choice mal eeec ytm h arxntto nldsobtl an orbitals includes frequencies, notation Matsubara approximate matrix bosonic to The is system. idea reference main The self-energy [14]. bee the BHM previously the [13]—has to electrons applied correlated strongly for o calculations concept of proof early systems. calculate CAOL in systems metric key such experi a to also in connection and time necessary first a the spectra, single-particle for lobes and Mott detuning discuss We the of temperatures), finite included. relevant experimentally are at fluctuation scale (also quantum length particular, finite In a explore on model. quantum to the limit a of and thermodynamic employ physics the the we 7], for Here 3, method [2, many-body [8]. systems case cavity one-dimensional few [1] the calculations of Existing mean-field diagonalization on exact methods. focused has numerical work to theoretical amenable photons related n-atceparameters one-particle H ˆ Ω[ H o rasi nlzduigtevrainlcluster variational the using analyzed is arrays ton i h aitoa lse prah(C)itoue first (VCA)—introduced Approach Cluster Variational The cor- strongly of model generic a represents (1) Hamiltonian JC lsfo uslclzdectn oplrtn to polaritons to excitons quasilocalized from cles ˆ s igepril xiainsetai h Mott the in spectra excitation Single-particle es. xssu otemperatures to up exists Σ nfiielnt cls hs igasi one in diagrams Phase scales. length finite on otgp ncnrs oBs-ubr models, Bose-Hubbard to contrast In gap. Mott ω 0 = ] = = 0 stecvt htneeg,and energy, photon cavity the is ǫ − g g |↑ . t st fW¨urzburg, Germany of rsity , ( ,3 2, a h i X Ω i,j h↑ i i † ω , Σ n ee .Littlewood B. Peter and i en h rn oeta.Tilself-energies Trial potential. grand the being Σ 0 a a i /g olw rmagnrlvrainlprinciple variational general a from follows i | i † + a fteifiiesse yta fafinite a of that by system infinite the of ~ r htnceto n niiainop- annihilation and creation photon are j [23], 20,USA 22203, 1 = ω + 0 ξ oi-ih Systems Solid-Light a X c ) i † k i T/g a = B n i L H ˆ { + n h atc osatt one. to constant lattice the and c = i t JC & ie)aeprmtie ythe by parametrized are sites) c g ǫ , h Σ − ( N 0 |↑ N c . ˆ ω , 03 p µ Σ = i p i h↓ i N H /L implying , ˆ = ˆ 0 c µ , p 2 JC i P use , , = ∆ αβ | h atom-photon The . c a } i ( i ( iω ftereference- the of a + g i † ω n a steui of unit the as |↓ 0 ) i h opti- The . − i + h↑ i µ |↑ ǫ 7.We [7]. defines i i h↑ i | effect ment a om- i † (1) i ) | t d n . ) s f , , 2 system Hamiltonian, i.e. Σ = Σ(ξc). For bosons [14], ∆ = 0 0.6 MI (n = 2) (a) 1D, T = 0, −1 c −1 Σ c L = 2 Ω=Ω + Tr ln G0  Tr ln (G ) . (2) 0.5 c − − Lc = 4 c c 0.4 L = 6 Here, , G , and Σ are the grand potential, Green’s func- g SF c Ω L = 8

µ / c tion, and self energy of an isolated cluster, and G0 is the non- 0.3 III c Σ MI (n = 1) interacting (g =0) Green’s function. G0, G , and are eval- 0.2 II uated at bosonic Matsubara frequencies iωn, and traces in- I clude frequency summation. The stationary solution is given 0.1 SF c 0 by ∂Ω/∂ξ = 0. Traces can be evaluated exactly using only 0 0.04 0.08 0.12 0.16 0.2 the poles of the Green’sfunctionbut not their weights [13, 14]. L = 1 x 1 −1 Σ (b) 2D, T = 0, ∆ = 0 c The poles ωm of G0 are obtained from a bosonic for- 0.5 Lc = 2 x 1 − L = 2 x 2 mulation of the Q-matrix method [15]. At temperature T > 0, SF c 0.4 Lc = 8 the required matrix diagonalization restricts Lc and Tmax. For g simplicity, we restrict ourselves to a single variational param- 0.3 III c c µ / II eter ξ = ω0. MI (n = 1) The full quantum dynamics are taken into account exactly 0.2 I on the length scale of the cluster Lc, andevenfor Lc =1 VCA 0.1 results are beyond the mean-field solution [1, 14]. The present SF 0 formulation cannot describe the properties of the SF phase, as 0 0.01 0.02 0.03 0.04 0.05 0.06 ˆ † the required symmetry-breaking term Hψ = ψ i(ai + ai ) (c) T = 0, ∆ = 0 P 0.5 1D, Lc = 8 (z = 2) cannot be cast into a single-particle operator. However, this 2D, L = 6 (z = 3) 0.4 c does not affect the accuracy of the phase boundary of the MI. 2D, L = 8 (z = 4) g c 0.3

The numerical effort is very moderate as compared to, e.g., µ / SF the density matrix renormalization group [8], and the VCA 0.2 MI (n = 1) provides T = 0 and T > 0 static and dynamic properties in 0.1 SF 0 one and two dimensions. 0 0.04 0.08 0.12 0.16 0.2 The T = 0 phase diagram in one (chain) and two dimen- t / g sions (square lattice) is shown in Fig. 1. There exists a series of Mott lobes with integer polariton density n = 0, 1,... FIG. 1: (Color online) Ground-state phase diagram: (a) 1D, (b) 2D int and (c) different coordination numbers z. Lines are guides to the eye. and compressibility κ ∂n/∂µ = 0 [1]. In contrast to re- cently proposed photonic≡ MI phases [4], a†a fluctuates even h i for constant, integer n. Inside the lobes, wherethe VCA yields and two dimensions. However, quantum fluctuation effects a solution, the system has an energygap Eg,p (Eg,h) foradding e e diminish quickly with increasing coordination number z (Fig- a particle (hole) equal to the vertical distance of µ from the ure 1(c)). In particular, the reentrant behavior with increas- upper (lower) phase boundary [24]. The spacing of the points ing t/g [16] and the cusplike tip indicative of the Berezinskii- 0 √ √ where adjacent lobes touch at µn = ω + g( n n + 1) Kosterlitz-Thouless transition [17] exist only for z =2. t/g = 0 decreases− quickly with increasing n, in contrast to Spectral properties play a key role in understanding con- the BHM where µ = Un. Convergence with cluster size L n c densed matter systems, but are notoriously difficult to cal- is surprisingly fast. The 1D data agree well with exact results, culate accurately. In cavities, the occupation and spectra although the VCA slightly underestimates t∗ (the value of t at can be directly monitored through luminescence spectra, with the lobe tip, t∗/g =0.2 in [8]). Figure 1(b) represents the first angular emission translating into momentum k, much more accurate (non-mean-field) phase diagram in 2D, arguably the straightforward than in CAOL. most important case for experimental realizations. The VCA yields the single-particle spectral function For t/g = 0, the MI states are ψn , where the n- Qi | ii A(k,ω) = π−1Im G(k,ω), shown in Fig. 2 for the pa- polariton eigenstate of Hˆ JC (the ,n branch in [1]) is a su- − |− i rameters marked in Fig. 1. Important for detection in ex- perposition of photonic ( ,n , with n photons) and excitonic periment, the MI state is characterized by cosinelike parti- states ( ,n 1 ), |↓ i |↑ − i cle and hole bands, separated by the Mott gap Eg, minimal at k = 0, which decreases with increasing t/g, and even- ψn = η(∆,n) ,n + φ(∆,n) ,n 1 . (3) | i |↓ i |↑ − i tually closes in the SF phase. The phase boundaries µ± Lobes with n> 1 are much smaller due to the effective polari- are related to the particle/hole dispersions εp/h(k = 0,t) via ton repulsion decreasing with n, and we focus on the n = 1 µ+(t) = εp(0,t) and µ−(t) = εh(0,t), and the Mott gap is case for which quantum effects are strongest, and which can given by Eg(t)= εp(0,t) εh(0,t) [18]. − be easily initialized experimentally. For t t∗, particles are much lighter than holes (e.g., ≪ Figures 1(a) and (b) show significant deviations from the tmp 0.33, tmh 0.81 in Fig. 2(c); obtained from parabolic parabolic lobes predicted by mean-field theory in both one fits).≈ The particle/hole≈ bandwidth scales almost linearly with 3

1.00 (a) 1D, I (b) 1D, II (c) 1D, III 0.5 0.80 (a) 1D, T = 0, ∆ / g = 2 0.60 0 Lc = 2 g

) / g 0.40 L = 4

µ -0.5 c µ /

- MI (n = 1) 0.20 Lc = 6 ω ( 0.00 -1 -0.20 -1.5 -π 0 π -π 0 π -π 0 π 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.80 0.68 (d) 2D, I (e) 2D, II (f) 2D, III 0.60 0.66 (b) 1D, T = 0, ∆ / g = -2 0.40 0.64 g ) / g

µ 0.20 µ / - MI (n = 1)

ω 0.62 ( 0.00 -0.20 0.6

π π π π π π π π π 0.58 (0,0) ( , ) ( ,0) (0,0) (0,0) ( , ) ( ,0) (0,0) (0,0) ( , ) ( ,0) (0,0) 0 0.005 0.01 0.015 0.02 0.025 k k k t / g 4 (c) ∆= 2g (d) ∆ = -2g 3 0.1 FIG. 2: Single-particle spectra at T = 0 in 1D [(a)–(c), Lc = 6] and × 2D [(d)–(f), Lc = 2 2]. I–III refer to the marks in Fig. 1, i.e., (a) 2 ) / g t/g = 0.16, (b) 0.08, (c) 0.01, (d) 0.05, (e) 0.02, (f) 0.005. µ - 1 0 ω ( 0 t, confirming the weakly interacting Bose gas picture for -0.1 -π 0 π -π 0 π the superfluid fraction of particles/holes doped into the MI k k [17]. The hole bandwidth in both one and two dimensions is about zt. For the BHM model, the ratio of the bandwidths FIG. 3: (Color online) Phase diagram for detuning (a) ∆/g = −2, (b) ∆/g = 2, and spectra for (c) ∆/g = 2, t/g = 0.7, and (d) is Wp/Wh =2 for the MI with n =1 due to the effective par- ∆/g = −2, t/g = 0.0.018 in 1D at T = 0. ticle (hole) hopping tp = (n + 1)t (th = nt) [14, 19]. For the JCM, the matrix elements for hopping of one particle or hole in the MI have to be evaluated using the dressed states ψn , | i and we find Wp/Wh = 3/2+ √2 2.91 in one dimension, found for ∆ = 0 does not hold. The excitonic MI with ≈ in good agreement with Figs. 2(a)–(c). The bandwidth ratio is n =1 is approximately given by i , 0 i, whereas we have Q |↑ i fairly independent of respectively because the inter- i , 1 i for the photonic case. Finally, the incoherent fea- t/U t/g Q |↓ i action energy of a single particle or hole is the same at every tures in Fig. 3(c) may originate from the finite cluster size Lc, site. In the BHM, emergent particle-hole symmetry leads to and should be addressed using other exact methods. ∗ mp mh 0 for t t [20]. This behavior is also seen Solid state quantum devices will operate at finite tempera- ≈ → → in the JCM (mp/mh 1.3 for t =0.16 in 1D), but the VCA ture, much higher than in CAOL (estimated at nK). This leads does not permit a detailed≈ analysis near t∗. to the important question of the stability of the MI state at In contrast to CAOL, the parameters ǫ and ω0 and hence T > 0, which will ultimately determine their technological the detuning ∆ = ǫ ω0 can easily be changed experimen- usefulness. Strictly speaking, there is no true MI at T > 0 tally. This permits to− tune the character of the bosonic parti- due to thermal fluctuations. However, there exist regions cles. The states ψn are excitonlike for ω0 ǫ, photonlike where fluctuations are small enough for the system to be- | i ≫ for ω0 ǫ, and polaritonlike for ω0 ǫ. The 1D phase have like a MI for experimental purposes [21]. We determine diagrams≪ for ∆/g = 2 are shown in Fig.≈ 3(a),(b). The exci- the region of existence of the MI using the stringent criterion ± −4 tonic system (∆/g =2) with small photon-mediated hopping ∆n = n nint 10 , corresponding to a “worst case sce- | − |≤ exhibits a large n = 1 Mott lobe, whereas the latter is very nario” since experimentally the MI will survive as long as the small in the photonic system with small exciton-mediated in- Mott gap g is large compared to thermal fluctuations. ∼ teraction. Reentrant behavior due to quantum fluctuations is Results for ∆=0 are shown for one and two dimensions seen for the excitonic case, but is absent in the photonic case in Fig. 4. Finite-size effects diminish quickly with increasing for large enough Lc. temperature (see Fig. 4(a)). The size of the n = 1 Mott lobe Single-particle spectra in one dimension near the lobe tips is significantly reduced with increasing temperature, and ∆n are shown in Fig. 3(c),(d). The excitonic system shows a very exceeds 10−4 at T ∗/g 0.03 in both one and two dimen- ≈ large ratio of particle and hole bandwidths (Wp/Wh 7), sions. This value is consistent with the onset of deviations of ≈ whereas Wp/Wh 2.1 (very close to the BHM, since ψ1 n from 1 in the exact atomic-limit results, and agrees with the , 1 ) for the photonic≈ case. These values result from the| dif-i≈ few-cavity results of [7]. Lobes with n> 1 disappear at much |↓ i ferent admixture of the states , 0 and , 1 to ψ1 depend- lower temperatures (not shown). In contrast to T = 0, Mott |↑ i |↓ i | i ing on ∆. In particular, the approximate relation Wh zt lobes with n and n 1 do not touch at t/g = 0, but are sur- ≈ ± 4

In summary, motivated by theoretical predictions and on- (a) 1D, ∆ = 0 T/g = 0.00 (L = 6) 0.5 c going experimental advances, we have studied polariton Mott T/g = 0.01 (Lc = 2) T/g = 0.01 (Lc = 4) phases by means of a versatile quantum many-body approach. 0.4 T/g = 0.02 (L = 2) c Phase diagrams, single-particle spectra and finite tempera-

g T/g = 0.03 (L = 2) 0.3 c ture effects have been related to the known features of Bose- µ / Hubbard models and possible experiments, and the novel 0.2 physics emerging from detuning has been explored. 0.1 This work was supported by the FWF Schr¨odinger Fellow- 0 0 0.04 0.08 0.12 0.16 0.2 ship No. J2583, the DFG research unit FOR538, and by NSF T/g = 0 Award No. DMR-0502047. (b) 2D, ∆ = 0 0.5 T/g = 0.005 (L = 1 x 1) T/g = 0.01 c T/g = 0.02 0.4 T/g = 0.03 g 0.3 µ / [1] A. D. Greentree, C. Tahan, J. H. Cole, and L. C. L. Hollenberg, 0.2 Nat. Phys. 2, 856 (2006). [2] M. J. Hartmann, F. G. S. L. Brand˜ao, and M. B. Plenio, Nat. 0.1 Phys. 2, 849 (2006). [3] D. G. Angelakis, M. F. Santos, and S. Bose, Phys. Rev. A 76, 0 0 0.01 0.02 0.03 0.04 0.05 0.06 031805(R) (2007). t / g [4] M. J. Hartmann and M. B. Plenio, Phys. Rev. Lett. 99, 103601 (2007); Phys. Rev. Lett. 100, 070602 (2008). FIG. 4: (Color online) Mott-like regions for different . T [5] N. Na, S. Utsunomiya, L. Tian, and Y. Yamamoto, Phys. Rev. A 77, 031803(R) (2008). [6] S.-C. Lei and R.-K. Lee, Phys. Rev. A 77, 033827 (2008). rounded by the normal fluid (NF) phase for sufficiently small [7] M. I. Makin, J. H. Cole, C. Tahan, L. C. L. Hollenberg, and t/g, and a transition NF-SF occurs at larger t/g. However, A. D. Greentree, Phys. Rev. A (accepted). while the VCA can distinguish between the MI and the NF, [8] D. Rossini and R. Fazio, Phys. Rev. Lett. 99, 186401 (2007). the phase boundary to the SF at T > 0 cannot be determined [9] M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, A. Sen, accurately. The gap in A(k,ω) increases with T [21]. Parti- and U. Sen, Adv. Phys. 56, 243 (2007). [10] D. E. Chang, V. Gritsev, G. Morigi, V. Vuletic, M. D. Lukin, cle/hole bands are still well defined at T > 0, and spectra (not and E. A. Demler, arXiv:0712.1817 (unpublished). shown) look very similar to T = 0. Excitation spectra of the [11] J. Koch, T. M. Yu, J. Gambetta, A. A. Houck, D. I. Schus- NF and SF deserve detailed future studies. ter, J. Majer, A. Blais, M. H. Devoret, S. M. Girvin, and R. J. The quasi-MI lobes in Fig. 4 are dominated by the point Schoelkopf, Phys. Rev. A 76, 042319 (2007). (t = 0,µ 0.3) where Eg is maximal and by the T = 0 [12] C. W. Lai, N. Y. Kim, S. Utsunomiya, G. Roumpos, H. Deng, lobe tip (µ∗≈,t∗). The MI ise destroyed by thermal excitation of M. D. Fraser, T. Byrnes, P. Recher, N. Kumada, T. Fujisawa, particles or holes, and therefore survives longest near µ =0.3 and Y. Yamamoto, Nature 450, 529 (2007). [13] M. Potthoff, M. Aichhorn, and C. Dahnken, Phys. Rev. Lett. 91, where the energy cost is largest. Besides, T ∗ is determined 206402 (2003); M. Potthoff, Eur. Phys. J B 36, 335 (2003). by the size of the lobe at T = 0, t = 0. Hence, the system [14] W. Koller and N. Dupuis, J. Phys.: Condens. Matter 18, 9525 at T > 0 is dominated by atomic-limit physics (yielding the (2006). same T ∗ in 1D and 2D) and the T = 0 fixed point (quantum [15] M. Aichhorn, E. Arrigoni, M. Potthoff, and W. Hanke, Phys. critical point) of the MI-SF transition [17]. For the 1D BHM Rev. B 74, 235117 (2006). [17], the VCA yields T ∗/U 0.059, lower than T ∗/U =0.2 [16] T. D. K¨uhner and H. Monien, Phys. Rev. B 58, R14741 (1998). found in [22] using a different≈ criterion [25]. The ratio of the [17] M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. Fisher, critical temperatures for the 1D JCM and 1D BHM is 0.56, Phys. Rev. B 40, 546 (1989). [18] N. Elstner and H. Monien, Phys. Rev. B 59, 12184 (1999). close to the ratio of the t =0, T =0 Mott gaps (0.59). [19] D. van Oosten, P. van der Straten, and H. T. C. Stoof, Phys. Rev. Solid state cavity-QED systems offer the possibility of large A 63, 053601 (2001). g. Taking a conservative estimate g = 1010 Hz [1] gives for [20] B. Capogrosso-Sansone, N. V. Prokof’ev, and B. V. Svistunov, the effective temperature of the quantum system T ∗ & 14 Phys. Rev. B 75, 134302 (2007). mK. Even taken as the actual temperature, this is well ac- [21] D. van Oosten, P. van der Straten, and H. T. C. Stoof, Phys. Rev. cessible experimentally. The hopping rate 1/t has to be fast A 67, 033606 (2003). [22] F. Gerbier, Phys. Rev. Lett. 99, 120405 (2007). enough to permit equilibration before photon loss and de- [23] As in [1]; ω0 is not an absolute frequency, so that this choice phasing set in. Indeed, taking t/g = 0.01 (inside the Mott does not affect results. −1 lobe in Fig. 4(b)), we obtain a realistic t & 0.14 ps. Fi- [24] For t/g = 0 and density n, a particle/hole corresponds to one nally, T ∗ for the n = 1 MI is enhanced by detuning ∆ > 0 site having occupation n + 1 (n − 1). (Fig. 3(a)), which additionally increases t∗ and µ∗ , estimated [25] There are no well-defined plateaus in the density even below ∗ as (µ∗,t∗) (0.08, 0.02) meV in 2D for T =∆=0| | . T /U = 0.2 in [22]. ≈