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:2104.13483

Particle Number Conservation and Block Structures in Matrix Product States

Markus Bachmayr, Michael Götte, Pfeffer, Max

Published Apr 27, 2021• 0 citations
arXivHTMLPDF

Abstract

The eigenvectors of the particle number operator in second quantization are characterized by the block sparsity of their matrix product state representations. This is shown to generalize to other classes of operators. Imposing block sparsity yields a scheme for conserving the particle number that is commonly used in applications in physics. Operations on such block structures, their rank truncation, and implications for numerical algorithms are discussed. Explicit and rank-reduced matrix product operator representations of one- and two-particle operators are constructed that operate only on the non-zero blocks of matrix product states.

DOI
10.48550/arxiv.2104.13483

Related papers

Particle number conservation and block structures in matrix product states

Markus Bachmayr, M Götte, Max J. Pfeffer

Abstract The eigenvectors of the particle number operator in second quantization are characterized by the block sparsity of their matrix product state representations. This is shown to generalize to other classes of operators. Imposing block sparsity yields a scheme for conserving the particle number that is commonly used in applications in physics. Operations on such block structures, their rank truncation, and implications for numerical algorithms are discussed. Explicit and rank-reduced matrix product operator representations of one- and two-particle operators are constructed that operate only on the non-zero blocks of matrix product states.

CALCOLO•May 25, 2022•6 citations

Low-rank eigenvalue solvers for block-sparse matrix product states

Markus Bachmayr, Sebastian Krämer, Max J. Pfeffer

We consider an iterative eigensolver for Schrödinger equations that constructs low-rank approximations of eigenfunctions with accuracy-adapted ranks, with particular focus on fermionic Schrödinger equations in second-quantized form and on matrix product state approximations enforcing particle number conservation. We provide a complete analysis of a solver based on preconditioned inverse iteration combined with rank truncation and propose a generalization to subspace iteration for the joint approximation of several eigenspaces. The practical performance of the method is illustrated by numerical tests for several model problems.

arXiv:2604.16118•Apr 17, 2026•0 citations

Particle number conservation in quantum many-body simulations with matrix product operators

Dominik Muth

Incorporating conservation laws explicitly into matrix product states (MPSs) has proven to make numerical simulations of quantum many-body systems much less resource consuming. We will discuss here to what extent this concept can be used in simulation where the dynamically evolving entities are matrix product operators (MPOs). Quite counter-intuitively the expectation of gaining in speed by sacrificing information about all but a single symmetry sector is not in all cases fulfilled. It turns out that often the entanglement imposed by the global constraint of fixed particle number is the limiting factor.

arXiv:1104.1611•Journal of Statistical Mechanics Theory and Experiment•Nov 25, 2011•13 citations

Beyond Sparsity: Quantum Block Encoding for Dense Matrices via Hierarchically Low Rank Compression

Kun Tang, Jun Lai

While quantum algorithms for solving large scale systems of linear equations offer potentially exponential speedups, their application has largely been confined to sparse matrices. This work extends the scope of these algorithms to a broad class of structured dense matrices arise in potential theory, covariance modeling, and computational physics, namely, hierarchically block separable (HBS) matrices. We develop two distinct methods to make these systems amenable to quantum solvers. The first is a pre-processing approach that transforms the dense matrix into a larger but sparse format. The second is a direct block encoding scheme that recursively constructs the necessary oracles from the HBS structure. We provide a detailed complexity analysis and rigorous error bounds for both methods. Numerical experiments are presented to validate the effectiveness of our approaches.

quant-pharXiv:2602.09745•Feb 11, 2026•0 citations

Space Reduction in Matrix Product State

Lihua Wang, Kwang S. Kim

We reconstruct a matrix product state (MPS) in reduced spaces using density matrix. This scheme applies to a MPS built on a blocked quantum lattice. Each block contains $N$ physical sites that have a local space of rank $R$. The simulation in the original spaces of rank $R^N$ is used to construct density matrices for every block. They are diagonalized and only the eigenvectors corresponding to significant diagonal elements are used to transform the original spaces to smaller ones and to reconstruct the MPS in those smaller spaces accordingly. Simulations in the reduced spaces are used to reliably extrapolate the result in unreduced spaces. Moreover, to obtain a required accuracy, the ratio of the reduced space rank over the original decreases quickly with $N$. The reduced space has a saturated rank to obtain a demanded accuracy when $N\rightarrow \infty$.

arXiv:1803.04844•Mar 13, 2018•0 citations

Block-diagonalization in second quantization

A. Tarantelli, L. S. Cederbaum

