Exploring local quantum many-body relaxation by atoms in optical superlattices

Effective light cone and digital quantum simulation of interacting bosons

The speed limit of information propagation is one of the most fundamental features in non-equilibrium physics. The region of information propagation by finite-time dynamics is approximately restricted inside the effective light cone that is formulated by the Lieb-Robinson bound. To date, extensive studies have been conducted to identify the shape of effective light cones in most experimentally relevant many-body systems. However, the Lieb-Robinson bound in the interacting boson systems, one of the most ubiquitous quantum systems in nature, has remained a critical open problem for a long time. This study reveals a tight effective light cone to limit the information propagation in interacting bosons, where the shape of the effective light cone depends on the spatial dimension. To achieve it, we prove that the speed for bosons to clump together is finite, which in turn leads to the error guarantee of the boson number truncation at each site. Furthermore, we applied the method to provide a provably efficient algorithm for simulating the interacting boson systems. The results of this study settle the notoriously challenging problem and provide the foundation for elucidating the complexity of many-body boson systems.

Adiabatic processes like isothermal processes

The objective of this work is to show that adiabatic processes can be very similar to isothermal ones. First, we show that the criteria for the compatibility of linear-response theory with the second law of thermodynamics for thermally isolated systems are the same as those for systems performing isothermal processes. Motivated by such a result, we explore the thermodynamic consequences of the time-average excess work, observing an unexpected existence of a well-defined relaxation time for thermally isolated systems that obeys the second law of thermodynamics. This is justified by recognizing that such systems, in the usual sense, present random relaxation time, which can be "averaged" by taking the time average of the relaxation function. Such a proceeding is very similar to what happens in isothermal processes, where a stochastic average must be done on the relaxation function to have a well-defined relaxation time. In the end, we analyze the Landau-Zener model from this new point of view, discussing the construction of slowly-varying processes from linear-response theory and observing negative entropy production rates for nonmonotonic and rapid protocols.

Observation of Thermalization and Information Scrambling in a Superconducting Quantum Processor

Understanding various phenomena in nonequilibrium dynamics of closed quantum many-body systems, such as quantum thermalization, information scrambling, and nonergodic dynamics, is crucial for modern physics. Using a ladder-type superconducting quantum processor, we perform analog quantum simulations of both the XX-ladder model and the one-dimensional XX model. By measuring the dynamics of local observables, entanglement entropy, and tripartite mutual information, we signal quantum thermalization and information scrambling in the XX ladder. In contrast, we show that the XX chain, as free fermions on a one-dimensional lattice, fails to thermalize to the Gibbs ensemble, and local information does not scramble in the integrable channel. Our experiments reveal ergodicity and scrambling in the controllable qubit ladder, and open the door to further investigations on the thermodynamics and chaos in quantum many-body systems.

Lieb-Robinson Bound and Almost-Linear Light Cone in Interacting Boson Systems

In this work, we investigate how quickly local perturbations propagate in interacting boson systems with Bose-Hubbard-type Hamiltonians. In general, these systems have unbounded local energies, and arbitrarily fast information propagation may occur. We focus on a specific but experimentally natural situation in which the number of bosons at any one site in the unperturbed initial state is approximately limited. We rigorously prove the existence of an almost-linear information-propagation light cone, thus establishing a Lieb-Robinson bound: the wave front grows at most as t log^{2}(t). We prove the clustering theorem for gapped ground states and study the time complexity of classically simulating one-dimensional quench dynamics, a topic of great practical interest.

Observation of Strong and Weak Thermalization in a Superconducting Quantum Processor

We experimentally study the ergodic dynamics of a 1D array of 12 superconducting qubits with a transverse field, and identify the regimes of strong and weak thermalization with different initial states. We observe convergence of the local observable to its thermal expectation value in the strong-thermalizaion regime. For weak thermalization, the dynamics of local observable exhibits an oscillation around the thermal value, which can only be attained by the time average. We also demonstrate that the entanglement entropy and concurrence can characterize the regimes of strong and weak thermalization. Our work provides an essential step toward a generic understanding of thermalization in quantum systems.

