In order to understand transport phenomena in a quasi-classical regime the Boltzmann transport equation (BTE) is one of the most frequently used tools. Therein, the key quantity is a collision integral – the quantity that encapsulates the properties of the medium under consideration. Usually the result of this integral is approximated by one single parameter, the relaxation time. However it leaves one wondering if such situation is sufficient: for instance, if the dispersion of bosons is non-linear then what will be the influence of this non-linearity on the BTE. Here we give a fully analytic solution of the collisions’ integral for 1D bosonic gases with non-linear dispersion and far out of equilibrium. Our analytic result is given in terms of the Lerch transcendent function and it has been obtained for the case of two sub-systems (one dragging another), by taking a maximum-entropy displaced Bose–Einstein ansatz for their distributions. Currently, there are numerous experiments performed far away from equilibrium, where distributions are massively shifted and our result will serve as a main building block to derive distributions of bosons, and later linear and nonlinear transport coefficients, in such regimes.
In order to understand transport phenomena in a quasi-classical regime the Boltzmann transport equation (BTE) is one of the most frequently used tools. Therein, the key quantity is a collision integral – the quantity that encapsulates the properties of the medium under consideration. Usually the result of this integral is approximated by one single parameter, the relaxation time. However it leaves one wondering if such situation is sufficient: for instance, if the dispersion of bosons is non-linear then what will be the influence of this non-linearity on the BTE. Here we give a fully analytic solution of the collisions’ integral for 1D bosonic gases with non-linear dispersion and far out of equilibrium. Our analytic result is given in terms of the Lerch transcendent function and it has been obtained for the case of two sub-systems (one dragging another), by taking a maximum-entropy displaced Bose–Einstein ansatz for their distributions. Currently, there are numerous experiments performed far away from equilibrium, where distributions are massively shifted and our result will serve as a main building block to derive distributions of bosons, and later linear and nonlinear transport coefficients, in such regimes.
We study, in the framework of the Boltzmann-Bloch-Nordheim equation (BBN), the kinetic properties of a boson gas above the Bose-Einstein transition temperature Tc. The BBN equation is solved numerically within an algorithm that we test with exact analytical results for the collision rate of a homogeneous system in thermal equilibrium. In the classical regime $(T\ensuremath{\gg}{T}_{c}),$ the relaxation time of a quadrupolar deformation in momentum space is proportional to the mean free collision time ${\ensuremath{\tau}}_{\mathrm{relax}}\ensuremath{\propto}{\ensuremath{\tau}}_{\mathrm{coll}}\ensuremath{\propto}{T}^{\ensuremath{-}1/2}.$ Approaching the critical temperature, quantum statistic effects become dominant, and the collision rate increases dramatically. Nevertheless, this does not affect the relaxation properties of the gas that depend only on the spontaneous collision term in BBN. The relaxation time becomes ${\ensuremath{\tau}}_{\mathrm{relax}}\ensuremath{\propto}(T\ensuremath{-}{T}_{c}{)}^{\ensuremath{-}1/2},$ exhibiting a critical slowing down. These phenomena can be experimentally confirmed looking at the damping properties of collective motions induced on trapped atoms. The possibility to observe a transition from collisionless (zero-sound) to hydrodynamic (such as first and second-sound) is finally discussed.
This thesis examines two aspects of the non-equilibrium dynamics of Bose-Einstein condensates (BECs). Our particular interests are the relaxation of near-integrable systems, and the inter-particle correlations present in BECs, both in equilibrium and following an abrupt change in the Hamiltonian. In the first part of this thesis we investigate the stability of solitary wave solutions to a system of two coupled one-dimensional non-linear Schrodinger equations (CNLSE). This system has an integrable point known as the Manakov[1] model which has analytic dark-bright soliton solutions[2]. We break the integrability of this model by varying the inter-species interaction strength, and employ variational and numerical methods in order to find dark-bright solitary wave solutions, away from the integrable point. The time evolution of these states is calculated via numerical integration of the CNLSE, which allows us to assess the stability of these solutions in isolation and to quantify their robustness against collisions with a dark soliton. We find that there is a broad region of the parameter space in which solutions for black-bright (stationary) solitary waves can be found. In this domain, the integrability breaking is revealed only during collisions. Prior to a collision, the numerical solutions are indistinguishable from true solitons, however, the collisional stability of these solitary waves is significantly affected by the extent to which integrability is broken. We find that there is a smooth transition between the non-dispersive particle-like nature of soliton interactions for integrable systems, and the destructive collisions which take place when the system is perturbed far from the integrable point. In the case of moving (grey-bright) solitary waves, we find a more restricted region of the parameter space which admits solitary wave solutions. In this domain we are able to find approximately stable solitary wave solutions using a variational ansatz. The collisional stability of these grey-bright solitary waves is also found to depend upon the extent to which integrability is broken. We conclude that the observation of stable, long-lived dark-bright solitary waves is not sufficient to indicate near-integrability of the two-component system. However, near-integrability may inferred if dark-bright solitary waves are observed to survive multiple collisions. These results are not only of theoretical interest, but may also inform future experiments which explore the behaviour of solitary waves in two-component Bose-Einstein condensates. Following this we investigate the properties of magnetic solitary wave solutions to the CNLSE. We find that the analytic ‘magnetic soliton’ solutions derived in reference [3] are dynamically unstable and relax to non-stationary magnetic solitary wave states over a short timescale. We investigate the collisional dynamics of these magnetic solitary waves and find that these excitations are remarkably robust. In the second part of this thesis we calculate the fluctuations in the number difference between two sub-regions of a Bose-Einstein condensate. This study is motivated by recent experiments performed by the Truscott group at the Australian National University (ANU), which indicate sub- Poissonian fluctuations even for condensates with large thermal depletion. This result cannot easily be reconciled with existing literature, which predicts strongly super-Poissonian fluctuations below the critical temperature (Tc). We develop one-dimensional, and three-dimensional, models to describe these experiments and employ the Bogoliubov formalism to explore the quantum fluctuations of the condensate, which have not previously been studied in detail. We calculate the relative number fluctuations between two halves of the condensate, and also between multiple pairs of spatial bins. For harmonically confined condensates at zero temperature, quantum correlations result in a suppression of the relative number fluctuations between the two halves of the condensate below the shot-noise level. However, the presence of even a small thermal component results in strongly super-Poissonian fluctuations below Tc; this is in agreement with the existing theoretical and experimental results. In addition, we find that the fluctuations of the harmonically confined condensate are non-uniform, with the relative number fluctuations being peaked towards the edges of the condensate. The recent studies of the Truscott group measure the fluctuations after free expansion. We therefore calculate the propagation of correlations following a sudden removal of the harmonic confinement. We find that free expansion results in a strong decrease in the relative number fluctuations. Condensates which exhibit super-Poissonian fluctuations when harmonically confined may therefore display number- squeezing after expansion. A second set of experiments by the Truscott group measured the density fluctuations present in small ‘atom laser’ pulses of atoms out-coupled from their BEC. We also model this procedure and find that immediately after the out-coupling process, the number fluctuations of the out-coupled atomic pulse are approximately Poissonian, irrespective of the correlations present in the atomic reservoir. However, the competing effects of interactions between the atomic pulse and the condensate reservoir, and free expansion of the pulse after out-coupling, can result in the measured fluctuations being either number squeezed or super-Poissonian. Our results will guide future experiments planned by the Truscott group.
The Navier-Stokes equations are paradigmatic equations describing hydrodynamics of an interacting system with microscopic interactions encoded in transport coefficients. In this work we show how the Navier-Stokes equations arise from the microscopic dynamics of nearly integrable 1D quantum many-body systems. We build upon the recently developed hydrodynamics of integrable models to study the effective Boltzmann equation with collision integral taking into account the non-integrable interactions. We compute the transport coefficients and find that the resulting Navier-Stokes equations have two regimes, which differ in the viscous properties of the fluid. We illustrate the method by computing the transport coefficients for an experimentally relevant case of coupled 1D cold-atomic gases.
The Boltzmann equation is a powerful theoretical tool for modeling the collective dynamics of quantum many-body systems subject to external perturbations. Analysis of the equation gives access to linear response properties including collective modes and transport coefficients, but often proves intractable due to computational costs associated with multidimensional integrals describing collision processes. Here, we present a method to resolve this bottleneck, enabling the study of a broad class of many-body systems that appear in fundamental science contexts and technological applications. Specifically, we demonstrate that a Gaussian mixture model can accurately represent equilibrium distribution functions, thereby allowing efficient evaluation of collision integrals. Inspired by cold atom experiments, we apply this method to investigate the collective behavior of a quantum Bose-Fermi mixture of cold atoms in a cigar-shaped trap, a system that is particularly challenging to analyze. We focus on monopole and quadrupole collective modes above the Bose-Einstein transition temperature, and find a rich phenomenology that spans interference effects between bosonic and fermionic collective modes, dampening of these modes, and the emergence of hydrodynamics in various parameter regimes. These effects are readily verifiable experimentally. Published by the American Physical Society 2024
We treat the non-equilibrium evolution of an open one-particle statistical system, subject to a potential and to an external “heat bath” (hb) with negligible dissipation. For the classical equilibrium Boltzmann distribution, Wc,eq, a non-equilibrium three-term hierarchy for moments fulfills Hermiticity, which allows one to justify an approximate long-time thermalization. That gives partial dynamical support to Boltzmann’s Wc,eq, out of the set of classical stationary distributions, Wc;st, also investigated here, for which neither Hermiticity nor that thermalization hold, in general. For closed classical many-particle systems without hb (by using Wc,eq), the long-time approximate thermalization for three-term hierarchies is justified and yields an approximate Lyapunov function and an arrow of time. The largest part of the work treats an open quantum one-particle system through the non-equilibrium Wigner function, W. Weq for a repulsive finite square well is reported. W’s (< 0 in various cases) are assumed to be quasi-definite functionals regarding their dependences on momentum (q). That yields orthogonal polynomials, HQ,n(q), for Weq (and for stationary Wst), non-equilibrium moments, Wn, of W and hierarchies. For the first excited state of the harmonic oscillator, its stationary Wst is a quasi-definite functional, and the orthogonal polynomials and three-term hierarchy are studied. In general, the non-equilibrium quantum hierarchies (associated with Weq) for the Wn’s are not three-term ones. As an illustration, we outline a non-equilibrium four-term hierarchy and its solution in terms of generalized operator continued fractions. Such structures also allow one to formulate long-time approximations, but make it more difficult to justify thermalization. For large thermal and de Broglie wavelengths, the dominant Weq and a non-equilibrium equation for W are reported: the non-equilibrium hierarchy could plausibly be a three-term one and possibly not far from Gaussian, and thermalization could possibly be justified.
In Chapter 17, we derived two-fluid hydrodynamic equations that include damping related to transport coefficients. Our entire analysis was based on the coupled ZNG equations for the condensate and in the thermal cloud. These involved a generalized GP equation for the condensate and a kinetic equation for the thermal atoms. A crucial role is played by the C 12 collision term in the kinetic equation, which describes the interactions between atoms in the condensate and in the thermal cloud. Our analysis of the deviation from the diffusive local equilibrium solution of the kinetic equation was based on the Chapman–Enskog approach, extensively developed for classical gases and first applied to Bose-condensed gases by Kirkpatrick and Dorfman (1983, 1985a). This approach required a careful treatment of the novel feature relating to the C 12 collisions both in the kinetic equation describing the thermal atoms and also in the source term Γ 12 in the generalized GP equation for the condensate. Using the Chapman– Enskog approach to solve the kinetic equation for a trapped Bose gas, we obtained explicit expressions for the function ψ( p, r , t ) that describes the deviation from diffusive local equilibrium, as defined by (17.25) and (17.39). This deviation can be related to various transport coefficients, as discussed in Chapter 17. These transport coefficients are determined by the solutions of the three integral equations (17.40)–(17.42) for the three contributions to the deviation function ψ( p, r , t ) in (17.39). In Section 18.1, we will solve these integral equations and obtain explicit expressions for the thermal conductivity k , the shear viscosity η and the four second viscosity coefficients ζ i .