Abstract The iterative methods to diagonalize matrices and many-body Hamiltonians can be reformulated as flows of Hamiltonians towards diagonalization driven by unitary transformations that preserve the spectrum. After a comparative overview of the various types of discrete flows (Jacobi, QR-algorithm) and differential flows (Toda, Wegner, White) that have been introduced in the past, we focus on the random XXZ chain with random fields in order to determine the best closed flow within a given subspace of running Hamiltonians. For the special case of the free-fermion random XX chain with random fields, the flow coincides with the Toda differential flow for tridiagonal matrices which is related to the classical integrable Toda chain and which can be seen as the continuous analog of the discrete QR-algorithm. For the random XXZ chain with random fields that displays a many-body-localization transition, the present differential flow should be an interesting alternative to compare with the discrete flow that has been proposed recently to study the many-body-localization properties in a model of interacting fermions (Rademaker and Ortuno 2016 Phys. Rev. Lett. 116 , 010404 ).
Abstract The iterative methods to diagonalize matrices and many-body Hamiltonians can be reformulated as flows of Hamiltonians towards diagonalization driven by unitary transformations that preserve the spectrum. After a comparative overview of the various types of discrete flows (Jacobi, QR-algorithm) and differential flows (Toda, Wegner, White) that have been introduced in the past, we focus on the random XXZ chain with random fields in order to determine the best closed flow within a given subspace of running Hamiltonians. For the special case of the free-fermion random XX chain with random fields, the flow coincides with the Toda differential flow for tridiagonal matrices which is related to the classical integrable Toda chain and which can be seen as the continuous analog of the discrete QR-algorithm. For the random XXZ chain with random fields that displays a many-body-localization transition, the present differential flow should be an interesting alternative to compare with the discrete flow that has been proposed recently to study the many-body-localization properties in a model of interacting fermions (Rademaker and Ortuno 2016 Phys. Rev. Lett. 116 , 010404 ).
Continuous unitary transformations, or flow equations, provide a powerful non-perturbative framework for diagonalizing quantum many-body Hamiltonians. This work presents a first-principles proof that the complex, high-dimensional dynamics generated by a Wegner-type flow equation can, under specific and well-defined conditions, rigorously collapse to an exactly integrable second-order ordinary differential equation. We demonstrate that for a one-dimensional spinless fermion model in the non-interacting limit (Δ₀=0), the flow equations for the coupling constants can be derived exactly, revealing a hidden mathematical structure. We term this phenomenon "Emergent Integrability," where the diagonalization flow itself possesses an integrable submanifold, even if the initial Hamiltonian is not integrable. This discovery, established through rigorous and unabridged derivations, reveals that the process of diagonalization can be fundamentally simpler than the initial problem, shedding new light on the intrinsic mathematical structure of renormalization group-like transformations.
This paper focuses on developing a new algorithm (En-DMRG) for obtaining highly excited states of disordered and strongly interacting spin systems. This line of research is stimulated by the growing interest in understanding the many-body localized (MBL) state, which emerges in disordered interacting systems at high energy densities through a disorder-driven dynamic phase transition. By studying the one-dimensional Heisenberg spin chain (with up to 72 spins) in a random field, the authors demonstrate the accuracy of the method in obtaining energy eigenstates and the corresponding statistical results of the quantum system in the MBL phase. The method presented in this work has the potential to be applicable in a variety of quantum statistical physics problems, possibly including the MBL phase in higher dimensions.
Abstract The Vosk–Altman strong disorder renormalization for the unitary dynamics of various random quantum spin chains is reformulated as follows: the local degree of freedom characterized by the highest eigenfrequency Ω can be considered as a high-frequency-Floquet-periodic-driving for the neighboring slower degrees of freedom. Then the two first orders of the high-frequency expansion for the effective Floquet Hamiltonian can be used to generate the emergent local integrals of motion (LIOMs) and to derive the renormalization rules for the effective dynamics of the remaining degrees of freedom. The flow for this effective Floquet Hamiltonian is equivalent to the RSRG-X procedure to construct the whole set of eigenstates that generalizes the Fisher RSRG procedure constructing the ground state. This general framework is applied to the random-transverse-field XXZ spin chain in its many-body-localized phase, in order to derive the renormalization rules associated to the elimination of the biggest transverse field and to the elimination of the biggest coupling respectively.
The spectral properties of the Heisenberg spin-1/2 chain with random fields are analyzed in light of recent works on the renormalization-group flow of the Anderson model in infinite dimension. We reconstruct the β function of the order parameter from the numerical data, and observe that it may not admit a one-parameter scaling form and a simple Wilson-Fisher fixed point. Rather, it appears to be more compatible with a two-parameter, Berezinskii–Kosterlitz-Thouless-like flow with a line of fixed points (the many-body localized phase) terminating at the localization transition critical point. We argue that this renormalization group framework provides a more coherent and intuitive explanation of numerical data, up to the system sizes available with the present technology.
We benchmark a variant of the recently introduced density matrix renormalization group (DMRG)-X algorithm against exact results for the localized random field XX chain. We find that the eigenstates obtained via DMRG-X exhibit a highly accurate l-bit description for system sizes much bigger than the direct, many-body, exact diagonalization in the spin variables is able to access. We take advantage of the underlying free fermion description of the XX model to accurately test the strengths and limitations of this algorithm for large system sizes. We discuss the theoretical constraints on the performance of the algorithm from the entanglement properties of the eigenstates, and its actual performance at different values of disorder. A small but significant improvement to the algorithm is also presented, which helps significantly with convergence. We find that, at high entanglement, DMRG-X shows a bias towards eigenstates with low entanglement, but can be improved with increased bond dimension. This result suggests that one must be careful when applying the algorithm for interacting many-body localized spin models near a transition.This article is part of the themed issue 'Breakdown of ergodicity in quantum systems: from solids to synthetic matter'.
This thesis is focused on many-body localization (MBL) and the development of algorithms using the tensor networks representation for many-body localized systems. Chapter 1 is a detailed discussion on the MBL phenomenon. Chapter 2 gives an introduction to the numerical techniques that are frequently used in studying many-body quantum systems. Chapter 3 studies the shift-and-invert matrix product state method for obtaining highly excited states of MBL systems. Chapter 4 studies an efficient full diagonalization method using Wegner unitary gates. Chapter 5 studies the distribution of the entanglement entropy in the MBL-ergodic transition region. Chapter 6 studies the the influence of pseudo-spin SU (2) symmetry on the MBL phase. Chapter 7 is a side project that studies the influence of the duality-twisted boundary conditions on the critical twist defect chains. The MBL-ergodic transition is a dynamical phase transition which happens at finite energy density for a disordered and isolated many-body interacting system. It is a natural extension of Anderson localization where the localized phase survives the inter-particle interaction. In Chapter 1, we describe its relation to quantum thermalization, typical lattice models that harbor a MBL phase, MBL system’s phenomenological description, the characteristics of the MBL phase, experimental studies of MBL phenomenon in 1D and higher dimensions, the latest developments and the open questions on the nature of the MBL-ergodic transition. In Chapter 2, we provide a short but self-contained discussion of basic numerical methods that are frequently used in studying MBL physics. We mainly focus on the exact diagonalization methods and the matrix product state (MPS) representation, and explain their advantages and disadvantages. In Chapter 3, because obtaining the interior eigenstates of MBL Hamiltonians is a key step to studying the MBL phases and the transitions, we develop algorithms that can capture individual excited states to high fidelity. Using the fact that eigenstates of many-body localized systems have area-law entanglement and can therefore be efficiently represented as a matrix product state (MPS), we designed two algorithms that can generate excited states in the MPS representation, and use them to test the basic properties of MBL in the regime of large one-dimensional random field Heisenberg chains that were previously inaccessible due to the limitations of exact diagonalization (ED). In Chapter 4, we investigate an efficient way to diagonlize a fully many-body localized (FMBL) Hamiltonian. Using the Wegner flow technique for subsystems, we constructively build unitary gates that can greedily transform a local part of the FMBL Hamiltonian with the aim of lowering its average energy variance. Applying the Wegner unitary successively can make the system “flow” towards the diagonalized form. We find that the performance of this method is mainly controlled by the length of each individual Wegner unitary gate, and gets better with increasing length of the gate. We compare the performance of this constructive method with previous methods based on quasi-Newtonian optimization methods, and find that they have similar performance, while our constructive method is far more efficient. We also comment the possibility of constructing a tensor network state that resembles a multi-scale entanglement renormalization ansatz (MERA) using this constructive method. In Chapter 5, using strong subadditivity (SSA) theorem, we develop two order parameters – cut-averaged entanglement entropy (CAEE) and its slope (SCAEE), and show that they can be used to directly identify the volume or area law scaling in single eigenstates. We study the distribution of the SCAEE over disorder realizations. This distribution appears Gaussian at weak disorder, while on the MBL side, the distribution of the SCAEE is peaked at zero slope and has an exponential tail. In the critical regime we find that the distribution of the SCAEE is bimodal both over multiple disorder realizations as well as for single disorder realizations. The variance of the SCAEE distribution in the transition region seems to grow with system size when considered over disorder realizations. Our system sizes are too small to pin down its maximal value, but they are consistent with (among other possibilities) the maximal variance possible which would lead to half the states having zero SCAEE and half having maximal SCAEE. This scenario would lead to an entanglement entropy at the transition which scales as a volume law with half its thermal value. In Chapter 6, we consider the strongly disordered one-dimensional Hubbard model with spin disorder, which, as we argue, cannot be categorized as a conventional MBL phase. Though the spin rotation SU (2) symmetry is broken by the spin disorder, this model still preserves the pseudo-spin SU (2) symmetry under periodic boundary conditions. On the theoretical side, using the pseudo-spin algebra we show that a significant number of the excited states at any disorder strength have logarithmic correction to their von Neumann entanglement entropy, which violates the area-law entanglement for a typical MBL phase. A common feature of this group of excited states is that they usually have far more double occupancies than single occupancies. On the numerical side, we studied the time evolution after quantum quench from two extreme cases of product states – all single occupancies at quarter filling, and all double occupancies at half filling. We find in the former case entanglement entropy behaves as in a typical MBL system, while the latter is clearly delocalized, which suggests the existence of incomplete set of local integrals of motion. With the above evidence, it is convincing that this model provides a playground for studying a non-ergodic, non-MBL phase, and its relation to continuous SU (2) symmetry. Chapter 7 deviates from the common scene of MBL physics. In this chapter, we consider a series of one-dimensional critical twist defect models with various boundary conditions. These twist defects are embedded in the background of a generalization of Kitaev Z 2 toric code model, and can also be understood as a discrete Z k gauge theory in its deconfined phase. We analytically construct the duality transformation of these critical models, and use numerical methods to extract information on the central charges, sound velocities, and conformal dimensions. With the numerical data, we show the relation between the CFT contents of the critical chains with periodic boundary conditions and the critical chains with duality-twisted boundary conditions.