Strongly correlated fermionic many-body systems are ubiquitous in nature. Their theoretical description poses challenging problems which are further complicated when imbalances in, e.g., the particle numbers of the involved species or their masses are introduced. In this thesis, a number of different approaches is developed and applied in order to obtain predictions for physical observables of such systems that mutually support and confirm each other. In a first step, analytically well-founded mean-field analyses are carried through. One- and three-dimensional ultracold Fermi gases with spin and mass imbalance as well as Gross-Neveu and NJL-type relativistic models at finite baryon chemical potential are investigated with respect to their analytic properties in general and the occurrence of spontaneous breaking of translational invariance in particular. Based on these studies, further methods are devised or adapted allowing for investigations also beyond the mean-field approximation. Lattice Monte Carlo simulations with imaginary imbalance parameters are employed to surmount the infamous sign problem and compute the equation of state of the respective unitary Fermi gases. Moreover, in-medium two-body analyses are used to confirm and explain the characteristics of inhomogeneously ordered phases. Finally, functional RG methods are applied to the unitary Fermi gas with spin and mass imbalance. Besides quantitatively competitive predictions for critical temperatures for the superfluid state, strong hints on the stability of inhomogeneous phases with respect to order parameter fluctuations in the regime of large mass imbalance are obtained. Combining the findings from these different theoretical studies suggests the possibility to find such phases in experiments presently in preparation.
Strongly correlated fermionic many-body systems are ubiquitous in nature. Their theoretical description poses challenging problems which are further complicated when imbalances in, e.g., the particle numbers of the involved species or their masses are introduced. In this thesis, a number of different approaches is developed and applied in order to obtain predictions for physical observables of such systems that mutually support and confirm each other. In a first step, analytically well-founded mean-field analyses are carried through. One- and three-dimensional ultracold Fermi gases with spin and mass imbalance as well as Gross-Neveu and NJL-type relativistic models at finite baryon chemical potential are investigated with respect to their analytic properties in general and the occurrence of spontaneous breaking of translational invariance in particular. Based on these studies, further methods are devised or adapted allowing for investigations also beyond the mean-field approximation. Lattice Monte Carlo simulations with imaginary imbalance parameters are employed to surmount the infamous sign problem and compute the equation of state of the respective unitary Fermi gases. Moreover, in-medium two-body analyses are used to confirm and explain the characteristics of inhomogeneously ordered phases. Finally, functional RG methods are applied to the unitary Fermi gas with spin and mass imbalance. Besides quantitatively competitive predictions for critical temperatures for the superfluid state, strong hints on the stability of inhomogeneous phases with respect to order parameter fluctuations in the regime of large mass imbalance are obtained. Combining the findings from these different theoretical studies suggests the possibility to find such phases in experiments presently in preparation.
Strongly coupled quantum matter displays a rich phenomenology including phase transitions and often unexpected collective behavior. Remarkable advances in experiments with ultracold Fermi gases allow us to gain deep insight into these intriguing systems. Their theoretical description, however, is often challenging as exact analytic solutions are available only in a few special cases, and approximate techniques such as mean-field or perturbation theory are of limited use. Numerical treatment with Monte Carlo (MC) methods has led to profound success in this regard. Unfortunately, for many systems - and especially for asymmetric quantum gases - the infamous sign problem slows progress due to an exponentially scaling of the computational effort with inceasing system size. In this thesis, we set out to explore the rich physics of two-component Fermi gases in the presence of finite spin polarization and/or mass imbalance. To surmount an arising sign problem, we learn from methodological advances made in the field of quantum chromodynamics and further develop these lattice approaches in the context of nonrelativistic Fermi gases. An extensive overview of the numerical methods is presented, including several toy problems to detail the capabilities and shortcomings of the developed approaches. With these tools in hand, we perform extensive benchmarks of the hybrid Monte Carlo method with imaginary asymmetries (iHMC) and the complex Langevin (CL) method, which is based on a complex version of stochastic quantization. Both approaches are shown to yield excellent results for the ground-state energy equation of state of mass-imbalanced Fermi gases in one spatial dimension. Due to its great versatility, the CL method is subsequently employed to study pairing in one-dimensional Fermi gases, for which suitable two-body correlations are computed, revealing unexpected pairing patterns for spin- and mass-imbalanced systems. Another major system of interest in this thesis is the paradigmatic unitary Fermi gas which is investigated at finite temperature and spin polarization. A precise determination of the density equation of state in the normal phase enables us to explore a broad range of thermodynamic properties. We infer valuable information on the finite-temperature phase diagram, such as a flat phase boundary of the normal-to-superfluid transition near the balanced limit and indications for the absence of an extensive pseudogap phase above this transition. The presented results provide experimentally testable ab initio predictions for a range of previously inaccessible thermodynamic quantities.
We study the phase diagram of mass- and spin-imbalanced unitary Fermi gases, in search for the emergence of spatially inhomogeneous phases. To account for fluctuation effects beyond the mean-field approximation, we employ renormalization group techniques. We thus obtain estimates for critical values of the temperature, mass and spin imbalance, above which the system is in the normal phase. In the unpolarized, equal-mass limit, our result for the critical temperature is in accordance with state-of-the-art Monte Carlo calculations. In addition, we estimate the location of regions in the phase diagram where inhomogeneous phases are likely to exist. We show that an intriguing relation exists between the general structure of the many-body phase diagram and the binding energies of the underlying two-body bound-state problem, which further supports our findings. Our results suggest that inhomogeneous condensates form for mass imbalances ${m}_{\ensuremath{\downarrow}}/{m}_{\ensuremath{\uparrow}}\ensuremath{\gtrsim}3$. The extent of the inhomogeneous phase in parameter space increases with increasing mass imbalance.
Fermi gases in strongly coupled regimes are inherently challenging for many-body methods. Although progress has been made analytically, quantitative results require ab initio numerical approaches, such as Monte-Carlo (MC) calculations. However, mass-imbalanced and spin-imbalanced gases are not accessible to MC calculations due to the infamous sign problem. For finite spin imbalance, the problem can be circumvented using imaginary polarizations and analytic continuation, and large parts of the phase diagram then become accessible. We propose to apply this strategy to the mass-imbalanced case, which opens up the possibility to study the associated phase diagram with MC calculations. We perform a first mean-field analysis which suggests that zero-temperature studies, as well as detecting a potential (tri)critical point, are feasible.
Two-component Fermi gases model the behavior of many systems in different fields of physics, and one of their interesting features is that they condense into superfluids at low temperatures. Ultracold atoms experiments represent one realization of such Fermi gases, and their great flexibility sparked active research into their phase structure. In the literature, there are studies of the phase structure using functional methods, mean-field approximations, and other approaches. In the present work, we aim to perform ab-initio Monte-Carlo simulations of the system to probe their phase structure for inhomogeneous phases in the presence of spin imbalance. Such simulations generally require a bosonization of the theoretical description of the system that rewrites the theory in terms of an auxiliary bosonic field. For this auxiliary field, many possible choices achieve this, and previous studies often use a field that corresponds to a density of fermions. In the present work, we develop a novel approach to this problem by bosonizing the system in terms of the so-called pairing field which corresponds to the superfluid order parameter. Even in the absence of mass- or spin-imbalance, this approach exhibits a sign problem that presents a problem to many Monte-Carlo methods. To circumvent this sign problem, we base our simulation on the Complex Langevin method. The central question of this thesis is if a simulation of a pairing-field-based formalism using the Complex-Langevin method is suitable to study the phase structure of two-component Fermi gases in the presence of spin imbalance. To this end, we develop a lattice theory based on the pairing field by discretizing the continuous Hamiltonian and performing a rigorous derivation of the path integral from there. We pay special attention to the derivation of lattice derivative operators and rescale the theory to be dimensionless. Beyond that, we develop an efficient notation for lattice theories that makes their handling arguably easier than that of continuum theories. With the obtained theory, we derive the Langevin equation for the system together with the expressions we need to sample observables and develop a numerical simulation on that basis. We use the simulation to study 0+1-dimensional systems as a proof-of-concept. Specifically, we calculate density equations of state and two-point functions for a range of dimensionless couplings and compare them to exact solutions. Additionally, we compare our obtained results for the two-point functions to general discussions of the analytic properties of correlation functions to gain additional insight into the qualitative and numerical behavior of our simulation. We found that our results for both density equations of state and two-point functions are in excellent agreement with the exact solutions. Going forward, we plan to study systems in d>0 spatial dimensions, beginning with simulations of 1+1-dimensional systems. In d=3 dimensions the pairing field's correspondence to the superfluid order parameter may allow us to efficiently study the spontaneous breaking of the U(1) symmetry of the system. On that basis, we can probe the phase diagram of the system for inhomogeneous phases at finite spin imbalance in future studies.
We study the phase diagram of strongly interacting fermions in one dimension at finite temperature, with mass and spin imbalance. By including the possibility of the existence of a spatially inhomogeneous ground state, we find regions where spatially varying superfluid phases are favored over homogeneous phases. We obtain estimates for critical values of the temperature, mass, and spin imbalance, above which these phases disappear. Finally, we show that an intriguing relation exists between the general structure of the phase diagram and the binding energies of the underlying two-body bound-state problem. Although our present study of the many-body aspects is based on a mean-field approximation, it may still yield phenomenologically important insights into the dynamics of strongly interacting fermions in one dimension. Most importantly, however, we expect that our present mean-field study can be considered worthwhile from a field-theoretical point of view as it may help to better understand the mechanisms underlying the emergence of inhomogeneous phases in higher dimensions as well as the general structure of the associated phase diagrams.
Fermionic matter is ubiquitous in nature, from the electrons in metals and semiconductors or the neutrons in the inner crust of neutron stars, to gases of fermionic atoms, like 40K or 6Li that can be created and studied under laboratory conditions. It is especially interesting to study these systems at very low temperatures, where we enter the world of quantum mechanical phenomena. Due to the Fermi-Dirac statistics, a dilute system of spin-polarised fermions exhibits no interactions and can be viewed as an ideal Fermi gas. However, interactions play a crucial role for fermions of several spin species. This thesis addresses several questions concerning interacting Fermi gases, in particular the transition between the normal and the superfluid phase and dynamical properties at higher temperatures. First we will look at the unitary Fermi gas: a two-component system of fermions interacting with divergent scattering length. This system is particularly interesting as it exhibits universal behaviour. Due to the strong interactions perturbation theory is inapplicable and no exact theoretical description is available. I will describe the Determinant Diagrammatic Monte Carlo algorithm with which the unitary Fermi gas can be studied from first principles. This algorithm fails in the presence of a spin imbalance (unequal number of particles in the two components) due to a sign problem. I will show how to apply reweighting techniques to generalise the algorithm to the imbalanced case, and present results for the critical temperature and other thermodynamic observables at the critical point, namely the chemical potential, the energy per particle and the contact density. These are the first numerical results for the imbalanced unitary Fermi gas at finite temperature. I will also show how temperatures beyond the critical point can be accessed and present results for the equation of state and the temperature dependence of the contact density. At sufficiently high temperatures a semiclassical description captures all relevant physical features of the system. The dynamics of an interacting Fermi gas can then be studied via a numerical simulation of the Boltzmann equation. I will describe such a numerical setup and apply it to study the collision of two spin-polarised fermionic clouds. When the two components are separated in an elongated harmonic trap and then released, they collide and for sufficiently strong interactions can bounce off each other several times. I will discuss the different types of the qualitative behaviour, show how they can be interpreted in terms of the equilibrium properties of the system, and explain how they relate to the coupling between different excitation modes. I will also demonstrate how transport coefficients, for instance the spin drag, can be extracted from the numerical data.