Skip to content

Relativistic Fields, Spinors, and the Ising Continuum Limit

The previous page produced a surprise: a two-dimensional Ising model written in terms of ordinary commuting spins contains a natural fermionic field. The field is not inserted by hand. It is the point-split product of an order operator and a disorder operator, and its fermionic sign comes from transporting the disorder cut past the order insertion.

This page translates that observation into continuum field theory. We first classify scalar, vector and spinor fields by their rotation laws. Then we explicitly project the square-lattice corner equation onto its two critical modes, obtaining a Majorana equation with m=4(K−Kc)/am=4(K-K_c)/a to leading scaling order in the conventions below. This limit assumes a∣k∣≪1a|k|\ll1 and ∣K−Kc∣≪1|K-K_c|\ll1 in the bulk, away from additional coincident insertions. The same calculation identifies the critical part of the Ising energy with the fermion mass operator and explains the logarithmic specific heat.

Required background. Order, disorder, and the Ising fermion supplies the four corner directions, transported cuts and exact local propagation relation.

Helpful background. Critical propagators and the upper critical dimension supplies the scaling interpretation of a correlation length; two-dimensional disorder lines supplies the underlying mixed-correlator sign proof. The rotation review below is self-contained.

A continuum field is not defined only by its equation of motion. It is also defined by how it transforms under spacetime symmetries. In Euclidean field theory, the local rotation group is SO(d)\mathrm{SO}(d), or more precisely its double cover Spin⁡(d)\operatorname{Spin}(d) when spinors are present. A field with components ΦA(x)\Phi_A(x) transforms as

ΦA′(x′)=D(R)ABΦB(x),x′=Rx,\Phi'_A(x')=D(R)_{AB}\Phi_B(x), \qquad x'=Rx,

where D(R)D(R) is a finite-dimensional representation carried by the field indices. A scalar has D(R)=1D(R)=1. A vector has D(R)=RD(R)=R. A spinor has a matrix S(R)S(R) satisfying

S(R)−1γiS(R)=Rijγj.S(R)^{-1}\gamma_i S(R)=R_{ij}\gamma_j.

This is the first difference between writing down a field and writing down a collection of components. A two-component object is not automatically a spinor; it is a spinor only if rotations act by the spin representation.

A free scalar field is the simplest example. In Euclidean signature one may write

S0[ϕ]=12∫ddx [(∂iϕ)2+m2ϕ2],S_0[\phi] ={1\over2}\int d^d x\, \left[(\partial_i\phi)^2+m^2\phi^2\right],

so the classical equation is

(−∂2+m2)ϕ=0.(-\partial^2+m^2)\phi=0.

After analytic continuation to Lorentzian signature this becomes the Klein–Gordon equation. The field carries no internal rotation index; all angular momentum of a scalar excitation is orbital.

A vector field carries one rotation index:

Ai′(x′)=RijAj(x).A_i'(x')=R_{ij}A_j(x).

A tensor field carries several such indices. For example,

Bij′(x′)=RikRjℓBkℓ(x).B_{ij}'(x')=R_{ik}R_{j\ell}B_{k\ell}(x).

These transformation laws are local. They do not yet say which field equations are correct or how many physical polarizations propagate. Those are dynamical questions.

A massive vector field can be described by the Euclidean Proca action

S[A]=∫ddx [14FijFij+m22AiAi],Fij=∂iAj−∂jAi.S[A]=\int d^d x\, \left[ {1\over4}F_{ij}F_{ij}+{m^2\over2}A_iA_i \right], \qquad F_{ij}=\partial_iA_j-\partial_jA_i.

Varying AiA_i gives

−∂jFji+m2Ai=0.-\partial_jF_{ji}+m^2A_i=0.

Expanding FjiF_{ji},

∂jFji=∂2Ai−∂i∂jAj,\partial_jF_{ji} =\partial^2A_i-\partial_i\partial_jA_j,

so the equation is

(−∂2+m2)Ai+∂i∂jAj=0.(-\partial^2+m^2)A_i+\partial_i\partial_jA_j=0.

Taking a divergence gives

m2∂iAi=0.m^2\partial_iA_i=0.

Therefore, for m≠0m\ne0,

∂iAi=0.\boxed{\partial_iA_i=0.}

The vector equation plus this constraint is equivalent to a Klein–Gordon equation for the transverse components:

(−∂2+m2)Ai=0,∂iAi=0.(-\partial^2+m^2)A_i=0, \qquad \partial_iA_i=0.

In momentum space, the decomposition is especially transparent:

Ai(k)=Ai⊥(k)+Ai∥(k),kiAi⊥(k)=0,A_i(k)=A_i^{\perp}(k)+A_i^{\parallel}(k), \qquad k_iA_i^{\perp}(k)=0,

with transverse projector for nonzero Euclidean momentum kk,

Pij⊥(k)=δij−kikjk2.P_{ij}^{\perp}(k)=\delta_{ij}-{k_i k_j\over k^2}.

