QCD-like phase diagram with Efimov trimers and Cooper pairs in resonantly interacting SU(3) Fermi gases

TL;DR

This paper explores the phase diagram of resonantly interacting SU(3) Fermi gases, revealing the interplay of Efimov trimers and color superfluidity, with implications for simulating dense QCD-like matter.

Contribution

It introduces a finite-range interaction model to study Efimov trimers and color superfluidity, providing the first phase diagram analysis including these effects in SU(3) Fermi gases.

Findings

Identification of a phase diagram with trimer, normal, and superfluid phases

Critical temperature estimates for color superfluidity

Demonstration of QCD-like phase behavior in ultracold gases

Abstract

We investigate color superfluidity and trimer formation in resonantly interacting SU(3) Fermi gases with a finite interaction range. The finite range is crucial to avoid the Thomas collapse and treat the Efimov effect occurring in this system. Using the Skorniakov-Ter-Martirosian (STM) equation with medium effects, we show the effects of the atomic Fermi distribution on the Efimov trimer energy at finite temperature. We show the critical temperature of color superfluidity within the non-selfconsistent $T$-matrix approximation (TMA). In this way, we can provide a first insight into the phase diagram as a function of the temperature $T$ and the chemical potential $\mu$. This phase diagram consists of trimer, normal, and color-superfluid phases, and is similar to that of quantum chromodynamics (QCD) at finite density and temperature, indicating a possibility to simulate the properties of such extremely dense matter in the laboratory.

Equations

$$H$$

$$\frac{m}{4\pi a}=-\frac{\nu }{g^2}+\sum _{\mathbf{p}}\frac{m}{{\mathbf{p}}^2},$$

$$R_{\ast }=\frac{4\pi }{m^2{g}^2}.$$

$$G(p)=\frac{1}{i{\omega }_n-{\xi }_{\mathbf{p}}^{\mathrm{F}}-\Sigma (p)},$$

$$\Sigma (p)=2g^{2}T\sum _{Q}D(Q)G_0(Q-p),$$

$$D(Q)=\frac{1}{i{\nu }_{n^{\prime }}-{\xi }_{\mathbf{Q}}^{\mathrm{B}}-\Xi (Q)},$$

$$\Xi (Q)$$

$$=$$

$$n=3T\sum _{p}G(p)-6T\sum _{Q}D(Q).$$

$$\frac{m^2\mu R_{\ast }}{2\pi }+\sum _{\mathbf{p}}\left[\frac{1}{2{\xi }_{\mathbf{p}}^{\mathrm{F}}}\tanh \left(\frac{{\xi }_{\mathbf{p}}^{\mathrm{F}}}{2T_{\mathrm{c}}}\right)-\frac{m}{{\mathbf{p}}^2}\right]=0.$$

$$T_3^{\mathrm{M}}(p,p^{\prime };P)$$

$$\left[-R_{\ast }\kappa (\mathbf{p}{)}^2+4\pi \sum _{\mathbf{k}}\left\{\frac{F(\mathbf{k},\mathbf{p})}{{\mathbf{k}}^2+\kappa (\mathbf{p}{)}^2}-\frac{1}{{\mathbf{k}}^2}\right\}\right]L\left(\mathbf{p}\right)=-8\pi \sum _{\mathbf{q}}\frac{F(\mathbf{q},\mathbf{p})L(\mathbf{q}+\mathbf{p}/2)}{{\mathbf{q}}^2+\kappa (\mathbf{p}{)}^2},$$

$$F(\mathbf{p},\mathbf{q})=1-f({\xi }_{\mathbf{p}+\mathbf{q}/2}^{\mathrm{F}})-f({\xi }_{\mathbf{p}-\mathbf{q}/2}^{\mathrm{F}})$$

$$F(\mathbf{p},\mathbf{q})=\theta ({\xi }_{\mathbf{p}+\mathbf{q}/2}^{\mathrm{F}})\theta ({\xi }_{-\mathbf{p}+\mathbf{q}/2}^{\mathrm{F}}),$$

$$E_3^{\mathrm{M}}=\frac{1}{m{R}_{\ast }^2}X\left(R_{\ast }/{\lambda }_T,\mu /T\right),$$

$$A(\mathbf{p},\omega )=-\frac{1}{\pi }\mathrm{Im}G(\mathbf{p},i{\omega }_n\to \omega +i\delta )$$

$$\rho (\omega )=\sum _{\mathbf{p}}A(\mathbf{p},\omega ).$$

