Tod M. Wright

Nonequilibrium dynamics of one-dimensional hard-core anyons following a quench: complete relaxation of one-body observables.

We demonstrate the role of interactions in driving the relaxation of an isolated integrable quantum system following a sudden quench. We consider a family of integrable hard-core lattice anyon models that continuously interpolates between noninteracting spinless fermions and strongly interacting hard-core bosons. A generalized Jordan-Wigner transformation maps the entire family to noninteracting fermions. We find that, aside from the singular free-fermion limit, the entire single-particle density matrix and, therefore, all one-body observables relax to the predictions of the generalized Gibbs ensemble (GGE). This demonstrates that, in the presence of interactions, correlations between particles in the many-body wave function provide the effective dissipation required to drive the relaxation of all one-body observables to the GGE. This relaxation does not depend on translational invariance or the tracing out of any spatial domain of the system.

Single-particle and many-body analyses of a quasiperiodic integrable system after a quench

In general, isolated integrable quantum systems have been found to relax to an apparent equilibrium state in which the expectation values of few-body observables are described by the generalized Gibbs ensemble. However, recent work has shown that relaxation to such a generalized statistical ensemble can be precluded by localization in a quasiperiodic lattice system. Here we undertake complementary single-particle and many-body analyses of noninteracting spinless fermions and hard-core bosons within the Aubry-Andre model to gain insight into this phenomenon. Our investigations span both the localized and delocalized regimes of the quasiperiodic system, as well as the critical point separating the two. Considering first the case of spinless fermions, we study the dynamics of the momentum distribution function and characterize the effects of real-space and momentum-space localization on the relevant single-particle wave functions and correlation functions. We show that although some observables do not relax in the delocalized and localized regimes, the observables that do relax in these regimes do so in a manner consistent with a recently proposed Gaussian equilibration scenario, whereas relaxation at the critical point has a more exotic character. We also construct various statistical ensembles from the many-body eigenstates of the fermionic and bosonic Hamiltonians and study the effect of localization on their properties.

A coordinate Bethe ansatz approach to the calculation of equilibrium and nonequilibrium correlations of the one-dimensional Bose gas

We use the coordinate Bethe ansatz to exactly calculate matrix elements between eigenstates of the Lieb-Liniger model of one-dimensional bosons interacting via a two-body delta-potential. We investigate the static correlation functions of the zero-temperature ground state and their dependence on interaction strength, and analyze the effects of system size in the crossover from few-body to mesoscopic regimes for up to seven particles. We also obtain time-dependent nonequilibrium correlation functions for five particles following quenches of the interaction strength from two distinct initial states. One quench is from the non-interacting ground state and the other from a correlated ground state near the strongly interacting Tonks-Girardeau regime. The final interaction strength and conserved energy are chosen to be the same for both quenches. The integrability of the model highly constrains its dynamics, and we demonstrate that the time-averaged correlation functions following quenches from these two distinct initial conditions are both nonthermal and moreover distinct from one another.

Dynamical thermalization and vortex formation in stirred two-dimensional Bose-Einstein condensates

We present a quantum-mechanical treatment of the mechanical stirring of Bose-Einstein condensates using classical field techniques. In our approach the condensate and excited modes are described using a Hamiltonian classical field method in which the atom number and (rotating frame) energy are strictly conserved. We simulate a T=0 quasi-two-dimensional condensate perturbed by a rotating anisotropic trapping potential. Vacuum fluctuations in the initial state provide an irreducible mechanism for breaking the initial symmetries of the condensate and seeding the subsequent dynamical instability. Highly turbulent motion develops and we quantify the emergence of a rotating thermal component that provides the dissipation necessary for the nucleation and motional damping of vortices in the condensate. Vortex lattice formation is not observed, rather the vortices assemble into a spatially disordered vortex liquid state. We discuss methods we have developed to identify the condensate in the presence of an irregular distribution of vortices, determine the thermodynamic parameters of the thermal component, and extract damping rates from the classical field trajectories.

Relaxation dynamics of the Lieb-Liniger gas following an interaction quench: a coordinate Bethe-ansatz analysis

We investigate the relaxation dynamics of the integrable Lieb-Liniger model of contact-interacting bosons in one dimension following a sudden quench of the collisional interaction strength. The system is initially prepared in its noninteracting ground state and the interaction strength is then abruptly switched to a positive value, corresponding to repulsive interactions between the bosons. We calculate equal-time correlation functions of the nonequilibrium Bose field for small systems of up to five particles via symbolic evaluation of coordinate Bethe-ansatz expressions for operator matrix elements between Lieb-Liniger eigenstates. We characterize the relaxation of the system by comparing the time-evolving correlation functions following the quench to the equilibrium correlations predicted by the diagonal ensemble and relate the behavior of these correlations to that of the quantum fidelity between the many-body wave function and the initial state of the system. Our results for the asymptotic scaling of local second-order correlations with increasing interaction strength agree with the predictions of recent generalized thermodynamic Bethe-ansatz calculations. By contrast, third-order correlations obtained within our approach exhibit a markedly different power-law dependence on the interaction strength as the Tonks-Girardeau limit of infinitely strong interactions is approached.