Transverse and longitudinal decomposition of a vector field

A vector field decomposes into components transverse and longitudinal to the momentum. The massive Proca equation imposes ∂iAi=0\partial_iA_i=0 on shell, leaving the spin-one polarizations.

This is a useful warning for later gauge theory. A vector field has a vector index, but not every component of that index represents an independent physical polarization. Gauge redundancy, constraints, and equations of motion all matter. In the Ising problem the analogous lesson is that the four corner labels of the lattice fermion are not four independent particles. They are a convenient lattice package whose long-distance content is a smaller spinor field.

Spinors are not ordinary tensors. For an active counterclockwise rotation, fix

Rθ=(cos⁡θ−sin⁡θsin⁡θcos⁡θ),x′=Rθx.R_\theta=\begin{pmatrix}\cos\theta&-\sin\theta\\\sin\theta&\cos\theta\end{pmatrix}, \qquad x'=R_\theta x.

In the chosen gamma basis, the spin transformation is

Ψ=(uv),Ψ′(x′)=SθΨ(x),Sθ=exp⁡ ⁣(−θ2γ1γ2)=(e−iθ/200eiθ/2).\begin{gathered} \Psi=\begin{pmatrix}u\\v\end{pmatrix},\qquad \Psi'(x')=S_\theta\Psi(x),\\ S_\theta=\exp\!\left(-\frac\theta2\gamma_1\gamma_2\right) =\begin{pmatrix}e^{-i\theta/2}&0\\0&e^{i\theta/2}\end{pmatrix}. \end{gathered}

Indeed, γ1γ2=i diag(1,−1)\gamma_1\gamma_2=i\,\mathrm{diag}(1,-1) and direct multiplication gives Sθ−1γiSθ=(Rθ)ijγjS_\theta^{-1}\gamma_iS_\theta=(R_\theta)_{ij}\gamma_j. Together with ∂i′=(Rθ)ij∂j\partial'_i=(R_\theta)_{ij}\partial_j, this makes the Dirac equation covariant. The same rotor is derived in the one-plane Spin example.

Thus u′(x′)=e−iθ/2u(x)u'(x')=e^{-i\theta/2}u(x) and v′(x′)=eiθ/2v(x)v'(x')=e^{i\theta/2}v(x). These signs depend on the complete coordinate-and-field convention; checking only a full turn would not detect their reversal. In particular,

S2π(u)=−u,S2π(v)=−v.S_{2\pi}(u)=-u, \qquad S_{2\pi}(v)=-v.

A spinor changes sign under a full 2π2\pi rotation. Only after a 4π4\pi rotation does it return to itself.

A counterclockwise active rotation gives opposite half-angle phases to u and v, with a minus sign after a full turn

For x′=Rθxx'=R_\theta x and γ1=σ1\gamma_1=\sigma_1, γ2=σ2\gamma_2=\sigma_2, the upper component acquires e−iθ/2e^{-i\theta/2} and the lower component eiθ/2e^{i\theta/2}. Both acquire −1-1 under a full turn. The coordinate rotation is schematic; its orientation and the component phases use the same convention.

The Euclidean two-dimensional Dirac equation is a first-order equation,

(γi∂i−m)Ψ=0.(\gamma_i\partial_i-m)\Psi=0.

With the gamma matrices stated above and Ψ=(u,v)T\Psi=(u,v)^T, this is equivalent to

(∂1+i∂2)u=mv,(∂1−i∂2)v=mu.\boxed{ (\partial_1+i\partial_2)u=m v, \qquad (\partial_1-i\partial_2)v=m u. }

At m=0m=0 these equations split into holomorphic and antiholomorphic pieces:

∂ˉu=0,∂v=0.\bar\partial u=0, \qquad \partial v=0.

Thus one component depends only on zz, and the other depends only on zˉ\bar z. This is the continuum seed of the Ising conformal field theory: the critical fermion separates into left-moving and right-moving Majorana components.

There is no contradiction between the first-order Dirac equation and the second-order Klein–Gordon equation. Squaring the operator gives

(γi∂i+m)(γj∂j−m)Ψ=(∂2−m2)Ψ,(\gamma_i\partial_i+m)(\gamma_j\partial_j-m)\Psi =(\partial^2-m^2)\Psi,

With these fixed conventions, this identity is exact. Multiplying by an overall minus sign gives the positive Euclidean operator −∂2+m2-\partial^2+m^2. Each component obeys the associated second-order equation, while the first-order equation also relates the components.

Recall the point-split order–disorder composite from the previous page. Put a spin at an original-lattice site xx and a disorder endpoint at a neighboring dual site x+eax+e_a, where the four half-lattice vectors have angles

θa=π4+π2(a−1),a=1,2,3,4.\theta_a={\pi\over4}+{\pi\over2}(a-1), \qquad a=1,2,3,4.

The four corner fields are

