Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2409.02362 · Journal of Physics A Mathematical and Theoretical

Bundled matrix product states represent low-energy excitations faithfully

Thomas E. Baker, Negar Seif

Published Sep 3, 2024• 0 citations
arXivHTMLPDF

Abstract

Abstract We consider a set of density matrices. All of which are written in the same orbital basis, but the orbital basis size is less than the total Hilbert space size. We ask how each density matrix is related to each of the others by establishing a norm between density matrices based on the truncation error in a partial trace for a small set of orbitals. We find that states with large energy differences must have large differences in their density matrices. Small energy differences are divided into two groups, one where two density matrices have small differences and another where they are very different, as is the case of symmetry. We extend these ideas to a bundle of matrix product states and show that bond dimension of the wavefunction ansatz for two states with large energy differences are larger. Meanwhile, low energy differences can have nearly the same bond dimensions for similar states.

DOI
10.1088/1751-8121/ad770f

Related papers

Theory of bundled matrix product states and low-lying excitations in quantum many-body systems

Negar Seif

This thesis develops a framework for understanding when different excited states of local quantum Hamiltonians can be efficiently represented in a common basis, one containing fewer states than the full Hilbert space. In the analytical part, we introduce metrics that quantify the ''closenes'' of excitations by comparing their truncated density matrices in a fixed single-particle basis, and relate these metrics to their energy differences. We show that large energy separations necessarily imply large differences between density matrices, while small energy differences fall into two categories: one in which the density matrices differ only slightly, and another in which they remain very different, as can occur in the presence of symmetry. We also extend these ideas to bundled matrix product states (MPS), in which multiple excitations are encoded simultaneously, and study how the required bond dimension grows as excitations with different energies are included.

•Jan 1, 2026•0 citations

Encoding the structure of many-body localization with matrix product operators

David Pekker, Bryan K. Clark

Anderson insulators are noninteracting disordered systems which have localized single-particle eigenstates. The interacting analog of Anderson insulators are the many-body localized (MBL) phases. The spectrum of the many-body eigenstates of an Anderson insulator is efficiently represented as a set of product states over the single-particle modes. We show that product states over matrix product operators of small bond dimension is the corresponding efficient description of the spectrum of an MBL insulator. In this language all of the many-body eigenstates are encoded by matrix product states (i.e., density matrix renormalization group wave functions) consisting of only two sets of low bond dimension matrices per site: the ${G}_{i}$ matrices corresponding to the local ground state on site $i$ and the ${E}_{i}$ matrices corresponding to the local excited state. All ${2}^{n}$ eigenstates can be generated from all possible combinations of these sets of matrices.

arXiv:1410.2224•Physical review. B./Physical review. B•Jan 10, 2017•42 citations

Geometry of density matrices. IV. The relationship between density matrices and densities

John E. Harriman

The relationship between densities and density matrices is explored in the case of a finite-basis-set expansion. The space of one-electron density matrices can be divided into two orthogonal subspaces with elements in one of them in one-to-one correspondence with densities. A component in the other does not contribute to the density. The set of densities is convex but there may be densities which cannot be obtained from a density matrix. The matrix of a local potential has a component only in the first subspace, and any such matrix can be obtained from a local potential. It is possible for Hamiltonian matrices differing by the matrix of a local potential to have a common ground-state eigenvector, so a Hohenberg-Kohn theorem cannot always be established. When it can, the explicit local potential with a given ground-state density can be formally obtained when appropriate conditions are satisfied. The details of the decomposition of the space of matrices and of subsequent developments depend on linear-dependency relationships among basis-set products, and are thus basis-set dependent.

Physical Review A•Feb 1, 1983•74 citations

Reduced Density Matrices for Valence Bond Wavefunctions. II. Symmetry and Reduction Procedures for Square Planar Nets

Jong H. Choi, R. P. Messmer, Walter R. Thorson