$$T_3^{\mathrm{M}}(p,p^{\prime };P)=-2g^{2}T\sum _q{G}_0(P-p-q)G_0(q)D(P-q)T_3^{\mathrm{M}}(q,p^{\prime };P),$$

$$T\sum _{i{\Omega }_n}=-{\oint }_C\frac{dx}{2\pi i}f(x),$$

$$T_3^{\mathrm{M}}(p,p^{\prime };P)$$

$$=$$

$$T_3^{\mathrm{M}}(p,p^{\prime };P)$$

$$D(\mathbf{q},{\xi }_{\mathbf{p}+\mathbf{q}}^{\mathrm{F}}+{\xi }_{\mathbf{p}}^{\mathrm{F}})\simeq D(\mathbf{q},E_3^{\mathrm{M}}-3\mu -{\xi }_{\mathbf{q}}^{\mathrm{F}}),$$

$$T_3^{\mathrm{M}}\left(\left\{\mathbf{q},E_3^{\mathrm{M}}-3\mu -{\xi }_{\mathbf{p}}^{\mathrm{F}}-{\xi }_{\mathbf{p}+\mathbf{q}}^{\mathrm{F}}\right\},p^{\prime };P\right)\simeq T_3^{\mathrm{M}}\left(\left\{\mathbf{q},{\xi }_{\mathbf{q}}^{\mathrm{F}}\right\},p^{\prime };P\right).$$

$$T_3^{\mathrm{M}}(p,p^{\prime };P)$$

$$D(\mathbf{q},E_3^{\mathrm{M}}-3\mu -{\xi }_{\mathbf{q}}^{\mathrm{F}})=\frac{4\pi }{m{g}^2}{\left[-R_{\ast }\kappa (\mathbf{q}{)}^2+4\pi \sum _{\mathbf{k}}\left\{\frac{F(\mathbf{k},\mathbf{q})}{{\mathbf{k}}^2+\kappa (\mathbf{q}{)}^2}-\frac{1}{{\mathbf{k}}^2}\right\}\right]}^{-1},$$

$$L(\mathbf{q})=\frac{m{g}^2}{4\pi }D(\mathbf{q},E_3^{\mathrm{M}}-3\mu -{\xi }_{\mathbf{q}}^{\mathrm{F}})T_3^{\mathrm{M}}(\{\mathbf{q},{\xi }_{\mathbf{q}}^{\mathrm{F}}\},p^{\prime };P).$$

$$\left[-R_{\ast }\kappa (\mathbf{p}{)}^2+4\pi \sum _{\mathbf{k}}\left\{\frac{F(\mathbf{k},\mathbf{p})}{{\mathbf{k}}^2+\kappa (\mathbf{p}{)}^2}-\frac{1}{{\mathbf{k}}^2}\right\}\right]L(\mathbf{p})$$

$$=8\pi \sum _{\mathbf{q}}\frac{1-f({\xi }_{\mathbf{p}+\mathbf{q}}^{\mathrm{F}})-f({\xi }_{\mathbf{q}}^{\mathrm{F}})}{m{E}_3^{\mathrm{M}}-{\mathbf{p}}^2-{\mathbf{q}}^2-\mathbf{p}\cdot \mathbf{q}}L(\mathbf{q}).$$

$$\left[-R_{\ast }\kappa (\mathbf{p}{)}^2+4\pi \sum _{\mathbf{k}}\left\{\frac{1}{{\mathbf{k}}^2+\kappa (\mathbf{p}{)}^2}-\frac{1}{{\mathbf{k}}^2}\right\}\right]L\left(\mathbf{p}\right)=-8\pi \sum _{\mathbf{q}}\frac{L(\mathbf{q}+\mathbf{p}/2)}{{\mathbf{q}}^2+\kappa (\mathbf{p}{)}^2},$$

Full text

I Introduction

The concept of quantum simulation has opened new possibilities for exploring the properties of novel materials as well as exotic matter in extreme conditions Bloch ; Georgescu ; Gross . Recently, ultracold atoms have been used as quantum simulators of strongly correlated systems thanks to the tunability of their physical parameters such as interparticle interaction Dalfovo ; Bloch2 ; Giorgini ; Chin . For example, ultracold atomic Fermi gases loaded into an optical lattice can realize the Hubbard model, which is relevant to high-$T_{\rm c}$ cuprate superconductors Gross ; Bloch3 ; Kohl ; Jordens ; Mazurenko . The antiferromagnetic behavior which plays an important role for the superconducting mechanism of high-$T_{\rm c}$ cuprates has been observed in this atomic system Gross ; Mazurenko . As another example, strongly interacting two-component Fermi gases can be used as quantum simulators of dilute neutron matter Gezerlis ; vanWyk ; Strinati . Indeed, the observed thermodynamic quantities of homogeneous Fermi gases near the unitarity limit Navon ; Horikoshi ; Tajima2 ; HorikoshiE quantitatively reproduce the equation of state of neutron matter obtained by numerical simulations APR which is crucial for understanding the interior of a neutron star Oertel ; Baym .