χa(x)=σxμx+ea.\chi_a(x)=\sigma_x\mu_{x+e_a}.

Because moving the disorder endpoint once around the spin crosses the branch cut once, the corner field is antiperiodic:

χa+4(x)=−χa(x).\chi_{a+4}(x)=-\chi_a(x).

This antiperiodicity is exactly the lattice form of spinorial behavior. A scalar function of a continuous direction would have integer angular harmonics, whereas an antiperiodic function has half-integer harmonics. Here, however, the direction is sampled at only four corners. The corner-label space is therefore four-dimensional, with four independent Fourier characters. A convenient set of representatives is

s∈{+12,−12,+32,−32}(s mod 4).s\in\left\{+{1\over2},-{1\over2},+{3\over2},-{3\over2}\right\} \qquad (s\ \mathrm{mod}\ 4).

A complete finite Fourier expansion is

fs(x)=14∑a=14e−isθaχa(x),χa(x)=∑s=±1/2, ±3/2eisθafs(x).\begin{gathered} f_s(x)=\frac14\sum_{a=1}^4e^{-is\theta_a}\chi_a(x),\\ \chi_a(x)=\sum_{s=\pm1/2,\,\pm3/2}e^{is\theta_a}f_s(x). \end{gathered}

Here fsf_s is the coefficient called usu_s on the previous page; the temporary name keeps it distinct from the continuum components u,vu,v. The angles are lifted continuously, so θa+4=θa+2π\theta_{a+4}=\theta_a+2\pi. Advancing a corner frame multiplies its character by eisθe^{is\theta}, whereas rotating the field and its argument actively transforms the coefficient by e−isθe^{-is\theta}. For example, a lattice quarter-turn gives χa′(Rπ/2x)=χa−1(x)\chi'_a(R_{\pi/2}x)=\chi_{a-1}(x) with the cut transported, and hence fs′(Rπ/2x)=e−isπ/2fs(x)f'_s(R_{\pi/2}x)=e^{-is\pi/2}f_s(x). This is the same convention as SθS_\theta above.

The critical projection retains f+1/2,f−1/2f_{+1/2},f_{-1/2}. The other two characters have a nonzero lattice kernel at criticality and can be eliminated in a long-wavelength expansion. They are not extra independent continuum particles or the beginning of an unlimited corner-harmonic tower. Inspect the figure’s original sites, dual face centers and the two successive reductions before following the matrix calculation.

Original spins and dual-face disorder insertions form four corner fields whose critical pair gives the continuum spinor and its massive equations

The four Ising corner fields split into a critical s=±1/2s=\pm1/2 pair and a gapped s=±3/2s=\pm3/2 pair. With the lifted angles and active rotation convention stated here, the critical coefficients give (u,v)(u,v) up to one common normalization, and the leading mass is m=4(K−Kc)/am=4(K-K_c)/a. In the figure, jj labels the corners and aa is lattice spacing. The geometry is schematic; the continuum step requires a∣k∣≪1a|k|\ll1 and ∣K−Kc∣≪1|K-K_c|\ll1.

Projecting the four-component lattice equation

Section titled “Projecting the four-component lattice equation”

Let ℓ\ell denote the lattice spacing in this calculation (ℓ=a\ell=a in the surrounding text), to distinguish it from the corner index aa. Define the antiperiodic shift TT by (Tχ)a=χa+1(T\chi)_a=\chi_{a+1}, including χ5=−χ1\chi_5=-\chi_1. The exact previous-page relation is

χa(x)−cosh⁡(2K)χa+1(x)+sinh⁡(2K)χa+2(x+δa)=0.\chi_a(x)-\cosh(2K)\chi_{a+1}(x) +\sinh(2K)\chi_{a+2}(x+\delta_a)=0.

This is a correlator equation away from spectator contacts. With inverse Fourier phase eik⋅xe^{ik\cdot x}, write it as A(k,K)χ(k)=0A(k,K)\chi(k)=0, where

A(k,K)=I−cosh⁡(2K)T+sinh⁡(2K)E(k)T2,E(k)=diag⁡(eik⋅δ1,…,eik⋅δ4),(δ1,δ2,δ3,δ4)=ℓ((0,1),(−1,0),(0,−1),(1,0)).\begin{gathered} A(k,K)=I-\cosh(2K)T+\sinh(2K)E(k)T^2,\\ E(k)=\operatorname{diag}(e^{ik\cdot\delta_1},\ldots,e^{ik\cdot\delta_4}),\\ (\delta_1,\delta_2,\delta_3,\delta_4) =\ell\big((0,1),(-1,0),(0,-1),(1,0)\big). \end{gathered}

Use the orthonormal matrix Bas=12eisθaB_{as}=\tfrac12e^{is\theta_a} with columns ordered s=(1/2,−1/2,3/2,−3/2)s=(1/2,-1/2,3/2,-3/2). Then B†B=IB^\dagger B=I and B†χ=2fB^\dagger\chi=2f. At zero momentum and critical coupling,

