Skip to content

Modular Berry Transport and Holonomy

Modular Berry transport compares an isolated modular spectral sector along a smooth family of states or regions. The comparison is not unique inside a degenerate sector: every change of frame that commutes with the modular generator is a zero-mode gauge transformation. Consequently, a closed-path matrix is meaningful only up to conjugation; its conjugacy class, eigenvalues, and class functions are the observables.

Required background. Modular spectra and spectral measures supply spectral projectors, direct-integral language, and the distinction between point and continuous spectrum.

Helpful background. State perturbations and susceptibility supply differentiable faithful families and their regulator controls.

The chapter’s differentiable-family map locates projected modular transport alongside state and shape response, while the independent validity gates separate the obligations a regulated holonomy must pass. In particular, a successful finite-dimensional calculation does not by itself settle the continuum spectral or algebra-identification gates.

Zero modes and constant-rank spectral blocks

Section titled “Zero modes and constant-rank spectral blocks”

Begin with a faithful finite regulator and a smooth Hermitian modular generator

K(λ)=U(λ)D(λ)U(λ)†.K(\lambda)=U(\lambda)D(\lambda)U(\lambda)^\dagger.

The relevant zero modes are operators in the kernel of the modular Liouvillian

LK(X):=[K,X].\mathcal L_K(X):=[K,X].

If K=∑ακαPαK=\sum_\alpha \kappa_\alpha P_\alpha, where the sum is over distinct eigenvalues, then

P0K(X)=∑αPαXPα.P_0^K(X)=\sum_\alpha P_\alpha X P_\alpha.

Thus a multiplicity-mαm_\alpha eigenspace contributes the full matrix algebra MmαM_{m_\alpha}, not merely mαm_\alpha independent phases. Equivalently,

P0K(X)=lim⁡T→∞12T∫−TTeisKXe−isK dsP_0^K(X)=\lim_{T\to\infty}\frac{1}{2T} \int_{-T}^{T}e^{isK}Xe^{-isK}\,ds

in the finite model, and in a continuum setting only when the mean ergodic limit exists in the stated topology. “Zero mode” here means zero Liouvillian frequency; it does not mean a zero eigenvalue of KK.

To transport one cluster, choose a spectral projector P(λ)P(\lambda) of fixed rank mm separated from the rest of the spectrum by a gap γ>0\gamma>0. Smoothness of the corresponding Riesz projector follows while that contour remains in the resolvent set Kato 1952, pp. 323–327. Internal eigenvalue coincidences are harmless if the whole cluster is retained. A crossing through the cluster boundary, or any closure of γ\gamma, removes the hypothesis needed for this connection.

Let Q(λ):Cm→HQ(\lambda):\mathbb C^m\to\mathcal H be an orthonormal frame for the chosen block, so Q†Q=1mQ^\dagger Q=\mathbf 1_m and QQ†=PQQ^\dagger=P. We use the anti-Hermitian frame connection

A:=Q†dQ,A†=−A,\mathcal A:=Q^\dagger dQ, \qquad \mathcal A^\dagger=-\mathcal A,

and orient parallel transport along increasing path parameter by

W[C]:=Pexp⁡ ⁣(−∮CA).W[C]:=\mathcal P\exp\!\left(-\oint_C\mathcal A\right).

For a complete diagonalizer, P0D(U†dU)P_0^D(U^\dagger dU) is the direct sum of these block connections. The superscript matters: U†dUU^\dagger dU is in the DD-frame, whereas P0KP_0^K acts in the physical KK-frame. This separates commuting frame motion from horizontal motion. The modular construction likewise defines zero modes by commutation with HmodH_{\rm mod} and treats their local frame freedom as gauge Czech et al. 2018, pp. 2–3, eqs. (4) and (12)–(13).

Under a change of block frame Q↦QgQ\mapsto Qg, g(λ)∈U(m)g(\lambda)\in U(m),

A↦g†Ag+g†dg,W[C]↦g(λ0)†W[C]g(λ0)\mathcal A\mapsto g^\dagger\mathcal A g+g^\dagger dg, \qquad W[C]\mapsto g(\lambda_0)^\dagger W[C]g(\lambda_0)

for a single-valued frame around a closed loop based at λ0\lambda_0. This is the Wilczek–Zee non-Abelian connection for a degenerate block Wilczek and Zee 1984, pp. 2111–2113. The curvature is

