Focus on Out-of-Equilibrium Dynamics in Strongly Interacting One-Dimensional Systems

In the past few years, there have been significant advances in understanding out-of-equilibrium dynamics in strongly interacting many-particle quantum systems. This is the case for 1D dynamics, where experimental advances—both with ultracold atomic gases and with solid state systems—have been accompanied by advances in theoretical methods, both analytical and numerical. This 'focus on' collection brings together 17 new papers, which together give a representative overview of the recent advances.

We use a broadband microbolometer array to measure the full three-dimensional (3D) terahertz (THz) intensity profile emitted from a two-color femtosecond plasma and subsequently focused in a geometry useful for nonlinear spectroscopic investigations. Away from the immediate focal region we observe a sharp, conical intensity profile resembling a donut, and in the focal region the beam collapses to a central, Lorentz-shaped profile. The Lorentzian intensity profile in the focal region can be explained by considering the frequency-dependent spot size derived from measurements of the Gouy phase shift in the focal region, and the transition from the donut profile to a central peak is consistent with propagation of a Bessel–Gauss beam, as shown by simulations based on a recent transient photocurrent model (You et al 2012 Phys. Rev. Lett. 109 183902). We combine our measurements to the first full 3D visualization of the conical THz emission from the two-color plasma.

We detail the experimental observation of the non-equilibrium many-body phenomenon prethermalization. We study the dynamics of a rapidly and coherently split one-dimensional Bose gas. An analysis based on the use of full quantum mechanical probability distributions of matter wave interference contrast reveals that the system evolves toward a quasi-steady state. This state, which can be characterized by an effective temperature, is not the final thermal equilibrium state. We compare the evolution of the system to an integrable Tomonaga–Luttinger liquid model, and show that the system dephases to a prethermalized state rather than undergoing thermalization toward a final thermal equilibrium state.

We study the non-equilibrium dynamics of an impurity in a harmonic trap that is kicked with a well-defined quasi-momentum, and interacts with a bath of free fermions or interacting bosons in a one-dimensional lattice configuration. Using numerical and analytical techniques we investigate the full dynamics beyond linear response, which allows us to quantitatively characterize states of the impurity in the bath for different parameter regimes. These vary from a tightly bound molecular state in a strongly interacting limit to a polaron (dressed impurity) and a free particle for weak interactions, with composite behaviour in the intermediate regime. These dynamics and different parameter regimes should be readily realizable in systems of cold atoms in optical lattices.

In the study of trapped two-component Bose gases, a widely used dynamical protocol is to start from the ground state of a one-component condensate and then switch half the atoms into another hyperfine state. The slightly different intra-component and inter-component interactions can then lead to highly non-trivial dynamics, especially in the density mismatch between the two components, commonly referred to as 'spin' density. We study and classify the possible subsequent dynamics, over a wide variety of parameters spanned by the trap strength and by the inter- to intra-component interaction ratio. A stability analysis suited to the trapped situation provides us with a framework to explain the various types of dynamics in different regimes.

Theoretical understanding of strongly correlated systems in one spatial dimension has been greatly advanced by the density-matrix renormalization group (DMRG) algorithm, which is a variational approach using a class of entanglement-restricted states called matrix product states (MPSs). However, DMRG suffers from inherent accuracy restrictions when multiple states are involved due to multi-state targeting and also the approximate representation of the Hamiltonian and other operators. By formulating the variational approach of DMRG explicitly for MPSs one can avoid errors inherent in the multi-state targeting approach. Furthermore, by using the matrix product operator (MPO) formalism, one can exactly represent the Hamiltonian and other operators relevant for the calculation. The MPO approach allows one-dimensional (1D) Hamiltonians to be templated using a small set of finite state automaton rules without reference to the particular microscopic degrees of freedom. We present two algorithms which take advantage of these properties: eMPS to find excited states of 1D Hamiltonians and tMPS for the time evolution of a generic time-dependent 1D Hamiltonian. We properly account for time-ordering of the propagator such that the error does not depend on the rate of change of the Hamiltonian. Our algorithms are applicable to microscopic degrees of freedom of any variety, and do not require Hamiltonian-specialized implementation. We benchmark our algorithms with a case study of the Ising model, where the critical point is located using entanglement measures. We then study the dynamics of this model under a time-dependent quench of the transverse field through the critical point. Finally, we present studies of a dipolar, or long-range Ising model, again using entanglement measures to find the critical point and study the dynamics of a time-dependent quench through the critical point.