B†A(0,Kc)B=diag⁡(0,0,2−2i,2+2i).B^\dagger A(0,K_c)B =\operatorname{diag}(0,0,2-2i,2+2i).

Write δK=K−Kc\delta K=K-K_c, k±=k1±ik2k_\pm=k_1\pm ik_2, and keep first order in δK\delta K and ℓk\ell k. The critical 2×22\times2 block is

AL(k,K)=(2e3iπ/4δKiℓ2e−iπ/4k−iℓ2eiπ/4k+2e−3iπ/4δK)+O(δK2, ∣δK∣ℓ∣k∣, ℓ2∣k∣2).\begin{aligned} A_{\rm L}(k,K) &=\begin{pmatrix} 2e^{3i\pi/4}\delta K & \tfrac{i\ell}{2}e^{-i\pi/4}k_-\\ \tfrac{i\ell}{2}e^{i\pi/4}k_+ & 2e^{-3i\pi/4}\delta K \end{pmatrix}\\ &\quad+O(\delta K^2,\,|\delta K|\ell|k|,\,\ell^2|k|^2). \end{aligned}

These coefficients can be checked directly from B†TB=diag⁡(eisπ/2)B^\dagger T B=\operatorname{diag}(e^{is\pi/2}), cosh⁡(2K)=2+2δK+O(δK2)\cosh(2K)=\sqrt2+2\delta K+O(\delta K^2) and sinh⁡(2K)=1+22δK+O(δK2)\sinh(2K)=1+2\sqrt2\delta K+O(\delta K^2). The derivative term comes from E(k)=I+idiag⁡(k⋅δa)+O(ℓ2∣k∣2)E(k)=I+i\operatorname{diag}(k\cdot\delta_a)+O(\ell^2|k|^2).

Eliminating the heavy block gives the Schur complement ALL−ALHAHH−1AHLA_{\rm LL}-A_{\rm LH}A_{\rm HH}^{-1}A_{\rm HL}. The heavy eigenvalues remain nonzero, and both mixed blocks vanish at k=0k=0. Their correction therefore starts at O(ℓ2∣k∣2)O(\ell^2|k|^2) and cannot alter the displayed first-order coefficients.

Replace ikiik_i by ∂i\partial_i and multiply the two equations, respectively, by (2/ℓ)eiπ/4(2/\ell)e^{i\pi/4} and (2/ℓ)e−iπ/4(2/\ell)e^{-i\pi/4}. The phases cancel:

(−m∂1−i∂2∂1+i∂2−m)(f+1/2f−1/2)=0,m=4δKℓ.\begin{gathered} \begin{pmatrix} -m & \partial_1-i\partial_2\\ \partial_1+i\partial_2 & -m \end{pmatrix} \begin{pmatrix}f_{+1/2}\\f_{-1/2}\end{pmatrix}=0,\\ m=\frac{4\delta K}{\ell}. \end{gathered}

The omitted operator terms have order δK2/ℓ+∣δK∣∣k∣+ℓ∣k∣2\delta K^2/\ell+|\delta K||k|+\ell|k|^2. In the scaling limit ℓ→0\ell\to0 with mm and physical momentum fixed, δK=mℓ/4\delta K=m\ell/4 and this remainder vanishes. The lattice spacing is now again denoted by aa.

The continuum components may be normalized as

(uv)=Zψa−1/2(f+1/2f−1/2).\begin{pmatrix}u\\v\end{pmatrix} =Z_\psi a^{-1/2} \begin{pmatrix}f_{+1/2}\\f_{-1/2}\end{pmatrix}.

The common dimensionless factor ZψZ_\psi fixes the chosen two-point amplitude; the homogeneous propagation equation does not determine it. No relative field rephasing is needed with these lifted angles. This completes the promised differential and mass normalization. Polyakov 1987, § 10.3.1, pp. 276–278, Eqs. (10.60)–(10.69) gives the lattice equation and soft-mode argument; the explicit matrix calculation above fixes our spacing and phase dictionary. His same-argument transformation u(x)↦eiθ/2u(Rθx)u(x)\mapsto e^{i\theta/2}u(R_\theta x) corresponds to the inverse of the active transformation used here.

The two low modes describe one Lorentzian Majorana field, with two real spinor components, rather than a charged complex Dirac field. In the Euclidean functional integral use independent Grassmann variables u,vu,v and

SM=12∫d2x [u(∂1+i∂2)u−v(∂1−i∂2)v−2muv],A=C(γi∂i−m),C=(01−10).\begin{gathered} S_{\rm M}=\frac12\int d^2x\, \big[u(\partial_1+i\partial_2)u -v(\partial_1-i\partial_2)v-2muv\big],\\ \mathcal A=C(\gamma_i\partial_i-m),\qquad C=\begin{pmatrix}0&1\\-1&0\end{pmatrix}. \end{gathered}

