Skip to content

Transverse-Field Ising Model

The one-dimensional transverse-field Ising model is a chain of spin-1/21/2 degrees of freedom in which ferromagnetic Ising exchange competes with a noncommuting transverse field, yielding an exactly solvable quantum phase transition between an ordered phase and a quantum paramagnet.

This dossier fixes the model record, conventions, solution-status claims, observables, limiting cases, and validation targets. The Transverse-Field Ising Model teaching article owns the detailed physical development and quasiparticle results. Jordan–Wigner Transformation owns the nonlocal spin-to-fermion derivation, and Quantum Phase Transitions owns the general scaling framework.

Unless a variant is named explicitly, this dossier uses the following model.

FieldBaseline choice
geometrya one-dimensional chain of LL sites
local degree of freedomone spin-1/21/2 space C2\mathbb C^2 per site
operatorsdimensionless Pauli matrices σjα\sigma_j^\alpha with eigenvalues ±1\pm1
interactionnearest-neighbor Ising exchange along zz
fielduniform transverse field along xx
couplingsJ>0J>0 and h≥0h\ge0
dimensionless controlg=h/Jg=h/J when J≠0J\ne0
boundariesopen or periodic, stated explicitly
energy offsetnone
longitudinal fieldabsent
disorder and long-range termsabsent
equilibrium limitL→∞L\to\infty after the finite-chain problem is defined

For bulk formulas, the periodic chain with L≥3L\ge3 is the default unless stated otherwise. Small periodic chains require extra care: for L=2L=2, writing both j=1j=1 and j=2j=2 terms in a nearest-neighbor sum counts the same undirected bond twice. The numerical contracts below therefore specify their bond lists rather than relying on an ambiguous shorthand.

The sign of hh is inessential for this isolated model because a global π\pi rotation about zz sends σjx↦−σjx\sigma_j^x\mapsto-\sigma_j^x while leaving every Ising bond unchanged. The restriction h≥0h\ge0 removes this unitary redundancy. The sign of JJ has more substantial boundary consequences and is discussed under variants.

The local and many-body Hilbert spaces are

Hj≃C2,HL=⨂j=1LHj,\mathcal H_j \simeq \mathbb C^2, \qquad \mathcal H_L = \bigotimes_{j=1}^{L}\mathcal H_j,

so

dim⁡HL=2L.\dim\mathcal H_L = 2^L.

The computational basis is the simultaneous eigenbasis of all σjz\sigma_j^z:

σz∣0⟩=+∣0⟩,σz∣1⟩=−∣1⟩.\begin{aligned} \sigma^z\lvert0\rangle &= +\lvert0\rangle, \\ \sigma^z\lvert1\rangle &= -\lvert1\rangle. \end{aligned}

A basis vector can therefore be labeled by a bit string

∣s1s2⋯sL⟩,sj∈{0,1}.\lvert s_1s_2\cdots s_L\rangle, \qquad s_j\in\{0,1\}.

The exchange term is diagonal in this basis. The transverse-field operator σjx\sigma_j^x flips bit jj, so the full Hamiltonian connects bit strings that differ at one site. This gives an immediate sparse-matrix construction with at most L+1L+1 nonzero entries per row before duplicate contributions are combined.

The Spin-1/21/2 Chain dossier owns the generic tensor-product ordering, local-operator embedding, and family-wide matrix checks. Here they are specialized to one Ising coupling and one transverse field.

For open boundary conditions,

HO=−J∑j=1L−1σjzσj+1z−h∑j=1Lσjx.\begin{aligned} H_{\mathrm O} ={}& -J \sum_{j=1}^{L-1} \sigma_j^z\sigma_{j+1}^z \\ &- h \sum_{j=1}^{L} \sigma_j^x. \end{aligned}

For periodic boundary conditions,

HP=−J∑j=1Lσjzσj+1z−h∑j=1Lσjx,\begin{aligned} H_{\mathrm P} ={}& -J \sum_{j=1}^{L} \sigma_j^z\sigma_{j+1}^z \\ &- h \sum_{j=1}^{L} \sigma_j^x, \end{aligned}

with

