Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Nonequilibrium dynamics of a one-dimensional Bose gas via the coordinate Bethe ansatz

Jan C. Zill

Published Jan 30, 2017• 0 citations
Publisher

Abstract

This thesis is concerned with non-equilibrium phenomena in interacting quantum many-body systems. Specifically, we investigate the time-evolution and relaxation dynamics of Bose gases in a highly restricted geometry, which constrains the dynamics to one spatial dimension. This leads to a description of the system in terms of a simple model Hamiltonian which permits exact many-body quantum mechanical solutions due to its integrability. In the first part of this thesis, a computational method is developed to obtain experimentally relevant correlation functions in the framework of the coordinate Bethe ansatz. We employ this method to compute exact ground-state correlation functions of the Lieb-Liniger gas for up to seven particles covering the whole regime of repulsive interactions. We also investigate the dynamics of the system after an instantaneous change of the interaction strength. This quantum quench deposits large amounts of energy that cannot be dissipated due to the system being closed, and so the dynamics far from equilibrium are probed. We prepare the system in two different initial states, and quench to the same final interaction strength. The latter is determined in such a way that the added energy due to the quench is the same for both scenarios. Conventional statistical mechanics predicts the same relaxed state, but due to the integrability of the system all considered correlation functions of the relaxed states differ from each other and also from the thermal ones. We then investigate the dynamics and relaxed state for a quench from zero to repulsive interactions in more detail, focussing on the mechanism of relaxation and the involved time-scales. We find that local correlation functions relax on time-scales determined by the interaction strength, in contrast to non-local correlation functions, whose relaxation time-scale is proportional to the system size. Next, we employ the same methodology to study the one-dimensional Bose gas with attractive interactions. In this case many-body bound states are permissible solutions of the Lieb-Liniger model. We compare exact ground-state correlation functions of up to seven particles to their corresponding mean-field solution. The latter displays a quantum phase transition at a critical interaction strength, marking the transition from a uniform-density state to a localized bright soliton. Our exact results agree remarkably well with the corresponding mean-field solution past the critical point. We also investigate the dynamics following an interaction strength quench, starting again from the non-interacting ground state. Bound states strongly influence correlation functions for all post-quench interaction strengths, and local correlation functions are largely increased compared to their initial value. In the last part of this thesis, we investigate the behavior of the one-dimensional Bose gas under periodic driving of the interaction strength. To this end, we extend the coordinate Bethe-ansatz formalism by employing Floquet theory to obtain solutions for the full time- dependent Hamiltonian. This realizes an extension to non-integrable systems in a way that allows controlled breaking of integrability. We compare the dynamics of the system following an interaction quench to a fixed final interaction strength to that of a system with periodically modulated post-quench interaction strength. For fast driving, the system evolves according to the time-averaged Hamiltonian and correlation functions are nearly identical to those of the undriven system. The response of the system qualitatively changes at resonances of the time-averaged Hamiltonian, where we observe energy absorption and a marked change in correlation functions. However, the system does not absorb energy indefinitely.

DOI
10.14264/uql.2017.182

Related papers

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

Jan C. Zill, Tod M. Wright, Karen V. Kheruntsyan +2 more

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 noninteracting 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.

arXiv:1601.00434•New Journal of Physics•Apr 13, 2016•34 citations

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

Jan C. Zill, Tod M. Wright, Karen V. Kheruntsyan +2 more

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.

arXiv:1407.4998•Physical Review A•Feb 12, 2015•43 citations

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

Jan Zill, Tod M. Wright, Karen V. Kheruntsyan +2 more

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.

arXiv:1705.09168•SciPost Physics•Feb 25, 2018•14 citations

A study of one-dimensional quantum gases

Andrew G. Sykes