Entanglement-Ergodic Quantum Systems Equilibrate Exponentially Well

One of the outstanding problems in nonequilibrium physics is to precisely understand when and how physically relevant observables in many-body systems equilibrate under unitary time evolution. General equilibration results show that equilibration is generic provided that the initial state has overlap with sufficiently many energy levels. But results not referring to typicality which show that natural initial states actually fulfill this condition are lacking. In this work, we present stringent results for equilibration for systems in which Rényi entanglement entropies in energy eigenstates with finite energy density are extensive for at least some, not necessarily connected, subsystems. Our results reverse the logic of common arguments, in that we derive equilibration from a weak condition akin to the eigenstate thermalization hypothesis, which is usually attributed to thermalization in systems that are assumed to equilibrate in the first place. We put the findings into the context of studies of many-body localization and many-body scars.

Exact relaxation dynamics and quantum information scrambling in multiply quenched harmonic chains

The quantum dynamics of isolated systems under quench condition exhibits a variety of interesting features depending on the integrable/chaotic nature of system. We study the exact dynamics of trivially integrable system of harmonic chains under a multiple quench protocol. Out of time ordered correlator of two Hermitian operators at large time displays scrambling in the thermodynamic limit. In this limit, the entanglement entropy and the central component of momentum distribution both saturate to a steady-state value. We also show that reduced density matrix assumes the diagonal form long after multiple quenches for large system size. These exact results involving infinite-dimensional Hilbert space are indicative of dynamical equilibration for a trivially integrable harmonic chain.

Thermalization and Heating Dynamics in Open Generic Many-Body Systems

The last decade has witnessed remarkable progress in our understanding of thermalization in isolated quantum systems. Combining the eigenstate thermalization hypothesis with quantum measurement theory, we extend the framework of quantum thermalization to open many-body systems. A generic many-body system subject to continuous observation is shown to thermalize at a single trajectory level. We show that the nonunitary nature of quantum measurement causes several unique thermalization mechanisms that are unseen in isolated systems. We present numerical evidence for our findings by applying our theory to specific models that can be experimentally realized in atom-cavity systems and with quantum gas microscopy. Our theory provides a general method to determine an effective temperature of quantum many-body systems subject to the Lindblad master equation and thus should be applicable to noisy dynamics or dissipative systems coupled to nonthermal Markovian environments as well as continuously monitored systems. Our work provides yet another insight into why thermodynamics emerges so universally.

Quantum Relaxation and Metastability of Lattice Bosons with Cavity-Induced Long-Range Interactions

The coupling of cold atoms to the radiation field within a high-finesse optical resonator, an optical cavity, induces long-range interactions which can compete with an underlying optical lattice. The interplay between short- and long-range interactions gives rise to new phases of matter including supersolidity (SS) and density waves (DW), and interesting quantum dynamics. Here it is shown that for hard-core bosons in one dimension the ground state phase diagram and the quantum relaxation after sudden quenches can be calculated exactly in the thermodynamic limit. Remanent DW order is observed for quenches from a DW ground state into the superfluid (SF) phase below a dynamical transition line. After sufficiently strong SF to DW quenches beyond a static metastability line DW order emerges on top of remanent SF order, giving rise to a dynamically generated SS state. Our method to handle infinite- and short-range interactions in the infinite system size limit opens a way to solve exactly other Hamiltonians with infinite- and short-range interactions as well.

Analytic model of thermalization: Quantum emulation of classical cellular automata

We introduce a method of quantum emulation of a classical reversible cellular automaton. By applying this method to a chaotic cellular automaton, the obtained quantum many-body system thermalizes while all the energy eigenstates and eigenvalues are solvable. These explicit solutions allow us to verify the validity of some scenarios of thermalization to this system. We find that two leading scenarios, the eigenstate thermalization hypothesis scenario and the large effective dimension scenario, do not explain thermalization in this model.

