We propose a model for addressing the superfluidity of two different Fermi species confined in a bilayer geometry of square optical lattices. The fermions are assumed to be molecules with interlayer s-wave interactions, whose dipole moments are oriented perpendicularly to the layers. Using functional integral techniques we investigate the BCS-like state induced in the bilayer at finite temperatures. In particular, we determine the critical temperature as a function of the coupling strength between molecules in different layers and of the interlayer spacing. By means of Ginzburg-Landau theory we calculate the superfluid density. We also study the dimerized BEC phase through the Berezinskii-Kosterlitz-Thouless transition, where the effective mass leads to identify the crossover from BCS to BEC regimes. The possibility of tuning the effective mass as a direct consequence of the lattice confinement, allows us to suggest a range of values of the interlayer spacing, which would enable observing this superfluidity within current experimental conditions.
We perform a variational quantum Monte Carlo simulation of the transition from a Bardeen-Cooper-Schrieffer superfluid (BCS) to a Bose-Einstein condensate (BEC) at zero temperature. The model Hamiltonian involves an attractive short range two body interaction and the atoms number $2N =330$ is chosen so that, in the non-interacting limit, the ground state function corresponds to a closed shell configuration. The system is then characterized by the s-wave scattering length $a$ of the two-particle collisions in the gas, which is varied from negative to positive values, and the Fermi wave number $k_F$. Based on an extensive analysis of the s-wave two-body problem, one parameter variational many-body wave functions are proposed to describe the ground state of the interacting Fermi gas from BCS to BEC states. We exploit properties of antisymmetrized many-body functions to develop efficient techniques that permit variational calculations for a large number of particles. It is shown that a virial relation between the energy per particle and the trapping energy is approximately valid for $-0.1<1/k_Fa<3.4$. The influence of the harmonic trap and the interaction potential as exhibited in two-body correlation functions is also analyzed.
