The description of quantum many-body systems is an exceedingly cumbersome challenge since the dimensions of their Hilbert spaces scale exponentially with the number of constituents. Whereas mean-field approaches are expected to provide accurate descriptions for high dimensional lattices, numerical approaches based on the density matrix renormalization group have very successfully been applied to one-dimensional chains. Yet for two-dimensional lattices there is still a pressing need for efficient descriptions as these are expected to host some very intriguing but still illusive phenomena in quantum many-body physics. Here we derive an exact equation of motion for the reduced density matrices of individual subsystems of quantum many-body systems of any lattice dimension and arbitrary system size. Our projection operator based approach is thus highly efficient and not restricted to unitary dynamics. We show its excellent accuracy for moderate size one-dimensional and small two-dimensional systems where very accurate or exact numerical solutions are available.
The description of quantum many-body systems is an exceedingly cumbersome challenge since the dimensions of their Hilbert spaces scale exponentially with the number of constituents. Whereas mean-field approaches are expected to provide accurate descriptions for high dimensional lattices, numerical approaches based on the density matrix renormalization group have very successfully been applied to one-dimensional chains. Yet for two-dimensional lattices there is still a pressing need for efficient descriptions as these are expected to host some very intriguing but still illusive phenomena in quantum many-body physics. Here we derive an exact equation of motion for the reduced density matrices of individual subsystems of quantum many-body systems of any lattice dimension and arbitrary system size. Our projection operator based approach is thus highly efficient and not restricted to unitary dynamics. We show its excellent accuracy for moderate size one-dimensional and small two-dimensional systems where very accurate or exact numerical solutions are available.
We derive an exact equation of motion for the reduced density matrices of individual subsystems of quantum many-body systems of any lattice dimension and arbitrary system size. Our projection operator based theory yields a highly efficient analytical and numerical approach. Besides its practical use it provides an interpretation and systematic extension of mean-field approaches and an adaption of open quantum systems theory to settings where a dynamically evolving environment has to be taken into account. We show its high accuracy for two significant classes of complex quantum many-body dynamics, unitary evolutions of nonequilibrium states in closed and stationary states in driven-dissipative systems.
Abstract The physical properties of a quantum many-body system can, in principle, be determined by diagonalizing the respective Hamiltonian, but the dimensions of its matrix representation scale exponentially with the number of degrees of freedom. Hence, only small systems that are described through simple models can be tackled via exact diagonalization. To overcome this limitation, numerical methods based on the renormalization group paradigm that restrict the quantum many-body problem to a manageable subspace of the exponentially large full Hilbert space have been put forth. A striking example is the density-matrix renormalization group (DMRG), which has become the reference numerical method to obtain the low-energy properties of one-dimensional quantum systems with short-range interactions. Here, we provide a pedagogical introduction to DMRG, presenting both its original formulation and its modern tensor-network-based version. This colloquium sets itself apart from previous contributions in two ways. First, didactic code implementations are provided to bridge the gap between conceptual and practical understanding. Second, a concise and self-contained introduction to the tensor-network methods employed in the modern version of DMRG is given, thus allowing the reader to effortlessly cross the deep chasm between the two formulations of DMRG without having to explore the broad literature on tensor networks. We expect this pedagogical review to find wide readership among students and researchers who are taking their first steps in numerical simulations via DMRG. Graphic abstract
The analytical, experimental, and even numerical investigation of quantum many-body systems is pushed to its limits rather rapidly as the dimension of the Hilbert space scales exponentially with the number of subsystems. As a result, tasks such as the exact diagonalization of general quantum many-body Hamiltonians and the analysis of quantum experiments are limited to small system sizes. In this thesis, we address two issues related to the exponential scaling of the Hilbert space dimension: (i) The approximation of the ground state energy of many-body systems with a lower bound and (ii) efficient strategies for quantum state tomography experiments.
Understanding the non-equilibrium quantum dynamics of many-body systems is one of the most challenging problems in modern theoretical physics. While numerous approximate and exact solutions exist for systems in equilibrium, examples of non-equilibrium dynamics of many-body systems that allow reliable theoretical analysis are few and far between. In this paper, we discuss a broad class of time-dependent interacting systems subject to external linear and parabolic potentials, for which the many-body Schrödinger equation can be solved using a scaling transformation. We demonstrate that scaling solutions exist for both local and non-local interactions, and derive appropriate self-consistency equations. We apply this approach to several specific experimentally relevant examples of interacting bosons in one and two dimensions. As an intriguing result, we find that weakly and strongly interacting Bose gases expanding from a parabolic trap can exhibit very similar dynamics.
This thesis contributes to the field of strongly correlated electron systems with studies in two distinct fields thereof: the specific nature of correlations between electrons in one dimension and quantum quenches in quantum impurity problems. In general, strongly correlated systems are characterized in that their physical behaviour needs to be described in terms of a many-body description, i.e. interactions correlate all particles in a complex way. The challenge is that the Hilbert space in a many-body theory is exponentially large in the number of particles. Thus, when no analytic solution is available - which is typically the case - it is necessary to find a way to somehow circumvent the problem of such huge Hilbert spaces. Therefore, the connection between the two studies comes from our numerical treatment: they are tackled by the density matrix renormalization group (DMRG) [1] and the numerical renormalization group (NRG) [2], respectively, both based on matrix product states. The first project presented in this thesis addresses the problem of numerically finding the dominant correlations in quantum lattice models in an unbiased way, i.e. without using prior knowledge of the model at hand. A useful concept for this task is the correlation density matrix (CDM) [3] which contains all correlations between two clusters of lattice sites. We show how to extract from the CDM, a survey of the relative strengths of the system’s correlations in different symmetry sectors as well as detailed information on the operators carrying long-range correlations and the spatial dependence of their correlation functions. We demonstrate this by a DMRG study of a one-dimensional spinless extended Hubbard model [4], while emphasizing that the proposed analysis of the CDM is not restricted to one dimension. The second project presented in this thesis is motivated by two phenomena under ongoing experimental and theoretical investigation in the context of quantum impurity models: optical absorption involving a Kondo exciton [5, 6, 7] and population switching in quantum dots [8, 9, 10, 11, 12, 13, 14, 15]. It turns out that both phenomena rely on the various manifestations of Anderson orthogonality (AO) [16], which describes the fact that the response of the Fermi sea to a quantum quench (i.e. an abrupt change of some property of the impurity or quantum dot) is a change of the scattering phase shifts of all the single-particle wave functions, therefore drastically changing the system. In this context, we demonstrate that NRG, a highly accurate method for quantum impurity models, allows for the calculation of all static and dynamic quantities related to AO and present an extensive NRG study for population switching in quantum dots.