Varying this action gives the two displayed component equations. Since CT=−CC^T=-C and CγiC\gamma_i is symmetric, integration by parts makes A\mathcal A antisymmetric with boundary conditions that remove its surface term. With a fixed Grassmann measure orientation the Gaussian integral is proportional to Pf⁡A\operatorname{Pf}\mathcal A; its sign is not specified by taking an arbitrary square root of a determinant. No pointwise condition v=u∗v=u^* is imposed. This action describes the local bulk scaling theory; a global partition function also requires the spin structure and sector combination inherited from the lattice boundaries.

The Kramers–Wannier relation is

sinh⁡2K sinh⁡2K∗=1.\sinh 2K\,\sinh 2K^*=1.

At the self-dual point,

K∗=K=Kc,sinh⁡2Kc=1.K^*=K=K_c, \qquad \sinh 2K_c=1.

Expanding around criticality gives

K∗−Kc=−(K−Kc)+O((K−Kc)2).K^*-K_c=-(K-K_c)+O((K-K_c)^2).

Therefore the continuum mass changes sign under duality:

m(K∗)=−m(K)+O(am2).m(K^*)=-m(K)+O(a m^2).

With the field phases fixed above, m>0m>0 is the ordered side K>KcK>K_c and m<0m<0 is the disordered side. Duality exchanges long-range order of the spin and disorder variables in the thermodynamic limit. The free-fermion correlation length is ξf=1/∣m∣\xi_{\rm f}=1/|m| in the continuum normalization. Other Ising channels can have different correlation-length amplitudes, but their critical divergence has the same exponent:

ξa∝1∣K−Kc∣.\frac{\xi}{a}\propto\frac1{|K-K_c|}.

Thus the two-dimensional Ising correlation-length exponent is

ν=1.\nu=1.

This agrees with the exact solution and is much sharper than mean-field theory. The deeper reason is that the critical Ising model is not an interacting scalar theory at long distances; it is a free massless Majorana fermion plus nontrivial spin and disorder fields.

The lattice energy is the operator conjugate to the coupling. For the nearest-neighbor square-lattice model,

Z(K)=∑{σ}exp⁡ ⁣(K∑⟨ij⟩σiσj),Z(K)=\sum_{\{\sigma\}}\exp\!\left(K\sum_{\langle ij\rangle}\sigma_i\sigma_j\right),

define the bond energy

ϵij=σiσj.\epsilon_{ij}=\sigma_i\sigma_j.

The dimensionless free energy F=−log⁡ZF=-\log Z satisfies

−∂F∂K=⟨∑⟨ij⟩ϵij⟩,-{\partial F\over\partial K} =\left\langle\sum_{\langle ij\rangle}\epsilon_{ij}\right\rangle,

and

−∂2F∂K2=∑⟨ij⟩,⟨kℓ⟩⟨ϵijϵkℓ⟩c.-{\partial^2 F\over\partial K^2} =\sum_{\langle ij\rangle,\langle k\ell\rangle} \left\langle \epsilon_{ij}\epsilon_{k\ell}\right\rangle_c.

Up to conventional powers of the temperature, this is the specific heat. Translational invariance turns it into an integral of the connected energy–energy correlator:

C∝∑x⟨ϵ(0)ϵ(x)⟩c.C\propto \sum_x \langle\epsilon(0)\epsilon(x)\rangle_c.

The lattice bond has a nonzero identity contribution. Its critical part has an expansion ϵij=c0+c1a ϵcont(x)+⋯\epsilon_{ij}=c_0+c_1a\,\epsilon_{\rm cont}(x)+\cdots, with nonuniversal coefficients; the connected correlator removes c0c_0. In the action convention above,

ϵcont(x)=∂LM∂m=−u(x)v(x)=iψˉ(x)ψ(x),ψ=u,ψˉ=−iv.\epsilon_{\rm cont}(x) =\frac{\partial\mathcal L_{\rm M}}{\partial m} =-u(x)v(x)=i\bar\psi(x)\psi(x), \qquad \psi=u,\quad\bar\psi=-iv.

Here ψˉ\bar\psi is the independent opposite-chirality Grassmann field, not the complex conjugate of ψ\psi. Anticommutation gives i(−iv)u=vu=−uvi(-iv)u=vu=-uv, so the factor of ii is consistent with the already fixed u,vu,v phases. The sign of c1c_1 relates this mass derivative to the particular lattice bond convention; it does not affect the connected two-point scaling.

Equivalently, the energy operator is the operator multiplying the fermion mass in the action. Since the Majorana field has scaling dimension 1/21/2, the bilinear has scaling dimension

Δϵ=1.\Delta_\epsilon=1.

At criticality,

⟨ϵ(x)ϵ(0)⟩c∼1∣x∣2Δϵ=1∣x∣2.\langle\epsilon(x)\epsilon(0)\rangle_c\sim {1\over |x|^{2\Delta_\epsilon}} ={1\over |x|^2}.