Many-body physics in the classical-field description of a degenerate Bose gas

The classical-field formalism has been widely applied in the calculation of normal correlation functions, and the characterization of condensation, in finite-temperature Bose gases. Here we discuss the extension of this method to the calculation of more general correlations, including the so-called anomalous correlations of the field, without recourse to symmetry-breaking assumptions. Our method is based on the introduction of U(1)-symmetric classical-field variables analogous to the modified quantum ladder operators of number-conserving approaches to the degenerate Bose gas, and allows us to rigorously quantify the anomalous and non-Gaussian character of the field fluctuations. We compare our results for anomalous correlation functions with the predictions of mean-field theories, and demonstrate that the nonlinear classical-field dynamics incorporate a full description of many-body processes which modify the effective mean-field potentials experienced by condensate and noncondensate atoms. We discuss the role of these processes in shaping the condensate mode, and thereby demonstrate the consistency of the Penrose-Onsager definition of the condensate orbital in the classical-field equilibrium. We consider the contribution of various noncondensate-field correlations to the overall suppression of density fluctuations and interactions in the field, and demonstrate the distinct roles of phase and density fluctuations in the transition of the field to the normal phase.

Finite-temperature dynamics of a single vortex in a Bose-Einstein condensate: Equilibrium precession and rotational symmetry breaking

We consider a finite-temperature Bose-Einstein condensate in a quasi-two-dimensional trap containing a single precessing vortex. We find that such a configuration arises naturally as an ergodic equilibrium of the projected Gross-Pitaevskii equation, when constrained to a finite conserved angular momentum. In an isotropic trapping potential, the condensation of the classical field into an off-axis vortex state breaks the rotational symmetry of the system. We present a methodology to identify the condensate and the Goldstone mode associated with the broken rotational symmetry in the classical-field model. We also examine the variation in vortex trajectories and thermodynamic parameters of the field as the energy of the microcanonical field simulation is varied.

Quantum quench dynamics of the attractive one-dimensional Bose gas via the coordinate Bethe ansatz

We use the coordinate Bethe ansatz to study the Lieb-Liniger model of a one-dimensional gas of bosons on a finite-sized ring interacting via an attractive delta-function potential. We calculate zero-temperature correlation functions for seven particles in the vicinity of the crossover to a localized solitonic state and study the dynamics of a system of four particles quenched to attractive interactions from the ideal-gas ground state. We determine the time evolution of correlation functions, as well as their temporal averages, and discuss the role of bound states in shaping the postquench correlations and relaxation dynamics.

Formation of Bose-Einstein condensates

The problem of understanding how a coherent, macroscopic Bose-Einstein condensate (BEC) emerges from the cooling of a thermal Bose gas has attracted significant theoretical and experimental interest over several decades. The pioneering achievement of BEC in weakly-interacting dilute atomic gases in 1995 was followed by a number of experimental studies examining the growth of the BEC number, as well as the development of its coherence. More recently there has been interest in connecting such experiments to universal aspects of nonequilibrium phase transitions, in terms of both static and dynamical critical exponents. Here, the spontaneous formation of topological structures such as vortices and solitons in quenched cold-atom experiments has enabled the verification of the Kibble-Zurek mechanism predicting the density of topological defects in continuous phase transitions, first proposed in the context of the evolution of the early universe. This chapter reviews progress in the understanding of BEC formation, and discusses open questions and future research directions in the dynamics of phase transitions in quantum gases.

C-Field Methods for Non-Equilibrium Bose Gases

We review c-field methods for simulating the non-equilibrium dynamics of degenerate Bose gases beyond the mean-field Gross-Pitaevskii approximation. We describe three separate approaches that utilise similar numerical methods, but have distinct regimes of validity. Systems at finite temperature can be treated with either the closed-system projected Gross-Pitaevskii equation (PGPE), or the open-system stochastic projected Gross-Pitaevskii equation (SPGPE). These are both applicable in quantum degenerate regimes in which thermal fluctuations are significant. At low or zero temperature, the truncated Wigner projected Gross-Pitaevskii equation (TWPGPE) allows for the simulation of systems in which spontaneous collision processes seeded by quantum fluctuations are important. We describe the regimes of validity of each of these methods, and discuss their relationships to one another, and to other simulation techniques for the dynamics of Bose gases. The utility of the SPGPE formalism in modelling non-equilibrium Bose gases is illustrated by its application to the dynamics of spontaneous vortex formation in the growth of a Bose-Einstein condensate.