Dynamics Of Quantum Many-body Systems Research Articles

Information scrambling -- a quantum thermodynamic perspective

Abstract Recent advances in quantum information science have shed light on the intricate dynamics of quantum many-body systems, for which quantum information scrambling is a perfect example. Motivated by considerations of the thermodynamics of quantum information, this perspective aims at synthesizing key findings from several pivotal studies, exploring various aspects of quantum scrambling. We consider quantifiers such as the Out-of-Time-Ordered Correlator (OTOC), the quantum Mutual Information, and the Tripartite Mutual Information (TMI), their connections to thermodynamics, and their role in understanding chaotic versus integrable quantum systems. Focusing on representative examples, we cover a range of topics, including the thermodynamics of quantum information scrambling, and the scrambling dynamics in quantum gravity models such as the Sachdev-Ye-Kitaev (SYK) model. By examining these diverse approaches, we highlight the multifaceted nature of quantum information scrambling and its significance in understanding the fundamental aspects of quantum many-body dynamics at the intersection of quantum mechanics and thermodynamics.

Full counting statistics as probe of measurement-induced transitions in the quantum Ising chain

Non-equilibrium dynamics of many-body quantum systems under the effect of measurement protocols is attracting an increasing amount of attention. It has been recently revealed that measurements may induce different non-equilibrium regimes and an abrupt change in the scaling-law of the bipartite entanglement entropy. However, our understanding of how these regimes appear{, how they affect the statistics of local quantities and}, finally whether they survive in the thermodynamic limit, is much less established. Here we investigate measurement-induced phase transitions in the Quantum Ising chain coupled to a monitoring environment. In particular we show that local projective measurements induce a quantitative modification of the out-of-equilibrium probability distribution function of the local magnetization. Starting from a GHZ state, the relaxation of the paramagnetic and the ferromagnetic order is analysed. In particular we describe how the probability distributions associated to them show different behaviour depending on the measurement rate.

The singularities of the rate function of quantum coherent work in one-dimensional transverse field Ising model

Quantum coherence will undoubtedly play a fundamental role in understanding of the dynamics of quantum many-body systems, thereby to reveal its genuine contribution is of great importance. In this paper, we specialize our discussions to the one-dimensional transverse field quantum Ising model initialized in the coherent Gibbs state. After quenching the strength of the transverse field, the effects of quantum coherence are studied by the rate function of quantum work distribution. We find that quantum coherence not only recovers the quantum phase transition destroyed by thermal fluctuations, but also generates some entirely new singularities both in the static state and dynamics. It can be manifested that these singularities are rooted in spin flips causing the sudden change of the domain boundaries of spin polarization. This work sheds new light on the fundamental connection between quantum critical phenomena and quantum coherence.

Entanglement Barrier and its Symmetry Resolution: Theory and Experimental Observation

An efficient method to extract higher-order density-matrix functionals and access operator entanglement is introduced and tested on experimental data, improving the characterization of nonequilibrium dynamics of many-body quantum systems.

Phase-space stochastic quantum hydrodynamics for interacting Bose gases

Hydrodynamic theories offer successful approaches that are capable of simulating the otherwise difficult-to-compute dynamics of quantum many-body systems. In this work we derive, within the positive-P phase-space formalism, a new stochastic hydrodynamic method for the description of interacting Bose gases. It goes beyond existing hydrodynamic approaches, such as superfluid hydrodynamics or generalized hydrodynamics, in its capacity to simulate the full quantum dynamics of these systems: it possesses the ability to compute non-equilibrium quantum correlations, even for short-wavelength phenomena. Using this description, we derive a linearized stochastic hydrodynamic scheme which is able to simulate such non-equilibrium situations for longer times than the full positive-P approach, at the expense of approximating the treatment of quantum fluctuations, and show that this linearized scheme can be directly connected with existing Bogoliubov approaches. Furthermore, we go on to demonstrate the usefulness and advantages of this formalism by exploring the correlations that arise in a quantum shock wave scenario and comparing its predictions to other established quantum many-body approaches.

Probing Transport and Slow Relaxation in the Mass-Imbalanced Fermi-Hubbard Model

