A fundamental longstanding problem in studying spin models is the efficient and accurate numerical simulation of the long-time behavior of larger systems. The exponential growth of the Hilbert space and the entanglement accumulation at long times pose major challenges for current methods. To address these issues, we employ the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) framework to simulate the many-body spin dynamics of the Heisenberg model in various settings, including the Ising and XYZ limits with different interaction ranges and random couplings. Benchmarks with analytical and exact numerical approaches show that ML-MCTDH accurately captures the time evolution of one- and two-body observables in both one- and two-dimensional lattices. A comparison with the discrete truncated Wigner approximation (DTWA) highlights that ML-MCTDH is particularly well-suited for handling anisotropic models and provides more reliable results for two-point observables across all tested cases. The behavior of the corresponding entanglement dynamics is analyzed to reveal the complexity of the quantum states. Our findings indicate that the rate of entanglement growth strongly depends on the interaction range and the presence of disorder. This particular relationship is then used to examine the convergence behavior of ML-MCTDH. Our results indicate that the multilayer structure of ML-MCTDH is a promising numerical framework for handling the dynamics of generic many-body spin systems.
A fundamental longstanding problem in studying spin models is the efficient and accurate numerical simulation of the long-time behavior of larger systems. The exponential growth of the Hilbert space and the entanglement accumulation at long times pose major challenges for current methods. To address these issues, we employ the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) framework to simulate the many-body spin dynamics of the Heisenberg model in various settings, including the Ising and XYZ limits with different interaction ranges and random couplings. Benchmarks with analytical and exact numerical approaches show that ML-MCTDH accurately captures the time evolution of one- and two-body observables in both one- and two-dimensional lattices. A comparison with the discrete truncated Wigner approximation (DTWA) highlights that ML-MCTDH is particularly well-suited for handling anisotropic models and provides more reliable results for two-point observables across all tested cases. The behavior of the corresponding entanglement dynamics is analyzed to reveal the complexity of the quantum states. Our findings indicate that the rate of entanglement growth strongly depends on the interaction range and the presence of disorder. This particular relationship is then used to examine the convergence behavior of ML-MCTDH. Our results indicate that the multilayer structure of ML-MCTDH is a promising numerical framework for handling the dynamics of generic many-body spin systems.
Numerical simulations of quantum spin models are crucial for a profound understanding of many-body phenomena in a variety of research areas in physics. An outstanding problem is the availability of methods to tackle systems that violate area laws of entanglement entropy. Such scenarios cover a wide range of compelling physical situations including disordered quantum spin systems among others. In this paper, we employ a numerical technique referred to as multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) to evaluate the ground state of several disordered spin models. ML-MCTDH has previously been used to study problems of high-dimensional quantum dynamics in molecular and ultracold physics but is here applied to study spin systems. We exploit the inherent flexibility of the method to present results in one and two spatial dimensions and treat challenging setups that incorporate long-range interactions as well as disorder. Our results suggest that the hierarchical multilayering inherent to ML-MCTDH allows to tackle a wide range of quantum many-body problems such as spin dynamics of varying dimensionality.
Multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) theory is a rigorous and powerful method to simulate quantum dynamics in complex many-body systems. This approach extends the original MCTDH theory of Meyer, Manthe, and Cederbaum to include dynamically contracted layers in a recursive way, within which the equations of motion are determined from the Dirac-Frenkel variational principle. This paper presents the general derivation of the theory and analyzes the important features that make the ML-MCTDH method numerically efficient. Furthermore, we discuss the generalization of the theory to treat many-body identical particles (fermions or bosons) as well as calculating energy eigenstates via the improved relaxation method.
The dynamics of a two-level system coupled to a spin bath is investigated via the numerically exact multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) theory. Consistent with the previous work on linear response approximation [N. Makri, J. Phys. Chem. B 103, 2823 (1999)], it is demonstrated numerically that this spin-spin-bath model can be mapped onto the well-known spin-boson model if the system-bath coupling strength obeys an appropriate scaling behavior. This linear response mapping, however, may require many bath spin degrees of freedom to represent the practical continuum limit. To clarify the discrepancies resulted from different approximate treatments of this model, the population dynamics of the central two-level system has been investigated near the transition boundary between the coherent and incoherent motions via the ML-MCTDH method. It is found that increasing temperature favors quantum coherence in the nonadiabatic limit of this model, which corroborates the prediction in the previous work [J. Shao and P. Hanggi, Phys. Rev. Lett. 81, 5710 (1998)] based on the non-interacting blip approximation (NIBA). However, the coherent-incoherent boundary obtained by the exact ML-MCTDH simulation is slightly different from the approximate NIBA results. Quantum dynamics in other physical regimes are also discussed.
The improved relaxation approach of Meyer and co-workers is extended by employing the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) theory. Within this approach the trial wave function is expressed by a recursive, multilayered expansion that is the same as in a normal ML-MCTDH simulation. Calculation of an energy eigenstate then proceeds by iteratively diagonalizing the Boltzmann operator in the top layer representation using a Lanczos/Arnoldi method while relaxing the single particle functions of all layers using the ML-MCTDH imaginary time propagation. Numerical iteration is stopped until the relative energy and/or density of the target state meet the convergence criteria. The application of this multilayer improved relaxation method is illustrated by computing the energy splitting for the lowest pair of eigenstates as well as the spin observables and entanglement entropy of the spin-boson model. A large number of degrees of freedom, up to 10,000 bath modes, can be treated, and thus the appropriate phase transition at a critical system-bath coupling strength may be predicted via appropriate extrapolations.
In this Colloquium, the wavefunction-based Multiconfigurational Time-Dependent Hartree approaches to the dynamics of indistinguishable particles (MCTDH-F for Fermions and MCTDH-B for Bosons) are reviewed. MCTDH-B and MCTDH-F or, together, MCTDH-X are methods for describing correlated quantum systems of identical particles by solving the time-dependent Schrödinger equation from first principles. MCTDH-X is used to accurately model the dynamics of real-world quantum many-body systems in atomic, molecular, and optical physics. The key feature of these approaches is the time-dependence and optimization of the single-particle states employed for the construction of a many-body basis set, which yields nonlinear working equations. We briefly describe the historical developments that have lead to the formulation of the MCTDH-X methods and motivate the necessity for wavefunction-based approaches. We sketch the derivation of the unified MCTDH-F and MCTDH-B equations of motion for complete and also specific restricted configuration spaces. The strengths and limitations of the MCTDH-X approach are assessed via benchmarks against an exactly solvable model and via convergence checks. We highlight some applications to instructive and experimentally-realized quantum many-body systems: the dynamics of atoms in Bose-Einstein condensates in magneto-optical and optical traps and of electrons in atoms and molecules. We discuss the current development and frontiers in the field of MCTDH-X: theories and numerical methods for indistinguishable particles, for mixtures of multiple species of indistinguishable particles, the inclusion of nuclear motion for the nonadia-batic dynamics of atomic and molecular systems, the so-called multilayer generalizations to the MCTDH-F and MCTDH-B methods, and the time-dependent orbital-adaptive coupled cluster theory are discussed.
Abstract Numerical techniques to efficiently model out-of-equilibrium dynamics in interacting quantum many-body systems are key for advancing our capability to harness and understand complex quantum matter. Here we propose a new numerical approach which we refer to as generalized discrete truncated Wigner approximation (GDTWA). It is based on a discrete semi-classical phase space sampling and allows to investigate quantum dynamics in lattice spin systems with arbitrary S ≥ 1/2. We show that the GDTWA can accurately simulate dynamics of large ensembles in arbitrary dimensions. We apply it for S > 1/2 spin-models with dipolar long-range interactions, a scenario arising in recent experiments with magnetic atoms. We show that the method can capture beyond mean-field effects, not only at short times, but it also can correctly reproduce long time quantum-thermalization dynamics. We benchmark the method with exact diagonalization in small systems, with perturbation theory for short times, and with analytical predictions made for models which feature quantum-thermalization at long times. We apply our method to study dynamics in large S > 1/2 spin-models and compute experimentally accessible observables such as Zeeman level populations, contrast of spin coherence, spin squeezing, and entanglement quantified by single-spin Renyi entropies. We reveal that large S systems can feature larger entanglement than corresponding S = 1/2 systems. Our analyses demonstrate that the GDTWA can be a powerful tool for modeling complex spin dynamics in regimes where other state-of-the art numerical methods fail.