---
title: Finite-Entanglement Scaling in Tensor Networks
url: https://www.emergentmind.com/topics/finite-entanglement-scaling
type: topic
---

# Finite-Entanglement Scaling in Tensor Networks

Searching arXiv for papers on finite-entanglement scaling and closely related developments.
Finite-entanglement scaling is the scaling framework in which the finite bond dimension $D$ or $\chi$ of a matrix product state (MPS), continuous MPS, or related tensor-network ansatz is treated as the control parameter that induces an effective infrared cutoff. In a critical system, finite bond dimension bounds the bipartite entanglement entropy, generates an emergent correlation length $\xi_D \sim D^\kappa$, and perturbs the state away from the conformal fixed point while preserving universal scaling relations in terms of the induced length scale. First established for infinite translationally invariant critical spin chains, the framework has since been used to extract central charges, scaling dimensions, and critical exponents, to distinguish finite-size and finite-entanglement regimes, and to formulate finite-correlation-length scaling in higher-dimensional tensor networks [0712.1976, 1204.3934, 1401.7654, 2411.01009, 2305.02091].

## 1. Foundational definition and scaling laws

The starting point is the incompatibility between critical entanglement growth and the finite Schmidt rank of an MPS. For an infinite, translationally invariant critical chain, a block of $L$ contiguous sites has entanglement entropy
$$
S(L)\simeq (c/3)\,\log L \qquad (L\gg 1),
$$
where $c$ is the central charge. By contrast, an MPS of bond dimension $\chi$ can never exceed
$$
S_{\max}=\log \chi.
$$
Calabrese–Cardy scaling away from criticality gives
$$
S(\xi)\simeq (c/6)\,\log \xi,
$$
with $\xi$ the correlation length. At the true critical point, a finite-$\chi$ MPS exhibits saturation of the half-chain entropy according to
$$
S_\chi \simeq \frac{c}{6}\,\log \xi_\chi,
$$
where $\xi_\chi$ is the emergent MPS correlation length. Numerical evidence then supports
$$
\xi_\chi \sim \chi^\kappa,
$$
and hence
$$
S_\chi \simeq \frac{c\kappa}{6}\,\log \chi.
$$
For the Ising universality class, $c=1/2$ and $\kappa=2$ give $S_\chi \simeq (1/6)\log \chi$; for the Heisenberg chain, the reported value is $\kappa\sim 1.37$ [0712.1976].

A complementary formulation uses $D$ rather than $\chi$ and emphasizes conformal invariance. At a conformally invariant critical point, one observes the universal scaling laws
$$
\xi_D \propto D^\kappa,\qquad
S_D \simeq \frac{c}{6}\,\log \xi_D + \mathrm{const}
       = \frac{c\kappa}{6}\,\log D + \mathrm{const},
$$
together with the analytic relation
$$
\kappa = \frac{6}{\,c\,\bigl(\sqrt{12/c}+1\bigr)}.
$$
This relation yields $\kappa\approx 2.03$ for the critical quantum Ising chain and $\kappa\approx 1.34$ for the Heisenberg chain [1102.4290, 1204.3934, 2306.08163].

| Model or universality class | Reported $c$ | Reported $\kappa$ |
|---|---:|---:|
| Quantum Ising | $1/2$ | $2$, $\approx 2.00\pm0.03$, or $\approx 2.03$ |
| Spin-$1/2$ Heisenberg | $1$ | $\sim 1.37$ or $\approx 1.34$ |
| ANNNI floating phase | $\approx 1.0$ | $\approx 1.34$ |

These laws motivate the term “finite-entanglement scaling”: the role played by physical size in finite-size scaling is taken over by the entanglement support permitted by the variational class [0712.1976, 1102.4290].

## 2. Effective correlation length and the transfer-matrix perspective