Constraints in the dynamics of quantum many-body systems can dramatically alter transport properties and relaxation timescales even in the absence of static disorder. Here, we report on the observation of such constrained dynamics arising from the distinct mobility of two species in the one-dimensional mass-imbalanced Fermi-Hubbard model, realized with ultracold ytterbium atoms in a state-dependent optical lattice. By displacing the trap potential and monitoring the subsequent dynamical response of the system, we identify suppressed transport and slow relaxation with a strong dependence on the mass imbalance and interspecies interaction strength, consistent with eventual thermalization for long times. Our observations demonstrate the potential for quantum simulators to provide insights into unconventional relaxation dynamics arising from constraints.

Tunable transport in the mass-imbalanced Fermi-Hubbard model

The late-time dynamics of quantum many-body systems is organized in distinct dynamical universality classes, characterized by their conservation laws and thus by their emergent hydrodynamic transport. Here, we study transport in the one-dimensional Hubbard model with different masses of the two fermionic species. To this end, we develop a quantum Boltzmann approach valid in the limit of weak interactions. We explore the crossover from ballistic to diffusive transport, whose timescale strongly depends on the mass ratio of the two species. For timescales accessible with matrix product operators, we find excellent agreement between these numerically exact results and the quantum Boltzmann equation, even for intermediate interactions. We investigate two scenarios which have been recently studied with ultracold atom experiments. First, in the presence of a tilt, the quantum Boltzmann equation predicts that transport is significantly slowed down and becomes subdiffusive, consistent with previous studies. Second, we study transport probed by displacing a harmonic confinement potential and find good quantitative agreement with recent experimental data [N. Darkwah Oppong et al., arXiv:2011.12411]. Our results demonstrate that the quantum Boltzmann equation is a useful tool to study complex non-equilibrium states in inhomogeneous potentials, as often probed with synthetic quantum systems.

Unified entropies and quantum speed limits for nonunitary dynamics

We discuss a class of quantum speed limits (QSLs) based on unified quantum ($\alpha,\mu$)-entropy for nonunitary physical processes. The bounds depend on both the Schatten speed and the smallest eigenvalue of the evolved state, and the two-parametric unified entropy. We specialize these results to quantum channels and non-Hermitian evolutions. In the first case, the QSL depends on the Kraus operators related to the quantum channel, while in the second case the speed limit is recast in terms of the non-Hermitian Hamiltonian. To illustrate these findings, we consider a single-qubit state evolving under the (i) amplitude damping channel, and (ii) the nonunitary dynamics generated by a parity-time-reversal symmetric non-Hermitian Hamiltonian. The QSL is nonzero at earlier times, while it becomes loose as the smallest eigenvalue of the evolved state approaches zero. Furthermore, we investigate the interplay between unified entropies and speed limits for the reduced dynamics of quantum many-body systems. The unified entropy is upper bounded by the quantum fluctuations of the Hamiltonian of the system, while the QSL is nonzero when entanglement is created by the nonunitary evolution. Finally, we apply these results to the XXZ model and verify that the QSL asymptotically decreases as the system size increases. Our results find applications to nonequilibrium thermodynamics, quantum metrology, and equilibration of many-body systems.

Tunable transport in the mass-imbalanced Fermi Hubbard model

The late-time dynamics of quantum many-body systems is organized in distinct dynamical universality classes, characterized by their conservation laws and thus by their emergent hydrodynamic transport. Here, we study transport in the one-dimensional Hubbard model with different masses of the two fermionic species. To this end, we develop a quantum Boltzmann approach valid in the limit of weak interactions. We explore the crossover from ballistic to diffusive transport, whose timescale strongly depends on the mass ratio of the two species. For timescales accessible with matrix product operators, we find excellent agreement between these numerically exact results and the quantum Boltzmann equation, even for intermediate interactions. We investigate two scenarios which have been recently studied with ultracold atom experiments. First, in the presence of a tilt, the quantum Boltzmann equation predicts that transport is significantly slowed down and becomes subdiffusive, consistent with previous studies. Second, we study transport probed by displacing a harmonic confinement potential and find good quantitative agreement with recent experimental data from [Oppong et al., arXiv:2011.12411]. Our results demonstrate that the quantum Boltzmann equation is a useful tool to study complex non-equilibrium states in inhomogeneous potentials, as often probed with synthetic quantum systems.