A systematic procedure is described for calculating the reduced density matrices of first and second order for a resonating valence bond wavefunction on square planar nets with Born-von Karman boundary conditions of finite period. The density matrices explicitly display the electron correlation in such a wavefunction, which is a model ground state for a two-dimensional Mott insulator. Because of the extremely large number of contributing valence bond structures, it is essential to utilize their symmetry under the space group for minimum computation. For the totally symmetric singlet state described, density matrix eigenvectors form basis vectors for irreducible representations. The irreducible representations spanned by pair states on the square plane lattice are listed, and a set of irreducible basis vectors for the two-particle density matrix is defined which clearly displays the nature of electron pair correlations in such a system (correlated geminals). The apparatus has been employed in a calculation described in the paper accompanying this one (succeeding article).

The Journal of Chemical Physics•Jul 1, 1972•3 citations

Two-electron wavefunctions are matrix product states with bond dimension three

Gero Friesecke, Benedikt R. Graswald

We prove the statement in the title for a suitable (wavefunction-dependent) choice of the underlying orbitals and show that 3 is optimal. Thus, for two-electron systems, the quantum chemistry density matrix renormalization group (QC-DMRG) method with bond dimension 3 combined with fermionic mode optimization exactly recovers the full configuration-interaction (FCI) energy.

arXiv:2109.10091•Journal of Mathematical Physics•Sep 1, 2022•1 citation

Correlation density matrices for one-dimensional quantum chains based on the density matrix renormalization group

Wolfgang Münder, Andreas Weichselbaum, Andreas Holzner +2 more

A useful concept for finding numerically the dominant correlations of a given ground state in an interacting quantum lattice system in an unbiased way is the correlation density matrix (CDM). For two disjoint, separated clusters, it is defined to be the density matrix of their union minus the direct product of their individual density matrices and contains all the correlations between the two clusters. We show how to extract from the CDM a survey of the relative strengths of the system's correlations in different symmetry sectors and the nature of their decay with distance (power law or exponential), as well as detailed information on the operators carrying long-range correlations and the spatial dependence of their correlation functions. To achieve this goal, we introduce a new method of analysing the CDM, termed the dominant operator basis (DOB) method, which identifies in an unbiased fashion a small set of operators for each cluster that serve as a basis for the dominant correlations of the system. We illustrate this method by analysing the CDM for a spinless extended Hubbard model that features a competition between charge density correlations and pairing correlations, and show that the DOB method successfully identifies their relative strengths and dominant correlators. To calculate the ground state of this model, we use the density matrix renormalization group, formulated in terms of a variational matrix product state (MPS) approach within which subsequent determination of the CDM is very straightforward. In an extended appendix, we give a detailed tutorial introduction to our variational MPS approach for ground state calculations for one-dimensional quantum chain models. We present in detail how MPSs overcome the problem of large Hilbert space dimensions in these models and describe all the techniques needed for handling them in practice.