The induced length scale is defined from the transfer matrix of the tensor-network state. In an infinite MPS, the leading transfer-matrix eigenvalues determine a finite correlation length,
$$
\xi_\chi = \frac{1}{\ln(\lambda_0/\lambda_1)},
$$
or, in equivalent conventions,
$$
\xi_\chi = -\,\frac{1}{\log|\lambda_1/\lambda_0|}.
$$
For iPEPS, the analogous expression is
$$
\xi(D,\chi) = -\,1\Big/\ln\!\Bigl|\tfrac{\lambda_2}{\lambda_1}\Bigr|.
$$
In all of these settings, finite bond dimension produces a nonzero transfer-matrix gap and therefore a finite effective correlation length even at criticality [0712.1976, 2411.01009].

A major refinement is to use the induced correlation length itself as the scaling variable. In the notation of translation-invariant (c)MPS, one defines
$$
\mu_2(D) = -1/\mathrm{Re}\,\lambda_2,
$$
with $\lambda_2$ the transfer-matrix eigenvalue with largest nonzero real part. The finite-entanglement-scaling hypothesis is then written directly as
$$
S(D)\simeq \frac{c}{6}\log \mu_2(D)+\mathrm{const},
$$
and, for a two-point function of an operator of scaling dimension $\Delta$,
$$
G(x;D)\simeq \mu_2(D)^{-2\Delta}\,g\!\bigl(x/\mu_2(D)\bigr).
$$
This choice is motivated by several theoretical arguments and by the numerical observation that scaling with $\mu_2(D)$ produces significantly smaller error bars than scaling with $\log D$ [1401.7654].

The same transfer-matrix viewpoint underlies higher-dimensional generalizations. In iPEPS, one builds the row-to-row transfer matrix $T$, approximates its leading eigenvector by an iMPS of bond dimension $\chi$ using CTMRG, and extracts $\lambda_1,\lambda_2$ to obtain $\xi$. In this formulation, one tunes $\chi\gg D$ so that the effective correlation length is dominated by $D$ [2411.01009].

## 3. Empirical diagnostics and extraction of conformal data

The original finite-entanglement-scaling analysis assembled five independent tests of the scaling hypothesis. First, direct measurement of $\xi_\chi$ gives $\xi_\chi\propto \chi^\kappa$ with a single exponent $\kappa$. Second, at criticality the order parameter obeys
$$
M_\chi \sim \chi^{-\beta\kappa/\nu}.
$$
For the Ising model, with $\beta=1/8$ and $\nu=1$, fitting $M_\chi$ against $\chi$ yields $\kappa\approx 2.03\pm 0.02$. Third, the pseudo-critical point $\lambda_\chi^*$ shifts according to
$$
t^* \equiv |\lambda_\chi^*-\lambda_c|/\lambda_c \sim \chi^{-\kappa/\nu},
$$
giving $\kappa\approx 2.1\pm 0.1$ in the Ising case. Fourth, the entropy of a finite block satisfies $S(L)\sim (c/3)\log L$ for $L\ll \xi_\chi$ and saturates for $L\gg \xi_\chi$, with the saturation length obeying $L_{\rm sat}\sim \chi^\kappa$. Fifth, critical exponents can be extracted by data collapse, for example by plotting $\chi^{\beta\kappa/\nu}M$ versus $\chi^{\kappa/\nu}(\lambda-\lambda_c)$ and tuning $\kappa$ until the curves collapse onto a universal curve [0712.1976].

The algorithmic setting in which these tests were first carried out was the infinite time evolving block decimation algorithm. In the two-site unit-cell ansatz,
$$
|\Psi\rangle= \cdots \Gamma^A\Lambda\Gamma^B\Lambda\Gamma^A\Lambda\Gamma^B\cdots,
$$
the diagonal matrix $\Lambda$ contains the Schmidt coefficients, and the half-chain entropy is
$$
S_\chi = -\sum_i \Lambda_{ii}^2 \log \Lambda_{ii}^2.
$$
The iTEBD workflow consists of starting from a random translationally invariant MPS of bond dimension $\chi$, applying a Suzuki–Trotter decomposition of $e^{-\tau H}$ into two-site gates with small time step $\epsilon$, truncating after each gate to the $\chi$ largest singular values, and iterating $\tau\to\infty$ until energy, entropy, and local observables converge. The results were cross-checked against standard infinite- and finite-DMRG runs, with agreement to $\sim 10^{-8}$ in energy and $10^{-4}$ in entropy away from the thin critical window [0712.1976].