The analog simulation of quantum chromodynamics (QCD) Greiner ; Fukushima , where quarks with three colors strongly interact with each other, can be regarded as one of the next important challenges for cold-atom physics Maeda ; Cirac ; Ohara . While numerous efforts have been made to explore the phase diagram of finite-density QCD where various phenomena such as color superconductivity Alford have been theoretically proposed, they could not be confirmed by high-energy experiments due to the extreme densities, nor by exact numerical simulations due to the sign problem. On the other hand, three-component fermion systems have been experimentally realized by using mixtures of fermionic atoms in three different internal states Ohara ; Ottenstein ; Huckans ; Williams ; Spiegelhalder ; Lompe ; Wenz ; Nakajima1 ; Nakajima2 . These systems reach the quantum degenerate regime around $T\simeq 0.3T_{\rm F}$ Ottenstein where $T_{\rm F}$ is the Fermi degenerate temperature of non-interacting atoms. In this regime, the existence of color superfluidity similar to color superconductivity in QCD has been theoretically examined in three-component Fermi gases He ; Rapp ; Paananen ; Bedaque ; Kantian ; Catelani ; Martikainen ; Ozawa . Although there are several differences between these atomic systems and QCD, they constitute a good starting point to study the strong-coupling effects in three-component fermion systems.

In most of the previous works on color superfluidity He ; Rapp ; Paananen ; Bedaque ; Kantian ; Catelani ; Martikainen ; Ozawa , zero-range interactions have been used to describe the attraction between the three components. However, such zero-range interactions are known to lead to a collapse of the system Thomas ; Blume in connection with the existence of Efimov trimers Efimov1 ; Efimov2 ; Braaten ; Naidon2 ; Greene ; Williams ; Nakajima1 ; Nakajima2 . A finite range of interactions is therefore necessary to prevent the collapse and properly treat Efimov trimers. The effect of finite range and the Efimov trimers on the phase diagram of three-component fermionic system was studied in Ref. Nishida . However, this study was limited to zero temperature.

In this paper, we investigate in-medium Efimov trimers and color superfluidity in resonantly interacting SU(3) Fermi gases with finite interaction range and finite temperature. By including medium effects in the Skorniakov-Ter-Martirosian (STM) equation STM which exactly describes Efimov physics in a three-body problem Naidon3 , we discuss how the trimer energy is changed at finite temperature. To explore many-body physics related to color superfluidity, we employ the non-selfconsistent $T$-matrix approximation (TMA) Pieri ; Perali1 ; Perali2 ; Tsuchiya ; Tsuchiya2 ; Watanabe ; Perali3 ; Mueller1 ; Palestini ; Ota ; Tajima3 , which can correctly describe the Bardeen-Cooper-Schrieffer to Bose-Einstein-condensation (BCS-BEC) crossover Eagles ; Leggett ; NSR ; SadeMelo ; Ohashi ; Chen in a strongly interacting two-component 40K Regal and 6Li Zwierlein Fermi gas. Combining the results of these two approaches, we obtain the phase diagram with respect to the temperature $T$ and the chemical potential $\mu$.

This paper is organized as follows. In Sec. II, we explain the Hamiltonian for resonant SU(3) Fermi gases and our framework of in-medium STM equation and non-selfconsistent TMA. In Sec. III, we show our numerical results of the in-medium Efimov trimer energy and discuss the finite temperature behavior of them, as well as the color-superfluid instability. Finally, we summarize this paper in Sec. IV. For simplicity, we take $\hbar=k_{\rm B}=1$ and the system volume $V$ is taken to be unity.

II Formalism

II.1 Hamiltonian

We consider a symmetric three-component fermionic system. We model the interaction between two different components by a two-channel Feshbach resonance model Chin ; Nishida consisting of an open channel and a closed channel described by a single diatomic molecular state. We neglect the interaction within the open channel which corresponds to the limit of closed-channel dominated resonances. Furthermore, the coupling between open channel and molecular states is taken as a contact-type coupling. The three-component fermionic system is therefore described by the following Hamiltonian,