Systematic Construction of Counterexamples to the Eigenstate Thermalization Hypothesis

We propose a general method to embed target states into the middle of the energy spectrum of a many-body Hamiltonian as its energy eigenstates. Employing this method, we construct a translationally invariant local Hamiltonian with no local conserved quantities, which does not satisfy the eigenstate thermalization hypothesis. The absence of eigenstate thermalization for target states is analytically proved and numerically demonstrated. In addition, numerical calculations of two concrete models also show that all the energy eigenstates except for the target states have the property of eigenstate thermalization, from which we argue that our models thermalize after a quench even though they do not satisfy the eigenstate thermalization hypothesis.

Eigenstate thermalization in the two-dimensional transverse field Ising model. II. Off-diagonal matrix elements of observables

We study the matrix elements of few-body observables, focusing on the off-diagonal ones, in the eigenstates of the two-dimensional transverse field Ising model. By resolving all symmetries, we relate the onset of quantum chaos to the structure of the matrix elements. In particular, we show that a general result of the theory of random matrices, namely, the value 2 of the ratio of variances (diagonal to off-diagonal) of the matrix elements of Hermitian operators, occurs in the quantum chaotic regime. Furthermore, we explore the behavior of the off-diagonal matrix elements of observables as a function of the eigenstate energy differences and show that it is in accordance with the eigenstate thermalization hypothesis ansatz.

Typical Relaxation of Isolated Many-Body Systems Which Do Not Thermalize

We consider isolated many-body quantum systems which do not thermalize; i.e., expectation values approach an (approximately) steady longtime limit which disagrees with the microcanonical prediction of equilibrium statistical mechanics. A general analytical theory is worked out for the typical temporal relaxation behavior in such cases. The main prerequisites are initial conditions which appreciably populate many energy levels and do not give rise to significant spatial inhomogeneities on macroscopic scales. The theory explains very well the experimental and numerical findings in a trapped-ion quantum simulator exhibiting many-body localization, in ultracold atomic gases, and in integrable hard-core boson and XXZ models.

Edge mode dynamics of quenched topological wires

The fermionic and Majorana edge mode dynamics of various topological systems are compared, after a sudden global quench of the Hamiltonian parameters takes place. Attention is focused on the regimes where the survival probability of an edge state has oscillations either due to critical or off-critical quenches. The nature of the wave functions and the overlaps between the eigenstates of different points in parameter space determine the various types of behaviors, and the distinction due to the Majorana nature of the excitations plays a lesser role. Performing a sequence of quenches, it is shown that the edge states, including Majorana modes, may be switched off and on. Also, the generation of Majoranas due to quenching from a trivial phase is discussed.

Typical fast thermalization processes in closed many-body systems

The lack of knowledge about the detailed many-particle motion on the microscopic scale is a key issue in any theoretical description of a macroscopic experiment. For systems at or close to thermal equilibrium, statistical mechanics provides a very successful general framework to cope with this problem. However, far from equilibrium, only very few quantitative and comparably universal results are known. Here a quantum mechanical prediction of this type is derived and verified against various experimental and numerical data from the literature. It quantitatively describes the entire temporal relaxation towards thermal equilibrium for a large class (in a mathematically precisely defined sense) of closed many-body systems, whose initial state may be arbitrarily far from equilibrium.

The relaxation of closed macroscopic systems towards thermal equilibrium is an ubiquitous experimental fact, but very difficult to characterize theoretically. Here, the author establishes a quantitative description of such relaxation under arbitrary typical conditions, capturing well experimental data.

Fate of Majorana fermions and Chern numbers after a quantum quench