The unitary matrix that brings a Hermitian H into block-diagonal form can be uniquely determined under very simple and transparent conditions. Int his work the block-diagonalization problem is investigated in the framework of the second quantization formalism. Starting with an operator Ḩ which in any n-particle Fock space has a well-defined matrix representation an attemt was made to answer the question whether the transformation matrices T which can be separately given in the various n-particle spaces can be considered matrix representations of the same operator. T̂. Interestingly, the very important result was reached that the block-diagonalization operator T̂ exists and is unique. As a particular example, attention was concentrated on the case of an operator Ĥ given by a one-particle operator. In this case the block-diagonalization operator can be constructed and given in explicit from. This approach is applied to the theory of Green’s functions where the block-diagonalization of the Hamiltonian has interesting consequences that are illustrated in some details.

Journal of Mathematical Physics•Apr 1, 1990•3 citations
The eigenvectors of the particle number operator in second quantization are characterized by the block sparsity of their matrix product state representations. This is shown to generalize to other classes of operators. Imposing block sparsity yields a scheme for conserving the particle number that is commonly used in applications in physics. Operations on such block structures, their rank truncation, and implications for numerical algorithms are discussed. Explicit and rank-reduced matrix product operator representations of one- and two-particle operators are constructed that operate only on the non-zero blocks of matrix product states.
Abstract The eigenvectors of the particle number operator in second quantization are characterized by the block sparsity of their matrix product state representations. This is shown to generalize to other classes of operators. Imposing block sparsity yields a scheme for conserving the particle number that is commonly used in applications in physics. Operations on such block structures, their rank truncation, and implications for numerical algorithms are discussed. Explicit and rank-reduced matrix product operator representations of one- and two-particle operators are constructed that operate only on the non-zero blocks of matrix product states.
We consider an iterative eigensolver for Schrödinger equations that constructs low-rank approximations of eigenfunctions with accuracy-adapted ranks, with particular focus on fermionic Schrödinger equations in second-quantized form and on matrix product state approximations enforcing particle number conservation. We provide a complete analysis of a solver based on preconditioned inverse iteration combined with rank truncation and propose a generalization to subspace iteration for the joint approximation of several eigenspaces. The practical performance of the method is illustrated by numerical tests for several model problems.
Incorporating conservation laws explicitly into matrix product states (MPSs) has proven to make numerical simulations of quantum many-body systems much less resource consuming. We will discuss here to what extent this concept can be used in simulation where the dynamically evolving entities are matrix product operators (MPOs). Quite counter-intuitively the expectation of gaining in speed by sacrificing information about all but a single symmetry sector is not in all cases fulfilled. It turns out that often the entanglement imposed by the global constraint of fixed particle number is the limiting factor.
While quantum algorithms for solving large scale systems of linear equations offer potentially exponential speedups, their application has largely been confined to sparse matrices. This work extends the scope of these algorithms to a broad class of structured dense matrices arise in potential theory, covariance modeling, and computational physics, namely, hierarchically block separable (HBS) matrices. We develop two distinct methods to make these systems amenable to quantum solvers. The first is a pre-processing approach that transforms the dense matrix into a larger but sparse format. The second is a direct block encoding scheme that recursively constructs the necessary oracles from the HBS structure. We provide a detailed complexity analysis and rigorous error bounds for both methods. Numerical experiments are presented to validate the effectiveness of our approaches.
We reconstruct a matrix product state (MPS) in reduced spaces using density matrix. This scheme applies to a MPS built on a blocked quantum lattice. Each block contains NNN physical sites that have a local space of rank RRR. The simulation in the original spaces of rank RNR^NRN is used to construct density matrices for every block. They are diagonalized and only the eigenvectors corresponding to significant diagonal elements are used to transform the original spaces to smaller ones and to reconstruct the MPS in those smaller spaces accordingly. Simulations in the reduced spaces are used to reliably extrapolate the result in unreduced spaces. Moreover, to obtain a required accuracy, the ratio of the reduced space rank over the original decreases quickly with NNN. The reduced space has a saturated rank to obtain a demanded accuracy when N→∞N\rightarrow \inftyN→∞.
The unitary matrix that brings a Hermitian H into block-diagonal form can be uniquely determined under very simple and transparent conditions. Int his work the block-diagonalization problem is investigated in the framework of the second quantization formalism. Starting with an operator Ḩ which in any n-particle Fock space has a well-defined matrix representation an attemt was made to answer the question whether the transformation matrices T which can be separately given in the various n-particle spaces can be considered matrix representations of the same operator. T̂. Interestingly, the very important result was reached that the block-diagonalization operator T̂ exists and is unique. As a particular example, attention was concentrated on the case of an operator Ĥ given by a one-particle operator. In this case the block-diagonalization operator can be constructed and given in explicit from. This approach is applied to the theory of Green’s functions where the block-diagonalization of the Hamiltonian has interesting consequences that are illustrated in some details.