A classic example of a quantum quench concerns the release of an interacting Bose gas from an optical lattice. The local properties of quenches such as this have been extensively studied; however, the global properties of these nonequilibrium quantum systems have received far less attention. Here we study several aspects of global nonequilibrium behavior by calculating the amount of work done by the quench as measured through the work distribution function. Using Bethe ansatz techniques, we determine the Loschmidt amplitude and work distribution function of the Lieb-Liniger gas after it is released from an optical lattice. We find the average work and its universal edge exponents from which we determine the long-time decay of the Loschmidt echo and highlight striking differences caused by the interactions as well as changes in the geometry of the system. We extend our calculation to the attractive regime of the model and show that the system exhibits properties similar to the super-Tonks-Girardeau gas. Finally, we examine the prominent role played by bound states in the work distribution and show that, with low probability, they allow for work to be extracted from the quench.
A classic example of a quantum quench concerns the release of an interacting Bose gas from an optical lattice. The local properties of quenches such as this have been extensively studied; however, the global properties of these nonequilibrium quantum systems have received far less attention. Here we study several aspects of global nonequilibrium behavior by calculating the amount of work done by the quench as measured through the work distribution function. Using Bethe ansatz techniques, we determine the Loschmidt amplitude and work distribution function of the Lieb-Liniger gas after it is released from an optical lattice. We find the average work and its universal edge exponents from which we determine the long-time decay of the Loschmidt echo and highlight striking differences caused by the interactions as well as changes in the geometry of the system. We extend our calculation to the attractive regime of the model and show that the system exhibits properties similar to the super-Tonks-Girardeau gas. Finally, we examine the prominent role played by bound states in the work distribution and show that, with low probability, they allow for work to be extracted from the quench.
We summarize in the present work exact results obtained for Tonks-Girardeau gases on one-dimensional optical lattices both for the ground state and nonequilibrium dynamics. On the theoretical side, impenetrable bosons offer the opportunity to study strongly interacting systems in one-dimensional lattices exactly, by means of the Jordan-Wigner transformation, and hence contribute to the topic of strong correlations at the center of interest in both condensed matter physics and quantum gases. This motivation is further enhanced by recent experimental realizations of such systems with ultracold atoms. After having shown their universal properties in equilibrium, we concentrate on their nonequilibrium dynamics. It will be shown that, starting from a pure Fock state, quasi-long-range correlations develop dynamically and lead to the formation of quasicondensates with a momentum determined by the underlying lattice. We expect this effect to be relevant for atom lasers with full control of the wavelength. Then, we will show that the free evolution of an initially confined Tonks-Girardeau gas leads to a momentum distribution that approaches at long times that of the equivalent fermionic system, giving rise to a bosonic gas with a Fermi edge, and hence a fermionization that can only be obtained out of equilibrium. Remarkably, although the momentum distribution function of the Tonks-Girardeau gas becomes equal to the one of the fermions, no loss in coherence is observed in the system, as reflected by a large occupation of eigenstates of the one-particle density matrix.
The sine-Gordon---equivalently, the massive Thirring---Hamiltonian is ubiquitous in low-dimensional physics, with applications that range from cold atom and strongly correlated systems to quantum impurities. We study here its nonequilibrium dynamics using the quantum quench protocol---following the system as it evolves under the sine-Gordon Hamiltonian from initial Mott-type states with large potential barriers. By means of the Bethe ansatz, we calculate exactly the Loschmidt amplitude, the fidelity, and work distribution characterizing these quenches for different values of the interaction strength. Some universal features are noted as well as an interesting duality relating quenches in different parameter regimes of the model.
We develop a general method for calculating the characteristic function of the work distribution of the Tonks-Girardeau (TG) gas confined to an arbitrary time-dependent trapping potential. We use the determinant representation of the many-body wave function to characterise the nonequilibrium thermodynamics of the TG gas and obtain exact and computationally tractable expressions in terms of Fredholm determinants for the mean work, the work probability distribution function, the nonadiabaticity parameter, and the Loschmidt amplitude. When applied to a harmonically trapped TG gas, our results for the mean work and the nonadiabaticity parameter reduce to those derived previously using an alternative approach. We next propose to use periodic modulation of the trap frequency in order to drive the system to highly non-equilibrium states by taking advantage of the phenomenon of parametric resonance. Under such driving protocol, the nonadiabaticity parameter may reach large values, which indicates a large amount of work being done on the system, as compared to sudden quench protocols considered previously. This scenario is realizable in ultracold atom experiments, aiding fundamental understanding of all thermodynamic properties of the system, in addition to facilitating possible optimization of high-efficiency quantum heat engine or quantum refrigeration cycles.
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.
This dissertation serves as a summary of my Ph.D. work numerically studying equilibrium and non-equilibrium properties of strongly-interacting one-dimensional (1D) boson systems. This work is motivated by the fact that 1D systems are realizable and highly controllable with ultracold atoms in optical lattice and atom chip experiments. We apply a recent worm algorithmic Monte Carlo approach developed for 1D continuous systems to study their equilibrium properties, both with and without an underlying lattice. We also apply an exact lattice approach based on the Bose-Fermi mapping to check our Monte Carlo results in the Tonks-Girardeau limit, and more importantly, to study far-from-equilibrium expansion dynamics of the systems.
We first study the scaling of one-particle correlations of the harmonically trapped Lieb-Liniger gas with changing temperature and interaction strength. Based on the universal behaviors of the density and momentum profiles, we are able to determine the effective parameters needed to fully characterize the system. We also find that the Tonks-Girardeau limit at low temperatures is the ideal regime for the experimental observation of the 1/k4 momentum tail.
An extra periodic lattice can drive the transition from superfluid to Mott insulator states. Exact and complete phase diagrams for such transitions are available only in the weak interacting and deep lattice limit, in which the system can be described using one-band Bose-Hubbard model. Beyond this limit, we use the worm algorithm in continuous space to map out the phase diagrams at various interaction strengths. We compare our phase diagrams with one-band Bose-Hubbard predictions and identify the regime where the one-band description breaks down. We introduce an inverse confined scattering solution to obtain effective Hubbard parameters, with which the Bose-Hubbard model provides correct results for strong interactions and deep lattices at unit filling.
In addition to the equilibrium properties, we also study the expansion dynamics of ultracold atoms in the hard-core limit. Experimentally, this is usually achieved by turning off confining potentials and letting atoms expand in optical lattices. Theoretical studies from initial ground states predicted the occurrence of fermionization of the momentum distribution after long expansion times. In addition, quasicondensation at finite momenta emerges when expanding from Mott insulating domains. Here, we develop a finite-temperature extension of the lattice approach for dynamics. We find the dynamical ferminoization of the momentum distributions at all temperatures. For expansion from initial Mott domains, we observe enhanced correlations reminiscent of dynamical quasicondensation. Surprisingly, we find the systems appear to cool down during the melting of the Mott domains. We use an emergent local Hamiltonian to understand these emergent phenomena.
We study the statistics of large deviations of the intensive work done in an interaction quench of a one-dimensional Bose gas with a large number N of particles, system size L, and fixed density. We consider the case in which the system is initially prepared in the noninteracting ground state and a repulsive interaction is suddenly turned on. For large deviations of the work below its mean value, we show that the large-deviation principle holds by means of the quench action approach. Using the latter, we compute exactly the so-called rate function and study its properties analytically. In particular, we find that fluctuations close to the mean value of the work exhibit a marked non-Gaussian behavior, even though their probability is always exponentially suppressed below it as L increases. Deviations larger than the mean value exhibit an algebraic decay whose exponent cannot be determined directly by large-deviation theory. Exploiting the exact Bethe ansatz representation of the eigenstates of the Hamiltonian, we calculate this exponent for vanishing particle density. Our approach can be straightforwardly generalized to quantum quenches in other interacting integrable systems.