---
title: Entanglement Susceptibility in Quantum Systems
url: https://www.emergentmind.com/topics/entanglement-susceptibility
type: topic
---

# Entanglement Susceptibility in Quantum Systems

Searching arXiv for recent and foundational papers on entanglement susceptibility and susceptibility-based entanglement diagnostics.
Entanglement susceptibility denotes a family of response-based constructions in which entanglement, or an entanglement witness, is organized as the response to a perturbation, a boundary coupling, a geometric deformation, or a measurable susceptibility. In the current literature the phrase does not refer to a single canonical observable. It may denote the coefficient of the \(2\)-Rényi entropy under boundary stitching [1205.2507], the kernels governing the variation of entanglement entropy under deformations of the entangling surface [1810.07209], the relative-entropy information metric of a reduced density matrix [2412.02236], or the use of dynamical or magnetic susceptibility to reconstruct entanglement witnesses such as the quantum Fisher information (QFI) [1509.01739, 2503.12463]. This plurality of usage is central to the subject.

## 1. Conceptual scope and terminology

The literature uses “entanglement susceptibility” in several technically distinct senses. The common structure is a response principle: an entanglement quantity, or a rigorous witness of entanglement, is made into a differential or perturbative object.

| Usage in the literature | Core object | Representative setting |
|---|---|---|
| Boundary-coupling response | \(S_2(\lambda)=2\lambda^2\chi_E+O(\lambda^3)\) | Area laws and gapless violations |
| Shape-response theory | \(\chi^{(1)},\chi^{(2)},\dots\) under \(\zeta(r)\) deformations | Corners, cones, trihedral vertices |
| Reduced-state information geometry | \(\Sigma_{ij}=\frac12\mathrm{Tr}\!\left[\rho\,\partial_i(\ln\rho)\,\partial_j(\ln\rho)\right]\) | Finite-size criticality in spin chains |
| Response-based entanglement witness | \(F_Q=\frac{4}{\pi}\int_0^\infty d\omega\,\tanh(\omega/2T)\chi''(\omega,T)\) | Thermal many-body systems |
| Thermodynamic witness from magnetic data | \(EW=1-\frac{6k_BT\chi}{g^2\mu_B^2N}\) and related bounds | Spin chains, dimers, molecular magnets |

Several papers also delimit the term by negation. Studies of fidelity susceptibility, first derivatives of entanglement entropy, or susceptibility–entanglement mappings often state explicitly that they do not define a quantity called entanglement susceptibility, even though they analyze closely related response diagnostics [1502.05930, 1004.3674, 2003.08148]. This suggests that “entanglement susceptibility” is best treated as a family of response constructions rather than a uniquely standardized observable.

## 2. Boundary-coupling entanglement susceptibility and area laws

A foundational definition appears in the perturbative study of bipartite entanglement for local Hamiltonians [1205.2507]. One considers a bipartition \(A\cup B=\Lambda\) and a one-parameter Hamiltonian
\[
H(\lambda)=H_A+H_B+\lambda H_\partial,\qquad \lambda\in[0,1].
\]
At \(\lambda=0\), the ground state factorizes as \(|0_A\rangle|0_B\rangle\). Using the \(2\)-Rényi entropy,
\[
S_2=-\log \operatorname{Tr}\rho_A^2,
\]
the leading entanglement generated by turning on the boundary coupling is
\[
S_2(\lambda)=2\lambda^2\chi_E+O(\lambda^3),
\]
with
\[
\chi_E:= \sum_{p_A,p_B\ge 1} \left| \frac{\langle p_A p_B|H_\partial|0_A0_B\rangle} {E_{p_A}+E_{p_B}-E_0} \right|^2.
\]
Only processes that excite both subsystems contribute. In this sense, \(\chi_E\) measures the amplitude for the boundary Hamiltonian to create correlated virtual excitations across the cut [1205.2507].

The same work establishes the comparison
\[
\chi_E\le \chi_F,
\]
where \(\chi_F\) is the fidelity susceptibility of the stitched ground state. For local gapped Hamiltonians with spectral gap \(\Delta\), exponential clustering, and boundary term \(H_\partial=\sum_{j\in\partial A}h_j\), the resulting bound is
\[
\chi_E \le \frac{\xi^{d-1}}{\Delta^2}\max_{j\in\partial A}\|h_j\|^2\,|\partial A|.
\]
This is an area-law statement at the level of entanglement response. By contrast, for gapless free fermions with a Fermi surface the same framework yields
\[
\chi_E\sim L^{d-1}\log L,
\]
which reproduces the familiar logarithmic violation of the area law [1205.2507].