Subsequent work turned these observations into a systematic procedure for extracting conformal data. In translationally invariant (c)MPS and MPS simulations, one performs linear regression of $S_{\rm hi}(D)$ versus $\log\mu_2(D)$ to obtain $c/6$, and of $\log G_O(s\,\mu_2(D);D)$ versus $\log\mu_2(D)$ to obtain $-2\Delta_O$. Applied to the Lieb–Liniger model, the massless relativistic boson, and the critical quantum Ising chain, this yielded central charges within a few percent and scaling dimensions accurate to $\sim 0.1\%$ for primary fields [1401.7654].

## 4. Finite-entanglement versus finite-size scaling, and unified crossover theory

In a critical one-dimensional system with periodic boundary conditions, finite size and finite entanglement generate two distinct artificial gaps. The finite-size gap is
$$
\Delta E_N \sim \xi_N^{-1}\sim 2\pi x_1/N,
$$
where $x_1$ is the smallest scaling dimension of the CFT. The finite-entanglement gap is
$$
\Delta E_D \sim \xi_D^{-1}\propto D^{-\kappa}.
$$
This produces two regimes in the $(N,D)$ plane. Finite-size scaling (FSS) corresponds to
$$
\xi_D \gg \xi_N,
$$
so that the MPS is accurate enough to resolve the periodic-chain gap and standard CFT formulas such as
$$
E_0(N)-E_0(\infty)= -v_f\pi c/(6N^2)+O(N^{-4})
$$
apply. Finite-entanglement scaling (FES) is the complementary regime
$$
\xi_N \gg \xi_D \qquad \Leftrightarrow \qquad N\gg \xi_D\sim D^\kappa,
$$
in which the MPS is effectively blind to the chain’s finite size or boundary conditions. A sharp crossover is observed when $\xi_D\simeq N$, and only the FSS regime maintains a finite overlap with the exact periodic-boundary-condition ground state in the large-$N$ limit [1204.3934].

This distinction resolves a common misconception: finite-$D$ simulations can yield correct universal local quantities while failing to reproduce global finite-$N$ data. In the FES regime, local observables and local scaling laws remain meaningful, whereas boundary-condition-sensitive quantities require remaining in the FSS regime and extrapolating $D\to\infty$ afterwards [1204.3934].

A broader unification includes equilibrium criticality, finite bond dimension, finite system size, and linear driving. With $g=Rt$, the four competing length scales are the equilibrium correlation length $\xi\sim |g|^{-\nu}$, the finite-bond-dimension length $\xi_D\sim D^\kappa$, the system size $L$, and the finite-driving length $\xi_R\sim R^{-1/r}$ with $r\equiv z+1/\nu$. The scaling ansatz for the half-chain entanglement entropy is
$$
S(g,D^{-1},L^{-1},R) = (c/6)\,\ln b + (c/6)\,S\bigl(g\,b^{1/\nu},\,D^{-1}b^{1/\kappa},\,L^{-1}b,\,R\,b^r\bigr).
$$
Choosing $b$ so that one scaled variable is of order unity yields the leading forms for equilibrium scaling, FES, FSS, and finite-time scaling. In this unified picture, whichever of $\xi$, $\xi_D$, $L$, or $\xi_R$ is shortest controls the leading logarithmic growth of $S$, while the longer scales enter as analytic subleading corrections [1806.10290].

## 5. Generalizations beyond the original one-dimensional setting