The interplay between the confinement potential and the electron–electron interactions causes reconstructions of quantum Hall edges. We study the consequences of this edge reconstruction for the relaxation of hot electrons injected into integer quantum Hall edge states. In translationally invariant edges, the relaxation of hot electrons is governed by three-body collisions, which are sensitive to the electron dispersion and thus to reconstruction effects. We show that the relaxation rates are significantly altered in different reconstruction scenarios.

Continuous phase transitions occur in a wide range of physical systems, and provide a context for the study of non-equilibrium dynamics and the formation of topological defects. The Kibble–Zurek (KZ) mechanism predicts the scaling of the resulting density of defects as a function of the quench rate through a critical point, and this can provide an estimate of the critical exponents of a phase transition. In this work, we extend our previous study of the miscible–immiscible phase transition of a binary Bose–Einstein condensate (BEC) composed of two hyperfine states in which the spin dynamics are confined to one dimension (Sabbatini et al 2011 Phys. Rev. Lett. 107 230402). The transition is engineered by controlling a Hamiltonian quench of the coupling amplitude of the two hyperfine states and results in the formation of a random pattern of spatial domains. Using the numerical truncated Wigner phase-space method, we show that in a ring BEC the number of domains formed in the phase transitions scales as predicted by the KZ theory. We also consider the same experiment performed with a harmonically trapped BEC, and investigate how the density inhomogeneity modifies the dynamics of the phase transition and the KZ scaling law for the number of domains. We then make use of the symmetry between inhomogeneous phase transitions in anisotropic systems, and an inhomogeneous quench in a homogeneous system, to engineer coupling quenches that allow us to quantify several aspects of inhomogeneous phase transitions. In particular, we quantify the effect of causality in the propagation of the phase transition front on the resulting formation of domain walls and find indications that the density of defects is determined during the impulse to adiabatic transition after the crossing of the critical point.

We present a general formalism describing the propagation of an arbitrary multi-particle wave packet in a one-dimensional multimode chiral channel coupled to an ensemble of emitters which are distributed at arbitrary positions. The formalism is based on a direct and exact resummation of a diagrammatic series for the multi-particle scattering matrix. It is complementary to the Bethe Ansatz approach and to approaches based on equations of motion, and it reveals a simple and transparent structure of scattering states. In particular, we demonstrate how this formalism works on various examples, including the scattering of one- and two-photon states off two- and three-level emitters, off an array of emitters, as well as the scattering of coherent light. We argue that this formalism can be used constructively for the study of the scattering of an arbitrary initial photonic state off emitters with an arbitrary degree of complexity.

This paper investigates numerically a phenomenon which can be used to transport a single q-bit down a J₁–J₂ Heisenberg spin chain using a quantum adiabatic process. The motivation for investigating such processes comes from the idea that this method of transport could potentially be used as a means of sending data to various parts of a quantum computer made of artificial spins, and that this method could take advantage of the easily prepared ground state at the so-called Majumdar–Ghosh point. We examine several annealing protocols for this process and find similar results for all of them. The annealing process works well up to a critical frustration threshold. There is also a brief section examining what other models this protocol could be used for, examining its use in the XXZ and XYZ models.

We investigate the relation between thermalization following a quantum quench and many-body localization in quasi-particle space in terms of the long-time full distribution function of physical observables. In particular, expanding on our recent work (Canovi et al 2011 Phys. Rev. B 83 094431), we focus on the long-time behavior of an integrable XXZ chain subject to an integrability-breaking perturbation. After a characterization of the breaking of integrability and the associated localization/delocalization transition using the level spacing statistics and the properties of the eigenstates, we study the effect of integrability breaking on the asymptotic state after a quantum quench of the anisotropy parameter, looking at the behavior of the full probability distribution of the transverse and longitudinal magnetization of a subsystem. We compare the resulting distributions with those obtained in equilibrium at an effective temperature set by the initial energy. We find that, while the long-time distribution functions appear to always agree qualitatively with the equilibrium ones, quantitative agreement is obtained only when integrability is fully broken and the relevant eigenstates are diffusive in quasi-particle space.