arXiv:0910.0753•New Journal of Physics•Jul 28, 2010•12 citations
Abstract We consider a set of density matrices. All of which are written in the same orbital basis, but the orbital basis size is less than the total Hilbert space size. We ask how each density matrix is related to each of the others by establishing a norm between density matrices based on the truncation error in a partial trace for a small set of orbitals. We find that states with large energy differences must have large differences in their density matrices. Small energy differences are divided into two groups, one where two density matrices have small differences and another where they are very different, as is the case of symmetry. We extend these ideas to a bundle of matrix product states and show that bond dimension of the wavefunction ansatz for two states with large energy differences are larger. Meanwhile, low energy differences can have nearly the same bond dimensions for similar states.
This thesis develops a framework for understanding when different excited states of local quantum Hamiltonians can be efficiently represented in a common basis, one containing fewer states than the full Hilbert space. In the analytical part, we introduce metrics that quantify the ''closenes'' of excitations by comparing their truncated density matrices in a fixed single-particle basis, and relate these metrics to their energy differences. We show that large energy separations necessarily imply large differences between density matrices, while small energy differences fall into two categories: one in which the density matrices differ only slightly, and another in which they remain very different, as can occur in the presence of symmetry. We also extend these ideas to bundled matrix product states (MPS), in which multiple excitations are encoded simultaneously, and study how the required bond dimension grows as excitations with different energies are included.
Anderson insulators are noninteracting disordered systems which have localized single-particle eigenstates. The interacting analog of Anderson insulators are the many-body localized (MBL) phases. The spectrum of the many-body eigenstates of an Anderson insulator is efficiently represented as a set of product states over the single-particle modes. We show that product states over matrix product operators of small bond dimension is the corresponding efficient description of the spectrum of an MBL insulator. In this language all of the many-body eigenstates are encoded by matrix product states (i.e., density matrix renormalization group wave functions) consisting of only two sets of low bond dimension matrices per site: the Gi{G}_{i}Gi​ matrices corresponding to the local ground state on site iii and the Ei{E}_{i}Ei​ matrices corresponding to the local excited state. All 2n{2}^{n}2n eigenstates can be generated from all possible combinations of these sets of matrices.
The relationship between densities and density matrices is explored in the case of a finite-basis-set expansion. The space of one-electron density matrices can be divided into two orthogonal subspaces with elements in one of them in one-to-one correspondence with densities. A component in the other does not contribute to the density. The set of densities is convex but there may be densities which cannot be obtained from a density matrix. The matrix of a local potential has a component only in the first subspace, and any such matrix can be obtained from a local potential. It is possible for Hamiltonian matrices differing by the matrix of a local potential to have a common ground-state eigenvector, so a Hohenberg-Kohn theorem cannot always be established. When it can, the explicit local potential with a given ground-state density can be formally obtained when appropriate conditions are satisfied. The details of the decomposition of the space of matrices and of subsequent developments depend on linear-dependency relationships among basis-set products, and are thus basis-set dependent.
A systematic procedure is described for calculating the reduced density matrices of first and second order for a resonating valence bond wavefunction on square planar nets with Born-von Karman boundary conditions of finite period. The density matrices explicitly display the electron correlation in such a wavefunction, which is a model ground state for a two-dimensional Mott insulator. Because of the extremely large number of contributing valence bond structures, it is essential to utilize their symmetry under the space group for minimum computation. For the totally symmetric singlet state described, density matrix eigenvectors form basis vectors for irreducible representations. The irreducible representations spanned by pair states on the square plane lattice are listed, and a set of irreducible basis vectors for the two-particle density matrix is defined which clearly displays the nature of electron pair correlations in such a system (correlated geminals). The apparatus has been employed in a calculation described in the paper accompanying this one (succeeding article).
A useful concept for finding numerically the dominant correlations of a given ground state in an interacting quantum lattice system in an unbiased way is the correlation density matrix (CDM). For two disjoint, separated clusters, it is defined to be the density matrix of their union minus the direct product of their individual density matrices and contains all the correlations between the two clusters. We show how to extract from the CDM a survey of the relative strengths of the system's correlations in different symmetry sectors and the nature of their decay with distance (power law or exponential), as well as detailed information on the operators carrying long-range correlations and the spatial dependence of their correlation functions. To achieve this goal, we introduce a new method of analysing the CDM, termed the dominant operator basis (DOB) method, which identifies in an unbiased fashion a small set of operators for each cluster that serve as a basis for the dominant correlations of the system. We illustrate this method by analysing the CDM for a spinless extended Hubbard model that features a competition between charge density correlations and pairing correlations, and show that the DOB method successfully identifies their relative strengths and dominant correlators. To calculate the ground state of this model, we use the density matrix renormalization group, formulated in terms of a variational matrix product state (MPS) approach within which subsequent determination of the CDM is very straightforward. In an extended appendix, we give a detailed tutorial introduction to our variational MPS approach for ground state calculations for one-dimensional quantum chain models. We present in detail how MPSs overcome the problem of large Hilbert space dimensions in these models and describe all the techniques needed for handling them in practice.