By using the diffusion Monte Carlo method we calculate the one-and two-body density matrix of an interacting Fermi gas at T = 0 in the BCS-BEC crossover. Results for the momentum distribution of the atoms, as obtained from the Fourier transform of the one-body density matrix, are reported as a function of the interaction strength. Off-diagonal long-range order in the system is investigated through the asymptotic behavior of the two-body density matrix. The condensate fraction of fermionic pairs is calculated in the unitary limit and on both sides of the BCS-BEC crossover.
PACS numbers:The physics of the crossover from Bardeen-CooperSchrieffer (BCS) superfluids to molecular Bose-Einstein condensates (BEC) in ultracold Fermi gases near a Feshbach resonance is a very exciting field that has recently attracted a lot of interest, both from the experimental [1,2] and the theoretical side [3]. An important experimental achievement is the observation of a condensate of pairs of fermionic atoms on the side of the Feshbach resonance where no stable molecules would exist in vacuum [4,5]. Although the interpretation of the experiment is not straightforward, as it involves an out-of-equilibrium projection technique of fermionic pairs onto bound molecules [6], it is believed that these results strongly support the existence of a superfluid order parameter in the strongly correlated regime on the BCS side of the resonance [5].The occurrence of off-diagonal long-range order (ODLRO) in interacting systems of bosons and fermions was investigated by C.N. Yang in terms of the asymptotic behavior of the one-and two-body density matrix [7]. In the case of a two-component Fermi gas with N ↑ spin-up and N ↓ spin-down particles, the one-body density matrix (OBDM) for spin-up particles, defined asdoes not possess any eigenvalue of order N ↑ . This behavior implies for homogeneous systems the asymptotic condition ρ 1 (rIn the above expression ψ † ↑ (r) (ψ ↑ (r)) denote the creation (annihilation) operator of spin-up particles. The same result holds for spin-down particles. ODLRO may occur instead in the two-body density matrix (TBDM), that is defined asFor an unpolarized gas with N ↑ = N ↓ = N/2, if ρ 2 has an eigenvalue of the order of the total number of particles N , the TBDM can be written as a spectral decomposition separating the largest eigenvalue,2 containing only eigenvalues of order one. The parameter α ≤ 1 in Eq. (3) is interpreted as the condensate fraction of pairs, in a similar way as the condensate fraction of single atoms is derived from the OBDM.The spectral decomposition (3) yields for homogeneous systems the following asymptotic behavior of the TBDMThe complex function ϕ is proportional to the order parameter ψ ↑ (r 1 )ψ ↓ (r 2 ) = αN/2ϕ(|r 1 − r 2 |), whose appearance distinguishes the superfluid state of the Fermi gas. Equation (4) should be contrasted with the behavior of Bose systems with ODLRO, where ρ 1 has an eigenvalue of order N [8], and consequently the largest eigenvalue of ρ 2 is of the order of N 2 . In thi...