Skip to content

Entanglement Entropy of Free Fields

Free scalar and fermion theories are the principal controlled benchmarks for subregion entropy. Their ground and thermal states are Gaussian, so correlation matrices, replica determinants, and—in conformal limits—symmetry methods can be compared on the same regulated problem. The comparison reveals which terms are universal and which depend on mass, state, geometry, boundary conditions, or cutoff.

Required background. Review Gaussian fields and sources, heat kernels and spectral determinants, and Gaussian correlation-matrix entropy. Helpful background. Replica branched geometries supplies an independent formulation.

For a free scalar on a spatial lattice,

H=12pTp+12qTKq,X=12K1/2,P=12K1/2H=\frac12 p^Tp+\frac12q^TKq, \qquad X=\frac12K^{-1/2}, \qquad P=\frac12K^{1/2}

in the vacuum. Restricting XX and PP to a set of sites gives the symplectic spectrum and entropy. A free, number-conserving fermion is similarly determined by the occupied-mode projector restricted to AA. These are exact statements about the discretized theory and provide small-system tests for replica or Monte Carlo implementations.

The interval result of Holzhey, Larsen, and Wilczek 1994, pp. 443–467 gives, for a two-dimensional CFT vacuum on the infinite line,

SA=c3logϵ+s1,S_A=\frac{c}{3}\log\frac{\ell}{\epsilon}+s_1,

where cc is the central charge and s1s_1 is cutoff dependent. For a free real massless scalar the zero mode must be specified before applying this expression in finite periodic volume. For a free Dirac fermion no scalar zero-mode ambiguity appears, although spin structure and boundary conditions still matter.

In higher dimensions a smooth entangling surface has a leading divergence proportional to its area. Srednicki 1993, pp. 666–669 provides the classic spherical-lattice scalar example. The coefficient of the power divergence depends on the regulator; subleading logarithms or finite combinations can be more robust under stated geometric and renormalization conditions.

The structural map places Entanglement Entropy of Free Fields on the route from a regulated subsystem to integer moments, spectral checks, analytic continuation, and a continuum claim.

A regulated subsystem yields integer density-matrix moments by spectral or replica routes, while the von Neumann limit additionally requires analytic and growth assumptions.

The spectral and replica routes must agree on matched integer moments at fixed regulator. Continuation from those moments to n=1n=1 is logically separate and must state its analytic domain, branch, growth conditions, and order of limits. Schematic.

A mass introduces a correlation length ξm1\xi\sim m^{-1}. For regions much smaller than ξ\xi, the answer approaches the critical form plus controlled mass corrections. For regions much larger than ξ\xi, correlations across a smooth boundary are localized within O(ξ)O(\xi), so the entropy is dominated by boundary terms. Thermal or excited states can add extensive contributions; one must not label those terms vacuum entanglement without subtracting or otherwise controlling the state dependence.

Corners and conical features generate additional logarithms. Physical boundaries introduce boundary-condition data. A curved entangling surface allows local intrinsic- and extrinsic-curvature invariants. Free-field coefficients are excellent benchmarks for these structures, but they are not automatically the coefficients of an interacting theory.

Choose a massive scalar interval or half-space and compute its entropy at several lattice spacings by two routes:

  1. restrict the vacuum covariance and diagonalize its symplectic spectrum;
  2. evaluate the corresponding integer replica determinant with identical boundary conditions and cutoff geometry.

At n=2n=2 compare the same Rényi entropy before any continuation. Then vary aa at fixed physical mass, region, and volume. Agreement at one spacing tests implementation; common scaling over several spacings tests the continuum interpretation. Keep the determinant normalization and zero-mode prescription visible so the two methods do not merely share the same hidden error.

Approach the massless periodic scalar in two ways: at fixed finite volume, and after taking the volume large. The constant mode has frequency mm, so its covariance becomes singular as m0m\to0. If two prescriptions yield different constants or additional logarithms, those pieces are infrared dependent. Do not absorb them into the ultraviolet area-law coefficient. Zero Modes, Boundaries, and Infrared Sensitivity treats this separation in detail.

Before exporting this calculation, use the validity map to check normalization, infrared data, spectral or continuation control, and matched continuum scaling independently.

A regulated entropy claim passes normalization and sewing, infrared control, spectral and continuation checks, and matched continuum scaling; each missing step causes a distinct failure.

Normalization and sewing establish the intended integer moment; zero-mode and boundary control establish the infrared state; spectral and continuation checks control n1n\to1; geometry matching and a scaling window establish the continuum target. Omitting any stage licenses only a weaker conclusion. Schematic.

  • Holzhey, Christoph, Finn Larsen, and Frank Wilczek. “Geometric and Renormalized Entropy in Conformal Field Theory.” Nuclear Physics B 424 (1994): 443–467. arXiv; DOI.
  • Srednicki, Mark. “Entropy and Area.” Physical Review Letters 71 (1993): 666–669. arXiv; DOI.
  • Peschel, Ingo. “Calculation of Reduced Density Matrices from Correlation Functions.” Journal of Physics A: Mathematical and General 36 (2003): L205–L208. arXiv; DOI.