[TABLE]

where $\xi_{\bm{p}}^{\rm F}=\frac{\bm{p}^{2}}{2m}-\mu$ and $\xi_{\bm{q}}^{\rm B}=\frac{\bm{q}^{2}}{4m}+\nu-2\mu$ are the kinetic energies of a fermion with mass $m$ and the diatomic molecule, respectively, and $\mu$ is the fermionic chemical potential. $c_{\bm{p},j}$ and $b_{\bm{q},ij}$ are annihilation operators of a Fermi atom with the internal state $j=1,2,3$ and a diatomic molecule made of $i-j$ atomic pair. The energy of diatomic molecules $\nu$ and the atom-dimer Feshbach coupling $g$ are related to the scattering length $a$ and the range parameter $R_{*}$ as follows,

[TABLE]

[TABLE]

We note that the effective range $r_{\rm e}$ of this two-channel interaction is here negative and is associated with $R_{*}=-\frac{1}{2}r_{\rm e}$. In this paper, we focus on the unitarity limit $1/a=0$.

II.2 non-selfconsistent $T$-matrix approximation

We first employ the non-selfconsistent $T$-matrix approximation within the Matsubara formalism Pieri ; Perali1 ; Perali2 ; Tsuchiya ; Tsuchiya2 ; Watanabe ; Perali3 ; Mueller1 ; Palestini ; Ota ; Tajima3 to determine thermodynamic properties such as critical temperature $T_{\rm c}$ of color superfluidity. In this framework, the thermal Green’s function of dressed atoms is given by

[TABLE]

where we used the four-momentum notation $p=(\bm{p},i\omega_{n})$ and $\omega_{n}=(2n+1)\pi T$ is the fermionic Matsubara frequency. $\Sigma(p)$ is the fermionic self-energy. In the case of just two fermions in vacuum, the self-energy is given by the diagram shown in Fig. 1(a). Retaining only this diagram, we obtain Tajima3

[TABLE]

where $Q=(\bm{Q},i\nu_{n^{\prime}})$ and $\nu_{n^{\prime}}=2\pi n^{\prime}T$ is the bosonic Matsubara frequency. $G_{0}(p)=1/\left(i\omega_{n}-\xi_{\bm{p}}^{\rm F}\right)$ and $D(Q)$ are the in-medium Green’s functions of a non-interacting fermion and a dressed molecule, respectively. We note that although this equation (6) as the same form as that in vacuum, here $G_{0}$ and $D$ contain the medium effects. We also note that the factor $2$ in Eq. (6) comes from the degree of freedom with respect to internal states.

Similarly, the thermal Green’s function of dressed molecules with the ultraviolet renormalization is given by

[TABLE]

where $\Xi(Q)$ is the bosonic self-energy diagrammatically shown in Fig. 1(b). Here again, we take the vacuum form

[TABLE]

where $f(\xi)=1/(e^{\xi/T}+1)$ is the Fermi-Dirac distribution function. The chemical potential $\mu/\varepsilon_{\rm F}$ with fixed number density $n$ [where $\varepsilon_{\rm F}=(2\pi^{2}n)^{2/3}/(2m)$ is the Fermi energy of an ideal Fermi gas] is obtained by solving the number equation,

[TABLE]

In addition, we obtain the critical temperature $T_{\rm c}$ from the Thouless criterion $D^{-1}(\bm{Q}=\bm{0},i\nu_{n}=0)=0$ Thouless , which gives ($1/a=0$)

[TABLE]

II.3 in-medium Skorniakov-Ter-Martirosian equation

To determine the trimer energy $E_{3}^{\rm M}$ in the medium, we consider the three-body $T$-matrix equation Brodsky ; Iskin which is diagrammatically shown in Fig. 2. In the vacuum case, it gives the so-called Skorniakov-Ter-Martirosian (STM) equation which exactly describes Efimov physics in three-body problems. Using the thermal Green’s functions within the Matsubara formalism, we obtain the in-medium three-body $T$-matrix

[TABLE]

We note that $p$, $p^{\prime}$, and $P$ are the four-momenta of incoming fermion, outcoming fermion, and the center-of-the-mass, respectively (we suppress the index of internal states for simplicity). The factor $2$ in the second term of r. h. s. of Eq. (13) comes from the degree of freedom with respect to internal states. $G_{0}(q)$ with $q=(\bm{q},i\Omega_{n})$ indicates intermediate states of fermions which are integrated. By including Pauli-blocking effects on $G_{0}(p)$ (see Appendix A), we obtain the in-medium STM equation in the unitarity limit $(1/a\rightarrow 0)$,

