Boltzmann equation simulation for a trapped Fermi gas of atoms

The dynamics of an interacting Fermi gas of atoms at sufficiently high temperatures can be efficiently studied via a numerical simulation of the Boltzmann equation. In this work we describe in detail the setup we used recently to study the oscillations of two spin-polarised fermionic clouds in a trap. We focus here on the evaluation of interparticle interactions. We compare different ways of choosing the phase space coordinates of a pair of atoms after a successful collision and demonstrate that the exact microscopic setup has no influence on the macroscopic outcome

Thermodynamic and hydrodynamic behaviour of interacting Fermi gases

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.This work was supported by the German Academic Exchange Service (DAAD), the Engineering and Physical Sciences Research Council (EPSRC) and the Cambridge European Trust. I am grateful to the Department of Applied Mathematics and Theoretical Physics (DAMTP), Wolfson College, the Cambridge Philosophical society, the ECT* in Trento, Durham University, the ENS in Paris and the Institute for Nuclear Theory at the University of Washington for providing me with travel grants and local financial support during conferences and meetings. This work has made use of the resources provided by the University of Cambridge High Performance Computing Service (http://www.hpc.cam.ac.uk/)

Analog Simulation of Weyl Particles with Cold Atoms

We study theoretically, numerically, and experimentally the relaxation of a collisionless gas in a quadrupole trap after a momentum kick. The non-separability of the potential enables a quasi thermalization of the single particle distribution function even in the absence of interactions. Suprinsingly, the dynamics features an effective decoupling between the strong trapping axis and the weak trapping plane. The energy delivered during the kick is redistributed according to the symmetries of the system and satisfies the Virial theorem, allowing for the prediction of the final temperatures. We show that this behaviour is formally equivalent to the relaxation of massless relativistic Weyl fermions after a sudden displacement from the center of a harmonic trap

Globe-hopping

We consider versions of the grasshopper problem (Goulko & Kent 2017 Proc. R. Soc. A473, 20170494) on the circle and the sphere, which are relevant to Bell inequalities. For a circle of circumference 2π, we show that for unconstrained lawns of any length and arbitrary jump lengths, the supremum of the probability for the grasshopper’s jump to stay on the lawn is one. For antipodal lawns, which by definition contain precisely one of each pair of opposite points and have length π, we show this is true except when the jump length ϕ is of the form π(p/q) with p, q coprime and p odd. For these jump lengths, we show the optimal probability is 1 − 1/q and construct optimal lawns. For a pair of antipodal lawns, we show that the optimal probability of jumping from one onto the other is 1 − 1/q for p, q coprime, p odd and q even, and one in all other cases. For an antipodal lawn on the sphere, it is known (Kent & Pitalúa-García 2014 Phys. Rev. A90, 062124) that if ϕ = π/q, where q∈N, then the optimal retention probability of 1 − 1/q for the grasshopper’s jump is provided by a hemispherical lawn. We show that in all other cases where 0 < ϕ < π/2, hemispherical lawns are not optimal, disproving the hemispherical colouring maximality hypotheses (Kent & Pitalúa-García 2014 Phys. Rev. A90, 062124). We discuss the implications for Bell experiments and related cryptographic tests