Symmetry Classification and Universality in Non-Hermitian Many-Body Quantum Chaos by the Sachdev-Ye-Kitaev Model

Spectral correlations are a powerful tool to study the dynamics of quantum many-body systems. For Hermitian Hamiltonians, quantum chaotic motion is related to random matrix theory spectral correlations. Based on recent progress in the application of spectral analysis to non-Hermitian quantum systems, we show that local level statistics, which probe the dynamics around the Heisenberg time, of a non-Hermitian $q$-body Sachdev-Ye-Kitev (nHSYK) model with $N$ Majorana fermions, and its chiral and complex-fermion extensions, are also well described by random matrix theory for $q > 2$, while for $q = 2$, they are given by the equivalent of Poisson statistics. For that comparison, we combine exact diagonalization numerical techniques with analytical results obtained for some of the random matrix spectral observables. Moreover, depending on $q$ and $N$, we identify $19$ out of the $38$ non-Hermitian universality classes in the nHSYK model, including those corresponding to the tenfold way. In particular, we realize explicitly $14$ out of the $15$ universality classes corresponding to non-pseudo-Hermitian Hamiltonians that involve universal bulk correlations of classes ${\rm AI}^\dagger$ and ${\rm AII}^\dagger$, beyond the Ginibre ensembles. These results provide strong evidence of striking universal features in non-unitary many-body quantum chaos, which in all cases can be captured by nHSYK models with $q > 2$.

Entanglement Dynamics From Random Product States: Deviation From Maximal Entanglement

We study the entanglement dynamics of quantum many-body systems and prove the following: (I) For any geometrically local Hamiltonian on a lattice, starting from a random product state the entanglement entropy is bounded away from the maximum entropy at all times with high probability. (II) In a spin-glass model with random all-to-all interactions, starting from any product state the average entanglement entropy is bounded away from the maximum entropy at all times. We also extend these results to any unitary evolution with charge conservation and to the Sachdev-Ye-Kitaev model. Our results highlight the difference between the entanglement generated by (chaotic) Hamiltonian dynamics and that of random states, for the latter is nearly maximal.

Noise-Induced Entanglement Transition in One-Dimensional Random Quantum Circuits

A random quantum circuit is a minimally structured model to study entanglement dynamics of many-body quantum systems. We consider a one-dimensional quantum circuit with noisy Haar-random unitary gates using density matrix operator and tensor contraction methods. It is shown that the entanglement evolution of the random quantum circuits is properly characterized by the logarithmic entanglement negativity. By performing exact numerical calculations, we find that, as the physical error rate is decreased below a critical value p c ≈ 0.056, the logarithmic entanglement negativity changes from the area law to the volume law, giving rise to an entanglement transition. The critical exponent of the correlation length can be determined from the finite-size scaling analysis, revealing the universal dynamic property of the noisy intermediate-scale quantum devices.

NISQ Algorithm for Hamiltonian simulation via truncated Taylor series

Simulating the dynamics of many-body quantum systems is believed to be one of the first fields that quantum computers can show a quantum advantage over classical computers. Noisy intermediate-scale quantum (NISQ) algorithms aim at effectively using the currently available quantum hardware. For quantum simulation, various types of NISQ algorithms have been proposed with individual advantages as well as challenges. In this work, we propose a new algorithm, truncated Taylor quantum simulator (TQS), that shares the advantages of existing algorithms and alleviates some of the shortcomings. Our algorithm does not have any classical-quantum feedback loop and bypasses the barren plateau problem by construction. The classical part in our hybrid quantum-classical algorithm corresponds to a quadratically constrained quadratic program (QCQP) with a single quadratic equality constraint, which admits a semidefinite relaxation. The QCQP based classical optimization was recently introduced as the classical step in quantum assisted eigensolver (QAE), a NISQ algorithm for the Hamiltonian ground state problem. Thus, our work provides a conceptual unification between the NISQ algorithms for the Hamiltonian ground state problem and the Hamiltonian simulation. We recover differential equation-based NISQ algorithms for Hamiltonian simulation such as quantum assisted simulator (QAS) and variational quantum simulator (VQS) as particular cases of our algorithm. We test our algorithm on some toy examples on current cloud quantum computers. We also provide a systematic approach to improve the accuracy of our algorithm.

