Motivated by a recent experiment Revelle et al., [Phys. Rev. Lett. 117, 235301 (2016)] that characterized the one- to three-dimensional crossover in a spin-imbalanced ultracold gas of $^{6}\mathrm{Li}$ atoms trapped in a two-dimensional array of tunnel-coupled tubes, we calculate the phase diagram for this system by using Hartree-Fock Bogoliubov-de Gennes mean-field theory and compare the results with experimental data. Mean-field theory predicts fully-spin-polarized normal, partially-spin-polarized normal, spin-polarized superfluid, and spin-balanced superfluid phases in a homogeneous system. We use the local density approximation to obtain density profiles of the gas in a harmonic trap. We compare these calculations with experimental measurements in Revelle et al. as well as previously unpublished data. Our calculations qualitatively agree with experimentally measured densities and coordinates of the phase boundaries in the trap and quantitatively agree with experimental measurements at moderate-to-large polarizations. Our calculations also reproduce the experimentally observed universal scaling of the phase boundaries for different scattering lengths at a fixed value of scaled intertube tunneling. However, our calculations have quantitative differences with experimental measurements at low polarization and fail to capture important features of the one- to three-dimensional crossover observed in experiments. These suggest the important role of physics beyond mean-field theory in the experiments. We expect that our numerical results will aid future experiments in narrowing the search for the Fulde-Ferrell-Larkin-Ovchinnikov phase.
Motivated by a recent experiment Revelle et al., [Phys. Rev. Lett. 117, 235301 (2016)] that characterized the one- to three-dimensional crossover in a spin-imbalanced ultracold gas of 6Li atoms trapped in a two-dimensional array of tunnel-coupled tubes, we calculate the phase diagram for this system by using Hartree-Fock Bogoliubov-de Gennes mean-field theory and compare the results with experimental data. Mean-field theory predicts fully-spin-polarized normal, partially-spin-polarized normal, spin-polarized superfluid, and spin-balanced superfluid phases in a homogeneous system. We use the local density approximation to obtain density profiles of the gas in a harmonic trap. We compare these calculations with experimental measurements in Revelle et al. as well as previously unpublished data. Our calculations qualitatively agree with experimentally measured densities and coordinates of the phase boundaries in the trap and quantitatively agree with experimental measurements at moderate-to-large polarizations. Our calculations also reproduce the experimentally observed universal scaling of the phase boundaries for different scattering lengths at a fixed value of scaled intertube tunneling. However, our calculations have quantitative differences with experimental measurements at low polarization and fail to capture important features of the one- to three-dimensional crossover observed in experiments. These suggest the important role of physics beyond mean-field theory in the experiments. We expect that our numerical results will aid future experiments in narrowing the search for the Fulde-Ferrell-Larkin-Ovchinnikov phase.
We present a systematic study of quantum phases in a one-dimensional spin-polarized Fermi gas. Three comparative theoretical methods are used to explore the phase diagram at zero temperature: the mean-field theory with either an order parameter in a single-plane-wave form or a self-consistently determined order parameter using the Bogoliubov--de Gennes equations, as well as the exact Bethe ansatz method. We find that a spatially inhomogeneous Fulde-Ferrell-Larkin-Ovchinnikov phase, which lies between the fully paired Bardeen-Cooper-Schrieffer (BCS) state and the fully polarized normal state, dominates most of the phase diagram of a uniform gas. The phase transition from the BCS state to the Fulde-Ferrell-Larkin-Ovchinnikov phase is of second order, and therefore there are no phase separation states in one-dimensional homogeneous polarized gases. This is in sharp contrast to the three-dimensional situation, where a phase separation regime is predicted to occupy a very large space in the phase diagram. We conjecture that the prediction of the dominance of the phase separation phases in three dimension could be an artifact of the non-self-consistent mean-field approximation, which is heavily used in the study of three-dimensional polarized Fermi gases. We consider also the effect of a harmonic trapping potential on the phase diagram, and find that in this case the trap generally leads to phase separation, in accord with the experimental observations for a trapped gas in three dimensions. We finally investigate the local fermionic density of states of the Fulde-Ferrell-Larkin-Ovchinnikov ansatz. A two-energy-gap structure appears, which could be used as an experimental probe of the Fulde-Ferrell-Larkin-Ovchinnikov states.
We present a theory for a lattice array of weakly coupled one-dimensional ultracold attractive Fermi gases (1D tubes) with spin imbalance, where strong intratube quantum fluctuations invalidate mean-field theory. We first construct an effective field theory, which treats spin-charge mixing exactly, based on the Bethe ansatz solution of the 1D single tube problem. We show that the 1D Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) state is a two-component Luttinger liquid, and its elementary excitations are fractional states carrying both charge and spin. We analyze the instability of the 1D FFLO state against intertube tunneling by renormalization group analysis, and find that it flows into either a polarized Fermi liquid or a FFLO superfluid, depending on the magnitude of interaction strength and spin imbalance. We obtain the phase diagram of the quasi-1D system and further determine the scaling of the superfluid transition temperature with intertube coupling.
We investigate the finite-temperature properties of an ultracold atomic Fermi gas with spin population imbalance in a highly elongated harmonic trap. Previous studies at zero temperature showed that the gas stays in an exotic spatially inhomogeneous Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superfluid state at the trap center; while moving to the edge, the system changes into either a nonpolarized Bardeen-Cooper-Schrieffer superfluid (P<Pc) or a fully polarized normal gas (P>Pc), depending on the smallness of the spin polarization P, relative to a critical value Pc. In this work, we show how these two phase-separation phases evolve with increasing temperature, and thereby construct a finite-temperature phase diagram. For typical interactions, we find that the exotic FFLO phase survives below one-tenth of Fermi degeneracy temperature, which seems to be accessible in the current experiment. The density profile, equation of state, and specific heat of the polarized system have been calculated and discussed in detail. Our results are useful for the ongoing experiment at Rice University on the search for FFLO states in quasi-one-dimensional polarized Fermi gases.
We present an extensive study of two-dimensional Larkin-Ovchinnikov (LO) superfluidity in a spin-imbalanced two-component atomic Fermi gas. In the context of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase, we explore a wide and generic class of pairing gap functions with explicit spatial dependency. The mean-field theory of such phases is applied through the Bogoliubov-de Gennes equations in which the pairing gap can be determined self-consistently. To systematically explore the configuration space we consider both the canonical and grand canonical ensembles where we control the polarization or chemical potentials of the system, respectively. The mean-field calculations enable us to understand the nature of the phase transitions in the fully paired Bardeen-Cooper-Schrieffer (BCS) state, exotic LO phase, and partially polarized free Fermi gas (NPP). The order of the phase transitions has been examined and, in particular, we find a weak first-order phase transition between the exotic inhomogeneous LO phase and the BCS phase. In comparison to the three-dimensional case, where the phase diagram is dominated by a generic separation phase, we predict a broader parameter space for the spatially inhomogeneous LO phase. By computing the superfluid density of the LO phase at different polarization, we show how the superfluidity of the system is suppressed with increasing spin polarization.
The Fulde-Ferrell Larkin-Ovchinnikov (FFLO) type superconductor is an exotic form of matter known as a supersolid because it simultaneously supports superfluid and magnetic order. In both condensed matter and ultracold atomic gases, this phase has yet to be conclusively observed despite long interest in both theory and experiment. In one-dimension (1D), unlike expectations for higher dimensions, the FFLO phase is found in a large region of the phase diagram [1]. The effect of quantum and thermal fluctuations, however, are expected to be reduced in higher dimensions. These considerations motivated the proposal to search for FFLO near the 1D-3D dimensional crossover [2], which we have identified and characterized [3]. We confine a spin-imbalanced Fermi gas of lithium 6 to an array of 1D tubes using a 2D optical lattice. We bring the system into the dimensional crossover by increasing the inter-tube tunneling rate and using a Feshbach resonance to tune interactions. The periodicity of these domain walls, which depends on the magnitude of the polarization, is a definitive signature of the FFLO phase. In this work, we develop new experimental and analytical methods to extract periodic signals that, if they indeed are domain walls, lay the foundation for full confirmation of the FFLO phase. [1] Y. A. Liao, A. S. C. Rittner, T. Paprotta, W. Li, G. B. Partridge, R. G. Hulet, S. K. Baur, and E. J. Mueller, “Spin-imbalance in a one-dimensional Fermi gas,” Nature 467, 567 (2010). [2] M. M. Parish, S. K. Baur, E. J. Mueller, and D. A. Huse, “Quasi-One- Dimensional Polarized Fermi Superfluids,” Phys. Rev. Lett. 99, 250403 (2007). [3] M. C. Revelle, J. A. Fry, B. A. Olsen, and R. G. Hulet, “1D to 3D Crossover of a Spin-Imbalanced Fermi Gas,” Physical Review Letters 117 (2016).
Ultracold atomic gases have emerged as an ideal laboratory system to emulate many-body physics in an unprecedentedly controllable manner. Numerous many-body quantum states and phases have been experimentally explored and characterized using the ultracold atomic gases, offering new insights into many exciting physics ranging from condensed matters to cosmology. In this thesis, we will present a systematic numerical study of a novel experimental system, population imbalanced two-component ultracold Fermi gases. We explore the phase diagram of this system in both 3D and 1D especially focusing on the exotic Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase, which is characterized by a spatially oscillating order parameter. In 3D, we solve for the stationary states of trapped imbalanced Fermi gases in a wide range of parameter space with a home-made parallel eigen-solver for Bogoliubov-de Gennes (BdG) equations. Our results show that there exists a metastable state with a FFLO type oscillating order parameter. In 1D, we simulate the dynamical expansion of the population imbalanced Fermi gases from the trap. A numerically quasi-exact scheme, time-evolving block decimation (TEBD), is introduced for the comparative studies with the solution of the time-dependent BdG equation. Our results predict that the existence of FFLO states will leave conspicuous signatures in the density profiles during the expansion. For further understanding of the interplay between the population imbalance and two-body pairing interaction between two spin components, we also study the spin transport properties through trapped ultracold Fermi gases. The preliminary results will be discussed.