σL+1α=σ1α.\sigma_{L+1}^{\alpha} = \sigma_1^{\alpha}.

The open chain has L−1L-1 exchange bonds; an ordinary ring has LL. Both JJ and hh have units of energy. For J≠0J\ne0, the overall scale and dimensionless competition are separated by

H=JH^(g),g=hJ.H = J\widehat H(g), \qquad g = \frac{h}{J}.

This dossier never uses σα\sigma^\alpha and a spin operator interchangeably. With dimensionless spin operators

sjα=12σjα,s_j^\alpha = \frac12\sigma_j^\alpha,

the same Hamiltonian becomes

H=−4J∑⟨j,k⟩sjzskz−2h∑jsjx.H = -4J \sum_{\langle j,k\rangle} s_j^z s_k^z - 2h \sum_j s_j^x.

If dimensionful angular-momentum operators Sjα=ℏsjαS_j^\alpha=\hbar s_j^\alpha are used, the coupling constants acquire the corresponding powers of ℏ−1\hbar^{-1}. A critical ratio quoted without the operator normalization is not a complete statement.

Some authors write the exchange in terms of domain-wall projectors,

nj,j+1DW=1−σjzσj+1z2.n_{j,j+1}^{\mathrm{DW}} = \frac{1-\sigma_j^z\sigma_{j+1}^z}{2}.

Then

−Jσjzσj+1z=−J+2Jnj,j+1DW.-J\sigma_j^z\sigma_{j+1}^z = -J + 2J n_{j,j+1}^{\mathrm{DW}}.

The projector form contains an explicit constant −J-J per bond. The baseline Hamiltonian above retains that constant; no later recentering is implied.

Define the global spin-flip operator

P=∏j=1Lσjx.P = \prod_{j=1}^{L} \sigma_j^x.

Because operators on different sites commute,

P2=I,P^2 = \mathbb I,

and conjugation gives

PσjzP=−σjz,PσjxP=σjx.\begin{aligned} P\sigma_j^zP &= -\sigma_j^z, \\ P\sigma_j^xP &= \sigma_j^x. \end{aligned}

Every Ising bond contains two σz\sigma^z factors, so

[H,P]=0.[H,P] = 0.

The Hilbert space splits into parity sectors P=±1P=\pm1. For L≥1L\ge1, each sector has dimension

dim⁡H+=dim⁡H−=2L−1.\dim\mathcal H_{+} = \dim\mathcal H_{-} = 2^{L-1}.

The longitudinal magnetization is odd under PP, while the energy, transverse magnetization, and even-zz correlators are parity even.

The uniform periodic chain is invariant under one-site translations and reflections. The uniform open chain retains reflection about its midpoint but not translation symmetry. Disorder, nonuniform fields, boundary pinning, or a defect bond can remove these symmetries without necessarily removing Ising parity.

In the computational basis, HH is a real symmetric matrix, so complex conjugation KK is an antiunitary symmetry of that matrix representation. This should not be confused with physical spin time reversal. Physical time reversal sends every spin component to its negative; therefore a nonzero transverse field changes sign under it. The exchange-only model is time-reversal invariant, whereas the field selects an oriented external magnetic field.

Neither total zz magnetization nor total xx magnetization commutes with the generic Hamiltonian:

[H,∑jσjz]≠0,[H,∑jσjx]≠0\left[ H, \sum_j\sigma_j^z \right] \ne 0, \qquad \left[ H, \sum_j\sigma_j^x \right] \ne 0

when Jh≠0Jh\ne0. Ising parity, rather than a continuous spin component, is the model’s internal conserved quantum number.

RegimeGround-state or low-energy pictureExact statement
h=0h=0, J>0J>0two aligned zz configurationsclassical commuting Ising limit
J=0J=0, h>0h>0unique product state polarized along +x+xindependent-site limit
0<g<10<g<1ordered thermodynamic phasetwo symmetry-related pure phases; finite chains retain parity eigenstates
g=1g=1scale-invariant critical pointgap closes linearly; Ising universality class
g>1g>1quantum paramagnetunique bulk phase adiabatically connected to the +x+x product state
h<0h<0reversed transverse polarizationunitarily equivalent to h>0h>0
J<0J<0antiferromagnetic Ising exchangeequivalent to J>0J>0 on a bipartite chain, subject to boundary frustration
longitudinal field nonzeroexplicit Ising-symmetry breakinggeneric lattice model is no longer quadratic-fermion integrable

