{'fraction_non_alphanumeric': 0.037069587871559055, 'fraction_numerical': 0.008882890235098402, 'mean_word_length': 4.462962962962963, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Quantum noise in a transversely pumped cavity Bose--Hubbard model

We investigate the quantum measurement noise effects on the dynamics of an atomic Bose lattice gas inside an optical resonator. We describe the dynamics by means of a hybrid model consisting of a Bose--Hubbard Hamiltonian for the atoms and a Heisenberg--Langevin equation for the lossy cavity field mode. We assume that the atoms are prepared initially in the ground state of the lattice Hamiltonian and then start to interact with the cavity mode. We show that the cavity field fluctuations originating from the dissipative outcoupling of photons from the resonator lead to vastly different effects in the different possible ground state phases, i.e., the superfluid, the supersolid, the Mott- and the charge-density-wave phases. In the former two phases with the presence of a superfluid wavefunction, the quantum measurement noise appears as a driving term leading to excess noise depletion of the ground state. The time scale for the system to leave the ground scale is determined analytically. For the latter two incompressible phases, the quantum noise results in the fluctuation of the chemical potential. We derive an analytical expression for the corresponding broadening of the quasiparticle resonances.


I. INTRODUCTION
Cavity Quantum Electrodynamics (CQED) is devoted to studying the interaction of the electromagnetic field with atoms under the best possible control of circumstances. The radiation field is tailored by resonators in order to select spatially and spectrally one or a few relevant field modes. The atomic positions are controlled as well; either by well-defined trajectories across the resonator, or by trapping the atoms in well-defined positions [1,2]. The control over the components allows for accessing the coherent quantum dynamics of the coupled atom-field system. The first milestone has been the demonstration of the vacuum Rabi splitting [3], which is the benchmark of strong coupling between the induced electric dipole of the atom and the cavity mode. The possibility of observing coherent processes, such as the Rabi oscillation, is limited by the spontaneous photon scattering into modes other than the cavity mode. The larger the electric dipole coupling strength g with respect to the atomic spontaneous emission rate γ and the cavity mode linewidth κ, the shorter is the dynamical timescale needed to resolve the Jaynes-Cummings spectrum. Various nonlinear quantum effects, e.g., the photon blockade [4], or the two-photon gateway [5] have been observed in the strong coupling regime of cavity QED.
With the use of ultracold atoms in CQED experiments, the magnitude of the collective coupling of an atomic ensemble to the cavity mode can be significantly enhanced to higher values than those characterising the loss rates. This allows for applying large detuning between all laser excitations and the internal atomic resonances, and thereby the atomic scattering loss processes can be significantly suppressed. As a consequence, the atom-cavity dynamics can be controlled on much longer time scales. This opportunity opened up the way to a new regime of cavity QED experiments, where the spatial motion of an atom cloud couples coherently to the dynamics of the cavity field mode [6]. The effective Hamiltonian describing the system includes characteristic frequencies well below the single-atom coupling strength g, e.g., the so-called recoil frequency ω R = k 2 /(2m), where k is the cavity mode wavenumber and m is the atomic mass. These experiments revealed, for example, a Dicke-type superradiant phase transition of an atomic superfluid in the cavity [7][8][9], and demonstrated optomechanical strong coupling between vibration and field intensity [10]. The limitation of the coherent dynamics on the time scales longer than the inverse of the recoil frequency originates from the spontaneous photon scattering again, however, (i) the spontaneous emission from atoms into free-space modes is strongly suppressed, (ii) the photon loss from the cavity mode into external modes is weakly coupled into the dynamics in a rather indirect way.
A very new generation of CQED experiments [11,12] introduces a time scale which is even longer than the inverse recoil frequency. Ultracold atoms trapped in optical lattices sustained by the cavity can tunnel between adjacent sites. The many-body quantum state of the atoms in the lattice sites becomes relevant to the dynamics of the cavity field amplitude [13][14][15]. The tunneling time depends on the depth of the trapping potential which is a novel control parameter in the cavity QED system. The effective Hamiltonian corresponds to the family of bosonic Hubbard-type lattice models extended to include the cavity field mode [16]. In these Hamiltonian systems exotic new phases of lattice bosons appear due to the cavity mediated global-range interactions [17][18][19][20][21][22][23][24][25][26][27][28]. In fermionic lattices, cavity-induced topologically nontrivial [29][30][31] states can be generated. All these research topics stimulate nowadays a significant theoretical and experimental interest in atomic lattice gases integrated in CQED systems.
Despite both the pronounced theoretical and experimental interest, very little effort has been made to study the limitations of a Hamiltonian approach. In a very recent paper Chiacchio and Nunnenkamp studied the time evolution of the density matrix with integrating the mas-FIG. 1. Illustration of the coupled cavity Bose-Hubbard model setup. An atomic cloud is loaded into a square optical lattice, which is inside a single-mode high-Q Fabry-Pérot resonator. The period of the cavity mode is approximately equal to that of the optical lattice. The cavity is pumped by the side through the light scattered by the atoms from the laser drive. The system is open, photons leak out from the cavity, resulting in heating and decoherence.
ter equation [32]. They found that in the bad cavity limit, which is close to the experimental situation of Ref. [11], the steady state is an infinite temperature state. In the experiments, however, they find a more or less coherent evolution during the time of the measurements.
In this paper we investigate the time limitation on the coherence in cavity Bose-Hubbard models. We consider the cavity photon loss as the dominant dissipative process. Since the outcoupled photons can be directed to a photo-detector, the fluctuations associated with the loss process can be equally well conceived as the measurement-induced backaction on the quantum system. We derive simple analytical formulae in terms of experimentally measurable quantities to quantify the time limitation of the Hamiltonian description. In Section II, we introduce one example of a cavity Bose-Hubbard model and derive the Heisenberg-Langevin equations that take into account the fluctuations of the cavity mode. In Sec. III, we derive an effective model for the lattice bosons by adiabatically eliminating the cavity field. In Sec. IV, we calculate and compare the effects of dissipation noise in the (i) superfluid-type and (ii) in the Mott-type phases. Finally, we summarise the results in Sec. V.

II. CAVITY BOSE-HUBBARD MODEL
An atomic degenerate Bose gas trapped in an optical lattice inside a high-finesse optical cavity gives rise to a system with competing short-and long-range interactions. The lattice dynamics of the ultracold atoms including on-site collisional interactions and tunneling between adjacent sites corresponds to the usual Bose-Hubbardmodel [33,34]. The optical lattice depth set by the power of its driving laser controls the strength of these short-range effects. On top of the bosonic lattice model, there is a long-range atom-atom interaction mediated by the cavity field in photon scattering processes. In order to be specific, we consider the geometry of the experimental setup in Ref. [11], illustrated in Fig 1. A sample of bosonic atoms is placed inside a two-dimensional optical lattice, which overlaps with a mode of an optical cavity. A single mode is selected of which the wavelength is very close to that of the optical lattice. The atoms are illuminated from the side by means of a far-detuned laser source which is close to resonance with the cavity mode. The atoms scatter photons coherently between the laser and the commonly coupled cavity mode [6]. The collective coupling of the atoms to the cavity mode leads to an indirect and infinite-range interaction between the atoms, which can be controlled independently from other parameters of the lattice dynamics by means of a small detuning of the drive frequency from the cavity resonance. The effective Hamiltonian of the system reads, see Ref. [11], where The first term is the standard Bose-Hubbard model describing the atomic dynamics inside the optical lattice. The atomic annihilation and creation operators at site i are b i and b † i , respectively. The tunneling is characterised by J, and the strength of the on-site interaction is denoted by U s . The second term represents the oscillator of the cavity field, with a being the photon annihilation operator, and ∆ C = ω L −ω C the cavity detuning, i.e., the frequency difference of the pumping laser and the empty cavity. The interaction term H I is a driving of the cavity mode which depends on the atom density n j = b † j b j , and as a key element, on a phase factor (−1) j . The phase comes from that the atoms scatter the laser light into the cavity with an amplitude depending on the sign of the cavity mode function at the given site. This sign alternates between adjacent sites of the optical lattice, discriminating "even" and "odd" sites in a checker-board pattern. The strength of the interaction is modelled by the parameter η, which is proportional to the pumping amplitude of the laser and inversely proportional to cavity detuning ∆ C .
Because of the very large detuning of the driving laser with respect to all atomic resonances, the atoms scatter photons only in a coherent manner and the spontaneous emission can be neglected. However, one must take into account that the photons leak out from the cavity through the mirrors and couple to the free-space modes. The loss process can be treated within the Markov approximation with an intensity loss rate of 2κ. The Heisenberg equations of motion of the operators are complemented by a damping term and a corresponding Langevin-type quantum noise term. The Heisenberg-Langevin equations are where ∂ t denotes the time derivative, ∆n = j (−1) j n j is the imbalance between the particle numbers on the even and odd sites. It can be seen from the form of this term that the atoms act collectively on the cavity mode amplitude. In Eq. (2b), the sum for ε is over the 4 neighbouring sites of the site j. The noise term in Eq. (2a) is a delta-correlated white noise with zero mean value: ξ(t) = 0. Furthermore, let us assume that the temperature is very low compared to optical frequencies, that is, normal-ordered correlations are zero. In particular, ξ † (t )ξ(t) = 0. All other correlations are determined by the bosonic commutation relations, in second order. The goal of this paper is to reveal the dynamical consequences of this noise term during the initial short evolution time.

III. ADIABATIC DYNAMICS OF THE ATOMS
Since the time scale of the cavity field relaxation is the far shortest, i.e. κ, ∆ C J, U s , the cavity mode can be slaved to the slow atomic lattice dynamics. Performing adiabatic elimination of the cavity dynamics, it results in a renormalisation of the parameters describing the atomic motion in the dynamical lattice. Furthermore, damping and decoherence are introduced into the atomic time evolution. These latter effects are modeled also by a Langevin-like noise, now as part of the atomic motion, which, in general, is not a white noise.
Let us start by integrating out the fast cavity field equations of motion (2a), with keeping the slow atomic operators constant, where the first term is the adiabatic steady state of the photon field, and the noise is which is a white noise filtered through the cavity mode. The commutation relation can be derived from that of the original Eq. (3), Combining Eqs. (4) with Eq. (2b), we arrive at the adiabatic dynamics where the last noise term includes the product of atomic operators and the noise R(t) = Σ(t) + Σ † (t), being a selfadjoint operator. The mean value of R(t) is zero, and the second-order correlations are evaluated similarly to Eq. (5) with the following result, When the photon decay κ −1 is much shorter than the other timescales of the problem, we can approximate Eqs. (5) and (7) by a delta-correlated noise, with This approximation corresponds to the broad bandwidth reservoir assumption used for Markovian decay. We note that the appearance of fluctuations due to the lossy cavity mode is accompanied by non-adiabatic drift terms, the so-called cavity cooling or heating, depending on the sign of the detuning ∆ C Ref. [35]. However, this is beyond the adiabatic approximation and can be safely neglected in the limit κ ω R .

IV. MANIFESTATION OF THE MULTIPLICATIVE NOISE IN THE DIFFERENT PHASES
In this section we analyze the effects of the noise created by the leakage of cavity photons on the atomic dynamics. We assume that this effect is small and that the noise can be considered as an additional effect on top of the quantum fluctuations arising from the interactions within the many-body system.
The study has to be done separately for the distinct cases corresponding to the different possible thermodynamic phases [18][19][20][21][22][23][24]. In fact, the model allows for four different ground state phases depending on the preservation or breaking of two different symmetries. Firstly, the global U(1) symmetry associated with particle number conservation, the breaking of which corresponds to the presence of off-diagonal long-range order. Secondly, the Z 2 symmetry of even and odd lattice sites, which the atom-cavity interaction Eq. (1d) can break spontaneously. When the Z 2 symmetry is not broken, all averages on the even and odd lattice sites are equal. When the Z 2 symmetry is broken, the expectation values of operators on the even and on the odd sublattice can be different. The four phases are: the superfluid phase with broken U(1) and unbroken Z 2 symmetry, the supersolid phase with both symmetries broken, the Mott phase with both phases unbroken, and finally, the charge-density wave (CDW) phase with a broken Z 2 and unbroken U(1) symmetry.
The broken U(1) symmetry phases emerge when the on-site atom-atom interaction is small compared to tunneling, U s J. In this weakly-interacting limit the equations of motion of the atomic field operators can be truncated at second order in the usual perturbative way, and can be eventually dealt with Bogoliubov transformation [36]. In the other limit, when U s J, one goes to the strongly interacting phases, where U(1) is unbroken. Here the Mott phase and the CDW phase are described after a canonical transformation to the low energy degrees of freedom [37,38].

A. Superfluid and supersolid phases
In the weakly interacting limit, where U s J, Bose-Einstein condensation occurs for sufficiently low temperatures. The wavefunction of the condensate is proportional to the expectation value of the atomic ladder operators, where the coefficient is the square root of the average atomic density of the condensate. Close to zero temperature almost all the atoms are condensed, and the condensate density can be approximated by the total density. It is convenient to separate the mean values from the fluctuations, where n = N/M is the total density of atoms, i.e. the total number of atoms N divided by the total number of sites M . Furthermore, µ plays the role of the chemical potential. When inserting Eq. (10) into Eq. (6) and collecting the terms proportional to √ n, we arrive at the Gross-Pitaevskii equation (GPE), We introduced the dimensionless parameter Y 2 = −2 ∆ C η 2 N/(∆ 2 C + κ 2 ) characterising the power of the pumping laser and the scaled on-site interaction strength g = nU s . Furthermore, we introduced the supersolid order parameter Θ, which is the difference between the even and odd site condensate densities, which are normalised to unity in agreement with Eq. (10), Equations (11), (12), (13) form a closed set for the chemical potential µ, the supersolid order parameter Θ and the condensate wavefunction β j . They have a remarkably simple solution both in the superfluid and in the supersolid phases. In both phases, where z is the coordination number of the lattice. For the square lattice, z = 4. Furthermore, by introducing β e = √ 1 + Θ for the even sites and β o = √ 1 − Θ for the odd sites, the solution becomes β j = β e , for j even, β o , for j odd.
Finally, the following equation has to be fulfilled: from which either Θ = 0, or µ = 2g − Y 2 . The superfluid phase is characterised by Θ = 0, i.e., with a homogeneous condensate β e = β o = 1. In this phase, from Eq. (14), we get µ = g − zJ. In the supersolid phase µ = 2g − Y 2 , and from Eq. (14), we arrive to The critical pumping power separating the two phases is at The chemical potential and the order parameters are plotted in Fig. 2.
Using again the substitution Eq. (10) in Eq. (6) but now keeping the terms linear in the atomic operators and the noise, we get the equation of motion for the fluctuations, We restrict ourselves to the low energy excitations and neglect the wavenumber dependence of the fluctuations. That is, we keep only the two relevant modes of the symmetry breaking: the fluctuations of the ladder operators of the even and odd sites. The 2 coupled equations become These Bogoliubov equations are constant-coefficient inhomogeneous linear differential equations with a Langevin-type noise added as driving. Another 2 equations describing the time evolution of δb † e and δb † o have to be added to get a closed set of equations for the fluctuations. These latter two equations are obtained by taking the hermitian conjugates of Eqs. (19). First, let us gather the fluctuation operators to the formal vector w = (δb e , δb o , δb † e , δb † o ). Then, the four Bogoliubov equations become where the 4 × 4 coefficient matrix M 4×4 and the noise vector Ξ can be read off from Eqs. (19). The solution is obtained by finding the normal modes, i.e., diagonalising the coefficient matrix. The eigenvalues of M 4×4 come in ± pairs, and 2 of the 4 eigenvalues are identically zero in the entire Bose condensed phase [36]. The normal mode with zero eigenvalue, which we call the zero mode, corresponds to phase fluctuations of the condensate. The corresponding eigenvector is proportional to the condensate wavefunction, and therefore, the zero mode does not describe any of the fluctuations orthogonal to the con-densate. Furthermore, its operator is anti-hermitian and is decoupled from the dynamics of the other normal mode with eigenvalues ±λ. This latter excitation describes density waves with a period of 2 lattice sites. At the transition point between the superfluid and the supersolid phases, this excitation also becomes soft, therefore, we call it as the soft mode. In other words, we introduce δb z and δb s , corresponding to the zero and soft modes, respectively, as (δb e , δb o ) = γ δb z + ε δb s , where γ = (β e , β o )/ √ 2, and ε = (β o , −β e )/ √ 2, orthogonal vectors. These two vectors span the 2-dimensional space, from which, only the direction orthogonal to the condensate, i.e., the soft mode is relevant for us. The soft mode decouples from the purely anti-hermitian zero mode, and its Bogoliubov equation is obtained directly from Eq. (20) simply by taking the appropriate submatrix of M 4×4 and subvector of Ξ, with the reduced vector v = (δb s , δb † s ) T and noise Ξ (t) = (ζ, −ζ) T . The 2 × 2 coefficient matrix reads while the effective noise operator is Note, that this noise operator inherits the properties of R(t), hence it is Hermitian and its correlation is determined by Eq. (8).
The noise leads to a diffusion-like process, that depletes the superfluid or supersolid ground states by heating the condensate atoms into the orthogonal fluctuation mode. In the following we calculate the rate by which the atoms leave the Bose condensed state. The matrix M is diagonalised by its real right and left eigenvectors, M r (k) = ω k r (k) and (l (k) ) T M = ω k (l (k) ) T . Their scalar product is conveniently normalised, (l (k) , r (k) ) = 1. Multiplying Eq. (22) from the left with (l (k) ) T , we obtain the equation of motion of the normal modes ρ k = (l (k) , v) that reads where Q k (t) = (l (k) , Ξ (t)) is the projection of the noise vector to each mode. The normal mode frequencies (eigenvalues of the coefficient matrix M ) are ω 1,2 = ± −1 λ, with The eigenvalues are plotted in Fig. 3. The corresponding normal modes form a Hermitian adjoint pair ρ 2 = ρ † 1 . Their correlation function can be directly calculated from the formal solution of Eq. (25), The fluctuation number is expressed with the normal modes The summation over the eigenvector components gives a factor of −1. According to Eq. (27), the system obeys an exponential relaxation law towards the steady state when the eigenvalue λ has an imaginary part, such as in Ref. [39]. This finding is similar to what was found in Ref. [32] in the bad-cavity limit and for the case of vanishing on-site interactions. For short times, t λ −1 the incoherent population, Eq. (28) builds up linearly in time, hence the rate characterising the diffusion is The right-hand side of Eq. (29) is written as a product of two factors. The first one sets the dimension and the order of magnitude of the diffusion, while the second factor is dimensionless and is on the order of unity. This latter quantity is plotted in Fig 4. In the superfluid phase it is proportional to Y 2 , while in the supersolid phase it goes to zero according to 1/Y 2 . The first factor has the dimension of s −1 , and is combined from 3 different frequency scales: the photon escape rate (κ), the magnitude of the cavity detuning (|∆ C |) and the optical-lattice tunneling rate (zJ/ ). In the experiments, the first two of the three frequencies are usually chosen to be close to each other, while zJ is much smaller then the other two. Furthermore, the number of lattice sites illuminated by the cavity mode is on the order of a thousand (M ∼ 10 3 ). Therefore, we estimate 2κzJ/(M |∆ C |) ∼ 10 −2 zJ -10 −3 zJ. Clearly, this magnitude of the diffusion process allows for considering the tunneling dynamics to be coherent on a long enough time scale. Notably, some of the technical noise sources can be more restrictive in this respect. However, one must be aware of that, with the cavity detuning ∆ C going to zero, the photon outcoupling process can generate a large number of photons within the cavity. These results are in complete analogy to the model without the additional optical lattice [40]; the noise leads to the departure of the system from its quantum ground state, and for long times it relaxes exponentially to a steady state with an excess noise depletion of the condensate [41] in the cavity-cooling regime, while it relaxes exponentially to the infinite-temperature state in the cavity-heating case.  29), by which the atoms leave the superfluid and supersolid ground states. More precisely, we plot the dimensionless quantity Y 2 zJ/(g − µ) 2 as a function Y 2 . The diffusion rate exhibits a cusp at the transition point between the superfluid and supersolid phases.