Therefore the specific heat diverges logarithmically:

C∼∫aξ2πr dr 1r2∼log⁡ξa∼log⁡1∣K−Kc∣.C\sim \int_a^\xi 2\pi r\,dr\,{1\over r^2} \sim \log {\xi\over a} \sim \log {1\over |K-K_c|}.

The same result appears in momentum space as the one-loop fermion bubble

C∼∫∣m∣1/ad2kk2∼log⁡1∣m∣a.C\sim \int_{|m|}^{1/a}{d^2k\over k^2} \sim \log {1\over |m|a}.

Energy operator, fermion bilinear, and the logarithmic specific heat

After subtracting its identity contribution and fixing the field amplitude, the critical part of the Ising bond energy is proportional to the fermion mass operator iψˉψi\bar\psi\psi. Its two-point function scales as 1/∣x∣21/|x|^2, so the integrated energy–energy correlator gives the logarithmic specific-heat singularity. The diagram represents the scaling argument; it does not fix its lattice amplitude.

Once the correct scaling field is known, the thermodynamic singularity follows from its connected two-point function. The exact lattice solution fixes the nonuniversal amplitude and regular terms; the continuum argument fixes the logarithm. The equivalent fermion-integral calculation appears in Polyakov 1987, § 10.3.1, p. 278, Eq. (10.70) and the following specific-heat formula. These estimates assume a≪ξa\ll\xi and a sample much larger than ξ\xi; at criticality in a finite sample the sample size supplies the infrared cutoff.

From lattice transport to critical thermodynamics

Section titled “From lattice transport to critical thermodynamics”

A continuum field is classified not only by its equation of motion but also by its rotation representation. Scalars transform trivially, vectors transform by the ordinary rotation matrix, and spinors transform by the double cover. The essential spinor fact is S(2π)=−1S(2\pi)=-1.

The Ising order–disorder composite has precisely this sign. Its four corner components are antiperiodic under a full rotation of the disorder endpoint around the order insertion. Their four independent corner characters may be represented by spins ±1/2\pm1/2 and ±3/2\pm3/2 modulo 44. In the scaling limit the leading ±1/2\pm1/2 combinations form a real two-component fermion. The local lattice transfer relation becomes the massive two-dimensional Dirac equation

(∂1+i∂2)u=mv,(∂1−i∂2)v=mu,m=4(K−Kc)ato leading scaling order.(\partial_1+i\partial_2)u=m v, \qquad (\partial_1-i\partial_2)v=m u, \qquad m=\frac{4(K-K_c)}a\quad\text{to leading scaling order}.

At K=KcK=K_c, the mass vanishes and the two components become holomorphic and antiholomorphic Majorana fields. The energy operator is the fermion mass bilinear iψˉψi\bar\psi\psi, whose scaling dimension is 11. Hence the specific heat is logarithmically singular:

C∼log⁡1∣K−Kc∣.C\sim \log {1\over |K-K_c|}.

The four Ising corner fields are not four independent continuum fermions, nor do four samples support an unlimited independent harmonic tower. They form a four-dimensional lattice-label space; the long-distance spin ±1/2\pm1/2 projection gives one Majorana fermion, while derivative descendants account for higher continuum-spin corrections.

The statement χa+4=−χa\chi_{a+4}=-\chi_a is not ordinary anticommutation of microscopic spins. It is monodromy of mixed order–disorder correlation functions. The original Ising variables are still commuting classical spins.

The sign of the continuum mass depends on phase conventions, but its oddness under Kramers–Wannier duality is physical. The two signs encode the ordered and disordered phases.

The Feynman-bubble picture of the specific heat should be understood as a continuum scaling argument. The lattice cutoff aa and the normalization of ϵ\epsilon determine nonuniversal constants, while the logarithmic divergence is universal.

Derive the transversality condition for a massive vector field from the Proca equation

−∂jFji+m2Ai=0,Fij=∂iAj−∂jAi.-\partial_jF_{ji}+m^2A_i=0, \qquad F_{ij}=\partial_iA_j-\partial_jA_i.

Then show that the transverse momentum-space projector is

Pij⊥(k)=δij−kikjk2.P_{ij}^{\perp}(k)=\delta_{ij}-{k_i k_j\over k^2}.
Solution

Taking a divergence of the Proca equation gives

−∂i∂jFji+m2∂iAi=0.-\partial_i\partial_jF_{ji}+m^2\partial_iA_i=0.

The first term vanishes because ∂i∂j\partial_i\partial_j is symmetric in i,ji,j, while FjiF_{ji} is antisymmetric:

∂i∂jFji=0.\partial_i\partial_jF_{ji}=0.

For m≠0m\ne0, this implies

∂iAi=0.\partial_iA_i=0.

In momentum space, the longitudinal part of a vector is proportional to kik_i. The projector onto this direction is

Pij∥=kikjk2.P_{ij}^{\parallel}={k_i k_j\over k^2}.

