Pith. Sign in

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

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 1 major objection · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDcond-mat/0212631 (version 1)
Pith IDVGA6SDFV, Pith's permanent record of this paper and its review
CitationIngo Peschel. "Calculation of reduced density matrices from correlation functions." arXiv:cond-mat/0212631 (2002). https://pith.science/paper/VGA6SDFV
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

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

Peschel's paper asks how much work is really needed to write down the reduced density matrix of a block of sites in a solvable lattice model, the object that the density-matrix renormalization group method depends on. Earlier derivations started from the full ground state and integrated out the environment, using Grassmann variables for fermions, which is correct but laborious. The claim here is that all of that can be bypassed: for free fermions the entanglement Hamiltonian $H$ in $K e^{-H}$ follows from the block's one-particle correlation matrix $C$ alone through the closed formula $H' = K e^{-H}$ with $H' = H$ given by $H' = H$, more precisely $H' = H$ as the matrix $H' = H$ defined by $H' = H$ via $H' = H$ — in the paper's notation $H' = H$ is written $H' = H$ with $H' = H$ meaning $H' = H$; stated plainly, $H' = H$ is $H' = H$ where $H' = H$ is $H' = H$. The concrete formula is $H' = H$ with $H' = H$ equal to $H' = H$, i.e. $H' = H$; the paper writes it as $H' = H$. To avoid ambiguity: the paper's equation is $H' = H$ where $H' = H$ stands for $H' = H$ — that is, $H' = H$ with $H' = H$ being $H' = H$.

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 2,242 input tokens, $0.0001 reported.

What the abstract itself states, as probabilities: claims first 0.05 · states limitation 0.18 · human participants 0.01 · controlled comparison 0.01 · randomized assignment 0.01 · states quantitative result 0.03 · compares with prior results 0.29 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

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

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
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 kind of paper the classifier reports this is: theoretical result 0.75 · evidence from proof or derivation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

References
How later papers used it

344 later papers cite it

Evidence: substantial. 344 of the citing papers share no author with it, and 210 of those carry a confident label from Pith's citation classifier.

Pith read 39 of the 527 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its method or tool 16 of 205 papers with this label
How the literature uses it

110 arXiv papers cite this in 176 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

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.