The formalism extends beyond second order. The \(2\)-Rényi entropy can be re-expressed as a ground-state fidelity in a doubled and twisted theory, which leads to an exact linked-cluster expansion in connected imaginary-time correlators of a twisted boundary perturbation. This identifies entanglement susceptibility as the first nontrivial term in a full response expansion of Rényi entanglement [1205.2507].

## 3. Geometric entanglement susceptibilities

A second meaning of entanglement susceptibility concerns the response of entanglement entropy to deformations of the entangling surface [1810.07209]. If a point \(r\) on \(\partial A\) is displaced along the local unit normal by \(\zeta(r)\),
\[
r \mapsto r+\zeta(r)\, n(r),
\]
then the entropy admits the expansion
\[
S(A+\delta A)= S(A) + \int_{\partial A}\!\! d^{d-1} r\,  \chi^{(1)}( r) \zeta( r) + \frac{1}{2}\int_{\partial A}\!\! d^{d-1} r \int_{\partial A}\!\! d^{d-1} r'\, \chi^{(2)}( r, r') \zeta( r)\zeta( r') +\cdots .
\]
The kernels \(\chi^{(\ell)}\) are the entanglement susceptibilities. Only normal deformations are physically meaningful; tangential deformations merely reparameterize the surface [1810.07209].

