Multi-component ultracold Fermi gases in optical lattices realize the multi-component Fermi-Hubbard model, in which pairing is renormalized by induced (particle-hole screening) interactions between components. Weak-coupling mean-field theory gives explicit exponential corrections to s- and d-wave gaps and mean-field transition temperatures, but it fails precisely where the physics is richest: across the BCS-BEC crossover, near half filling in two dimensions where the screening interaction becomes singular, in the strongly correlated regime relevant to d-wave pairing, and for the true (Berezinskii-Kosterlitz-Thouless) transition temperature. The open problem is a controlled, quantitative many-body treatment of induced-interaction effects on pairing in multi-component lattice Fermi gases valid in these regimes, including the interplay with antiferromagnetic and Mott-insulating orders, which existing partial results (two-component lattice calculations of the Gorkov-Melik-Barkhudarov correction, continuum crossover treatments, and attractive-Hubbard quantum Monte Carlo) leave unresolved.
Background
In a dilute two-component Fermi gas, Bardeen-Cooper-Schrieffer (BCS) theory predicts a superfluid transition temperature set by an exponentially small pairing gap. Gorkov and Melik-Barkhudarov showed that particle-hole fluctuations screen the attractive interaction between the two pairing components, reducing by a factor ; these 'induced interactions' are the leading beyond-mean-field correction in the weak-coupling limit. In a mixture of fermion species (hyperfine states of Li, or the nuclear-spin states of alkaline-earth atoms such as Yb and Sr that realize an approximate SU() symmetry, under which the spin states interact identically), the extra components contribute additional screening loops that can either suppress or enhance pairing. An optical lattice - a standing-wave light field creating a periodic potential for the atoms - reduces the problem to the Fermi-Hubbard model with on-site couplings between components and hopping . On a lattice the pairing gap acquires an exponential induced-interaction correction , where is the density of states and is an angular average of the Lindhard screening function (the standard density response of a Fermi sea); near half filling (one particle per site) in two dimensions, Umklapp processes (momentum transfers that include a reciprocal-lattice vector) make the Fermi surface square-shaped (nested), and diverges, so the correction is no longer perturbatively small. Repulsive on-site interaction with a weak nearest-neighbor attraction favors d-wave pairing (gap changing sign under rotation, as in cuprate superconductors), for which mean-field expressions exist near half filling; but at strong coupling the same lattice realizes Mott-insulating (interaction-localized) and antiferromagnetic (AF) phases that mean field ignores. Strong attraction, tunable by Feshbach resonances (magnetic-field-tuned scattering resonances), drives the BCS-BEC crossover from weakly bound Cooper pairs to tightly bound molecules. Finally, in two dimensions the true transition is of Berezinskii-Kosterlitz-Thouless (BKT) type - vortex-antivortex unbinding at below the mean-field - which the lattice mean-field calculation does not evaluate. Heiselberg's 2009 analysis (arXiv:0905.2107) computed gaps, , and induced interactions to leading orders for multi-component lattice gases, finding, for example, that components enhance the 2D gap by in the dilute regime while near half filling the screening becomes singular; the present problem is everything that analysis explicitly could not do.
Problem Statement
Obtain a controlled, quantitative many-body treatment of induced (screening) interactions and their effect on superfluid pairing in multi-component Fermi gases in optical lattices - i.e., in the multi-component Fermi-Hubbard model with components - beyond the weak-coupling mean-field regime in which current results are derived. This means a theory, with a controlled approximation whose validity is established in the regime it treats, not a re-derivation at the same weak-coupling mean-field level. The single overarching objective is a theory of induced-interaction renormalization of pairing valid where the existing perturbative mean-field framework fails, and its linked scope comprises: (i) the behavior of the induced interaction and the resulting gap and transition-temperature renormalization from weak coupling through the BCS-BEC crossover (strong attraction, up to the lattice unitary threshold where two-body bound states appear) - as motivating context, the source notes that for the continuum uniform gas Nozieres-Schmitt-Rink-type estimates (fluctuation corrections beyond mean field) reduce pairing by only ~20% whereas experiments indicate a factor ~2-3; (ii) the non-perturbative regime near half filling in two dimensions, where Umklapp-enhanced screening makes the induced interaction comparable to or larger than the on-site coupling and the exponential weak-coupling correction breaks down; (iii) the correction of the near-half-filling d-wave mean-field gap by strong on-site correlations, antiferromagnetic exchange, and Mott physics, including how the induced-interaction enhancement or suppression of d-wave pairing survives correlations; and (iv) the true two-dimensional transition temperature (vortex-unbinding BKT transition) including induced interactions and lattice effects, not only the mean-field . Each item is a regime or observable of the same induced-interaction renormalization that the weak-coupling mean-field baseline explicitly identifies as beyond its validity; they are not independent questions, and controlled progress in any of them advances the overarching objective. Relative terms are determinate by field convention: 'weak coupling' means and induced interactions small compared to on-site couplings, 'strong coupling' the approach to or of order several, and 'near half filling' the filling range around per site where the 2D screening function diverges.
References
- H. Heiselberg, Induced interactions and superfluidity in optical lattices with multi-component Fermi gases, arXiv:0905.2107 (2009), https://arxiv.org/abs/0905.2107 (primary source of the open problem)
- D.-H. Kim, P. Törmä, J.-P. Martikainen, Induced Interactions for Ultracold Fermi Gases in Optical Lattices, Phys. Rev. Lett. 102, 245301 (2009), doi:10.1103/PhysRevLett.102.245301, https://doi.org/10.1103/PhysRevLett.102.245301
- J.-P. Martikainen, J. J. Kinnunen, P. Törmä, C. J. Pethick, Induced Interactions and the Superfluid Transition Temperature in a Three-Component Fermi Gas, Phys. Rev. Lett. 103, 260403 (2009), doi:10.1103/PhysRevLett.103.260403, https://doi.org/10.1103/PhysRevLett.103.260403
- H. Heiselberg, C. J. Pethick, H. Smith, L. Viverit, Influence of Induced Interactions on the Superfluid Transition in Dilute Fermi Gases, Phys. Rev. Lett. 85, 2418 (2000), doi:10.1103/PhysRevLett.85.2418, https://doi.org/10.1103/PhysRevLett.85.2418
- L. Pisani, A. Perali, P. Pieri, G. C. Strinati, Entanglement between pairing and screening in the Gorkov-Melik-Barkhudarov correction to the critical temperature throughout the BCS-BEC crossover, Phys. Rev. B 97, 014528 (2018), doi:10.1103/PhysRevB.97.014528, https://doi.org/10.1103/PhysRevB.97.014528
- R. A. Fontenele, N. C. Costa, R. R. dos Santos, T. Paiva, Two-dimensional attractive Hubbard model and the BCS-BEC crossover, Phys. Rev. B 105, 184502 (2022), doi:10.1103/PhysRevB.105.184502, https://doi.org/10.1103/PhysRevB.105.184502
- T. Hazra, N. Verma, M. Randeria, Bounds on the Superconducting Transition Temperature: Applications to Twisted Bilayer Graphene and Cold Atoms, Phys. Rev. X 9, 031049 (2019), doi:10.1103/PhysRevX.9.031049, https://doi.org/10.1103/PhysRevX.9.031049
- M. A. Cazalilla, A. M. Rey, Ultracold Fermi gases with emergent SU(N) symmetry, Rep. Prog. Phys. 77, 124401 (2014), doi:10.1088/0034-4885/77/12/124401, https://doi.org/10.1088/0034-4885/77/12/124401
- C. Feng, E. Ibarra-García-Padilla, K. R. A. Hazzard, R. T. Scalettar, S. Zhang, E. Vitali, Metal-insulator transition and quantum magnetism in the SU(3) Fermi-Hubbard model, Phys. Rev. Research 5, 043267 (2023), doi:10.1103/PhysRevResearch.5.043267, https://doi.org/10.1103/PhysRevResearch.5.043267