Finite-entanglement scaling was quickly applied to models with more intricate phase structure. In the one-dimensional ANNNI model, iMPS finite-entanglement scaling was used to study second-order transitions and the supposedly infinite-order Kosterlitz–Thouless transition. The analysis used the correlation length, entanglement entropy, and the second derivative of the energy with respect to the external field, together with an off-critical ansatz
$$
\log\xi_D
=\kappa\log D-\log\bigl[D^\kappa/\xi_\infty + e^{-\alpha D^\kappa/\xi_\infty}\bigr]+\mathrm{const}.
$$
At $\kappa=-0.25$, the FM–PM transition gave $h_c\approx 0.536$, $\nu\approx 1.02$, and $c\approx 0.5$. At $\kappa=-0.75$ and $h=0.40$, the floating phase was found to be critical with $c\approx 1.0$ and measured correlation-length exponent $\kappa\approx 1.34$. At the floating–paramagnetic transition, $\xi_D\propto D^\kappa$ held up to $D\approx 80$, beyond which $\log\xi_D$ bent downward, and fitting the crossover ansatz yielded $h_c\approx 0.42\text{–}0.44$ together with an essential singularity in $\xi_\infty^{-1}$ [1102.4290].

In two-dimensional iPEPS, the framework reappears as finite-correlation-length scaling. For the disorder parameter
$$
\tilde m_g(R_x,R_y)=\Bigl\langle \prod_{i\in R}U_i(g)\Bigr\rangle,
$$
one writes, in the critical region,
$$
\tilde m_g \sim \exp\!\bigl[-\,\xi\,\alpha_{\rm peri}(y,\xi)-\xi^2\,\alpha_{\rm area}(y,\xi)\bigr].
$$
Dimensional analysis gives
$$
\alpha_{\rm peri}(y,\xi)=\alpha_{\rm peri}(y)+\xi^{-1}f_0\!\bigl(y\,\xi^{d+z-\Delta_y}\bigr),\qquad
\alpha_{\rm area}(y,\xi)=\alpha_{\rm area}(y)+\xi^{-2}f_1\!\bigl(y\,\xi^{d+z-\Delta_y}\bigr).
$$
For the $(2+1)$-dimensional transverse-field Ising model, the extrapolated values were $J_c\approx 0.32847$, $\nu=0.62997$, $\alpha_{\rm peri,c}=0.0382$, and $\alpha_{\rm area,c}\approx 0$, implying a pure perimeter law at criticality. Data collapse used the scaling variable
$$
x=(J-J_c)\,\xi^{1/\nu}
$$
and the rescaled observables $(\alpha_{\rm peri}-\alpha_{\rm peri,c})\xi$ and $\alpha_{\rm area}\xi^2$ [2411.01009].

A further extension addresses two-dimensional metals with a Fermi surface. For the optimal tensor-network state approximation to a metallic state, the entanglement entropy of a contractible region of linear size $L$ is proposed to scale as
$$
S_{TNS}(L,D)\simeq \alpha\,L\,\log[\,\xi\,f(L/\xi)\,],
$$
with
$$
\alpha = \frac{1}{24\pi}\oint_{\partial R}\oint_{FS}|dS_x\!\cdot dS_k|.
$$
The scaling function satisfies $f(x)\sim \Lambda x$ for $x\ll 1$ and $f(x)\to \Lambda^{-1}$ for $x\gg 1$, reproducing respectively the Widom-type $L\log L$ growth and area-law saturation. In the Gaussian-fermion PEPS study of a spinless circular Fermi surface at half filling, the reported scaling exponent was $\kappa\simeq 1.074$, and moderate bond dimensions $D\approx 16\ldots 64$ already produced correlation lengths of order $10\ldots 10^2$ lattice sites. The same work also emphasized generalized Lieb–Schultz–Mattis–Oshikawa–Hastings constraints at fractional filling and their implications for charge-conserving PEPS [2305.02091]. This suggests that finite-entanglement scaling is not confined to one-dimensional conformal critical points, although the universal objects being scaled differ substantially.

## 6. Effective-field-theory and boundary-CFT interpretations