F=dA+A∧A=12Fij dλi∧dλj.\mathcal F=d\mathcal A+\mathcal A\wedge\mathcal A =\frac12\mathcal F_{ij}\,d\lambda^i\wedge d\lambda^j.

For a sufficiently small loop of linear size ℓ\ell in a smooth patch, define its antisymmetric area tensor by Σij=∫Σdλi∧dλj\Sigma^{ij}=\int_\Sigma d\lambda^i\wedge d\lambda^j. Then

W=1−12FijΣij+O(ℓ3).W=\mathbf1-\frac12\mathcal F_{ij}\Sigma^{ij}+O(\ell^3).

A Gaussian modular eigenspace around a loop

Section titled “A Gaussian modular eigenspace around a loop”

Here is a reproducible regulated QFT application. Take four fermionic wave-packet modes, grouped as a two-level orbital index with two identical flavors, and the faithful Gaussian state

ρ(n)=Z−1exp⁡ ⁣[−c†h(n)c],h(n)=ε n⋅σ⊗12,ε>0.\rho(\boldsymbol n)=Z^{-1} \exp\!\left[-c^\dagger h(\boldsymbol n)c\right], \qquad h(\boldsymbol n)=\varepsilon\,\boldsymbol n\cdot\boldsymbol\sigma\otimes\mathbf1_2, \qquad \varepsilon>0.

For a fermionic Gaussian state the one-particle correlation matrix is C=(1+eh)−1C=(\mathbf1+e^h)^{-1} Peschel 2003, pp. L205–L207, eqs. (3)–(7). The one-particle eigenvalues ±ε\pm\varepsilon are separated by 2ε2\varepsilon and each has multiplicity two. We transport the +ε+\varepsilon block along the latitude

n(ϕ)=(sin⁡θ0cos⁡ϕ,sin⁡θ0sin⁡ϕ,cos⁡θ0),0≤ϕ≤2π.\boldsymbol n(\phi) =(\sin\theta_0\cos\phi,\sin\theta_0\sin\phi,\cos\theta_0), \qquad 0\le\phi\le2\pi.

Choose

Q(ϕ)=(cos⁡(θ0/2)eiϕsin⁡(θ0/2))⊗12.Q(\phi)= \begin{pmatrix} \cos(\theta_0/2)\\ e^{i\phi}\sin(\theta_0/2) \end{pmatrix}\otimes\mathbf1_2.

The figure separates the spectral condition from the gauge statement. Inspect the nonzero 2ε2\varepsilon gap all along the latitude, then follow how changing the U(2)U(2) frame conjugates the closed-loop Wilson map without changing its spectrum.

A latitude loop carries an isolated rank-two plus-epsilon modular block separated by a gap two epsilon; a periodic eigenframe defines a horizontal frame and a Wilson map whose conjugacy data are invariant, while gap closure removes intrinsic spectral selection of the former block.

Along n(ϕ)\boldsymbol n(\phi) at θ0=π/3\theta_0=\pi/3, the +ε+\varepsilon one-particle block has rank two and remains separated from the −ε-\varepsilon block by 2ε>02\varepsilon>0. A periodic eigenframe QQ defines the horizontal frame Q~=QW\widetilde Q=QW, and a periodic frame change conjugates W[C]W[C], so only conjugacy data are frame independent. When ε→0\varepsilon\to0, the former rank-two projector is no longer intrinsically selected; transporting a chosen limiting subbundle would require extra structure. The diagram is schematic and not to scale.

Then

Aϕ=isin⁡2 ⁣(θ02)12,W[C]=e−iπ(1−cos⁡θ0)12.\mathcal A_\phi =i\sin^2\!\left(\frac{\theta_0}{2}\right)\mathbf1_2, \qquad W[C]=e^{-i\pi(1-\cos\theta_0)}\mathbf1_2.

At θ0=π/3\theta_0=\pi/3, the solid angle is Ω=π\Omega=\pi and W=−i12W=-i\mathbf1_2. This is a holonomy of an isolated one-particle modular block; the induced many-body Gaussian transport follows by second quantization, but it should not be confused with a finite-rank eigenspace of a continuum local-QFT modular operator.

