We use numerical simulation to study ultracold, quantum degenerate, atomic gases. In the first part of the thesis we study the effects of disorder, introduced via a bichromatic optical lattice, in one and two dimensional systems. We employ the Aubry-Andr\'{e} model and use time-dependent numerical simulations to investigate the disorder dependent transition to strong localisation present in the model. Weak s-wave interactions are added to the model and we observe the interaction between localisation and interaction induced self-trapping. We then add a tilted lattice potential to the model. In the homogeneous model this induces Bloch oscillations. While one might expect that a strong enough force will break the strong localisation or self-trapping, within the bounds of the single-band model, the trapping effect of the Bloch oscillations reinforces both of the other effects leading to increased confinement, albeit lacking the clear single frequency oscillation signature of pure Bloch oscillations. Along with the two dimensional bichromatic optical lattice we add a term to the Hamiltonian equivalent to that of a uniform external magnetic field on charged particles. Since the experimental realisation of this model would employ neutral atoms, the magnetic field is synthetic, the equivalent effect being produced by an appropriate set of lasers and magnetic fields. We show that in the ballistic regime (weak bichromatic disorder) the system displays positive magnetoresistance. Conversely in the strong localisation regime the system exhibits negative magnetoresistance. In the latter part of the thesis we use density functional theory to calculate the ground-state density of a harmonically trapped dipolar Fermi gas. We then use these to calculate the lowest energy collective mode oscillation frequencies under the hydrodynamic approximation. We find that increasing the strength of the dipoles has the effect of increasing the mode frequencies. The increase saturates for large dipole strengths. We verify this analytically and show that such is due to the local nature of the two dimensional energy functional and not dependent on the specific equation of state. We employ an average density approximation to construct an energy functional for the inhomogeneous, 2D degenerate Fermi gas. The ground-state densities for a cylindrically symmetric harmonic trap are compared to the Kohn-Sham results, showing extremely good agreement in the tail region and good agreement with the exact ground-state energy. We then do the same for higher order polynomial traps and obtain improved agreement for higher degree.
We use numerical simulation to study ultracold, quantum degenerate, atomic gases. In the first part of the thesis we study the effects of disorder, introduced via a bichromatic optical lattice, in one and two dimensional systems. We employ the Aubry-Andr'{e} model and use time-dependent numerical simulations to investigate the disorder dependent transition to strong localisation present in the model. Weak s-wave interactions are added to the model and we observe the interaction between localisation and interaction induced self-trapping. We then add a tilted lattice potential to the model. In the homogeneous model this induces Bloch oscillations. While one might expect that a strong enough force will break the strong localisation or self-trapping, within the bounds of the single-band model, the trapping effect of the Bloch oscillations reinforces both of the other effects leading to increased confinement, albeit lacking the clear single frequency oscillation signature of pure Bloch oscillations. Along with the two dimensional bichromatic optical lattice we add a term to the Hamiltonian equivalent to that of a uniform external magnetic field on charged particles. Since the experimental realisation of this model would employ neutral atoms, the magnetic field is synthetic, the equivalent effect being produced by an appropriate set of lasers and magnetic fields. We show that in the ballistic regime (weak bichromatic disorder) the system displays positive magnetoresistance. Conversely in the strong localisation regime the system exhibits negative magnetoresistance. In the latter part of the thesis we use density functional theory to calculate the ground-state density of a harmonically trapped dipolar Fermi gas. We then use these to calculate the lowest energy collective mode oscillation frequencies under the hydrodynamic approximation. We find that increasing the strength of the dipoles has the effect of increasing the mode frequencies. The increase saturates for large dipole strengths. We verify this analytically and show that such is due to the local nature of the two dimensional energy functional and not dependent on the specific equation of state. We employ an average density approximation to construct an energy functional for the inhomogeneous, 2D degenerate Fermi gas. The ground-state densities for a cylindrically symmetric harmonic trap are compared to the Kohn-Sham results, showing extremely good agreement in the tail region and good agreement with the exact ground-state energy. We then do the same for higher order polynomial traps and obtain improved agreement for higher degree.
In this thesis we study the effects of the dipole-dipole interaction on the physics of ultracold quantum gases, both bosonic and fermionic, within the theoretical framework provided by the mean-field regime. This kind of interaction takes place in ultracold atomic gases (for instance 52Cr or 164Dy) due to their atomic magnetic dipole moment, and in ultracold molecular gases due to the magnetic or electric dipole moment. \n \nIn the case of quantum gases of bosonic atoms, or Bose-Einstein condensates, the dipole-dipole interaction can be studied within mean-field approximation using the Gross-Pitaevskii equation, which now contains a new non-linear term due to the dipole-dipole interaction. We investigate, on the one hand dipolar condensates confined in harmonic traps, and on the other dipolar condensates confined in toroidal traps. In the harmonic geometry, our focus is on the study of the ground state and the quantized vortex state, where the density profile is characterized as well as some properties leading to the process of vortex formation, such as the critical frequency and the energy barrier that has to be overcome to bring the vortex from the surface to the centre of the gas. We finish the study of dipolar condensates in harmonic traps by dynamically simulating the precession frequency of an off-center vortex in a non-rotating condensate. \n \nIn the toroidal geometry the dipolar effects are strongly magnified when the polarization axis of the dipoles is perpendicular to the trap symmetry axis. In this case, the anisotropic structure of the density can be understood as the response of the system to the double-well effective potential along the ring. We have studied the dynamics of this system when the initial number of atoms in the left and right wells is imbalanced, predicting Josephson and self-trapping oscillations depending on the initial condition. This has led us to name this new system as Self-induced Josephson Junction. We have studied in detail the self-trapping regime and we have seen that the particle flux inversion is closely related to the crossing of vortices across the Josephson junctions. This result opens the door to establishing a more direct connection between the phase-slip regime, widely addressed in superfluid helium, and the self-trapping regime of condensates. \n \nIn the case of quantum gases of fermionic dipolar particles, we have studied how the radial quadrupole mode allows one to distinguish between hydrodynamic and collisionless regimes. We have analytically calculated the frequency of this mode in the mean-field approximation, generalizing the results from the Thomas-Fermi approximation for trapped ideal Fermi gases. On the one hand, we observe that the frequency in the hydrodynamic regime is smaller than in non-dipolar Fermi gases, while in the collisionless regime the frequency is larger or smaller than that corresponding to the non-interacting system depending on the geometry of the harmonic trap. On the other hand, we predict that reducing the trap deformation (aspect ratio) an observable jump in the frequency of the radial quadrupole mode would take place, which would correspond to the transition between the collisionless and hydrodynamic regimes, for instance when the gas undergoes the transition to the superfluid state.
This thesis contains theoretical research on ultracold quantum gases in spatially periodic potentials, featuring high-frequency periodic driving, long-range interactions, or both. The largest part features deep potentials where the behaviour of the gas is well-described by quantum lattice models. The periodic driving is then integrated out to obtain effective time-independent descriptions. One project investigates emergent long-range interactions in a stationary, weak, spatially periodic potential, where a lattice theory is not appropriate.In two short introductory chapters, the topics of ultracold atoms, optical lattice potentials, periodic driving, and long-range interactions, are sketched, without attempting to give a complete overview. Some experimentally relevant length and energy scales are given, but the main focus is on deriving and constructing theoretical descriptions of ultracold gases in various spatially and/or temporally periodic potentials.The first project presented, is focused on how the single-particle spectrum of a Bose gas in a non-separable two-dimensional square lattice is affected by high-frequency periodic driving; the most striking conclusion is that under suitable circumstances, it develops two inequivalent minima, leading to finite-momentum Bose-Einstein condensation. Perturbative calculations indicate that local interactions induce spontaneous time-reversal symmetry breaking (TRSB) in such a system.The second project investigates the interplay between kinetic frustration and long-range interactions in fermionic gases. Both a mean-field approximation and exact diagonalisations predict that such a system, studied in the more specific realisation of a weakly interacting dipolar fermionic gas in a 2D triangular lattice, also leads to spontaneous TRSB. Furthermore, a density wave could form at quarter filling, where the Fermi surface is perfectly nested. Perhaps more interesting yet, a spatially inhomogeneous TRSB pattern is predicted, confined to the low-density sublattice that emerges in the density wave.The third project revolves around the question of supersolidity in the presence of a gauge field. Applying the Bogolyubov approximation to a variation of the extended Bose-Hubbard model, indicates that combining an artificial staggered magnetic field in a 2D square lattice with nearest-neighbour density-density interactions, not only leads to a supersolid with staggered vortices, but also induces an inhomogeneous distribution of the associated currents around the elementary plaquette.In the fourth project, a one-dimensional Bose gas with strong local interactions in a weak lattice at incommensurate densities is shown to feature excitations corresponding to excess or deficit particles. The excitations interact repulsively at long distances, in spite of the fact that the underlying atoms themselves do not. As a consequence, the incommensurability of the density with the lattice can drive a transition to a density wave and even a supersolid.The four above-mentioned research projects combine bosonic and fermionic gases, weak and strong interactions, perturbative and mean-field approximations, effective field theories and exact diagonalisations. The main overall conclusion is that long-range interactions and high-frequency periodic driving lead to a very diverse range of fascinating phenomena in ultracold lattice gases.
In this thesis we address relevant problems of the physics of quantum disordered systems from a numerical and theoretical point of view, with specific attention to the connection of our findings with ultracold atomic gases experiments. We concentrate on two main issues: the interplay between localization and interaction in disordered systems and the problem of localization in correlated random potentials. The first problem is investigated considering the expansion of a weakly interacting Bose gas in a bichromtic optical lattice. We observe that interaction has a destructive effect on the disorder-induced localization and leads to a subdiffusive expansion of the atomic gas. By comparing three characteristic energy scales of the system one can identify three different spreading regimes: weak chaos, strong chaos and self-trapping. The spreading behaviour in these regimes is predicted theoretically and verified numerically. We also interpreted existing experimental data on the basis of our findings and showed that there is a qualitative agreement between our numerical simulations and experiments. The second problem is investigated proposing a new model of correlated disorder that can be implemented experimentally using ultracold dipolar gases. We show that this model is characterized by the presence of both short and long range correlations. We study the localization properties of the model and highlight the role played by short and long range correlations in the determination of those properties. In particular we show that when short-range correlations are dominant, extended states can appear in the spectrum. The effect of long-range correlations is instead to restore localization over the whole spectrum and lead to counterintuitive behaviours of the localization length. More precisely, depending on the localization regime they can enhance or reduce the localization length at the centre of the band.
The experimental achievement of Bose-Einstein condensation and Fermi degeneracy with ultracold gases boosted tremendous progresses both in theoretical methods and in the development of new experimental tools. Among them, intriguing possibilities have been opened by the implementation of optical lattices: periodic potentials for neutral atoms created by interfering laser beams. Degenerate gases in optical lattices can be forced in highly anisotropic traps, reducing the effective dimensionality of the system. From a fundamental point of view, the behavior of matter waves in reduced dimensions sheds light on the intimate properties of interparticle interactions. Furthermore, such reduced-dimensional systems can be engineered to quantum-simulate fasci- nating solid state systems, like bidimensional crystals, in a clean and controllable environment. Motivated by the exciting perspectives of this field, we devote this Thesis to the theoretical study of two systems where matter waves propagate in reduced dimensions. The long-range and anisotropic character of the dipole-dipole interaction critically affects the behavior of dipolar quantum gases. The continuous experimental progresses in this flourishing field might lead very soon to the creation of degenerate dipolar gases in optical potentials. In the first part of this Thesis, we investigate the emergence of a single dipolar-induced resonance in the two-body scattering process in quasi-one dimensional geometries. We develop a two-channel approach to describe such a resonance in a highly elongated cigar-shaped harmonic trap, which approximates the single site of a quasi-one-dimensional optical lattice. At this stage, we develop a novel atom-dimer extended Bose-Hubbard model for dipolar bosons in this quasi-one-dimensional optical lattice. Hence we investigate the T = 0 phase diagram of the model by exact diagonalization of a small-sized system, highlighting the effects of the dipolar-induced resonance on the many-body behavior in the lattice. In the second part of the Thesis, we present a general scheme to realize cold-atom quantum simulators of bidimensional atomic crystals, based on the possibility to independently trap two different atomic species. The first one constitutes a two-dimensional matter wave which interacts only with the atoms of the second species, deeply trapped around the nodes of a two-dimensional optical lattice. By introducing a general analytic approach, we investigate the matter-wave transport properties. We propose some illustrative applications to both Bravais (square, triangular) and non-Bravais (graphene, kagomeÌ ) lattices, studying both ideal periodic systems and experimental- sized, eventually disordered, ones. The features of the artificial atomic crystal critically depend on the two-body interspecies interaction strength, which is shown to be widely tunable via 0D-2D mixed-dimensional resonances. Keywords: matter waves, reduced dimensions, dipolar-induced resonances, mixed-dimensional resonances, extended Bose-Hubbard model, atomic artificial crystals.
The study of collective behaviour in many-body systems often explores fundamentally new ideas absent from the mere constituents of such a system. A paradigmatic model for these studies is the spin-1/2 XXZ chain and its fermionic equivalent. This thesis can be broadly divided into the study of two fundamental aspects of this model. Firstly, we discuss localisation phenomena in one dimensional lattices as often experimentally realised in cold atom systems. Secondly, we investigate how disorder and symmetry influence heat transport in spin chains. More specifically, in the first part we consider a system of non-interacting fermions in one dimension subject to a single-particle potential consisting of a strong optical lattice, a harmonic trap, and uncorrelated on-site disorder. We investigate a global inhomogeneous quantum quench and present numerical and analytical results for static and dynamical properties. We show that the approach to the non-thermal equilibrium state is extremely slow and that it implies a sensitivity to disorder parametrically stronger than that expected from Anderson localisation. We also consider the above system in a strong non-uniform electric field. In the non-interacting case, due to Wannier-Stark localisation, the single-particle wave functions are exponentially localised without quenched disorder. We show that this system remains localised in the presence of nearest-neighbour interactions and exhibits physics analogous to models of conventional many-body localisation. The second part explores the hydrodynamics of the disordered XYZ spin chain. Using time-evolving block decimation on open chains of up to 400 spins attached to thermal baths, we probe the energy transport of this system. Our principal findings are as follows. For weak disorder there is a stable diffusive region that persists up to a critical disorder strength that depends on the XY anisotropy. Then, for disorder strengths above this critical value energy transport becomes increasingly subdiffusive.