# Dynamical Cluster Quantum Monte Carlo Study of the Single Particle Spectra of Strongly Interacting Fermion Gases

## Abstract

We study the single-particle spectral function of resonantly-interacting fermions in the unitary regime, as described by the three-dimensional attractive Hubbard model in the dilute limit. Our approach, based on the Dynamical Cluster Approximation and the Maximum Entropy Method, shows the emergence of a gap with decreasing temperature, as reported in recent cold-atom photoemission experiments, for coupling values that span the BEC-BCS crossover. By comparing the behavior of the spectral function to that of the imaginary time dynamical pairing susceptibility, we attribute the development of the gap to the formation of local bound atom pairs.

###### pacs:

One of the most exciting recent developments in correlated systems has been the observation of pairing and superfluidity of ultracold atomic fermions interacting via an s-wave Feshbach resonance (1); (2); (3). The precise tunability of the Feshbach resonance mediated interaction allows the exploration of superfluidity for coupling extending from the BCS regime, with loosely-bound Cooper pairs, to the BEC regime, with local pairs at strong attraction. Thus such experiments provide a highly controllable setting to study phenomena that appears in a range of other systems, including neutron stars (4), quark matter (5), and correlated electronic systems such as the high temperature superconductors (6).

While the low-temperature properties of ultracold fermion gases are relatively well understood and essentially described by the BEC-BCS crossover wavefunction (2), the finite temperature behavior still presents several mysteries, particularly in the unitary regime, where the s-wave scattering length between the two species of fermion diverges. At strong couplings preformed pairs are predicted to exist above the superfluid transition temperature (6), possibly yielding a pseudogap in the density of states (7). But, how the pseudogap phenomena is reflected in experimental probes, such as radio frequency spectroscopy (8); (9), is still an open question.

To shed light on the spectral properties of resonantly-interacting fermion gases, we compute the single particle spectral function as a function of frequency and wavevector . The spectral function has recently been measured in experiments with trapped K atoms via a novel momentum-resolved generalization (10) of radio frequency spectroscopy (8); (9). Various analytical methods have been developed to compute in the BEC-BCS crossover based on the ladder approximation (11); (12), self-consistent many-body theories (13); (15); (14) and large- expansions (16). These approaches are hampered by the lack of a small parameter in the unitary regime.

In this letter we combine the Dynamical Cluster Approximation (17); (18) and the Maximum Entropy Method (19) to obtain the spectral function in the thermodynamic limit without uncontrollable approximations (see also Ref. (20)). We observe a strong suppression of low-energy spectral weight in conjunction with the expected onset of pairing correlations. Our simulations also capture the concomitant formation of a double-peak structure in with decreasing temperature in the strong-coupling regime. At unitarity, we find a momentum distribution consistent with the predictions of Tan (21), and an estimate of the universal parameters and characterizing the unitary Fermi gas (22).

The Dynamical Cluster Approximation (DCA), unlike conventional QMC, allows the extraction of the single-particle spectral function in the thermodynamic limit. One disadvantage of the DCA in the present setting is that it applies to lattice models, while cold atoms possess the itinerant-fermion dispersion . However, we address this by focusing on the dilute limit of the three-dimensional cubic Hubbard model

(1) |

where annihilates a spin- fermion at site , the angle brackets restrict the sum over nearest neighbor sites with the matrix element , and is the attractive on-site interaction. Here, the chemical potential is tuned to keep the filling at the small value , where is half filling. The Fermi energy associated with this filling is . In this dilute limit, the hopping dispersion can be approximated by its low expansion , thus describing free fermions if we take .

By tuning close to , the interaction strength where a two-body bound state forms in vacuum (23), we can capture the unitary regime in cold Fermi gases. We expect that interaction effects are insensitive to the underlying lattice when the scattering length

(2) |

exceeds the lattice spacing. We find, as discussed below, that Eq. (1) works well in the BCS and unitary regime, but begins to deviate from the free-fermion results in the deep-BEC regime of tightly-bound local pairs where lattice effects cannot be neglected.

The DCA covers the lattice with clusters of sites using periodic boundary conditions. Correlations up to length are considered exactly while long range correlations are treated at the mean-field level (17); (18). We use the Hirsch-Fye quantum Monte Carlo (QMC) algorithm (24) as the quantum solver (25), and the maximum entropy method (19) to analytically continue the imaginary-time spectra to real frequencies. In the DCA the momentum resolution of the self energy is determined by a set of -points corresponding to the cluster reciprocal space(18). We use a Fourier Transform interpolation scheme to approximate on other k-points in the Brillouin zone. We use a Betts cluster(26) with . Our Hirsch-Fye QMC algorithm is efficient as long as is below the bandwidth , allowing us to study the BEC-BCS crossover but not the asymptotic BEC limit.

