Describing a quantum impurity coupled to one or more noninteracting fermionic reservoirs is a paradigmatic problem in quantum many-body physics. While historically the focus has been on the equilibrium properties of the impurity-reservoir system, recent experiments with mesoscopic and cold-atomic systems enabled studies of highly nonequilibrium impurity models, which require novel theoretical techniques. We propose an approach to analyze impurity dynamics based on the matrix-product state (MPS) representation of the Feynman-Vernon influence functional (IF). The efficiency of such a MPS representation rests on the moderate value of the temporal entanglement (TE) entropy of the IF, viewed as a fictitious ``wave function'' in the time domain. We obtain explicit expressions of this wave function for a family of one-dimensional reservoirs, and analyze the scaling of TE with the evolution time for different reservoir's initial states. While for initial states with short-range correlations we find temporal area-law scaling, Fermi-sea-type initial states yield logarithmic scaling with time, closely related to the real-space entanglement scaling in critical $1d$ systems. Furthermore, we describe an efficient algorithm for converting the explicit form of general reservoirs' IF to MPS form. Once the IF is encoded by a MPS, arbitrary temporal correlation functions of the interacting impurity can be efficiently computed, irrespective of the form of impurity interactions and bath geometry. The approach introduced here can be applied to a number of experimental setups, including highly nonequilibrium transport via quantum dots and real-time formation of impurity-reservoir correlations, as well as in nonequilibrium dynamical mean-field theory.
Describing a quantum impurity coupled to one or more noninteracting fermionic reservoirs is a paradigmatic problem in quantum many-body physics. While historically the focus has been on the equilibrium properties of the impurity-reservoir system, recent experiments with mesoscopic and cold-atomic systems enabled studies of highly nonequilibrium impurity models, which require novel theoretical techniques. We propose an approach to analyze impurity dynamics based on the matrix-product state (MPS) representation of the Feynman-Vernon influence functional (IF). The efficiency of such a MPS representation rests on the moderate value of the temporal entanglement (TE) entropy of the IF, viewed as a fictitious ``wave function'' in the time domain. We obtain explicit expressions of this wave function for a family of one-dimensional reservoirs, and analyze the scaling of TE with the evolution time for different reservoir's initial states. While for initial states with short-range correlations we find temporal area-law scaling, Fermi-sea-type initial states yield logarithmic scaling with time, closely related to the real-space entanglement scaling in critical 1d systems. Furthermore, we describe an efficient algorithm for converting the explicit form of general reservoirs' IF to MPS form. Once the IF is encoded by a MPS, arbitrary temporal correlation functions of the interacting impurity can be efficiently computed, irrespective of the form of impurity interactions and bath geometry. The approach introduced here can be applied to a number of experimental setups, including highly nonequilibrium transport via quantum dots and real-time formation of impurity-reservoir correlations, as well as in nonequilibrium dynamical mean-field theory.
We introduce an efficient method to simulate the dynamics of an interacting quantum impurity coupled to noninteracting fermionic reservoirs. Viewing the impurity as an open quantum system, we describe the reservoirs by their Feynman-Vernon influence functionals (IFs). The IFs are represented as matrix-product states in the temporal domain, which enables an efficient computation of the dynamics for arbitrary interactions. We apply our method to study quantum quenches and transport in an Anderson impurity model, including highly nonequilibrium setups, and find a favorable performance compared to state-of-the-art methods. The computational resources required for an accurate computation of the dynamics scale polynomially with evolution time, indicating that a broad class of out-of-equilibrium quantum impurity problems are efficiently solvable. This approach will provide additional insights into the dynamical properties of mesoscopic devices and correlated materials.
Feynman–Vernon influence functional (IF) was originally introduced to describe the effect of a quantum environment on the dynamics of an open quantum system. We apply the IF approach to describe quantum many-body dynamics in isolated spin systems, viewing the system as an environment for its local subsystems. While the IF can be computed exactly only in certain many-body models, it generally satisfies a self-consistency equation, provided the system, or an ensemble of systems, are translationally invariant. We view the IF as a fictitious wavefunction in the temporal domain, and approximate it using matrix-product states (MPS). This approach is efficient provided the temporal entanglement of the IF is sufficiently low. We illustrate the broad applicability of the IF approach by analyzing several models that exhibit a range of dynamical behaviors, from thermalizing to many-body localized. In particular, we study the non-equilibrium dynamics in the quantum Ising model in both Floquet and Hamiltonian settings. We find that temporal entanglement entropy may be significantly lower compared to the spatial entanglement and analyze the IF in the continuous-time limit. We simulate the thermodynamic-limit evolution of local observables in various regimes, including the relaxation of impurities embedded in an infinite-temperature chain, and the long-lived oscillatory dynamics of the magnetization associated with the confinement of excitations. Furthermore, by incorporating disorder-averaging into the formalism, we analyze discrete time-crystalline response using the IF of a bond-disordered kicked Ising chain. In this case, we find that the temporal entanglement entropy scales logarithmically with evolution time. The IF approach offers a new lens on many-body non-equilibrium phenomena, both in ergodic and non-ergodic regimes, connecting the theory of open quantum systems to quantum statistical physics.
Feynman-Vernon influence functional (IF) was originally introduced to describe the effect of a quantum environment on the dynamics of an open quantum system. We apply the IF approach to describe quantum many-body dynamics in isolated spin systems, viewing the system as an environment for its local subsystems. While the IF can be computed exactly only in certain many-body models, it generally satisfies a self-consistency equation, provided the system, or an ensemble of systems, are translationally invariant. We view the IF as a fictitious wavefunction in the temporal domain, and approximate it using matrix-product states (MPS). This approach is efficient provided the temporal entanglement of the IF is sufficiently low. We illustrate the versatility of the IF approach by analyzing several models that exhibit a range of dynamical behaviors, from thermalizing to many-body localized. In particular, we study the non-equilibrium dynamics in the quantum Ising model in both Floquet and Hamiltonian settings. We find that temporal entanglement entropy may be significantly lower compared to the spatial entanglement and analyze the IF in the continuous-time limit. We simulate the thermodynamic-limit evolution of local observables in various regimes, including the relaxation of impurities embedded in an infinite-temperature chain, and the long-lived oscillatory dynamics of the magnetization associated with the confinement of excitations. By incorporating disorder-averaging into the formalism, we analyze discrete time-crystalline response using the IF method. In this case, we find that the temporal entanglement entropy scales logarithmically with evolution time. The IF approach offers a new lens on many-body non-equilibrium phenomena, both in ergodic and non-ergodic regimes, connecting the theory of open quantum systems theory to quantum statistical physics.
The Grassmann time-evolving matrix product operator (GTEMPO) method, which represents the Feynman-Vernon influence functional as a temporal matrix product state, has been shown to be a flexible and potentially scalable solution for fermionic quantum impurity problems. In this work, we extend GTEMPO to solve fermionic impurity problems in the Nambu formalism, in which the impurity is coupled to a superconducting bath. A key insight is that by employing the Bogoliubov transformation for the superconducting bath, one could obtain the analytic expression of the Feynman-Vernon influence functional in a similar form to the case of a normal bath, after which the core algorithms of GTEMPO can be straightforwardly adapted. We demonstrate the accuracy of our method by benchmarking it against exact diagonalization in several exactly solvable cases, and against the continuous-time quantum Monte Carlo method using converged dynamical mean field theory (DMFT) iterations on the imaginary contour in the non-integrable case. In all cases, we perform both imaginary- and real-time calculations to illustrate the flexibility of our method. These results illustrate that our method could be potentially useful as an impurity solver in DMFT as well as its non-equilibrium extension for fermionic impurity problems in the Nambu formalism.
Recent experimental and theoretical works have made much progress toward understanding nonequilibrium phenomena in thermalizing systems, which act as thermal baths for their small subsystems, and many-body localized systems, which fail to do so.The description of time evolution in many-body systems is generally challenging due to the dynamical generation of quantum entanglement.In this work, we introduce an approach to study quantum many-body dynamics, inspired by the Feynman-Vernon influence functional.Focusing on a family of interacting, Floquet spin chains, we consider a Keldysh path-integral description of the dynamics.The central object in our approach is the influence matrix, which describes the effect of the system on the dynamics of a local subsystem.For translationally invariant models, we formulate a self-consistency equation for the influence matrix.For certain special values of the model parameters, we obtain an exact solution which represents a perfect dephaser (PD).Physically, a PD corresponds to a many-body system that acts as a perfectly Markovian bath on itself: at each period, it measures every spin.For the models considered here, we establish that PD points include dual-unitary circuits investigated in recent works.In the vicinity of PD points, the system is not perfectly Markovian, but rather acts as a bath with a short memory time.In this case, we demonstrate that the self-consistency equation can be solved using matrix-product states (MPS) methods, as the influence matrix temporal entanglement is low.A combination of analytical insights and MPS computations allows us to characterize the structure of the influence matrix in terms of an effective "statistical-mechanics" description.We finally illustrate the predictive power of this description by analytically computing how quickly an embedded impurity spin thermalizes.The influence matrix approach formulated here provides an intuitive view of the quantum many-body dynamics problem, opening a path to constructing models of thermalizing dynamics that are solvable or can be efficiently treated by MPS-based methods and to further characterizing quantum ergodicity or lack thereof.
Experimental advances have made it possible to realize and control quantum many-body systems, allowing the experimental study of non-equilibrium phenomena. Due to the Hilbert space growing exponentially with the number of particles, advanced numerical techniques are usually required to understand many-body effects at the theoretical level. In this thesis, we address the question of how to efficiently simulate the non-equilibrium dynamics of many-body systems using matrix product states (MPS). In particular, we consider quantum systems in interaction with macroscopic environments, as they appear in Anderson impurity problems, quantum thermodynamics, or open quantum systems. We develop a low-entanglement representation of the environment with only short-range interactions, perfectly suitable for simulations with MPS. We further show that an interleaved ordering of tight-binding chains can significantly reduce the creation of entanglement, as opposed to an intuitive implementation of the Hamiltonian's geometry. Our approach allows long-time simulations and an analysis of the environmental dynamics, in which, after a sudden quench, we find clear signatures of many-body effects. We employ our techniques to compute spectral functions of impurity models, with possible applications to dynamical mean-field theory impurity solvers, and to calculate dissipation in the non-equilibrium Anderson model with explicit time-dependence, as relevant for ongoing experiments with oscillating tip atomic force microscopes.