Numerical and analytic results have been used to characterize quantum transport in spin chains, showing the existence of both ballistic and diffusive motion. Experiments have shown that heat transfer is surprisingly always diffusive. The scattering from phonons and impurities have been postulated to be the two factors critical in causing the diffusive transport. In this work, we evaluate the transport process by incorporating a bath of phonons and impurities in order to understand the role played by each of the factors. While methods like time-dependent density matrix renormalization group (tDMRG) can be used to simulate isolated spin chains, the coupling with phonons make simulations significantly more challenging. The recently developed multisite tensor network path integral (MS-TNPI) method builds a framework for simulating the dynamics in extended open quantum systems by combining ideas from tDMRG and Feynman-Vernon influence functional. This MS-TNPI is used to characterize dynamics in open, extended quantum systems. Simulations are done with the commonly used sub-Ohmic, Ohmic and super-Ohmic spectral densities describing the phononic bath. We show that while the transport in presence of impurities eventually becomes diffusive, the exact details are dependent on the specifics of the interactions and amount of impurities. In contrast, the presence of a bath makes the transport diffusive irrespective of the parameters characterizing the bath.
Numerical and analytic results have been used to characterize quantum transport in spin chains, showing the existence of both ballistic and diffusive motion. Experiments have shown that heat transfer is surprisingly always diffusive. The scattering from phonons and impurities have been postulated to be the two factors critical in causing the diffusive transport. In this work, we evaluate the transport process by incorporating a bath of phonons and impurities in order to understand the role played by each of the factors. While methods like time-dependent density matrix renormalization group (tDMRG) can be used to simulate isolated spin chains, the coupling with phonons make simulations significantly more challenging. The recently developed multisite tensor network path integral (MS-TNPI) method builds a framework for simulating the dynamics in extended open quantum systems by combining ideas from tDMRG and Feynman-Vernon influence functional. This MS-TNPI is used to characterize dynamics in open, extended quantum systems. Simulations are done with the commonly used sub-Ohmic, Ohmic and super-Ohmic spectral densities describing the phononic bath. We show that while the transport in presence of impurities eventually becomes diffusive, the exact details are dependent on the specifics of the interactions and amount of impurities. In contrast, the presence of a bath makes the transport diffusive irrespective of the parameters characterizing the bath.
Understanding the transport properties of quantum many-body systems is a central challenge in condensed matter and statistical physics. Theoretical studies usually rely on two main approaches: Dynamics of linear-response functions in closed systems and boundary-driven dynamics governed by Markovian master equations for open systems. While the equivalence of their dynamical behavior has been explored in recent studies, a systematic comparison of the transport coefficients obtained from these two classes of approaches remains a largely open question. Here, we address this gap by comparing and contrasting the dc diffusion constant Ddc according to the two approaches, focusing on the specific example of magnetization transport in the spin-1/2 XXZ chain. Using exact numerical simulations for finite system sizes, we find (i) a clear mismatch between the two Ddc and (ii) a strong dependence of Ddc on the system-bath coupling strength for the open system, where neither (i) nor (ii) tend to vanish in the thermodynamic limit. These findings suggest limitations of the open-system approach to transport coefficients. To gain insight into the origin of (i) and (ii), we go beyond Ddc and analyze the full time dependence of the diffusion coefficient D(t) in the open system. In this way, we find that both (i) and (ii) vanish up to a finite time scale. While this time scale gradually increases with system size and tends to be macroscopic in the thermodynamic limit, this increase is still slow compared to the increase of the time to reach the steady state, where (i) and (ii) do not vanish. This observation can be seen as a wrong, yet unavoidable order of limits of long times first and large system sizes afterwards.
Based on the references [P1–P6], we investigate in this cumulative PhD thesis the dynamics and transport in closed and open interacting spin chain models. We focus on (i) the connection between the open and the closed system dynamics [P1, P2], (ii) the influence of various system parameters on the transport behavior in closed and open systems [P3, P5], and in addtion (iii) the dynamics and transport in closed models [P4, P6]. As for (i), we explore the connection of the dynamics in open systems, modeled by the Lindblad equation, and the correlation functions in the closed system [P1]. Building on this connection for magnetization transport, we study whether the time evolution of the open quantum system can be described in terms of classical correlation functions generated by the Hamiltonian equations of motion for real vectors [P2]. Second, we study in (ii) the influence of different parameters on the diffusion constants in closed and open systems [P3, P5]. In Ref. [P3], we study the impact of integrability-breaking perturbations in both closed and open spin chain systems. While we find a continuous change of the diffusion constant with perturbation strength in the closed system, in the open system, in the weak-perturbation limit, we indicate the possibility of a diverging diffusion constant. Also, we investigate the influence of the system-bath coupling on the transport behavior in boundary-driven open systems [P5]. We observe that the diffusion constants depend on the system-bath coupling and that in the nonequilibrium steady state, they inevitably retain finite-size effects, limiting the computation of long-time transport coefficients. Additional for (iii), in Ref. [P4], we first examine the impact of perturbations on the dynamics in closed spin chain models. We employ projector-operator techniques to understand, on perturbative grounds, the nontrivial damping of the standard correlation function. Second, we also discuss a dynamical quantum typicality approach for studying the energy transport in bipartite systems at low temperatures [P6].
This thesis is concerned with open quantum systems, and more specifically, quantum impurity models. By this we mean a small local quantum system in contact with a macroscopic non-interacting environment. This can be used to model individual impurities in metals and quantum information systems where the influence from the surrounding environment is not negligible. The numerical renormalization group (NRG) is the traditional method to study quantum impurity models. However its application is limited when dealing with real-time dynamics and bosonic systems. In recent years some of NRG's techniques have been introduced to the density matrix renormalization group method (DMRG), which itself is the most powerful numerical method to study one-dimensional quantum systems. The resulting method shows great potential, and this thesis explores and extends the power of the NRG+DMRG combination in treating open quantum systems. We focus mainly on two types of problem. The first is an open quantum system with a time-dependent Hamiltonian, which for example could be the theoretical description of various problems encountered in qubit manipulation. We combine NRG discretization and adaptive time-dependent DMRG (t-DMRG) to study the dissipative Landau-Zener problem. We also use this method to study the quantum decoherence process and the dynamical properties of the telegraph noise model. The results show that the NRG and t-DMRG combination is a fast, accurate and versatile method for such problems. The second type of problem we study involves the quantum critical properties of one- and two-bath spin-boson models. Unlike fermion and spin models, models with bosons are difficult to treat numerically as each boson basis has an infinite number of dimensions. By introducing the optimal boson basis into the variational matrix product state method (VMPS), which is a variant of DMRG, we can deal with an effective local boson basis as large as 1010. This is the crucial improvement over NRG which can only deal with at most a few dozen local basis states. With this powerful tool we have settled a controversy about the nature of the quantum phase transition in a sub-Ohmic spin-boson model regarding quantum to classical mapping. There, NRG fails to yield the right result due to the highly truncated local boson basis. We also explore the phase diagram of the two-bath spin-boson model and find a new critical phase. We demonstrate that NRG+VMPS with optimal boson basis represents a powerful setting to study quantum impurity models with a bosonic environment.
Using the adaptive time-dependent density matrix renormalization group method, we numerically study the spin dynamics and transport in one-dimensional spin-1/2 systems at zero temperature. Instead of computing transport coefficients from linear-response theory, we study the real-time evolution of the magnetization starting from spatially inhomogeneous initial states. In particular, we are able to analyze systems far away from equilibrium with this setup. By computing the time dependence of the variance of the magnetization, we can distinguish diffusive from ballistic regimes, depending on model parameters. For the example of the anisotropic spin-1/2 chain and at half filling, we find the expected ballistic behavior in the easy-plane phase, while in the massive regime the dynamics of the magnetization is diffusive. Our approach allows us to tune the deviation of the initial state from the ground state and the qualitative behavior of the dynamics turns out to be valid even for highly perturbed initial states in the case of easy-plane exchange anisotropies. We further cover two examples of nonintegrable models, the frustrated chain and the two-leg spin ladder, and we encounter diffusive transport in all massive phases. In the former system, our results indicate ballistic behavior in the critical phase. We propose that the study of the time dependence of the spatial variance of particle densities could be instrumental in the characterization of the expansion of ultracold atoms in optical lattices as well.
High-temperature spin transport in integrable quantum spin chains exhibits a rich dynamical phase diagram, including ballistic, superdiffusive, and diffusive regimes. While integrability is known to survive in static and periodically driven systems, its fate in the complete absence of time-translation symmetry, particularly in interacting random quantum circuits, has remained unclear. Here we construct integrable random quantum circuits built from inhomogeneous XXZ R-matrices. Remarkably, integrability is preserved for arbitrary sequences of gate layers, ranging from quasiperiodic to fully random, thereby explicitly breaking both continuous and discrete time-translation symmetry. Using large-scale time-dependent density-matrix renormalization group simulations at infinite temperature and half filling, we map out the resulting spin-transport phase diagram and identify ballistic, superdiffusive, and diffusive regimes controlled by the spectral parameters of the R-matrices. The spatiotemporal structure of spin correlations within each regime depends sensitively on the inhomogeneity, exhibiting spatial asymmetry and sharp peak structures tied to near-degenerate quasiparticle velocities. To account for these findings, we develop a generalized hydrodynamics framework adapted to time-dependent integrable circuits, yielding Euler-scale predictions for correlation functions, Drude weights, and diffusion bounds. This approach identifies the quasiparticles governing transport and quantitatively captures both the scaling exponents and fine structures of the correlation profiles observed numerically. Our results demonstrate that exact Yang-Baxter integrability is compatible with stochastic quantum dynamics and establish generalized hydrodynamics as a predictive framework for transport in time-dependent integrable systems.