In this work we formulate an efficient method for the description of fully many-body localized systems in weak contact with thermal environments at temperature T. The key idea is to exploit the representation of the system in terms of quasi-local integrals of motion (l-bits) to efficiently derive the generator for the quantum master equation in Born-Markov approximation. We, moreover, show how to compute the steady state of this equation efficiently by using quantum-jump Monte-Carlo techniques as well as by deriving approximate kinetic equations of motion. As an example, we consider a one-dimensional disordered extended Hubbard model for spinless fermions, for which we derive the l-bit representation approximately by employing a recently proposed method valid in the limit of strong disorder and weak interactions. Coupling the system to a global thermal bath, we study the transport between two leads with different chemical potentials at both of its ends. We find that the temperature-dependent current is captured by an interaction-dependent version of Mott's law for variable range hopping, where transport is enhanced/lowered depending on whether the interactions are attractive or repulsive, respectively. We interpret these results in terms of spatio-energetic correlations between the l-bits.
In this work we formulate an efficient method for the description of fully many-body localized systems in weak contact with thermal environments at temperature T. The key idea is to exploit the representation of the system in terms of quasi-local integrals of motion (l-bits) to efficiently derive the generator for the quantum master equation in Born-Markov approximation. We, moreover, show how to compute the steady state of this equation efficiently by using quantum-jump Monte-Carlo techniques as well as by deriving approximate kinetic equations of motion. As an example, we consider a one-dimensional disordered extended Hubbard model for spinless fermions, for which we derive the l-bit representation approximately by employing a recently proposed method valid in the limit of strong disorder and weak interactions. Coupling the system to a global thermal bath, we study the transport between two leads with different chemical potentials at both of its ends. We find that the temperature-dependent current is captured by an interaction-dependent version of Mott's law for variable range hopping, where transport is enhanced/lowered depending on whether the interactions are attractive or repulsive, respectively. We interpret these results in terms of spatio-energetic correlations between the l-bits.
Abstract In this work we formulate an efficient method for the description of fully many-body localized systems in weak contact with thermal environments at temperature T. The key idea is to exploit the representation of the system in terms of quasi-local integrals of motion (l-bits) to efficiently derive the generator for the quantum master equation in Born–Markov approximation. We, moreover, show how to compute the steady state of this equation efficiently by using quantum-jump Monte-Carlo techniques as well as by deriving approximate kinetic equations of motion. As an example, we consider a one-dimensional disordered extended Hubbard model for spinless fermions, for which we derive the l-bit representation approximately by employing a recently proposed method valid in the limit of strong disorder and weak interactions. Coupling the system to a global thermal bath, we study the transport between two leads with different chemical potentials at both of its ends. We find that the temperature-dependent current is captured by an interaction-dependent version of Mott’s law for variable range hopping, where transport is enhanced/lowered depending on whether the interactions are attractive or repulsive, respectively. We interpret these results in terms of spatio-energetic correlations between the l-bits.
We investigate heat transport in a spin-1/2 Heisenberg chain, coupled locally to independent thermal baths of different temperature. The analysis is carried out within the framework of the theory of open systems by means of appropriate quantum master equations. The standard microscopic derivation of the weak-coupling Lindblad equation in the secular approximation is considered, and shown to be inadequate for the description of stationary nonequilibrium properties like a nonvanishing energy current. Furthermore, we derive an alternative master equation that is capable of describing a stationary energy current and, at the same time, leads to a completely positive dynamical map. This paves the way for efficient numerical investigations of heat transport in larger systems based on Monte Carlo wave function techniques.
The interplay between interactions and quenched disorder can result in rich dynamical quantum phenomena far from equilibrium, particularly when many-body localization prevents the system from full thermalization. With the aim of tackling this interesting regime, here we develop a semianalytical flow equation approach to study the time evolution of strongly disordered interacting quantum systems. We apply this technique to a prototype model of interacting spinless fermions in a random on-site potential in both one and two dimensions. Key results include (i) an explicit construction of the local integrals of motion that characterize the many-body localized phase in one dimension, ultimately connecting the microscopic model to phenomenological descriptions, (ii) calculation of these quantities in two dimensions, and (iii) an investigation of the real-time dynamics in the localized phase which reveals the crucial role of l-bit interactions for enhancing dephasing and relaxation.
In various fields from quantum physics to biology, the open quantum dynamics of a system consisting of interacting subsystems emphasizes its fundamental functionality. The local approach, deriving a dissipator in a master equation by ignoring the inter-subsystem interaction, has been widely used to describe the reduced dynamics due to its robustness to keep the positivity of a density operator. However, one critique is that a stationary state obtained by the approach in the limit of weak system-environment coupling is written in the form of the Gibbs state for the partial Hamiltonian by excluding the inter-subsystem interaction from the total one of the relevant system. As an alternative, the global approach, deriving a dissipator with including the inter-subsystem interaction, under the Born--Markov and secular approximations has attracted much attention, and there is debate concerning its violation of positivity in the short-term region and/or limited parameter region for the Bohr frequencies of the subsystems. In this study, we present a formalism that leads to the time-convolutionless (time-local) master equation obtained by extending the global approach beyond the Born-Markov and secular approximations. We apply it to the excitation energy transfer between interacting sites in which only the terminal site weakly interacts with a bosonic environment of finite temperature in a manner beyond the rotating-wave approximation. We find that the formulation (1) gives the short-time behavior while preserving positivity, (2) shows the oscillatory features that the secular approximation would obscure, and (3) leads to a stationary state very near to the Gibbs state for the total Hamiltonian of the relevant system.
In various fields from quantum physics to biology, the open quantum dynamics of a system consisting of interacting subsystems emphasizes its fundamental functionality. The local approach, deriving a dissipator in a master equation by ignoring the intersubsystem interaction, has been widely used to describe the reduced dynamics due to its robustness to keep the positivity of a density operator. However, one critique is that a stationary state obtained by the approach in the limit of weak system-environment coupling is written in the form of the Gibbs state for the partial Hamiltonian by excluding the intersubsystem interaction from the total one of the relevant system. As an alternative, the global approach, deriving a dissipator with including the intersubsystem interaction, under the Born-Markov and secular approximations has attracted much attention, and there is debate concerning its violation of positivity in the short-time region and/or limited parameter region for the Bohr frequencies of the subsystems. In this paper, we present a formalism that leads to the time-convolutionless (time-local) master equation obtained by extending the global approach beyond the Born-Markov and secular approximations. We apply it to the excitation energy transfer between interacting sites in which only the terminal site weakly interacts with a bosonic environment of finite temperature in a manner beyond the rotating-wave approximation. We find that the formulation (1) gives the short-time behavior while preserving positivity, (2) shows the oscillatory features that the secular approximation would obscure, and (3) leads to a stationary state very near to the Gibbs state for the total Hamiltonian of the relevant system.
Despite the exponentially large amount of information required in the quantum description of many-body systems, finite size numerical simulations have paved the way for recent progress in the understanding of quantum phases of matter both in and out-of-equilibrium. \n \nIn this spirit, using exact diagonalization for a small number of spins ~ 20, first, we investigate the fate of both ergodic quantum many-body systems that thermalize at long times and systems with quenched disorder that exhibit many-body localization (MBL) and fail to thermalize (due to an extensive set of local conserved degrees of freedom) under periodic driving. Ergodic systems always delocalize in energy space and heat up to infinite temperature, for both local driving, in which the time dependence of the Hamiltonian is restricted to a few contiguous degrees of freedom, and global driving. On the other hand, MBL systems remain localized at finite energy density for local driving, while in the case of global driving these systems can either heat up to infinite temperature or remain localized depending on the parameters of the drive. Underlying the latter is the emergence of an effective time independent MBL Hamiltonian describing the long time dynamics of the system. These numerical predictions have recently been verified experimentally in a one dimensional quasi-disordered optical lattice of interacting fermions wherein unitary time evolution can be experimentally probed. \n \nSubsequently, we study the stability of MBL systems with respect to another perturbation, in this case by thermal inclusions, which remains poorly understood. We consider a simple model to address this question: a two level system interacting with 𝑁 localized bits subject to random fields. In the large 𝑁 limit, the localized bits are thermalized for any nonzero ratio of interactions to disorder strength. Moreover, at finite 𝑁, in the transition region, the single-site eigenstate entanglement entropies exhibit bimodal distributions, so that localized bits are either "on" (strongly entangled) or "off" (weakly entangled) in eigenstates. The clusters of "on" bits vary significantly between eigenstates of the same sample, which provides evidence for a heterogenous discontinuous transition out of the localized phase in single site observables. These results also imply the MBL phase is unstable in systems with short-range interactions and quenched randomness in dimensions 𝑑 that are high but finite. \n \nThe second part of this thesis is motivated by the recent breakthroughs of artificial intelligence in domains such as computer vision involving highly dimensional data of complexity not unlike the one found in the description of many-body systems. Using Monte Carlo simulations to generate synthetic data, we find that kernel methods, a class of machine learning algorithms, can learn the mathematical \nform of physical discriminators associated with phase transitions, such as order parameters and Hamiltonian constraints, for a set of two dimensional spin models: the ferromagnetic Ising model, a conserved-order-parameter Ising model, and the Ising gauge theory. \n \nOur work shows that the combination of numerical simulations with the ability of machine learning algorithms for automatic feature detection hold promise to become standard tools in the exploration and characterization of different phases of matter in condensed matter and other many-body systems.