The numerical simulation of two-dimensional quantum many-body systems away from equilibrium constitutes a major challenge for all known computational methods. We investigate the utility of tree tensor network (TTN) states to solve the dynamics of the quantum Ising model in two dimensions. Within the perturbative regime of small transverse fields, TTNs faithfully reproduce analytically known, but nontrivial and physically interesting dynamics of magnetic domains, for lattices up to 16 × 16 sites. Limitations of the method related to the rapid growth of entanglement entropy are explored within more general, paradigmatic quench settings. We provide and discuss comprehensive benchmarks regarding the benefit of GPU acceleration and the impact of using local operator sums on the performance.
The numerical simulation of two-dimensional quantum many-body systems away from equilibrium constitutes a major challenge for all known computational methods. We investigate the utility of tree tensor network (TTN) states to solve the dynamics of the quantum Ising model in two dimensions. Within the perturbative regime of small transverse fields, TTNs faithfully reproduce analytically known, but nontrivial and physically interesting dynamics of magnetic domains, for lattices up to 16 × 16 sites. Limitations of the method related to the rapid growth of entanglement entropy are explored within more general, paradigmatic quench settings. We provide and discuss comprehensive benchmarks regarding the benefit of GPU acceleration and the impact of using local operator sums on the performance.
The discovery of dynamical quantum phase transitions (DQPTs) has fundamentally challenged the traditional view that phase transitions only occur in thermal equilibrium. Experimental platforms and 1D numerical methods, like matrix product states (MPS), have made great progress. However, exploring true 2D DQPTs remains difficult due to finite-size limitations and the geometric biases of quasi-1D cylinder mappings. Here, we bypass these limitations by deploying a tree tensor network (TTN) approach. This allows us to directly compute the quench dynamics of the transverse-field Ising model (TFIM) on an open 2D square lattice. Because the TTN architecture naturally mirrors 2D lattice connectivity, we can extract the global Loschmidt echo. Our simulations reveal that while deep quenches yield standard DQPTs, quenching within the ferromagnetic phase produces an anomalous dynamical response. In this regime, the rate function exhibits sharp non-analytic peaks even as the macroscopic order parameter maintains its initial sign. This decoupled behavior strongly indicates that local spin excitations drive 2D DQPTs, rather than the macroscopic domain-wall motions seen in 1D chains. These results provide a quantitative numerical baseline for understanding non-equilibrium quantum matter in higher dimensions.
The quantum dynamics of many-qubit systems is an outstanding problem that has recently driven significant advances in both numerical methods and programmable quantum processing units. In this work, we employ a comprehensive toolbox of state-of-the-art numerical approaches to classically simulate the dynamics of the two-dimensional transverse field Ising model. Our methods include three different tensor network techniques -- matrix product states, tree-tensor networks, and two-dimensional tensor-networks under the belief propagation approximation -- as well as time-dependent variational Monte Carlo with Neural Quantum States. We focus on two paradigmatic dynamical protocols: (i) quantum annealing through a critical point and (ii) post-quench dynamics. Our extensive results show the quantitative predictions of various state-of-the-art numerical methods providing a benchmark for future numerical investigations and experimental studies with the aim to push the limitations on classical and QPUs. In particular, our work connects classical simulability to different regimes associated with quantum dynamics in Rydberg arrays - namely, quasi-adiabatic dynamics, the Kibble-Zurek mechanism, and quantum quenches.
We investigate the disordered spin-21Heisenberg model in two dimensions and employ tree tensor networks (TTNs) with a physics-informed structural optimization of the tree layout, to simulate dynamics in the many-body localization problem. We find that TTNs are able to capture two-dimensional entanglement patterns more effectively than matrix product states (MPS) while being more efficient to contract than projected entangled pair states (PEPS) to probe larger systems and longer times. Structural optimization of the trees based on time evolution of the entanglement in the system allows to keep the necessary bond dimensions low and to maximally exploit the increased expressiveness of TTNs over MPS. In this way, we achieve more accurate results in all considered parameter regimes both below and above the ergodicity-to-localization crossover at a comparable compute-time cost.
The efficient numerical simulation of nonequilibrium real-time evolution in isolated quantum matter constitutes a key challenge for current computational methods. This holds in particular in the regime of two spatial dimensions, whose experimental exploration is currently pursued with strong efforts in quantum simulators. In this work we present a versatile and efficient machine learning inspired approach based on a recently introduced artificial neural network encoding of quantum many-body wave functions. We identify and resolve key challenges for the simulation of time evolution, which previously imposed significant limitations on the accurate description of large systems and long-time dynamics. As a concrete example, we study the dynamics of the paradigmatic two-dimensional transverse-field Ising model, as recently also realized experimentally in systems of Rydberg atoms. Calculating the nonequilibrium real-time evolution across a broad range of parameters, we, for instance, observe collapse and revival oscillations of ferromagnetic order and demonstrate that the reached timescales are comparable to or exceed the capabilities of state-of-the-art tensor network methods.
In this thesis I present contributions to the development, analysis and application of tensor network state methods for numerically exact quantum dynamics in one and two-dimensional lattice systems. The setting of numerically exact quantum dynamics is introduced in Chapter 2. This includes a discussion of exact diagonalization approaches and massively parallel implementations thereof as well as a brief introduction of tensor network states. In Chapter 3, I perform a detailed analysis of the performance of n-ary tree tensor network states for simulating the dynamics of two-dimensional lattices. This constitutes the first application of this class of tensor network to dynamics in two spatial dimensions, a long-standing challenge, and the method is found to perform on par with existing state-of-the-art approaches. Chapter 4 showcases the efficacy of a novel tensor network format I developed, tailored to electron-phonon coupled problems in their single-electron sector, through an application to the Holstein model. The applicability of the approach to a broad range of parameters of the model allows to reveal the strong influence of the spread of the electron distribution on the initial state of the phonons at the site where the electron is introduced, for which a simple physical picture is offered. I depart from method development in Chapter 5 and analyse the prospects of using tensor network states evolved using the time-dependent variational principle as an approximate approach to determine asymptotic transport properties with a finite, moderate computational effort. The method is shown to not yield the correct asymptotics in a clean, non-integrable system and can thus not be expected to work in generic systems, outside of finely tuned parameter regimes of certain models. Chapters 6 and 7 are concerned with studies of spin transport in long-range interacting systems using tensor network state methods. For the clean case, discussed in Chapter 6, we find that for sufficiently short-ranged interactions, the spreading of the bulk of the excitation is diffusive and thus dominated by the local part of the interaction, while the tail of the excitation decays with a powerlaw that is twice as large as the powerlaw of the interaction. Similarly, in the disordered case, analysed in Chapter 7, we find subdiffusive transport of spin and sub-linear growth of entanglement entropy. This behaviour is in agreement with the behaviour of systems with local interactions at intermediate disorder strength, but provides evidence against the phenomelogical Griffith picture of rare, strongly disordered insulating regions. We generalize the latter to long-ranged interactions and show that it predicts to diffusion, in contrast to the local case where it results in subdiffusive behaviour.
Tensor network states (TNSs) offer an efficient representation for the ground state of quantum many-body systems and play an important role in the simulations of them. Numerous TNSs have been proposed in the past few decades. However, due to the high cost of TNS for two-dimensional systems, a balance between the encoded entanglement and computational complexity of TNS is yet to be reached. In this work we introduce a new tree tensor network-based (TTN-based) TNS dubbed a ``fully augmented tree tensor network (FATTN) by releasing the constraint in augmented tree tensor network (ATTN). When disentanglers are augmented in the physical layer of TTN, FATTN can provide more entanglement than TTN and ATTN. At the same time, FATTN maintains the scaling of computational cost with bond dimension in TTN and ATTN. Benchmark results on the ground-state energy for the transverse Ising model are provided to demonstrate the improvement of accuracy of FATTN over TTN and ATTN. Moreover, FATTN is quite flexible, which can be constructed as an interpolation between TTN and the multiscale entanglement renormalization ansatz to reach a balance between the encoded entanglement and the computational cost.