In this thesis we study the physics of quantum many-body systems confined to one-dimensional geometries. The work was motivated by the recent success of experimentalists in developing atom traps, capable of restricting the motion of the individual atoms to a single spatial dimension. Specifically, we look at aspects of the one-dimensional Bose gas including; excitation spectrum, correlation functions, and dynamical behaviour. In Chapter \ref{ch:excitation1D} we consider the Lieb-Liniger model of interacting bosons in one-dimension. We numerically solve the equations arising from the Bethe ansatz solution for the exact many-body wave function in a finite-size system of up to twenty particles for attractive interactions. We discuss novel features of the solutions, including deviations from the well-known string solutions due to finite size effects. We present excited state string solutions in the limit of strong interactions and discuss their physical interpretation, as well as the characteristics of the quantum phase transition that occurs as a function of interaction strength in the mean-field limit. Our results are compared to those obtained via exact diagonalization of the Hamiltonian in a truncated basis. In Chapter \ref{ch:g2} we analytically calculate the spatial nonlocal pair correlation function for an interacting uniform one dimensional Bose gas at finite temperature and propose an experimental method to measure nonlocal correlations. Our results span six different physical realms, including the weakly and strongly interacting regimes. We show explicitly that the characteristic correlation lengths are given by one of four length scales: the thermal de Broglie wavelength, the mean interparticle separation, the healing length, or the phase coherence length. In all regimes, we identify the profound role of interactions and find that under certain conditions the pair correlation may develop a global maximum at a finite interparticle separation due to the competition between repulsive interactions and thermal effects. In Chapter \ref{ch:casimirdrag} we study the drag force below the critical velocity for obstacles moving in a superfluid. The absence of drag is well established in the context of the mean-field Gross-Pitaevskii theory. We calculate the next order correction due to quantum and thermal fluctuations and find a non-zero force acting on a delta-function impurity moving through a quasi-one-dimensional Bose-Einstein condensate at all subcritical velocities and at all temperatures. The force occurs due to an imbalance in the Doppler shifts of reflected quantum fluctuations from either side of the impurity. Our calculation is based on a consistent extension of Bogoliubov theory to second order in the interaction strength, and finds new analytic solutions to the Bogoliubov-de Gennes equations for a gray soliton. In Chapter \ref{ch:solitons} we study the effect of quantum noise on the stability of a soliton. We find the soliton solutions exactly define the reflectionless potentials of the Bogoliubov-de Gennes equations. This results in complete stability of the solitons in a purely one dimensional system. We look at the modifications to the density profile of a black soliton due to quantum fluctuations.

•Mar 1, 2010•0 citations

Quench dynamics of one dimensional multi-component quantum gases

Huijie Guan

This thesis studies the quench dynamics of strongly correlated quantum systems described by one dimensional integrable Hamiltonians. We develop the Yudson approach for such systems in terms of contour integrals allowing the expansion of arbitrary states in terms of the Bethe Ansatz eigenstates which in turn provides the means to calculate the time evolution of such observables as densities or noise correlations from any initial states. As a motivation to the present work, I present the state of art ultracold atom techniques and fundamental questions related to the quench dynamics. Then various Bethe ansatz solutions are discussed. This is followed by an introduction to the Yudson approach where advantages and difficulties are listed. As applications of the Yudson approach, I studied the nonequilibrium dynamics of the Lieb-Liniger model,(bosonic) Gaudin-Yang model, quenched from a Mott state with a super fluid Hamiltonian. Integration contours are specified for various models and different interactions. It is shown that the Yudson approach incorporates all free states and bound states into the contour, separating them apart sheds light on the validity of the String hypothesis. The result is affirmative for Lieb-Linger model, but not for the Gaudin-Yang model. Our calculations for the density and correlation shows that for interacting system, if the pre-quench state has negligible overlap among the particles, the system retains this feature after the quench. Particularly, normalized noise function c(z; -z) at the origin shows different stages. Shortly after the quench, the sign are different for attractive and repulsive interaction. Then they both quickly approaches the value where possibility to find both particles equals zero. Then the value increases gradually for attractive models, while c(0; 0) remains small for repulsive systems.

•Jan 1, 2018•0 citations