At h=0h=0, the ferromagnetic ground configurations are

∣↑↑⋯↑⟩and∣↓↓⋯↓⟩.\lvert\uparrow\uparrow\cdots\uparrow\rangle \quad\text{and}\quad \lvert\downarrow\downarrow\cdots\downarrow\rangle.

For a periodic chain their common energy is

E0=−JL,E_0 = -JL,

and for an open chain it is −J(L−1)-J(L-1). A flipped domain costs 2J2J at each domain wall. Periodic boundary conditions force the total number of domain walls to be even.

At J=0J=0 and h>0h>0, the unique ground state is

∣→→⋯→⟩,\lvert\rightarrow\rightarrow\cdots\rightarrow\rangle,

where

σx∣→⟩=+∣→⟩.\sigma^x\lvert\rightarrow\rangle = +\lvert\rightarrow\rangle.

Its energy is

E0=−hL.E_0 = -hL.

The first excited manifold contains one site polarized along −x-x and lies 2h2h above the ground state.

Conjugating every even site by σx\sigma^x reverses σjz\sigma_j^z on that sublattice while leaving σjx\sigma_j^x unchanged. It therefore maps J<0J<0 to J>0J>0 on an open chain and on an even periodic ring. An odd periodic ring is frustrated and cannot be mapped to an unfrustrated ferromagnetic ring without leaving a twisted bond.

Adding

Vz=−λ∑jσjzV_z = -\lambda \sum_j\sigma_j^z

breaks PP explicitly because PVzP=−VzPV_zP=-V_z. It also destroys the generic free-fermion closure of the lattice Hamiltonian. The resulting model may still have controlled limits, continuum descriptions, or special integrable field-theory regimes, but those do not make the generic finite lattice model exactly solvable.

The uniform one-dimensional model with nearest-neighbor exchange and no longitudinal field is exactly reducible to independent fermionic quasiparticles. The standard route is

spins→Jordan–Wignerquadratic fermions→Fouriermomentum pairs→Bogoliubovnormal modes.\text{spins} \xrightarrow{\text{Jordan–Wigner}} \text{quadratic fermions} \xrightarrow{\text{Fourier}} \text{momentum pairs} \xrightarrow{\text{Bogoliubov}} \text{normal modes}.

The phrase “exactly solvable” is quantity dependent.

QuantityStatus in the baseline chain
finite-size spectrumexact after boundary and parity sectors are fixed
bulk quasiparticle dispersionelementary closed form
thermodynamic ground-state energyone-dimensional integral
gap and critical pointexact
spontaneous order parameterexact in the thermodynamic ordered phase
partition functionexact from independent quasiparticles
local spin correlationsexact determinant or Pfaffian representations; asymptotics require further analysis
entanglement entropyexact correlation-matrix construction for fermionic Gaussian states
generic real-time Gaussian observablesreducible to mode evolution and Wick contractions
model with longitudinal fieldnot generically a quadratic-fermion problem
disordered or long-range variantssolvability must be reassessed model by model

Exact diagonalization of a finite 2L2^L matrix is not, by itself, an analytic integrability claim. Conversely, exact quasiparticle energies do not make every nonlocal spin correlation a one-line calculation.

For the baseline convention and the thermodynamic chain, the positive quasiparticle dispersion is

ε(q)=2J1+g2−2gcos⁡q,q∈[−π,π).\varepsilon(q) = 2J \sqrt{ 1+g^2-2g\cos q }, \qquad q\in[-\pi,\pi).

Its minimum gives the bulk gap

Δbulk=2J∣1−g∣.\Delta_{\mathrm{bulk}} = 2J\lvert1-g\rvert.

The gap closes at

gc=1,g_c = 1,

or h=Jh=J in this Pauli convention. At criticality,