B. Mott and CDW phases
When the strength of the on-site repulsion is much larger than the amplitude of the hopping (U J), multiple occupations of a site become energetically penalised. For commensurate fillings, i.e., when the number of atoms is integer times the number of lattice sites, each lattice site becomes populated with exactly the same number of atoms. In the same time, particle number fluctuations become suppressed at each individual site. Depending on the strength of the transverse laser driving, the system can be in a Mott state or in a CDW state. For weak pumping, the Mott state is realised, where each site has the same number of particles. In contrast, when the pumping strength is large, the CDW state is the ground state, where all even sites have the same occupation and all odd sites have the same occupation, but these two are different. Both in the Mott-and CDW phases, the low energy excitations are no longer the Bogoliubov quasiparticles but rather particles and holes of a strongly correlated system. These quasiparticle excitations have to be introduced separately for each phase. To be specific, we are going to study the excitations over the Mott state with one particle per site and over the CDW state with one particle per two sites. The generalisation of the theory to other Mott-and CDW states is straightforward and the conclusions do not change qualitatively. For a more transparent presentation, we also completely neglect the effect of the kinetic energy, which can be reintroduced with the help of perturbation theory, as was done for the Mott-superfluid transition in Refs. [37,38].
Without tunneling, the Hamiltonian (1) becomes In the one particle per site Mott lobe, we follow the route of Ref. [37] and truncate the single site Hilbert space for the three lowest occupied states, |0 j , |1 j , |2 j . We introduce three Schwinger boson operators creating these states, The state |vac is a fictitious vacuum state, not present in our original Hilbert space. In order to exclude the unphysical states, these Schwinger bosons must share a single excitation at each site, therefore we have the constraint In the present case, the ground state is a tensor product state of the singly occupied sites, That is, in the ground state, the occupation of the t 1 boson is 1 at each site, while the other bosons have zero occupations. Low energy excited states have these other bosonic excitations mixed in with some amplitude much smaller than unity, while the t 1,j bosons having still an occupation very close to unity. Thus, by using Eq. (32), the Hamiltonian can be written as Equation (34) is quadratic in the excitations, i.e., both the hole-type excitations created by t † 0,j and the particlelike excitations of t † 2,j . In the Mott phase, the effect of the photon field is through renormalising the (chemical) potential. Integrating out the photon field, we still have Eqs. (4), but now with To lowest order in the single particle per site Mott phase, we can approximate ∆n ≈ j (−1) j = 0. That is, in the Mott phase there is no classical part of the cavity field a ss = 0. Therefore, the chemical potential is renormalised only by a noise, Therefore, in the Mott phase, for zero tunneling, the effect of cavity decay is to make the chemical potential noisy. As the Mott phase is gapped, there is no effect of the noise on the ground state while we stay away from the phase boundary. The fluctuation of the chemical potential affects, on the other hand, the particle-and hole-like excitations of the Mott phase: the corresponding quasiparticle resonances get broadened by an amount of When a small enough tunneling is introduced, the situation can be more complicated, as the noise can dress also the quasiparticles. This effect must be small, as it has to vanish for zero tunneling. In Ref. [32] the authors found anomalous and normal diffusion towards the steady state similar to the standard optical lattice Mott insulator [42]. The present analysis was performed for the first Mott sate, i.e., the one with n = 1 particles per site. This analysis can simply be generalised to any of the Mott states with n particles per site by keeping the most relevant three states, namely |n − 1 j , |n j , and |n + 1 j at each site. Equations through (33)- (36) has to be modified, although in a straightforward manner. Nevertheless, the final result for the noisy chemical potential, Eq. (37), holds for each Mott state.
In order to generalise the previous analysis further to the CDW phase, we need to introduce two sublattices in a checkerboard setting. One sublattice has n e particles on each site, while on the other sublattice each site has n o atoms. We shall refer to these sublattices as the even and odd sublattice, respectively. A unit cell now contains two neighbouring sites, each from a different sublattice. The Hamiltonian reads The sum goes over the unit cells indexed by m. The number of unit cells is half the number of sites. For simplicity we work with the n e = 1 and n o = 0 CDW phase. As the occupation number of each site on the even sublattice is unity, we use the same 3 states and the same 3 Schwinger bosons as in Eqs. (31). On the odd sublattice, we have no atoms at all in the ground state, and we use only two Schwinger bosons corresponding to the empty and to the singly occupied sites. That is, if site j is the even site in unit cell m, |1 m,e = t † 1,m |vac , |0 m,e = t † 0,m |vac , |2 m,e = t † 2,m |vac .
Where the photon field is adiabatically eliminated in order to reach the final result. In this phase, ∆n = N , since the odd sites are empty, and ∆µ is negative (note, that ∆ C is negative). In fact, this CDW phase is stable while the energy of the s 1,m particle excitations on the odd sites is positive, i.e., while |∆µ| > µ. Below that, the system relaxes to the n = 1 Mott phase. Notice, that the noise term randomly pulls ∆µ, and the boundary of the phase gets smoothed out.