A more recent line of work interprets finite bond dimension as an actual relevant deformation of the critical theory rather than merely a numerical cutoff. In one formulation,
$$
H_{\rm eff}=H_*+\sum_g c_g\,\xi^{\Delta_g-2}O_g+\cdots,
$$
where the most relevant operator allowed by the residual MPS symmetry controls the deformation. In another,
$$
H_{\rm eff}=H_{\rm CFT}+g\!\int\!dx\,\mathcal O_{\rm pert}(x),
$$
with $g_0\sim \xi(\chi)^{-(D-\Delta)}$ at the lattice scale. Under RG,
$$
\frac{dg(\ell)}{d\ell}=(D-\Delta)\,g(\ell),
$$
and the induced gap forms at a scale $\ell_*=\ln\xi$. This picture is explicitly compatible with the standard law $\xi(\chi)\sim \chi^\kappa$ and with the central-charge relation
$$
\kappa(c)=\frac{6}{c\bigl(1+\sqrt{12/c}\bigr)}.
$$
The finite-entanglement perturbation therefore has a field-theoretic interpretation without invalidating the original scaling hypothesis [2306.08163, 2411.03954].

An important consequence is a boundary-CFT interpretation of the entanglement Hamiltonian. The reduced density matrix can be written as
$$
\rho_A \propto e^{-2\pi H_E},
$$
and the spectrum of $H_E$ is that of a BCFT on a strip of width $W=\ln\xi$, with two boundaries: an entanglement boundary fixed by the bipartition and a physical boundary determined by the relevant deformation selected by the MPS symmetry. In the critical Ising chain, a generic MPS without $\mathbb Z_2$ symmetry activates the $\sigma$ field and yields the $\sigma$ tower, whereas a $\mathbb Z_2$-symmetric MPS forbids $\sigma$, activates $\varepsilon$, and yields the identity and $\varepsilon$ towers. Analogous constructions were demonstrated for the three-state Potts model and the free compact boson, including symmetry-engineered conformal boundary conditions and the recovery of the boson radius from the entanglement spectrum [2306.08163].

A related development is the proposal of a renormalization-group self-congruent point. The claim is that the running coupling at the induced infrared scale,
$$
g_*(\chi)=g_0(\chi)\,\xi(\chi)^{D-\Delta},
$$
approaches a $\chi$-independent constant. The transfer-matrix spectrum then becomes scale-invariant but differs from the unperturbed BCFT values. For the Ising model, the reported universal transfer-matrix gap ratios are
$$
\{R_1,R_2,R_3,\dots\}=\{0,\;0.1769,\;0.8130,\;1,\dots\},
$$
rather than the free–free BCFT values
$$
\{0,\;0.25,\;0.75,\;1,\dots\}.
$$
The paper argues that this does not alter the validity of finite-entanglement scaling because the self-congruent point lies within the scaling regime: bulk observables continue to scale with $\xi(\chi)$, while boundary-spectrum data are shifted in a universal, $\chi$-independent way [2411.03954].

Several caveats recur across the literature. Corrections of order $O(1/\log D)$ can bias estimates of $c$ and $\kappa$ if the accessible bond dimensions are too small. The off-critical ansatz used in the ANNNI analysis is phenomenological and could in principle be refined. Weakly gapped regions, especially near Kosterlitz–Thouless transitions, may require sufficiently large $D$ before deviations from critical scaling become visible. For finite systems with periodic boundary conditions, faithful extraction of global finite-$N$ quantities still requires remaining in the finite-size-scaling regime rather than the finite-entanglement regime. In cMPS-based conformal-data extraction, the extraction of OPE coefficients remains to be implemented [1102.4290, 1204.3934, 1401.7654].

Finite-entanglement scaling is therefore best understood as a family of scaling theories centered on an induced correlation length. In its original MPS form, it provides a thermodynamic-limit analogue of finite-size scaling. In its transfer-matrix formulation, it turns the emergent length scale into the primary observable. In its more recent effective-field-theory and BCFT formulations, it interprets finite bond dimension as a symmetry-constrained relevant perturbation with a nontrivial boundary spectrum. Across these formulations, the core statement remains the same: finite entanglement imposes a controlled infrared scale that can be used, rather than merely tolerated, to extract universal critical information [0712.1976, 1401.7654, 2306.08163, 2411.03954].

Source: https://www.emergentmind.com/topics/finite-entanglement-scaling