[TABLE]

where $\kappa(\bm{q})^{2}=\frac{3}{4}\bm{q}^{2}-mE_{3}^{\rm M}$, $L(\bm{p})$ is the function defined by Eq. (33) in Appendix A, and

[TABLE]

is the statistical factor associated with the Fermi-Dirac distribution of atoms given by $f(\xi_{\bm{p}}^{\rm F})=1/\left(e^{\xi_{\bm{p}}^{\rm F}/T}+1\right)$. One can find that the ordinary STM equation in a three-body problem is recovered by setting $F(\bm{p},\bm{q})=1$. We note that Eq. (15) is qualitatively consistent with the previous works Niemann ; Tajima1 . We note however that Ref. Niemann is restricted to $T=0$ with the choice

[TABLE]

instead of Eq. (15). Here $\theta(x)$ denotes the Heaviside step function. This choice allows the possibility of positive-energy solution of STM equation, called Cooper triples. In contrast, our result Eq. (15) (see Appendix A) does not allow such solutions and we focus on the region where $E_{3}^{\rm M}\leq 0$. We briefly note that similar in-medium three-body equations were employed in nuclear physics Beyer2 ; Beyer ; Barbieri ; Beyer3 ; Kuhrts ; Pepin ; Kvinikhidze ; Mattiello .

III Results

Figure 3 shows the ground-state trimer energy $E_{3}^{\rm M}$ with the medium effects, which can be obtained numerically from Eq. (14) as a dimensionless function

[TABLE]

where $\lambda_{T}=\sqrt{\frac{2\pi}{mT}}$ is the thermal de Broglie wavelength. Physically, the range parameter $R_{*}$ gives the typical size of the Efimov trimer Petrov . In this regard, the ratio between $R_{*}$ and $\lambda_{T}$ represents how trimer states are affected by finite temperature effects. In the case of $R_{*}\ll\lambda_{T}$ and \mu/T\ \raise 1.29167pt\hbox{<}\kern-8.00003pt\lower 3.01385pt\hbox{\sim}\ 0, $E_{3}^{\rm M}$ approaches the vacuum limit given by $E_{3}^{\rm V}=-0.01385/(mR_{*}^{2})$ which is close to that of a universal trimer Gogolin , since the trimer size is small enough compared to the typical thermal length scale. In addition, the ratio $\mu/T$ is associated with the fugacity $z=e^{\mu/T}$ and represents Pauli-blocking effects due to the atomic Fermi distribution, which plays a significant role when $\mu>0$. The absolute value of $E_{3}^{\rm M}$ is greatly reduced in the Fermi degenerate region. Finally, $E_{3}^{\rm M}$ disappears in the region where $\mu/T$ and $R_{*}/\lambda_{T}$ are relatively large. The physical interpretation of these effects is that Fermi atoms from the medium weaken the Efimov attraction between three atoms forming a trimer state. This phenomenon is somewhat similar to the Gor’kov-Melik-Barkhudarov (GMB) correction in weak-coupling superconductors for which the size of Cooper pairs is large GMB and the pairing interaction is screened by the medium Floerchinger ; Yu ; Ruan ; Pisani1 ; Pisani2 .

To see how these effects would appear in actual experiments at given temperatures and densities, we plot in Fig. 4 (a) the typical temperature dependence of trimer energy $E_{3}^{\rm M}$ at different range parameters according to the density equation of state obtained from the non-selfconsistent $T$-matrix approximation (TMA). By numerically solving the number equation Eq. (10), we obtain the temperature-dependent chemical potential $\mu$ as shown in Fig. 4 (b). We then use $\mu$ as an input for in-medium STM equation given by Eq. (14). This gives $E_{3}^{\rm M}$ showing in Fig. 4 (a). $E_{3}^{\rm M}$ has a peak structure which can be understood from two effects, that is, the evolution of the Fermi chemical potential and the decrease of $\lambda_{T}$ compared to the trimer size. In the high-temperature limit, we can neglect the interaction effects and $\mu$ is given by the number equation of ideal gases $n=3\sum_{\bm{p}}e^{-\xi_{\bm{p}}^{\rm F}/T}$. This gives approximately $\mu\simeq-\frac{3}{2}T{\rm ln}T$ and for large temperatures $E_{3}^{\rm M}$ reproduces the vacuum result $E_{3}^{\rm V}$ due to the large negative $\mu$. Decreasing the temperature makes the density and temperature effects more visible which weakens the trimer ($E_{\rm 3}^{\rm M}$ increases). However, at very low temperature $\mu$ becomes almost constant while $\lambda_{T}$ increases as $1/\sqrt{T}$ with decreasing temperature. As a result, it becomes larger and larger compared to trimer size $R_{*}$ (which is fixed in this figure), which suppresses the temperature effects and the trimer strengthens ($E_{\rm 3}^{\rm M}$ decreases). This decrease becomes sharper with increasing $R_{*}$ as shown in Fig. 4 (a). We note that our calculation is stopped at $T=T_{\rm c}$, where Eq. (10) is invalid below $T_{\rm c}$ due to the existence of the superfluid gap.