We have obtained the single-particle spectral function , with the retarded Green’s function, across the BEC-BCS crossover and for a wide range of temperatures. Our displayed results are along the line connecting the center of the first Brillouin zone to the corner , with referring to the wavevector magnitude. The spectral function along other directions is qualitatively equivalent. We focus on the large scattering length limit with , where we approximate the Fermi wavelength by the itinerant-fermion formula . Our simulations are in the normal state, meaning we cannot fully describe the broken symmetry state below the superfluid transition at . Above , however, expect our results to be quantitatively accurate. Previous work (23) has demonstrated that Eq. (1) yields reliable results for the transition temperature of itinerant superfluid Fermi gases. By taking the extrapolated transition temperature from Ref. (23) we get at unitarity coupling and filling .

The most striking feature of our results, shown in Fig. 1, is the strong suppression of the spectral weight near zero frequency with increasing attraction towards the BEC limit. Fig. 1 displays results for four coupling values at the temperature . In all plots, the behavior at high- is an artifact of the underlying lattice that we shall ignore; this is reasonable since experiments only probe the occupied part of the spectral function. Fig. 1a, in the weak-coupling BCS regime , shows an essentially free-particle dispersion for and . The narrowing of the linewidth as the Fermi surface is approached from is clearly seen, consistent with the system being a Fermi liquid in this regime. With increasing coupling a two branch structure begins to form at the Fermi surface. While this behavior is barely visible at the unitary point, , shown in Fig. 1 (b), it increases in magnitude with increasing interaction strength, as seen in panels (c) and (d). In Fig. 1(d), the deep BEC limit, the two branches of the spectral function are clearly visible, and qualitatively consistent with the recent cold-atom photoemission experiments (10).

To examine the onset of the pseudogap with decreasing temperature, we plot at fixed as a function of energy in Fig. 2. For each coupling we have chosen the value of momentum for which the distance between the two spectral branches is the shortest. We find the pseudogap, with decreasing , only for coupling values beyond the unitary point. In the BEC regime, as seen in Fig. 2, the high temperature spectral functions show a Gaussian shape at high temperatures that develops a pseudogap at lower temperatures. The pseudogap separates two asymmetric peaks, reminiscent of recent analytic results (12). We define the pseudogap as one-half the minimum peak-to-peak distance between the two spectral function peaks; our values for at the lowest temperature (characterizing the evolution of the pseudogap across the BEC-BCS crossover) are plotted in Fig. 3(b).

Even though our simulation does not explore the superfluid state, it is able to capture the onset of pairing correlations and estimate the universal parameters and characterizing the unitary Fermi gas (22). For example, the internal energy per particle can be expressed as

(3) |

near unitarity (22); (27). We quantify from our extracted values of , that is adjusted to always keep filling . Fig. 3(a) displays as a function of coupling. The chemical potential is positive and close to in the BCS regime and reduces monotonically with increasing attraction towards the BEC limit of local pairing. We find at unitarity, yielding our prediction for the universal parameter , which is consistent with recent experimental and theoretical results, see e.g. Refs. (2); (28).

We also calculate the momentum distribution

(4) |

with the Fermi function. The momentum distribution has been previously measured in the BEC-BCS crossover (29). Fig. 3(c) displays our results for showing the expected behavior: a sharp step in in the BCS regime that broadens with increasing attraction towards the BEC limit. To gain additional insight, we fit these curves to to test whether the Tan relation (21)

(5) |

governing the large behavior of the momentum occupation, holds for our system. Although our system mimics a free-fermion dispersion at small momenta, at large lattice effects become inescapable, and it is not clear if the Tan relation will be valid since it probes the extremely short distance behavior of interacting fermions. Despite this, we find that, at unitarity, the optimal fitting parameters are , and . The and are reasonably small and the fit, shown as a solid curve for the unitary plot in Fig. 3(c), is seen to be quite accurate, showing that our system does exhibit the universal properties of a unitary Fermi gas. From the theoretical relation , with governing the interaction dependence of the energy per particle near unitary (21) (Eq. (3)), we extract the prediction . In the BCS and BEC regimes, the best-fit values of and are larger, so that we do not find a good fit to the Tan formula. In particular in the deep BEC region, one expects significant occupation of high momentum states. While this certainly occurs in our system, the finite size of our momentum-space Brillouin zone restricts the possible occupied momenta, preventing our system from truly accessing the deep BEC limit.

It is natural to assume that the pseudogap observed in our single-particle spectra arises from the formation of s-wave pairing correlations. To test this assumption, we study the local s-wave pair susceptibility at spatial distance and imaginary time . Setting and , we have

(6) |

where is the s-wave pair operator and the average occurs over system sites. This quantity, , probes local pairing correlations with the longest distance along the imaginary time direction.