Therefore the projector onto the orthogonal subspace is

Pij⊥=δij−Pij∥=δij−kikjk2.P_{ij}^{\perp}=\delta_{ij}-P_{ij}^{\parallel} =\delta_{ij}-{k_i k_j\over k^2}.

It obeys

kiPij⊥=kj−k2kjk2=0,k_iP_{ij}^{\perp}=k_j-{k^2k_j\over k^2}=0,

so P⊥AP^{\perp}A is transverse.

Exercise 2: Half-angle rotation and corner characters

Section titled “Exercise 2: Half-angle rotation and corner characters”

Let a two-dimensional spinor transform under rotations by

Sθ(u)=e−iθ/2u,Sθ(v)=eiθ/2v.S_\theta(u)=e^{-i\theta/2}u, \qquad S_\theta(v)=e^{i\theta/2}v.

Use the active convention Ψ′(Rθx)=SθΨ(x)\Psi'(R_\theta x)=S_\theta\Psi(x). Check Sθ−1γ1Sθ=cos⁡θ γ1−sin⁡θ γ2S_\theta^{-1}\gamma_1S_\theta=\cos\theta\,\gamma_1-\sin\theta\,\gamma_2, then show that a 2π2\pi rotation gives a minus sign. Explain why a continuously defined antiperiodic angular field has half-integer harmonics, and why four corner directions contain only four independent characters, labeled by s=±1/2,±3/2s=\pm1/2,\pm3/2 modulo 44.

Solution

Direct multiplication gives

Sθ−1γ1Sθ=(0eiθe−iθ0)=cos⁡θ γ1−sin⁡θ γ2.S_\theta^{-1}\gamma_1S_\theta =\begin{pmatrix}0&e^{i\theta}\\e^{-i\theta}&0\end{pmatrix} =\cos\theta\,\gamma_1-\sin\theta\,\gamma_2.

For θ=2π\theta=2\pi,

eiθ/2=eiπ=−1,e−iθ/2=e−iπ=−1.e^{i\theta/2}=e^{i\pi}=-1, \qquad e^{-i\theta/2}=e^{-i\pi}=-1.

Thus

S2π(u)=−u,S2π(v)=−v.S_{2\pi}(u)=-u, \qquad S_{2\pi}(v)=-v.

Now suppose a function of the corner angle has an angular harmonic eisθe^{is\theta}. Under θ↦θ+2π\theta\mapsto\theta+2\pi it transforms as

eis(θ+2π)=e2πiseisθ.e^{is(\theta+2\pi)}=e^{2\pi i s}e^{is\theta}.

Antiperiodicity requires

e2πis=−1.e^{2\pi i s}=-1.

Therefore

s∈12+Z.s\in {1\over2}+\mathbb Z.

Thus a continuously defined antiperiodic angular field has half-integer harmonics. On the four corner angles θa=θ0+πa/2\theta_a=\theta_0+\pi a/2, however,

ei(s+4)θa=e4iθ0eisθa,e^{i(s+4)\theta_a}=e^{4i\theta_0}e^{is\theta_a},

so ss and s+4s+4 give the same corner character up to an aa-independent phase. There are only four independent characters. One may choose the representatives

s=+12, −12, +32, −32(mod 4).s=+{1\over2},\ -{1\over2},\ +{3\over2},\ -{3\over2} \qquad (\mathrm{mod}\ 4).

Using

γ1=(0110),γ2=(0−ii0),Ψ=(uv),\gamma_1=\begin{pmatrix}0&1\\1&0\end{pmatrix}, \qquad \gamma_2=\begin{pmatrix}0&-i\\i&0\end{pmatrix}, \qquad \Psi=\begin{pmatrix}u\\v\end{pmatrix},

show that

(γi∂i−m)Ψ=0(\gamma_i\partial_i-m)\Psi=0

is equivalent to

(∂1+i∂2)u=mv,(∂1−i∂2)v=mu.(\partial_1+i\partial_2)u=m v, \qquad (\partial_1-i\partial_2)v=m u.

What happens at m=0m=0?

Solution

Compute

γ1∂1Ψ=(∂1v∂1u),γ2∂2Ψ=(−i∂2vi∂2u).\gamma_1\partial_1\Psi =\begin{pmatrix}\partial_1 v\\ \partial_1 u\end{pmatrix}, \qquad \gamma_2\partial_2\Psi =\begin{pmatrix}-i\partial_2 v\\ i\partial_2 u\end{pmatrix}.

Therefore

γi∂iΨ=((∂1−i∂2)v(∂1+i∂2)u).\gamma_i\partial_i\Psi =\begin{pmatrix} (\partial_1-i\partial_2)v\\ (\partial_1+i\partial_2)u \end{pmatrix}.

The equation γi∂iΨ=mΨ\gamma_i\partial_i\Psi=m\Psi gives