From these results, we get a first insight of the phase diagram with respect to the chemical potential $\mu$ and temperature $T$, as shown in Fig. 5. We expect a trimer phase (TP) where all of the atoms are bound into ground-state Efimov trimers, a color superfluid phase (CSF) where all three kinds of pairs are condensed, and a normal phase (NP) where the atoms form neither trimers nor condensed pairs. The boundary between TP and NP is estimated by the curve where $E_{3}^{\rm M}=0$. Threshold curve $E_{3}^{\rm M}=0$ does not mean any phase transition between trimer and normal phases but it should be a good indication of how the trimer character disappears in the high-density region of this system. CSF is defined by the region below $T_{\rm c}$. We note that $T_{\rm c}$ approaches [math] in low-density (zero-range) limit at $\mu mR_{*}^{2}=0$ Tajima1 . We also note that the high density limit of the critical temperature $T_{\rm c}^{\rm HDL}=0.133T_{\rm F}$ is simply obtained from the Bose-Einstein condensation (BEC) temperature of diatomic molecules in the presence of thermal-excited fermions and $\mu$ approaches [math] in this limit Tajima1 . This indicates that the system undergoes a crossover from unitary Cooper pairs to BEC of closed channel molecules. This behavior is specific to the narrow-resonance two-channel model used in this work. Interestingly, Fig. 6 is similar to the phenomenological phase diagram of QCD consisting of the hadron phase (analogue of TP), the color superconducting phase (analogue of CSF), and the deconfined quark phase (analogue of NP) Fukushima . However, while the phase transition at $T=T_{\rm c}$ between CSF and NP is of the second order, that of color superconductivity is of the first order due to the gauge coupling Alford . Moreover, the conjectured BEC-BCS crossover AbukiNishida in QCD with increasing the chemical potential is opposite to the BCS-BEC crossover found in this model at $1/a=0$. We stress again this is a particularity of the narrow-resonance two-channel model.

Although our in-medium STM equation cannot be justified near $T=T_{\rm c}$ where the Bose-Einstein distribution of diatomic molecules increases at the low energy region, one can find a quantum-phase-transition-like behavior around $\mu mR_{*}^{2}=0$. The inset of Fig. 5 shows the magnification around $\mu mR_{*}^{2}=0$, where the two curves we calculated cross each other. In reality, the region near this point is expected to be dominated by strong multi-body correlation due to the competition between trimer formation and color superfluidity Pepin , which cannot be captured by our treatment.

Finally, we look at pairing fluctuations above $T_{\rm c}$. Indeed, it is known that pairing fluctuations become strong near $T_{\rm c}$ in two-component Fermi gases near the unitarity limit.

Figure 6 shows the single-particle spectral function

[TABLE]

obtained from the analytic continuation of $G(p)$ given by Eq. (5) ($\delta$ is an infinitesimally small positive number), at $\mu mR_{*}^{2}=0.0298$ and $TmR_{*}^{2}=0.0302$ which is just above $T_{\rm c}$ indicated in the inset of Fig. 5. One can see that the atomic dispersion has a gap structure near $\omega=0$ even in the absence of the superfluid gap. This excitation gap in the normal phase originates from strong pairing fluctuations (preformed Cooper pair) and is called pseudogap, which has been extensively discussed for various strongly correlated quantum systems such as high-$T_{\rm c}$ superconductors PRL80.149 ; RMP79.353 , ultracold Fermi gases Pieri ; Perali1 ; Perali2 ; Tsuchiya ; Tsuchiya2 ; Watanabe ; Perali3 ; Mueller1 ; Palestini ; Ota ; Tajima3 ; Tajima4 ; Mueller2 ; Stewart ; Gaebler ; Feld ; Magierski ; Wlazlowski ; Sagi , color superconductivity Kitazawa ; He2 , and nuclear matter Schnell ; Bozek ; Huang . This single-particle excitation property is accessible by photo-emission spectrum measurement in cold atom systems Stewart ; Gaebler ; Sagi ; Feld . In principle, such experiments could also observe many-body effects associated with in-medium Efimov trimers. However, treating such many-body effects theoretically would require a self-consistent approach including both two-body and three-body correlations in the self-energy.