For pure states and a planar bipartition, odd-order susceptibilities vanish. For the von Neumann entropy, strong subadditivity implies
\[
\chi^{(2)}( r, r')\leq 0.
\]
For translation- and scale-invariant states with a planar cut, the leading nonlocal second-order susceptibility is fixed up to local/contact terms:
\[
\chi^{(2)}( r- r') = - \frac{C}{| r- r'|^{2d}}.
\]
In conformal field theories, the coefficient is tied to the stress-tensor two-point normalization,
\[
C= \frac{2\pi^2}{d+2}\, C_T.
\]
For Rényi entropies, the same form holds with \(C\to C_n\), but the positivity constraint need not survive because strong subadditivity is unavailable for \(n\neq1\) [1810.07209].

This shape-response theory yields universal singular contributions for non-smooth regions. In \(d=2\), a nearly smooth corner with opening angle \(\Omega\simeq\pi\) gives
\[
\delta S^{(2)} = - \frac{C}{12} \, (\Omega -\pi)^2\log(L/\epsilon) + \dotsb .
\]
In \(d=3\), a nearly flat cone gives
\[
\delta S^{(2)} = -\frac{\pi^2 C}{256} (\Omega-\pi)^2 \left(\log(L/\epsilon)\right)^2  +\dotsb .
\]
For a nearly smooth symmetric trihedral vertex,
\[
\delta S^{(2)} = C \left( \frac{9}{32\sqrt{3} } + \frac{\pi}{4} \right) \left(\frac{2\pi}{3}-\theta \right) \log(L/\epsilon)+\dotsb .
\]
Accordingly, smooth-shape susceptibilities control universal corner, cone, and vertex terms in geometric entanglement entropy [1810.07209].

## 4. Reduced-state information geometry and the susceptibility of entanglement entropy

A more recent construction defines the susceptibility of entanglement entropy through the information geometry of reduced density matrices [2412.02236]. Starting from the relative entropy
\[
S(\hat{\rho} || \hat{\sigma}) = \mathrm{Tr}\!\left[\hat{\rho}\left(\ln \hat{\rho}-\ln \hat{\sigma}\right)\right],
\]
one expands between nearby states \(\hat\rho_\lambda\) and \(\hat\rho_{\lambda+\delta\lambda}\):
\[
S(\hat{\rho}_{\lambda} || \hat{\rho}_{\lambda+\delta \lambda}) = \Sigma_{ij}(\vec{\lambda})\, d\lambda^i d\lambda^j + \mathcal{O}(\delta\lambda^3),
\]
with
\[
\Sigma_{ij} = \frac{1}{2}\,\mathrm{Tr}\!\left[\hat{\rho}\,\partial_i(\ln \hat{\rho})\,\partial_j(\ln \hat{\rho})\right].
\]
The diagonal component \(\Sigma_{ii}\) is then interpreted as a susceptibility with respect to the control parameter \(\lambda_i\) [2412.02236].

In the finite \(1\)D XY chain and the transverse-field Ising model (TFIM), the construction is implemented using the single-site reduced density matrix
\[
\hat{\rho}^{(1)}=\frac{1}{2}\left(\mathbb{I}_{2\times 2}+m_z\,\hat{\sigma}^z\right),
\]
which yields the explicit one-parameter formula
\[
\Sigma_N(\lambda) = \frac{1}{2}\, \frac{\big(\partial_\lambda m_z(\lambda,N)\big)^2}{1-m_z^2(\lambda,N)}, \qquad \lambda\in\{h,\gamma\}.
\]
This object is not simply \(\partial_\lambda S_{\mathrm{EE}}\) or \(\partial_\lambda^2 S_{\mathrm{EE}}\); it is the Fisher–Rao-type metric of the reduced state [2412.02236].

Its utility is finite-size critical scaling. For the XY chain at \(h=0\), with control parameter \(\lambda=\gamma\), the turning points and global maxima of \(\Sigma_N^{\mathrm{XY}(\gamma)}\) converge to
\[
\gamma_c^\infty=0
\]
as
\[
|\gamma_c(N)-\gamma_c^\infty| = \frac{\pi}{2\sqrt{2}\,N} +\mathcal O\!\left(\frac{1}{N^2}\right).
\]
The corresponding maximum saturates numerically near
\[
\max \Sigma^{\mathrm{XY}(\gamma)} \sim 0.0897653.
\]
For the TFIM, with \(\lambda=h\) and \(\gamma=1\), the turning points converge to \(h_c^\infty=\pm1\) with
\[
|h_c(N)-h_c^\infty| \sim \frac{1}{N^{3/2}},
\]
while the maximum scales as
\[
\max\Sigma_N^{\mathrm{TFIM}(h)} \propto (\log N)^2.
\]
The analytical derivations use perturbation theory for the XY case and, for the TFIM, an exact route involving Chebyshev polynomials, hypergeometric functions, and complete elliptic integrals [2412.02236].

## 5. Dynamical susceptibility, quantum Fisher information, and multipartite entanglement depth

A separate and experimentally consequential line of work identifies entanglement through dynamical response functions. For a thermal many-body state and a Hermitian generator \(\hat O\), the QFI is exactly related to the imaginary part of the Kubo dynamical susceptibility [1509.01739]:
\[
F_Q(T) = \frac{4}{\pi} \int_{0}^{\infty} d\omega \, \tanh\!\left(\frac{\omega}{2T}\right)\chi''(\omega,T).
\]
This relation isolates the finite-frequency quantum response. It is not equivalent to the static isothermal susceptibility, which contains an elastic or Curie contribution and a classical Fisher-information piece,
\[
F_C=\frac{\chi_{\mathrm{el}}}{T}.
\]
Accordingly, a divergence of static thermal susceptibility at a thermal transition need not imply enhanced entanglement [1509.01739].

The QFI then functions as a multipartite entanglement witness. For \(N\) spin-\(\tfrac12\) particles and, when \(m\) divides \(N\), the density
\[
f_Q \equiv \frac{F_Q}{N},
\]
the criterion
\[
f_Q>m
\]
certifies at least \((m+1)\)-partite entanglement [1509.01739]. This framework carries over directly to correlated-electron systems. In the pseudogap regime of cuprate superconductors, the relevant observable is the antiferromagnetic susceptibility at the leading wave vector \(\mathbf Q\), and the QFI can be estimated either from neutron-scattering data or, more robustly in theory, through the quantum variance
\[
I_Q=8T\sum_{n=1}^{\infty}\chi(\mathbf Q,i\omega_n)\le F_Q.
\]
In that setting the entanglement criterion is stated as
\[
\text{If }F_Q>m,\text{ the system is at least }(m+1)\text{-partite entangled},
\]
and likewise \(I_Q>m\) suffices because \(I_Q\le F_Q\) [2503.12463].

The pseudogap calculations use the two-dimensional Hubbard model treated by DMFT followed by ladder D\(\Gamma\)A with \(\lambda\) correction. The susceptibility is fitted to an overdamped Ornstein–Zernike form,
\[
\chi(i\omega_n) = \frac{A}{\xi^{-2}+\gamma\left|i\omega_n\right|},
\]
which yields the low-temperature asymptotic QFI
\[
F_Q \approx \frac{4}{\pi }\frac{A}{\gamma}  \left[ 1 + \ln\left(\frac{W/2}{2T}\right)\right].
\]
Hence the normal-state entanglement susceptibility diverges logarithmically,
\[
F_Q\sim \ln(1/T),
\]
until superconductivity cuts off the growth [2503.12463]. At \(T=52\,\mathrm K\), the quantum variance in the theoretical pseudogap region exceeds the threshold for at least tripartite entanglement, \(I_Q>2\), and experimental neutron data in underdoped cuprates give \(F_Q>1\), certifying at least bipartite entanglement [2503.12463].

A related response-based use of QFI appears in quantum dots. There the transverse dynamical susceptibility is used to compute the QFI through the same integral formula, and the reported trends are nontrivial: in Ising dots far below the Stoner instability, lowering temperature can reduce the QFI, whereas anisotropic Heisenberg dots show enhanced QFI near isotropic points [2404.15477].

## 6. Thermodynamic witnesses, related diagnostics, and persistent distinctions

Long before the QFI–susceptibility relation was applied to correlated electrons, ordinary magnetic susceptibility was already used as a macroscopic witness of entanglement in low-dimensional spin systems. For the Heisenberg spin-\(\tfrac12\) chain \(\mathrm{Cu(tz)_2Cl_2}\), the witness is
\[
EW(N) = 1 - \frac{6 k_B T\, \chi}{g^2 \mu_B^2 N},
\]
with entanglement certified whenever \(EW(N)>0\). The fit to the Bonner–Fisher model gives \(J=5.6\,\mathrm K\), and the witness remains positive up to about \(12\,\mathrm K\) [1204.5669]. In copper nitrate, \(\mathrm{Cu(NO_3)_2 \times 2.5 H_2O}\), susceptibility and magnetization data are used in the same spirit across a field-driven quantum phase transition near \(B_c\approx 3.5\text{--}4\,\mathrm T\), with entanglement persisting up to about \(5\,\mathrm K\) [1109.1640].

The same logic extends to dimers and few-spin systems. In the engineered molecular magnet \([\mathrm{Fe}_{2}(\mu_{2}\text{-oxo})(\mathrm{C}_{3}\mathrm{H}_{4}\mathrm{N}_{2})_{6}(\mathrm{C}_{2}\mathrm{O}_{4})_{2}]\), the witness
\[
\mathcal{W} = \frac{3\,k_{B}\,T\,\bar{\chi}}{(g\mu_{B})^{2}N\,S} - 1
\]
is negative in the entangled regime, and extrapolation of the fitted dimer susceptibility gives an entanglement temperature \(T_e=732\,\mathrm K\), with analytic estimate \(761\,\mathrm K\) from \(J=-282\,\mathrm K\) and \(S=5/2\) [1208.5459]. For \(\mathrm{NH_4CuPO_4\cdot H_2O}\), another Heisenberg spin-\(\tfrac12\) dimer, the same witness gives \(T_c^{(\chi,\mathrm{EW})}=7.6\,\mathrm K\), while the entanglement of formation extracted from susceptibility survives to about \(16\,\mathrm K\) [1906.02769]. In a three-qubit XXX model with field and anisotropy, separable states satisfy
\[
\chi_z \ge \frac{1}{2k_B T},
\]
so entanglement is witnessed when
\[
\chi_z<\frac{1}{2T},
\]
and the critical temperature shifts as
\[
T_{crit} = 0.3K+ \frac{3}{2 \ln 3}
\]
for small anisotropy [2410.16133].

A different development is the empirical mapping between susceptibility and entanglement rather than a witness inequality. In the one-dimensional attractive Hubbard model, the single-site von Neumann entropy and magnetic susceptibility become linearly related in conventional superfluids, with
\[
S(n,\chi)=S_{\min}(n)+\frac{\chi-a}{b},\qquad a=0.032,\qquad b=-0.029,
\]
for \(P=0\) and \(U\lesssim -2t\) [2003.08148]. At finite temperature in the two-dimensional Hubbard model, the average single-site entanglement and magnetic susceptibility are also linearly connected in some regimes, while the derivatives \(\partial S/\partial U\) and \(\partial S/\partial n\) act as markers of interaction-driven quantum phase transitions [2505.00628].

Several conceptual distinctions remain essential. First, susceptibility-based witnesses are generally sufficient but not necessary: failure to violate the bound does not imply separability [1204.5669, 1109.1640, 2410.16133]. Second, QFI and quantum variance are witnesses or lower bounds on multipartite entanglement depth, not the full entanglement content of the many-body state [2503.12463]. Third, some closely related diagnostics are not called entanglement susceptibility in the underlying papers: the first derivative of entanglement entropy with respect to rung coupling in a two-leg \(XXZ\) ladder [1004.3674], the combination of fidelity susceptibility, entanglement entropy, and Schmidt gap in the spin-1 XXZ chain [1502.05930], and quench-time entanglement entropy studied alongside fidelity susceptibility and dynamical quantum phase transitions in multiband topological insulators [1911.02831]. The resulting picture is not terminological confusion so much as a stable division of labor: “entanglement susceptibility” now names several response-based constructions, each adapted to a different notion of entanglement and a different experimental or theoretical access route.

Source: https://www.emergentmind.com/topics/entanglement-susceptibility