(∂1−i∂2)v=mu,(∂1+i∂2)u=mv.(\partial_1-i\partial_2)v=m u, \qquad (\partial_1+i\partial_2)u=m v.

These are the desired equations, just written in the opposite order.

At m=0m=0 they become

(∂1+i∂2)u=0,(∂1−i∂2)v=0.(\partial_1+i\partial_2)u=0, \qquad (\partial_1-i\partial_2)v=0.

Since

∂ˉ=12(∂1+i∂2),∂=12(∂1−i∂2),\bar\partial={1\over2}(\partial_1+i\partial_2), \qquad \partial={1\over2}(\partial_1-i\partial_2),

we get

∂ˉu=0,∂v=0.\bar\partial u=0, \qquad \partial v=0.

Thus uu is holomorphic and vv is antiholomorphic.

Use the Kramers–Wannier relation

sinh⁡2K sinh⁡2K∗=1\sinh 2K\,\sinh 2K^*=1

to show that near the self-dual point,

K∗−Kc=−(K−Kc)+O((K−Kc)2).K^*-K_c=-(K-K_c)+O((K-K_c)^2).

Why does this imply that the Ising fermion mass is duality-odd?

Solution

Set

K=Kc+δK,K∗=Kc+δK∗.K=K_c+\delta K, \qquad K^*=K_c+\delta K^*.

At criticality,

sinh⁡2Kc=1.\sinh 2K_c=1.

Expanding to first order,

sinh⁡2K=1+2cosh⁡2Kc δK+O(δK2),\sinh 2K =1+2\cosh 2K_c\,\delta K+O(\delta K^2),

and similarly

sinh⁡2K∗=1+2cosh⁡2Kc δK∗+O((δK∗)2).\sinh 2K^* =1+2\cosh 2K_c\,\delta K^*+O((\delta K^*)^2).

Multiplying and using sinh⁡2Ksinh⁡2K∗=1\sinh 2K\sinh 2K^*=1 gives

1+2cosh⁡2Kc(δK+δK∗)+O(δK2)=1.1+2\cosh 2K_c(\delta K+\delta K^*)+O(\delta K^2)=1.

Therefore

δK∗=−δK+O(δK2),\delta K^*=-\delta K+O(\delta K^2),

or

K∗−Kc=−(K−Kc)+O((K−Kc)2).K^*-K_c=-(K-K_c)+O((K-K_c)^2).

The continuum mass is proportional to the relevant perturbation away from criticality:

m∝K−Kc.m\propto K-K_c.

Since duality reverses K−KcK-K_c, it reverses mm. Thus the mass term is duality-odd.

Assume the energy operator at the critical Ising point has scaling dimension Δϵ=1\Delta_\epsilon=1. Show that the integrated energy–energy correlator diverges logarithmically. Then express the answer in terms of the correlation length ξ\xi and the mass mm.

Solution

At criticality, an operator of scaling dimension Δϵ\Delta_\epsilon has two-point function

⟨ϵ(x)ϵ(0)⟩c∼1∣x∣2Δϵ.\langle\epsilon(x)\epsilon(0)\rangle_c\sim {1\over |x|^{2\Delta_\epsilon}}.

For Δϵ=1\Delta_\epsilon=1,

⟨ϵ(x)ϵ(0)⟩c∼1∣x∣2.\langle\epsilon(x)\epsilon(0)\rangle_c\sim {1\over |x|^2}.

The specific heat is proportional to the integral of this correlator. With a UV cutoff aa and an IR cutoff ξ\xi,

C∼∫aξ2πr dr 1r2=2π∫aξdrr=2πlog⁡ξa,C\sim \int_a^\xi 2\pi r\,dr\,{1\over r^2} =2\pi\int_a^\xi {dr\over r} =2\pi\log {\xi\over a},

up to the normalization of ϵ\epsilon.

The mass gap cuts off correlations at

ξ∼1∣m∣.\xi\sim {1\over |m|}.

Since m∝K−Kcm\propto K-K_c,

C∼log⁡1∣m∣a∼log⁡1∣K−Kc∣C\sim \log {1\over |m|a} \sim \log {1\over |K-K_c|}

plus nonuniversal regular terms.

  • A. M. Polyakov, Gauge Fields and Strings. Contemporary Concepts in Physics, Vol. 3. Harwood Academic Publishers, 1987, § 10.3.1, doi:10.1201/9780203755082.
  • P. Di Francesco, P. Mathieu, and D. Sénéchal, Conformal Field Theory. Graduate Texts in Contemporary Physics. Springer, 1997, Chapters 4 and 12.
  • B. M. McCoy and T. T. Wu, The Two-Dimensional Ising Model. Harvard University Press, 1973.
  • T. D. Schultz, D. C. Mattis, and E. H. Lieb, “Two-Dimensional Ising Model as a Soluble Problem of Many Fermions,” Reviews of Modern Physics 36, 856–871 (1964), doi:10.1103/RevModPhys.36.856.

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