The pseudogap can also be seen in the single-particle density of states $\rho(\omega)$, which is defined by

[TABLE]

It is shown in Fig. 7 (a) at $T=T_{\rm c}$. This quantity clearly shows the pseudogap effect as a dip structure around $\omega=0$. The pseudogap disappears away from $\mu=T=0$, that is, in the high-density (large range parameter) limit as found in the case of two-component Fermi gases with negative effective range Tajima3 . To characterize this many-body phenomenon, we introduce the pseudogap size $\Delta_{\rm pg}$ defined as the half width of the dip Watanabe [see the inset of Fig. 7 (b)]. One can find that $\Delta_{\rm pg}$ grows when $\mu mR_{*}^{2}$ approaches [math] and reaches a maximum value $\Delta_{\rm pg}\simeq 0.55\varepsilon_{\rm F}$ at $\mu mR_{*}^{2}=0$. This enhancement of $\Delta_{\rm pg}$ indicates that many-body effects associated with pairing fluctuations are important in the region around $\mu=T=0$ in the phase diagram of Fig. 5. This confirms the expected competition between formation of Cooper pairs and that of Efimov trimers. This competition would occur around $E_{3}^{\rm M}=0$, which is shown as the vertical dashed line in Fig. 7(b).

IV Summary

In this paper, we have investigated some of the strong-coupling effects occurring in resonantly interacting SU(3) Fermi gases with a finite interaction range, namely, the in-medium Efimov trimer and the critical temperature $T_{\rm c}$ of color superfluidity.

The trimer formation is weakened by the medium effects, which consist of thermal agitation and Fermi pressure due to Pauli exclusion. The trimer is affected by thermal agitation when the thermal de Broglie wavelength is comparable to the trimer size given by the range of interaction. The Pauli-blocking effects are significant when the chemical potential $\mu$ becomes positive. As in the case of the Gor’kov-Melik-Barkhudarov corrections in weakly coupled superconductors where the pairing interaction forming loosely bound Cooper pairs is screened by electrons from the medium, the medium effects in three-component Fermi gases become stronger when the trimer size or the atomic Fermi sphere becomes larger. We have shown how these effects would appear in actual experiments varying temperature at fixed number density.

Finally, we have investigated the phase diagram with respect to the chemical potential $\mu$ and temperature $T$. Our calculations indicate the existence of three phases: trimers phase, normal phase, and color superfluidity. Interestingly, the obtained phase diagram is analogous to the phenomenological QCD phase diagram which consists of hadron, deconfined quark, and color superconducting phases. We emphasize that the finite interaction range plays an important role to obtain such a QCD-like phase diagram in this atomic system.

Near $\mu=T=0$ in our phase diagram, the system is expected to be dominated by strong two-body and three-body correlations resulting from the competition between trimer formation and color superfluidity. This idea is supported by our calculation of the in-medium trimer energy and the single-particle spectral function which exhibits strong pairing fluctuations near the color superfluid phase transition. A self-consistent treatment of two-body and three-body correlations is required to understand this interesting regime, which is left as a future problem. Our analysis does not exclude the possibility of trimer superfluidity due to a possible residual attraction between the trimers Endo3 . Such a state would be similar to the $p$-wave superfluidity in a Bose-Fermi mixture Kinnunen3 . Although we consider a resonant interaction in this paper, the phase diagram would quantitatively change when tuning the scattering length. In particular, the color superfluid phase (or molecular BEC) would start at a lower chemical potential in the case of a finite scattering length Nishida .

Acknowledgements.

H. T. thanks T. Hatsuda, Y. Nishida, and G. Baym for useful discussion. H. T. was supported by a Grant-in-Aid for JSPS fellows (No.17J03975). P. N. was supported by the RIKEN Incentive Research Project. This work was partially supported by RIKEN iTHEMS program.

Appendix A Derivation of the in-medium STM equation

