11 Quantum Criticality in the Biased Dicke Model
Hanjie Zhu( 朱汉杰 ) ,Guofeng Zhang( 张国锋 ) ∗ ,Heng Fan( 范桁 ) Key Laboratory of Micro-Nano Measurement-Manipulation and Physics (Ministry of Education), School of Physics and Nuclear Energy Engineering;
State Key Laboratory of Software Development Environment,
Beihang University, Xueyuan Road No. 37, Beijing 100191, China State Key Laboratory of Low-Dimensional Quantum Physics, Tsinghua University, Beijing 100084, China Key Laboratory of Quantum Information, University of Science and Technology of China, Chinese Academy of Sciences, Hefei 230026, China Beijing National Laboratory for Condensed Matter Physics, and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
The biased Dicke model describes a system of biased two-level atoms coupled to a bosonic field, and is expected to produce new phenomena that are not present in the original Dicke model. In this paper, we study the critical properties of the biased Dicke model in the classical oscillator limits. For the finite-biased case in this limit , We present analytical results demonstrating that the excitation energy does not vanish for arbitrary coupling. This indicates that the second order phase transition is avoided in the biased Dicke model, which contrasts to the original Dicke model. We also analyze the squeezing and the entanglement in the ground state, and find that a finite bias will strongly modify their behaviors in the vicinity of the critical coupling point.
The Dicke model, which considers a system of a single-mode bosonic field coupled to N two-level atoms, plays a key role as a model illustrating the collective and coherent effects of many atoms in quantum optics . Numerous efforts have been paid to understand its properties over the past few decades, which have resulted in a wide variety of phenomena . One of these interesting phenomena is that the Dicke model undergoes an equilibrium phase transition in the classical limit as the coupling between the atoms and bosonic field reaches a specific value. This type of phase ∗ Correspondence and requests for materials should be addressed to G.Z.( [email protected]) transition is known as the superradiant phase transition (SPT) and has been intensely discussed . The early work on the SPT in the Dicke model mostly considers the classical spin (CS) limit, where the number of atoms N tends to infinity. In the classical oscillator (CO) limit, where the ratio of the atomic transition frequency to the bosonic field frequency approaches infinity, the Dicke model also experiences a SPT even at a finite N. This situation had been largely overlooked, and some efforts have been devoted to understand it
10, 12 . This is due to the fact that the Dicke model is mainly realized in the cavity quantum electrodynamics (QED) systems, where the coupling between the atoms and cavity field is weak compared to the atomic and cavity frequencies. Therefore, the achievement of the critical coupling requires using large atomic ensembles. Recently, the superconducting circuit-QED systems have achieved the ultrastrong coupling between the qubit and oscillator . This enables us to study the strong coupling effects including the SPT even for a single qubit. The other advantage of the circuit-QED systems is that the qubit parameters can be easily adjusted by the applied bias current, gate voltage, and microwave fields . However, this controllability produces an additional bias term in the Hamiltonian describing the qubit, and thus leads to the biased Dicke model in the circuit-QED systems. While most previous studies have focused on the zero-biased case, one can expect that this additional bias term will produce new phenomena that are not present in the original Dicke model. Several approaches have been developed to describe the behavior of this biased system . It has been shown in recent paper that the bias terms will create additional coupling to the environment and increase the fragility of non-classical states . The bias term will also smear out the SPT in the CS limit . The system that we consider is composed of (cid:1840) biased qubits coupled to a single-mode bosonic field. The biased Dicke model describing this system is given by (ℏ = 1) (cid:1834) = (cid:1834) (cid:3044) + (cid:1834) (cid:3033) + (cid:1834) (cid:3036)(cid:3041)(cid:3047) , ( 1 ) where (cid:1834) (cid:3044) = (cid:3533) (cid:3436)Ω2 (cid:2026) (cid:3053)(cid:3036) + (cid:2013)2 (cid:2026) (cid:3051)(cid:3036) (cid:3440) (cid:3015)(cid:3036)(cid:2880)(cid:2869) , (cid:1834) (cid:3033) = (cid:2033)(cid:1853) (cid:2878) (cid:1853), (cid:1834) (cid:3036)(cid:3041)(cid:3047) = (cid:2019) (cid:3533) (cid:2026) (cid:3051)(cid:3036)(cid:3015)(cid:3036)(cid:2880)(cid:2869) ((cid:1853) + (cid:1853) (cid:2878) ). ( 2 ) Here (cid:1853) and (cid:1853) (cid:2878) denote the annihilation and creation operators for the bosonic field with frequency (cid:2033) respectively, and {(cid:2026) (cid:3036)(cid:3038) ; (cid:1861) = (cid:1876), (cid:1877), (cid:1878)} are Pauli matrices for the k-th qubit with transition frequency Ω and bias (cid:2013) . The field couples to the atoms uniformly with the coupling strength (cid:2019) . This system is equivalent to a pseudospin of length (cid:1840)/2 coupled to a single mode field, and its Hamiltonian can be written as (cid:1834) = (cid:2033)(cid:1853) (cid:2878) (cid:1853) + Ω2 (cid:1836) (cid:3053) + (cid:2013)2 (cid:1836) (cid:3051) + (cid:2019)(cid:1836) (cid:3051) ((cid:1853) + (cid:1853) (cid:2878) ), ( 3 ) where the angular momentum operators (cid:1836) (cid:3036) = ∑ (cid:2026) (cid:3036)(cid:3038)(cid:3015)(cid:3038)(cid:2880)(cid:2869) , (cid:1861) = (cid:1876), (cid:1877), (cid:1878) . The biased Dicke model can be realized easily in the system of the flux qubits coupled to a quantum oscillator . In the basis of clockwise and anticlockwise qubit persistent currents, the Hamiltonian of the flux qubit can be written as (cid:1834) (cid:3044)(cid:3029) =((cid:2007)(cid:2026) (cid:3053) + (cid:2013)(cid:2026) (cid:3051) )/2 . Here (cid:2007) is the tunnel splitting and (cid:2013) = 2(cid:1835) (cid:3043) (Φ (cid:3032)(cid:3051) − 3Φ (cid:2868) /2) is the energy bias, where (cid:1835) (cid:3043) is the persistent current in the qubit, Φ (cid:3032)(cid:3051) is the external flux threading the qubit loop and Φ (cid:2868) is the flux quantum. The interaction between the qubit and the oscillator can be described by (cid:2019)(cid:2026) (cid:3051) ((cid:1853) + (cid:1853) (cid:2878) ) , where (cid:2019) = (cid:1839)(cid:1835) (cid:3043) (cid:1835) (cid:2868) is the coupling strength, (cid:1835) (cid:2868) = (cid:3493)(cid:2033)/(4(cid:1838)) is the measure for zero-point current fluctuations, (cid:1838) is the inductance of the wire. Then the system Hamiltonian reads (cid:1834) = (cid:2033)(cid:1853) (cid:2878) (cid:1853) + 12 ((cid:2007)(cid:2026) (cid:3053) + (cid:2013)(cid:2026) (cid:3051) ) + (cid:2019)(cid:2026) (cid:3051) ((cid:1853) + (cid:1853) (cid:2878) ), ( 4 ) and is equivalent to the single qubit case of the biased Dicke model. Mean-field (MF) theory results
First we discuss the properties of the system using the mean-field method. Similar to the non-biased spin-boson system , the mean-field ansatz for the ground-state wavefunction of the biased Dicke model is a product of a spin coherent state |(cid:2016)(cid:1767) and a boson coherent state |(cid:2009)(cid:1767) , which reads |(cid:2032) (cid:3014)(cid:3007) (cid:1767) = |(cid:2016)(cid:1767)⨂|(cid:2009)(cid:1767).
The values of the spin inversion and the mean photon number can be written as: 〈(cid:1836) (cid:3053) 〉 = −(cid:1840)cos(cid:2016), 〈(cid:1853) (cid:2878) (cid:1853)〉 = (cid:2009) (cid:2870) . ( 5 ) Then the energy functional is given as (cid:1831)((cid:2009), (cid:2016)) = (cid:2033)(cid:2009) (cid:2870) + 12 (cid:1840)(cid:2013)sin(cid:2016) − (cid:2007)2 (cid:1840)cos(cid:2016) + 2(cid:2009)(cid:2019)(cid:1840)sin(cid:2016). ( 6 )
By minimizing the energy with respect to α and θ respectively, we acquire the relations (cid:2009) = − (cid:2019)(cid:2033) (cid:1840)sin(cid:2016), ( 7 ) − 2(cid:2019) (cid:2870) (cid:2033) (cid:1840)sin(cid:2016)cos(cid:2016) + (cid:2013)2 cos(cid:2016) + (cid:2007)2 sin(cid:2016) = 0. ( 8 ) The second equation is identical to the semiclassical calculation in Ref. [11]. The Eq. (8) can be expressed as (cid:3436)(cid:2018) (cid:2870) − 1cosθ(cid:3440) sin(cid:2016) = (cid:2013) (cid:4593) . ( 9 ) Here we introduce the rescaled coupling constant (cid:2018) = 2(cid:2019)√(cid:1840) √(cid:2033)(cid:2007)⁄ and the rescaled bias (cid:2013) (cid:4593) = (cid:2013)/Ω . For the (cid:2013)′ = 0 case, our model returns to the original Dicke model, and the mean-field results describes a second-order transition (see Ref. [10]). When (cid:2018) < (cid:2018) (cid:3030) = 1 , the solution of the Eqs. (7, 8) is (cid:2009) = 0 and (cid:2016) = 0 . Thus the system is in the normal phase. Above the (cid:2018) (cid:3030) the Eqs. (7, 8) possess nonzero solution and both the field and the atoms acquire macroscopic occupations. This situation corresponds to the superradiant phase. For (cid:2013)′ ≠ 0 , this equation can be reduced to a quartic equation, and therefore its solution is not directly applicable to investigate the system. Here we focus on the positive (cid:2013) (cid:4593) case, and make some general statements about this solution. When (cid:2018) ≤ 1 , the left hand side of Eq. (9) is a monotonically decreasing function and therefore Eq. (9) only has one negative solution (cid:2016) (cid:2879) [Fig. 1 (c)]. For (cid:2018) > 1 , the left hand side of Eq. (9) has a local maximum in (cid:2016) > 0 and a local minimum in (cid:2016) < 0 [Fig. 1 (b)]. If the coupling constant (cid:2018) exceeds a certain critical value, there can be one negative solution and two positive solutions [Fig. 1 (a)]. In this case, it is easy to verify that the ground state is obtained when using the negative solution (cid:2016) (cid:2879) . This negative solution changes continuously with the coupling constant even when (cid:2018) goes across the critical value. For the positive solutions, the larger one is dynamically stable and is corresponded to the high-energy state, while the smaller one is dynamically unstable. Here we call this dynamically stable solution as (cid:2016) (cid:2878) . From the above discussion we can give an intuitive illustration of the effect of the bias term (cid:2013)(cid:1836) (cid:3051) . It acts as a symmetry breaking field and creates an asymmetric effective potential. For the positive bias, it lowers the energy functional at (cid:2016) (cid:2879) , thus the negative solution is always corresponded to the ground state.
Quantum corrections to the mean-field theory for the CO limit
We begin by deriving the effective model for the classical oscillator (CO) limit, where the ratio of the atomic transition frequency to the bosonic field frequency approaches infinity. According to the MF results, for nonzero (cid:2013) both the field and the atoms acquire macroscopic occupations, therefore we need to consider the fluctuations above the mean-field ground state |(cid:2016)(cid:1767)⨂|(cid:2009)(cid:1767) . To do this, we shift the field operator and rotate the spin operator as follows: (cid:1853) → (cid:1853) + (cid:2009); (cid:1836) (cid:3053) → (cid:1836) (cid:3053) cos(cid:2016) + (cid:1836) (cid:3051) sin(cid:2016). ( 10 )
Here the parameters (cid:2009) and (cid:2016) are corresponded to the mean-field state |(cid:2016)(cid:1767)⨂|(cid:2009)(cid:1767) , which are the solutions of Eqs. (7, 8). Making these transformations, the Hamiltonian of Dicke model becomes (cid:1834) = (cid:2033)(cid:1853) (cid:2878) (cid:1853) + (cid:2033)(cid:2009)((cid:1853) (cid:2878) + (cid:1853)) + (cid:2033)(cid:2009) (cid:2870) + Ω (cid:4593) (cid:3053) + (cid:2019)cos(cid:2016)(cid:1836) (cid:3051) ((cid:1853) (cid:2878) + (cid:1853)) − (cid:2019)sin(cid:2016)(cid:1836) (cid:3053) ((cid:1853) (cid:2878) + (cid:1853)), ( 11 ) where Ω (cid:4593) = Ωcos(cid:2016) − (cid:2013)sin(cid:2016) − 4(cid:2009)(cid:2019)sin(cid:2016) . In this new Hamiltonian, the mean-field ground state is |−(cid:1862)(cid:1767)⨂|0(cid:1767) thus no macroscopic occupation exists. For the Ω/(cid:2033) → ∞ limit, the low-energy part of the Hilbert space would be confined in the subspace ℋ ↓ . By using a unitary transformation (cid:1847)((cid:2010), (cid:2011)) = (cid:1857)(cid:1876)(cid:1868)(cid:3419)(cid:1861)(cid:1836) (cid:3052) (cid:4670)(cid:2010)((cid:1853) (cid:2878) + (cid:1853)) + (cid:2011)((cid:1853) (cid:2878) + (cid:1853)) (cid:2870) (cid:4671)(cid:3423) , we can decouple the ℋ ↓ and ℋ ↑ subspaces up to second order in (cid:2019)/Ω . Here (cid:2010)Ω (cid:4593) + (cid:2019)cos(cid:2016) = 0, (cid:2011)Ω (cid:4593) − 2(cid:2010)(cid:2019)sin(cid:2016) = 0. ( 12 ) Then we can obtain the effective low-energy Hamiltonian by projecting onto ℋ ↓ , which reads (cid:1834) (cid:3004)(cid:3016) = (cid:1766)−(cid:1862)|(cid:1847) (cid:2878) (cid:1834)(cid:1847)|−(cid:1862)(cid:1767) = (cid:2033)(cid:1853) (cid:2878) (cid:1853) + (cid:1840)(cid:2010)(cid:2019)cos(cid:2016)((cid:1853) (cid:2878) + (cid:1853)) (cid:2870) + (cid:2033)(cid:2009) (cid:2870) − (cid:1840) Ω (cid:4593) ( 13 ) With squeeze operator (cid:1845)((cid:2022)) , this Hamiltonian can be easily diagonalized to give (cid:1834) (cid:3004)(cid:3016)(cid:4593) = (cid:2033) (cid:4593) (cid:1853) (cid:2878) (cid:1853) + (cid:1829) . Here (cid:1829) is a constant we do not care about, and (cid:2022) = 14 ln (cid:4678)1 − 4(cid:1840)(cid:2019) (cid:2870) cos (cid:2870) (cid:2016)(cid:2033)Ω (cid:4593) (cid:4679), (cid:2033) (cid:4593) = (cid:2033)(cid:3496)1 − 4(cid:1840)(cid:2019) (cid:2870) cos (cid:2870) (cid:2016)(cid:2033)Ω (cid:4593) . ( 14 ) Now we have obtained the excitation energy (cid:2033) (cid:4593) . For (cid:2013)′ = 0 , the excitation energy is found to be (cid:2033) (cid:4593) = (cid:2033)√1 − (cid:2018) (cid:2870) and vanishes at (cid:2018) = 1 . This vanishing energy scale locates the QPT in the original Dicke model . For (cid:2013)′ > 0 , the low-energy properties are very different from the (cid:2013)′ = 0 case. When the coupling constant (cid:2018) is below the critical value, the Eq. (9) only has one negative solution (cid:2016) (cid:2879) , and this negative solution gives the low-energy sector of (cid:1834) . Above the critical coupling value, the Eq. (9) produces three solutions: one of them is an unstable stationary point, while the other two solutions (cid:2016) (cid:2879) and (cid:2016) (cid:2878) lead to two spectrums. For finite value of (cid:2013)′ , these two spectrums are not identical. From the numerical results we know that Ω (cid:4593) ((cid:2016) (cid:2879) ) > Ω (cid:4593) ((cid:2016) (cid:2878) ) for positive (cid:2013)′ , therefore the spectrum which corresponds to the positive solution is being lifted, while the other spectrum is being lower. This can be shown in Fig. (2). In this figure we plot the lowest 30 energy levels compare to the ground state as a function of the coupling constant (cid:2018) . Below the critical point, the low-lying energy levels are similar to that of a harmonic oscillator. Above the critical point, the energy spectrum turns into two equally spaced sets. The energy gap between these two sets is of the (cid:1841)((cid:2013)) . For the Ω/(cid:2033) → ∞ limit, this gap is much larger than the excitation energy (cid:2033) (cid:4593) . Therefore we can reasonably regard the lower one as the low-energy spectrum. It is obvious that the rescaled excitation energy (cid:2033) (cid:4593) /(cid:2033) only depends on the rescaled coupling constant (cid:2018) and the rescaled bias (cid:2013) (cid:4593) . In Fig. 3 we show the convergence to this analytical value when the CO limit is approached. This confirms that (cid:2033) (cid:4593) is the real excitation energy in the CO limit.
Crucially, we note that the excitation energy remains finite for arbitrary (cid:2018) in the presents of finite bias (cid:2013)′ . According to the Eq. (14), the excitation energy vanishes only if tan(cid:2016) = √(cid:2013)′ (cid:3119) , and can be never fulfilled for the negative solution. We show in Fig. 4 the rescaled excitation energy as a function of the rescaled coupling constant (cid:2018) and the rescaled bias (cid:2013) (cid:4593) . As in this figure we see that the excitation energy becomes zero only at (cid:2013)′ = 0 and (cid:2018) = 1 , which indicates the QPT. This scenario is strongly suppressed as (cid:2013) (cid:4593) is increased, and there is no longer any sign of a critical point. Thus we conclude that the characteristic energy scale of fluctuations above the ground state does not vanish for arbitrary coupling, and the second order phase transition cannot occur for finite bias. For the positive solution, the relation tan(cid:2016) = √(cid:2013)′ (cid:3119) can be satisfied when the coupling reaches the critical value, and it seems that the excitation energy would vanish in this case. However, we note that the positive solution coincides with the unstable stationary point when this relation is satisfied, and therefore it does not correspond to a real spectrum. For finite field frequency, the coupling between two spectrums becomes much stronger as the coupling approaches the critical value. By considering this coupling, the eigenstates of the high-energy sector are given by the superposition of two different spectrums. Therefore the excitation energy of the high-energy spectrum does not become zero even in the vicinity of the critical coupling. From previous discussion we see that the excitation energy vanishes only at (cid:2013)′ = 0 and (cid:2018) = 1 , demonstrating the existence of the QPT. We now discuss the critical behavior of the system as bias tends to zero. For (cid:2018) = 1 , the excitation energy vanishes as (cid:2013)′ → 0 from either direction. In the (cid:2013)′ → 0 (cid:2878) limit , the Eq. (9) can be reduced to (cid:2016) = −(2(cid:2013)′) (cid:2869)(cid:2871) . ( 15 ) Then the excitation energy (cid:2033) (cid:4593) can be shown to vanish as (cid:2033) (cid:4593) ~(cid:2033)(cid:2013)′ (cid:2869)(cid:2871) . ( 16 ) Meanwhile, we identify the oscillator variance
Δ((cid:1853) + (cid:1853) (cid:2878) ) as the characteristic length scale. Our later calculations will show that this length diverges as |(cid:2013) (cid:4593) | (cid:2879)(cid:2869)/(cid:2874) , from which we find that (cid:1878) = 2 and (cid:1874) = 1/6 . Here (cid:1874) is the critical exponent and (cid:1878) is the dynamic critical exponent. The role of A-square term
In the minimal-coupling Hamiltonian, a term containing the square of the vector potential
A(cid:3553) (cid:2870) is present. Although this term can be neglected in most cases, it can prevent the SPT in cavity QED systems . We now consider the critical behavior of the biased Dicke model which contains the
A(cid:3553) (cid:2870) term. With the
A(cid:3553) (cid:2870) term included, an additional term κ (cid:2868) (a (cid:2878) + a) (cid:2870) is added in the biased Dicke Hamiltonian. The parameter λ and κ (cid:2868) are not independent of each other, and are related by κ (cid:2868) =α (cid:2868) Nλ (cid:2870) /Ω . Here α (cid:2868) is an independent parameter decided by the field and the atoms. Then the Eq. (9) should be modified as (cid:4678) (cid:2018) (cid:2870) (cid:2009) (cid:2018) (cid:2870) − 1cosθ(cid:4679) sin(cid:2016) = (cid:2013) (cid:4593) ( 17 ) The left hand side of Eq. (17) is a monotonically decreasing function for arbitrary coupling κ since α (cid:2868) > 1 is always satisfied according to the Thomas-Reiche-Kuhn sum rule. Thus the critical condition cannot be reached if the A(cid:3553) (cid:2870) term is not neglected, which is known as the no-go theorem . The effective low-energy Hamiltonian can also be calculated by repeating the same step in the previous section, which reads (cid:1834) (cid:3004)(cid:3016) = (cid:2033)(cid:1853) (cid:2878) (cid:1853) + ((cid:1840)(cid:2010)(cid:2019)cos(cid:2016) + (cid:2018) (cid:2868) )((cid:1853) (cid:2878) + (cid:1853)) (cid:2870) + (cid:1829) (cid:2869) . ( 18 ) Here C (cid:2869) is a constant we do not care about. Then the excitation energy is given as (cid:2033) (cid:4593) = (cid:2033)(cid:3496)1 − 4(cid:1840)(cid:2019) (cid:2870) cos (cid:2870) (cid:2016)(cid:2033)Ω (cid:4593) + 4(cid:2018) (cid:2868) (cid:2033) = (cid:2033)(cid:3496)1 − ΩΩ (cid:4593) (cid:2018) (cid:2870) cos (cid:2870) (cid:2016) + (cid:2009) (cid:2868) (cid:2018) (cid:2870) . ( 19 ) We note that the excitation energy can never vanish for arbitrary α (cid:2868) provided ε (cid:4593) ≠ 0 . Noticing that parameter θ here is the negative solution to the Eq. (9) if we replaced κ with an effective coupling constant κ/(cid:3493)1 + α (cid:2868) κ (cid:2870) . Since the excitation energy of the biased Dicke model remains finite for any coupling, we have (cid:4593) (cid:2018) (cid:2870) (cid:2009) (cid:2018) (cid:2870) cos (cid:2870) (cid:2016) > 0. ( 20 ) With this inequility we find ω (cid:4593) ω⁄ > (cid:3493)1 − (1 + α (cid:2868) κ (cid:2870) ) + α (cid:2868) κ (cid:2870) = 0 , which indicates that in the CO limit the SPT cannot occur in the finite biased case regardless of the A(cid:3553) (cid:2870) term.
Ground state wave function and its squeezing and entanglement
We now consider the ground state wave function of the system in CO limit. In contrast to the doubly degenerate ground state in the original Dicke model, for finite bias the ground state remains nondegenerate even above the critical coupling. According to our previous discussion, the ground state wave function in the CO limit is given as |(cid:2032)((cid:2016) (cid:2879) )(cid:1767) = (cid:1857) (cid:2869)(cid:2870)(cid:3036)(cid:3087) (cid:3127) (cid:3011) (cid:3300) (cid:1830)((cid:2009))(cid:1847)((cid:2010), (cid:2011))(cid:1845)((cid:2022))|0(cid:1767)|(cid:1862), −(cid:1862)(cid:1767), ( 21 ) where (cid:1830)((cid:2009)) is the displacement operator. In this expression the operator exp(cid:3435)(cid:1861)(cid:2016) (cid:2879) (cid:1836) (cid:3052) /2(cid:3439) (cid:1830)((cid:2009)) comes from the mean-field results, (cid:1847)((cid:2010), (cid:2011)) is corresponded to the spin-oscillator correlation, while the operator (cid:1845)((cid:2022)) squeezes the bosonic field and produces the field fluctuations. It is interesting that no operator except (cid:1847)((cid:2010), (cid:2011)) in this expression is responsible for the spin fluctuations. Since the operator (cid:1847)((cid:2010), (cid:2011)) corresponds to the quantum corrections for the finite (cid:2033) case, we can expect that the spin fluctuations vanish in the CO limit. In fact, the spin variance
Δ(cid:1836) (cid:3051) = 〈(cid:1836) (cid:3051)(cid:2870) 〉 − 〈(cid:1836) (cid:3051) 〉 (cid:2870) for small (cid:2033)/Ω is of order (cid:1841)((cid:2010) (cid:2870) ) , or equivalently (cid:1841)((cid:2033) (cid:2870) /Ω (cid:2870) ) . The vanishing of Δ(cid:1836) (cid:3051) as (cid:2033)/Ω → 0 shows that the spin fluctuations are strongly suppressed in the CO limit. In order to discuss the fluctuations and squeezing of the field, we use the variance (Δx) (cid:2870) and (Δp) (cid:2870) of the field position operator x = (a (cid:2878) + a)/√2ω and the momentum operator p = i(a (cid:2878) − a)(cid:3493)ω/2 . Using the results from the Eq. (21) the variances are obtain as (Δ(cid:1876)) (cid:2870) = 12(cid:2033) (cid:1857) −2(cid:2022) , (Δ(cid:1868)) (cid:2870) = (cid:2033)2 (cid:1857) + (cid:1841) ( (cid:2010) ) . ( 22 ) Thus, we see that the field becomes a squeezed coherent state |α, ξ(cid:1767) in the CO limit. In Fig. 5 (left panel) we plot the parameter ξ as a function of the rescaled coupling constant κ . The parameter ξ diverges at the QPT point ( ε′ = 0, κ = 1 ), which shows that the momentum variance is being strongly squeezed in the vicinity of the QPT point, and (Δx) (cid:2870) becomes strongly antisqueezed. This result is different from the CS limit case, where only a slight squeezing of the momentum variance presents as the system approaches the QPT point . The field fluctuations, which can be described by (Δx) (cid:2870) , diverge only at the QPT point. Near the QPT point, the parameter ξ behaves as e (cid:2870)(cid:2974) ∝ (ε (cid:4593) ) (cid:2879)(cid:2869)/(cid:2871) . Therefore both the fluctuations and the squeezing drop rapidly as the rescaled bias ε (cid:4593) is increased from zero. This is agreed with our previous results, where the system experiences QPT only at the point ε′ = 0 and κ = 1 . In addition to the fluctuations and the squeezing, we now study the spin-field entanglement. Here we use the entropy S = −Tr(ρlnρ) to quantify the entanglement, which is calculated with the reduced spin density matrix ρ . For ε (cid:4593) = 0 case, the ℤ (cid:2870) symmetry in the original Dicke model will leads to symmetrized ground-state wave function |ψ(cid:1767) = (|ψ(θ (cid:2879) )(cid:1767) ± |ψ(θ (cid:2878) )(cid:1767))/√2 when the system is in the superradiant phase. Thus the ground states obtain non-zero entanglement for κ > 1 in the CO limit. This entanglement will be destroyed in the present of finite bias due to the absence of ℤ (cid:2870) symmetry and the ground state becomes nondegenerate. For ε (cid:4593) ≠ 0 case, the ground state wave function can be written as |ψ(θ (cid:2879) )(cid:1767) = U(β, γ)|α, ξ(cid:1767)|θ (cid:2879) (cid:1767) . It is obvious that the entanglement is created completely by the operator U(β, γ) . Based on our previous discussion, the parameters β and γ become zero as ω/Ω → 0 , which indicates that the spin-field entanglement vanish in the CO limit. In Fig. 5 (right panel) we plot the entropy as a function of the rescaled coupling constant κ for different ω/Ω . The vanishing of entropy as ω/Ω → 0 shows the suppression of entanglement in the CO limit. Different from the CO limit, the spin fluctuations are not suppressed in the CS limit. The effective model for the CS limit can be obtained by performing the Holstein-Primakoff transformation (cid:1836) (cid:2878) = 2(cid:1854) (cid:2878) (cid:3493)2(cid:1862) − (cid:1854) (cid:2878) (cid:1854), (cid:1836) (cid:2879) = 2(cid:3493)2(cid:1862) − (cid:1854) (cid:2878) (cid:1854)(cid:1854), (cid:1836) (cid:3053) = 2((cid:1854) (cid:2878) (cid:1854) − (cid:1862)), ( 23 ) which is given as (cid:1834) (cid:3004)(cid:3020) = (cid:2033)(cid:1853) (cid:2878) (cid:1853) + Ω (cid:4593) (cid:1854) (cid:2878) (cid:1854) + √(cid:1840)(cid:2019)cos(cid:2016)((cid:1854) (cid:2878) + (cid:1854))((cid:1853) (cid:2878) + (cid:1853)) − (cid:1840)Ω (cid:4593) . ( 24 ) This Hamiltonian is bilinear in the bosonic operators and can be simply diagonalized by the Bogoliubov transformation. Here we use the position operator and the momentum operator of bosonic mode (cid:1854) in Eq. (24) to describe the squeezing in the spin. Then the variances in the ground state can be calculate, which are given as (Δ(cid:1876) (cid:3028) ) (cid:2870) = − cos (cid:2026) + 12(cid:2033) + sin (cid:2026), (Δ(cid:1876) (cid:3029) ) (cid:2870) = − sin (cid:2026) + 12(cid:2033) + cos (cid:2026), (Δ(cid:1868) (cid:3028) ) (cid:2870) = (cid:2033) − (cid:2026) + (cid:2033) + (cid:2026), (Δ(cid:1868) (cid:3029) ) (cid:2870) = (cid:2033) − (cid:2026) + (cid:2033) + (cid:2026), ( 25 ) where the subscripts denote the bosonic modes (cid:1853) and (cid:1854) respectively. In Fig. 6, we show the analytical values of these variances as a function of κ for different biases. We see that a finite bias will diminish the sharp increase of (Δ(cid:1876) (cid:3028) ) (cid:2870) and (Δ(cid:1876) (cid:3029) ) (cid:2870) as κ approaches QPT point in the zero biased case. The slight squeezing of (Δ(cid:1868) (cid:3028) ) (cid:2870) and (Δ(cid:1868) (cid:3029) ) (cid:2870) is also weakened when bias is increased. However, the bias has little effect on the squeezing as the coupling constant κ is far from the QPT point, and the variances stay largely constant for different bias. Conclusion
We have analyzed the properties of the biased Dicke model in the CO limit. For finite bias, the mean-field results show that the ground state remains nondegenerate even above the critical coupling. The low-energy effective Hamiltonians are found for the CO limit. We find that the low-energy spectrum turns into two equally spaced spectrums when the coupling is above the critical value, and the energy gap between these two spectrums is much larger than the excitation energy, thus we can reasonably regard the lower one as the low-energy spectrum. For finite bias, the excitation energy remains finite for arbitrary coupling. This indicates that the second order phase transition is avoided for finite bias in the CO limit. We then discuss the excitation energy of the biased Dicke model in the presence of the (cid:2157)(cid:3553) (cid:2870) term. The resulted effective low-energy Hamiltonian is given and we show that, the SPT cannot occur for the finite biased case regardless of the (cid:2157)(cid:3553) (cid:2870) term. The results also demonstrate that the (cid:2157)(cid:3553) (cid:2870) term prevents the system from reaching the critical coupling, thus can be compensated by other interaction terms. However, a finite bias will always break the ℤ (cid:2870) symmetry and destroy the superradiant phase completely. Thus the SPT is hard to recover unless we manage to restore the symmetry. We have also calculated the ground state and discuss its squeezing in both limits. The results show that the bias term will suppress the squeezing strongly in the vicinity of the QPT point, and has little effect as the coupling is far from the QPT point. Moreover, the squeezing properties are quite different in two classical limits. In the CO limit, the momentum variance of the field is being strongly squeezed in the vicinity of the QPT point, in comparison with the slight squeezing of the momentum variance in the CS limit. The entanglement of the spins and the field will behave dramatically different for finite bias. Instead of remaining a finite value , the entanglement vanishes as we approach the CO limit. Thus the entanglement of the spins and the field will be suppressed in the presence of the bias term. [1] Dicke, R. H. Coherence in Spontaneous Radiation Processes.
Phys. Rev. (1), 99-110 (1954). [2] Scully, M. O. and Zubairy M. S. in Quantum Optics 1 st edn, Ch. 6, 194-196 (Cambridge University Press, 1997). [3] Hepp, K. and Lieb, E. H. On the superradiant phase transition for molecules in a quantized radiation field: the Dicke maser model.
Ann. Phys. (2), 360-404 (1973). [4] Emary, C. and Brandes, T. Quantum chaos triggered by precursors of a quantum phase transition: The Dicke model.
Phys. Rev. Lett. , 044101 (2003). [5] Emary, C. and Brandes, T. Chaos and the quantum phase transition in the Dicke model,
Phys. Rev. E. (6), 066203 (2003). [6] Dimer, F., Estienne, B., Parkins, A. S. and Carmichael, H. J. Proposed realization of the Dicke-model quantum phase transition in an optical cavity QED system.
Phys. Rev. A. , 013804 (2007). [7] Bastidas, V. M., Emary, C., Regler, B. and Brandes, T.
Nonequilibrium Quantum Phase Transitions in the Dicke Model.
Phys. Rev. Lett. , 043003 (2012). [8] Baumann, K. et al ., Dicke quantum phase transition with a superfluid gas in an optical cavity.
Nature .7293, 1301-1306 (2010). [9]
Baumann, K., Mottl, R., Brennecke, F. and Esslinger, T. Exploring Symmetry Breaking at the Dicke Quantum Phase Transition.
Phys. Rev. Lett. , 140402 (2011). [10]
Bakemeier, L., Alvermann, A. and Fehske, H. Quantum phase transition in the Dicke model with critical and non-critical entanglement.
Phys. Rev. A. , 043821 (2012). [11] Ashhab, S. and Franco N, Qubit-oscillator systems in the ultrastrong-coupling regime and their potential for preparing nonclassical states.
Phys. Rev. A. , 042311(2010). [12] Ashhab, S. Superradiance transition in a system with a single qubit and a single oscillator.
Phys. Rev. A. , 013826 (2013). [13] Wallraff, A. et al ., Strong coupling of a single photon to a superconducting qubit using circuit quantum electrodynamics.
Nature , 162-167 (2004). [14]
Devoret, M. H. et al ., Circuit ‐ QED: How strong can the coupling between a Josephson junction atom and a transmission line resonator be?
Annalen Der Physik .10-11, 767–779 (2007). [15] Niemczyk, T. et al ., Circuit quantum electrodynamics in the ultrastrong-coupling regime.
Nature Physics , 772-776 (2010). [16] Scalari, G. et al., Ultrastrong coupling of the cyclotron transition of a 2D electron gas to a THz metamaterial.
Science
Forn-Díaz, P., Lisenfeld, J., Marcos, D., García-Ripoll, J. J., Solano, E., Harmans, C. J. P. M. and Mooij, J. E. Observation of the Bloch-Siegert shift in a qubit-oscillator system in the ultrastrong coupling regime.
Phys. Rev. Lett. , 237001 (2010). [18]
Fedorov, A., Feofanov, A. K., Macha, P., Forn-Díaz, P., Harmans, C. J. P. M. and Mooij, J. E. Strong Coupling of a Quantum Oscillator to a Flux Qubit at Its Symmetry Point.
Phys. Rev. Lett. , 060503 (2010). [19]
Xiang, Z. L. et al ., Hybrid quantum circuits: Superconducting circuits interacting with other quantum systems.
Rev. Mod. Phys . , 623-653 (2013). [20] Hausinger, J. and Grifoni, M. Qubit-oscillator system: an analytical treatment of the ultra-strong coupling regime.
Phys. Rev. A. , 062320 (2010). [21] Zhang, Y. Y., Chen, Q. H. and Zhao, Y. Generalized rotating-wave approximation to biased qubit-oscillator systems.
Phys. Rev. A. , 033827 (2013). [22] Emary, C. and Brandes, T. Phase transitions in generalized spin-boson (Dicke) models.
Phys. Rev. A. , 053804 (2004). [23] Genway, S. et al ., Generalized Dicke Nonequilibrium Dynamics in Trapped Ions.
Phys. Rev. Lett. , 023603 (2014). [24]
Porras, D. et al ., Quantum Simulation of the Cooperative Jahn-Teller Transition in 1D Ion Crystals.
Phys. Rev. Lett. , 235701 (2012). [25]
Rza ż ewski, K. et al., Phase Transitions, Two-Level Atoms, and the A Term.
Phys. Rev. Lett. , 432 (1975). [26] Nataf, P. and Ciuti, C. No-go theorem for superradiant quantum phase transitions in cavity QED and counter-example in circuit QED.
Nat. Commun . , 72 (2010). [27] Viehmann, O. et al ., Superradiant Phase Transitions and the Standard Description of Circuit QED.
Phys. Rev. Lett. , 113602 (2011). [28]
Vukics, A. et al., Elimination of the A-Square Problem from Cavity QED.
Phys. Rev. Lett. , 073601 (2014).
Figure captions
Figure 1: Graphical solutions of Eq. (9) for different values of coupling (cid:2018) . The (cid:1853), (cid:1854), (cid:1855) lines represent the left side in Eq. (9), while the horizontal line is the right side in the equation.
Figure 2: The rescaled energy levels of the first 30 excited states relative to the ground state energy as functions of the coupling strength (cid:2018) . Here we take
Ω (cid:2033)⁄ = 60 , (cid:2013) (cid:4593) = 0.11 with (cid:1840) = 5 . Figure 3: The rescaled excitation energy as functions of the coupling strength κ . The dashed line represents the analytical result in Eq. (14), while the a, b, c lines are the excitation energies of the first excited state relative to the ground state. Here we take (cid:2013) (cid:4593) = 0.11 with (cid:1840) = 5 . Figure 4: The excitation energies as functions of the coupling strength κ and the bias ε′ in the CO limit.
Figure 5: (Left panel) The squeezing parameter ξ of the ground state as a function of coupling κ for different biases in the CO limit. (Right panel) The entanglement entropy S of the ground state as a function of coupling κ for different Ω/(cid:2033) . Both panels start from the curve for (cid:1840) = 5 . Figure 6: The squeezing variances of the field (left panel) and the spin (right panel) of the ground state for different biased cases in the CS limit. Here solid lines denote the variances of the position operators, whereas dashed lines correspond to the variances of the momentum operators. The Hamiltonian is on scaled resonance: (cid:2033) = Ω = 1 . Acknowledgments
This work is supported by the National Natural Science Foundation of China (Grant No. 11574022 and 11174024) and the Open Fund of IPOC (BUPT) grants Nos. IPOC2013B007, also supported by the Open Project Program of State Key Laboratory of Low-Dimensional Quantum Physics (Tsinghua University) grants Nos. KF201407, and Beijing Higher Education (Young Elite Teacher Project) YETP 1141.
Author contributions
G.Z. conceived and designed the research. G.Z. and H.Z. performed analysis and wrote the manuscript. H.Z. prepared all the figures. H.F. reviewed the paper and gave some valuable suggestion.
Additional information
Competing financial interests: The authors declare no competing financial interests. κ b: = 1.2 κ a: = 1.5 κ a b c ε Ω/ θ θ - θ - θ - θ + κ ( E - E ) / n0 ω a: Ω ω / = 200 b: Ω ω / = 60 c: Ω ω / = 20 κ ( E - E ) / ω a bc ω ’ / κ ε’ a: 0 ε΄ =b: 0.01 ε΄ =c: 0.05 ε΄ =d: 0.1 ε΄ = abcd κξ S a: Ω ω / = 10b:
Ω ω / = 20c:
Ω ω / = 40 abc κ Δ κ a: 0 ε΄ =b: 0.01 ε΄ =c: 0.05 ε΄ = abc κ Δ a: 0 ε΄ =b: 0.01 ε΄ =c: 0.05 ε΄ ==