Confined ultracold atoms in optical lattices provide a versatile platform for simulating lattice models of strongly correlated quantum systems, where pairing phenomena and superfluid phases can be explored under controlled conditions. While the crossover between the Bardeen-Cooper-Schrieffer (BCS) phase and the Bose-Einstein condensation (BEC) is well understood in homogeneous systems, spatial confinement breaks translational symmetry and reshapes correlation patterns, making the BCS-BEC identification in trapped geometries challenging and allowing unconventional phases to emerge with no direct analog in homogeneous systems. Here we present a characterization of the BCS-BEC crossover in harmonically confined one-dimensional Fermi-Hubbard chains. Our analysis combines Density Matrix Renormalization Group (DMRG) simulations and entanglement-based diagnostics with effective models describing the formation of tightly bound fermion pairs. This combined approach enables a detailed understanding of how the interplay between interactions and confinement reshapes the crossover, leading to insulating regions coexisting with persistent superfluid correlations. Within this framework, we further introduce conditioned correlation functions whose power-law decay allows a clear distinction between BCS-like and BEC-like regimes. The consistency between the effective descriptions and the numerical DMRG results yields a unified picture of the crossover in harmonically confined geometries.
The BCS-BEC (Bose-Einstein condensation) crossover in a lattice is a powerful paradigm that describes how a superconductor deviates from the Bardeen-Cooper-Schrieffer physics as the attractive interaction increases. Optical lattices loaded with binary mixtures of cold atoms allow one to access this phenomenon experimentally in a clean and controlled way. We show that, however, the possibility to study this phenomenon in actual cold-atoms experiments is limited by the effect of the trapping potential. Real-space dynamical mean-field theory calculations show indeed that interactions and the confining potential conspire to pack the fermions in the center of the trap, which approaches a band insulator when the attraction becomes sizeable. Interestingly, the energy gap is spatially more homogeneous than the superfluid condensate order parameter. We show how this physics reflects in several observables, and we propose an alternative strategy to disentangle the effect of the harmonic potential and measure the intrinsic properties resulting from the interaction strength.
This work studies zero-temperature Fermi gases with short-ranged and attractive interactions with focus on the dimensionality aspects by confining the system from three to two dimensions. Interacting Fermi gases are superfluids characterized as being strongly correlated if the scattering length is much larger than the other scales present. In the many-body context, this category of systems requires a careful treatment since mean-field approaches are not expected to hold under the conditions which define the strong pairing. For this reason, we employ quantum Monte Carlo methods to study the dimensionality effects in two-component Fermi gases in the Bardeen-Cooper-Schrieffer to Bose-Einstein condensate (BCS-BEC) crossover. To be able to investigate such aspects, we consider a three-dimensional gas with periodic boundary conditions in the xy plane and a confining potential in the z direction capable of reducing the dimensionality of the system, which defines a dimensional BCS-BEC crossover. We calculate the equation of state as the energy per particle as a function of the dimensionality and as a function of the interaction parameter, which are characterized by the oscillator length in the z direction and the scattering length respectively. The greater confinement in the system lowers the interparticle distances, which increases the effects of the attractive interactions between antiparallel-spin fermions.We also calculate the density profile in the z direction over the crossover for different confinements and different interaction parameters to quantify their influence on the spatial distribution of the particles. This study connects the perspectives between theoretical and experimental works in low-dimensional ultracold Fermi gases by utilizing an approach that is similar to what it is realized in modern experiments.
The continuous crossover between a Bardeen-Cooper-Schrieffer (BCS)-type superfluid of fermion pairs and a Bose-Einstein condensate (BEC) of tightly bound bosonic molecules can be attributed to the spontaneous breaking of global
U(1) gauge symmetry which underlies both quantum condensation phenomena. Recently much attention has been paid to this problem, since Feshbach resonances allow for an experimental implementation of crossover physics in cold fermion gases. The strong interactions close to resonance call for an analysis beyond Mean Field Theory. We develop a systematic functional integral approach for the description of this phenomenon. Starting from a Yukawa-type atom-molecule model, a symmetry analysis allows to both construct the equation of state and to classify the thermodynamic phases in a unified way. The onset of superfluidity is signalled by the emergence of a massless Goldstone mode associated with the broken continuous U(1) symmetry. Beyond Mean Field, we include fluctuations of the molecule field self-consistently via the solution of suitable Schwinger-Dyson equations. The phase diagram is computed, and a variety of universal features are established. A new form of crossover from an exactly solvable narrow resonance limit to broad resonances or pointlike interactions is found. At low temperature our results agree well with quantum Monte Carlo simulations and recent experiments. Our approach is further developed in the frame of functional renormalization group equations. While the effective bosonic theory in the BEC regime shows the characteristics of a Bogoliubov theory for small temperatures, the phase transition is of second order.
Recent experiments with ultracold fermionic atoms in optical lattices have provided a tuneable and clean realization of the attractive Hubbard model (AHM). In view of this, several physical properties may be thoroughly studied across the crossover between weak [Bardeen-Cooper-Schrieffer, (BCS)] and strong [Bose-Einstein condensation (BEC)] couplings. Here we report on extensive determinant Quantum Monte Carlo (DQMC) studies of the AHM on a square lattice from which several different quantities have been calculated and should be useful as a road map to experiments. We have obtained a detailed phase diagram for the critical superconducting temperature
Tc in terms of the band filling $\ensuremath{\langle}n\ensuremath{\rangle}$ and interaction strength
U from which we pinpoint a somewhat wide region $|U|/t\ensuremath{\approx}5\ifmmode\pm\else\textpm\fi{}1$ (
t is the hopping amplitude) and $\ensuremath{\langle}n\ensuremath{\rangle}\ensuremath{\approx}0.79\ifmmode\pm\else\textpm\fi{}0.09$ leading to a maximum ${T}_{c}\ensuremath{\approx}0.16t$. Two additional temperature scales, namely, pairing
Tp and degeneracy
Td have been highlighted: The former sets the scale for pair formation (believed to be closely related to the scale for the gap of spin excitations in cuprates), whereas the latter sets the scale for dominant quantum effects. Our DQMC data for the distribution of doubly occupied sites for the momentum distribution function, and for the quasiparticle weight show distinctive features on both sides of the BCS-BEC crossover, being also suggestive of an underlying crossover between Fermi- and non-Fermi liquid behaviors.
We investigate the crossover from the Bardeen-Cooper-Schrieffer (BCS) state of weakly bound Cooper pairs to the Bose-Einstein condensate (BEC) of strongly bound molecular dimers in a gas of ultracold atoms loaded on a two-dimensional optical lattice. By using the mean-field BCS equations of the emerging Hubbard model and the concept of off-diagonal long-range order for fermions we calculate analytically and numerically the pair binding energy, the energy gap, and the condensate fraction of Cooper pairs as a function of interaction strength and filling factor of atoms in the lattice at zero temperature.
Recent experiments with ultracold fermionic atoms in optical lattices have provided a tuneable and clean realization of the attractive Hubbard model (AHM). In view of this, several physical properties may be thoroughly studied across the crossover between weak (Bardeen-Cooper-Schrieffer, BCS) and strong (Bose-Einstein condensation, BEC) couplings. Here we report on extensive determinant Quantum Monte Carlo (DQMC) studies of the AHM on a square lattice, from which several different quantities have been calculated and should be useful as a roadmap to experiments. We have obtained a detailed phase diagram for the critical superconducting temperature,
Tc, in terms of the band filling, $\ave{n}$, and interaction strength,
U, from which we pinpoint a somewhat wide region
∣U∣/t≈5±1 (
t is the hopping amplitude) and $\ave{n} \approx 0.79 \pm 0.09$ leading to a maximum
Tc≈0.16t. Two additional temperature scales, namely pairing,
Tp, and degeneracy,
Td, have been highlighted: the former sets the scale for pair formation (believed to be closely related to the scale for the gap of spin excitations in cuprates), while the latter sets the scale for dominant quantum effects. Our DQMC data for the distribution of doubly occupied sites, for the momentum distribution function, and for the quasiparticle weight show distinctive features on both sides of the BCS-BEC crossover, being also suggestive of an underlying crossover between Fermi- and non-Fermi liquid behaviors.