ε(q)=4J∣sin⁡q2∣.\varepsilon(q) = 4J \left| \sin\frac q2 \right|.

For small physical wave number k=q/ak=q/a,

ε∼2Ja∣k∣,\varepsilon \sim 2Ja\lvert k\rvert,

so the emergent velocity is

v=2Jaℏ.v = \frac{2Ja}{\hbar}.

The thermodynamic ground-state energy density is

e0(g)=−Jπ∫0π1+g2−2gcos⁡q dq.e_0(g) = -\frac{J}{\pi} \int_0^\pi \sqrt{ 1+g^2-2g\cos q } \, \mathrm dq.

Useful exact critical data are

z=1,ν=1,β=18,η=14,c=12.z=1, \qquad \nu=1, \qquad \beta=\frac18, \qquad \eta=\frac14, \qquad c=\frac12.

Here β\beta is the order-parameter exponent, not inverse temperature, and cc is the central charge of the critical continuum theory. The spontaneous longitudinal magnetization in a selected thermodynamic pure phase is

mz=(1−g2)1/8,0≤g<1.m_z = \left(1-g^2\right)^{1/8}, \qquad 0\le g<1.

It vanishes for g≥1g\ge1. This nonzero value belongs to a symmetry-broken thermodynamic state. The exact parity eigenstate of any finite chain has ⟨Mz⟩=0\langle M_z\rangle=0.

The Jordan–Wigner string turns the spin boundary condition into a parity-dependent fermion boundary condition. For a periodic spin chain, the allowed momenta therefore depend on the eigenvalue of PP. One must diagonalize the correct momentum grid in each parity sector and then impose the sector’s fermion-parity constraint.

Three finite-size energy differences are often called “the gap”:

  1. the splitting between the lowest even- and odd-parity states;
  2. the lowest excitation energy within a fixed parity sector;
  3. the bulk quasiparticle gap inferred after taking L→∞L\to\infty.

They need not agree at finite LL. In the ordered phase, the first can be exponentially small even while the bulk excitation gap remains nonzero. At criticality, the low-lying spacings scale as L−1L^{-1} and depend on boundary conditions and conformal sector.

At h=0h=0 on a periodic ring, a local spin flip creates two domain walls and costs 4J4J. The dispersion formula approaches a single-quasiparticle energy 2J2J as g→0g\to0, but parity and boundary constraints decide which combinations are physical in a specified finite spin sector. This is a standard example of why the bulk dispersion alone is not a complete finite-size spectrum.

See Boundary Conditions on Lattices and Symmetry Sectors for the general bookkeeping.

The ground-state energy, excitation energies, and sector-resolved gaps diagnose limiting states, criticality, and implementation errors. A complete finite-size report states both boundary conditions and the symmetry sectors being compared.

Define the intensive operators

Mz=1L∑j=1Lσjz,Mx=1L∑j=1Lσjx.M_z = \frac1L \sum_{j=1}^{L} \sigma_j^z, \qquad M_x = \frac1L \sum_{j=1}^{L} \sigma_j^x.

MzM_z is odd under PP and is the order parameter. MxM_x is parity even and measures polarization along the applied field. In a finite parity eigenstate,

⟨Mz⟩=0,\langle M_z\rangle = 0,

so ordered behavior is instead exposed by ⟨Mz2⟩\langle M_z^2\rangle, long-distance correlations, a small symmetry-breaking source, or the near-degeneracy of opposite-parity states.

Equal-time correlators include

Czz(r)=⟨σjzσj+rz⟩C_{zz}(r) = \langle \sigma_j^z\sigma_{j+r}^z \rangle

and the connected function

Czzc(r)=Czz(r)−⟨σjz⟩⟨σj+rz⟩.C_{zz}^{\mathrm c}(r) = C_{zz}(r) - \langle\sigma_j^z\rangle \langle\sigma_{j+r}^z\rangle.

For a translation-invariant ring, one useful convention for the static structure factor is

Szz(q)=1L∑j,keiq(j−k)⟨σjzσkz⟩.S_{zz}(q) = \frac1L \sum_{j,k} e^{iq(j-k)} \langle \sigma_j^z\sigma_k^z \rangle.