In Fig. 4, we plot as function of temperature (left panel) and coupling (right panel). On panel (a), we can see that vanishes when at weak coupling. This signals that, in the weak coupling region, atoms can not form stable local pairs; although, of course, Cooper pairing correlations will set in at very much lower temperatures, outside the range of our computational results. When is large enough, towards the BEC limit, shows a finite -dependent intercept in the zero temperature limit, showing the establishment of local pairs. Fig. 4(b) shows the coupling dependence of at low temperatures. While this quantity is small at weak coupling in the BCS regime, near unitarity it undergoes a rapid increase as local pairing is established.

To summarize, we have demonstrated that the Dynamical Cluster Quantum Monte Carlo, along with the Maximum Entropy Method, may be used to make quantitatively-reliable predictions for the single-particle spectral function of unitary Fermi gases. Our spectral function results capture features seen in recent photoemission experiments (10). We find the onset of a pairing pseudogap with decreasing temperature (above the expected superfluid transition temperature) for coupling values beyond the unitary point in the BEC regime. Since our method treats the short range correlations of the system exactly, it holds particular promise for obtaining quantitatively accurate predictions for radio frequency spectra and photoemission data in the strongly-interacting unitary regime of superfluid Fermi gases. Future work will extend this method to incorporate the effect of a background trapping potential, a necessary step to quantitatively understand experimental results.

###### Acknowledgements.

We acknowledge useful discussions with Dana Browne, Ehsan Khatami, Karlis Mikelsons, Unjong Yu, and Zhaoxin Xu. This research was supported by an allocation of computing time from the Ohio Supercomputer Center, by the NSF OISE-0952300 and the Louisiana Board of Regents, under grant No. LEQSF (2008-11)-RD-A-10.### References

- V. Gurarie and L. Radzihovsky, Ann. of Phys. 322, 2 (2007).
- S. Giorgini, L.P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 80, 1215 (2008).
- W. Ketterle and M. Zwierlein, in Proc. Int. Sch. Phys. “Enrico Fermi” 164, ed. by M. Inguscio, W. Ketterle and C. Salomon, (IOS Press, 2008).
- A. Gezerlis, and J. Carlson, Phys. Rev. C 77, 032801(R) (2008).
- M.G. Alford et al, Rev. Mod. Phys. 80, 1455 (2008).
- Q. Chen et al, Phys. Rep. 412, 1 (2005).
- M. Randeria, in Proc. Int. Sch. Phys. “Enrico Fermi”, 136, ed. by G. Iadonisi, J. R. Schrieffer, and M. L. Chiafalo, (IOS Press, 1998).
- C. Chin et al, Science 305, 1128 (2004).
- C.H. Schunck et al, Nature 454, 739 (2008).
- J.T. Stewart, J.P. Gaebler, and D.S. Jin, Nature 454, 744 (2008).
- P. Massignan, G.M. Bruun, H.T.C. Stoof, Phys. Rev. A 77, 031601 (2008).
- S. Tsuchiya, R. Watanabe, and Y. Ohashi, Phys. Rev. A 80, 033613 (2009).
- R. Haussmann et al, Phys. Rev. A 75, 023610 (2007).
- R. Haussmann, M. Punk, and W. Zwerger, Phys. Rev. A 80, 063612 (2009).
- Q. Chen and K. Levin, Phys. Rev. Lett. 102, 190402 (2009).
- M. Veillette et al, Phys. Rev. A 78, 033614 (2008).
- M. H. Hettler et al, Phys. Rev. B 58, R7475 (1998); M. H. Hettler et al, Phys. Rev. B 61, 12739 (2000).
- Th. Maier et al, Rev. Mod. Phys. 77, 1027 (2005).
- M. Jarrell and J.E. Gubernatis, Phys. Rep. 269, 133 (1996).
- P. Magierski et al, Phys. Rev. Lett. 103, 210403 (2009).
- S. Tan, Ann. Phys. 323, 2952 (2008); Ann. Phys. 323, 2971 (2008).
- A. Bulgac and G.F. Bertsch, Phys. Rev. Lett. 94, 070401 (2005).
- E. Burovski et al, Phys. Rev. Lett. 96, 160402 (2006).
- J. E. Hirsch and R. M. Fye, Phys. Rev. Lett. 56, 2521, (1986).
- M. Jarrell et al, Phys. Rev. B 64, 195130 (2001).
- D. D. Betts, H. Q. Lin and J. S. Flynn, Can. J. Phys. 77, 353 (1999).
- J. Carlson et al, Phys. Rev. Lett. 91, 050401 (2003).
- M.Y. Veillette, D.E. Sheehy, and L. Radzihovsky, Phys. Rev. A 75, 043614 (2007).
- C.A. Regal, et al Phys. Rev. Lett. 95, 250404 (2005).