Eigenstate fluctuation theorem in the short- and long-time regimes.

The canonical ensemble plays a crucial role in statistical mechanics in and out of equilibrium. For example, the standard derivation of the fluctuation theorem relies on the assumption that the initial state of the heat bath is the canonical ensemble. On the other hand, the recent progress in the foundation of statistical mechanics has revealed that a thermal equilibrium state is not necessarily described by the canonical ensemble but can be a quantum pure state or even a single energy eigenstate, as formulated by the eigenstate thermalization hypothesis (ETH). Then a question raised is how these two pictures, the canonical ensemble and a single energy eigenstate as a thermal equilibrium state, are compatible in the fluctuation theorem. In this paper, we theoretically and numerically show that the fluctuation theorem holds in both of the long- and short-time regimes, even when the initial state of the bath is a single energy eigenstate of a many-body system. Our proof of the fluctuation theorem in the long-time regime is based on the ETH, while it was previously shown in the short-time regime on the basis of the Lieb-Robinson bound and the ETH [Phys. Rev. Lett. 119, 100601 (2017)0031-900710.1103/PhysRevLett.119.100601]. The proofs for these time regimes are theoretically independent and complementary, implying the fluctuation theorem in the entire time domain. We also perform a systematic numerical simulation of hard-core bosons by exact diagonalization and verify the fluctuation theorem in both of the time regimes by focusing on the finite-size scaling. Our results contribute to the understanding of the mechanism that the fluctuation theorem emerges from unitary dynamics of quantum many-body systems and can be tested by experiments with, e.g., ultracold atoms.

Growth of entanglement entropy under local projective measurements

Non-equilibrium dynamics of many-body quantum systems under the effect of measurement protocols is attracting an increasing amount of attention. It has been recently revealed that measurements may induce an abrupt change in the scaling-law of the bipartite entanglement entropy, thus suggesting the existence of different non-equilibrium regimes. However, our understanding of how these regimes appear and whether they survive in the thermodynamic limit is much less established. Here we investigate these questions on a one-dimensional quadratic fermionic model: this allows us to reach system sizes relevant in the thermodynamic sense. We show that local projective measurements induce a qualitative modification of the time-growth of the entanglement entropy which changes from linear to logarithmic. However, in the stationary regime, the logarithmic behavior of the entanglement entropy do not survive in the thermodynamic limit and, for any finite value of the measurement rate, we numerically show the existence of a single area-law phase for the entanglement entropy. Finally, exploiting the quasi-particle picture, we further support our results analysing the fluctuations of the stationary entanglement entropy and its scaling behavior.

Real time evolution with neural-network quantum states

A promising application of neural-network quantum states is to describe the time dynamics of many-body quantum systems. To realize this idea, we employ neural-network quantum states to approximate the implicit midpoint rule method, which preserves the symplectic form of Hamiltonian dynamics. We ensure that our complex-valued neural networks are holomorphic functions, and exploit this property to efficiently compute gradients. Application to the transverse-field Ising model on a one- and two-dimensional lattice exhibits an accuracy comparable to the stochastic configuration method proposed in [Carleo and Troyer, Science 355, 602-606 (2017)], but does not require computing the (pseudo-)inverse of a matrix.

Non-equilibrium dynamics of the open quantum O(n)-model with non-Markovian noise: exact results

The collective and purely relaxational dynamics of quantum many-body systems after a quench at temperature T = 0, from a disordered state to various phases is studied through the exact solution of the quantum Langevin equation of the spherical and the O(n)-model in the limit n → ∞. The stationary state of the quantum dynamics is shown to be a non-equilibrium state. The quantum spherical and the quantum O(n)-model for n → ∞ are in the same dynamical universality class. The long-time behaviour of single-time and two-time correlation and response functions is analysed and the universal exponents which characterise quantum coarsening and quantum ageing are derived. The importance of the non-Markovian long-time memory of the quantum noise is elucidated by comparing it with an effective Markovian noise having the same scaling behaviour and with the case of non-equilibrium classical dynamics.