Its normalization must be stated when comparing codes or references. Correlation Functions Overview and Structure Factors own the general conventions.

For a bond set E\mathcal E,

NDW=∑(j,k)∈E1−σjzσkz2.N_{\mathrm{DW}} = \sum_{(j,k)\in\mathcal E} \frac{ 1-\sigma_j^z\sigma_k^z }{2}.

This operator counts antiparallel bonds in a computational-basis configuration and measures their quantum expectation in a superposition. It is directly tied to the exchange energy.

For a contiguous region AA, the von Neumann entropy is

SA=−Tr⁡(ρAln⁡ρA).S_A = -\operatorname{Tr} \left( \rho_A\ln\rho_A \right).

Away from criticality it saturates with subsystem size in the ground state of the gapped chain. At the critical point it grows logarithmically, with a coefficient controlled by c=1/2c=1/2 and by the boundary geometry. Entanglement and Criticality owns the general scaling laws.

Exchange and field favor incompatible local descriptions. The exchange term prefers definite σz\sigma^z alignment, while the field term continually mixes those configurations through spin flips. Varying gg changes the ground state even at zero temperature because it changes quantum fluctuations rather than thermal populations.

For a finite chain with h>0h>0, the exact ground state has definite parity. Deep in the ordered regime it approaches an even superposition of the two aligned configurations, while the lowest odd state becomes nearly degenerate. Spontaneous symmetry breaking requires an order of limits:

lim⁡λ→0+lim⁡L→∞⟨Mz⟩λ,L≠0,\lim_{\lambda\to0^+} \lim_{L\to\infty} \langle M_z\rangle_{\lambda,L} \ne 0,

where λ\lambda is a longitudinal source. Reversing the limits leaves the finite-system parity symmetry intact. Spontaneous Symmetry Breaking develops this distinction.

At small gg, low-energy excitations are naturally described as fluctuating domain walls. The exact fermionic quasiparticles are the delocalized normal modes of those defects. Their statistics are a property of the nonlocal solution variables; the original Hilbert space remains a tensor product of spins.

At g=1g=1, the correlation length and correlation time diverge, the gap closes, and long-distance observables lose sensitivity to microscopic details. The critical point belongs to the two-dimensional classical Ising universality class and is described by a c=1/2c=1/2 conformal field theory in the scaling limit.

No finite-temperature ordered phase in one dimension

Section titled “No finite-temperature ordered phase in one dimension”

For the short-range one-dimensional chain, every nonzero temperature introduces a finite density of domain walls. The finite-temperature correlation length can be large at low temperature, but true long-range Ising order does not persist for T>0T>0. This does not contradict the zero-temperature quantum critical point.

Minimal Worked Example: Two-Site Open Chain

Section titled “Minimal Worked Example: Two-Site Open Chain”

For L=2L=2 with one open bond,

H=−Jσ1zσ2z−h(σ1x+σ2x).H = -J\sigma_1^z\sigma_2^z - h \left( \sigma_1^x+\sigma_2^x \right).

Introduce parity-adapted states

∣F±⟩=∣00⟩±∣11⟩2,∣A±⟩=∣01⟩±∣10⟩2.\begin{aligned} \lvert F_{\pm}\rangle &= \frac{ \lvert00\rangle \pm \lvert11\rangle }{\sqrt2}, \\ \lvert A_{\pm}\rangle &= \frac{ \lvert01\rangle \pm \lvert10\rangle }{\sqrt2}. \end{aligned}

The P=+1P=+1 block in the basis {∣F+⟩,∣A+⟩}\{\lvert F_+\rangle,\lvert A_+\rangle\} is

H+=(−J−2h−2hJ).H_+ = \begin{pmatrix} -J & -2h \\ -2h & J \end{pmatrix}.

The two negative-parity states do not mix:

H∣F−⟩=−J∣F−⟩,H∣A−⟩=+J∣A−⟩.\begin{aligned} H\lvert F_-\rangle &= -J\lvert F_-\rangle, \\ H\lvert A_-\rangle &= +J\lvert A_-\rangle. \end{aligned}

