Skip to content

Shape Deformations and Modular Perturbation Theory

A shape deformation changes the region, its causal domain, and generally its regulator boundary. Modular generators for two different regions cannot be differentiated until their algebras and domains have been identified. Once that choice is explicit, horizon stress fluxes, pullback terms, displacement operators, and contacts organize the response.

Required background. Ward-identity contact terms supply the distributional corrections, modular Hamiltonian domains supply the unbounded-generator conventions, and entanglement shape dependence supplies geometric variations.

Helpful background. Nonlocal modular generators explain why perturbative locality need not persist.

The overview’s differentiable-family map makes clear why a moving region first needs a common algebra identification, while its independent validity gates apply separately to operator domains, contacts, and regulator motion. The preceding chapter’s governed discussion of smooth modes, cusps, and fixed cutoffs supplies the relevant geometric figure; duplicating it here would not add a new relationship.

Let the reference entangling surface have embedding Xμ(y)X^\mu(\mathbf y) and deform it by

Xμ(y)⟼Xμ(y)+ϵζμ(y).X^\mu(\mathbf y) \longmapsto X^\mu(\mathbf y)+\epsilon\zeta^\mu(\mathbf y).

A tangential ζμ\zeta^\mu is locally a reparametrization. Its normal part changes the surface, the spatial region AϵA_\epsilon, and the causal domain D(Aϵ)D(A_\epsilon). Choose a diffeomorphism or unitary pullback

αϵ:A(Aϵ)⟶A(A0).\alpha_\epsilon: \mathcal A(A_\epsilon)\longrightarrow\mathcal A(A_0).

Only αϵ(KAϵ)\alpha_\epsilon(K_{A_\epsilon}) is a family on one representation. Changing the pullback shifts a commutator-like connection term, although invariant matrix elements agree when operator motion and contacts are transformed consistently.

For a smooth deformation of the vacuum half-space, Faulkner, Leigh, Parrikar, and Wang 2016, §1.1, eqs. (9)–(14) find, in their inward-pointing and null-coordinate convention,

αϵ(KA)=KA0−2πϵ∫H+ζ+T+++2πϵ∫H−ζ−T−−+ϵδKpullback+O(ϵ2).\alpha_\epsilon(K_A) =K_{A_0} -2\pi\epsilon\int_{\mathcal H_+}\zeta^+T_{++} +2\pi\epsilon\int_{\mathcal H_-}\zeta^-T_{--} +\epsilon\delta K_{\rm pullback} +O(\epsilon^2).

The source’s eq. (14) gives δKpullback\delta K_{\rm pullback} explicitly in its continuation convention. It must not be dropped or rewritten with a guessed factor of ii. Reversing the null normals, taking an outward deformation, or using inverse modular flow changes signs. The formula is an equality of regulated operators or controlled matrix elements; it is not an operator-norm expansion of arbitrary continuum modular Hamiltonians.

An exact CFT-vacuum family makes the first application particularly clean. Use the global metric convention (+−−−)(+---) and let

xμ=λξμ+yμ,ξμ=(1,1,0,…,0),x^\mu=\lambda\xi^\mu+y^\mu, \qquad \xi^\mu=(1,1,0,\ldots,0),

parametrize a null plane. The affine coordinate is λ\lambda and y∈Rd−2\mathbf y\in\mathbb R^{d-2} labels its generators. Let the future horizon start at λ=γ(y)\lambda=\gamma(\mathbf y). With Tλλ=TμνξμξνT_{\lambda\lambda}=T_{\mu\nu}\xi^\mu\xi^\nu, the vacuum-subtracted modular Hamiltonian on that horizon is

Hγ=2π∫dd−2y∫γ(y)∞dλ [λ−γ(y)]Tλλ(λ,y).H_\gamma =2\pi\int d^{d-2}y \int_{\gamma(\mathbf y)}^\infty d\lambda\, [\lambda-\gamma(\mathbf y)] T_{\lambda\lambda}(\lambda,\mathbf y).

Casini, Testé, and Torroba 2017, §1, eqs. (1.5)–(1.7), and §3.2 establish this local null-horizon expression and explain its normalization. Away from the null surface the generator need not be local.

Choose the compact smooth profile

γϵ(y)=ϵfL(y),\gamma_\epsilon(\mathbf y)=\epsilon f_L(\mathbf y),

with

