We present a numerical study of the micro-dynamical roots of dissipation in two colliding mesoscopic clouds of point-like fermions as a function of the scattering length and of temperature approaching full quantum degeneracy. This study, which is motivated by current experiments on ultracold gaseous mixtures of fermionic atoms inside magnetic traps, combines the solution of the coupled Vlasov-Landau equations for the Wigner distribution functions with a locally adaptive importance-sampling technique for handling collisional interactions. The results illustrate the consequences of genuinely quantum collisional phenomena, and in particular the role of Pauli blocking in the transition to hydrodynamic behaviour. We also compare the computed quantum collision rate as a function of temperature in the weak-coupling case with theoretical results assuming that equilibrium distributions determine the quantum collision integral.
We present a numerical study of the micro-dynamical roots of dissipation in two colliding mesoscopic clouds of point-like fermions as a function of the scattering length and of temperature approaching full quantum degeneracy. This study, which is motivated by current experiments on ultracold gaseous mixtures of fermionic atoms inside magnetic traps, combines the solution of the coupled Vlasov-Landau equations for the Wigner distribution functions with a locally adaptive importance-sampling technique for handling collisional interactions. The results illustrate the consequences of genuinely quantum collisional phenomena, and in particular the role of Pauli blocking in the transition to hydrodynamic behaviour. We also compare the computed quantum collision rate as a function of temperature in the weak-coupling case with theoretical results assuming that equilibrium distributions determine the quantum collision integral.
Experiments at JILA on the collisions between two oscillating spin-polarized components of a Fermi gas of 40K atoms have shown that this setup is an important tool for studying the dynamics of quantum gases. We have developed advanced numerical methods to solve the Vlasov–Landau equations for a binary mixture of atomic Fermi gases inside a harmonic trap at very low temperatures. In particular, a locally adaptive importance-sampling technique allows us to handle collisional interactions several orders of magnitude more efficiently than in standard Monte Carlo calculations. These methods have been applied to study the dynamical and transport properties of a fermion mixture in the presence of mean-field interactions and of collisions between its components for system parameters ranging from the collisionless to the hydrodynamic regime.
We present a theoretical study of the dynamical behavior of a gas made of ultracold fermionic atoms, which during their motions can collide with a much smaller number of thermal bosonic impurities. The atoms are confined inside harmonic traps and the interactions between the two species are treated as due to s-wave scattering with a negative scattering length modeling the $^{40}\mathrm{K}\penalty1000-\hskip0pt^{87}\mathrm{Rb}$ fermion-boson system. We set the fermions into motion by giving a small shift to their trap center and examine two alternative types of initial conditions, referring to (i) a close-to-equilibrium situation in which the two species are at the same temperature (well below the Fermi temperature and well above the Bose-Einstein condensation temperature); and (ii) a far-from-equilibrium case in which the impurities are given a Boltzmann distribution of momenta while the fermions are at very low temperatures. The dynamics of the gas is evaluated by the numerical solution of the Vlasov-Landau equations for the one-body distribution functions, supported by some analytical results on the collisional properties of a fermion gas. We find that the trapped gaseous mixture is close to the collisionless regime for values of the parameters corresponding to current experiments, but can be driven towards a collisional regime even without increasing the strength of the interactions, either by going over to heavier impurity masses or by matching the width of the momentum distribution of the impurities to the Fermi momentum of the fermion gas.
We study the transition from the collisionless to the hydrodynamic regime in a two-component spin-polarized mixture of 40 K atoms by exciting its dipolar oscillation modes inside harmonic traps. The time evolution of the mixture is described by the Vlasov-Landau equations and numerically solved with a fully three-dimensional concurrent code. We observe a master/slave behaviour of the oscillation frequencies depending on the dipolar mode that is excited. Regardless of the initial conditions, the transition to hydrodynamics is found to shift to lower values of the collision rate as the temperature decreases.
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.
We investigate how collisional interactions between the condensate and the thermal cloud influence the distinct dynamical regimes (Josephson plasma, phase-slip-induced dissipative regime and macroscopic quantum self-trapping) emerging in ultracold atomic Josephson junctions at non-zero subcritical temperatures. Specifically, we discuss how the self-consistent dynamical inclusion of collisional processes facilitating the exchange of particles between the condensate and the thermal cloud impacts both the condensate and the thermal currents, demonstrating that their relative importance depends on the system's dynamical regime. Our study is performed within the full context of the Zaremba-Nikuni-Griffin (ZNG) formalism, which couples a dissipative Gross-Pitaevskii equation for the condensate dynamics to a quantum Boltzmann equation with collisional terms for the thermal cloud. In the Josephson plasma oscillation and vortex-induced dissipative regimes, collisions markedly alter dynamics at intermediate-to-high temperatures, amplifying damping in the condensate imbalance mode and inducing measurable frequency shifts. In the self-trapping regime, collisions destabilize the system even at low temperatures, prompting a transition to Josephson-like dynamics on a temperature-dependent timescale. Our results show the interplay between coherence, dissipation, and thermal effects in a Bose-Einstein condensate at finite temperature, providing a framework for tailoring Josephson junction dynamics in experimentally accessible regimes.