We study the dynamics of two types of pairs of excitations which are bound despite their strong repulsive interaction. One corresponds to doubly occupied sites in one-dimensional Bose–Hubbard systems, the so-called doublons. The other is pairs of neighboring excited spins in anisotropic Heisenberg spin-1/2 chains. We investigate the possibility of decay of the bound pairs due to resonant scattering by a defect or due to collisions of the pairs. We find that the amplitudes of the corresponding transitions are very small. This is the result of destructive quantum interference and explains the stability of the bound pairs.

It is analytically shown that the asymptotic correlations following a quantum quench in exactly solvable models can sometimes look essentially thermal provided the initial coupling between the system eigenmodes induces a large gap. We study this phenomenon using simple models, which also illustrate the relationship between the entanglement spectrum of the initial state and the generalized Gibbs ensemble describing the long-time correlations after the quench. We also show that the effective temperature characterizing the correlations is not related to the energy fluctuations after the quench, and therefore does not have thermodynamic meaning. The latter observation implies a breakdown of the fluctuation–dissipation theorem.

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.

Single-particle momentum spectra for a dynamically evolving one-dimensional Bose gas are analysed in the semi-classical wave limit. Representing one of the simplest correlation functions, these provide information on a possible universal scaling behaviour. Motivated by the previously discovered connection between (quasi-) topological field configurations, strong wave turbulence and non-thermal fixed points of quantum field dynamics, soliton formation is studied with respect to the appearance of transient power-law spectra. A random-soliton model is developed for describing the spectra analytically, and the analogies and differences between the emerging power laws and those found in a field theory approach to strong wave turbulence are discussed. The results open a new perspective on solitary wave dynamics from the point of view of critical phenomena far from thermal equilibrium and the possibility of studying this dynamics by experiment without the need for detecting solitons in situ.

Stochastic exclusion processes play an integral role in the physics of non-equilibrium statistical mechanics. These models are Markovian processes, described by a classical master equation. In this paper, a quantum mechanical version of a stochastic hopping process in one dimension is formulated in terms of a quantum master equation. This allows the investigation of coherent and stochastic evolution in the same formal framework. The focus lies on the non-equilibrium steady state. Two stochastic model systems are considered: the totally asymmetric exclusion process and the fully symmetric exclusion process. The steady-state transport properties of these models are compared to the case with additional coherent evolution, generated by the XX Hamiltonian.

We analyze the recently developed folding algorithm (Bañuls et al 2009 Phys. Rev. Lett. 102 240603) for simulating the dynamics of infinite quantum spin chains and we relate its performance to the kind of entanglement produced under the evolution of product states. We benchmark the accomplishments of this technique with respect to alternative strategies using Ising Hamiltonians with transverse and parallel fields, as well as XY models. Also, we evaluate its capability of finding ground and thermal equilibrium states.

We study under-barrier tunneling for a pair of energetically bound bosonic atoms in an optical lattice with a barrier. We identify the conditions under which this exotic molecule tunnels as a point particle with the coordinate given by the bound pair center of mass and discuss the atomic co-tunneling beyond this regime. In particular, we quantitatively analyze resonantly enhanced co-tunneling, where two interacting atoms penetrate the barrier with higher probability than a single atom.

We study the relaxation dynamics of the one-dimensional Tomonaga–Luttinger model after an interaction quench, paying particular attention to the momentum dependence of the two-particle interaction. Several potentials of different analytical forms are investigated that all lead to universal Luttinger liquid (LL) physics in equilibrium. The steady-state fermionic momentum distribution shows universal behavior in the sense of the LL phenomenology. For generic regular potentials, the large time decay of the momentum distribution function toward the steady-state value is characterized by a power law with a universal exponent that depends only on the potential at zero momentum transfer. The commonly employed ad hoc procedure fails to give this exponent. In addition to quenches from zero to positive interactions, we also consider the abrupt changes in the interaction between two arbitrary values. Additionally, we discuss the appearance of a factor of two between the steady-state momentum distribution function and that obtained in equilibrium at equal two-particle interaction.