Therefore the complete spectrum is

spec⁡(H)={−J2+4h2,−J,+J,+J2+4h2}.\operatorname{spec}(H) = \left\{ -\sqrt{J^2+4h^2}, -J, +J, +\sqrt{J^2+4h^2} \right\}.

At J=1J=1 and h=0.7h=0.7, the sorted dimensionless energies are

EJ∈{−1.720465053409,−1,+1,+1.720465053409}.\frac EJ \in \left\{ -1.720465053409, -1, +1, +1.720465053409 \right\}.

This example tests operator normalization, basis ordering, parity, Hermiticity, and the open-bond convention without requiring the many-mode exact solution.

MB-B001 fixes the complete small-chain validation data. Its stronger assembly case uses L=4L=4, periodic boundaries, J=1J=1, and h=0.7h=0.7. Required fingerprints include

dim⁡H=16,E0J=−4.563856065203,Tr⁡H=0,116Tr⁡H2=4(J2+h2)=5.96,[H,P]=0.\begin{aligned} \dim\mathcal H &= 16, \\ \frac{E_0}{J} &= -4.563856065203, \\ \operatorname{Tr}H &= 0, \\ \frac1{16} \operatorname{Tr}H^2 &= 4(J^2+h^2) = 5.96, \\ [H,P] &= 0. \end{aligned}

The benchmark page owns the complete level-and-degeneracy table. A code should pass the L=2L=2 open-chain test before the periodic L=4L=4 test, because a wrong two-site ring can accidentally conceal double-counted bonds.

The second moment follows without diagonalization. Distinct nonidentity Pauli strings are orthogonal under the Hilbert–Schmidt inner product, so for an ordinary LL-site ring,

12LTr⁡H2=L(J2+h2).\frac1{2^L} \operatorname{Tr}H^2 = L(J^2+h^2).

This trace identity checks all matrix entries, not merely the lowest eigenvalue.

MB-B009 uses the critical open chain h=J>0h=J>0. Its exact global gap is

ΔL=4Jsin⁡(π4L+2).\Delta_L = 4J \sin \left( \frac{\pi}{4L+2} \right).

For example,

Δ8J=0.369073437853208.\frac{\Delta_8}{J} = 0.369073437853208.

Writing ℓ=L+12\ell=L+\tfrac12 gives

ℓΔLπJ=sin⁡xx,x=π4ℓ,\frac{ \ell\Delta_L }{ \pi J } = \frac{\sin x}{x}, \qquad x = \frac{\pi}{4\ell},

and hence

ℓΔLπJ=1−π296ℓ2+O(ℓ−4).\frac{ \ell\Delta_L }{ \pi J } = 1 - \frac{\pi^2}{96\ell^2} + O(\ell^{-4}).

This contract tests the dynamical exponent z=1z=1 and the leading boundary-aware correction. It is a global gap, not a gap restricted to one parity sector.

The artifact catalog reserves

notebooks/many-body-statistical/transverse_field_ising_ed.ipynb

for the spin-bit Hamiltonian and complete small-chain spectrum. The notebook is planned rather than silently treated as published. Until it passes the release gates on Reproducible Notebooks, the versioned benchmark contracts are the authoritative numerical targets.

VariantWhat changesCanonical route
classical Ising chainremove the noncommuting transverse field and study thermal configurationsIsing Chain model card
antiferromagnetic chainchange the exchange sign and track frustration or sublattice rotationTransverse-Field Ising Model
longitudinal-field chainbreak Ising parity and generic free-fermion integrabilityQuantum Phase Transitions
XY chainadd unequal xx- and yy-exchange couplingsExact Solutions Preview
random Ising chainmake bonds or fields site dependentseparate disorder and strong-disorder methods are required
long-range Ising modelreplace nearest-neighbor exchange by algebraically decaying couplingsuniversality and causal structure can change
higher-dimensional Ising modelchange the bond graphno generic Jordan–Wigner free-fermion solution
driven or quenched chainspecify an initial state and time-dependent protocolnonequilibrium observables require their own contract