Exact time evolution of space- and time-dependent correlation functions after an interaction quench in the one-dimensional Bose gas

Jorn Mossel, Jean-Sébastien Caux

We consider the non-equilibrium dynamics of the interacting Lieb–Liniger gas after instantaneously switching the interactions off. The subsequent time evolution of the space- and time-dependent correlation functions is computed exactly. Different relaxation behavior is observed for different correlation functions. The long time average is compared with the predictions of several statistical ensembles. The generalized Gibbs ensemble restricted to a fixed number of particles is shown to give correct results at large times for all length scales.

arXiv:1201.1885•New Journal of Physics•Jul 3, 2012•104 citations
This thesis is concerned with non-equilibrium phenomena in interacting quantum many-body systems. Specifically, we investigate the time-evolution and relaxation dynamics of Bose gases in a highly restricted geometry, which constrains the dynamics to one spatial dimension. This leads to a description of the system in terms of a simple model Hamiltonian which permits exact many-body quantum mechanical solutions due to its integrability. In the first part of this thesis, a computational method is developed to obtain experimentally relevant correlation functions in the framework of the coordinate Bethe ansatz. We employ this method to compute exact ground-state correlation functions of the Lieb-Liniger gas for up to seven particles covering the whole regime of repulsive interactions. We also investigate the dynamics of the system after an instantaneous change of the interaction strength. This quantum quench deposits large amounts of energy that cannot be dissipated due to the system being closed, and so the dynamics far from equilibrium are probed. We prepare the system in two different initial states, and quench to the same final interaction strength. The latter is determined in such a way that the added energy due to the quench is the same for both scenarios. Conventional statistical mechanics predicts the same relaxed state, but due to the integrability of the system all considered correlation functions of the relaxed states differ from each other and also from the thermal ones. We then investigate the dynamics and relaxed state for a quench from zero to repulsive interactions in more detail, focussing on the mechanism of relaxation and the involved time-scales. We find that local correlation functions relax on time-scales determined by the interaction strength, in contrast to non-local correlation functions, whose relaxation time-scale is proportional to the system size. Next, we employ the same methodology to study the one-dimensional Bose gas with attractive interactions. In this case many-body bound states are permissible solutions of the Lieb-Liniger model. We compare exact ground-state correlation functions of up to seven particles to their corresponding mean-field solution. The latter displays a quantum phase transition at a critical interaction strength, marking the transition from a uniform-density state to a localized bright soliton. Our exact results agree remarkably well with the corresponding mean-field solution past the critical point. We also investigate the dynamics following an interaction strength quench, starting again from the non-interacting ground state. Bound states strongly influence correlation functions for all post-quench interaction strengths, and local correlation functions are largely increased compared to their initial value. In the last part of this thesis, we investigate the behavior of the one-dimensional Bose gas under periodic driving of the interaction strength. To this end, we extend the coordinate Bethe-ansatz formalism by employing Floquet theory to obtain solutions for the full time- dependent Hamiltonian. This realizes an extension to non-integrable systems in a way that allows controlled breaking of integrability. We compare the dynamics of the system following an interaction quench to a fixed final interaction strength to that of a system with periodically modulated post-quench interaction strength. For fast driving, the system evolves according to the time-averaged Hamiltonian and correlation functions are nearly identical to those of the undriven system. The response of the system qualitatively changes at resonances of the time-averaged Hamiltonian, where we observe energy absorption and a marked change in correlation functions. However, the system does not absorb energy indefinitely.
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 noninteracting 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.
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.
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.
In this thesis we study the physics of quantum many-body systems confined to one-dimensional geometries. The work was motivated by the recent success of experimentalists in developing atom traps, capable of restricting the motion of the individual atoms to a single spatial dimension. Specifically, we look at aspects of the one-dimensional Bose gas including; excitation spectrum, correlation functions, and dynamical behaviour. In Chapter \ref{ch:excitation1D} we consider the Lieb-Liniger model of interacting bosons in one-dimension. We numerically solve the equations arising from the Bethe ansatz solution for the exact many-body wave function in a finite-size system of up to twenty particles for attractive interactions. We discuss novel features of the solutions, including deviations from the well-known string solutions due to finite size effects. We present excited state string solutions in the limit of strong interactions and discuss their physical interpretation, as well as the characteristics of the quantum phase transition that occurs as a function of interaction strength in the mean-field limit. Our results are compared to those obtained via exact diagonalization of the Hamiltonian in a truncated basis. In Chapter \ref{ch:g2} we analytically calculate the spatial nonlocal pair correlation function for an interacting uniform one dimensional Bose gas at finite temperature and propose an experimental method to measure nonlocal correlations. Our results span six different physical realms, including the weakly and strongly interacting regimes. We show explicitly that the characteristic correlation lengths are given by one of four length scales: the thermal de Broglie wavelength, the mean interparticle separation, the healing length, or the phase coherence length. In all regimes, we identify the profound role of interactions and find that under certain conditions the pair correlation may develop a global maximum at a finite interparticle separation due to the competition between repulsive interactions and thermal effects. In Chapter \ref{ch:casimirdrag} we study the drag force below the critical velocity for obstacles moving in a superfluid. The absence of drag is well established in the context of the mean-field Gross-Pitaevskii theory. We calculate the next order correction due to quantum and thermal fluctuations and find a non-zero force acting on a delta-function impurity moving through a quasi-one-dimensional Bose-Einstein condensate at all subcritical velocities and at all temperatures. The force occurs due to an imbalance in the Doppler shifts of reflected quantum fluctuations from either side of the impurity. Our calculation is based on a consistent extension of Bogoliubov theory to second order in the interaction strength, and finds new analytic solutions to the Bogoliubov-de Gennes equations for a gray soliton. In Chapter \ref{ch:solitons} we study the effect of quantum noise on the stability of a soliton. We find the soliton solutions exactly define the reflectionless potentials of the Bogoliubov-de Gennes equations. This results in complete stability of the solitons in a purely one dimensional system. We look at the modifications to the density profile of a black soliton due to quantum fluctuations.
This thesis studies the quench dynamics of strongly correlated quantum systems described by one dimensional integrable Hamiltonians. We develop the Yudson approach for such systems in terms of contour integrals allowing the expansion of arbitrary states in terms of the Bethe Ansatz eigenstates which in turn provides the means to calculate the time evolution of such observables as densities or noise correlations from any initial states. As a motivation to the present work, I present the state of art ultracold atom techniques and fundamental questions related to the quench dynamics. Then various Bethe ansatz solutions are discussed. This is followed by an introduction to the Yudson approach where advantages and difficulties are listed. As applications of the Yudson approach, I studied the nonequilibrium dynamics of the Lieb-Liniger model,(bosonic) Gaudin-Yang model, quenched from a Mott state with a super fluid Hamiltonian. Integration contours are specified for various models and different interactions. It is shown that the Yudson approach incorporates all free states and bound states into the contour, separating them apart sheds light on the validity of the String hypothesis. The result is affirmative for Lieb-Linger model, but not for the Gaudin-Yang model. Our calculations for the density and correlation shows that for interacting system, if the pre-quench state has negligible overlap among the particles, the system retains this feature after the quench. Particularly, normalized noise function c(z; -z) at the origin shows different stages. Shortly after the quench, the sign are different for attractive and repulsive interaction. Then they both quickly approaches the value where possibility to find both particles equals zero. Then the value increases gradually for attractive models, while c(0; 0) remains small for repulsive systems.
We consider the non-equilibrium dynamics of the interacting Lieb–Liniger gas after instantaneously switching the interactions off. The subsequent time evolution of the space- and time-dependent correlation functions is computed exactly. Different relaxation behavior is observed for different correlation functions. The long time average is compared with the predictions of several statistical ensembles. The generalized Gibbs ensemble restricted to a fixed number of particles is shown to give correct results at large times for all length scales.