A basis-independent discretization uses overlap matrices Mj=Q(ϕj)†Q(ϕj+1)M_j=Q(\phi_j)^\dagger Q(\phi_{j+1}) and their polar unitaries Lj=Mj(Mj†Mj)−1/2L_j=M_j(M_j^\dagger M_j)^{-1/2}. The links must be invertible, which holds for a sufficiently fine mesh on this smooth constant-rank loop. With ϕj=2πj/N\phi_j=2\pi j/N, define WN=LN−1†⋯L0†W_N=L_{N-1}^\dagger\cdots L_0^\dagger. For this loop,

Mj=[cos⁡2 ⁣θ02+e2πi/Nsin⁡2 ⁣θ02]12.M_j=\left[ \cos^2\!\frac{\theta_0}{2} +e^{2\pi i/N}\sin^2\!\frac{\theta_0}{2} \right]\mathbf1_2.

The convergence check at θ0=π/3\theta_0=\pi/3 is:

NN∥WN−(−i12)∥2\lVert W_N-(-i\mathbf1_2)\rVert_2
86.295×10−26.295\times10^{-2}
161.529×10−21.529\times10^{-2}
323.794×10−33.794\times10^{-3}
649.468×10−49.468\times10^{-4}

The controls are the spectral gap 2ε2\varepsilon, fixed rank two, loop orientation, polar-decomposition convention, and step number NN. The observed factor-of-four reduction under N↦2NN\mapsto2N is the expected O(N−2)O(N^{-2}) error for this smooth symmetric loop.

Replace the frame at every sampled point by QjgjQ_jg_j with arbitrary gj∈U(2)g_j\in U(2) and gN=g0g_N=g_0. Each link transforms as

Lj↦gj†Ljgj+1,L_j\mapsto g_j^\dagger L_jg_{j+1},

so the product telescopes:

WN↦g0†WNg0.W_N\mapsto g_0^\dagger W_Ng_0.

For a generic noncentral holonomy, individual entries change while the eigenvalues, tr⁡WN\operatorname{tr}W_N, and det⁡WN\det W_N do not. In the latitude benchmark above, WNW_N is proportional to 12\mathbf1_2, so conjugation leaves even its entries unchanged; the calculation verifies the covariance law but is not an entrywise gauge-sensitivity test. A result that changes by more than conjugation under the gjg_j test is basis contamination, not modular holonomy. If gN≠g0g_N\ne g_0, one has changed the endpoint identification; the missing transition function must be included before comparing closed-loop answers.

For a family of regions, the local algebras differ. One must first give embeddings or isomorphisms into a common reference algebra; changing those identifications changes the connection. In local QFT the modular spectrum is commonly continuous, so an eigenframe may be replaced by a direct integral, a spectral window, or an ergodic projection defined only in a weak topology. The inverse of ad⁡K\operatorname{ad}_K on the nonzero-frequency sector also requires a named domain.

A continuum claim must therefore state the algebra identification, spectral window or projector, topology of P0P_0, constant-rank or gap condition, common domain, path-discretization convergence, and regulator order. The special geometric interpretation in vacuum two-dimensional CFT and AdS3_3 is established under conformal and holographic hypotheses Czech et al. 2018, pp. 3–5; it is not a theorem that generic modular holonomy equals spacetime curvature.

Calling an eigenvalue of KK a zero mode. The zero-mode projector acts on operators and selects the kernel of [K,⋅][K,\cdot]. Degeneracy therefore enlarges the gauge algebra.

Mixing the physical and diagonal frames. Apply P0DP_0^D to U†dUU^\dagger dU, or conjugate both objects back to the KK-frame. Mixing them produces a frame-dependent formula.

Continuing through a closing gap. Internal degeneracy is allowed inside a retained block. A rank change or a collision with the excluded spectrum invalidates the same smooth bundle.

Let K=diag⁡(κ,κ,ν)K=\operatorname{diag}(\kappa,\kappa,\nu) with κ≠ν\kappa\ne\nu. Compute P0K(X)P_0^K(X) for a general 3×33\times3 matrix and identify the gauge group of the κ\kappa block.

Solution

The distinct-eigenvalue projectors are Pκ=diag⁡(1,1,0)P_\kappa=\operatorname{diag}(1,1,0) and Pν=diag⁡(0,0,1)P_\nu=\operatorname{diag}(0,0,1). Hence