The Ising Chain Hamiltonian card is the compact formula locator. This dossier is the convention-complete record. The teaching article remains the place for derivation and interpretation.

  • Calling the exchange-only diagonal Hamiltonian and the transverse-field quantum model the same problem without qualification.
  • Replacing Pauli matrices by sα=σα/2s^\alpha=\sigma^\alpha/2 without changing the couplings.
  • Writing a periodic nearest-neighbor sum for L=2L=2 and counting the same physical bond twice.
  • Quoting gc=1g_c=1 without stating the Hamiltonian and operator normalization.
  • Assuming a finite parity eigenstate has nonzero longitudinal magnetization.
  • Confusing the exponentially small parity-partner splitting with the nonzero bulk quasiparticle gap in the ordered phase.
  • Using one fermionic momentum grid for both parity sectors of a periodic spin chain.
  • Interpreting exact quasiparticle energies as proof that every spin correlation has an elementary closed form.
  • Treating a longitudinal-field perturbation as though the quadratic solution survived unchanged.
  • Inferring a finite-temperature ordered phase from a large but finite low-temperature correlation length.

Show directly that P=∏jσjxP=\prod_j\sigma_j^x commutes with the open-chain Hamiltonian. Determine the parity of MzM_z and MxM_x.

Solution

On one site,

σxσzσx=−σz,σxσxσx=σx.\sigma^x\sigma^z\sigma^x = -\sigma^z, \qquad \sigma^x\sigma^x\sigma^x = \sigma^x.

Conjugation by PP therefore changes the sign of each σjz\sigma_j^z. A bond contains two such factors:

P(σjzσj+1z)P=σjzσj+1z.P \left( \sigma_j^z\sigma_{j+1}^z \right) P = \sigma_j^z\sigma_{j+1}^z.

Each σjx\sigma_j^x is also invariant, so PHP=HPHP=H and [H,P]=0[H,P]=0. The same conjugation gives

PMzP=−Mz,PMxP=Mx.PM_zP = -M_z, \qquad PM_xP = M_x.

Thus MzM_z is parity odd and MxM_x is parity even.

An article writes

H′=−Js∑jsjzsj+1z−hs∑jsjx,H' = -J_s \sum_j s_j^zs_{j+1}^z - h_s \sum_j s_j^x,

with sα=σα/2s^\alpha=\sigma^\alpha/2. Express JsJ_s and hsh_s in terms of this dossier’s JJ and hh. At what ratio hs/Jsh_s/J_s does the critical point occur?

Solution

Using

σjzσj+1z=4sjzsj+1z,σjx=2sjx,\sigma_j^z\sigma_{j+1}^z = 4s_j^zs_{j+1}^z, \qquad \sigma_j^x = 2s_j^x,

gives

Js=4J,hs=2h.J_s = 4J, \qquad h_s = 2h.

The baseline critical point is h=Jh=J, so

hsJs∣c=2J4J=12.\left. \frac{h_s}{J_s} \right|_{\mathrm c} = \frac{2J}{4J} = \frac12.

The apparent factor-of-two shift is entirely a normalization change.

Use the parity-adapted basis to derive the complete L=2L=2 open-chain spectrum. Which sector contains the ground state when h>0h>0?

Solution

The positive-parity block is

H+=(−J−2h−2hJ),H_+ = \begin{pmatrix} -J & -2h \\ -2h & J \end{pmatrix},

whose characteristic polynomial is

det⁡(EI−H+)=E2−J2−4h2.\det(E\mathbb I-H_+) = E^2-J^2-4h^2.

Its eigenvalues are ±J2+4h2\pm\sqrt{J^2+4h^2}. The two negative-parity states ∣F−⟩\lvert F_-\rangle and ∣A−⟩\lvert A_-\rangle have energies −J-J and +J+J. Therefore

spec⁡(H)={−J2+4h2,−J,+J,+J2+4h2}.\operatorname{spec}(H) = \left\{ -\sqrt{J^2+4h^2}, -J, +J, +\sqrt{J^2+4h^2} \right\}.

For h>0h>0,

J2+4h2>∣J∣,\sqrt{J^2+4h^2} > \lvert J\rvert,

