A generic scenario of a crossover between weak-coupling and strong-coupling limits of a many-body quantum system is studied. Powerful numerical methods of investigating the strongly correlated regime are developed and applied to answering a number of important questions of the strongly correlated fermions physics. The internal structure of an exciton in an insulating medium is studied throughout the crossover between the (analytically solvable) strong- and weak-coupling limits. In the intermediate regime, accurate results for the exciton ground state are obtained and used to establish the limits of applicability of the corresponding limiting cases for the first time. A novel combination of the diagrammatic determinant Monte Carlo and the worm algorithm is developed and used for investigating the unitarity regime of the BEC-BCS crossover via simulating the dilute Fermi-Hubbard model. Superfluid phase transition and thermodynamic properties of the normal state are characterized for both lattice fermions and continuum unitary Fermi gas. The relation to current experiments with trapped ultracold quantum gases is discussed.
A generic scenario of a crossover between weak-coupling and strong-coupling limits of a many-body quantum system is studied. Powerful numerical methods of investigating the strongly correlated regime are developed and applied to answering a number of important questions of the strongly correlated fermions physics. The internal structure of an exciton in an insulating medium is studied throughout the crossover between the (analytically solvable) strong- and weak-coupling limits. In the intermediate regime, accurate results for the exciton ground state are obtained and used to establish the limits of applicability of the corresponding limiting cases for the first time. A novel combination of the diagrammatic determinant Monte Carlo and the worm algorithm is developed and used for investigating the unitarity regime of the BEC-BCS crossover via simulating the dilute Fermi-Hubbard model. Superfluid phase transition and thermodynamic properties of the normal state are characterized for both lattice fermions and continuum unitary Fermi gas. The relation to current experiments with trapped ultracold quantum gases is discussed.
Ultracold atomic gases may be the ultimate quantum simulator.These isolated systems have the lowest temperatures in the observable universe, and their properties and interactions can be precisely and accurately tuned across a full spectrum of behaviors, from few-body physics to highly-correlated many-body effects.The ability to impose potentials on and tune interactions within ultracold gases to mimic complex systems mean they could become a theorist's playground.One of their great strengths, however, is also one of the largest obstacles to this dream: isolation.This thesis touches on both of these themes.First, methods to characterize phases and quantum critical points, and to construct finite temperature phase diagrams using experimentally accessible observables in the Bose Hubbard model are discussed.Then, the transition from a weakly to a strongly interacting Bose-Fermi mixture in the continuum is analyzed using zero temperature numerical techniques.Real materials can be emulated by ultracold atomic gases loaded into optical lattice potentials.We discuss the characteristics of a single boson species trapped in an optical lattice (described by the Bose Hubbard model) and the hallmarks of the quantum critical region that separates the superfluid and the Mott insulator ground states.We propose a method to map the quantum critical region using the single, experimentally accessible, local quantity R, the ratio of compressibility to local number fluctuations.The procedure to map a phase diagram with R is easily generalized to inhomogeneous systems and generic many-body Hamiltonians.We illustrate it here using quantum Monte Carlo simulations of the 2D Bose Hubbard model.Secondly, we investigate the transition from a degenerate Fermi gas weakly coupled List of Tables Table Page 2.1 Techniques to Probe Atoms in Optical Lattices . . . . . .
Ultracold atomic Fermi gases have been a popular topic of research, with attention being paid recently to two-dimensional (2D) gases. In this work, we perform T=0 ab initio diffusion Monte Carlo calculations for a strongly interacting two-component Fermi gas confined to two dimensions. We first go over finite-size systems and the connection to the thermodynamic limit. After that, we illustrate pertinent 2D scattering physics and properties of the wave function. We then show energy results for the strong-coupling crossover, in between the Bose-Einstein condensation (BEC) and Bardeen-Cooper-Schrieffer (BCS) regimes. Our energy results for the BEC-BCS crossover are parametrized to produce an equation of state, which is used to determine Tan's contact. We carry out a detailed comparison with other microscopic results. Finally, we calculate the pairing gap for a range of interaction strengths in the strong coupling regime, following from variationally optimized many-body wave functions.
We will present recent experimental work on the thermodynamics of strongly interacting Fermi gases. We have developed a general method to probe with high precision the Equation of State (EoS) of locally homogeneous ultracold gases. This allows stringent tests of recent many-body theories. First, we focus on the finite-temperature EoS of the unpolarized unitary gas. Precise thermometry is provided by adding to the Fermi gas of 6Li a trace of bosonic 7Li. We show that the low-temperature properties of the strongly interacting normal phase are well described by Fermi liquid theory and we localize the superfluid transition. Second, we address the zero-temperature EoS of the spin-polarized system in the BEC-BCS crossover. Surprisingly, despite strong interactions, the polarized phase behaves as a mixture of two ideal gases: a Fermi gas of majority atoms and a non-interacting gas of dressed quasi-particles, the Fermi polarons. We also obtain the equation of state of the superfluid state as a function of interaction strength, that we compare to Monte Carlo simulations and to the Lee-Huang-Yang corrections for low-density bosonic and fermionic superfluids.
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.
We study the system of a dilute gas of fermions in 3-dimensions, with attractive interactions tuned to the unitarity point, using the non-perturbative Restricted Path Integral Monte Carlo (R-PIMC) method. The pairing and superfluid properties of this system are calculated at finite temperature. The total energy at very low temperature from our results agrees closely with that of previous ground-state Quantum Monte Carlo calculations. We identify the temperature T � � 0.70ǫF below which pairing correlations develop, and estimate the critical temperature for the superfluid transition Tc � 0.25ǫF from a finite size scaling analysis of the superfluid density. PACS numbers: The recent experiments, starting with generating a degenerate gas of cold atoms[3], the use of Feshbach resonance to tune the effective interaction between the fermions[4, 5], the measurements of the gap to identify a pairing scale[6] and the measurements of vortices[7] to identify the superfluid state, have all given an impetus to the study of superfluidity in the BCS-BEC crossover regime. The nature of superconductivity or superfluidity in a many-particle system with an increasing pairing attraction was shown [1, 2] to interpolate smoothly between the BCS regime for weak attraction between the fermions, and the Bose-Einstein condensation (BEC) regime for strong coupling. In the weak coupling regime, the Cooper pair size is much larger than the interparticle spacing, and the simultaneous pairing and condensation of fermions is well described by the BCS theory. In the strong coupling regime, fermions form strongly bound bosonic molecules at a pairing temperature scale T ∗ that is considerably higher than the BEC temperature Tc. In a two species system of fermions with an attractive interaction between them, the unitary point is defined by the divergence of the zero energy s-channel scattering length, as, for two particles in free space. The unitary point is interesting since it describes a strongly correlated system with universal properties[8]. In general, there are two length scales that define the system: the interparticle spacing ∝ n −1/3 , and the scattering length as which contains information about the interaction potential. At the unitary point however, since as diverges, all the properties of the system are described by a single length scale k −1