In the sequence of quenches to either nontopological phases or other topological phases, we study the stability of Majorana fermions at the edges of a two-dimensional topological superconductor with spin-orbit coupling and in the presence of a Zeeman term. Both instantaneous and slow quenches are considered. In the case of instantaneous quenches, the Majorana modes generally decay, but for a finite system there is a revival time that scales to infinity as the system size grows. Exceptions to this decaying behavior are found in some cases due to the presence of edge states with the same momentum in the final state. Quenches to a topological Z(2) phase reveal some robustness of the Majorana fermions in the sense that even though the survival probability of the Majorana state is small, it does not vanish. If the pairing is not aligned with the spin-orbit Rashba coupling, it is found that the Majorana fermions are fairly robust with a finite survival probability. It is also shown that the Chern number remains invariant after the quench, until the propagation of the mode along the transverse direction reaches the middle point, beyond which the Chern number fluctuates between increasing values. The effect of varying the rate of change in slow quenches is also analyzed. It is found that the defect production is nonuniversal and does not follow the Kibble-Zurek scaling with the quench rate, as obtained before for other systems with topological edge states.

Quench dynamics of one-dimensional interacting bosons in a disordered potential: elastic dephasing and critical speeding-up of thermalization

The dynamics of interacting bosons in one dimension following the sudden switching on of a weak disordered potential is investigated. On time scales before quasiparticles scatter (prethermalized regime), the dephasing from random elastic forward scattering causes all correlations to decay exponentially fast, but the system remains far from thermal equilibrium. For longer times, the combined effect of disorder and interactions gives rise to inelastic scattering and to thermalization. A novel quantum kinetic equation accounting for both disorder and interactions is employed to study the dynamics. Thermalization turns out to be most effective close to the superfluid-Bose-glass critical point where nonlinearities become more and more important. The numerically obtained thermalization times are found to agree well with analytic estimates.

Expansion dynamics of interacting bosons in homogeneous lattices in one and two dimensions

We experimentally and numerically investigate the expansion of initially localized ultracold bosons in homogeneous one- and two-dimensional optical lattices. We find that both dimensionality and interaction strength crucially influence these nonequilibrium dynamics. While the atoms expand ballistically in all integrable limits, deviations from these limits dramatically suppress the expansion and lead to the appearance of almost bimodal cloud shapes, indicating diffusive dynamics in the center surrounded by ballistic wings. For strongly interacting bosons, we observe a dimensional crossover of the dynamics from ballistic in the one-dimensional hard-core case to diffusive in two dimensions, as well as a similar crossover when higher occupancies are introduced into the system.

Relaxation dynamics of disordered spin chains: localization and the existence of a stationary state

We study the unitary relaxation dynamics of disordered spin chains following a sudden quench of the Hamiltonian. We give analytical arguments, corroborated by specific numerical examples, to show that the existence of a stationary state depends crucially on the spectral and localization properties of the final Hamiltonian, and not on the initial state. We test these ideas on integrable one-dimensional models of the Ising or XY class, but argue more generally on their validity for more complex (nonintegrable) models.

Constructing the generalized Gibbs ensemble after a quantum quench

Using a numerical renormalization group based on exploiting an underlying exactly solvable nonrelativistic theory, we study the out-of-equilibrium dynamics of a 1D Bose gas (as described by the Lieb-Liniger model) released from a parabolic trap. Our method allows us to track the postquench dynamics of the gas all the way to infinite time. We also exhibit a general construction, applicable to all integrable models, of the thermodynamic ensemble that has been suggested to govern this dynamics, the generalized Gibbs ensemble. We compare the predictions of equilibration from this ensemble against the long time dynamics observed using our method.

Quantum entanglement in random physical states