so the ground state is the negative eigenvalue of the P=+1P=+1 block.

Prove that every computational-basis configuration on a periodic ring has an even number of domain walls. Why does this constrain the h=0h=0 excitation energy?

Solution

Assign zj=±1z_j=\pm1 to site jj. A domain wall occurs when zjzj+1=−1z_jz_{j+1}=-1. Around a ring,

∏j=1Lzjzj+1=(∏j=1Lzj)2=1.\prod_{j=1}^{L} z_jz_{j+1} = \left( \prod_{j=1}^{L}z_j \right)^2 = 1.

If NDWN_{\mathrm{DW}} bonds have product −1-1, the same product is (−1)NDW(-1)^{N_{\mathrm{DW}}}. Hence NDWN_{\mathrm{DW}} is even. Each wall changes one bond energy from −J-J to +J+J and costs 2J2J, so the lowest domain-wall excitation of the periodic h=0h=0 spin chain contains two walls and costs 4J4J.

5. Extract the critical finite-size exponent

Section titled “5. Extract the critical finite-size exponent”

Starting from

ΔL=4Jsin⁡(π4L+2),\Delta_L = 4J \sin \left( \frac{\pi}{4L+2} \right),

derive its leading large-LL behavior and identify zz.

Solution

Let ℓ=L+12\ell=L+\tfrac12. Then

ΔL=4Jsin⁡(π4ℓ).\Delta_L = 4J \sin \left( \frac{\pi}{4\ell} \right).

Using sin⁡x=x−x3/6+O(x5)\sin x=x-x^3/6+O(x^5) gives

ΔL=πJℓ[1−π296ℓ2+O(ℓ−4)]=πJL+O(L−2).\begin{aligned} \Delta_L &= \frac{\pi J}{\ell} \left[ 1 - \frac{\pi^2}{96\ell^2} + O(\ell^{-4}) \right] \\ &= \frac{\pi J}{L} + O(L^{-2}). \end{aligned}

At a quantum critical point, a characteristic gap scales as ΔL∝L−z\Delta_L\propto L^{-z}. Therefore z=1z=1.

  1. E. Lieb, T. Schultz, and D. Mattis, “Two soluble models of an antiferromagnetic chain,” Annals of Physics 16, 407–466 (1961), doi:10.1016/0003-4916(61)90115-4.
  2. P. Pfeuty, “The one-dimensional Ising model with a transverse field,” Annals of Physics 57, 79–90 (1970), doi:10.1016/0003-4916(70)90270-8.
  3. E. Barouch and B. M. McCoy, “Statistical mechanics of the XY model. II. Spin-correlation functions,” Physical Review A 3, 786–804 (1971), doi:10.1103/PhysRevA.3.786.
  4. J. B. Kogut, “An introduction to lattice gauge theory and spin systems,” Reviews of Modern Physics 51, 659–713 (1979), doi:10.1103/RevModPhys.51.659.
  5. S. Sachdev, Quantum Phase Transitions, 2nd ed., Cambridge University Press (2011), doi:10.1017/CBO9780511973765.
  6. B. K. Chakrabarti, A. Dutta, and P. Sen, Quantum Ising Phases and Transitions in Transverse Ising Models, Springer (1996), doi:10.1007/978-3-540-49865-0.
  7. A. Dutta et al., Quantum Phase Transitions in Transverse Field Spin Models: From Statistical Physics to Quantum Information, Cambridge University Press (2015), doi:10.1017/CBO9781107706057.
  8. G. Vidal, J. I. Latorre, E. Rico, and A. Kitaev, “Entanglement in quantum critical phenomena,” Physical Review Letters 90, 227902 (2003), doi:10.1103/PhysRevLett.90.227902.
  9. P. Calabrese and J. Cardy, “Entanglement entropy and quantum field theory,” Journal of Statistical Mechanics P06002 (2004), doi:10.1088/1742-5468/2004/06/P06002.
  10. S. Suzuki, J.-i. Inoue, and B. K. Chakrabarti, Quantum Ising Phases and Transitions in Transverse Ising Models, 2nd ed., Springer (2013), doi:10.1007/978-3-642-33039-1.