It is generally assumed that a condensate of paired fermions at equilibrium is characterized by a macroscopic wavefunction with a well-defined, immutable phase. In reality, all systems have a finite size and are prepared at non-zero temperature; the
condensate has then a finite coherence time, even when the system is isolated in its evolution and the particle number $N$ is fixed. The loss of phase memory is due to interactions of the condensate with the excited modes that constitute a dephasing environment. This fundamental effect, crucial for applications using the condensate of pairs macroscopic coherence, was scarcely studied. We link the coherence time to the condensate phase dynamics, and we show with a microscopic theory that the time derivative of the condensate phase operator $hat{theta}_0$ is proportional to a chemical potential operator that we construct including both the pair-breaking and pair-motion excitation branches. In a single realization of energy $E$, $hat{theta}_0$ evolves at long times as $-2mu_{rm mc}(E)t/hbar$ where $mu_{rm mc}(E)$ is the microcanonical chemical potential; energy fluctuations from one realization to the other then lead to a ballistic spreading of the phase and to a Gaussian decay of the temporal coherence function with a characteristic time $propto N^{1/2}$. In the absence of energy fluctuations, the coherence time scales as $N$ due to the diffusive motion of $hat{theta}_0$. We propose a method to measure the coherence time with ultracold atoms, which we predict to be tens of milliseconds for the canonical ensemble unitary Fermi gas.
Due to atomic interactions and dispersion in the total atom number, the order parameter of a pair-condensed Fermi gas experiences a collapse in a time that we derive microscopically. As in the bosonic case, this blurring time depends on the derivativ
e of the gas chemical potential with respect to the atom number and on the variance of that atom number. The result is obtained first using linearized time-dependent Bogoliubov-de Gennes equations, then in the Random Phase Approximation, and then it is generalized to beyond mean field. In this framework, we construct and compare two phase operators for the paired fermionic field: The first one, issued from our study of the dynamics, is the infinitesimal generator of adiabatic translations in the total number of pairs. The second one is the phase operator of the amplitude of the field of pairs on the condensate mode. We explain that these two operators differ due to the dependence of the condensate wave function on the atom number.
We propose and analyze a scheme to entangle the collective spin states of two spatially separated bimodal Bose-Einstein condensates. Using a four-mode approximation for the atomic field, we show that elastic collisions in a state-dependent potential
simultaneously create spin-squeezing in each condensate and entangle the collective spins of the two condensates. We investigate mostly analytically the non-local quantum correlations that arise in this system at short times and show that Einstein-Podolsky-Rosen (EPR) entanglement is generated between the condensates. At long times we point out macroscopic entangled states and explain their structure. The scheme can be implemented with condensates in state-dependent microwave potentials on an atom chip.