Most states in the Hilbert space are maximally entangled. This fact has proven useful to investigate--among other things--the foundations of statistical mechanics. Unfortunately, most states in the Hilbert space of a quantum many-body system are not physically accessible. We define physical ensembles of states acting on random factorized states by a circuit of length k of random and independent unitaries with local support. We study the typicality of entanglement by means of the purity of the reduced state. We find that for a time k=O(1), the typical purity obeys the area law. Thus, the upper bounds for area law are actually saturated, on average, with a variance that goes to zero for large systems. Similarly, we prove that by means of local evolution a subsystem of linear dimensions L is typically entangled with a volume law when the time scales with the size of the subsystem. Moreover, we show that for large values of k the reduced state becomes very close to the completely mixed state.

Biorthonormal transfer-matrix renormalization-group method for non-Hermitian matrices

A biorthonormal transfer-matrix renormalization-group (BTMRG) method for non-Hermitian matrices is presented. This BTMRG produces a dual set of biorthonormal bases to construct the renormalized transfer matrix with only half the dimensions of the matrix of a conventional transfer-matrix renormalization group (TMRG). We show that under generic conditions, such biorthonormal bases always exist. Based on a special E·S·E scheme (where S and E represent the system and environment blocks, respectively, and the two dots in between represent two additional physical sites), the BTMRG method can achieve zero truncation of any reduced state in describing both current left and right Perron states so as to reach a high degree of efficiency and accuracy. We believe that the BTMRG constitutes a more powerful and robust tool than conventional TMRG for non-Hermitian matrices and that it would allow us to better understand the collective behaviors and emerging phenomena of strongly correlated many-body systems. We also show that this scheme is particularly adapted to the calculation of the two-site correlation function of a one-dimensional quantum or two-dimensional classical lattice model.

Absence of thermalization in nonintegrable systems

We establish a link between unitary relaxation dynamics after a quench in closed many-body systems and the entanglement in the energy eigenbasis. We find that even if reduced states equilibrate, they can have memory on the initial conditions even in certain models that are far from integrable. We show that in such situations the equilibrium states are still described by a maximum entropy or generalized Gibbs ensemble, regardless of whether a model is integrable or not, thereby contributing to a recent debate. In addition, we discuss individual aspects of the thermalization process, comment on the role of Anderson localization, and collect and compare different notions of integrability.

Effect of rare fluctuations on the thermalization of isolated quantum systems

We consider the question of thermalization for isolated quantum systems after a sudden parameter change, a so-called quantum quench. In particular, we investigate the prerequisites for thermalization, focusing on the statistical properties of the time-averaged density matrix and of the expectation values of observables in the final eigenstates. We find that eigenstates, which are rare compared to the typical ones sampled by the microcanonical distribution, are responsible for the absence of thermalization of some infinite integrable models and play an important role for some nonintegrable systems of finite size, such as the Bose-Hubbard model. We stress the importance of finite size effects for the thermalization of isolated quantum systems and discuss two scenarios for thermalization.

Equivalence condition for the canonical and microcanonical ensembles in coupled spin systems

It is typically assumed, without justification, that a weak coupling between a system and a bath is a necessary condition for the equivalence of a canonical ensemble and a microcanonical ensemble. For instance, in a canonical ensemble, temperature emerges if the system and the bath are uncoupled or weakly coupled. We investigate the validity region of this weak-coupling approximation, using a coupled composite-spin system. Our results show that the spin coupling strength can be as large as the level spacing of the system, indicating that the weak-coupling approximation has a much wider region of validity than usually expected.

Thermalization in a coherently driven ensemble of two-level systems

We investigate the coherent quantum time evolution of a driven mesoscopic chain of two-level systems that interact via the van der Waals interaction in their excited state. The Hamiltonian is the sum of a classical lattice gas Hamiltonian and an off-diagonal driving term without classical counterpart. Starting from a product state we observe-beyond a certain interaction strength-thermalization of the system with respect to observables of the classical lattice gas. This transition can be studied experimentally with Rydberg atoms, ions, or polar molecules. We suggest how to experimentally determine the temperature of the thermal state which should allow for thermometry of the internal degrees of freedom of cold Rydberg gases whose external dynamics is frozen.