P0K(X)=PκXPκ+PνXPν=(X11X120X21X22000X33).P_0^K(X)=P_\kappa XP_\kappa+P_\nu XP_\nu =\begin{pmatrix} X_{11}&X_{12}&0\\ X_{21}&X_{22}&0\\ 0&0&X_{33} \end{pmatrix}.

The full 2×22\times2 block survives because every operator in it commutes with κ12\kappa\mathbf1_2. Orthonormal frames in that block are related by U(2)U(2), not merely U(1)2U(1)^2. The discarded matrix elements have Liouvillian frequencies ±(κ−ν)\pm(\kappa-\nu).

2. Reproduce the Gaussian latitude holonomy

Section titled “2. Reproduce the Gaussian latitude holonomy”

Starting from the frame Q(ϕ)Q(\phi) above, calculate Aϕ\mathcal A_\phi and W[C]W[C]. Reverse the loop orientation and state what changes.

Solution

Differentiation gives

∂ϕQ=(0ieiϕsin⁡(θ0/2))⊗12,\partial_\phi Q= \begin{pmatrix}0\\ i e^{i\phi}\sin(\theta_0/2)\end{pmatrix} \otimes\mathbf1_2,

so Q†∂ϕQ=isin⁡2(θ0/2)12Q^\dagger\partial_\phi Q=i\sin^2(\theta_0/2)\mathbf1_2. Path ordering is trivial, and our sign convention gives

W=exp⁡ ⁣[−2πisin⁡2 ⁣θ02]12=e−iπ(1−cos⁡θ0)12.W=\exp\!\left[-2\pi i\sin^2\!\frac{\theta_0}{2}\right]\mathbf1_2 =e^{-i\pi(1-\cos\theta_0)}\mathbf1_2.

Reversing the path changes dϕ→−dϕd\phi\to-d\phi, hence W→W−1=W†W\to W^{-1}=W^\dagger. At θ0=π/3\theta_0=\pi/3, −i12-i\mathbf1_2 becomes +i12+i\mathbf1_2.

3. Test gauge covariance and a gap closure

Section titled “3. Test gauge covariance and a gap closure”

Show directly that the discrete polar link transforms as Lj↦gj†Ljgj+1L_j\mapsto g_j^\dagger L_jg_{j+1}. Then set ε→0\varepsilon\to0 in the Gaussian model and explain exactly which conclusion fails.

Solution

The overlap becomes Mj′=gj†Mjgj+1M_j'=g_j^\dagger M_jg_{j+1}. Therefore

Mj′†Mj′=gj+1†(Mj†Mj)gj+1,M_j'^\dagger M_j'=g_{j+1}^\dagger(M_j^\dagger M_j)g_{j+1},

and functional calculus gives (Mj′†Mj′)−1/2=gj+1†(Mj†Mj)−1/2gj+1(M_j'^\dagger M_j')^{-1/2}=g_{j+1}^\dagger(M_j^\dagger M_j)^{-1/2}g_{j+1}. Multiplication yields Lj′=gj†Ljgj+1L_j'=g_j^\dagger L_jg_{j+1}, so a closed product changes only by conjugation.

At ε=0\varepsilon=0, h=0h=0 has one four-dimensional eigenspace. The former +ε+\varepsilon rank-two projector is no longer selected by the spectrum: the gap 2ε2\varepsilon has closed and any rank-two subspace is an extra choice. The previously computed matrix remains the limit along the chosen family, but it is not an intrinsic holonomy of the degenerate endpoint without additional structure.

  • Czech, Bartłomiej, Lampros Lamprou, Samuel McCandlish, and James Sully. “Modular Berry Connection for Entangled Subregions in AdS/CFT.” Physical Review Letters 120 (2018): 091601. DOI; Open PDF.
  • Kato, Tosio. “On the Perturbation Theory of Closed Linear Operators.” Journal of the Mathematical Society of Japan 4 (1952): 323–337. DOI.
  • Peschel, Ingo. “Calculation of Reduced Density Matrices from Correlation Functions.” Journal of Physics A: Mathematical and General 36 (2003): L205–L208. DOI; Open PDF.
  • Wilczek, Frank, and A. Zee. “Appearance of Gauge Structure in Simple Dynamical Systems.” Physical Review Letters 52 (1984): 2111–2114. DOI.

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