V. DISCUSSION AND SUMMARY
In this paper, we studied the dynamics of a lattice Bose gas, which apart from the external optical lattice, is immersed in the optical field of a single-mode high-Q Fabry-Pérot interferometer. The wavelength of the resonator is tuned close to that of the optical lattice in order to allow the lattice gas to form a grating. This atomic grating can scatter light into the cavity mode either with constructive interference, thus allowing a classical cavity field to build up, or scatter the light destructively thereby destroying the cavity field. There is a phase transition in between these two regimes separating two phases with different lattice periodicity.
We derived analytical expressions for the superfluidsupersolid phase transition in the weak-coupling limit, and an analytical formula for the excitation energy of the critical density wave. In the experimental realisation of the cavity Bose-Hubbard model, the optical resonator is an inherently open quantum system. Photons from the laser pump are scattered into the cavity, and finally leak out through the cavity mirrors. The corresponding dissipation process leads to a quantum noise, that can excite the system out from its ground states. It is a fundamental question what time limitation this intrinsic quantum process imposes on the validity of the ground state description, which is substantial to all the studies relying on a Hamiltonian description of the cavity Bose-Hubbard model. We calculated analytically the effects of photon loss dissipation in the superfluid and Mott phases. We showed that the dissipative quantum noise has substantially different effects in the two phases. In the superfluid phase, the noise heats the system out from its ground state with a time scale 2κzJ/(M |∆ C |). This time scale is a slow one except for very small detunings.
In the Mott-type phase the gap protects the population in the ground state against the photon field fluctuation noise and allows only a much slower approach towards equilibrium. However, the noise blurs also the phase boundaries and affects the excitation spectrum by the quasiparticle resonance broadening with an amount of 2 κ η 2 /(κ 2 + ∆ 2 C ).