Breakdown of thermalization in finite one-dimensional systems

We use quantum quenches to study the dynamics and thermalization of hard core bosons in finite one-dimensional lattices. We perform exact diagonalizations and find that, far away from integrability, few-body observables thermalize. We then study the breakdown of thermalization as one approaches an integrable point. This is found to be a smooth process in which the predictions of standard statistical mechanics continuously worsen as the system moves toward integrability. We establish a direct connection between the presence or absence of thermalization and the validity or failure of the eigenstate thermalization hypothesis, respectively.

Supersonic quantum communication

When locally exciting a quantum lattice model, the excitation will propagate through the lattice. This effect is responsible for a wealth of nonequilibrium phenomena, and has been exploited to transmit quantum information. It is a commonly expressed belief that for local Hamiltonians, any such propagation happens at a finite "speed of sound". Indeed, the Lieb-Robinson theorem states that in spin models, all effects caused by a perturbation are essentially limited to a causal cone. We show that for meaningful translationally invariant bosonic models with nearest-neighbor interactions (addressing the challenging aspect of an experimental realization) this belief is incorrect: We prove that one can encounter accelerating excitations under the natural dynamics that allow for reliable transmission of information faster than any finite speed of sound. It also implies that the simulation of dynamics of strongly correlated bosonic models may be much harder than that of spin chains even in the low-energy sector.

Time-dependent current-density functional theory for generalized open quantum systems

In this article, we prove the one-to-one correspondence between vector potentials and particle and current densities in the context of master equations with arbitrary memory kernels, therefore extending time-dependent current-density functional theory (TD-CDFT) to the domain of generalized many-body open quantum systems (OQS). We also analyse the issue of A-representability for the Kohn-Sham (KS) scheme proposed by D'Agosta and Di Ventra for Markovian OQS [Phys. Rev. Lett. 2007, 98, 226403] and discuss its domain of validity. We suggest ways to expand their scheme, but also propose a novel KS scheme where the auxiliary system is both closed and non-interacting. This scheme is tested numerically with a model system, and several considerations for the future development of functionals are indicated. Our results formalize the possibility of practising TD-CDFT in OQS, hence expanding the applicability of the theory to non-Hamiltonian evolutions.

Relaxation of antiferromagnetic order in spin-1/2 chains following a quantum quench

We study the unitary time evolution of antiferromagnetic order in anisotropic Heisenberg chains that are initially prepared in a pure quantum state far from equilibrium. Our analysis indicates that the antiferromagnetic order imprinted in the initial state vanishes exponentially. Depending on the anisotropy parameter, oscillatory or nonoscillatory relaxation dynamics is observed. Furthermore, the corresponding relaxation time exhibits a minimum at the critical point, in contrast to the usual notion of critical slowing down, from which a maximum is expected.

Dynamical simulations of charged soliton transport in conjugated polymers with the inclusion of electron-electron interactions

We present numerical studies of the transport dynamics of a charged soliton in conjugated polymers under the influence of an external time-dependent electric field. All relevant electron-phonon and electron-electron interactions are nearly fully taken into account by simulating the monomer displacements with classical molecular dynamics and evolving the wave function for the pi electrons by virtue of the adaptive time-dependent density matrix renormalization group simultaneously and nonadiabatically. It is found that after a smooth turn on of the external electric field the charged soliton is accelerated at first up to a stationary constant velocity as one entity consisting of both the charge and the lattice deformation. An Ohmic region (6 mV/A</=E(0)</=12 mV/A) where the stationary velocity increases linearly with the electric field strength is observed. The relationship between electron-electron interactions and charged soliton transport is also investigated in detail. We find that the dependence of the stationary velocity of a charged soliton on the on-site Coulomb interactions U and the nearest-neighbor interactions V is due to the extent of delocalization of the charged soliton defect.