Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 10568.
Published online 2017 September 5. doi:  10.1038/s41598-017-10686-3
PMCID: PMC5585414

Quantum mixed phases of a two-dimensional polarized degenerate Fermi gas in an optical cavity


The coupling of ultracold fermions to a high-finesse optical cavity can result in novel many-body phenomena, and has attracted significant interests at present. Here we consider a realization of the Fermi-Dicke model with controllable parameters, based on a two-dimensional polarized degenerate Fermi gas coupled to an optical cavity. We analytically investigate the ground-state properties of such system under the mean-field approximation. We find the system can exhibit a rich phase diagram depending on the fermion-photon coupling strength and the atomic resonant frequency. Contrasting to the bosonic counterpart, a first-order quantum phase transition between the superradiant phase and the normal phase featuring two Fermi surfaces can occur for the weak atomic resonant frequency, and there is a unique mixed phase where this normal phase and the superradiant phase coexist. The experimental detection of our results is also discussed.


Cavity quantum electrodynamics (QED) systems, which remarkably illuminate the fundamental interaction between light and matter, have emerged as a novel platform to explore the many-body physics, and promise potential applications in quantum information processing and quantum computing. Recently, the coupling of a Bose-Einstein condensate to a high-finesse optical cavity has been experimentally achieved18, thus opens a new avenue that combines cavity QED with ultracold atoms. In particular, the bosonic atoms occupying the same quantum state can interact identically with a single-mode quantized field, which can result in a strong collective matter-field interaction. This has led to the remarkable experimental observation35 of the celebrated second-order quantum phase transition from the normal phase to the superradiant (SR) phase predicted more than 40 years ago9, 10. Cavities moreover allow unconventional dynamical optical potentials which can induce a rich variety of strongly correlated many-body phenomena11.

While current cavity QED experiments have focused on bosons, there are surging interests in exploring the novel physics arising from the coupling of ultracold fermions to the optical cavity1225. Unlike bosons, fermions obey the Pauli exclusion principle, and two fermions with weak attractive interaction can form Cooper pairs which are responsible for superconductivity26. By coupling ultracold fermions to a high-finesse optical cavity, exotic phenomena have been predicted to arise. For example, recent studies on spinless fermions in the cavity-induced dynamical optical potential have revealed the crucial role of the Fermi statistics on the SR phase transition at moderate and high densities1719. It has been shown that the cavity-assisted spin-orbit coupling27, 28 can induce a topological SR phase20. Moreover, the cavity-induced artificial magnetic field22, chiral phases23, and non-trivial topological states24 have been reported. Interestingly, when fermions in an infinite lattice are gauge coupled to a cavity mode, a SR phase exhibiting a directed particle flow is found to arise above an infinitesimal pumping threshold25.

In this report, we consider a two-dimensional (2D) polarized degenerate Fermi gas coupled to a high-finesse optical cavity and realize a Fermi-Dicke model. Our setup relies on two Raman transitions induced by the quantized cavity field and two transverse pumping lasers, and allows for flexible controllability of the Hamiltonian parameters, including the fermion-photon coupling strength and the effective Zeeman field represented by the atomic resonant frequency. Based on this microscopic model, we study the ground state under the mean-field approximation and find several distinct properties compared to its bosonic counterpart9, 10. In particular, we predict a first-order quantum phase transition between the SR phase and the normal phase featuring two Fermi surfaces for the weak atomic resonant frequency. As such, the system exhibits a unique mixed phase where the normal and SR phases coexist. Finally, we discuss possible experiment observations of our results.


Model and Hamiltonian

Motivated by the experiment with bosonic atoms5, we instead consider an ensemble of ultracold four-level fermions coupled to a high-finesse optical cavity; see Fig. 1. The fermions are confined in a far-off-resonance optical trap in the yz plane [see Fig. 1(a)], and their motion in the x direction is frozen, thus effectively realizing a 2D system. The cavity mode is driven by a linearly polarized laser, while the fermions are pumped by two transverse lasers which are left- and right-circular polarized in the yz plane, respectively. We consider fermionic atoms with four internal levels, which contain two degenerate ground states (|↑right angle bracket and |↓right angle bracket) and two excited states (|1right angle bracket and |2right angle bracket) [see Fig. 1(b)]. As clearly illustrated in Fig. 1, in our setup two Raman processes can be induced from the quantized cavity field and the two transverse pumping lasers.

Figure 1
An ensemble of ultracold four-level fermions coupled to a high-finesse optical cavity. (a) Schematic of the proposed setup: the ultracold fermions (black online) are confined in a far-of-resonance optical trap (yellow online) in the yz plane along with ...