Scaling of temporal entanglement in proximity to integrability

Describing dynamics of quantum many-body systems is a formidable challenge due to rapid generation of quantum entanglement between remote degrees of freedom. A promising approach to tackle this challenge, which has been proposed recently, is to characterize the quantum dynamics of a many-body system and its properties as a bath via the Feynman-Vernon influence matrix (IM), which is an operator in the space of time trajectories of local degrees of freedom. Physical understanding of the general scaling of the IM's temporal entanglement and its relation to basic dynamical properties is highly incomplete to present day. In this Article, we analytically compute the exact IM for a family of integrable Floquet models - the transverse-field kicked Ising chain - finding a Bardeen-Cooper-Schrieffer-like "wavefunction" on the Schwinger-Keldysh contour with algebraically decaying correlations. We demonstrate that the IM exhibits area-law temporal entanglement scaling for all parameter values. Furthermore, the entanglement pattern of the IM reveals the system's phase diagram, exhibiting jumps across transitions between distinct Floquet phases. Near criticality, a non-trivial scaling behavior of temporal entanglement is found. The area-law temporal entanglement allows us to efficiently describe the effects of sizeable integrability-breaking perturbations for long evolution times by using matrix product state methods. This work shows that tensor network methods are efficient in describing the effect of non-interacting baths on open quantum systems, and provides a new approach to studying quantum many-body systems with weakly broken integrability.

Quantum aging and dynamical universality in the long-range O(N→∞) model.

Quantum quenches to or near criticality give rise to the phenomenon of aging, manifested by glassylike dynamics at short times and far from equilibrium. The recent surge of interest in the dynamics of quantum many-body systems has rejuvenated interest in this phenomenon. Motivated by the ubiquitous long-range interactions in emerging experimental platforms, it is vital to study quantum aging in such settings. In this paper, we investigate the dynamical universality and aging in the d-dimensional O(N) model with the long-range coupling 1/x^{d+σ} and in the mean-field limit N→∞ that allows an exact treatment. An immediate consequence of long-range coupling is the emergence of nonlinear light cones. We focus on the correlation and response functions, and identify a rich scaling behavior depending on how the corresponding space-time positions are located relative to each other, via a local light cone, and to the time of the quench via a global quench light cone. We determine the initial-slip exponent that governs the short-time dependence of two-point functions. We highlight the qualitative features of aging due to the long-range coupling, in particular in the region outside the light cones. As an important consequence of long-range coupling, the correlation function decays as 1/x^{d+σ} outside the quench light cone while increasing polynomially with the total time after quench. This is while, for short-time differences, the two-time response function "equilibrates" at all distances even outside this light cone. Our analytic findings are in excellent agreement with exact numerics, and provide a useful benchmark for modern experimental platforms with long-range interactions.

Variational classical networks for dynamics in interacting quantum matter

Dynamics in correlated quantum matter is a hard problem, as its exact solution generally involves a computational effort that grows exponentially with the number of constituents. While a remarkable progress has been witnessed in recent years for one-dimensional systems, much less has been achieved for interacting quantum models in higher dimensions, since they incorporate an additional layer of complexity. In this work, we employ a variational method that allows for an efficient and controlled computation of the dynamics of quantum many-body systems in one and higher dimensions. The approach presented here introduces a variational class of wavefunctions based on complex networks of classical spins akin to artificial neural networks, which can be constructed in a controlled fashion. We provide a detailed prescription for such constructions and illustrate their performance by studying quantum quenches in one- and two-dimensional models. In particular, we investigate the nonequilibrium dynamics of a genuinely interacting two-dimensional lattice gauge theory, the quantum link model, for which we have recently shown -- employing the technique discussed thoroughly in this paper -- that it features disorder-free localization dynamics [P. Karpov et al., Phys. Rev. Lett. 126, 130401 (2021)]. The present work not only supplies a framework to address purely theoretical questions but also could be used to provide a theoretical description of experiments in quantum simulators, which have recently seen an increased effort targeting two-dimensional geometries. Importantly, our method can be applied to any quantum many-body system with a well-defined classical limit.