fL(y)={exp⁡ ⁣(1−11−∣y∣2/L2),∣y∣<L,0,∣y∣≥L.f_L(\mathbf y)= \begin{cases} \displaystyle \exp\!\left(1-\frac{1}{1-\lvert\mathbf y\rvert^2/L^2}\right), &\lvert\mathbf y\rvert<L,\\[6pt] 0,&\lvert\mathbf y\rvert\geq L. \end{cases}

Here L>0L>0, ϵ\epsilon has dimensions of length, and max⁡fL=1\max f_L=1. Differentiating the exact formula gives the leading correction

δH=−2πϵ∫dd−2y fL(y)∫0∞dλ Tλλ(λ,y).\delta H =-2\pi\epsilon \int d^{d-2}y\,f_L(\mathbf y) \int_0^\infty d\lambda\, T_{\lambda\lambda}(\lambda,\mathbf y).

This is the requested compact deformation of a vacuum half-space. It is a transverse smearing of a half-sided null-energy operator, not a pointwise stress tensor.

The approximation has a reproducible matrix-element error. Taylor’s theorem applied on each null generator gives the exact remainder

R2=2π∫dd−2y∫0ϵfL(y)du [ϵfL(y)−u]Tλλ(u,y).R_2 =2\pi\int d^{d-2}y \int_0^{\epsilon f_L(\mathbf y)}du\, [\epsilon f_L(\mathbf y)-u] T_{\lambda\lambda}(u,\mathbf y).

At a fixed UV regulator, on a common finite-energy core, define

MΨ,Φ(y)=sup⁡∣u∣≤∣ϵ∣∣⟨Ψ∣Tλλ(u,y)∣Φ⟩∣.M_{\Psi,\Phi}(\mathbf y) =\sup_{\lvert u\rvert\leq\lvert\epsilon\rvert} \left\lvert \langle\Psi|T_{\lambda\lambda}(u,\mathbf y)|\Phi\rangle \right\rvert.

Then

∣⟨Ψ∣R2∣Φ⟩∣≤πϵ2∫dd−2y fL(y)2MΨ,Φ(y).\lvert\langle\Psi|R_2|\Phi\rangle\rvert \leq \pi\epsilon^2 \int d^{d-2}y\,f_L(\mathbf y)^2 M_{\Psi,\Phi}(\mathbf y).

In the regulated theory the shape-control parameter is

η=max⁡∣∇γϵ∣≤2.171∣ϵ∣L.\eta=\max\lvert\nabla\gamma_\epsilon\rvert \leq2.171\frac{\lvert\epsilon\rvert}{L}.

For example, ∣ϵ∣/L=10−2\lvert\epsilon\rvert/L=10^{-2} gives η<2.18×10−2\eta<2.18\times10^{-2}. The error ceiling is a statement about declared regulated matrix elements, not an operator norm. In the continuum, TλλT_{\lambda\lambda} is an operator-valued distribution, so the corresponding bound must be stated in a test-function or energy norm rather than with the pointwise MΨ,ΦM_{\Psi,\Phi}. Longitudinal light-ray integrals may be defined first with a smooth large-λ\lambda cutoff and then removed on the common domain.

In a replica description the twist operator is a codimension-two defect. Fix the displacement-operator sign by

∂μTtotμa(x)=−δΣ(x)Da(y)+derivative contacts,\partial_\mu T_{\rm tot}^{\mu a}(x) =-\delta_\Sigma(x)D^a(\mathbf y)+\text{derivative contacts},

where aa is normal to the defect. Billò, Gonçalves, Lauria, and Meineri 2016, §5.1, eqs. (5.21)–(5.24) derive the defect Ward identities and show how the displacement insertion represents shape response.

At second order, a separated-point correlator ⟨DaDb⟩\langle D^aD^b\rangle is not the whole Hessian. Varying the defect, surface measure, normal frame, stress tensor, and counterterms creates coincident contacts. Replica and horizon calculations agree only after matching normal orientation, analytic continuation in replica number, and the local counterterm scheme.

Cusp. Replace the bump by

γα(y)=α∣y1∣χ(y/L),\gamma_\alpha(\mathbf y) =\alpha\lvert y_1\rvert\chi(\mathbf y/L),

where α\alpha is a dimensionless cusp amplitude and χ\chi is smooth, compactly supported, and equals one near the origin. This α\alpha is distinct from the length-valued smooth-bump amplitude used above. The profile is continuous but not C1C^1, and distributionally

∂y12γα=2αδ(y1)+regular terms near the cusp.\partial_{y_1}^2\gamma_\alpha =2\alpha\delta(y_1)+\text{regular terms near the cusp}.

A smooth shape expansion that keeps only separated points misses this localized curvature/contact contribution. The general smooth-pullback remainder is therefore withdrawn. The stronger exact null-plane formula can still survive for a continuous cut on the light-ray domain; this does not promote the cusp to a smooth tangent.

Fixed-coordinate cutoff. Put a short-distance affine cutoff a>0a>0 on each null generator. A cutoff co-moving with the cut gives

Fco(γ)=∫γ+a∞dλ (λ−γ)T(λ),F_{\rm co}(\gamma) =\int_{\gamma+a}^\infty d\lambda\,(\lambda-\gamma)T(\lambda),

whereas a coordinate-fixed cutoff gives

Ffixed(γ)=∫a∞dλ (λ−γ)T(λ).F_{\rm fixed}(\gamma) =\int_a^\infty d\lambda\,(\lambda-\gamma)T(\lambda).

At γ=0\gamma=0,

Ffixed′(0)−Fco′(0)=aT(a).F_{\rm fixed}'(0)-F_{\rm co}'(0)=aT(a).

After transverse smearing, the apparent first-order modular response differs by 2πa∫fLT(a)2\pi a\int f_LT(a). This term cannot be discarded before controlling the near-cut stress matrix element and its counterterms. Thus a statement made with a fixed coordinate cutoff does not automatically describe a physically co-moving regulator.

Differentiating generators on different algebras. Supply the pullback and common quadratic-form domain first.

Keeping only one horizon flux. Future, past, and pullback terms are tied to the chosen Cauchy surface and identification.

Calling a cusp a small smooth deformation. Small amplitude does not bound derivatives. Smoothness and amplitude are independent control parameters.

  1. Show that a purely tangential deformation of a closed entangling surface is a reparametrization to first order.
Solution

Write ζμ=vi∂iXμ\zeta^\mu=v^i\partial_iX^\mu. Under the coordinate change yi↦yi+ϵviy^i\mapsto y^i+\epsilon v^i, the embedding becomes

Xμ(y+ϵv)=Xμ(y)+ϵvi∂iXμ+O(ϵ2).X^\mu(\mathbf y+\epsilon\mathbf v) =X^\mu(\mathbf y) +\epsilon v^i\partial_iX^\mu+O(\epsilon^2).

Thus the displaced points describe the same geometric surface to first order. A nonzero response from this direction alone signals an unmatched pullback, boundary, anomaly, or regulator term.

  1. Differentiate the exact null-cut modular Hamiltonian twice with respect to a one-parameter cut γs=sf\gamma_s=sf.
Solution

For one generator define

G(a)=∫a∞dλ (λ−a)T(λ).G(a)=\int_a^\infty d\lambda\,(\lambda-a)T(\lambda).

Leibniz differentiation gives

G′(a)=−∫a∞T(λ)dλ,G′′(a)=T(a).G'(a)=-\int_a^\infty T(\lambda)d\lambda, \qquad G''(a)=T(a).

Therefore

dHsfds∣0=−2π∫dd−2y f∫0∞dλ Tλλ,\left.\frac{dH_{sf}}{ds}\right|_0 =-2\pi\int d^{d-2}y\,f \int_0^\infty d\lambda\,T_{\lambda\lambda},

and

d2Hsfds2∣0=2π∫dd−2y f2Tλλ(0,y).\left.\frac{d^2H_{sf}}{ds^2}\right|_0 =2\pi\int d^{d-2}y\,f^2 T_{\lambda\lambda}(0,\mathbf y).

The second expression is distributional and requires the stated smearing and common domain.

  1. Derive the boundary mismatch between a co-moving affine cutoff and a coordinate-fixed cutoff.
Solution

For the co-moving prescription,

Fco′(0)=−aT(a)−∫a∞T(λ)dλ.F_{\rm co}'(0) =-aT(a)-\int_a^\infty T(\lambda)d\lambda.

For the fixed prescription only the weight changes, so

Ffixed′(0)=−∫a∞T(λ)dλ.F_{\rm fixed}'(0) =-\int_a^\infty T(\lambda)d\lambda.

Subtracting gives aT(a)aT(a). Multiplication by the modular normalization and transverse profile yields 2πa∫dd−2y f(y)T(a,y)2\pi a\int d^{d-2}y\,f(\mathbf y)T(a,\mathbf y). Its limit is controlled only if the near-boundary stress insertion and counterterms are controlled.

  • Billò, Marco, Vasco Gonçalves, Edoardo Lauria, and Marco Meineri. “Defects in Conformal Field Theory.” Journal of High Energy Physics 2016, no. 4 (2016): 091. DOI; arXiv.
  • Casini, Horacio, Eduardo Testé, and Gonzalo Torroba. “Modular Hamiltonians on the Null Plane and the Markov Property of the Vacuum State.” Journal of Physics A: Mathematical and Theoretical 50, no. 36 (2017): 364001. DOI; arXiv.
  • Faulkner, Thomas, Robert G. Leigh, Onkar Parrikar, and Huajia Wang. “Modular Hamiltonians for Deformed Half-Spaces and the Averaged Null Energy Condition.” Journal of High Energy Physics 2016, no. 9 (2016): 038. DOI; arXiv.

Original QFT.org content:CC BY 4.0, unless an item supplies different terms. Third-party material retains its own terms.