Since we calculate the ground-state trimer energy $E_{3}^{\rm M}$ with the medium effects, we set $P=(\bm{0},i\zeta_{l})$ in the three-body $T$-matrix equation given by Eq. (13) Brodsky ; Iskin [where $i\zeta_{l}=(2l+1)\pi T$ is the fermionic Matsubara frequency]. We obtain

[TABLE]

where we ignore the first term of r. h. s. of Eq. (13) which is negligible near the pole of $T_{3}^{\rm M}$. The summation of the Matsubara frequency $i\Omega_{n}$ [where $q=(\bm{q},i\Omega_{n})$] in Eq. (20) can be replaced by the contour integral with respect to an anticlockwise path $C$ FW enclosing the pole of the Fermi-Dirac distribution function $f(x)$, namely, $x=i\Omega_{n}$ as

[TABLE]

since ${\rm Res}f(x=i\Omega_{n})=-T$. We note that $C$ can be deformed to a clockwise path $C^{\prime}$ which encloses the poles of $G_{0}$, $D$, and $T_{3}^{\rm M}$, which give medium effects associated with the momentum distributions of atoms, molecules, and trimers, respectively. For simplicity, we consider only the pole of $G_{0}$ to incorporate the effects of the atomic Fermi-Dirac distribution. This approximation is justified in the high temperature regime where the fugacity $z=e^{\mu/T}$ is small. In this regime, the atomic Fermi distribution function is approximately given by $f(\xi_{\bm{p}}^{\rm F})\simeq e^{-\frac{p^{2}}{2mT}}z$, whereas the molecular and trimer distribution functions are $b(\xi_{\bm{Q}}^{\rm B})\simeq e^{-\frac{\bm{Q}^{2}}{4mT}}z^{2}$ and $f(\xi_{\bm{P}}^{\rm T})\simeq e^{-\left(\frac{\bm{P}^{2}}{6m}+E_{3}^{\rm M}\right)/T}z^{3}$, respectively [where $b(\xi)=1/(e^{\xi/T}-1)$ is the Bose-Einstein distribution function and $\xi_{\bm{P}}^{\rm T}=\bm{P}^{2}/(6m)-3\mu+E_{3}^{\rm M}$ is the kinetic energy of a trimer]. Using this approximation, we can analytically perform the energy integration as

[TABLE]

where $q_{1}=(\bm{q},i\zeta_{l}-i\omega_{n}-\xi_{\bm{p}+\bm{q}}^{\rm F})$ and $q_{2}=(\bm{q},\xi_{\bm{q}}^{\rm F})$. To obtain the in-medium STM equation, we perform the analytic continuations $i\omega_{n}\rightarrow\xi_{\bm{p}}^{\rm F}$ and $i\zeta_{l}\rightarrow E_{3}^{\rm M}-3\mu$ in Eq. (22). In this way, we obtain

[TABLE]

Furthermore, the first integrand in Eq. (27) gives dominant contributions near its pole, corresponding to $E_{3}^{\rm M}-3\mu-\xi_{\bm{p}}^{\rm F}-\xi_{\bm{q}}^{\rm F}-\xi_{\bm{q}+\bm{p}}^{\rm F}=0$. We can then approximate the arguments of $D$ and $T_{3}^{\rm M}$ as

[TABLE]

and

[TABLE]

By substituting Eqs. (28) and (29) into Eq. (27), we obtain

[TABLE]

where

[TABLE]

is obtained from the analytic continuation of Eq. (7). We note that $\kappa(\bm{q})^{2}=\frac{3}{4}\bm{q}^{2}-mE_{3}^{\rm M}$ and $F(\bm{k},\bm{q})$ is the statistical factor defined by Eq. (15). We introduce

[TABLE]

We note that although $L(\bm{q})$ implicitly depends on $p^{\prime}$ and $P$, they do not change the equation of $E_{3}^{\rm M}$. By using $L(\bm{q})$, Eq. (31) can be rewritten by

[TABLE]

Finally, we obtain the in-medium STM equation, that is, Eq. (14) by making the substitutions $\bm{q}\rightarrow\bm{q}+\bm{p}/2$ and $\bm{p}\rightarrow-\bm{p}$ in Eq. (34). We note that by taking the vacuum limit $\mu\rightarrow-\infty$ where $f(\xi_{\bm{p}}^{\rm F})\rightarrow 0$, Eq. (14) reproduces the ordinary STM equation of a three-body problem at $1/a=0$ Nishida given by

[TABLE]

which gives the ground-state trimer energy $E_{3}^{\rm V}=-0.01385/(mR_{*}^{2})$ in vacuum.