Interacting lattice Hamiltonians at high temperature generically give rise to energy transport governed by the classical diffusion equation; however, predicting the rate of diffusion requires numerical simulation of the microscopic quantum dynamics. For the purpose of predicting such transport properties, computational time evolution methods must be paired with schemes to control the growth of entanglement to tractably simulate for sufficiently long times. One such truncation scheme---dissipation-assisted operator evolution (DAOE)---controls entanglement by damping out components of operators with large Pauli weight. In this paper, we generalize DAOE to treat fermionic systems. Our method instead damps out components of operators with large fermionic weight. We investigate the performance of DAOE, the new fermionic DAOE (FDAOE), and another simulation method, density matrix truncation (DMT), in simulating energy transport in an interacting one-dimensional Majorana chain. The chain is found to have a diffusion coefficient scaling like interaction strength to the fourth power, contrary to naive expectations based on Fermi's Golden rule---but consistent with recent predictions based on the theory of weak integrability breaking. In the weak interaction regime where the fermionic nature of the system is most relevant, FDAOE is found to simulate the system more efficiently than DAOE.
Interacting lattice Hamiltonians at high temperature generically give rise to energy transport governed by the classical diffusion equation; however, predicting the rate of diffusion requires numerical simulation of the microscopic quantum dynamics. For the purpose of predicting such transport properties, computational time evolution methods must be paired with schemes to control the growth of entanglement to tractably simulate for sufficiently long times. One such truncation scheme---dissipation-assisted operator evolution (DAOE)---controls entanglement by damping out components of operators with large Pauli weight. In this paper, we generalize DAOE to treat fermionic systems. Our method instead damps out components of operators with large fermionic weight. We investigate the performance of DAOE, the new fermionic DAOE (FDAOE), and another simulation method, density matrix truncation (DMT), in simulating energy transport in an interacting one-dimensional Majorana chain. The chain is found to have a diffusion coefficient scaling like interaction strength to the fourth power, contrary to naive expectations based on Fermi's Golden rule---but consistent with recent predictions based on the theory of weak integrability breaking. In the weak interaction regime where the fermionic nature of the system is most relevant, FDAOE is found to simulate the system more efficiently than DAOE.
Charge and energy are expected to diffuse in interacting systems of fermions at finite temperatures, even in the absence of disorder, with the interactions inducing a crossover from the coherent and ballistic streaming of quasi-particles at early times, to incoherent diffusive behavior at late times. The relevant crossover timescales and the transport coefficients are both controlled by the strength of interactions. In this work we develop a numerical method to simulate such systems at high temperatures, applicable in a wide range of interaction strengths, by adapting Dissipation-assisted Operator Evolution (DAOE) to fermions. Our fermion DAOE, which approximates the exact dynamics by systematically discarding information from high n-point functions, is tailored to capture non-interacting dynamics exactly, thus providing a good starting point for the weakly interacting problem. Applying our method to a microscopic model of weakly interacting fermions, we numerically demonstrate that the crossover from ballistic to diffusive transport happens at a time tD∼1/Δ2 and that the diffusion constant similarly scales as D∼1/Δ2, where Δ is the interaction strength. We substantiate this scaling with a Fermi's golden rule calculation in the operator spreading picture, interpreting tD as the fermion-fermion scattering time and lifetime of the single-particle Green's function.
We introduce the dissipation-assisted operator evolution (DAOE) method for calculating transport properties of strongly interacting lattice systems in the high temperature regime. DAOE is based on evolving observables in the Heisenberg picture, and applying an artificial dissipation that reduces the weight on non-local operators. We represent the observable as a matrix product operator, and show that the dissipation leads to a decay of operator entanglement, allowing us to capture the dynamics to long times. We test this scheme by calculating spin and energy diffusion constants in a variety of physical models. By gradually weakening the dissipation, we are able to consistently extrapolate our results to the case of zero dissipation, thus estimating the physical diffusion constant with high precision.
In the absence of disorder and interactions, fermions move coherently and their associated charge and energy exhibit ballistic spreading, even at finite energy density. In the presence of weak interactions and a finite energy density, fermion-fermion scattering leads to a crossover between early-time ballistic and late-time diffusive transport. The relevant crossover timescales and the transport coefficients are both functions of interaction strength, but the question of determining the precise functional dependence is likely impossible to answer exactly. In this work we develop a numerical method (fDAOE) which is powerful enough to provide an approximate answer to this question, and which is consistent with perturbative arguments in the limit of very weak interactions. Our algorithm, which adapts the existing dissipation-assisted operator evolution (DAOE) to fermions, is applicable to systems of interacting fermions at high temperatures. The algorithm approximates the exact dynamics by systematically discarding information from high n -point functions, and is tailored to capture noninteracting dynamics exactly. Applying our method to a microscopic model of interacting fermions, we numerically determine crossover timescales and diffusion constants for a wide range of interaction strengths. In the limit of weak interaction strength ( Δ ), we demonstrate that the crossover from ballistic to diffusive transport happens at a time tD∼1/Δ2 and that the diffusion constant similarly scales as D∼1/Δ2 . We confirm that these scalings are consistent with a perturbative Fermi's golden rule calculation, and we provide a heuristic operator-spreading picture for the crossover between ballistic and diffusive transport. Published by the American Physical Society 2024
We consider the task of simulating the real-time dynamics of weakly interacting fermionic systems. In particular, we focus on computing the expectation value of a local observable A at time t. By analyzing the convergence of the perturbative expansion in the interaction strength λ for the Heisenberg-picture observable, we propose a polynomial-time algorithm for estimating this expectation value in the weakly interacting regime λ∣t∣2D+1=O(1), when the Hamiltonian is geometrically local on a D-dimensional lattice. Importantly, this condition is independent of the system size. If the goal is instead to approximate the time-evolved observable in normalized Frobenius norm, we extend the convergence regime to λ∣t∣=O(1) with quasi-polynomial runtime. When the non-interacting part exhibits Anderson localization, our polynomial-time algorithm can be extended up to λ∣t∣=O(1), modulo polylogarithmic factors. Our algorithm brings together ideas from continuous-time QMC, diagrammatic QMC, and Majorana Propagation, but with a new Heisenberg-picture operator-growth analysis that makes the sampling complexity rigorously controllable. This leads to provably efficient classical algorithms in regimes where the interaction is weak enough that the sampling variance remains bounded independently of system size. Together, these results identify broad regimes in which weak interactions, locality, and localization can be leveraged to make real-time fermionic dynamics classically tractable.
In interacting quantum systems, the single-particle Green's function is expected to decay in time due to the interaction induced decoherence of quasiparticles. In the limit of weak interaction strengths (Δ), a naive application of Fermi's Golden Rule (FGR) predicts an O(Δ2) quasiparticle decay rate. However, for 1d fermions on the lattice at T>0, this calculation gives a divergent result and the scaling of the quasiparticle lifetime with interaction strength remains an open question. In this work we propose a solution to this question: combining numerical simulations using the recently introduced dissipation-assisted operator evolution (DAOE) method, with non-perturbative diagrammatic re-summations, we predict a logarithmic enhancement of the quasiparticle decay rate τ−1∼Δ2logΔ−2. We argue that this effect is present in a wide variety of well-known weakly interacting quantum fermionic and bosonic systems, and even in some classical systems, provided the non-interacting limit has quasiparticles with a generic dispersion.
We use nonequilibrium steady states to study the effect of dissipation-assisted operator evolution (DAOE) on the scaling behavior of transport in one-dimensional spin chains. We consider three models in the XXZ family: the XXZ model with staggered anisotropy, which is chaotic; XXZ model with no external field and tunable interaction, which is Bethe-ansatz integrable and (in the zero interaction limit) free-fermion integrable; and the disordered XY model, which is free-fermion integrable and Anderson localized. We find evidence that DAOE's effect on transport is controlled by its effect on the system's conserved quantities. To the extent that DAOE preserves those symmetries, it preserves the scaling of the system's transport properties; to the extent it breaks those conserved quantities, it pushes the system towards diffusive scaling of transport.