REVIEW
Calculation of reduced density matrices from correlation functions
Condensed Matter · Submitted 30 Dec 2002 · 13 references
Pith's verdict: A short, honest methods note that gives the working route from correlation functions to entanglement spectra in free lattice models — the formula is older than the paper, but the unified derivation is the reason people cite it.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
For solvable lattice systems, the reduced density matrix of a subsystem is fixed by its ordinary two-point correlation functions, so the entanglement Hamiltonian can be read off by taking a logarithm of a single matrix.
The reading
What carries the argument
The load-bearing device is the identification between the correlation matrix of the block and the exponential form of its density matrix. Because $K e^{-H}$ with a quadratic $H$ generates Wick-compatible expectation values, matching the block's one-particle function forces the eigenvalues $H' = H$ of $C$ and $H' = H$ of $H$ to satisfy $H' = H$, which in matrix form is $H' = H$. The logarithm of a matrix built from $C$ is therefore the entanglement Hamiltonian itself, and the whole entanglement spectrum is one diagonalization away from data that can be read off ordinary correlation functions.
Extended reading notes
Core claim
The paper's central claim is that the reduced density matrix of a subsystem of a solvable fermionic or bosonic lattice system is completely determined by the subsystem's correlation functions, with no need to trace out the environment explicitly. For a free-fermion hopping Hamiltonian whose eigenstates are Slater determinants, the reduced density matrix takes the exponential form $K e^{-H}$ with $H = H$ quadratic in the subsystem's fermion operators, and the matrix $H$ is obtained from the block's one-particle correlation matrix $C_{ij} = K e^{-H}$ by the relation $H' = H$ — in the paper's notation $H' = H$ is written as $H' = H$ where $H' = H$ denotes $H' = H$. The eigenvalues of $C$ and of $H$ are tied by $H' = H$, so diagonalizing the $M imes M$ matrix $C$ for an $M$-site block gives the full entanglement spectrum. For systems with pairing, diagonalized by a Bogoliubov transformation, an additional anomalous correlation matrix is required, and the analogous relation is $H' = H$. For coupled harmonic oscillators the reduced density matrix is again Gaussian, and the single-particle entanglement energies follow from the eigenvalues $H' = H$ of the product of the position- and momentum-correlation matrices through $H' = H$.
Load-bearing premise
The argument assumes without deriving it here that the subsystem's reduced density matrix is the exponential of a quadratic fermion (or boson) operator; if that Gaussian form failed, matching the two-point correlations would not pin down the density matrix.
Editorial extensions
If this is right
- Every entanglement quantity of a block in a free-fermion or Gaussian-boson ground state — entropy, spectrum, and all Rényi entropies — can be computed by diagonalizing a single $M imes M$ correlation matrix instead of constructing and tracing out the full wavefunction.
- For models where the correlation functions have closed analytic forms, such as nearest-neighbour hopping on a chain or on the square lattice at half filling, only the diagonalization step remains between the model and its exact entanglement spectrum.
- The construction carries over to finite temperature, since the same correlation-matrix logic applies to thermal states of the same solvable systems.
- For systems with pairing, the anomalous correlations enter as a second matrix, and the combination $H' = H$ plays the role that $C$ alone plays for pure hopping.
- For a single oscillator, the entire reduced density matrix is fixed by the product of the position and momentum variances, which for a homogeneous system reduces to the product of the mean inverse frequency and the mean frequency of the normal modes.
Reading between the lines
- Because the formula returns a Gaussian density matrix for any input correlation matrix, it also defines a canonical Gaussian approximation to the reduced density matrix of an interacting system; whether that approximation captures the dominant entanglement spectrum is a testable question the paper does not pursue.
- The author notes that when the entanglement energies are large, the corresponding eigenvalues of $C$ sit exponentially close to 0 or 1, so the logarithm becomes numerically ill-conditioned exactly where the spectrum is deepest; practical implementations would do better to work with the eigenvalue relation directly in extended precision.
- The bosonic result makes the single-site entanglement spectrum of any Gaussian oscillator ground state a two-number object, which suggests a cheap diagnostic for when a harmonic lattice model is close to a critical point: the product of the two variances should diverge there.
Original abstract
It is shown that for solvable fermionic and bosonic lattice systems, the reduced density matrices can be determined from the properties of the correlation functions. This provides the simplest way to these quantities which are used in the density-matrix renormalization group method.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
It is shown that for solvable fermionic and bosonic lattice systems, the reduced density matrices can be determined from the properties of the correlation functions. This provides the simplest way to these quantities which are used in the density-matrix renormalization group method.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | It is shown that for solvable fermionic and bosonic lattice systems, the reduced density matrices can be determined from the properties of the correlation functions. | result 0.77 | 0.84 | 0.91 | asserted 1.00 |
| s1 | This provides the simplest way to these quantities which are used in the density-matrix renormalization group method. | interpretation 0.98 | 0.35 | 0.14 | asserted 1.00 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
References
How later papers used it
344 later papers cite it
16 later papers by other authors are listed below: 16 used its method or tool.
Uses its method or tool 16 of 205 papers with this label
-
Performed numerical simulations using this paper's real time approach on the lattice of a given section before taking the continuum limit.
Aspects of quantum information in finite density field theory (2020) · 1.00 · paraphrased: the licence does not allow quoting
-
Evaluated the mutual information using this paper's correlator matrix method generalized to mixed density matrices.
Signatures of quantum phases in a dissipative system (2023) · 1.00 · paraphrased: the licence does not allow quoting
-
Calculated the oscillator part of the reduced density matrix using methods of free theories originally introduced by this paper.
Momentum-Space Entanglement in Heisenberg Spin-Half Ladders (2014) · 1.00 · paraphrased: the licence does not allow quoting
-
Calculated the oscillator part using this paper's method.
Entanglement spectra between coupled Tomonaga-Luttinger liquids: Applications to ladder systems and topological phases (2013) · 1.00 · paraphrased: the licence does not allow quoting
-
Calculated the entanglement entropy of the excited eigenstates using this paper's formula, which is valid also for excited eigenstates.
Bipartite entanglement entropy of the excited states of the free fermions and the harmonic oscillators (2019) · 1.00 · paraphrased: the licence does not allow quoting
-
Adopted this paper's method within the Hartree-Fock approximation to compute the correlation matrix between sites restricted within a domain.
Mutual information and correlations across topological phase transitions in topologically ordered graphene zigzag nanoribbons (2023) · 1.00 · paraphrased: the licence does not allow quoting
-
Followed this paper's formalism to encode the reduced density matrix of the Gaussian subsystem in the block covariance matrix.
Observable signatures of exceptional points from left-right eigenstate distinction (2026) · 1.00 · paraphrased: the licence does not allow quoting
-
Used this paper's formalism to construct the reduced density matrix in confirming numerically that the entropy scales with the square root of time for a noisy 1D hopping model.
Quantum Entanglement Growth Under Random Unitary Dynamics (2016) · 1.00 · paraphrased: the licence does not allow quoting
-
Used this paper's correlation matrix method for Gaussian states to study quantum quenches from a massive theory to a scale-invariant theory.
On the Time Scaling of Entanglement in Integrable Scale-Invariant Theories (2021) · 1.00 · paraphrased: the licence does not allow quoting
-
Used this paper's formula to give the entanglement Hamiltonian matrix from the correlation matrix.
Hyperfine Structure of Quantum Entanglement (2023) · 0.94 · paraphrased: the licence does not allow quoting
-
Denoting the eigenvalues of this $4ℓ×4ℓ$ matrix by $γ_α(t), α=1,…, 4ℓ$, the entanglement entropy is given by [ref]
Quench dynamics of entanglement from crosscap states (2024) · Free fermion model · 0.97
-
This key property enables us to write down the entanglement Hamiltonian $K_Ω$ of the ground state $|ψ⟩$ for an arbitrary region $Ω$ [ref] [equation] where $P_Ω=ΩP Ω$ is the spatial restriction of the spectral projector $P$.
Generalized Real-space Chern Number Formula and Entanglement Hamiltonian (2022) · Electric Hall conductance · 0.86
-
Wrote the entanglement Hamiltonian of the free fermion state in a subregion according to this paper, in terms of the two-particle correlation matrix.
Conserved Quantities from Entanglement Hamiltonian (2021) · 0.97 · paraphrased: the licence does not allow quoting
-
The reduced density matrix can be written as [ref] : [equation] where the $[ γ γ̅ ]$ contains all the Majorana fermions of subsystem $1$.
Entanglement entropy in quantum spin chains with broken parity number symmetry (2021) · Vacuum state · 0.98
-
The many-body overlap between the initial and evolved Slater determinants is given by [ref] [equation] If ${λ_k}$ are the eigenvalues of $O_m$, this becomes [equation]
Krylov Complexity in Periodically Driven CFTs and Critical Fermions (2026) · K-complexity of the correlation matrix. · 0.95
-
If $|GS⟩$ is the ground state, the classification is based on the projector [ref] [equation] where $k$ is the $d$ dimensional momentum in the translationally invariant directions, $i,j$ label the further spatial coordinates, and $α,β$ the unit cell degrees of freedom.
Topological Classification of One-Dimensional Chiral Symmetric Interfaces (2024) · Introduction · 0.85
How the literature uses it
110 arXiv papers cite this in 176 sentences
Papers citing it · 62 sentences without a confident label
Uses its method or tool 68 papers
-
The EE logarithmic coefficient for this model is finally calculated numerically in a radial lattice, using the method described in [ref] where the entropy is written in terms of vacuum two points correlators.
Numerical determination of the entanglement entropy for a Maxwell field in the cylinder · Introduction · 0.99
-
In this paper we solve this problem by exploiting the relation between ([ref]) and ([ref]) with the expressions obtained by Chung and Peschel [ref] for $ρ_V$ in terms of correlators for free boson and fermion discrete systems.
Analytic results on the geometric entropy for free fields · Introduction · 0.85
-
According to it, if the reduced density matrix satisfies the Wick decomposition property, i.e. the $n$-point functions factor out into two-point functions (see [ref] ), the full reduced density matrix $ρ_X$, of dimension $2^|X|$, can be obtained from the two-point correlation matrix, whose dimension is $|X|$.
Entanglement of several blocks in fermionic chains · Entanglement entropy and Toeplitz determinants · 0.88
Cites it as background 18 papers
-
Although the Gibbs state is Gaussian and determined by a finite-temperature Majorana correlation matrix [ref] , the stabilizer entropy is not an ordinary Gaussian thermodynamic observable.
Hidden Conformal Boundary Data in Finite-Temperature Stabilizer Entropy · Introduction · 0.85
-
As it is well known [ref] , due to the fact that the $n$-point vacuum expectation value satisfies the Wick decomposition property, it is possible to derive the reduced density matrix from the two point correlation function.
On the M\"obius transformation in the entanglement entropy of fermionic chains · Finite range Hamiltonians, compact Riemann surfaces and entanglement entropy · 0.82
-
To our knowledge, this goes back to [ref] for generic subregions in fermionic systems, to [ref] for a single block of spins, and to [ref] for generic subregions in bosonic systems.
Aspects of Pseudo Entropy in Field Theories · Calculation of Pseudo Entropy in Free Theories · 0.86
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.