The considered system can be effectively described by the following 2D Hamiltonian


where HˆF is the Hamiltonian for the free four-level fermions, i.e.,


Here, ψˆi(r) and ψˆi(r) (i = 1, 2,  ↓ ,  ↑ ) are the field operators for the fermionic atoms with mass M, μ is the chemical potential, and ω i denotes the frequency of the ith eigenstate. The quantized cavity field and the driving laser are described by the Hamiltonian


where aˆ (aˆ) creates (annihilates) the quantized cavity field with frequency ω c, and ε (ω l) labels the magnitude (frequency) of the driving laser. The Hamiltonian HˆAR represents the interaction between fermions and the transverse pumping lasers, which under the rotating-wave approximation can be written as


with Ω1 and Ω2 (ω A and ω B) being the Rabi frequencies (frequencies) of the two lasers. Finally, Hamiltonian HˆAP describes the interaction between fermions and the quantized cavity fields, i.e.,


with g 1 and g 2 labeling the fermion-photon coupling strengths associated with the two pumping lasers, respectively.

It is more transparent to recast the time-dependent Hamiltonian (1) into a time-independent form. Introducing a unitary transformation Uˆ(t)=exp(iHˆ't), where


with ω l = (ω B+ω A)/2, we transform the Hamiltonian as Hˆ1=Uˆ(t)HˆT(t)Uˆ(t)+i[Uˆ(t)/t]Uˆ(t). Explicitly, we have


Here, ω˜=ωcωl is the effective cavity frequency, Δ1ω1 − ωA/2 (Δ2=ω2ωB/2) denotes the detuning from the excited state |1right angle bracket (|2right angle bracket), and ω˜= ωωA/2 (ω˜=ω+ωB/2) describes the effective frequency associated with the internal state |↑right angle bracket (|↓right angle bracket).

To proceed, we recall that in the experiments35, a weak driving (ε  0) and large detuning (|Δ1,2|{Ω1,2,g1,2,ω˜,ω0}) have been considered. Assuming similar scenarios here, this allows us to ignore the term ε(aˆ+aˆ) in the Hamiltonian (7), as well as adiabatically eliminate both excited states |1right angle bracket and |2right angle bracket29, 30. This way, we obtain


Further considering


the Hamiltonian (8) can be simplified into


Here, the factor 1/N, with N being the total atom number, is introduced so as to ensure a finite free energy per fermion in the thermodynamic limit9, 10. Furthermore, ω 0 = (ω˜ω˜)/2 is the effective resonant frequency between the ground states |↑right angle bracket and |↓right angle bracket, which in Eq. (10) acts as an effective Zeeman field. Without loss of generality, below we shall take ω0 > 0. In addition, the parameter η=Ng1Ω1/(2Δ1) =Ng2Ω2/(2Δ2) is the effective fermion-photon coupling strength, and ω=Nζ+ω˜ labels the atom-number dependent cavity frequency with ζ = |g1|21 = |g2|22. We remark that all the parameter choice here are motivated by the experimental considerations35.

The Hamiltonian (10) represents the paradigmatic Fermi-Dicke model31 describing the fermion-photon interaction. We emphasize that our setup allows flexible controllability of all Hamiltonian parameters. For example, both ω 0 and ω can be tuned by modifying the frequencies of the driving laser and the transverse pumping lasers, while η can be controlled via the Rabi frequencies of the transverse pumping lasers.

Ground-state properties

Our goal is to investigate the ground state of the Hamiltonian (10). To this end, it is more convenient to transform to the momentum space representation. Writing


where Cˆk,σ annihilates a fermion in the internal state σ with momentum k and S is the system size (hereafter S = 1 is set for convenience), we obtain


Here, ξkεk − μ, εkk2/2M is the kinetic energy, and n=KF2/(2π)=EFM/π is the 2D density of fermions with EF=KF2/(2M) being the Fermi energy and K F being the Fermi momentum.

In solving the ground state of the Hamiltonian (12), we will rely on the mean-field approximation, i.e., by replacing aˆ with its steady-state value. Specifically, we write down the Heisenberg-Langevin equation for the cavity field operator aˆ 32, 33, taking into account of the cavity decay with rate κ, i.e.,


For the quantum noise operator γˆin(t), we have γˆin(t)γˆin(t')=2κδ(tt') and γˆin(t)γˆin(t')=0. Since quantum noise usually occurs on a much shorter time scale than 1/κ 34, its average effect can be generically ignored on the time scale relevant for the steady state. When 1/κ is much shorter than the time scales of system dynamics, a steady-state solution3, 35 to Eq. (13) can be found, i.e.,


Motivated by the experiments3, 4 which shows that the mean-photon number aˆaˆ=|α|2 determines the SR properties, we henceforth refer to it as the SR order parameter.

By approximating aˆaˆ in Eq. (12) using Eq. (14), we obtain a quadratic Hamiltonian


Here, Ψˆk=(Cˆk,,Cˆk,)T is the standard Nambu spinor, and M k is the photon-number dependent BdG matrix, i.e.,


with η¯=η(α+α)/n. Now, the Hamiltonian (15) can be easily diagonalized as


where αˆk,± describe the fermionic Bogoliubov quasiparticles, whose energy is given by


with η¯2=4ω2η2|α|2/[n(ω2+κ2)]. Note that Ek can be either positive or negative depending on parameter choices. Regrouping terms with positive and negative energies using the Heaviside step function Θ(x), we recast Eq. (17) as


where E G is the ground-state energy expressed by


Despite the formal expression of the ground-state energy, Eq. (20) contains two variables, i.e., μ and |α|, which are to be determined from the particle number equation EG/∂μ = − n and the SR equation EG/∂(|α|) = 0. Explicitly, we have



where f(  −  Ek) = Θ(−Ek) − Ekδ(−Ek) and χ¯=η¯2+ω02. In deriving Eqs (21) and (22), we have used the identity Θ'(x) = δ(x) with δ(x) being the delta function.

Equations (2022) allow complete specification of the ground state: The ground-state energy can be obtained in a self-consistent manner by solving Eqs (21) and (22) for a fixed atom density n. It is important to bear in mind that solutions of above equations must be supplemented with a careful stability analysis, which we describe below.

Phase diagram

We now detail our analysis on the ground-state properties. We shall be interested in four quantities: the ground-state energy per fermion E¯G=EG/n, the chemical potential μ, the scaled mean-photon number |α¯|2, and the scaled polarization36 measuring the response to the effective Zeeman field defined by


As we shall see, their ground-state values depend crucially on the strength of the effective Zeeman field ω 0 and the effective fermion-photon coupling strength η, in particular, new features compared to their boson counterpart are found in the regime ω0 < EF. Note that all energies will be measured in units of E F hereupon.

In order to gain some intuition, let us first consider the simplest case η = 0 corresponding to a free Fermi gas. Transforming Eq. (20) into an integral over momenta as usual, we calculate the ground-state energy per fermion as


Because μω0 < 0 will entail a unphysical result E¯G=0 which excludes existence of real fermions37, 38, we shall limit our subsequent discussions in the regime μω0 > 0. There, the Fermi gas can exhibit different ground states depending on whether ω0 > μ or ω0 < μ.

When ω0 > μ, Θ(μ − ω0) = 0 and Eq. (24) becomes


Based on Eqs (21), (23) and (25), we obtain


Thus the Fermi gas is in a normal phase characterized by a full polarization and one Fermi surface with μ = 2EF, and we shall refer to it as the N-I phase.

In contrast, when ω0 < μ, Θ(μ − ω0) = 1 and Eq. (24) becomes


Based on Eqs (21), (23) and (25), we have


Hence the Fermi gas is only partially polarized and exhibits two Fermi surfaces, i.e., μEFω0 and μEF − ω0. We shall call it as the N-II phase. It is clear from Eqs (26) and (28) that, by tuning the parameter ω 0 across the critical point ω0c=EF 37, the Fermi gas undergoes a first-order transition between the N-I to N-II phases.

Now, suppose there exists a weak effective coupling between the photon and fermions, i.e., η is small. Still, the noninteracting terms in the Hamiltonian (12) dominate the system dynamics, such that |α|2 = 0 is expected to persist for the ground state and the Fermi gas remains in the N-I or N-II phase (depending on values of ω 0).

However, the ground-state properties of the Fermi gas can change drastically from both normal phases for a strong photon-fermion coupling when η is increased above a threshold. As we show now, when the interaction part of the Hamiltonian (12) dominates, the system can acquire a macroscopic collective excitation with |α|2  0, i.e., a quantum phase transition into the SR phase occurs when η is above some critical values. Integrating Eq. (20) over the momenta, we obtain


with μ±=μ±χ¯ and |α¯|2=|α|2/n being the scaled mean-photon number. We next analyze Eq. (29) when (i) μ < 0 and (ii) μ ≥ 0.

Consider first the case μ < 0, where Eq. (29) becomes


and Eqs (2123) take the following form




We can obtain two set of solutions, namely,




Obviously, while solutions (34) are associated with the normal phase, solutions (35) exhibit SR property. Furthermore, as mentioned earlier, the ground-state solutions must satisfy the stability condition as defined by 2E¯G/(|α¯|)2>0. Consequently, we find that the stability condition for the normal solution [see Eq. (34)] is η<ηc(1) with


Instead, the SR solution [see Eq. (35)] is only stable when η>ηc(1). Notice that since μ=μχ¯<0, this additionally requires η>ηc(2) with


Following from above analysis, we find that for ω0 ≥ EF [c.f. Fig. 2(a)] and thereby ηc(1)>ηc(2), the ground state of the Fermi gas is in the N-I phase for 0<η<ηc(1), but transits into the SR phase for η>ηc(1); see summarized results below


with corresponding relevant quantities given by




Figure 2
Comparison of the critical points ηc(1) and ηc(2) for (a) ω0 ≥ EF and (b) ω0 < EF. When ω0 ≥ EF, ηc(1)>ηc(2), and μ < 0 for both 0<η<ηc(1) and η>ηc(1). When ...

The ground-state properties described by Eqs (38)–(41) exhibit two features. First, the quantum phase transition between the N-I and SR phases is second order in nature: while the first-order derivative of E¯G is continuous with respect to η, its second order derivative develops discontinuity at the critical point ηc(1), as in the case with ultracold Bose atoms9, 10. Second, in contrast to the N-I phase, the SR phase features macroscopic collective excitations for both fermions and photons, i.e., we have both |α¯|0 and m¯0. Also note that while m¯η2, we see μ ∼ η−2. These analytical results are confirmed by our numerical calculations as illustrated by Fig. 3, where we plot E¯G, μ, |α¯|2, and m¯ as functions of η, respectively.

Figure 3
(a) The scaled ground-state energy E¯G/EF, (b) the chemical potential μ/EF, (c) the scaled mean-photon number |α¯|2, and (d) the scaled polarization m¯ as functions of the effective atom-photon coupling strength η/EF. Analytical ...

Next, turning to the case μ=μχ¯0, we find from Eq. (29) that


Imposing the stability condition 2E¯G/(|α¯|)2>0, and keeping in mind μ ≥ 0, we find that in the regimes ω0 < EF and 0<η<ηc(2), the relevant ground-state solutions are


Note that for ω0 < EF and η>ηc(2), we would have μ < 0 [c.f. Fig. 2(b)], where the ground-state solutions are of the form of Eq. (35). Thus in contrast to the case ω0 ≥ EF, the ground state of the Fermi gas with ω0 < EF exhibits the N-II phase for 0<η<ηc(2) and the SR phase for η>ηc(2). We summarize results below


We see that a quantum phase transition from the N-II phase to the SR phase occurs at the critical point ηc(2), which is first order in nature and is accompanied by a sudden change in |α¯|2 [see also red-dashed curves in Fig. 3(c)].

Interestingly, at the critical regime η=ηc(2), the scaled ground-state energies of the N-II and SR phases become equal, meaning both phases coexist. We shall therefore call it the N-II-SR mixed phase. In order to characterize this mixed phase, let x 0 stand for the fraction of the N-II part in the mixed phase, which can take arbitrary value in the regime 0 ≤ x0 ≤ 1. We find (see Methods Section for detailed derivation)





Since both the SR and N-II phases exhibit m¯0 [see Eqs (35) and (28)], the former due to the the macroscopic collective excitation of fermions while the later induced by the effective Zeeman field, it is naturally expected that both mechanisms contribute to the nonzero polarization of the N-II-SR mixed phase [see Eq. (48)]. In addition, we find two first-order quantum phase transitions, from the N-II-SR phase to the N-II phase and to the SR phase, respectively. This is different from the boson counterpart, as well as the previously discussed N-I phase, where the corresponding transitions are second order9, 10.

Collecting all above results, we plot in Fig. 4 the phase diagram in the entire parameter regimes of ω0/EF and η/EF. As predicted previously, while the quantum phase transition from the N-I phase to the SR phase is second order, a first-order transition occurs between the N-II and SR phases due to coexistence of both phases at the critical line. In addition, the phase diagram exhibits a tri-critical point (see the red dotted curves), where the quantum phase transition changes its character from the first to the second order.

Figure 4
Phase diagram as a function of the effective resonant frequency ω0/EF and the effective fermion-photon coupling strength η/EF. The atom-number dependent cavity frequency ω and the cavity decay rate κ are the same as those in Fig.  ...

Parameter estimation and possible experimental observation

We now provide an estimation of the relevant parameters taking 40K atom as the example. For the fermionic 40K atoms with the Fermi energy EF ∼ 0.73 MHz, the ground states with 2S1/2 are represented by |↑right angle bracket = |F = 9/2, m F = 9/2right angle bracket and |↓right angle bracket = |F = 9/2, m F = 7/2right angle bracket. For the excited states with 2P1/2, we choose |1right angle bracket = |F = 9/2, m F = 7/2right angle bracket and |2right angle bracket = |F = 9/2, m F = 9/2right angle bracket, with F and mF labeling the total angular momentum and magnetic quantum numbers, respectively.

Considering the optical properties of the 40K D1-line, we take the cavity length as 178 μm and the wavelengths of the transverse pumping lasers as 770 nm. This ensures that both g 1 and g 2 characterizing the fermion-photon coupling strength are at the order of MHz, such that the rotating-wave approximation leading to Eqs (4) and (5) is valid. For a cavity with a waist radius 27 μm and a finesse ~105, we estimate its decay rate κ is at the order of MHz. For the effective fermion-photon coupling strength η=Ng1Ω1/(2Δ1) =Ng2Ω2/(2Δ2), thanks to the prefactor N, a magnitude of MHz can be achieved simply by modifications of the Rabi frequencies of the transverse pumping lasers, despite the large detuning required by the adiabatic approximation leading to Eq. (8). Finally, the effective resonant frequency ω 0 = (ω˜ω˜)/2 and the atom-number dependent cavity frequency ω=Nζ+ω˜ can be easily controlled by tuning the frequencies of the driving and transverse pumping lasers. In the experiments3, both ω 0 and ω can be tuned from − GHz to GHz and even beyond.

We conclude this section by briefly discussing how to probe the predicted quantum phases and phase diagrams. As has been elaborated, the physics of the ground-state phases is determined by the scaled mean-photon number |α¯|2 and the scaled polarization m¯. Building on the development of the state of the art experimental techniques, m¯ can be measured by observing the different density distributions between the two-component Fermi gas39, 40, while α can be detected using the calibrated single-photon counting modules which allows for the in situ monitor of the intra-cavity light intensity3. We thus expect that our results to be feasible within the experimental capabilities.


In realistic experiments, there typically exists a shallow harmonic external confining potential, which can be modeled by V(r)=ω2r2/2, where ω is the harmonic trap frequency and r is the harmonic trap radius. In order to estimate its effect, we first use the local density approximation41 to obtain an effective chemical potential μ(r) = μ0 − V(r), where μ 0 is the chemical potential at the center of the harmonic trap and μ(r) determines the total density nnn. Then, the density distribution n(r) and the magnetization distribution m(r) can be solved from42


In Fig. 5, we plot the density distributions n(r)/n 0 in the radial direction for (a) ω0E˜F and (b) ω0<E˜F, respectively, where n0=ME˜F/(2π), RT=2E˜F/ω is the Thomas-Fermi radius, and E˜F=Nω is the Fermi energy in the trapped systems. As shown in Fig. 5(a), for ω0E˜F (fully polarized fermions), the density profile does not depend on the effective fermion-photon coupling strength η, which implies that the critical line η=ηc(1) for the second-order quantum phase transition from the N-I phase to the SR phase is unaffected by the trapping potential. However, for ω0<E˜F (partially polarized fermions), the density profile depends strongly on the effective fermion-photon coupling strength η, leading to the modification of the critical line for the first-order quantum phase transition from the N-II phase to the SR phase, as shown in Fig. 5(b). Moreover, the predicted N-II-SR mixed phase only exists near the center of the trap and the corresponding density profile has thus a jump discontinuity, which provides an experimentally observable signature. We also note that the presence of the harmonic trap renders the system to become of finite size, i.e., the density profile vanishes for size r/RT > rc/RT, where r c denotes the size of the trap, as shown in both Fig. 5(a) and 5(b).

Figure 5
The density distributions n(r)/n0 for (a) ω0E~F and (b) ω0<E~F along the radial direction r/RT. The atom-number dependent cavity frequency ω and the cavity decay rate κ are the same as those in Fig.  ...

In addition, when the decay rate κ, which has the order of MHz in experiments, is much larger than the external atomic degrees of freedom, the relatively slow atoms at times t ≫ 1/κ feel an average affect of steady-state photons, i.e., Eq. (14). Furthermore, when the system enters the SR phase with a large mean-photon number, the commutator [aˆ,aˆ]=1 can be ignored at the leading order results. Thus, it is reasonable to neglect the non-commutability between aˆ and aˆ and just set them to their steady-state averaged value3, 4, 1225. Based on the above two arguments, the employed mean-field approximation, which assumes to approximate the photon annihilation operator by its steady-state averaged value, provides a good description of the system within our context.

In summary, we have analytically investigated the ground-state properties of a 2D polarized degenerate Fermi gas coupled to a high-finesse optical cavity. By solving the photon-number dependent BdG equation, we have found rich quantum phases and phase diagrams, which depend crucially on the fermion-photon coupling strength and the atomic resonant frequency (effective Zeeman field). In particular, for the weak atomic resonant frequency, we have shown that there exists a mixed phase where the N-II and SR phases coexist. In addition, we have revealed a first-order quantum phase transition from the N-II phase to the SR phase, which contrasts to the familiar second-order transition for the case with bosons. Finally, we have presented a parameter estimation and have addressed briefly how to detect these predicted quantum phases and phase diagrams in experiments.


In this Section, we present a detailed derivation of Eqs (45)–(48). We denote by μNIISR ∈ [μNIIμSR] the chemical potential in the N-II-SR mixed phase, and let nNII (n SR) and μNII (μ SR) denote the atom density and the chemical potential in the N-II (SR) phase, respectively. In terms of the fraction 0 ≤ x0 ≤ 1 of the N-II phase, we have42


When x 0 = 1, we simply have the N-II phase [see Eq. (28)], and nNIIMμNII/π. Instead, when x 0 = 0, E¯GSR(μSR,ω0,η), μ SR, |α¯|SR, n SR, and m¯SR are given by






where A = 1 − ωη2/[EF(ω2κ2)].

Using the phase equilibrium condition E¯GNII(μNIISR,ω0)= E¯GSR(μNIISR,ω0,η) 37, 42, together with Eqs (28) and (51), we see that when μNIISRμNIIEF corresponding to x 0 = 1, the phase boundary between the N-II phase and the N-II-SR mixed phase is described by η=ηc(2). On the other hand, when μNIISRμSR = 2EF − 2ωη2/(ω2κ2) corresponding to x 0 = 0, the phase boundary between the N-II-SR mixed phase and the SR phase is also described by η=ηc(2). In addition, for the N-II-SR mixed phase, we can determine μNIISR from the relation E¯GNII(μNIISR,ω0)= E¯GSR(μNIISR,ω0,η), which gives


This is the same as Eq. (46). Substituting Eq. (56), nNII, and n SR into Eq. (50), we find x 0 can take arbitrary values ranging from 0 to 1.

Finally, let us prove that the ground state of the N-II-SR mixed phase is stable. Due to the existence of both N-II and SR phases, the scaled ground-state energy in this mixed phase is defined by37, 42


Thus we obtain



which are either less than or equal to zero, i.e., the ground state of the N-II-SR mixed phase is stable at η=ηc(2).

Substituting Eqs (28), (51), and (56) into Eqs (53) and (57), we derive Eqs (45)–(47). In addition, according to Eq. (50), we obtain Eq. (48).


This work is supported in part by the National Key R&D Program of China under Grant No. 2017YFA0304203; the NSFC under Grants No. 11674200, No. 11422433, No. 11604392, No. 11434007, and No. 61378049; the PCSIRT under Grant No. IRT13076; the FANEDD under Grant No. 201316; SFSSSP; OYTPSP; and SSCC.

Author Contributions

Author Contributions

Y.F., K.Z., J.F., F.M., G.C., and S.J. conceived the idea, Y.F., K.Z., and J.F. performed the calculation, M.F., G.C., and S.J. wrote the manuscript, F.M. and G.C. supervised the whole research project.


Competing Interests

The authors declare that they have no competing interests.


Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Contributor Information

Feng Mei, moc.361@iemluar.

Gang Chen, moc.361@179gnagnehc.


1. Brennecke F, et al. Cavity QED with a Bose-Einstein condensate. Nature. 2007;450:268–271. doi: 10.1038/nature06120. [PubMed] [Cross Ref]
2. Colombe Y, et al. Strong atom-field coupling for Bose-Einstein condensates in an optical cavity on a chip. Nature. 2007;450:272–276. doi: 10.1038/nature06331. [PubMed] [Cross Ref]
3. Baumann K, Guerlin C, Brennecke F, Esslinger T. Dicke quantum phase transition with a superfluid gas in an optical cavity. Nature. 2010;464:1301–1306. doi: 10.1038/nature09009. [PubMed] [Cross Ref]
4. Baumann K, Mottl R, Brennecke F, Esslinger T. Exploring Symmetry Breaking at the Dicke Quantum Phase Transition. Phys. Rev. Lett. 2011;107 doi: 10.1103/PhysRevLett.107.140402. [PubMed] [Cross Ref]
5. Baden MP, Arnold KJ, Grimsmo AL, Parkins S, Barrett MD. Realization of the Dicke Model Using Cavity-Assisted Raman Transitions. Phys. Rev. Lett. 2014;113 doi: 10.1103/PhysRevLett.113.020408. [PubMed] [Cross Ref]
6. Landig R, Brennecke F, Mottl K, Donner T, Esslinger T. Measuring the dynamic structure factor of a quantum gas undergoing a structural phase transition. Nat. Commun. 2015;6 doi: 10.1038/ncomms8046. [PMC free article] [PubMed] [Cross Ref]
7. Klinder J, Keßler H, Reza Bakhtiari M, Thorwart M, Hemmerichl A. Observation of a Superradiant Mott Insulator in the Dicke-Hubbard Model. Phys. Rev. Lett. 2015;115 doi: 10.1103/PhysRevLett.115.230403. [PubMed] [Cross Ref]
8. Landig R, et al. Quantum phases from competing short- and long-range interactions in an optical lattice. Nature. 2016;532:476–479. doi: 10.1038/nature17409. [PubMed] [Cross Ref]
9. Wang YK, Hioe FT. Phase Transition in the Dicke Model of Superradiance. Phys. Rev. A. 1973;7:831–836. doi: 10.1103/PhysRevA.7.831. [Cross Ref]
10. Hioe FT. Phase Transitions in Some Generalized Dicke Models of Superradiance. Phys. Rev. A. 1973;8:1440–1445. doi: 10.1103/PhysRevA.8.1440. [Cross Ref]
11. Ritsch H, Domokos P, Brennecke F, Esslinger T. Cold atoms in cavity-generated dynamical optical potentials. Rev. Mod. Phys. 2013;85:553–601. doi: 10.1103/RevModPhys.85.553. [Cross Ref]
12. Kanamoto R, Meystre P. Optomechanics of a Quantum-Degenerate Fermi Gas. Phys. Rev. Lett. 2010;104 doi: 10.1103/PhysRevLett.104.063601. [PubMed] [Cross Ref]
13. Sun Q, Hu X-H, Ji A-C, Liu WM. Dynamics of a degenerate Fermi gas in a one-dimensional optical lattice coupled to a cavity. Phys. Rev. A. 2011;83 doi: 10.1103/PhysRevA.83.043606. [Cross Ref]
14. Müller M, Strack P, Sachdev S. Quantum charge glasses of itinerant fermions with cavity-mediated long-range interactions. Phys. Rev. A. 2012;86 doi: 10.1103/PhysRevA.86.023604. [Cross Ref]
15. Padhi B, Ghosh S. Cavity Optomechanics with Synthetic Landau Levels of Ultracold Fermi Gas. Phys. Rev. Lett. 2013;111 doi: 10.1103/PhysRevLett.111.043603. [PubMed] [Cross Ref]
16. Guo X, Ren Z, Guo G, Peng J. Ultracold Fermi gas in a single-mode cavity: Cavity-mediated interaction and BCS-BEC evolution. Phys. Rev. A. 2012;86 doi: 10.1103/PhysRevA.86.053605. [Cross Ref]
17. Keeling J, Bhaseen MJ, Simons BD. Fermionic Superradiance in a Transversely Pumped Optical Cavity. Phys. Rev. Lett. 2014;112 doi: 10.1103/PhysRevLett.112.143002. [PubMed] [Cross Ref]
18. Piazza F, Strack P. Umklapp Superradiance with a Collisionless Quantum Degenerate Fermi Gas. Phys. Rev. Lett. 2014;112 doi: 10.1103/PhysRevLett.112.143003. [PubMed] [Cross Ref]
19. Chen Y, Yu Z, Zhai H. Superradiance of Degenerate Fermi Gases in a Cavity. Phys. Rev. Lett. 2014;112 doi: 10.1103/PhysRevLett.112.143004. [PubMed] [Cross Ref]
20. Pan J-S, Liu X-J, Zhang W, Yi W, Guo G-C. Topological Superradiant phases in a Degenerate Fermi Gas. Phys. Rev. Lett. 2015;115 doi: 10.1103/PhysRevLett.115.045303. [PubMed] [Cross Ref]
21. Chen Y, Zhai H, Yu Z. Superradiant phase transition of Fermi gases in a cavity across a Feshbach resonance. Phys. Rev. A. 2015;91 doi: 10.1103/PhysRevA.91.021602. [Cross Ref]
22. Kollath C, Sheikhan A, Wolff S, Brennecke F. Ultracold Fermions in a Cavity-Induced Artificial Magnetic Field. Phys. Rev. Lett. 2016;116 doi: 10.1103/PhysRevLett.116.060401. [PubMed] [Cross Ref]
23. Sheikhan A, Brennecke F, Kollath C. Cavity-induced chiral phases of fermionic quantum gases. Phys. Rev. A. 2016;93 doi: 10.1103/PhysRevA.93.043609. [Cross Ref]
24. Sheikhan A, Brennecke F, Kollath C. Cavity-induced generation of non-trivial topological states in a two-dimensional Fermi gas. Phys. Rev. A. 2016;94 doi: 10.1103/PhysRevA.94.061603. [Cross Ref]
25. Zheng W, Cooper NR. Superradiance Induced Particle Flow via Dynamical Gauge Coupling. Phys. Rev. Lett. 2016;117 doi: 10.1103/PhysRevLett.117.175302. [PubMed] [Cross Ref]
26. Cooper LN. Bound Electron Pairs in a Degenerate Fermi Gas. Phys. Rev. 1956;104:1189–1190. doi: 10.1103/PhysRev.104.1189. [Cross Ref]
27. Deng Y, Cheng J, Jing H, Yi S. Bose-Einstein Condensates with Cavity-Mediated Spin-Orbit Coupling. Phys. Rev. Lett. 2014;112 doi: 10.1103/PhysRevLett.112.143007. [PubMed] [Cross Ref]
28. Dong L, Zhou L, Wu B, Ramachandhran B, Pu H. Cavity-assisted dynamical spin-orbit coupling in cold atoms. Phys. Rev. A. 2014;89 doi: 10.1103/PhysRevA.89.011602. [Cross Ref]
29. Dimer F, Estienne B, Parkins AS, Carmichael HJ. Proposed realization of the Dicke-model quantum phase transition in an optical cavity QED system. Phys. Rev. A. 2007;75 doi: 10.1103/PhysRevA.75.013804. [Cross Ref]
30. Fan J, et al. Hidden continuous symmetry and Nambu-Goldstone mode in a two-mode Dicke model. Phys. Rev. A. 2014;89 doi: 10.1103/PhysRevA.89.023812. [Cross Ref]
31. Dicke RH. Coherence in Spontaneous Radiation Processes. Phys. Rev. 1954;93:99–110. doi: 10.1103/PhysRev.93.99. [Cross Ref]
32. Ford GW, Lewis JT, O’Connell RF. Quantum Langevin equation. Phys. Rev. A. 1988;37:4419–4428. doi: 10.1103/PhysRevA.37.4419. [PubMed] [Cross Ref]
33. Scully, M. O. & Zubariry, M. S. Quantum Optics (Cambridge University, 1997).
34. Larson J, Morigi G, Lewenstein M. Cold Fermi atomic gases in a pumped optical resonator. Phys. Rev. A. 2008;78 doi: 10.1103/PhysRevA.78.023815. [Cross Ref]
35. Nagy D, Szirmai G, Domokos P. Self-organization of a Bose-Einstein condensate in an optical cavity. Eur. Phys. J. D. 2008;48:127–137. doi: 10.1140/epjd/e2008-00074-6. [Cross Ref]
36. Sheehy DE, Radzihovsky L. BEC-BCS Crossover in “ Magnetized” Feshbach-Resonantly Paired Superfluids. Phys. Rev. Lett. 2006;96 doi: 10.1103/PhysRevLett.96.060401. [PubMed] [Cross Ref]
37. Sheehy DE, Radzihovsky L. BEC–BCS crossover, phase transitions and phase separation in polarized resonantly-paired superfluids. Ann. Phys. 2007;322:1790–1924. doi: 10.1016/j.aop.2006.09.009. [Cross Ref]
38. Sheehy DE. Fulde-Ferrell-Larkin-Ovchinnikov phase of two-dimensional imbalanced Fermi gases. Phys. Rev. A. 2015;92 doi: 10.1103/PhysRevA.92.053631. [Cross Ref]
39. Partridge GB, Li W, Kamar RI, Liao Y-a, Hulet RG. Pairing and Phase Separation in a Polarized Fermi Gas. Science. 2006;311:503–505. doi: 10.1126/science.1122876. [PubMed] [Cross Ref]
40. Zwierlein MW, Schirotzek A, Schunck CH, Ketterle W. Fermionic Superfluidity with Imbalanced Spin Populations. Science. 2006;311:492–496. doi: 10.1126/science.1122318. [PubMed] [Cross Ref]
41. Kohn W, Sham LJ. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965;140:A1133–A1138. doi: 10.1103/PhysRev.140.A1133. [Cross Ref]
42. He L, Zhuang P. Phase diagram of a cold polarized Fermi gas in two dimensions. Phys. Rev. A. 2008;78 doi: 10.1103/PhysRevA.78.033613. [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group