Physics 239a:
Where do quantum field theories come from?
Lecturer: McGreevy

Last updated: 2024/05/17, 12:49:47

1

0.1

Introductory remarks . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

4

1 QFT from springs

7

2 Ising spins

16

2.1

Quantum-to-classical correspondence . . . . . . . . . . . . . . . . . . . . . .

17

2.2

Transverse-Field Ising Model . . . . . . . . . . . . . . . . . . . . . . . . . . .

32

2.2.1

Duality between spin flips and domain walls . . . . . . . . . . . . . .

38

2.2.2

Mean field theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

41

2.2.3

Correlation functions and long-range order . . . . . . . . . . . . . . .

43

2.2.4

Physics of structure factors. . . . . . . . . . . . . . . . . . . . . . .

44

2.2.5

Solution of Ising chain in terms of Majorana fermions . . . . . . . . .

48

2.2.6

Continuum limit . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

54

2.2.7

Duality and Jordan-Wigner with periodic boundary conditions . . . .

58

2.2.8

Coherent state path integrals for fermions . . . . . . . . . . . . . . .

67

2.2.9

Scaling near the fixed-point theory of the Ising phase transition. . . .

74

2.2.10 Beyond the quantum Ising chain . . . . . . . . . . . . . . . . . . . . .

76

3 Bosons and Bosonization

77

3.0.1

Spins, bosons, fermions . . . . . . . . . . . . . . . . . . . . . . . . . .

77

3.0.2

Using the JW solution to learn about the bosons . . . . . . . . . . .

79

3.0.3

Continuum description . . . . . . . . . . . . . . . . . . . . . . . . . .

81

3.0.4

Fermi liquid in d = 1 + 1. . . . . . . . . . . . . . . . . . . . . . . . .

82

2

3.0.5
3.1

Bosonization, part 1: counting . . . . . . . . . . . . . . . . . . . . . .

86

Free boson CFT in 1+1d, aka Luttinger liquid theory . . . . . . . . . . . . .

89

3.1.1

Application 1: Bosonization dictionary . . . . . . . . . . . . . . . . .

99

3.1.2

Application 2: Briefly, what is a Luttinger liquid? . . . . . . . . . . . 106

3.1.3

Application 3: sine-Gordon and Thirring . . . . . . . . . . . . . . . . 111

3.1.4

XY transition from SF to Mott insulator . . . . . . . . . . . . . . . . 113

4 CFT
4.1

4.2

The stress tensor and conformal invariance (abstract CFT) . . . . . . . . . . 118
4.1.1

Geometric interpretation of conformal group . . . . . . . . . . . . . . 121

4.1.2

Infinite conformal algebra in D = 2. . . . . . . . . . . . . . . . . . . . 122

Radial quantization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 124
4.2.1

4.3

118

The Virasoro central charge in 2d CFT . . . . . . . . . . . . . . . . . 128

Back to general dimensions. . . . . . . . . . . . . . . . . . . . . . . . . . . . 129
4.3.1

Constraints on correlation functions from CFT . . . . . . . . . . . . . 129

4.3.2

Thermodynamics of a CFT . . . . . . . . . . . . . . . . . . . . . . . 131

4.3.3

Too few words about the Conformal Bootstrap . . . . . . . . . . . . . 132

5 D =2+1

134

5.1

Duality for the transverse-field Ising model in D = 2 + 1 . . . . . . . . . . . 134

5.2

Z2 lattice gauge theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 136

5.3

Phase diagram of Z2 gauge theory and topological order . . . . . . . . . . . 137

5.4

Superfluid and Mott Insulator in D = 2 + 1 . . . . . . . . . . . . . . . . . . 142

5.5

Compact electrodynamics in D = 2 + 1 . . . . . . . . . . . . . . . . . . . . . 144
3

0.1

Introductory remarks

I begin with some discussion of my goals for this course. This is a special topics course
directed at graduate students in theoretical physics; this includes high-energy theory andcondensed matter theory and maybe some other areas, too.
The subject of the course is regulated quantum field theory (QFT): we will study quantum
field theories which can be constructed by starting from systems with finitely many degrees
of freedom per unit volume, with local interactions between them. Often these degrees of
freedom will live on a lattice and perhaps could arise as a model of a solid. So you could saythat the goal is to understand how QFT arises from condensed matter systems. (However,I do not plan to discuss any chemistry or materials science.)
Why do I make a big deal about ‘regulated’ ? For one thing, this is how it comes to us(when it does) in condensed matter physics. For another, such a description is required if wewant to know what we’re talking about. For example, we need it if we want to know what
we’re talking about well enough to explain it to a computer. Many QFT problems are too
hard for our brains. A related but less precise point is that I would like to do what I canto erase the problematic perspective on QFT which ‘begins from a classical lagrangian andquantizes it’ etc, and leads to a term like ‘anomaly’.
The crux of many problems in physics is the correct choice of variables with which to label
the degrees of freedom. Often the best choice is very different from the obvious choice; aname for this phenomenon is ‘duality’. We will study many examples of it (Kramers-Wannier,Jordan-Wigner, bosonization, Wegner, particle-vortex, perhaps others). This word is dangerous (it is one of the forbidden words on my blackboard) because it is about ambiguities
in our (physics) language. I would like to reclaim it.
An important bias in deciding what is meant by ‘correct’ or ‘best’ in the previous paragraphis: we will be interested in low-energy and long-wavelength physics, near the groundstate.
For one thing, this is the aspect of the present subject which is like ‘elementary particle
physics’; the high-energy physics of these systems is of a very different nature and bears little
resemblance to the field often called ‘high-energy physics’ (for example, there is volume-law
entanglement).
Also, notice that above I said ‘a finite number of degrees of freedom per unit volume’.
This last is important, because we are going to be interested in the thermodynamic limit.
Questions about a finite amount of stuff (this is sometimes called ‘mesoscopics’) tend to bemuch harder.
An important goal for the course is demonstrating that many fancy phenomena precious
to particle physicists can emerge from very humble origins in the kinds of (completely welldefined) local quantum lattice models we will study. Here I have in mind: fermions, gauge
theory, photons, anyons, strings, topological solitons, CFT, and many other sources of wonder
I’m forgetting right now.
4

Coarse classification of translation-invariant states of matter
So the most basic question we can ask is: how many degrees of freedom are there at thelowest energies? There are essentially three possibilities:
1. None.
2. Some.
3. A lot.
A more informative tour through that list goes like this. First let me make the assumptionthat the system has translation invariance, so we can label the excitations by momentum.
(Relaxing this assumption is interesting but we won’t do it this quarter.)
1. None: Such a system has an energy gap (‘is gapped’): the energy difference ∆E =E1 − E0 between the first excited state and the groundstate is nonzero, even in thethermodynamic limit. Note that ∆E is almost always nonzero in finite volume. (Recall,for example, the spectrum of the electromagnetic field in a box of linear size L: En ∼n.) The crucial thing here (in contrast to the case of photons) is that this energy stays
Lfinite even as L → ∞.
The excitations of such a system are massive particles1 .
2. Some: An example of what I mean by ‘some’ is that the system can have excitations
which are massless particles, like the photon.
The lowest
p energy degrees of freedom occur at isolated points in momentum space:ω(k) = c ~k · ~k vanishes at ~k = 0.
In this category I also put the gapless fluctuations at a critical point; in that case, it’s
not necessarily true that ω ∼ k integer and those excitations are not necessarily particles.
But they are still at k = 02 .
3. A lot: What I mean by this is Fermi surfaces (but importantly, not just free fermions).
Let’s go through that list one more time more slowly. Let’sreconsider the case of gapped systems. Different gapped states
are different if we can’t deform the hamiltonian to get from one
to the other without closing the gap. You might be botheredthat it is hard to imagine checking that there is no way aroundthe wall. It is therefore important to find sharp characterizations of such states, like integer labels, which cannot changesmoothly.
1
2

Verstraete et al claims to have proved a version of this statement.
or some other isolated points in momentum space.

5

Even the lowest-energy (even below the gap) physics of gapped systems can be deeply
fascinating. Such a thing is a topological field theory: it is a theory of groundstates, and it
can provide a way to distinguish states of matter. For example, it may be that the numberof groundstates depends on the topology of the space on which we put the system. Thisphenomenon is called topological order. Another possibility is that even if the system in
infinite space has an energy gap, if we cut the space open, new stuff can happen; for example
there may be gapless edge modes. Both of these phenomena happen in quantum Hall systems.
A reason to think that an interface between the vacuum and a gappedstate of matter which is distinct from the trivial one might carry gaplessmodes is that the couplings in the hamiltonian are forced to pass throughthe wall where the gap closes. (In fact there are important exceptions
to this conclusion.)
An energy gap (and no topological order or special edge modes) should
probably be the generic expectation for what happens if you pile together
a bunch of degrees of freedom and couple them in some haphazard (buttranslation invariant!) way. At the very least this follows on generalgrounds of pessimism: if you generically got something interesting by
doing this, physics would be a lot easier (or more likely: we wouldn’t find it interesting
anymore). Gaplessness is something special that needs to be explained. Here is alist of some possible reasons for gaplessness (if you find another, you should write it down):
1. broken continuous symmetry (Goldstone bosons)
2. tuning to a critical point – notice that this requires some agent to do the tuning, andwill only occur on some subspace of the space of couplings of nonzero codimension.
3. continuous unbroken gauge invariance (photons)
4. Fermi surface (basically only in this case do we get gapless degrees of freedom at somelocus of dimension greater than one in momentum space)
5. edge of topological phase: non-onsite realization of symmetry, anomaly inflow.
6. CFT with no relevant operators
Each entry in this list is something to be understood. We will get to many of them. I havepromised various people that I would devote special attention to the second option, in thecase of 1+1 dimensions; it is also relevant to the last two options.
Another important axis along which to organize states of matter is by symmetry. Specifically, we can label them according to the symmetry group G that acts on their Hilbert
space and commutes with the Hamiltonian. You will notice that here I am speaking about
what are called global symmetries, that is, symmetries (not redundancies of our labelling,
like gauge transformations).
6

There are many refinements of this organization. We can ask how the symmetry G isrealized, in at least three senses:
1. most simply, what representations of the group appear in the system?
2. is the symmetry preserved by the groundstate? If not, this is called ‘spontaneous
symmetry breaking’.
3. is it ‘on-site’ ? Alternatively, is it ‘anomalous’ ? What are its anomaly coefficients? I’ll
postpone the explanation of these terms.
Some useful lattice models (and their symmetry), which I hope to discuss this quarter:• toric code (none), lattice gauge theory (none, necessarily)• ising (Z2 )• potts, clock (SN , ZN )• XY (spins, rotors) (U(1))• Heisenberg (SU(2))• O(n) (O(n))• lattice (tight-binding) fermions• lattice (tight-binding) bosons• Kitaev honeycomb
Resources: The material which comes closest to my present goals is a course taught
by Ashvin Vishwanath at Berkeley called ‘Demystifying quantum field theory’ and many
insights can be found in his notes. Some material from the books by Sachdev, Wen, Fradkin,
Polyakov will be useful.

1

QFT from springs

My purpose here is to remind us that QFT need not be regarded as a mystical beast thatbegins its life with infinitely many degrees of freedom near every point in space. In addition, the following discussion contains rudimentary versions of many of the ideas we’ll need.
Actually it contains even more of them than I thought when I decided to include it – someof them are simply missing from the textbook treatments of this system.
7

Consider a one-dimensional box of point masses, each connected by a spring to its neighbor.
This is a nice simple model of a crystalline solid (we can speak about its shortcomings later).
When in equilibrium, the masses form a regular one-dimensional crystal lattice (equally
spaced mass points). Now let qn denote the displacement of the nth mass from its equilibrium
position xn and let pn be the corresponding momentum. Assume there are N masses andimpose periodic boundary conditions: qn+N = qn . The equilibrium positions themselves arexn = na, n = 1, 2...N
where a is the lattice spacing. The Hamiltonian for the collection of masses is:

N  2
X
1
pn
2
+ κ (qn − qn−1 ) + λq4 .
H=
2m
2
n=1

(1)

I’ve include a token anharmonic term λq4 to remind us that we are leaving stuff out; forexample we might worry whether we could use this model to describe melting. Now set λ = 0
because we are going to study small deviations from q = 0.
This hamiltonian above describes a collection of coupled oscillators, with a matrix of spring
constants V = kab qa qb . If we diagonalize the matrix of spring constants, we will have a
description in terms of decoupled oscillators, called normal modes.
Since our system has (discrete) translation invariance, these modes are labelled by awavenumber k 3 :NN
1 X ikxn
1 X ikxn
e qn , pk = √e pn ,
qk = √
N n=1
N n=1(Notice that in the previous expressions I didn’t use boldface; that’s because this step isreally just classical physics. Note the awkward (but in field theory, inevitable) fact that
we’ll have (field) momentum operators pk labelled by a wavenumber aka momentum.) The
nice thing about the fourier kernel is that it diagonalizes the translation operator:
Teikx ≡ eik(x+a) = eika eikx .

Regulators: Because N is finite, k takes discrete values (1 = eikN a ); this is a longwavelength “IR” property. Because of the lattice structure, k is periodic (only eikan , n ∈ Z
appears): k ≡ k + 2π/a; this is a short-distance “UV” property. The range of k can be taken
to be
2π(N − 1)
0≤k≤
.
Na
3

The inverse transformation is:
2π/a

2π/a

1 X −ikxn
pn = √
e
pk .
N k>0

1 X −ikxn
qn = √
e
qk ,
N k>0

8

Because of the periodicity in k, we can equivalently label the set of wavenumbers by:ππ2πor − < k ≤ .
aaaThis range of independent values of the wavenumber in a lattice model is called the Brillouinzone. There is some convention for choosing a fundamental domain which prefers the lastone but I haven’t found a reason to care about this.
0<k≤Summary: Because the system is in a box (periodic), k-space is discrete. Because thesystem is on a lattice, k-space is periodic. There are N oscillator modes altogether.
The whole hamiltonian is a bunch of decoupled oscillators, labelled by these funny wave
numbers:X  pk p−k 12+ mωk qk q−k
H=
2m
2
k
where the frequency of the mode labelled k isrκ
|k|a
ωk ≡ 2
sin
.
m
2

(2)

Why might we care about this frequency? For one thing, consider the Heisenberg equationof motion for the deviation of one spring:pni∂t qn = [qn , H] =, i∂t pn = [pn , H]m
Combining these gives:mq̈n = −κ ((qn − qn−1 ) − (qn − qn+1 )) = −κ (2qn − qn−1 − qn+1 ) .
In terms of the fourier-mode operators:
mq̈k = −κ (2 − 2 cos ka) qk .
P
Plugging in a fourier ansatz in time qk (t) = ω e−iωt qk,ω turns this into an algebraic equation
 2 |k|awhich says ω 2 = ωk2 = 2κsin 2 for the allowed modes. We see that (the classical versionmof) this system describes waves:


k1/a
0 = ω 2 − ωk2 qk,ω ' ω 2 − vs2 k 2 qk,ω .
The result for small k is the fourier transform of the wave equation:∂t2 − vs2 ∂x2 q(x, t) = 0 .

(3)

vs is the speed of propagation of the waves, in this case the speed of sound. Comparing to
the dispersion relation (2), we have found
r
∂ωk
κ
vs =
|k→0 = a.
∂km
9The description we are about to give is a quantization of sound waves.
Notice that when k = 0, ωk = 0. We are going to have to treat this mode specially; there
is a lot of physics in it.
So far the fact that quantumly [qn , pn0 ] = i~δnn0 1 hasn’t mattered in our analysis (go backand check). For the Fourier modes, this implies the commutator
XX[qk , pk0 ] =Ukn Uk0 n0 [qn , p0n ] = i~1
Ukn Uk0 n = i~δk,−k0 1.
n,n0n(In the previous expression I called Ukn = √1N eikxn the unitary matrix realizing the discrete
Fourier kernel.)
To make the final step to decouple the modes with k and −k, introduce the annihilationand creation operatorsrr

1 ~mωk ~ †
For k 6= 0:
qk =ak + a−k , pk =ak − a†−k .
2mωki
2They satisfy[ak , a†k0 ] = δkk0 1.
In terms of these, the hamiltonian is
H0 =

X


~ωk

1
a†k ak +

2

k


+

p20
2m

– it is a sum of decoupled oscillators, and a free particle describing the center-of-mass. It
is worth putting together the final relation between the ‘position operator’ and the phonon
annihilation and creation operators:r~ X 1  ikx
1−ikx †qn =e ak + e
ak + √ q0
(4)
√
2m k
ωk
N
and the corresponding relation for its canonical conjugate momentumr
1
1 ~m X √  ikx−ikx †
ωk e ak − e
ak + √ p0 .
pn =
i
2 k
N
Notice that these expressions are formally identical to the formulae in a QFT textbook
expressing a scalar field in terms of creation and annihilation operators.
[End of Lecture 1]
The groundstate is obtained from|0i ⊗ |p0 = 0i,where ak |0i = 0, ∀ k, and p0 |p0 i = p0 |p0 i .
10A set of excited states isa†k |0iThis state has energy ~ωk above the groundstate. In a box of size L = N a, the smallest
energy excitation has k1 = 2π
and energy
L
1 L→∞
→ 0.
(5)
L
(Note that here I am taking L → ∞ to implement the thermodynamic limit of infinitely
many degrees of freedom; the lattice spacing can remain finite for this purpose – it is not acontinuum limit.) So according to our definition, this system is gapless. Why? Goldstone4 :the system has a symmetry under qn → qn +  for all n. If everyone moves to the left threefeet, none of the springs are stretched. This is the dance enacted by the k = 0 mode. Ifnearly everyone moves nearly three feet to the left, the springs will only be stretched a little;
hence the modes with small k have small ω.
∆E ∼Tower of States: Now I will say a few words about the zeromode, which is horribly
mistreated in all the textbook discussions of this system that I’ve seen (Le Bellac, AltlandSimons, and my own quantum mechanics lecture notes!). There is no potential at all for thismode – it drops out of the (qn − qn+1 )2 terms. It just has a kinetic term, which we can thinkof as the center-of-mass energy of the system. How much energy does it cost to excite thismode? Notice that if everyone moves to the left by a, the system comes back to itself (I amassuming that the masses are indistinguishable particles): |{qn }i ' |{qn + a}i. In terms of
the k = 0 mode, this is
!
N
N
X
√
1
1 X
−i0xn
qn e
'√
qn + N a , i .e. q0 ' q0 + N a.
q0 = √
N n=1
N n=1
This means that the wavefunction for the zeromode must satisfy√
2πZ
eip0 q0 = eip0 (q0 + N a) =⇒ p0 ∈ √Naand the first excited state has energy
p20
1 1|p0 = √2π =
Na
2m
2 Nm



2π
a

2
.

This is a victory at least in the sense that we expect the center of mass of the system tohave an intertial mass N m. Notice that the spacing of these states depends differently on
the parameters than that of the ones from the nonzero-k phonon states.
But actually this phenomenon is ubiquitous: it happens whenever we take a system which
breaks a continous symmetry (here: a solid breaks continuous translation invariance)5 and4I should probably give a review of Goldstone’s theorem here. The relevant statement for our purposes
is: if the groundstate is not invariant under a continuous symmetry of H, the spectrum is gapless. A usefulreference for a more specific statement is this paper.
5
The fine print:

11

put it in finite volume, i.e. depart from the thermodynamic limit. In particular, in finite
volume the zeromode associated with a conserved quantity (here the momentum) produces
a tower of states with a different level-spacing (as a function of system size L = N a) thanthe usual modes (5). (It is sometimes called the Anderson Tower of States in the study ofmagnetism or the Rotator spectrum in lattice gauge theory). In this case, both towers go like
1/N , but this is a coincidence. In other cases the tower from the zeromode is more closely
1∼ L1d ∼ N1 ) than the particle momentum tower (which always goes
spaced (it goes like volume
1
like L1 ∼ N 1/d(or maybe L12 )), so the tower of states from the zeromode is usually much closer
together, and in the thermodynamic limit L → ∞, they combine to form the degeneratevacua associated with spontaneous symmetry breaking. 6
Phonons. Back to the phonon excitations. We can specify basis states for this Hilbertspace nk1  nk2
a†k1
a†k2
· · · |0i = |{nk1 , nk2 , ...}iby a collection of occupation numbers nk , eigenvalues of the number operator for each normalmode.
We can make a state with two phonons:|k, k 0 i = a†k a†k0 |0i
and so on. As you can see, phonons are identical bosons: the state with several of them only
says how many there are with each wavenumber.
So you see that we have constructed an approximation to the Fock space of a (massless)scalar field from a system with finitely many degrees of freedom per unit volume (here,
length), and in fact finitely many degrees of freedom altogether, since we kept the IR regulator
L finite. It is worth pausing to appreciate this: we’ve been forced to discover a frameworkfor quantum systems in which particles can be created and annihilated, very different fromthe old-fashioned point of view where we have a fixed Hilbert space for each particle.
Many aspects of the above discussion are special to the fact that our hamiltonian was
quadratic in the operators. Certainly our ability to completely solve the system is. Notice
1. Actually it’s important that the order parameter doesn’t commute with the Hamiltonian; the exception
is ferromagnets, where the order parameter is the total spin itself, which is a conserved quantity andtherefore can be nonzero even in finite volume. So the tower is collapsed at zero in that case.
2. Actually, a one-dimensional mattress of oscillators will not spontaneously break continuous translation
symmetry even in infinite volume. This is a consequence of the Coleman-Mermin-Wagner theorem:
the positions of the atoms still fluctuate too much, even when there are infinitely many of them in arow; more than one dimension is required to have the crystal really sit still. You’ll see the effects ofthese fluctuations on the problem set when you study the Debye-Waller factors. This does not vitiate
our conclusions above at all.
6

The definitive discussion of this subject can be found in the last few pages of P. Anderson, Concepts in
Solids.

12

that the number of phonons of each momentum nk ≡ a†k ak is conserved for each k. But if
we add generic cubic and quartic
P terms in q (or if we couple our atoms to7 the photon field)
even the number of phonons k nk will no longer be a conserved quantity. So a descriptionof such particles which forced us to fix their number wouldn’t be so great. More generally, I
would like to emphasize that not every QFT can be usefully considered as nearly a bunch of
harmonic oscillators. Finding ways to think about the ones which can’t is an important job.
A nice example where we can see the importance of the tower of states and of the quantization of phonon number is the Mössbauer effect: when scattering high-energy photons offa solid, there is a finite amplitude for scattering events which produce zero phonons. Thismeans that all the momentum transfer goes to the center of mass mode, which has negligible
recoil as N → ∞, since its inertial mass is N m. This allows for very sharp absorption
lines, which if the atom were in free space would be washed out (i.e. broadened to a width2γ).) by the random recoils of the individual atoms (as depicted in the comic strip
Erecoil = (p2m
below).

I’ve made a big deal about the regulators here. One reason we care about them is if weremove them (N → ∞, a → 0) and ask badP
questions, we’ll get infinity. For example, we
could think about the vacuum energy E0 = 21 k ~ωk . There is physics in there (e.g. Casimirforces), but we will not discuss that now.

7

Note that it is possible to make a quadratic action for conserved particles, but this requires adding moredegrees of freedom – the required U(1) symmetry must act something like(q 1 , q 2 ) → (cos θq 1 , sin θq 2 ).
We can reorganize this as a complex field Φ = q 1 + iq 2 on which the symmetry acts by Φ → eiθ Φ. This is
illustrated on the problem set.

13

Continuum limit
The thing we have discovered is a quantum field theory.
2p+ V (q),At this point I will use the path-integral description of a 1d particle with H = 2m
a basic statement of which is the following formula for the propagator

hq|e

−iHt

Z q(t)=q|q0 i =

Rt

1 2[dq]ei 0 dt ( 2 q̇ −V (q)) .

q(0)=q0

Q τ
Here [dq] ≡ N M
l=1 dq(tl ) – the path integral measure is defined by a limiting procedure,and N is a normalization factor that always drops out of physical quantities.
Recall that the key step in the derivation of this statement is the evaluation of the propagator for an infinitesimal time:p2

hq2 |e−iH∆t |q1 i = hq2 |e−i∆t 2m e−iV (q) |q2 i + O(∆t2 ) .
An integral expression for this can be obtained by inserting resolutions of the identity

 Z
Z
2
dq|qihq|dp|pihp|
1=1 =in between the two exponentials.
Scalar field theory in one dimension
If we use the path integral description, some of these things (in particular the continuum, sound-wave limit) are more obvious-seeming. The path integral for our collection ofoscillators isZ
Z = [dq1 · · · dqN ]eiS[q]

R
P 1
P 1
2
2
with S[q] = dt
m
q̇
−
V
({q})
. V ({q}) =nnn 2 κ (qn+1 − qn ) . Now let’s try to
n 2take the continuum limit a → 0, N → ∞. Basically the only thing we need is to think of

qn = q(x = na) as defining a smooth function:[Note that the continuum field is often called φ(x) instead of q(x) for some reason. At least
the letters q(x) and φ(x) look similar.]
We now have(qn − qn−1 )2 ' a2 (∂x q)2 |x=na .
14

Now the path integral becomes:ZZ=[Dq]eiS[q]with Dq now representing an integral over all configurations q(t, x) (defined by this limit)
and
Z
Z
Z
Z

1
2
2
2
2
4
S[q] = dt dx µ (∂t q) − µvs (∂x q) − rq − uq − ... ≡ dt dxL2where I’ve introduced some parameters µ, vs , r, u determined from m, κ... in some ways that
we needn’t worry about right now (e.g. µ = m/a, the mass per unit lengthR of the chain). L is
the Lagrangian density whose integral over space is the Lagrangian L = dxL. The ellipses... represent terms of higher order in the Taylor expansion of q(x), which are suppressed bycorrespondingly more powers of a: an ∂xn q.
The equation of motion (stationary phase condition) is0=

δS= −µq̈ − µvs2 ∂x2 q − rq − 2uq 3 − ...
δq(x, t)From the phonon problem, we automatically found r = u = 0, and the equation of motion
is just the wave equation (3). This happened because of the symmetry qn → qn + . This is
the operation that translates the whole crystal, It guarantees low-energy phonons near k = 0
because it means q(x) can only appear in S via its derivatives.
The expressions for the field operators in terms of mode operators are nows Z
1  ikx~dk √e ak + e−ikx a†kq(x) =2µωkand its canonical conjugate momentum
r Z

√ 
1 ~µ
p(x) =
dk ωk eikx ak − e−ikx a†k .
i2(p(x) is the quantum operator associated with the field-momentum π above.) These equations are the same as (4), but with squinting. They are just the ones in introductory QFT
textbooks; the stray factors of µ arise because we didn’t ‘canonically normalize’ our fields√
and absorb the µs into the field, e.g. defining φ ≡ µq would get rid of them. Notice thatthe position along the chain x here is just a label on the fields (it was n), not a quantum
operator.
The field q is called a scalar field because it doesn’t have any indices decorating it. This is tobe distinguished from e.g. the Maxwell field, which is a vector field. (Note that vibrations of acrystal in three dimensions actually do involve vector indices. We will omit this complicationfrom our discussion.)15

2

Ising spins

The Ising model has many guises. There is this from statistical mechanics:XPZ=
e−K hjli sj sl .
{sj }There is this quantum spin system:
HT F IM = −JX

(gx Xj + gz Zj Zj+1 ) .

j

And there is this 2d conformal field theory:ZS[χ] = d2 z χ∂ z χ(6)which I first encountered on the worldsheet of a superstring. An important part of our job
is to understand the connections between these things. One thing they have in common is aZ2 symmetry, sj → −sj or Zj → −Zj or χ → −χ.
Let me say a few introductory words about quantum spin systems. This means that wehave a collection of two-state systems (aka qbits) Hj = span{| ↑j i, | ↓j i} distributed overspace and coupled somehow:O
H=
Hj , dim (H) = 2Njwhere N is the number of sites.
One qbit: To begin, consider just one two-state system. There are four independent hermitian operators acting on this Hilbert space. Besides the identity, there are the three Paulis,
which I will denote by X, Y, Z instead of σ x , σ y , σ z :






0 1
0 −i
1 0
x
y
z
X≡σ =
, Y≡σ =, Z≡σ =
1 0
i 0
0 −1
This notation (which comes to us from the quantum information community) makes theimportant information larger and is therefore better, especially for those of us with limited
eyesight.
[End of Lecture 2]
They satisfy
XY = iZ, XZ = −ZX, X2 = 1,
and all cyclic permutations X → Y → Z → X of these statements.
Multiple qbits: If we have more than one site, the paulis on different sites commute:[σj , σl ] = 0, j 6= li .e. Xj Zl = (−1)δjl Zl Xj ,where σj is any of the three paulis acting on Hj .
16

2.1

Quantum-to-classical correspondence[Kogut, Sachdev chapter 5, Goldenfeld §3.2]
Goals : Facilitate rapid passage between Lagrangian and Hamiltonian descriptions. Harness
your knowledge of stat mech.
1.
Important conclusions: correlation length ∼ gap

Let’s begin with the classical ising model in a (longitudinal) magnetic field:
XPPZ=
e−K hjli sj sl −h j sj .

(7)

{sj }

Here I am imagining we have classical spins sj = ±1 at each site of some graph, and hjli
denotes pairs of sites which share a link in the graph. You might be tempted to call K theinverse temperature, which is how we would interpret if we were doing classical stat mech;
resist the temptation.

17

One qbit from ising chain
First, let’s think about the case when the graph in (7) is just a chain:
X

Z1 =

−Hc

e

Hc = −K,{sl =±1}MτX

sl sl+1 − h

Mτ
X

l=1

sl

(8)

l=1

These ss are now just Mτ numbers, each ±1 – there are 2Mτ terms in thissum. (Notice that the field h breaks the s → −s symmetry of the summand.) The parameter
K > 0 is the ‘inverse temperature’ in the Boltzmann distribution; I put these words in quotes
because I want you to think of it as merely a parameter in the classical hamiltonian.
For definiteness let’s suppose the chain loops back on itself,
sl+Mτ = sl (periodic boundary conditions).
P
Q
Using the identity e l (...)l = l e(...)l ,
Z1 =

Mτ
XY

T1 (sl , sl+1 )T2 (sl ){sl } l=1

where
T1 (s1 , s2 ) ≡ eKs1 s2 ,T2 (s) ≡ ehs .

What are these objects? The conceptual leap is to think of T1 (s1 , s2 ) as a 2 × 2 matrix:
 K

e
e−K
T1 (s1 , s2 ) = −K
= hs1 |T1 |s2 i,eeK s1 s2which we can then regard as matrix elements of an operator T1 acting on a 2-state quantum
system (hence the boldface). And we have to think of T2 (s) as the diagonal elements of thesame kind of matrix:
 h

e
0
δs1 ,s2 T2 (s1 ) == hs1 |T2 |s2 i.
0 e−h s1 s2
So we have

Z1 = tr (T1 T2 )(T1 T2 ) · · · (T1 T2 ) = trTMτ
{z
}
|

(9)

Mτ times

where I’ve written

11

T ≡ T22 T1 T22 = T† = Tt
for convenience (so it’s symmetric). This object is the transfer matrix. What’s the traceover in (9)? It’s a single two-state system – a single qbit (or quantum spin) that we’veconstructed from this chain of classical two-valued variables.
18

Even if we didn’t care about quantum spins, this way of organizing the partition sum of
the Ising chain does the sum for us (since the trace is basis-independent, and so we might
as well evaluate it in the basis where T is diagonal):Mττ
Z1 = trTMτ = λM
+ + λ−

where λ± are the two eigenvalues of the transfer matrix, λ+ ≥ λ− :(p
2 cosh K
h→0
λ± = eK cosh h ± e2K sinh2 h + e−2K →2 sinh K.

(10)

In the thermodynamic limit, Mτ  1, one of them dominates the free energy Mτ !
λ−
τ
τ
∼ λM
1+
e−F = Z1 = λM
+
+ .
λ+

Now I command you to think of the transfer matrix as
T = e−∆τ H
the propagator in euclidean time (by an amount ∆τ ), where H is the quantum hamiltonian
operator for a single qbit (note the boldface to denote quantum operators). So what’s H?
To answer this, let’s rewrite the parts of the transfer matrix in terms of paulis, thinking ofs = ± as Z-eigenstates. For T2 , which is diagonal in the Z basis, this is easy:
T2 = ehZ .
To write T1 this way, stare at its matrix elements in the Z basis:

 K
e
e−K
hs1 |T1 |s2 i = −KeeK s1 s2and compare them to those of
eaX+b1 = eb eaX = eb (cosh a + X sinh a)which areaX+b1

hs1 |eb|s2 i = ecosh a sinh asinh a cosh a
s1 ,s2

So we want eb sinh a = e−K , eb cosh a = eK which is solved by
e−2K = tanh a .
So we want to identify
T1 T2 = eb1 +aX ehZ ≡ e−∆τ H
19

(11)

for small ∆τ . This requires that a, b, h scale like ∆τ , and so we can combine the exponents.
Assuming that ∆τ  1 (so that we can ignore the corrections to the leading term in the
Baker-Campbell-Hausdorff formula), the result is
H = E0 −∆X − h̄Z .
2bh2aHere E0 = ∆τ, h̄ = ∆τ, ∆ = ∆τ. (Note that it’s not surprising that the Hamiltonian for~ , since these operators span the set of
an isolated qbit is of the form H = d0 1 + d~ · σhermitian operators on a qbit; but the relation between the parameters that we’ve foundwill be important.)

To recap, let’s go backwards: consider the quantum system consisting of a single spin with
H = E0 − ∆2 X + h̄Z . Set h̄ = 0 for a moment. Then ∆ is the energy gap between the
groundstate and the first excited state (hence the name). The thermal partition function isX
hs|e−βH |si,
(12)
ZQ (T ) = tre−H/T =s=±in the Z basis, Z|si = s|si. I emphasize that T here is the temperature to which we aresubjecting our quantum spin; β = T1 is the length of the euclidean time circle. Break up the
euclidean time circle into Mτ intervals of size ∆τ = β/Mτ . Insert many resolutions of unity(this is called ‘Trotter decomposition’)X
hsMτ |e−∆τ H |sMτ −1 ihsMτ −1 |e−∆τ H |sMτ −2 i · · · hs1 |e−∆τ H |sMτ i .
ZQ =s1 ...sMτ

The RHS8 is the partition function of a classical Ising chain, Z1 in (8), with h = 0 and K
given by (11), which in the present variables is:


β∆
−2K
e
= tanh.
(13)
2Mτ
Notice that if our interest is in the quantum model with couplings E0 , ∆, we can use any
Mτ we want – there are many classical models we could use9 . For given Mτ , the couplings
we should choose are related by (13).
A quantum system with just a single spin (for any H not proportional to 1) clearly has aunique groundstate; this statement means the absence of a phase transition in the 1d Ising
chain.
8
9

RHS = right-hand side, LHS = left-hand side, BHS = both-hand side
If we include the Z term, we need to take ∆τ small enough so that we can writee−∆τ H = e∆τ 2 X e−∆τ (E0 −h̄Z) + O(∆τ 2 )∆20

More than one spin10
Let’s do this procedure again, supposing the graph in question is a cubic lattice with morethan one dimension, and let’s think of one of the directions as euclidean time, τ . We’ll end
up with more than one spin.
We’re going to rewrite the sum in (7) as a sum of products of (transfer) matrices. I will draw the pictures associated to a square lattice, but this is not a crucial limitation. Label points on the lattice by a vector ~n of integers; a unit vector in the time direction is. FirstP τ̌−Hrewrite the classical action Hc in Zc =e c as11X
Hc = −(Ks(~n + τ̌ )s(~n) + Kx s(~n + x̌)s(~n))
~
n

X
X 1
2
s(~n + x̌)s(~n)(s(~n + τ̌ ) − s(~n)) − 1 − Kx
= K
2
~
n
~
n
X
L(l + 1, l)= const +(14)

rows at fixed time, lwith12X
1 X
1
L(s, σ) = K(s(j) − σ(j))2 − Kx(s(j + 1)s(j) + σ(j + 1)σ(j)) .
22jjσ and s are the names for the spins on successive time slices, as in the figure at left.
The transfer matrix between successive time slices is a 2M ×
2 matrix:
hs|T|σi = Tsσ = e−L(s,σ) ,Min terms of which
Z=

X
{s}

−Hc

e

=

Mτ
XY

Ts(l,j),s(l+1,j) = trH TMτ .

{s} l=1

This is just as in the one-site case; the difference is that now the hilbert space has a twostate system for every site on a fixed-l slice of the lattice. I will call this “space”, and label
these sites by an index j. (NoteN that nothing we say in this discussion requires space to be
one-dimensional.) So H = j Hj , where each Hj is a two-state system.
10

This discussion comes from this paper of Fradkin and Susskind, and can be found in Kogut’s review
article.
11
The minus signs in red were flubbed in lecture.
R
12
Note that ‘L’ is for ‘Lagrangian’. In retrospect, I should have called Hc ≡ S = dτ L, as Kogut does.
‘S’ is for ‘action’.

21

The diagonal entries of Ts,σ come from contributions where s(l) = σ(l): they come with afactor of Ts=σ = eL(0 flips) withX
L(0 flips) = −Kx
σ(j + 1)σ(j).
jThe one-off-the-diagonal terms come fromσ(j) = s(j), except for one site where instead σ(j) = −s(j).
This gives a contributionX
1
1
L(1 flips) = K(1 − (−1))2 − Kx(σ(j + 1)σ(j) + s(j + 1)s(j)) .
|2
{z
} 2
j
=2K

Similarly,X
1
L(n flips) = 2nK − Kx(σ(j + 1)σ(j) + s(j + 1)s(j)) .
2
j
Now we need to figure out who is H, as defined by
T = e−∆τ H ' 1 − ∆τ H ;we want to consider ∆τ small and must choose Kx , K to make it so. We have to match thematrix elements hs|T|σi = Tsσ :P

T (0 flips)sσ= δsσ PeKx j s(j)s(j+1)' 1 − ∆τ H|0 flips−2K 12 Kx j (σ(j+1)σ(j)+s(j+1)s(j))
T (1 flip)sσ = ee' −∆τ H|1 flip
P
−2nK 21 Kx j (σ(j+1)σ(j)+s(j+1)s(j))' −∆τ H|n flips
T (n flips)sσ = ee(15)From the first line, we learn that Kx ∼ ∆τ ; from the second we learn e−2K ∼ ∆τ ; we’ll call theratio which we’ll keep finite g ≡ Kx−1 e−2K . To make τ continuous, we take K → ∞, Kx → 0,
holding g fixed. Then we see that the n-flip matrix elements go like e−nK ∼ (∆τ )n and canbe ignored – the hamlitonian only has 0- and 1-flip terms.
To reproduce (15), we must take!
HTFIM = −J

g

X

Xj +

j

X

Zj+1 Zj.

j

Here J is a constant with dimensions of energy that we pull out of ∆τ . The first termis the ‘one-flip’ term; the second is the ‘zero-flips’ term. The first term is a ‘transverse
magnetic field’ in the sense that it is transverse to the axis along which the neighboring
spins interact. So this is called the transverse field ising model. We’re going to understandit completely below. As we’ll see, it contains the universal physics of the 2d Ising model,
including Onsager’s solution. The word ‘universal’ requires some discussion.
[End of Lecture 3]
22Symmetry of the
Q transverse field quantum Ising model: HTFIM has a Z2 symmetry,generated by S = j Xj , which acts by
SZj = −Zj S, SXj = +Xj S, ∀j;
On Z eigenstates it acts as:
S|{sj }j i = |{−sj }j i .
It is a symmetry in the sense that:[HTFIM , S] = 0.
Notice that S2 =

Q

2
†
−1
j Xj = 1,, and S = S = S .

By ‘a Z2 symmetry,’ I mean that the symmetry group consists of two elements G = {1, S},
and they satisfy S2 = 1, just like the group {1, −1} under multiplication. This group is
G = Z2 . (For a bit of context, the group ZN is realized by the N th roots of unity, undermultiplication.)
The existence of this symmetry of the quantum model is a direct consequence of the fact thatthe summand of the classical system was invariant under the operation sj → −sj , ∀j. This
meant that the matrix elements of the transfer matrix satisfy Ts,s0 = T−s,−s0 which implies
the symmetry of H. (Note that symmetries of the classical action do not so immediately
imply symmetries of the associated quantum system if the system is not as well-regulated
as ours is. This is the phenomenon called ‘anomaly’.)23Quantum Ising in d space dimensions to classical ising in d + 1 dims
[Sachdev, 2d ed p. 75] Just to make sure it’s nailed down, let’s go backwards again. The
partition function of the quantum Ising model at temperature T is
1

ZQ (T ) = trNM
e− T HI = tr e−∆τ HI
j=1 HjMτ

The transfer matrix here e−∆τ HI is a 2M ×2M matrix. We’re going to take ∆τ → 0, Mτ → ∞,holding T1 = ∆Mτ fixed. Let’s use the usual13 ‘split-step’ trick of breaking up the noncommuting parts of H:e−∆τ HI ≡ Tx Tz + O(∆τ 2 ).
Tx ≡ eJg∆τ

P

j Xj

Tz ≡ eJ∆τ,P

j Zj Zj+1

Now insert a resolution of the identity in the Z-basis,X|{sj }ih{sj }|, Zj |{sj }i = sj |{sj }i,1=.

sj = ±1.

{sj }M
j=1

many many times, one between each pair of transfer operators; this turns the transfer operators into transfer matrices. The Tz bit is diagonal, by design:
Tz |{sj }i = eJ∆τP

j sj sj+1|{sj }i.

The Tx bit is off-diagonal, but only on a single spin at a time:Yh{s0j }|Tx |{sj }i =hs0j |eJg∆τ Xj |sj i
{z
}
|
j
2×2

Acting on a single spin at site j, this 2 × 2 matrix is just the one from the previous discussion
around (11):0hs0j |eJg∆τ Xj |sj i = e−b eKsj sj , e−2b =1sinh (2Jg∆τ ) , e−2K = tanh (Jg∆τ ) .
2

Notice that it wasn’t important to restrict to 1 + 1 dimensions here. The only difference isin the Tz bit, which gets replaced by a product over all neighbors in higher dimensions:h{s0j }|Tz |{sj }i = δs,s0 eJ∆τP

hjli sj slwhere hjli denotes nearest neighbors, and the innocent-looking δs,s0 sets the spins sj = s0j
equal for all sites.
13

By ‘usual’ I mean that this is just like in the path integral of a 1d particle, when we write∆τ2e−∆τ H = e− 2m p e−∆τ V (q) + O(∆τ 2 ).

24

Label the time slices by a variable l = 1...Mτ .
Z = tre− T1 HI=MτX Y

h{sj (l + 1)}|Tz Tx |{sj (l)}i{sj (l)} l=1

The sum on the RHS runs over the 2M Mτ values of sj (l) = ±1, which is the right set of things
to sum over in the d + 1-dimensional classical ising model. The weight in the partition sum
is



X

X 
τ
exp
Z = e|−bM
J∆τ sj (l)sj+1 (l) + Ksj (l)sj (l + 1)
{z }
{z
} |
{z
}
|
unimportant{sj (l)}j,lj,lspace deriv, from Tztime deriv, from Tx

constant=X

e−Hclassical ising

spins

except that the the couplings are a bit anisotropic: the couplings in the ‘space’ direction
Kx = J∆τ are not the same as the couplings in the ‘time’ direction, which satisfy e−2K =tanh (Jg∆τ ). (At the critical point K = Kc , this can be absorbed in a rescaling of spatial
directions, as we’ll see later.)25

Dictionary. So this establishes a mapping between classical systems in d + 1 dimensions
and quantum systems in d space dimensions. Here’s the dictionary:statistical mechanics in d + 1 dimensions

quantum system in d space dimensions

transfer matrix

euclidean-time propagator, e−∆τ H

statistical ‘temperature’(lattice-scale) coupling K

free energy in infinite volume

groundstate energy: e−F = Z = tre−βH → e−βE0periodicity of euclidean time Lτtemperature: β = T1 = ∆τ Mτ

statistical averages

groundstate expectation valuesof time-ordered operators

β→0Note that this correspondence between classical and quantum systems is not an isomorphism.
For one thing, we’ve seen that many classical systems are related to the same quantum
system, which does not care about the lattice spacing in time. There is a set of physicalquantities which agree between these different classical systems, called universal, which is
the information in the quantum system. More on this below.
Consequences for phase transitions and quantum phase transitions.
One immediate consequence is the following. Think about what happens at a phase transition of the classical problem. This means that the free energy F (K, ...) has some kind of
singularity at some value of the parameters, let’s suppose it’s the statistical temperature,i.e. the parameter we’ve been calling K. ‘Singularity’ means breakdown of the Taylor expansion, i.e. a disagreement between the actual behavior of the function and its Taylor series –a non-analyticity. First, this can only happen in the thermodynamic limit (at the very leastMτ → ∞), since otherwise there are only a finite number of terms in the partition sum andF is an analytic function of K (it’s a polynomial in e−K ).
An important dichotomy is between continuous phase transitions (also called second orderor higher) and first-order phase transitions; at the latter, ∂K F is discontinous at the transition, at the former it is not. This seems at first like an innocuous distinction, but think
about it from the point of view of the transfer matrix for a moment. In the thermodynamic
limit, Z = λ1 (K)Mτ , where λ1 (K) is the largest eigenvalue of T(K). How can this have asingularity in K? There are two possibilities:
1. λ1 (K) is itself a singular function of K. How can this happen? One way it can happen
is if there is a level-crossing where two completely unrelated eigenvectors switch whichis the smallest (while remaining separated from all the others).
26

This is a first-order transition. A distinctive feature of a first order transition is alatent heat: although the free energies of the two phases are equal at the transition(they have to be in order to exchange dominance there), their entropies (and hence
energies) are not: S ∝ ∂K F jumps across the transition.
2. The other possibility is that the eigenvalues of T have an accumulation point at K =Kc , so that we can no longer ignore the contributions from the other eigenvalues totrTMτ , even when Mτ = ∞. This is the exciting case of a continuous phase transition.
In this case the critical point Kc is really special.

Now translate those statements into statements about the corresponding quantum system.
Recall that T = e−∆τ H – eigenvectors of T are eigenvectors of H! Their eigenvalues are
related by
λa = e−∆τ Ea ,so the largest eigenvalue of the transfer matrix corresponds to the smallest eigenvalue of H:
the groundstate. The two cases described above are:
1. As the parameter in H varies, two completely orthogonal states switch which one is
the groundstate. This is a ‘first-order quantum phase transition’, but that name is a
bit grandiose for this boring phenomenon, because the states on the two sides of thetransition don’t need to know anything about each other, and there is no interesting
critical theory. For example, the third excited state need know nothing about thetransition.
2. At a continuous transition in F (K), the spectrum of T piles up at the top. This meansthat the spectrum of H is piling up at the bottom: the gap is closing. There is agapless state which describes the physics in a whole neighborhood of the critical point.

Using the quantum-to-classical dictionary, the groundstate energy of the TFIM at thetransition reproduces Onsager’s tour-de-force free energy calculation.
Another failure mode of this correspondence: there are some quantum systems which when
Trotterized produce a stat mech model with non-positive Boltzmann weights, i.e. e−Hc < 027for some configurations; this requires Hc to be complex. These models are less familiar! (Anexample where this happens is the spin- 12 chain.) The quantum phase transitions of such
quantum systems are not just ordinary finite-temperature transitions of familiar classical statmech systems. So for the collector of QFTs, there is something to be gained by studying
quantum phase transitions.
Correlation functions[Sachdev, 2d ed p. 69] First let’s construct correlation functions of spins in the classical
Ising chain, (8), using the transfer matrix. (We’ll study correlation functions in the TFIMin §2.2.3.) Let
1 X −Hce sl sl 0C(l, l0 ) ≡ hsl sl0 i =Z1{sl }l

By translation invariance, this is only a function of the difference C(l, l0 ) = C(l − l0 ). For
simplicity, set the external field h = 0. Also, assume that l0 > l (as we’ll see, this is timeordering of the correlation function). In terms of the transfer matrix, it is:
1  Mτ −l0
l0 −l
l
0
C(l − l ) = tr T
ZT ZT .
(16)
Z
Notice that there is only one operator Z = σ z here; it is the matrix
Zss0 = δss0 s .
All the information about the index l, l0 is encoded in the location in the trace.
Let’s evaluate this trace in the basis of T eigenstates. When h = 0, we have T = eK 1 +e−K X, so these are X eigenstates:
T| →i = λ+ | →i,

T| ←i = λ− | →i .

Here | →i ≡ √12 (| ↑i + | ↓i).
In this basis


0 1hα|Z|βi =, α, β =→ or ← .
1 0 αβSo the trace (aka path integral) has two terms: one where the system spends l0 − l steps inthe state | →i (and the rest in | ←i), and one where it spends l0 − l steps in the state | →i.
The result (if we take Mτ → ∞ holding fixed l0 − l) is
0

0

0

0

τ −l +l l −l
λ+ Mτ →∞λMτ −l +l λl−−l + λM
0−→ tanhl −l K .
C(l − l) = +MτMτ
λ+ + λ−0

You should think of the insertions as
sl = Z(τ ), τ = ∆τ l.
28

(17)

So what we’ve just computed isC(τ ) = hZ(τ )Z(0)i = tanhl K = e−|τ |/ξ(18)

where the correlation time ξ satisfies
11=
ln coth K .
ξ
∆τ

(19)

Notice that this is the same as our formula for the gap, ∆, in (13).14 This connection between
the correlation length in euclidean time and the energy gap is general and important.
For large K, ξ is much bigger than the lattice spacing:ξ K1 1 2K
' e  1.
∆τ
2This is the limit we had to take to make the euclidean time continuous.
Notice that if we had taken l < l0 instead, we would have found the same answer with l0 − lreplaced by l − l0 .
[End of Lecture 4]

14

Seeing this requires the following cool hyperbolic trig fact:
If e−2K = tanh X then e−2X = tanh K(20)(i.e. this equation is ‘self-dual’) which follows from algebra. Here (13) says X = T ∆Mτ = ∆τ ∆ while (19)says X = ∆τ /ξ. Actually this relation (20) can be made manifestly symmetric by writing it as
1 = sinh 2X sinh 2K .
(You may notice that this is the same combination that appears in the Kramers-Wannier self-duality condition.) I don’t know a slick way to show this, but if you just solve this quadratic equation for e−2K and boil
it enough, you’ll find tanh X.

29

Continuum scaling limit and universality[Sachdev, 2d ed §5.5.1, 5.5.2] Now we are going to grapple with the term ‘universal’. Let’s
think about the Ising chain some more. We’ll regard Mτ ∆τ as a physical quantity, the proper
length of the chain. We’d like to take a continuum limit, where Mτ → ∞ or ∆τ → 0 or
maybe both. Such a limit is useful if ξ  ∆τ . This decides how we should scale K, h in thelimit. More explicitly, here is the prescription: Hold fixed physical quantities (i.e. eliminate
the quantities on the RHS of these expressions in favor of those on the LHS):
the correlation length,
the length of the chain,
physical separations between operators,
the applied field in the quantum system,

ξ
Lτ
τ
h̄

1
' ∆τ e2K ,2= ∆τ Mτ ,= (l − l0 )∆τ,= h/∆τ.

(21)

while taking ∆τ → 0, K → ∞, Mτ → ∞.
What physics of the various chains will agree? Certainly only quantities that don’t depend
explicitly on the lattice spacing; such quantities are called universal.
15
Consider the thermal freep energy of the single quantum spin (12) : The energy spectrumof our spin is E± = E0 ± (∆/2)2 + h̄2 , which means


 q
2
2
F = −T log ZQ = E0 − T ln 2 cosh β (∆/2) + h̄(just evaluate the trace in the energy eigenbasis). In fact, this is just the behavior of theising chain partition function in the scaling limit (21), since, in the limit (10) becomes
r


q
2ξ
∆τ
λ± '
1 + 4h̄2 ξ 21±∆τ2ξ
and so in the scaling limit (21)
F ' Lτ 


K
−
∆τ}
| {z




q

1
Lτ
−2
2
− ln 2 cosh
ξ + 4h̄  ,Lτ
2cutoff-dependent vac. energy

which is the same (up to an additive constant) as the quantum formula under the previouslymade identifications T = L1τ , ξ −1 = ∆.
We can also use the quantum system to compute the correlation functions of the classicalchain in the scaling limit (17). They are time-ordered correlation functions:C(τ1 − τ2 ) = ZQ−1 tre−βH (θ(τ1 − τ2 )Z(τ1 )Z(τ2 ) + θ(τ2 − τ1 )Z(τ2 )Z(τ1 ))15[Sachdev, 1st ed p. 19, 2d ed p. 73]

30

whereZ(τ ) ≡ eHτ Ze−Hτ .
This time-ordering is just the fact that we had to decide whether l0 or l was bigger in (16).
For example, consider what happens to this when T → 0. Then (inserting 1 =in an energy eigenbasis H|ni = En |ni),X
C(τ )|T =0 =|h0|Z|ni|2 e−(En −E0 )|τ |P

n |nihn|,nwhere the |τ | is taking care of the time-ordering. This is a spectral representation of thecorrelator. For large τ , the contribution of |ni is exponentially suppressed by its energy, so
the sum is approximated well by the lowest energy state for which the matrix element isnonzero. Assuming this is the first excited state (which in our two-state system it has no
choice!), we haveτ →∞C(τ )|T =0 ' e−τ /ξ , ξ = 1/∆,
where ∆ is the energy gap.
In these senses, the quantum theory of a single qbit is the universal theory of the Ising
chain. For example, if we
P began with a chain that had in addition next-nearest-neighbor0interactions, ∆Hc = K
j s(j)s(j + 2), we could redo the procedure above. The scaling
limit would not be exactly the same; we would have to scale K 0 somehow (it would also
have to grow in the limit). But we would find the same 2-state quantum system, and when
expressed in terms of physical variables, the ∆τ -independent terms in F would be identical,
as would the form of the correlation functions, which ise−|τ |/ξ + e−(Lτ −|τ |)/ξ.
C(τ ) = hZ(τ )Z(0)i =1 + e−Lτ /ξ
(Note that in this expression we did not assume |τ |  Lτ as we did before in (18), to which
this reduces in that limit.)

31

2.2

Transverse-Field Ising Model

Whether or not you liked the derivation above of its relation to the euclidean statisticalmechanics Ising model, we are going to study the quantum system whose hamiltonian isX
X
HTFIM = −J g
Xj +Zj Zl  .
(22)
j

hjli

Some of the things we say next will be true in one or more spatial dimensions.
Notice that J has units of energy; we could choose units where it’s 1. In 1d (or on bipartite
lattices), the sign of J does not matter for determining what state of matter we realize: if
J < 0, we can relabel our operators: Z̃j = (−1)j Zj and turn an antiferromagnetic interaction
into a ferromagnetic one. So let’s assume g, J > 0.
This model is interesting because of the competition between the two terms: the Xj term
wants each spin (independently of any others) to be in the state | →ij which satisfies
1
Xj | →ij = | →ij . | →ij = √ (| ↑ij + | ↓ij ) .
2
In conflict with this are the desires of −Zj Zj+1 , which is made happy (i.e. smaller) by themore cooperative states | ↑j ↑j+1 i, or | ↓j ↓j+1 i. In fact, it would be just as happy about
any linear combination of these a| ↑j ↑j+1 i + b| ↓j ↓j+1 i and we’ll come back to this point.
Another model which looks like it might have some form of competition isXXπ
Hboring = cos θ
Zj + sin θ
Xj , θ ∈ [0, ]2jj
Why is this one boring? Notice that we can continuously interpolate between the states
enjoyed by these two terms: the groundstate of H1 = cos θZ + sin θX isθθ|θi = cos | ↑i + sin | ↓i
22
– as we vary θ from 0 to π/2 we just smoothly rotate from | ↑z i to | ↑x i.
How do we know the same thing can’t happen in the transverse-field Ising chain? Symmetry.
We’ve already seen that the Ising model has a G = Z2 symmetry which acts by Zj → SZj S† =−Zj , Xj → SXj S† = +Xj , where the unitary S commutes with HTFIM : SHTFIM S† =HT F IM . The difference with Hboring is that HTFIM has two phases in which G is realized
differently on the groundstate.
g = ∞ : First, let’s take g so big that we may ignore the ZZ ferromagnetic term, soX
Hg→∞ = −
Xj .
j

32

(The basic idea of this discussion will apply in any dimension, on any lattice.) Since allterms commute, the groundstate is the simultaneous groundstate of each term:
Xj |gsi = +|gsi, ∀j,=⇒ |gsi = ⊗j | →ij .

Notice that this state preserves the symmetry in the sense that S|gsi = |gsi. Such asymmetry-preserving groundstate is called a paramagnet.
g = 0 : Begin with g = 0.
H0 = −JX

Zj Zj+1

jhas groundstates|+i ≡ | ↑↑ · · · ↑i, |−i ≡ | ↓↓ · · · ↓i,or any linear combination. Note that the states |±i are not symmetric: S|±i = |∓i, and so
we are tempted to declare that the symmetry is broken by the groundstate.

You will notice, however, that the states|

1(|+i ± |−i)2±i ≡ √

are symmetric – they are S eigenstates, so S maps them to themselves up to a phase. It gets
worse: In fact, in finite volume (finite number of sites of our chain), with g 6= 0, |+i and |−iare not eigenstates, and |+ i is the groundstate. BUT:
1. The two states |+i and |−i only mix at order N in perturbation theory P
in g, since wehave to flip all N spins using the perturbing hamiltonian ∆H = −gJ j Xj to getfrom one to the other. The tunneling amplitude is thereforeN →∞T ∼ g N h − |X1 X2 · · · XN |+i → 0.
2. There’s a reason for the symbol I used to denote the symmetric states: at large N ,these ‘cat states’ are superpositions of macroscopically distinct quantum states. Suchthings don’t happen, because of decoherence: if even a single dust particle in the roommeasures the spin of a single one of the spins, it measures the value of the whole chain.
In general, this happens very rapidly.
P
3. Imagine we add a small symmetry-breaking perturbation: ∆H = − j hZj ; this splits
the degeneracy between |+i and |−i. If h > 0, |+i is for sure the groundstate. Consider
preparing the system with a tiny h > 0 and then setting h = 0 after it settles down.
If we do this to a finite system, N < ∞, it will be in an excited state of the h = 0
Hamiltonian, since |+i will not be stationary (it will have a nonzero amplitude to33

tunnel into |−i). But if we take the thermodynamic limit before taking h → 0, it will
stay in the state we put it in with the ‘training field’ h. So beware that there is asingularity of our expressions (with physical significance) that means that the limits
do not commute:
lim lim Z 6= lim lim Z.
N →∞ h→0h→0 N →∞The physical one is to take the thermodynamic limit first.
The conclusion of this brief discussion is that spontaneous symmetry breaking actually
happens in the N → ∞ limit. At finite N , |+i and |−i are approximate eigenstates whichbecome a better approximation as N → ∞.
This state of a Z2 -symmetric system which spontaneously breaks the Z2 symmetry is calleda ferromagnet.
So the crucial idea I want to convey here is that there must
be a sharp phase transition at some finite g: the situation cannot continuously vary from one unique, symmetric groundstate
S|gsg1 i = |gsg1 i to two symmetry-breaking groundstates:
S|gs± i = |gs∓ i. We’ll make this statement more precise when we discuss the notion of longrange order. First, let’s see what happens when we try to vary the coupling away from theextreme points.
g  1 An excited state of the paramagnet, deep in thephase,
P is achieved by flipping one spin. With H = H∞ =−gJ j Xj , this costs energy 2gJ above the groundstate.
There are N such states, labelled by which spin we flipped:|ni ≡ | → · · · → |{z}← → · · · i,

(H∞ − E0 ) |ni = 2gJ|ni, ∀nnth site

When g is not infinite, we can learn a lot from (1st order) degenerate perturbation theory
in the ferrmomagnetic term. The key information is the matrix elements of the perturbing
hamiltonian between the degenerate manifold of states. Using the fact that Zj | →i = | ←i,so,
Zj Zj+1 | →j ←j+1 i = | ←j →j+1 iX
hn ± 1|Zj Zj+1 |ni = 1,j

the ferromagnetic term hops the spin flip by one site. Within the degenerate subspace, it
acts asHeff |ni = −J (|n + 1i + |n − 1i) + (E0 + 2gJ) |ni.
It is a kinetic, or ‘hopping’ term for the spin flip.
34

Let’s see what this does to the spectrum. Assume periodic boundary conditions and N
sites total. Again this is a translation invariant problem (in fact the same one, basically),which we solve by Fourer transform:(
xj ≡ ja,
1 X −ikxje|ki,|ni ≡ √
k = 2πm, m = 1..NN j
Na
On the momentum states, we have(H − E0 ) |ki = (−2J cos ka + 2gJ) |ki.
The dispersion of these spinon particles is
k→0

(k) = 2J(g − cos ka) ∼ ∆ + J(ka)2(23)

with ∆ = 2J(g − 1) – there is an energy gap (noticethat ∆ does not depend on system size). So these are
k2
massive particles, with dispersion  = ∆ + 2M
+ ... where ∆ is the energy to create one atrest (notice that the rest energy is not related to its inertial mass M −1 = 2Ja2 ).
A particle at j is created by the creation operator Zj :|ni = Zn |gs∞ i.
And it is annihilated by the annihilation operator Zj – you can’t have two spin flips at thesame location! These particles are their own antiparticles.
The number of such particles is counted by the operator
is only conserved modulo two, however.

Pj (−Xj ). The number of particles

What happens as g gets smaller? The gap to creating a spin
flip at large g looks like 2J(g − 1). If we take this formula seriously, we predict that at g = 1 it costs zero energy to create
spin flips: they should condense in the vacuum. Condensing
spin flips means that the spins point in all directions, and thestate is paramagnetic. (We shouldn’t take it seriously because it’s just first order in perturbation theory, but it turns out to be exactly right.)
It’s possible to develop some more evidence for this picture and understanding of thephysics of the paramagnetic phase in the Ising chain by doing more perturbation theory, andincluding states with two spin flips. Notice that for a state with two spin-flip particles, thetotal momentum k no longer uniquely determines the energy, since the two spin-flips canhave a relative momentum; this means that there is a two-particle continuum of states, once
we have enough energy to make two spin flips. For more on this, see e.g. Sachdev (2d ed)§5.2.2. In particular the two spin-flips can form boundstates, which means the two-particle
continuum is actually slightly below 2∆.
[End of Lecture 5]35g  1 Now let’s consider excitations of the ferromagnet, about the state |+i = | ↑↑ · · · ↑i.
P
It is an eigenstate of H0 = −J j Zj Zj+1 and its (groundstate) energy is E0 = −JN . We
can make an excitation by flipping one spin:|··· ↑ ↑ ↑ · ↓ · ↑ ↑ ↑ ···iThis makes two bonds unhappy, and costs 2J + 2J = 4J. But once we make it there aremany such states: the hamiltonian is the same amount of unhappy if we also flip the nextone.
|··· ↑ ↑ ↑ · ↓ ↓ · ↑ ↑ ···i
The actual elementary excitation is a domain wall (or kink), which only costs 2J. The
domain wall should be regarded as living between the sites. It is not entirely a local object, since with periodic boundary conditions, we must make two, which can then moveindependently. To create two of them far apart, we must change the state of many spins.
At g = 0 the domain walls are localized in the sense that a domain wall at a fixed position isan energy eigenstate (just like the spinons
at g = ∞), with the same energy for any position.
P
But now the paramagnetic term − j gXj is a kinetic term for the domain walls:
Xj+1 | · · · ↑↑↑j · ↓j+1 ↓↓ · · · i = | · · · ↑↑↑j ↑j+1 · ↓j+2 ↓ · · · i .
|
{z
} |
{z
}
j̄

=|j̄+1i

Just like in our g  1 discussion, acting on a state with an even number of well-separated
domain walls(Heff − E0 ) |j̄i = −gJ (|j̄ + 1i + |j̄ − 1i) + 2J|j̄iwhere the diagonal term is the energy cost of one domain wall at rest. Again this is diagonalized in k-space with energy
one dwall (k) = 2J(1 − g cos ka)
Again, this calculation is almost ridiculously successful at predicting the location of the phase transition:g→1∆DW = 2J(1 − g) → 0.
Notice that although our discussion of the paramagnetic state g  1 can be applied in any
d ≥ 1, the physics of domain walls is very dimension-dependent.

36Interpretation of the stability of the SSB state in terms of domain walls:
If at finite N , with periodic boundary conditions, we prepare the system in the state |+i,tunneling to |−i requires creation of a pair of domain walls ∆E = 4J, which then move all
the way around the circle, giving the tunneling rateN Yhj̄ + 1|Heff |j̄i
j=1

∆E

∼

1
(gJ)N
∼ g N ∼ e−N log g .
NJ(For g < 1, log g1 > 0.) The tunneling rate goes like e−N – it is exponentially small in thesystem size.

37

2.2.1

Duality between spin flips and domain walls

The discussion we’ve just made of the small-g physics has a lot in common with the large-g
physics. More quantitatively, the dispersion relation one dwall (k) for a single domain wall
looks nearly the same as that of one spin flip (23). In fact they are mapped to each otherby the replacement
1(24)
g → , J → Jg.
g
Notice that this takes small g (weak coupling of domain walls, strong coupling of spin flips)
to large g (strong coupling of domain walls, weak coupling of spin flips).
In fact, there is a change of variables that (nearly) interchanges the two sides of the phase
diagram. Suppose the system is on an interval – open boundaries – the chain just stops at
j = 1 and j = N . (We do this to avoid the constraint of an even number of domain walls.)
We can specify a basis state in the Z-basis by the direction (up or down along Z) of the firstspin and the locations of domain walls.
Consider the operator, diagonal in this basis, which measures whether there is a domain
wall between j and j + 1:(
+1, if zj̄− 1 = zj̄+ 1
2
2= (−1)disagreement .
τj̄x ≡ Zj̄− 1 Zj̄+ 1 =2
2−1, if zj̄− 1 = −zj̄+ 1
2

2

Notice that τj̄2 = 1, τj̄† = τj̄ . Similarly, consider the operator that creates a domain wall at
j̄:
Y
Xj .
τj̄z ≡ Xj̄+ 1 Xj̄+ 3 · · · =
2

2

j>j̄

This operator flips all the spins to the right of the link in question (and fixes our reference
first spin). It, too, is hermitian and squares to one. Finally, notice thatτj̄z τj̄x0 = (−1)δj̄j̄0 τj̄x0 τj̄z
just like Z and X (since when j̄ = j̄ 0 , they contain a single Z and X at the same site). The
domain walls can be represented in terms of two-state systems living on the links.16
Notice that the inverse of the map from X, Z to τ x , τ z iszz
Xj = τj−
.
1τ
j+ 1
2

2

(The right hand side is an inefficient way to flip a single spin at j: namely, flip all the spinsright of j − 1, and then flip back all the spins to the right of j.)
16

Note that in lecture I reversed the names of τ x and τ z ; I think this way is a little better.

38

So the 1d TFIM hamiltonian in bulk is
HTFIM = −J

X

(gXj + Zj Zj+1 )

j

= −J

X


z
gτj̄z τj̄+1+ τj̄x .

(25)j̄This is the TFIM hamiltonian again with Z → τ z and X → τ x and the couplings mappedby (24).
This is in fact the same map as Kramers-Wannier duality (or rather it is mapped to it by
the quantum-to-classical map). As K&W argued, if there is a single phase transition it mustoccur at the self-dual point g = 1.
Notice that the paramagnetic (disordered) groundstate of the original system is a condensate of domain walls, in the following sense. The operator that creates a domain wall has an
expectation value:Yhτj̄z i = hgsg=∞ |τj̄z |gsg=∞ i = hgsg=∞ |Xj |gsg=∞ i = 1 ∀j̄ .
j>j̄(For g ∈ (1, ∞), this expectation value is less than one but nonzero, just like how |hZi|
decreases from 1 as g grows from zero.) Although there is a condensate, there is no order, inthe sense that an expectation value of X does not break any symmetry of HTFIM . (There isanother state where hτ x i = −1, namely the one where all the spins are pointing to the left.
But (at large g) it’s a very high-energy state.)
An important point (and the reason ‘duality’ is a dangerous word): the two sides of thephase diagram are not the same. On one side there are two groundstates related by thebroken symmetry, on the other side there is a unique symmetric groundstate. That’s how
we knew there had to be a phase transition! I will say more about this mismatch.
Open boundaries. Let us make sure we can reproduce the correct number of groundstatesin the two phases. To get this right, we have to be precise about the endpoint conditions.
Let’s study the case where we have N sites in a row; the first and last sites have only one
neighbor. The Hamiltonian is!
N−1
X(gXj + Zj Zj+1 ) + gXN .
HTFIM = −J
j=1

The duality map isx
Zj Zj+1 = τj+
1 , j = 1, 2...N − 1,2zz
Xj = τj−
,
1τ
j+ 1
2

j = 1...N2

In terms of the domain-wall variables, the hamiltonian is
!

N
−1 
X
1
x
z
z
HTFIM = −Jg
τj−
+ τj+
+ τNz − 1 τNz + 1 .
1
1τ
j+ 21
2
2
2
2
g
j=1
39

.

But now there are two special cases:• τNz + 1 = 1: this operator flips all the spins with j > N ; but there are no spins with j > N .
2
So it is the identity operator.
• τ 1x : this operator measures whether or not there is a domain wall between j = 1 and
2
j = 0, τ 1x = Z0 Z1 . But there is no spin at j = 0. One way to think about this is to put
2
a “ghost spin” at j = 0 which is always in the state Z0 = 1. So τ 1x = Z1 : it measures the2value of our reference spin.
Pxxx
At g = 0: Hg=0 = −J N
j=2 τj− 1 and the groundstate is τj− 1 = 1 for j = 2...N . But τ 1222
does not appear, so there are two degenerate groundstates, eigenstates of τ 1x with eigenvalue
2±, which are just |±i, the states with no domain walls: all the other spins agree with thexfirst one in a state where τj̄>1
= 1.

P
N −1 zzz. The first term requires agreementτ+ττAt g = ∞, Hg=∞ = −Jg
1
1
1
j=1
N− 2
j− 2 j+ 2zz
between τj− 1 = τj+ 1 for j = 1...N − 1, but the second term requires τNz − 1 = 1. So all the222others must be +1, too. This is the unique groundstate of the paramagnet.

Comment: In the notation of our lattice model (14), the critical curve between the twophases is atsinh(2Kx ) sinh(2K) = 1 .
This relationship follows by Kramers-Wannier duality. In the scaling limit, Kx ∼ ∆τ 
1, K  1. So the critical condition is
1 = sinh(2Kx ) sinh(2K) ' 2Kx

1
1
e2K
= J∆τ
= .
2
Jg∆τ
g

This determines the location of the quantum phase transition, in agreement with all theother stuff we’ve said about it.

40

2.2.2

Mean field theory

So far we have used perturbation theory about g = 0 and g = ∞ to delicately inch our way
toward the interior of the phase diagram of the TFIM in one dimension (22). Here we takea plunge and try to guess the groundstate for all g. The nice thing about trying to guess
the groundstate is the Rayleigh-Schrödinger hedge: the energy expectation in any state is
an upper bound for the groundstate energy; minimizing the energy within a class of guesses
is called a ‘variational approach’.
The name for the particular guess we’ll make is ‘mean field theory’, which means that
we completely ignore entanglement between different sites, and suppose that the state is a
product state|M F T i = |ψ1 i ⊗ |ψ2 i · · · |ψj i · · · ...
If we further assume translational invariance then the state at every site is the same and wehave one bloch sphere to minimize over for each g:θ −iϕ/2
θ iϕ/2| ←i .
|ňi = ⊗j | ↑ň ij = ⊗j cos e | →i + sin e22j(Here θ is the angle ň makes with the x axis, and ϕ is the azimuthal angle in the yz plane,from the z-axis.) To evaluate the energy expectation in this state, we only need to know
single-qbit expectations:h ↑ň |X| ↑ň i = cos θ,h ↑ň |Z| ↑ň i = sin θ cos ϕ.

So the energy expectation isE(θ, ϕ) = −N J sin2 θ cos2 ϕ + g cos θ .
This is extremized when ϕ = 0, π and when0 = ∂θ E = N J sin θ (2 cos θ − g) .
Notice that when θ = 0, the two solutions of ϕ are
the same, since the ϕ coordinate degenerates at thepole. The solutions at cos θ = g/2 only exist wheng/2 < 1. In that case they are minima (see the figure)since ∂θ2 E|cos θ=g/2 > 0, while ∂θ2 E|θ=0 = N J(g − 2)is negative for g < 2. (Notice that ϕ = π can be
included by allowing θ ∈ (−π, π], as in the figure.)
So mean field theory predicts a phase transition at
g = 2, from a state where hZj i = sin θ to one wherehZi = 0. It overestimates the range of the orderedphase because it leaves out fluctuations which tend
to destroy the order.
41

Let’s study the behavior near the transition, where θ is small. Then the energy can beapproximated by its Taylor expansion
g−2 2 1 4
θ + θ
E(θ) ' N J −2 +24(where I have set g = gc = 2 except in the crucial quadratic term). This has minima at√
hZj i = sin θ ' θ = ± gc − g .
(26)
The energy behaves like

EM F T (g) =32 4 (gc − g) ,g < gcg ≥ gc0,

Notice that ∂g E is continuous at the transition. (Recall that the groundstate energy ofthe quantum system is equal to the free energy of the corresponding stat mech system,so ∂g E ∝ ∂T F continuous is the same criterion for a continuous transition.) So mean
field theory (correctly) predicts a continuous quantum phase transition between the orderedphase and the disordered phase. The location of the transition is wrong (mean field theory
overestimates the size of the ordered region because it leaves out lots of order-destroying
fluctuations), and so are other properties, such as the exponent in (26), which should be 1/8
instead of 1/2.
[End of Lecture 6]42

2.2.3

Correlation functions and long-range order

A useful set of physical (observable! as we’ll discuss in a moment) quantities are the correlation functions of the spins:
C(r, r0 ) ≡ h0|Zr Zr0 |0i .
Notice that this is an equal-time correlator. If the groundstate is translation invariant, thisis only a function of the separation between the spins C(r, r0 ) = C(r − r0 ). (Notice that
I am using rs and js to label the positions interchangably, but r = ja, where a is thelattice spacing) We can take further advantage of translation symmetry by thinking about
the Fourier transform of this quantity, which is called the static structure factor:
X
S(q) ≡
e−iqr C(r) ;rthe sum over positions is over r ∈ a{1...N }, and (assuming periodic boundary conditions
again) the argument q ∈ N2πa {1...N }.
Correlations in the quantum Ising chain
At g = ∞, all the way in the paramagnetic phase,h0|Zj Zl |0i|g=∞ = δjl– different sites are totally uncorrelated. At finite but large g  1 we have instead:h0|Zj Zl |0i|g1|xj −xl |a'e−|xj −xl |/ξ(27)– the fluctuations of the spin flips communicate between the sites, but only over a shortrange, because the spin flips are massive.
In contrast, for g  1, in the ferromagnetic phase,h0|Zj Zl |0i|g1|xj −xl |a'N02 (g)(28)

where N02 (g) = 1 for g = 0. N0 < 1 for g > 1.
More simply, if we train the system appropriately (as we discussed in point 3), then we’ll
haveh0|Zj |0i = ±N0 (g).
N0 is the order parameter, the magnetization. Another way to state the issue about theorder of limits of h → 0, N → ∞: in a finite system the magnetization is zero:
lim lim hZi i = 0,

N →∞ h→0but

lim lim hZi i = N0 .

h→0 N →∞In this sense, the ferromagnet is said to have long range order.
43

OK but here’s the simple important point: There’s no smooth function of g that caninterpolate between (27) and (28). A virtue of the notion of long-range order is it gives
sharp distinctions between phases.

2.2.4

Physics of structure factors

Consider the limit where N → ∞ so that the allowed momenta fill in a continuous Brillouinzone q ∈ 2π
{1..N } → (0, π] ' (−π, π]. So the inverse Fourier transform isN
Z πdq
S(q)eiqR .
C(R) =
2π−π
When there is long range order,R→∞

C(R) ∼ N02
what does S(q) do? It means a singularity at q = 0:
S(q) = 2πδ(q)N02 + regular(29)

where ‘regular’ means terms which are smooth at q = 0. This singularity is a sharp signatureof the broken symmetry.
We can learn something by considering a spectral representation of S(q), obtained byresolving the identity by energy eigenstates |ni:X
|hn|Zq |0i|2 .
(30)
S(q) = h0|Zq Z−q |0i =nHere I used Zq ≡ √1N

P

re

−iqr

Zr = Z†−q .

A sum rule: C(R = 0) = h0|Zr Zr |0i = h0|0i = 1. But the (inverse) Fourier transform is
Z πdq
1 = C(R = 0) =S(q) .
−π 2πThis is some short-distance information which remembers that the objects we’re looking
at square to one. Given the spectral interpretation, this sum rule means that as we varyparameters, the spectral weight can move around, but it can’t go away. In particular, when
the system orders as in (29), the 2πN02 contribution at q = 0 has to come from somewhere.
Dynamical structure factor

C(r, t) ≡ h0|Zr (t)Z0 (0)|0i.
44This reduces to C(r) if we set t = 0. It is useful to think of Zr (t) in Heisenberg representation:
Zr (t) = e+iHt Zr e−iHt .

S(q, ω) =

X

e

−iqr

Z ∞
−∞

r

Notice that S(q) =structure factor is:R

dte+iωt C(r, t).

dωS(q, ω) (not S(ω = 0)). The spectral representation of the dynamical
X
S(q, ω) =|hn|Zq |0i|2 2πδ(ω + E0 − En ) .
(31)nSo (30) is obtained by integrating over ω 17 which simply removes the ω dependence.
In words: S(q, ω) is the number of states of momentum q and energy ω, weighted by (thesquare of) their overlap with the state engendered upon the groundstate by the spin operator.
This formula is crying out to be interpreted in terms of Fermi’s Golden Rule.
Before we do that let’s include finite temperature, too. Let
C(x, t; x0 t0 ) =1 −H/T
tre
Z(x, t)Z(x0 t0 )Z

where again Z are Heisenberg picture operators. Its fourier transform is
Z
Z ∞
0
0
d
dt eiω(t−t )−i~q·(~x−~x ) C(x, t; x0 t0 ).
S(q, ω) = d x−∞

Notice that the time integral here is over real time, not euclidean
time. To find the specPtral representation, again, insert an energy eigenbasis 1 = n |nihn| in between each time
evolution operator and each field, and do the space and time integrals. The result is:S(q, ω) =2π X −En0 /T2|hn0 |Z(q)|ni| δ (ω − En − En0 ) .
e
ZV n,n0

V is the volume of space.
Now let’s interpret this in terms of a transition rate in a scattering experiment, via Fermi’s
Golden Rule. The initial state is |n0 i with probability pn0 = e−En0 /T /Z given by the Boltzmann distribution. (In the T → 0 limit, we return to (31), where the initial state is thegroundstate with probability 1.) The perturbing hamiltonian is linear in the spin operator
Hpert ∝ Z(q). With these assumptions, the transition rate per unit time to the state |mi isexactly S(q, ω). This explains the connection between the dynamical structure factor andexperiments involving scattering of neutrons: a neutron couples to the spins via its magnetic
dipole moment, which couples linearly to the spin operator; ω, q are the change in energy andR1with the usual 2π
convention: d̄ωS(q, ω) = S(q). Thanks to Sid Parameswaran for alerting me to my
2π problem here. His notes for a related course can be found here. You should trust his signs over mine!
1745

momentum of the system. The transition probability (at zero temperature) is proportional
to something like:~ |0i ⊗ |neutron, initiali|2 .
|hn| ⊗ hneutron, final|~Sneutron · σ
The (known) neutron states allow us to determine q, ω. You could ask me how we do an~.
experiment which couples just to Z, and not all of σ
Given our understanding of the spectrum at large and small g, we have a prediction forthe behavior of S(q, ω): at large g, the lowest-energy excitation is a single spin flip, whichconveniently is created by Z, and so contributes to S(q, ω). But it only contributes if q andω are related in the right way, namely ω = q . This means that there will be a sharp peakalong this curve: for ω in the right range, S(q, ω) = Zδ(ω − q ). For larger ω, we access themulti-particle continuum. In fact, Z|0i only includes states with an odd number of particles,
so the continuum in S(q, ω) starts near 3∆. What happens at small g? At small g, the spin
flip is not stable – it can decay into a pair of domain walls, which are free to move apart.
Although these domain walls behave in many ways like particles, the operator Zj cannot
create a single domain wall. (Recall that the domain wall creation operator is τj̄x which ismade of a string of Xs.) So the spectral function of Z has only a continuous mess at smallg – a multiparticle continuum of domain walls, no sharp delta-function peak.
People18 have done this experiment on a system which is well-described by the TFIM
(CoNb2 O6 ). It can be approximated as a collection of decoupled chains of spins with nearestneighbor ferromagnetic Ising interactions. The transverse field term can be produced byapplying a (surprise) transverse magnetic field, and hence g can be varied. Indeed at large
g one sees a sharp single-particle peak, and at small g there is a featureless mess. Here is asummary of the results of the experiment (from Ashvin Vishwanath):

A further feature
P of these experiments which I must mention is: if we apply a longitudinal
field, ∆H =j hz Zj , we confine the domain walls, in the following sense. This term
18

Coldea et al, Science 327 (2010) 177.

46

means that a region of down-spins costs an energy proportional to its length. This is aninteraction energy between the domain walls at the edges of the region which is linear in
their separation, a constant attractive force between them. Their spectrum of boundstates
is extremely interesting and represents a physical appearance of E8 symmetry.

47

2.2.5

Solution of Ising chain in terms of Majorana fermions

Where does the Ising model come from? Just as the classical ising model arises inmany ways (lattice gas, double well potential....), quantum spin systems arise in many ways.
for example:
1. spinful fermions at half-filling, with charge dofs quenched.
2. Similarly, spinless fermions hopping around on a lattice also realize a spin system:
there are two states per site. How to relate fermion operators and spin operators? In
1+1 dimensions, we can answer this in complete explicitness, below.
3. If we view down-spin as vacuum and up-spin as the presence of a boson, we canview a hard-core boson system as a spin-1/2 system. here it is natural to preserve aU(1) ⊂ SU(2) symmetry which rotates the phase of X + iY. the conserved quantity is
Z = (2nB − 1) where nB is the boson number.
[End of Lecture 7]

Jordan-Wigner in 0+1 dimensions
Consider a single fermionic operator c satisfying
c2 = 0,{c, c† } = 1 .

This constructs the Hilbert space of a qbit as follows:c| →i = 0,c† | →i = | ←i, c† | ←i = 0, c| ←i = | →i .

The two states of the qbit just describe the presence or absence of the fermion in this state.
We can rewrite the Pauli operators as
Z = c + c† ,which flips the spin,Y=1c − c† ,i

X = 1 − 2c† c.
Here c† c counts
the number of fermions. Also note that the raising and lowering operators01
are σ + ≡= c† , σ − = c and the number operator c† c (whose eigenvalues are 0 and
0 01) is c† c = 12 X + 1.

48

The story is less simple with more than one fermion mode operator, because they mustsatisfy {cα , cβ } = 0, rather than just commuting, as the would if we just used the above
map for each mode.
[0904.2771, Fradkin] Let’s look at the TFIM again:X
HTFIM = −J
(gXj + Zj Zj+1 )jhas a phase transition between large-g and small-g phases.
(Disordered) large g: excitations are created by Zj – they are spin flips. The groundstateis a condensate of domain walls: hτ z i =6 0.
Q(Ordered) small g: excitations are created by the ‘disorder’ operator τj̄z = j>j̄ Xj – they
are domain walls. The groundstate is a condensate of spins hZj i =6 0, i.e. a ferromagnet.
So we understand what are the ‘correct variables’ (in the sense of the introductory discussion
to this course) at large and small g. I claim that the Correct Variables everywhere in thephase diagram are obtained by ”attaching a spin to a domain wall”. These words mean thefollowing: letYzχj ≡ Zj τj+
Xj 0
1 = Zj2χ̃j ≡j 0 >jYz
Yj τj+ 1 = −iZj
Xj 02
j 0 ≥j(32)

The first great virtue of this definition is that these operators agree with the creators of theelementary excitations in both regimes we’ve studied: When g  1, hZj i ' 1 and morezzstrongly, Zj = hZj i + small, so χj ' hZj iτj+
, the domain wall creation operator.
1 ' τ
j+ 12
2
z
Similarly, when g  1, τjz ' 1+ small, so χj ' Zj hτj+
1 i ' Zj , which is the spin flipper on2
the paramagnetic vacuum.
Now let us consider the algebra of these χs:• They are real: χ†j = χj , χ̃†j = χ̃j .
and• They are fermions:if i 6= j, χj χi + χi χj ≡ {χj , χi } = 0,{χ̃j , χ̃i } = 0, {χj , χ̃i } = 0.

(33)This is because the spin flip Zj in χj changes sign when it moves through the domain wall
created by χi . When they are at the same site:χ2j = 1 = χ̃2j .

In summary:{χi , χj } = 2δij , {χ̃i , χ̃j } = 2δij , .
49

Notice that (33) means that χi cares about χj even if |i − j|  1. Fermions are weird andnon-local!
Notice that these fermions don’t have a spin index. I should mention that this change ofvariables is called Jordan-Wigner.
Real fermion operators like this are called Majorana fermion operators. We can make morefamiliar-looking objects by making complex combinations:1(χj − iχ̃j )
2

cj ≡

=⇒

c†j =

1
(χj + iχ̃j )2

These satisfy the more familiar anticommutation relations:{ci , c†j } = δij , {ci , cj } = 0,{c†i , c†j } = 0, 2
and in particular, c†i = 0, like a good fermion creation operator should.
We can write HTFIM in terms of the fermion operators. We need to know how to write Xj
and Zj Zj+1 .
The operator which counts spin flips in the paramagnetic phase is†

Xj = −iχ̃j χj = −2c†j cj + 1 = (−1)cj cj .
(To get this we can use (32) and YZ = iX and (τ z )2 = 1.) Here c†j cj = nj measures thenumber of fermions at the site j and is either 0 or 1, since they are fermions. At each site
| →j i = |nj = 0i,| ←j i = |nj = 1i,.

The number of spin flips is the number of fermions.
The operator which counts domain walls is
Zj Zj+1 = iχ̃j+1 χj .
(Check: iχ̃j+1 χj = iYj+1

Qk≥j+2 Xk Zj

Ql≥j+1 Xl = iYj+1 Xj+1  Zj .)





|

HTFIM = −J



X

{z

=Zj+1

}

(iχ̃j+1 χj + giχj χ̃j )jis quadratic in these variables, for any g! Free at last! (It is quadratic in the cs, too, since
they are linear in the χs.)50

Comments:• Notice that the relation
Xj = 1 − 2c†j cj
is exactly implementing the simple idea that spinless fermions on a lattice produce
two-states per site which we can regard as spin up or spin down (in this case it’s up or
down along x): The states X = ±1 correspond to c† c = 0 and 1 respectively.
• Notice that the description in terms of majoranas is preferred over the complex fermionsbecause the phase rotation symmetry generated by the fermion number c† c is not asymmetry of HTFIM – in terms of the cs, it contains terms of the form cj cj+1 whichchange the total number of c fermions (by ±2). It is the hamiltonian for a superconductor, in which the continuous fermion number symmetry is broken down to a Z2
subgroup. Fermion number is still conserved mod two, and this is the Z2 symmetry of
the Ising model, which acts by Z → −Z.
• A useful thing to remember about majorana operators {χ, γ} = 0 is that (iχγ)† =−iγχ = +iχγ is hermitian.
• Another useful fact:†c(−1)c c = −c(34)which is true because the BHS only nonzero if the number is nonzero before the annihilation operator acts, in which case we get (−1)1 . Similarly (the conjugate equation),††(−1)c c c† = −c† , and (−1)c c c = c and so on.
• This procedure of “attaching spin to a domain wall” led to fermions. This maybe
isn’t so surprising in one dimension. But there are analogs of this procedure in higher
dimensions. In 2+1 dimensions, the analog is to attach charge to a vortex (or to
attach magnetic flux to charge). This leads to transmutation of statistics from bosons
to fermions and more generally to anyons and the fractional quantized Hall effect. In
3+1 dimensions, the analog seems to be attaching charge to a magnetic monopole toproduce a ‘dyon’.

51

Dual fermionsLetγj+ 1 = −χ̃j+1 , γ̃j+ 1 = χj .
2

2

We are dividing up the pair of majorana modes that we had previously associated with thesite j between the two neighboring links.
In terms of these H is
X

HTFIM = −Jiγ̃j̄ γj̄ + giγ̃j̄+1 γj̄.

j̄=j+ 21We have interchanged the form of the two terms. The two expressions are related by χ ↔
γ, J ↔ Jg, g ↔ 1/g.
Now let’s think about the two phases in terms of the fermion fields.
g  1 The disordered state, with a unique groundstate is governed by
X
X†
Hg→∞ = −Jgiχj χ̃j = −Jg(−1)cj cj .
jj

The groundstate of this is just the state |0i with no fermions, cj |0i = 0 for all j, on which†(−1)cj cj = 1.
g1
Hg→0 = −JXiχ̃j+1 χj .

jThis is a hopping term for the χ fermions, which will be diagonalized by some bandstructure,
which will then be partially filled by fermions. It’s easier just to use the dual fermion
variables, in terms of whichXXč† č
Hg→0 = −Jiγ̃j̄ γj̄ = −J(−1) j̄ j̄
j̄=j+ 21j̄=j+ 21
where čj̄ ≡ 21 γj̄ − iγ̃j̄ are the dual complex fermions. The groundstate of this is just thestate with no dual fermions: čj̄ |0̌i = 0, ∀j̄.
Again we should think about the endpoint conditions. Consider what happens for an open
chain. It’s useful to account for the states as follows: think of dividing up each site into a pairof sites (white and black, as in the figure, one from the lattice, one from the dual lattice) with
one majorana mode living at each. To help visualize, let’s say χ̃j = −γj− 1 lives at the white2
site (at j− 21 ) and χj = γ̃j+ 1 lives at the red site to its right (at j). Note that γ 1 = −χ̃1 , γ̃ 1 =2
2
2
χ0 . The g (para) term is a coupling between the white and red sites at fixed j. The ferro
term couples a red site at j − 21 to the next white site at j. This means that at small g, in an
open chain, one site at each end will be left out – will not appear in the hamiltonian at all.
52

This results in unpaired majorana
zeromodes at the ends of the chain.
Let a† ≡ 21 (iχ̃1 + χN ). The algebra {a, a† } = 1 must be represented on the groundstates. Thisrequires a pair of states(g  1)(g  1)|0i = |+i, a† |0i = |−i.
The SSB degeneracy of the ferromagnet is mapped by this fermionization map to a topologicaldegeneracy in terms of the fermions.
Splitting of the energy these two states is small because the modes are separated in space
and the bulk is gapped. The modes are separated by a distance much bigger than thecorrelation length, ξ ∼ ∆1 . The splitting comes from a term
∆H = a† a = iχ̃1 χN
and we (again) estimate that  ∼ e−N/ξ .
[End of Lecture 8]53

2.2.6

Continuum limit[Sachdev, p. 139 of 2nd ed, p. 48 of 1st ed]
We found in 2.1 that a quantum lattice model (with continuous time and discrete space)can be Trotterized into a lattice statistical mechanics model (with discrete (euclidean) timeand discrete space). We also found a reverse map, via the continuum scaling limit. In thissubsection we discuss the analogous scaling limit that takes a quantum lattice model to acontinuum quantum field theory (with continuous time and continuous space).
Label site j of the lattice by its position xj = ja, where a is the lattice spacing. L = N ais the length of the chain. Wavevectors lie in the interval k ∈ (− πa , πa ].
In terms of complex fermions
Xj = 1 − 2c†j cj ,Zj = −





Y
Y
†
(1 − 2c†i ci ) cj + c†j = − (−1)ci ci cj + c†j .
i>jand their Fourier modes ck ≡ √1N
HTFIM = J

X

i>j

P

−ikxj

, the TFIM hamiltonian is
2(g − cos ka)c†k ck − i sin ka c†−k c†k + c−k ck − g
j cj ek

Notice that the second set of terms violates fermion number conservation in units of two; thisis the same statement that the ferromagnetic term creates spin flips in pairs. The constantat the end is irrelevant unless we are after the Onsager free energy.
This Hamiltonian is quadratic in ck s, but not quite diagonal. The final solution for the
spectrum involves one more operation the fancy name for which is ‘Bogoliubov transformation’, which is the introduction of new (complex) mode operators which mix particles andholes:
γk = uk ck − ivk c†−kDemanding that the new variables satisfy canonical commutators {γk , γk†0 } = δk,k0 requires
uk = cos (φk /2) , vk = sin (φk /2). We fix the angles φk by demanding that the hamiltonianwithin terms of γk be diagonal – no γk γ−k terms. The resulting condition is tan φk = 21 (k)
(k)


p
P
1 (k) = 2J(g − cos ka), 2 (k) = −J sin ka, and H = k k γk† γk − 12 , with k = 21 + 22 .

The end result is that the exact single-particle (single γ) dispersion ispk = 2J 1 + g 2 − 2g cos ka .
The argument of the sqrt is positive for g ≥ 0. This is minimized at k = 0, which tells us the exact gap at all g:k ≥ 0 = 2J|1 − g| = ∆(g)54

which, ridiculously, is just what we got from 1st order perturbation theory on each side of
the transition.19
Comment on relation to Onsager. Notice that by the quantum-to-classical mapping,
this solution has all the universal information of Onsager’s solution. In particular, his exact
free energy can be obtained by computing the groundstate energy of the fermions. (See thehomework.)
The critical theory is scale invariant. At g → 1, the fermions at k = 0 become gapless:kg=gc =1=c|k|,where the speed of propagation is c = 2Ja. There is an emergent Lorentz symmetry, at leastin the spectrum, at the critical point. And the fact that the long-wavelength (k ∼ 0) havethe lowest energy allows a continuum description that forgets the lattice details.
Near g → gc = 1 (i.e. |g − gc |  gc ),v
2ug − gcg∼gc =1 u 2
k = cuk +ta| {z }
m2 →0a
There is a diverging length scale, which is ξ = m1 = |g−g.
c|Some useful language:• The correlation length critical exponent is ν in ξ ∼ |g − gc |−ν , the rate at which thecorrelation length diverges near the critical point. (Here ν = 1.)• The dynamical critical exponent is the power law in the dispersion at the critical point:
(k) ∝ k z .

(35)

Notice that if we rescale space and time likex → λx,t → λz t(36)

it will preserve the dispersion (35). (Here z = 1, as in a relativistic theory where timeand space can be boosted into each other.)We can expect the continuum description (at least the long-wavelength information) to beinvariant under the scale transformation (36).
19

Notice by the way that this ‘single-particle excitation’ in terms of the spins is a state with a single
domain wall, and therefore only makes sense in the thermodynamic limit; in finite volume, the spin hilbert
space only contains states with an even number of fermions. See §2.2.7 below for more on this point.

55

When ξ  a, define the continuum fermion field1
Ψ(xj ) ≡ √ cja√where the factor of a is designed to convert the kronecker deltas in{cj , c†j 0 } = δjj 0

into {Ψ(x), Ψ† (x0 )} = δ(x − x0 ).

The fourier transform is

e−ikx
dx √ Ψ(x) .
LLet’s plug this into HTFIM , focussing on (the lightest) modes with ka  1.
Z

ck =

Z

X
J
(g − cos ka)c†k ck(g − gc )dx Ψ† (x)Ψ(x)kHere themeans we are dropping terms of higher order in the expansion in powers of a∂x .
Z
X
c
−iJ
sin kac−k ckdx Ψ(x)† ∂x Ψ(x)†
2
k

So
H

c
2

Z

†

†


dx Ψ(x) ∂x Ψ(x) − Ψ(x)∂x Ψ(x) + ∆ZdxΨ† Ψ

with ∆ = 2J|g − 1|.

[Fradkin, 2d ed p. 119] To understand the significance of this result it will be best to use
the (fermion coherent state) path integral, which we’ll do below in §2.2.8. But we can learnsomething just by thinking about the Heisenberg equations of motion.

HTFIM = −J

X

(gXj + Zj Zj+1 )

j

= −iJ

X

(gχ(j)χ̃(j) − χ(j)χ̃(j + 1)) .

j

Notice that the Z symmetry action on the fermions isYYS=
Xj = iN
χ(j)χ̃(j).
jj

The Heisenberg eom are i∂t O = [H, O]:i∂t χ(j) = iJ (g χ̃(j) − χ̃(j + 1))
56

(37)

i∂t χ̃(j) = iJ (−gχ(j) + χ(j − 1))(38)(Free fields means linear equations of motion.) In the continuum,χ(j + 1) ' χ(xj ) + a∂x χ(xj ) + O(a2 ∂x2 )
so


1
1−g
∂t χ(x) ' −χ̃(x) − ∂x χ̃(x)
aJ
aJ

1−g
1
∂t χ̃(x) 'χ(x) − ∂x χ(x)
aJ
aJ

(39)

If we let χ± = 12 (χ ∓ χ̃) and rescale the speed of propagation into the time variable, t ≡aJx0 ,∂0 χ+ = ∂x χ+ + mχ−∂0 χ− = −∂x χ− − mχ+(40)with m = 1−g. When m → 0, at the critical point, these are chiral fermions:a(∂0 ∓ ∂x ) χ± = 0.
The χs are still real, so these are majorana fermions. In a language we will introduce soon,
this is an example of a 1+1d CFT with central charge c = ( 21 , 12 ).
Even a bit away from the critical point, we can reconstruct the relativistic symmetry, andorganize (40) into following Dirac equation
0 = iγ µ ∂µ χ + imχ,
Notice that a spin index has emerged to save the spin-statistics theorem. Here we havechosen a basis of gamma matrices where




0 1
0 1
0
1
γ =
, γ =
−1 0
1 0
(in the ± basis) which satisfy {γ µ , γ ν } = 2η µν as they should. So the Dirac equation above
is

 
m
∂0 + ∂xχ+0=i.
−∂0 + ∂xm
χ−(The Majorana condition I am imposing is χ̄ = χT γ 0 , aka χ† = χT , so χ± are real.)572.2.7

Duality and Jordan-Wigner with periodic boundary conditions

How do I define the creation operator for a single domain wall with periodic boundary
conditions on the spins? I can’t. The fact that I can only make an even number means thatthe sign of the operator is not well-defined – in fact you can’t measure the single-domain
wall operator. And since the fermions χ are made from this operator, you can’t measure asingle fermion operator either.
Although there is a sign ambiguity, we can still define the domain wall operators:
YXjτj̄x = Zj̄− 1 Zj̄+ 1 , τj̄z =22

N ≥j>j̄– in the second definition, we have made a choice: even though ZN +1 = Z1 , we’ve decided
to pick out the link where j̄ = N + 21 to stop the string of Xs. As we saw in the case of open
boundary conditions, this means that τ 1 and τN + 1 are special:22

τNz + 1 = 1
2

since the product is over the empty set.
z

τ1 =

N
Y

2

Xj = S

j=1

is the symmetry generator. This is significant because this operator commutes with HT F IM ,it means that eigenstates of τ 1z are energy eigenstates. And notice that it determines the2boundary conditions on the τ s in following sense: combining the previous two equations,

Q
τ 1z
if j Xj = S = 1z
2
1 = τN + 1 =Q2
−τ 1z
if j Xj = S = −12This operator τ 1z which says whether the τ s are periodic or antiperiodic is a dynamical
2
variable.

58

Gauge theory formulation of the 1d duality
There is a way to get around the obstacle of no-lone-domain wall on a circle, which generalizes to higher dimensions (but it is kind of overkill in this case). It is a reformulation of
the TFIM as a lattice gauge theory. This is worth doing because in duality transformations
in higher dimensions, the gauge theory formulation cannot be avoided. So this is a context
where we can introduce the concepts of gauge theory without lots of other complications.
We begin by enlarging our system, by adding, in addition to the real qbits on the sites of
the chain another set of fake qbits on the links:

enlarge H = ⊗j Hj ⊂ Hbig ≡ (⊗j Hj ) ⊗ ⊗j̄ Ȟj̄ .
Anticipating the answer, I will label the two states on a link j̄ as |τj̄z i = | ± 1i. So as for any
two-state system, associated with each link there are pauli operators
τj̄z |τj̄z i = τj̄z |τj̄z i, τj̄x |τj̄z i = | − τj̄z i , τj̄z τj̄x0 = (−1)δj̄j̄0 τj̄x τj̄z0 .
We want to pick out a subspace of this larger fake Hilbert space which is the real one.
One way to think about this is to declare that the following ‘gauge transformation’ is anequivalence:G : τj̄z → sj̄ τj̄z , sj̄ = ±1 .
(41)
Notice that this group G is big – we pick an element of Z2 (i.e. a sign) sj̄ = ±1 for eachlink of the lattice, its order is |Z2 |N = 2N . If a configuration of τ z s is equivalent to its image
under this map, we can just pick sj̄ = τj̄z and this means that any configuration of τ z s isequivalent to 1: we can just get rid of the τ s and we are back at our original Hilbert space.
This is called unitary gauge. Below we will implement this more concretely. But first we
must decide how the gauge redundancy G acts on our original variables. It will be useful tochoose(42)
G : Xj → sj− 1 Xj sj+ 1 , Zj → Zj .
22This is the transformation law of a link variable in lattice gauge theory. The site of theoriginal lattice is a link of the dual lattice. Let me use notation which emphasizes this:xz
Xj ≡ σj−, Zj ≡ σj−
.
1
1
,j+ 1
,j+ 1
2

2

2

2

The generator of the gauge transformation isz
x z
Gj̄ = σj̄−1,j̄ τj̄ σj̄,j̄+1 .

What I mean by ‘generator’ is that the action of the symmetry on states is |ψi → G|ψi andthe action on operators is
O → GOG† .
The operator Gj̄ generates the transformation where sj̄ = −1 and all the other sj̄ 0 = +1.
Check that this reproduces (41) and (42).
59The physical hilbert space is cut out of Hbig by projecting to eigenstates of these operatorswith eigenvalue 1:
Hphys = {|ψi ∈ Hbig such that Gj̄ |ψi = |ψi, ∀j} ⊂ Hbig .
This is just the statement that physical states don’t change under our made-up redundancy20 .
In order for this construction to make sense, our hamiltonian acting on the big hilbert spacemust be gauge invariant. That is, we must have
[H, Gj̄ ] = 0, ∀j̄.
This way, we can simultaneously diagonalize H and Gj̄ . A gauge invariant hamiltonian
acting on the big hilbert space is:
X

z
z
z x
z
Hbig = −J
σj̄−1,(43)j̄ σj̄,j̄+1 + gτj̄ σj̄,j̄+1 τj̄+1
j̄In unitary gauge, τ z = 1, we can erase the τ z s and this reduces to HTFIM . This is why I
chose it.
On the other hand, using the local constraintszzx
1 = Gj̄ ⇔ σj̄−1,j̄ σj̄,j̄+1 = τj̄

on Hphys ,we have
Hbig = −JXzτj̄x + gτj̄z σj̄,x j̄+1 τj̄+1
.

(44)

j̄

In this last expression, we still have the gauge freedom to play with – too many variables.
Let’s use it to get rid of the original σs. On an open chain, we can use the transformation
σj̄,z j̄+1 → sj̄ σj̄,z j̄+1 sj̄+1(45)

to set σ z = 1 for each j̄. In that case, the expression (44) precisely reduces to our dual
description of the TFIM in terms of domain wall operators, (25).
But with periodic boundary conditions, the object
W≡NY

σj̄,z j̄+1

j̄=1

is a sign which does not change under the transformation (45). It is called a Wilson line
operator. We can’t get rid of it. In fact it is the familiar object
YW=
Xj = S,j20

In lecture I called this constraint the ‘Gauss’ Law constraint’. This is a correct and useful terminology
which I have not explained yet. Sorry for the extra jargon. The reason it is correct is that if you really squint
at this expression you can recognize it as a lattice version of the divergence of the electric flux coming outof the site j̄. This will become clear in §5.2.

60

the Z2 symmetry charge.
You can see from (43) that if W = −1, we have APBCs on τ z . We can choose a gauge
where we get rid of all the σs except the first, σ z1 , 3 = W. (This is a manifestation of the
2 2fact that we had to label a configuration of the spins by the location of the domain walls
plus the value of one spin, which we are again choosing to be the first one.) Then the lastlink in the chain has an extra factor:
N − 32

Hdual = −JXzτj̄x + gτj̄z τj̄+1+ τNx − 1 + WgτN − 1 τ 1
2

j̄= 21

2

2[End of Lecture 9]

61

Gauge theory without metaphysics: Kitaev’s toric code in 1d
Since there seemed to be some discontent about the sudden appearance of gauge theory in
the previous discussion, I will try to make amends here by giving an example of whence thisconstruction can come. It is a 1d version of Kitaev’s toric code. The analogous construction
in higher dimensions is extremely important as a demonstration of the emergence of gauge
theory from a local lattice model. This silly reduction of it also does that, but gauge theory
in 1d is not quite as exciting as it is in higher dimensions.
Suppose I tell you that the following is a high-energy description of a lattice quantum
system: The Hilbert space isHbig ≡ (⊗j Hj ) ⊗ ⊗j̄ Ȟj̄
(nothing is fake, there is no redundancy, no nonsense, this is really the Hilbert space), andthe Hamiltonian is
HTC = HG + Hbig

P  z
z
as before, and the new term relative towhere Hbig = −J j̄ σj̄−1,j̄ σj̄,z j̄+1 + gτj̄z σj̄,x j̄+1 τj̄+1
the previous discussion isX
HG ≡ −Egiant
Gj̄
j̄
z
x z
with Gj̄ = σj̄−1,j̄ τj̄ σj̄,j̄+1 as before. Here Egiant is an energy scale which is a gajillion times
bigger than any other energy in the problem (Egiant  J  gJ).

A few relevant facts:[Hbig , Gj̄ ] = 0,

[Gj̄ , Gj̄ 0 ] = 0, ∀j̄, j̄ 0 .

This means that we can diagonalize Hbig by first diagonalizing HG (which is a good idea,
since its coefficient is so giant), and further, we can do that by diagonalizing each Gj̄ at thesame time.
So who are these Gj̄ s? Gj̄ is unitary (so should be thought of as enacting a transformationon Hbig ), and, in this example, G2j̄ = 1, it is an element of Z2 , so its eigenvalues are ±1. For
each Gj̄ , the eigenstate with eigenvalue +1 is the low-energy one.
Let us pause to note that this model actually has a local symmetry: the symmetry generatedby Gj̄ acts only on the dofs at or next to j̄, and this is a symmetry for any j̄. However, withthe sign of Egiant as indicated, all states on which this symmetry actually acts nontriviallyhave an enormous energy ∼ Egiant . This means that we may forget about them and focus
on the low-energy subspace of the Hilbert space
Hphys ≡ {|ψi ∈ Hbig such that Gj̄ |ψi = |ψi, ∀j} ⊂ Hbig .
You have seen this before, on the previous page: it is the Hilbert space of the Z2 gauge
theory. We have emerged a gauge theory.
62

Now let me explain better what I was saying about “choosing a gauge”. What I meant is
the following. States of Hphys correspond to orbits of the action of {Gj̄ } on Hbig .21 We canthink of the projection to the G = 1 subspace as choosing a representative of each orbit of
the action of (each) G.
Consider an (oversimplified) example of a single two-state system with G = σ z . In the
X-basis, this symmetry action has two orbits:1
1√ (| →i + | ←i) = | ↑i and √ (| →i − | ←i) = | ↓i22which are determined by the eigenvalue σ z =↑, ↓. If we know that we only want to keep thestate with eigenvalue +1 then we can label an orbit by the first term. In this case there is
only one orbit, and the projection is one-dimensional. We can describe this by saying that
we have used the gauge redundancy to set σ x = 1, and simply set σ x = 1 wherever it appears
in our Hamiltonian.
More generally, we will have some action of a (local in space) group g ∈ G of the form Ug
on our Hbig . This means that Hbig forms a (in general highly reducible) representation of G.
The orbits of G are
1 X
p
Ug |g0 i|G| g∈G
where |G| is the order of the group G, and |g0 i is a reference state (like a highest-weightstate in the theory of Lie groups). We may label the orbits by the reference state. And moreto the point, states that lie in the same orbit correspond to the same unit eigenvector, andcan be regarded as equivalent. This is the origin of the usual way of thinking about gauge
redundancy.
We’ll have more to say about the real toric code in d > 1 later (in §5.2). We have shownthat this model reduces at energies E  Egiant to the gauge theory description of the TFIM,
which we have solved by the Jordan-Wigner transformation. The situation is a little bit
different in 2d in that the model is not solvable away from the zero-correlation length limit.
21

Given a unitary operator on a Hilbert space and a state |ψ0 i in that space, what I mean by ‘the orbitof |ψ0 i under the action of U’ is the set of images under repeated action of U:
orbit of |ψ0 i under action of U ≡ {|ψ0 i, U|ψ0 i, U2 |ψ0 i...}(if U represents a finite group, this contains a finite number of elements; call this order(U)). The associatedeigenvector of U with eigenvalue 1 is1porder(U)order(U)X

Ul |ψ0 i.

l=0

It is just the average over the group. (For non-abelian groups, we must speak about the orbit of |ψ0 i under
the action of the whole group.)

63

Boundary conditions on fermions
How does this arise in terms of the fermions? Just as for domain walls, the sign of thefermion operator is also not well-defined. With periodic boundary conditions on the spins(Zj+N = Zj ), we can also define the fermions by Q
putting a branch cut at j = N , that is,
we’ll write the Jordan-Wigner formula as χj = Zj N ≥i>j Xi , where the string of Xi s stops
at N . This is a choice. The inverse map is then Y
Zj = c†j + cj
eiπni .
(46)
N ≥i>j

The ferromagnetic term in the hamiltonian, for j 6= N , is



 


Zj Zj+1 = c†j + cj eiπnj+1 c†j+1 + cj+1 = c†j + cj −c†j+1 + cj+1
where I used the identity (34). (This is as in our previous discussion.) But with periodic
boundary conditions on Z, the term at the branch cut is: Y
ZN Z1 = c†N + cN c†1 + c1
eiπnj .
N ≥j>1This last red factor includes all of the sites except the first:
YN ≥j>1

iπnj

e

=

N
Y

eiπnj eiπn1 = eiπN eiπn1 ≡ (−1)N eiπn1 .

j=1PHere N ≡ j c†j cj operator which counts the total number of fermions, and (−1)N is called
the ‘fermion parity’22 . This is a symmetry generator: [H, (−1)N ] = 0, since H only contains
terms with an even number of fermions. Returning to the outlying ferromagnetic coupling,
ZN Z1 = c†N + cN c†1 + c1 eiπn1 (−1)N
{z
}
|
†

=



  =c1 −c1c†N + cN c†1 − c1 (−1)N .

(47)

This is the same as the other terms, but for a factor of −(−1)N . This means that we getperiodic boundary conditions on c, i.e. cN +1 = +c1 if the total number of fermions is
odd and antiperiodic boundary conditions, cN +1 = −c1 , in sectors where N is even. Insummary23 :cN +1 = −c1 (−1)N.
2223

please don’t confuse the number operator N with the number of sites N – sorry about that
This fact is perhaps simpler in terms of the majoranas:
for j < N , − iχ̃j+1 χj = Zj Zj+1 ,

64

but

− iχ̃N χ1 = ZN Z1 (−S) .

So the full Hamiltonian is
N
−1 

 


X
X †
1 PBC on Z
†
†
†
†
HTFIM
=−
cj + cj cj+1 − cj+1 − cN + cN c1 − c1 (−1)N − 2gcj cj .
J
j=1
j
| {z }=N

We can diagonalize this by fourier transformation
1 X ikja
cj = √e ck
N kwhere the allowed values of k depend on the boundary conditions on c:2π{1...N }
Na 

2π 1
k∈
+ {1...N }Na 2

PBC, (−1)N = −1 :k∈APBC, (−1)N = +1 :(48)

And the Hamiltonian in each case is: XX †ck ck (2g − 2 cos ka) + (ck c−k i sin ka + h.c.) ≡
hk
H/J =
k

k

The dynamics determine the boundary conditions: is it less energy to have an odd numberof fermions or an even number? First notice that the dispersion has an inversion symmetryk → −k, so each k value is degenerate in energy with −k, unless it is fixed by this map.
The fixed points are k = 0, π/a (recall that ka = −π ' −π + 2π = π). These only occur for
PBC.
The twofold (in particular even) degeneracy for most k means that the occupation Nk6=−konly changes by two, and we don’t care about the associated hk . More algebraically, for
PBC,
X †
X †
ck ck + c†0 c0 + c†π cπN=
cj cj = 2jk6=0,π(for APBC we can simply omit the last term and conclude that N is even, so (−1)N = 1.)
So for determining (−1)N , we only care about hka=0 and hka=π .
hka=π = c†π cπ (2g + 2) > 0

empty in the groundstate, ∀g ≥ 0=⇒

As usual, the zero-momentum mode is the hero of the story:
hka=0 = c†0 c0 (2g − 2)which changes sign at g = gc = 1:
• For g > gc ,, h0 > 0 is a positive energy cost for filling the zeromode. So the groundstatefor g > gc has this odd mode empty and hence (−1)N = 1, and APBCs.

65

• For g < gc ,, h0 < 0 is an incentive for filling the zeromode, of which there is odd
number (namely, 1), so the system can lower its energy by having PBC and fillingthe k = 0 mode. In the thermodynamic limit L = N a → ∞, these two groundstates
become degenerate.

Comments on the splitting between the groundstates
Comments about the L-dependence of the splitting between these states:(1) We would like to identify them with (linear combinations of) the two ferromagnetic
groundstates, which become the fully-polarized (in Z) states |+i and |−i as g → 0. As
we’ve argued from many other points of view, these states are split by an amount which isexponentially small in L: ∆E ∼ e−L/ξ .
(2) At the critical point, where the fermions are massless, the energy splitting between the
groundstates with PBC and APBC (in the scaling limit, where a → 0), is an oft-used result
1 ~c(PBC is higher energy). The L-dependence is determined byin string theory, and is 16
Ldimensional analysis, since there is no other scale. Away from the critical point I don’t know
the answer at the moment.
(3) Note that these two sectors with different (−1)N do not mix via HTFIM , since (−1)N is
a symmetry generator; you can’t create a single fermion using HTFIM . There is no tunnelingbetween these vacua. So these two states must in fact be the cat states √12 (|+i ± |−i) which
are eigenvectors of the Z2 symmetry of the TFIM, which, after all, is
S=

N
Y
j=1

N
Y
†
Xj =
(−1)cj cj = (−1)N .
j=1(Recall that the symmetric cat state √12 (|+i + |−i) is the ground state in finite volume,so this must be the fermion state with antiperiodic boundary conditions.) Being forced tosuperpose states with different boundary conditions may be discomfiting, but it is part of
the resolution of the puzzle of an odd number of domain walls on a circle.

662.2.8

Coherent state path integrals for fermions[Shankar, Principles of QM, path integrals revisited. In this chapter of his great QM textbook, Shankar sneaks in lots of insights useful for modern condensed matter physics]
Consider a single fermion mode operator (as in our discussion of Jordan-Wigner in 0+1
dimensions);{c, c} = 0, {c† , c† } = 0, {c, c† } = 1 .
With a single mode, the general Hamiltonian is
H = c† c (ω0 − µ)(ω0 and µ are (redundant when there is only one mode) constants). We might be interested
in its thermal partition function
H
Z = tr e− T .
ω0 −µ(In this example, it happens to equal Z = 1 + e− T , as you can see by computing the tracein the eigenbasis of n = c† c. But never mind that; the one mode is a proxy for many, where
it’s not quite so easy to sum.) How do we trotterize this? That is, what is the correspondingclassical system? One answer is to use the 0d Jordan-Wigner map from §2.2.5. That routetakes us back where we came from. Here’s another answer. We can do the Trotterizing using
any resolution of the identity on H, so there can be many very-different-looking answers tothis question.
Let’s define coherent states for fermionic operators:c|ψi = ψ|ψi.

(49)

Here ψ is a c-number (not an operator), but acting twice with c we see that we must have
ψ 2 = 0. ψ is a grassmann number. These satisfy
ψ1 ψ2 = −ψ2 ψ1 , ψc = −cψ

(50)– they anticommute with each other and with fermionic operators, and commute with ordinary numbers and bosons. They seem weird but they are easy. We’ll need to consider
multiple grassmann numbers when we have more than one fermion mode, where {c1 , c2 } = 0
will require that they anticommute {ψ1 , ψ2 } = 0 (as in the definition (50)); note that we will
be simultaneously diagonalizing operators which anticommute.
It’s very easy to solve the equation (49):|ψi = |0i − ψ|1i
where as above |0i is the empty state (c|0i = 0 ) and |1i = c† |0i is the filled state. (Check:c|ψi = c|0i − cψ|1i = +ψc|1i = ψ|0i = ψ|ψi.)

67

Similarly, the left-eigenvector of the creation operator ishψ̄|c† = hψ̄|ψ̄,hψ̄| = h0| − h1|ψ̄ = h0| + ψ̄h1|.

Notice that these states are weird in that they are elements of an enlarged hilbert space withgrassmann coefficients (usually we just allow complex numbers). Also, ψ̄ is not the complex
conjugate of ψ and hψ̄| is not the adjoint of |ψi. Rather, their overlap ishψ̄|ψi = 1 + ψ̄ψ = eψ̄ψ .

Grassmann calculus summary. In the last expression we have seen an example of theamazing simplicity of Taylor’s theorem for grassmann functions:f (ψ) = f0 + f1 ψ .
Integration is just as easy and its the same as taking derivatives:ZZψdψ = 1,1dψ = 0.
With more than one grassmann we have to worry about the order:ZZ
1 = ψ̄ψdψdψ̄ = − ψ̄ψdψ̄dψ.
The only integral, really, is the gaussian integral:Ze−aψ̄ψ dψ̄dψ = a.
Many of these giveZe−ψ̄·A·ψ dψ̄dψ = det A.
 
A11 A12 · · ·
ψ1

..


. · · ·  ... Here ψ̄ · A · ψ ≡ ψ̄1 , · · · , ψ̄M A21
. One way to get this expression is
..
..
..
ψM
.
.
.
to change variables to diagonalize the matrix A.
Rhψ̄ψi ≡

ψ̄ψe−aψ̄ψ dψ̄dψ1R= − = −hψ ψ̄i.
−a
ψ̄ψ
a
e
dψ̄dψ

If for many grassman variables we use the action S =hψ̄i ψj i =Pi ai ψ̄i ψi (diagonalize A above) thenδij≡ hījiai

and Wick’s theorem here ishψ̄i ψ̄j ψk ψl i = hīlihj̄ki − hīkihj̄li.

68

(51)

Back to quantum mechanics: The resolution of 1 in this basis isZ
1 = dψ̄dψ e−ψ̄ψ |ψihψ̄|(52)

And if A is a bosonic operator (made of an even number of grassmann operators),Z
trA = dψ̄dψ e−ψ̄ψ h − ψ̄|A|ψi .
(Note the minus sign; it will lead to a deep statement.) So the partition function is:ZH
e− T
Z = dψ̄0 dψ0 e−ψ̄0 ψ0 h − ψ̄0 |
|ψ0 i
|{z}
=(1 − ∆τ H) · · · (1 − ∆τ H)

|

{z

}

M times

Now insert (52) in between each pair of Trotter factors to get
Z=

Z MY
−1

dψ̄l dψl e−ψ̄l ψl hψ̄l+1 |(1 − ∆τ H)|ψl i .

l=0

Because of the −ψ̄ in (52), to get this nice expression we had to define an extra letterψ̄M = −ψ̄0 , ψM = −ψ0(53)

so we could replace h − ψ̄0 | = hψ̄N |.

[End of Lecture 10]

Now we use the coherent state property to turn the matrix elements into grassmann-valued
functions:


∆τ →0
hψ̄l+1 | 1 − ∆τ H(c† , c) |ψl i = hψ̄l+1 | 1 − ∆τ H(ψ̄l+1 , ψl ) |ψl i = eψ̄l+1 ψl e−∆τ H(ψ̄l+1 ,ψl ) .
It was important that in H all cs were to the right of all c† s, i.e. that H was normal ordered.)
So we have
Z =

Z MY
−1

dψ̄l dψl e−ψ̄l ψl eψ̄l+1 ψl e−∆τ H(ψ̄l+1 ,ψl )
l=0



Z MY
−1




 ψ̄l+1 − ψ̄ldψ̄l dψl exp ψl − H(ψ̄l+1 , ψl )
∆τ  ∆τ| {z }
l=0=∂τ ψ̄
!
Z
Z 1/T
Z
'
[Dψ̄Dψ] expdτ ψ̄(τ ) (−∂τ − ω0 + µ) ψ(τ ) = [Dψ̄Dψ]e−S[ψ̄,ψ] . (54)
=

0

69

Points to note:• In the penultimate step we defined, as usual, continuum fields
ψ(τl = ∆τ l) ≡ ψl ,ψ̄(τl = ∆τ l) ≡ ψ̄l .

• We elided the difference H(ψ̄l+1 , ψl ) = H(ψ̄l , ψl ) + O(∆τ ) in the last expression. Thisdifference is usually negligible and sometimes helpful (an example where it’s helpful isthe discussion of the number density below).
• The APBCs (53) on ψ(τ + T1 ) = −ψ(τ ) mean that in its fourier representation24
XX
ψ(τ ) = Tψ(ω)e−iωn τ , ψ̄(τ ) = Tψ̄(ω)eiωn τ
(55)
n

n

the Matsubara frequencies
ωn = (2n + 1)πT, n ∈ Z
are half-integer multiples of πT .
• The measure [Dψ̄Dψ] is defined by this equation, just as in the bosonic path integral.
• The derivative of a grassmann function is also defined by this equation; note thatψl+1 − ψl is not ‘small’ in any sense.
• In the last step we integrated by parts, i.e. relabeled terms in the sum, so
X
X
X
X
X
X

ψ̄l (ψl − ψl−1 ) .
ψ̄l ψl = −ψ̄l0 ψl−1 −ψ̄l ψl =ψ̄l+1 − ψ̄l ψl =
ψ̄l+1 ψl −ll

l0 =l−1lllNote that no grassmanns were moved through each other in this process.
The punchline of this discussion for now (for purposes of our study of the scaling limit of
the TFIM) is that the euclidean action isZS[ψ̄, ψ] = dτ ψ̄∂τ ψ + H(ψ̄, ψ).
The first-order kinetic term we’ve found ψ̄∂τ ψ is sometimes called a ’Berry phase term’.
Note the funny-looking sign.
Continuum limit warning (about the red ' in (54)). The Berry phase term is actually
N
−1
X
l=0
24

ψ̄l+1 (ψl+1 − ψl ) = TXψ̄(ωn ) 1 − eiωn τ ψ(ωn )ωnψ̄ is still not the complex conjugate of ψ but the relative sign is convenient.

70

and in (54) we have kept only the leading nonzero term:
1 − eiωn τ → iωn τ.
Clearly this replacement is just fine ifωn τ  1for all ωn which matter. Which ωn contribute? I claim that if we use a reasonable H =Hquadratic + Hint , reasonable quantities like Z, hO† Oi, are dominated by ωn  τ −1 .
There’s more we can learn from what we’ve done here that I don’t want to pass up. Let’s
use this formalism to compute the fermion density at T = 0:
1 −H/T †trec c.
ZThis is an example where the annoying ∆τ s in the path integral not only matter, but are
extremely friendly to us.
hNi =Frequency space, T → 0.
Let’s change variables to frequency-space fields, which diagonalize S. The Jacobian is 1
(since fourier transform is unitary):YT →0dψ̄(ωn )dψ(ωn ) → Dψ̄(ω)Dψ(ω).
Dψ̄(τ )Dψ(τ ) =n

The partition function is!
ZXψ̄(ωn ) (iωn − ω0 + µ) ψ(ωn ) .
Z = Dψ̄(ω)Dψ(ω) exp TωnNotice that in the zero-temperature limit
Z
Z
X
dω
T
7→
≡ d̄ω.
2π
ωnR
P(This is the same fact as V k 7→ d̄d k in the thermodynamic limit.) So the zerotemperature partition function is
Z ∞

Z
T →0
Z =
Dψ̄(ω)Dψ(ω) expd̄ω ψ̄(ω) (iω − ω0 + µ) ψ(ω) .
−∞Using the gaussian-integral formula (51) you can see that the propagator for ψ ishψ̄(ω1 )ψ(ω2 )i =

δω1 ,ω2
T }
| {z

T →0

→ δ(ω1 −ω2 )

71

2π
.
iω1 − ω0 + µ

(56)

2π/T
In particular hψ̄(ω)ψ(ω)i = iω−ω. δ(ω = 0) = 1/T is the ‘volume’ of the time direction.
0 +µBack to the number density. Using the same strategy as above, we have
1
hNi =
Z

Z MY
−1+1 

dψ̄l dψl e−ψ̄l ψl

l=0

−1
 MY
hψ̄l+1 |(1 − ∆τ H(c† c))|ψl il=1hψ̄N +1 |c† c|ψN i
{z
}
|

.

=ψ̄N +1 ψN =ψ̄(τN +∆τ )ψ(τN )This formula has a built-in point-splitting of the operators!
Z
1
hNi =
Dψ̄Dψ e−S[ψ̄,ψ] ψ̄(τN + ∆τ )ψ(τN )
Z
Z ∞
eiω∆τ
d̄ω
=
= θ(µ − ω0 ).
iω − ω0 + µ
−∞

(57)

Which is the right answer: the mode is occupied in the groundstate only if ω0 < µ. In
the last step we used the fact that ∆τ > 0 to close the contour in the UHP; so we only
pick up the pole if it is in the UHP. NoticeR Λ dω that this quantity is very UV sensitive: if we∼ log Λ, the integral diverges logarithmically.
put a frequency cutoff on the integral,ω
For most calculations the ∆τ can be ignored, but here it told us the right way to treat thedivergence. 25
Bosons. This is a good place to mention quickly that we can do the same thing for bosons,using ordinary SHO (simple harmonic oscillator) coherent states. What I mean by ‘bosons’is a many-body system whose Hilbert space can be written as H = ⊗k Hk where k is a label(could be real space, could be momentum space) and Hk = span{|0ik , a†k |0ik , 21 a†kis the SHO Hilbert space.

2|0ik , ...}For each normal mode a, coherent states are
a|φi = φ|φi;†|φi = N eφa |0i.

The eigenbra of a is hφ|, with?

hφ|a = hφ|φ? , hφ| = h0|e+φ a N .
25

The calculation between the first and second lines of (57) is familiar to field theorists – it can be described
as a feynman diagram with one line between the two insertions. More prosaically, it isZXeiωn ∆τeiω∆τ(55) 2 X i(ωn −ωm )τ +iωn ∆τ
(56)
T →0
hψ̄(τN + ∆τ )ψ(τN )i = Tehψ̄(ωn )ψ(ωm )i = T→d̄ω.
iωn − ω0 + µiω − ω0 + µ
nm
m

72

Their overlap is26 :?

hφ1 |φ2 i = eφ1 φ2 .
2

If we choose N = e−|φ| /2 , they are normalized, but it is more convenient to set N = 1. The
overcompleteness relation on Hk isZ
dφdφ? −|φ|21k =e|φihφ|.
π
It will be convenient to arrange all our operators into sums of normal-ordered operators:: ak a†l :=: a†l ak := a†l ak
with all annihilation operators to the right of all creation operators. Coherent state expectation values of such operators can be built from the monomialshφ|Y  † MkYak(ak )Nk |φi =(φ?k )Mk (φk )Nk .
kk

Let the Hamiltonian be H = H({a†k }, {ak }) =: H :, normal ordered. By now you know
how to derive the path integral using this resolution of the identity:
Z = tre−H/T
Z
N
Y
PN
?
=
dφl e− l=0 (φl+1 (φl+1 −φl )−∆τ H(φl+1 φl ))ZφN +1 =φ0 l=0
R 1/T
?
?
'
[Dφ] e− 0 dτ (φ ∂τ φ+H(φ ,φ)) .

(58)

φ(0)=φ(1/T )Real time. If you are interested in real-time propagation, rather than euclidean time, justreplace the euclidean propagator e−τ H 7→ e−itH . The result, for example, for the amplitude
to propagate from one bose coherent state to another is
hφf , tf |e−itH

Z φ(tf )=φf|ψ0 , t0 i =

i

R tf

?

?

Dφ? Dφ e ~ t0 dt(i~φ ∂t φ−H(φ,φ )) .

φ(t0 )=φ0

Note the distinguishing feature of the Berry phase term that it produces a complex term in
the real-time action.

26

Check this by expanding the coherent states in the number basis and doing the integrals
Z
Z 2π
Z
n+n0
dφdφ? −φφ? n ? n0dθ i(n−n0 )θ ∞e
φ (φ ) =e
due−u u 2
π
2π
0
0
P
to get 1 = n |nihn|.

73

2.2.9

Scaling near the fixed-point theory of the Ising phase transition.

Armed with the fermion coherent state path integral, let’s go back to the continuum hamiltonian describing the TFIM near g = 1:
Z
Z
Z

c
†
†
†
H
dx Ψ(x) ∂x Ψ(x) − Ψ(x)∂x Ψ(x) + ∆ dxΨ Ψ ≡ dx h2with ∆ = 2J|g − 1|.
The euclidean-time action is S[Ψ, Ψ† ] =

R

dτ

R

dxL, with

L = Ψ̄∂τ Ψ + hc= Ψ̄∂τ Ψ +Ψ̄∂x Ψ̄ − Ψ∂x Ψ + ∆Ψ̄Ψ.
2(59)

Here I have simply replaced the field operators Ψ(x), Ψ† (x) with their (right and left, respectively) grassmann eigenvalues Ψ(x), Ψ̄. This lagrangian is more appealing in terms of
the self-conjugate variables ψ± :
Ψ = ψ+ + ψ− + i (ψ+ − ψ− ) , Ψ̄ = ψ+ + ψ− − i (ψ+ − ψ− )XL=
ψ± (∂τ ± i∂x ) ψ± + ∆ψ+ ψ− .
±If we let z ≡ x + iτ , and turn off the mass perturbation, this is the third avatar of the Isingmodel I mentioned at the beginning in (6).
Let’s analyze its behavior under scale transformations. In order to make S invariant, wemust scale time and space the same way (z = 1), and we must scale1

x → λx, t → λt, Ψ → λ− 2 Ψ.
So just by scaling, we can see that its correlators behave ashΨ(x)† Ψ(0)i ∼1.
x(Since it’s a gaussian theory, you can also calculate the coefficient.) More generally, thescaling dimension δ of an operator O(x) in O → λ−δ O in a scale-invariant QFT determines
its vacuum autocorrelation functions to behO(x)† O(0)i ∼

1
.
x2δ

(Later we will interpret δ as an eigenvalue of the dilatation operator.)
The mass perturbation ∆ will violate the scale invariance, sinceZZ+1dxdτ Ψ̄Ψ → λdxdτ Ψ̄Ψ .
74

As we make lengths and times bigger, it gets bigger – it’s a relevant perturbation, which determines the behavior at long distances. The scaling dimension of this relevant perturbation
determines the correlation length critical exponent: the correlation length is a length scaleξ = aλ at which the relevant coupling ∆ = λν ∆0 ∼ 1 has evolved to an order-one value (∆0
is its ‘bare’ value, at the lattice scale a, and ν is its scaling dimension). Eliminating λ = ξ/afrom these equations, we get
1
1 = (ξ/a)ν ∆0 =⇒ ξ ∼ ν ,∆0so ν, the scaling dimension of the relevant operator, is indeed the correlation length criticalexponent defined above.
Let’s ask whether there are any other relevant operators, in the sense that they scale tolarger values at larger λ (longer wavelengths). If there are more it means that our fixed
point is multicritical – we would have to tune more than one parameter to reach it. We’regoing to demand that the Ising symmetry is preserved; this is the fermion number symmetry,
Ψ(x) → −Ψ(x). So we should only consider operators made of an even number of fermions.
And all the other operators we can make out of an even number of Ψ are irrelevant, in thesense that they scale to smaller values at larger λ. The marginal-looking operator Ψ̄ΨΨ̄Ψ
vanishes by Fermi statistics. Other operators have more derivatives have no chance:Zdxdτ Ψ̄∂x2 Ψ ∼ λ−1Zdxdτ Ψ̄∂x ΨΨ̄∂x Ψ ∼ λ−2This means we have to tune only one parameter (g − gc ∼ ∆) to get to the critical theory.
(Since there are no classically-marginal operators, there is no need to do any perturbating
to determine the effects of turning them on a little bit – small quantum corrections cannot
overcome the order-1 amounts by which these operators are irrelevant.)
There are operators of dimension less than two, namely the operator Z itself; this is anoperator which creates a branch cut in the fermions, which after all are made from the domain
wall operators. To understand its dimension, it is useful to use a bosonized description, which
we will do in the next section. Its scaling dimension turns out to be 1/8. This is importantbecause it determines the rate at which the magnetization vanishes as g → gc from below.
It’s important that the Z2 symmetry forbids us from adding this (relevant) operator to theaction. So in the presence of the symmetry, the Ising fixed point has only one relevant
perturbation (the fermion mass).
Notice how different this critical theory is from the guess would have made from the mean
field theory above in §2.2.2: we would have guessed that the critical fluctuations should be
those of the mean field θ, which when I set g = gc becomes a massless free boson, with many
relevant operators which preserve the θ → −θ symmetry, since demanding scale invarianceofZ
S0 [θ] = dxdτ (∂τ θ)2 + (∂x θ)275

implies that θ is dimensionlesss. So all even powers θ2 , θ4 ... (in fact, all even functions, likecos θ) have a chance to be relevant. We’ll see that it’s not actually quite so far off. (And infact is the correct variables in d > 1, and the correct fixed point for d ≥ 3 (the upper criticaldimension).)
The nontrivial scaling of Z is an example of an anomalous dimension: the engineering
dimension is different from the correct scaling behavior. Dimensional analysis is inviolable.
What has happened is that powers of the short-distance scales are involved; that is: it is J 1/8 Z(or a−1/8 Z) that has finite correlation functions in the scaling limit. In relativistic QFT thisphenomenon is called ‘wavefunction renormalization’ – we must multiplicatively renormalize
our operators to get finite (not zero, not infinity) correlators in the long-wavelength limit.

2.2.10Beyond the quantum Ising chain

The Jordan-Wigner transformation is useful in spin chains with many other choices of Hamiltonian and symmetry group, as you’ll see on the problem set, and below.
Here we have a choice to make. Do we study systems which have more symmetry at eachsite or do we study ising spins in more dimensions? We’ll get to both eventually.
[End of Lecture 11]763

Bosons and Bosonization

I mentioned that another realization of two-state systems at each site comes from self-loathingbosons. Let’s see how that comes about in a model with U(1) symmetry.

3.0.1

Spins, bosons, fermionsConsider the following spin system (qbits at the sites of a lattice)µX
wX(Xi Xj + Yi Yj ) +Zj .
HXX = −
2
2 jhiji

On the homework, you Jordan-Wignered the 1d version of this27 and foundX  †
w cj+1 cj + h.c. + µc†j cj
HXX = −j(Note that µ is playing the role of the transverse-field term.) The fact that we have equal‘ferromagnetic’ couplings involving Xs and Ys means that this system preserves a continuous
U(1) symmetry, whose unitary implementation isY θUθ =ei 2 Zj .
j

It acts on the spin operators by
Xj 7→ Uθ Xj U†θ = cos θXj + sin θYj ,Yj 7→ Uθ Yj U†θ = cos θYj − sin θXj , Zj 7→ Zj .

(Please beware my signs.) It therefore commutes withof the fermions it acts
PHXX . In
Pterms†iθ
by c → Pe c: it is the total fermion particle number j nj ≡ j cj cj that is conserved as aresult, [ j nj , H] = 0 – not just mod two like when we had only Z2 symmetry.
The fermion description is nice because it is quadratic; the fermions are free in this subspaceof the phase diagram. If we add the nearest-neighbor Zi Zj term, we add
X
X∆H = −v
Zi Zj = −v
ni nj ,hiji

hiji

which is quartic in the fermion operators, and therefore not simple. A description in termsof bosons will allow us to understand this direction in the phase diagram.
27

using the (rotated X → Z relative to our previous one) map




Y
Y
†
Zj = 1 − 2c†j cj , Xj = − (1 − 2c†i ci ) cj + c†j = − (−1)ci ci cj + c†j , Yj = −iZj Xj
i>j

i>j

77

An explicit map from spin- 12 s to hardcore bosons is:
Xj = bj + b†j , Yj = −i bj − b†j ,Zj = 1 − 2b†j bj .

(60)

(The inverse is b = σ + ). Unlike Jordan-Wigner, we don’t have to attach any kind of string toanybody to make this work: before the hardcore projection, the bosons commute at different
sites[b̂i , b̂†j ] = δij .
The price of that niceness is that the map (60) relates the qbit hilbert space to the subspaceof the (much bigger!) boson Hilbert space where the occupation numbers b†j bj ≡ nbj = 0 or
1. And the objects in (60) only satisfy the Pauli algebra on this subspace (check it!). Worse,
the boson operators, projected to the no-double-occupation (b†j bj ≤ 1∀j) subspace are nolonger canonically conjugate. The projected bosonsb ≡ Pb̂P, b† ≡ Pb̂† P,

P ≡ |0ih0| + |1ih1|

instead satisfy[b, b† ] = 1 − 2b† b.

(61)

Check it by acting with the both-hand side on |ni with n = 0, 1. (61) in fact says that the
projected hardcore bosons at a single site anticommute to 1 . But for different sites they
commute:[bi , b†j ] = δij 1 − 2b†j bj .

(60) is not the only way to map spins to bosons. The bosons above are closely related to
Holstein-Primakoff bosons ([h, h† ] = 1), which satisfy
σ+ ≡

p
p
†
1
(X + iY) = 2S − h† hh , σ − = σ + = h† 2S − h† h,2

Z = 2(S − h† h), (62)with the (non-holonomic) constraint 0 ≤ h† h ≤ 2S. S is the spin – 2S + 1 is the number ofstates per site. These have the virtue that within the constrained subspace, the canonicalboson commutators [h, h† ] = 1 imply the Pauli spin algebra. The price is the sqrt factors.
A different map which has no sqrts or weird commutators, and makes the SU(2) action on
the spins manifest, is Schwinger bosons:
σ + = a† b,11 †(X − iY) = b† a, Z =a − b† .
2
2

We get spin S if we restrict to the sector with a† a + b† b = 2S (which is a nice holonomic
constraint that we can impose with a lagrange multiplier). If I don’t satisfy your thirst forlearning about this, see Dan Arovas’ notes here.

78

Persevering nevertheless with our hardcore projected bosons, let’s write the hamiltonian
above in these variables:

X †
X †
HXX = −w
bi bj + b†j bi − µbj bj .
j

hiji

P
P
A lovely hopping term and a chemical potential for the boson number, j nbj = j b†j bj .
We can make the projection to no-double-occupancy more palatable by implementing it
energetically (like we did the gauge constraint in the toric code construction above):

X †
X
X

HB = −wb̂i b̂j + b̂†j b̂i − µnbj + Unbj nbj − 1 . .
j

hiji

j

where B is for ‘big’ or ‘boson’. This last term is zero if nbj = 0, 1, but exacts a huge energetic
penalty ∆E = U if a site is occupied by two bosons. So the low-energy physics of HB is
that of HXX , with its slightly weird bosons.

3.0.2

Using the JW solution to learn about the bosons

The free fermion solution is great when it is available, but it doesn’t help much away fromthe solvable subspace of the parameters. Let’s see what it says about the bosons and then
we’ll move away from the solvable subspace.

HXX =

X

k c†k ck ,

k = −2w cos ka − µ.

k

The boson number is the same as the fermion number, whose groundstate expectation value
we can read off from the dispersion k by counting negative-energy states:µ ≤ −2w (I)
0,

1
†
†
µ
1
hbj bj i = (1 − hZj i) = hcj cj i = 1 − π cos−1 2w ,|µ| ≤ 2w (II)21,
µ ≥ 2w(III)
The groundstates at |µ| > 2w are fully-polarizedstates of the spins – the spins all point up or all
point down in the Z-basis (like g → ∞ of the TFIM).
This is a paramagnetic state where the U(1) symmetry is preserved, since it rotates the spins about theZ axis. There is no long-range order, because there is
no Z → −Z symmetry: the spins are told which way
to point. In terms of the bosons, µ < −2w is a state
where there are no bosons at all, a completely trivial insulator. The state at µ > 2w is moreinteresting. Here there is exactly one boson per site: this is a Mott insulator: the bosons canmove in principle, but do not because they would have to go near their despised neighbor.
79

These states have a gap. In terms of the fermions, this is because the band is completely
empty (µ < −2w) or completely filled (µ > 2w).

In the intermediate range of µ, we have a compressible superfluid state. ‘Compressible’ means that as we
6= 0)vary µ a little, the density changes a little, ( ∂n∂µas you can see from the figure. ‘Superfluid’ is a little
harder to explain. Note that although I say this isa superfluid, there is not in fact long-range order inthis case, because of we have continuous symmetry
in d < 2: the spins can’t decide in which direction in the XY-plane to point, so in fact
hbi = 12 hX + iYi = 0. This is a consequence of the Hohenberg-Mermin-Wagner-Coleman
theorem which I hope to discuss more later. One thing you can see clearly, however, is thatthe state is gapless, since
P there is a Fermi surface. In fact, a finite compressibility requires
gaplessness: since [H, j nj ] = 0, changing µ does not change the eigenvectors of H, only
the eigenvalues; so only if there is an unoccupied mode just above the chemical potential
can we change the density by varying µ a little.
In summary, we have:+|gs, Ii = | ↑↑↑ · · · ↑i =

=|

i

=|

i

+
|gs, IIi=+|gs, IIIi = | ↓↓↓ · · · ↓i =(63)

In region II we don’t have a simple description of the states in terms of the spin or thebosonic variables so far.
The points µ/2w = ±1 therefore describe Mott-superfluid and superfluid-Mott transitions(in terms of the fermions, they are band-insulator to metal transitions). The boson density
demonstrates that there is critical behavior: near the µ/2w = −1 transition, it is(0,µ < −2w

hnbj i '
µ 1/2
1+
,
µ>
∼ −2w
2w

80

3.0.3

Continuum description

The continuum description near the critical point µ + 2w ∼ 0
is


Z
2
†
†∇
Ψ + µF Ψ Ψ
H = dx Ψ
2mwhere the fermion chemical potential is µF = µ+2w. The U(1)symmetry acts as Ψ → eiθ Ψ. The inertial fermion mass here1
is mI = 2wa
2 . On the other hand, the ‘rest mass’ (the energy to create a fermion at rest) is
determined by µF (it’s equal to −µF for µF < 0).
As we did for the TFIM critical theory, let’s study the scaling of possible perturbations of
this critical theory. The euclidean action for the critical theory (µF = 0) isZS[Ψ̄, Ψ] =dτ dx Ψ̄∂τ Ψ + h
 
 
Z
∇2
Ψ .
(64)
=
dτ dx Ψ̄ ∂τ +2m
Notice that we have to scale time and space differently to make S invariant: (x, τ ) →(λx, λ2 τ ) so that both terms scale the same way. Thus, this fixed point has dynamical
1
exponentRz = 2. The fermion scaling is Ψ → λ− 2 Ψ (again). Notice that z 6= 1 makes themeasure dτ dx scale like λd+z = λ3 – it’s as if we were in higher dimensions. A marginal
operator has dimension d + z = 3.
The operator Ψ̄Ψ changes µF and has dimension 1; this is the relevant operator we haveto tune to get to the fixed point.
The operator Ψ̄∂x Ψ has dimension 2 and seems to be relevant. However: (1) it breaks theparity symmetry x → −x of the fixed point and so setting its coefficient to zero is technicallynatural. More importantly, (2) turning it on will merely move the minimum of the band
away from ka = π – it will not destroy the fact that the band has a minimum! So this justrenormalizes the actual critical value of µF , which was non-universal anyway.
There are no other relevant operators. The four-fermion operator Ψ† ΨΨ† Ψ = 0 again. Thismeans that the leading four-fermion operator coming from Zj Zj is actually Ψ† ∂x ΨΨ† ∂x Ψ
which is irrelevant.

81

3.0.4

Fermi liquid in d = 1 + 1.

Now consider what happens for µF > 0 (|µ| < 2w). The groundstate in this regime (returning
for a moment to the microscopic variables) isY †ck |Ωi, (ck |Ωi ≡ 0, ∀k)|k|≤kF

With k = −J cos ka, the Fermi surface is two points. For µF  w, we can use the quadratic
dispersion at the bottom√ of the band and the lowest-energy excitations above the groundstate
occur at k = ±kF = ± 2mµ. Let’s re-expand our continuum field to focus on these modes:
Ψ(x, τ ) = eikF x R(x, τ ) + e−ikF x L(x, τ )(65)where the names arise from the direction of propagation of the modes.
Z
Z

†
†
HF = dx R (−ivF ∂x ) R + L (ivF ∂x ) L ≡ dx hF
with vF = ∂k k = 4Ja is the Fermi velocity. To do a scaling analysis it will help to use theaction. The real-time Lagrangian isZLF =dx iL† ∂t L + iR† ∂t R − hF
ZZSF = i dt dx R† (∂t − vF ∂x ) R + L† (∂t + vF ∂x ) L(66)

Now there are no symmetry-preserving (U(1) and translations) relevant operators that are
not already present. Again the dimension of L, R is ∆L = ∆R = 21 (Invariance of LF under1(x, t) → λ(x, t) requires (L, R) → λ− 2 (L, R).)R† R) is relevant (∆ = 1), but merely shifts the location of the fermi• δL = µL† L (or P
surface (this is j Zj )• δL = L† ∂x L is marginal (∆ = 2), but merely renormalizes the Fermi velocity.
• δL = L† ∂x2 L is the first correction that restores the microscopic bandstructure. It andall the others are irrelevant.
†• δL
R has nonzero momentum – it violates translation invariance. This is
P = gL
j
(−1)
Z
j.
j

• δL = ∆eiϕ LR + h.c. explicitly breaks the U(1) P
symmetry and is therefore forbidden–
this is a superconducting perturbation. This is j eiϕ (Xj + iYy ) + h.c..

82• δL = vL† LR† R has dimension ∆ = P
2 and so is marginal. So we’ll have to think harderabout this one. This comes from v j λZj Zj+1 in the spin hamiltonian. It turns outthat this coupling is exactly marginal – i.e. there is a line of fixed points parametrized
by v. The parameter is called the ‘Luttinger liquid parameter’ and we will understandit through bosonization.
• δL = L† ∂x R∂x L† R has scaling dimension 3 and is therefore irrelevant. Anything withmore derivatives or more fermions is even more irrelevant.
Many features of this discussion extend to Fermi surfaces in d > 1 + 1, in particular thefact that the relevant operators merely shift kF and vF . This is the origin of Fermi liquid
theory. The analysis of the 4-fermion terms is even more interesting in that case.
[End of Lecture 12]83

Chiral anomaly
Notice that it seems that we have enhanced the particle-number symmetry to separate
conservation of R† R and L† L.
But notice, too, that we would be overcounting degrees of freedom if we imagined that R, L
had all possible momenta: they are connected by the bottom of the band. This seemingly
metaphysical consideration has an important consequence: if we apply an electric field, weconvert R excitations into L excitations. This means that the axial charge QA ≡ R† R − L† L
is not a conserved quantity. This is called the chiral anomaly.
See my penultimate lecture of Physics 215C for more on this. I reproduce part of the
discussion here:
The physics of the chiral anomaly[Polyakov, page 102; Kaplan 0912.2560 §2.1] Consider non-relativistic free (i.e. no 4-fermion
1 ~2
k . This isinteractions) fermions in 1+1 dimensions, e.g. with 1-particle dispersion ωk = 2m
a convenient description of our lattice model if µ  t, so we only care about the bottom of
the band, but it is not crucial. The groundstate of N such fermions is described by filling
the N lowest-energy single particle levels, up the Fermi momentum: |k| ≤ kF are filled. We
keep our infrared regulator so that the levels are discrete – put them in a box of length L,
so that kn = 2πn. (In Figure 1, the red circles are possible 1-particle states, and the greenLones are the occupied ones.) The lowest-energy excitations of this groundstate come fromtaking a fermion just below the Fermi level |k1 | <∼ kF and putting it just above |k2 | >∼ kF ;
the energy cost is
Ek1 −k2 =1
kF1(kF + k1 )2 −(kF − k2 )2 '(k1 − k2 )2m2mm

– we get relativistic dispersion with velocity vF = kmF . The fields near these Fermi points in
k-space satisfy the Dirac equation(ω − δk) ψL = 0, (ω + δk) ψR = 0.

84

It would therefore seem to imply a conserved axial
current – the number of left moving fermions minus
the number of right moving fermions. But the fields
ψL and ψR are not independent; with high-enoughenergy excitations, you reach the bottom of the band(near k = 0 here) and you can’t tell the difference.
This means that the numbers are not separately conserved.
We can do better in this 1+1d example and showthat the amount by which the axial current is violated
is given by the anomaly formula. Consider subjecting
our poor 1+1d free fermions to an electric field Ex (t)which is constant in space and slowly varies in time.
Suppose we gradually turn it on and then turn it off;
here gradually means slowly enough that the process
is adiabatic. Then each particle experiences a force∂t p = eEx and its net change in momentum isZ∆p = e dtEx (t).

Figure 1:

Green dots represent occupied 1-particle states. Top: In the
groundstate. Bottom: After applying
Ex (t).

This means that the electric field puts the fermions in a state where the Fermi surface k = kF
has shifted to the right by ∆p, as in the figure. Notice that the total number of fermions isof course the same – charge is conserved.
Now consider the point of view of the low-energy theory at the Fermi points. This theory
has the actionZ
S[ψ] = dxdtψ̄ (iγ µ ∂µ ) ψ ,(67)
where γ µ are 2 × 2 and the upper/lower component of ψ creates fermions near the left/rightFermi point. In the process above, we have added ∆NR right-moving particles and taken
away ∆NL left-moving particles, that is added ∆NL left-moving holes (aka anti-particles).
The axial charge of the state has changed by
Z
Z
Z
L
L
e
e
∆p
= ∆p = e dtEx (t) =dtdxEx =µν F µν
∆QA = ∆(NL − NR ) = 2
2π/L
π
π
π
2π
R
On the other hand, the LHS is ∆QA = ∂ µ JµA . We can infer a local version of this equationby letting E vary slowly in space as well, and we conclude that∂µ JAµ =eµν F µν .
2πThis agrees exactly with the anomaly equation in D = 1 + 1 produced by the continuum
QFT calculation of the anomaly using (67).

85

3.0.5

Bosonization, part 1: counting[Sachdev, chapter 20] Let’s quantize this systemon a circle of length L with APBC: L(x + L) =−L(x), R(x + L) = −R(x). The mode expansion is:
2πlxi
1 X
R(x) = √Rl e LL l∈Z+ 1
2

and a similar expression for L. We’ll focus on R fora while. The modes satisfy{Rl , R†l0 } = δll0 .
and the hamiltonian (minus the vacuum energy) is
2πvF Xl R†l Rl −E0
HR =| {z }
L
1l∈Z+ 2≡nRlso HR |gsi = 0. (The full H is H = HR + HL .) The groundstate has all modes with l > 0
empty and all modes with l < 0 filled. Notice that the Hamiltonian is symmetric under
the interchange l → −l: removing a particle with negative momentum −|l| (or energy) adds
momentum |l|. The charge operator isX
QR =: R†l Rl :∈ Z , QR |gsi = 0, [QR , H] = 0.
lLet’s count states. The thermal partition sum for the rightmover is
ZR (T ) = trHR e−HR /TY
Y X
1 2πvFR
e− T L |l|nl ≡
=
=

1+q

R

q |l|nl

2πvF(q ≡ e− T L )l∈Z+ 12 nR
l =0,1l∈Z+ 21 nR
l =0,1

Y

X


|l|

l∈Z+ 12

=

∞ 
Y

1

1 + q n− 2

2

(68)

n=1

This is an infinite-product representationof an elliptic theta function.
Another way to arrive at this function is to keep track of theparticle-hole excitations. If we take a particle out of the Fermi sea
from a filled level l0 and put it in an empty level l0 + n, we create an
excitation with momentum l0 + n − (l0 ) = n. n is an integer. Thechange in energy is∆ = vF ∆k =

2πvF
n.
L
86

Notice that it is independent of l0 ; this is because of the linear dispersion near the Fermi
surface. The two fermions (the particle and the hole) move at the same speed, and in thesame direction (because we are in 1d and they are both right-moving!) so it might be agood idea to think of them as a single particle. The operator which
this happen is
P makes†††
Rl0 +n Rl0 . It is a bosonic operator. More precisely, the object ρn ≡ l0 Rl0 +n Rl0 is a bosonicoperator which accomplishes the above operation and can act many times without giving) and the energy by 2πnv
zero. It raises the momentum by n (in units of 2π
L
L F
An arbitrary fermion state is|F i =

YR†l

Y

Rl0 |gsi.

l0 <0

l>0

If it has charge QR , we can rearrange this as|F i = (particle-hole excitations) |QR i
where |QR i is the lowest-energy state with charge QR . This state has only ls of the same
sign, and only the first QR of them. For QR > 0 it is|QR i = R†Q + 1 · · · R†3 R†1 |gsi.
R

2

2

2

(For QR < 0, remove the daggers and reverse the signs.) It has energy (for either sign of
QR )
|QR |
2πvF XπvF 2
Q .
E0 (QR ) =
l=
L
L R
1
l= 2

For simplicity, let’s consider particle-hole excitations
above the groundstate, |gsi = |QR = 0i, as in the figure at right. (To understand what happens for QR 6= 0,
just shift all the labels by QR .) The state in the figure
is|63331i =






R†− 1 R− 7
R†1 R− 5R†3 R− 3R†11 R− 1 |gsi
R†− 7 R− 9
2

2

2

2

2

2

2

2

2

It is important that we act first (first means its at theright of the queue) with the boson operator making
the largest jump – this guarantees that we never try toannihilate a fermion that is not there!
So we may label the excitation above |QR i by a set of
numbers rn indicating how many steps up the fermion
in box n = −l − 21 was pushed:
r1 ≥ r2 ≥ r3 ... ≥ rn? = 0.
87

2The corresponding state is:ρ†rn? · · · ρ†r2 ρ†r1 |QR i
This information can be represented by a Young diagram: For example, we represent thestate in the figure, which has {rn } = {6, 3, 3, 3, 1}, by:E, QR = 0 = ρ†1 ρ†3 ρ†3 ρ†3 ρ†6 |QR = 0i
Notice that some of the rn may be equal; the state still makes sense because these operators
are bosons.
We can see how the Young diagram encodes the state interms of the fermions, if we do anticommutator algebra
to cancel R†l and Rl whenever they both appear:E=
R†− 7 R− 9 · R†− 1 R− 7 · R†1 R− 5 · R†3 R− 3 · R†11 R− 1 |gsi
2
2
2
2
2
2
2
2
2
2


= R†11 R 3 † R 1 † R− 3 R− 5 R− 9 |gsi
2

2

2

2

2

2

Look at the momentum indices on the annihilation andcreation operators. You can see these very numbers are
the (±) the numbers of boxes to the right and below thediagonal of the diagram respectively.
What is the hamiltonian in terms these boson operators?
2πvF
πvF 2
QR +
H =
L
L
B

∞
X

nρ†n ρn .

n=l+ 12 = 21

Figure 2:

This is a chiral boson, since n > 0. Its partition sum is
B

R
Zboson(T ) = trHB e−H /T !
∞∞XYQ2R /2=q

The two ways of looking at the Young diagram which
determine the boson state and thefermion state, respectively.

X

q nmn

(69)

n=1 mn =ρ† ρn =0,1,..∞nQR =−∞

|

{z

1
1−q n

}

This is an infinite-sum representation of the same elliptic theta function. (For more about
the theta functions, if you must, see e.g. Polchinski vol. 1, page 215.)

88

3.1

Free boson CFT in 1+1d, aka Luttinger liquid theory

Now we take an apparent detour. One way to motivate what we’ll do is to imagine that oursystem tries to produce by spontaneous symmetry breaking a Goldstone boson for our U(1)symmetry. We will study the physics of the resulting scalar field – whose field space is thebroken group U(1), i.e. the circle – and see that the symmetry is in fact restored. We’ll learn
a lot of other things, too.
Consider the following (real-time) continuum action for a free massless scalar field in 1+1
dimensions:
Z
Z L
Z
T
2
2
dt
dx (∂0 φ) − (∂x φ) = 2T dxdt∂+ φ∂− φ .
(70)
S[φ] =
2
0
I have set the velocity of the bosons to c = 1 by rescaling t. Here x± ≡ t ± x are lightcone
coordinates; the derivatives are ∂± ≡ 12 (∂t ± ∂x ). Space is a circle: the point labelled x is
the same as the point labelled x + L.
We will assume that the field space of φ itself is periodic:
φ(x, t) ≡ φ(x, t) + 2π,∀x, t .

So the field space is a circle S 1 with (angular) coordinate φ. It can be useful to think of
the action (70) as describing the propagation of a string, since a field configuration describes
ana embedding of the real two dimensional space into the target space, which here is a circle.
This is easier to picture if we consider more than one (say N ) scalar fields, φa , which we cantake to be coordinates on an N -dimensional space with metric gab (φ). Then we can consider
the more general actionZa
S[φ ] = T dxdt∂+ φa ∂− φb gab (φ(x, t)) .
This is called a non-linear sigma model. Note that ‘T ’ is for ‘tension’28 .
[End of Lecture 13]
Let’s retreat now to the case where the target space is a circle. To get back the action (70),
we should set gφφ = 2. Notice that by changing coordinates on the field space, e.g. φ̃ ≡ Rφwith α independent of x, t, we change the prefactor in the action:
2T ∂+ φ∂− φ =28

2T∂+ φ̃∂− φ.
R2

As you might expect, there are many conventions in the literature for this quantity. The string theory1
literature is commonly written in terms of α0 ≡ 2πT. The literature on Luttinger liquids (e.g. in Sachdev’s
10book) uses K ≡ 2α = πT .

89√For example, by setting R2 = 2T it looks like we can make it disappear. But since φ iscompact, it returns in the periodicity:φ̃ ≡ φ̃ + 2πR.
So R is like the radius of the circle. This leads to a very annoying multiplicity of conventions.
Noticethat if φ is not compact, the path integral gets a very awkward infinite part fromR∞dφ
0 (with φ0 constant). But these are just the same issues that arise in the path−∞
integral of a particle on the real line; this integral will produce delta functions which enforce
conservation of momentum in the φ-space.
First a little bit of classical field theory. The equations of motion for φ are0=δS∝ ∂ µ ∂µ φ ∝ ∂+ ∂− φ
δφ(x, t)which is solved byφ(x, t) ≡ φL (x+ ) + φR (x− ) .
In euclidean time, φL,R depend (anti-)holomorphically on the complex coordinate z ≡ 12 (x +iτ ) and the machinery of complex analysis becomes useful.

Symmetries: Since S[φ] only depends on φ through its derivatives, there is a simple
symmetry φ → φ + . By the Nöther method29 the associated current is
jµ = T ∂µ φ .

(71)This symmetry is translations in the target space, and I will sometimes call the associatedconserved charge ‘momentum’.
There is another symmetry which is less obvious. It comes about because of the topologyof the target space. Since φ(x, t) ≡ φ(x, t) + 2πm, m ∈ Z describe the same point (itis a redundancy in our description, in fact a discrete gauge redundancy), we don’t needφ(x + L, t) = φ(x, t). It is enough to haveφ(x + L, t) = φ(x, t) + 2πm, m ∈ Z
The number m cannot change without the string breaking: it is a topological charge, awinding number:
Z L
1
x=L FTC 1
φ(x, t)|x=0 =dx∂x φ .
(72)
m=
2π
2π 0
29
In case you have forgotten, here’s the method: if the transformation with constant  is a symmetry, thenthe variation of the action with  = (x, t) must be proportional to ∂µ , and so that it vanishes ∀φ when is constant:
Z
Z
IBP
S[φ + (x)] − S[φ] = dxdt∂µ (x)j µ = − (x)∂µ j µ .

But if the equations of motion are obeyed, then the action is invariant under any variation, including thisone, for arbitrary (x). But this means that ∂µ j µ = 0, the current is conserved.

90

The associated current whose charge density is π1 ∂x φ (which integrates over space to the
topological charge) is
1 µν
1(∂x φ, −∂0 φ)µ = ∂ν φ.
j̃µ =
2π2πThis is conserved because of the equality of the mixed partials: µν ∂µ ∂ν = 0.
Let’s expand in normal modes: φ = φL + φR withrφL (t + x) = qL + (p + w)(t + x) − i| {z }1≡ 2T
pLr
φR (t − x) = qR + (p − w)(t − x) − i| {z }
1≡ 2T
pR

L X ρn in(t+x) 2πL ,e
4πT n6=0 n
L X ρ̃n in(t−x) 2πL ,e
4πT n6=0 n(73)

The factor of n1 is a convention whose origin you will appreciate below, as are the othernormalization factors. Real φ means ρ†n = ρ−n (If we didn’t put the i it would have been
−ρ−n ).
RL
Here q ≡ L1 0 dxφ(x, t) is the center-of-mass position of the string. The canonical momentum for φ is π(x, t) = T ∂0 φ(x, t) = T (∂+ φL + ∂− φR ).
QM. Now we’ll do quantum mechanics. Recall that a quantum mechanical particle on acircle has momentum quantized in units of integers over the period. Since φ is periodic, thewavefunction(al)s must be periodic in the center-of-mass coordinate q with period 2π, andthis means that the total (target-space) momentum must be an integer
Z L
Z L
(73)
Z 3 j = π0 ≡dxπ(x, t) = T
dx∂t φ = LT 2p00

So our conserved charges are quantized according to
p=

j
,
2LT

w

(73)(72) πm

=

L

,

j, m ∈ Z .

(Don’t confuse the target-space momentum j with the ‘worldsheet momentum’ n!)(Note that this theory is scale-free. We could use this freedom to choose units where
L = 2π.)Now I put the mode coefficients in boldface:r

1
φL (x ) = qL +pL x+ − i2T

L X ρn i 2π nx+
e L,
4πT n6=0 nr1
L X ρ̃n i 2π nx−
φR (x− ) = qR +pR x− − ie L,
2T
4πT n6=0 n
+

91

(74)

The nonzero canonical equal-time commutators are[φ(x), π(x0 )] = −iδ(x − x0 )which determines the commutators of the modes (this was the motivation for the weird
normalizations)[qL , pL ] = [qR , pR ] = i,[ρn , ρ†n0 ] = nδn,n0 , or [ρn , ρn0 ] = nδn+n0 ,and the same for the rightmovers with twiddles. This is one simple harmonic oscillator foreach n ≥ 1 (and each chirality); the funny normalization is conventional.
 2

Z

 1Z
π
2
H =
dx π(x)φ̇(x) − L =dx+ T (∂x φ)
2
T
∞
X

1
2
2
= L
pL + pR +π(ρ−n ρn + ρ̃−n ρ̃n ) + a}|4T {zn=12
π0+ T2 w2
2T

1
=
2L


 2
∞


X
j
2
+ T (2πm) + π
n Nn + Ñn + a
T
n=1

(75)

Here a is a (UV sensitive) constant which will not be important for us (it is very importantin string theory), which is the price we pay for writing the hamiltonian as a sum of normalordered terms – the modes with negative indices are to the right and they annihilate thevacuum:ρn |0i = 0, ρ̃n |0i = 0, for n > 0 .
Energy eigenstates can be labelled by a target-momentum j and a winding m. Notice that
there is an operator w whose eigenvalues are w. So when I write |0i above, I really shouldlabel a vacuum of the oscillator modes with p, w.
Nn ≡ n1 ρ−n ρn is the number operator; if we redefine an ≡[an , a†m ] = δnm and Nn = a†n an is the ordinary thing.

√ −1n ρn (n > 0), we have

Notice that (71) means that there are separately-conserved left-moving and right-moving
currents:(jL )µ = (jLz , jLz̄ )µ ≡ (j+ , 0)µ(jR )µ = (jRz , jRz̄ )µ ≡ (0, j− )µHere jL only depends on the modes ρn , and jR only depends on the modes ρ̃n :
r
2π
π X
+
+
ρn ei L nx
j+ = ∂+ φ = ∂+ φ(x ) = p + w +LT n6=0j− = ∂− φ = ∂− φ(x− ) = p − w +

92

r

2π
π X
−
ρ̃n ei L nx
LT n6=0

Here’s an Observation (T-duality): At large T(think of this as a large radius of the target space), themomentum modes are closely-spaced in energy, and exciting the winding modes is costly, since the string has atension, it costs energy-per-unit-length T to stretch it.
But the spectrum (75) is invariant under the operationm ↔ j,T ↔

1(2π)2 T

which takes the radius of the circle to its inverse andexchanges the momentum and winding modes. This is called T-duality. The required duality
map on the fields is
φL + φR ↔ φL − φR .
√(The variable R in the plot is R ≡ πT . )You probably won’t be surprised√ to learn that the value of T where the two spectra coincide– the self T-dual radius (R = 2) – is special. In fact there is an enhanced SU(2) × SU(2)symmetry at that point. We will understand this soon and there is a homework problem on
it. There are other special points, such as R = 2 where a winding mode w1 and a momentum
mode p2 become degenerate.
This is the free-fermion radius. The next collision (between w2√
and p3 ) occurs at R = 3 which is in fact a model with supersymmetry.
It will sometimes be useful to consider the euclidean-space formulation: it → τ . Then afunction of the lightcone coordinates (t ± x) becomes a function of the complex coordinatesw = (x + iτ ) , w̄ = (x − iτ )and the free scalar action isZ
S=T

d2 w ∂w φ∂¯w φ.

(Note that d2 w = 2dxdτ and ∂w = 12 (∂x − i∂τ ).) Note also that (following Polchinski’s
conventions) holomorphic = leftmoving.
Let’s think about who are the good operators in this model. The first thing we might try
is φ itself. Since the path integral is gaussian, we havehφ(z, z̄)φ(0, 0)i = −

1
z z̄
1
z
1
z̄
log 2 = −
log −
log
πT
a
πT
a πT
a

(76)

where a is a UV cutoff.
Three ways to get this answer (not really independent): (1) Use the equations of motion(i.e. the Ward identity). What I mean is this: the path integral is independent of choice of93integration variable. Therefore:
Z
0=

[dφ]

∂
(anything) .
∂φ(z, z̄)In particular,
Z
0=

[dφ]


∂
¯
e−S φ(0) = δ 2 (z, z̄) + 2T ∂ ∂hφ(z,z̄)φ(0)i .
∂φ(z, z̄)(77)

Two comments: (a) in the Hamiltonian formalism, the delta function comes about when thederivatives hit the time-ordering theta-function. (b) The eom is in fact an operator equation,¯
meaning that it (∂ ∂φ(z)= 0) is a true equation as long as there are no other operators
nearby. In (77) we have one other operator nearby, which produces one contact term. We’ll
come back to this point.
We conclude that2∂z ∂¯z hφ(z, z̄)φ(0, 0)i = − δ 2 (z, z̄) .
T
30This equation is solved by the math fact :∂ ∂¯ log |z|2 = 2πδ 2 (z, z̄) .

(78)(Clearly it’s true for z 6= 0. To see the delta function, integrate the both hand side.)(2) Fourier: use the momentum space propagator (beware my factors here):
Z
ip·x
2 e
hφ(z, z̄)φ(0, 0)i ∝ − d̄ p 2
pThis is just the Fourier space solution of the PDE (78).
(3) We can also arrive at this answer using operator algebra by decomposing φ = φ> + φ<into creation pieces and annihilation pieces.
Notice the virtue that this is a sum of holomorphic and antiholomorphic (actually meromorphic) pieces. Also notice that the cutoff dependence is just an additive constant. But
this is a weird propagator: it grows in magnitude as the points get farther apart. Notice alsothat we got a log because a massless scalar in 1+1d has scaling dimension zero. This meansthat φ is not such a great operator – its effects are not local.
The Hohenberg-Mermin-Wagner-Coleman Fact. This (logarithmic) growth of thecorrelations of a massless scalar field is what prevents true spontaneous breakdown of continuous symmetries in two dimensions and below (notice from the fourier representation thatthe growth is only faster in D < 2).
30

2

R In case you are worried about the factors of two: I am using Polchinski’s convention that d z = 2dxdt,
so d2 zδ 2 (z, z̄) = 1 means δ 2 (z, z̄) = 12 δ(x)δ(t).

94

If you look at Coleman’s paper linked above it says some nonsense about ‘there is no
massless scalar field theory in D = 2’; this is false. The thing there isn’t is Goldstonebosons, as his title correctly says. The field theory is fine, it’s just that the φ operator
itself is bad. And the fact that the fluctuations of φ (the supposed order parameter) are IR
divergent means that there is no long-range order.
However, there are other operators which are good. To begin, the derivatives of φ are good
operators. For example, the currents havehj+ (z)j+ (0)i =k.
z2

The power here says that j+ has scaling dimension 1. Actually, we can say somethingstronger: if we scale z and z̄ independently, it has scaling dimension (1, 0).
hj− (z̄)j− (0)i =kz̄ 2

has scaling dimension (0, 1).
(A current in a d-dimensional critical theory always has scaling dimension d: you have to
integrate j0 over space and get a pure number, the charge.)
Vertex operators. Since φ has logarithmic correlators, you might think that exponentiating it is a good idea. First let’s take advantage of the fact that the correlator (76) splits
into left and right to write φ(z, z̄) = φL (z) + φR (z̄):
hφL (z)φL (0)i = −

z
1
log ,
πT
a

hφR (z̄)φR (0)i = −

1
z̄
log ,
πT
a

hφL (z)φR (0)i = 0 .

(79)

Another set of good operators with definite scaling dimension arises from:
Vα,β (z, z̄) =: ei(αφL (z)+βφR (z̄)) : .
This is a composite operator which we have defined by normal-ordering. The normal ordering
prescription is: q, p, −, + , that is: positive-momentum modes (lowering operators) go on
the right, and p counts as a lowering operator, P
so in particular using the expansion (pleasebeware my factors here): φL (z) = qL + pL z + i n6=0 ρnn wn , we have
: eiαφL (z) :≡ eiαqL eiαpL z eiαPρn nn<0 n w

eiα

P

ρn n
n>0 n w

(I used the definition w ≡ e2πiz/L .)How should we think about this operator ? In the QM of a free particle, the operator eipx
inserts momentum p – it takes a momentum-space wavefunction ψ(p0 ) = hp0 |ψi and gives
hp0 |eipx |ψi = ψ(p0 + p).
95

It’s the same thing here, with one more twist.
In order for Vα,β to be well-defined under φ → φ + 2π, we’d better have α + β ∈ Z –momentum is quantized, just like for the particle (the center of mass is just a particle). Let’s
consider what the operator Vα,β does to a winding and momentum eigenstate |w, pi (withno oscillator excitations, ρn |p, wi = 0, n < 0):Vαβ (0)|w, pi = ei(α+β)q0 ei(α−β)φ̃0 eiαPn<0 ρneiαP

n>0 ρn|w, pi = eiαPn<0 ρn|w + α − β, p + α + βi
(80)
The monster in front here creates oscillator excitations; the important thing is that thewinding number has been incremented by α − β; this means that α − β must be an integer,
too. We conclude thatα + β ∈ Z, α − β ∈ Z(81)
so they can both be half-integer, or they can both be integers.
Using the gaussian result from homework one (or moving the annihilation operators to theright) their correlators areD0hVα,β (z, z̄)Vα0 ,β 0 (0, 0)i =
z

α2
πT

β2

.

(82)

z̄ πT

The zeromode prefactor D0 is:00

D0 = hei((α+α )qL +(β+β )qR ) i0 = δα+α0 δβ+β 0 .
This is charge conservation.
We conclude that the operator Vα,β has scaling dimension(hL , hR ) =1(α2 , β 2 ).
2πT

Notice the remarkable fact that the exponential of a dimension-zero operator manages tohave nonzero scaling dimension. This requires that the multiplicative prefactor depend on
the cutoff a to the appropriate power (and it is therefore nonuniversal). We could perform
a multiplicative renormalization of our operators V to remove this cutoff dependence fromthe correlators.
The values of α, β allowed by single-valuedness of φ and its wavefunctional are integers.
We see (at least) three special values of the parameter T :• The SU(2) radius: When 2πT = 1, the operators with (n, m) = 1 are marginal.
Also, the operators with (n, m) = (1, 0) and (n, m) = (0, 1) have the scaling behaviorof currents, and by holomorphicity are in fact conserved.
• The free fermion radius: when 2πT = 2, V1,0 has dimension ( 21 , 0). which is the
behavior of a left-moving free fermion.
96• The supersymmetric radius: when 2πT = 23 , V1,0 has dimension ( 32 , 0) and represents a supersymmetry current.
The expression for the correlator is easily generalized to multiple vertex operators:NYδPj αj ,0h Vαj (z)i = QN
αi ·αj2πT
j=1
i<j (zi − zj )Operator product expansion (OPE).
In QFT, if we act with a bunch of local operators close together (compared to their separation from other operators), their net effect is that of another local operator. If we have abasis of operators in our CFT, we may expand this net effect in terms of this basis. So far
this is true in any QFT. In a scale-invariant theory, we can choose a basis of operators ofdefinite scaling dimension31 ,Oa → λ−∆a Oa ,hOa† (z)Oa (0)i ∼

1
z 2∆a

.

We can organize this expansion (the operator product expansion) in order of increasing dimension, {Oa , ∆(a) ≤ ∆(b) if a < b}. Then the operators with the longest range of influence
come first:
OA (z)OB (0)z→0∼CAB1 (z)O1 (0) + CAB2 (z)O2 (0) + · · ·

=

cAB1
cAB2
O1 (0) + ∆ +∆ −∆2 O2 (0) + · · ·
∆
+∆
−∆
1
z A B
z A B

(83)

In general, there will be power-law singularities in the coefficients, produced by the collisionsof the operators; the powers are determined by scaling. But at some point in the expansion,
the effects start to vanish as z → 0, i.e. ∆A + ∆B − ∆n < 0 for ∆n big enough. The terms
in the OPE after this are regular terms and can often be ignored, for example because they
don’t effect correlation functions without more insertions nearby. In the future, we will be
writing expressions like1
∂φ(z)∂φ(0) ∼ 2 + regular terms .
zAfter this, “plus regular terms” will be implicit.
We saw earlier that free massless fermions in 1+1d have groundstate correlators that golike1hΨ† (x)Ψ(0)i ∝ − .
(84)
x
31

Some fine print: actually, this requires that the dilatation operator is hermitian. This fails to be the casein some non-unitary theories, for example in logarithmic CFTs where the dilatation operator involves Jordan
blocks, and where therefore the OPE involves also some logarithms (which shift under scale transformations).

97

1
In terms of the OPE, this can be expressed as ψR† (x, t)ψR (0) ∼ − t−x
or in euclidean coords1ψR† (z̄)ψR (0) ∼ .
z̄
The operators eiφ(x) satisfy
2nm: einφR (z) :: eimφR (0) :=: ei(n+m)φR (0) : z 2πT + ...

(85)

Note that the power of z agrees with our general OPE formula since n2 + m2 − (n + m)2 = −2nm.
When we set n = −m and πT = 1, this looks like the fermion expression (84).
In general the ... in (85) can contain further singular terms coming from the Taylor expansion of φL (z) = φL (0) + z∂φL (0) + · · · . That is, going more slowly, we have
2nm

: einφR (z) :: eimφR (0) : = : einφR (z)+imφR (0) : z 2πT2nm
2nm= : ei(n+m)φR (0) : z 2πT + in∂φR (0) : ei(n+m)φR (0) : z 2πT +1 + ... (86)Note that the bosonization identities which we will discuss further below are statements
about the singular parts of the OPE of these operators – the finite, less-universal pieces are
not part of the conversation.
[End of Lecture 14]983.1.1

Application 1: Bosonization dictionary[Fradkin 2d ed, page 133; Sachdev 2d ed, page 420] Let’s set T = π1 for a bit. The left andright U(1) currents in terms of bosons and fermions areΨ†L ΨL =

1
∂z φ,
2π

Ψ†R ΨR =

1 ¯
∂z φ.
2π

The fermions themselves are:
ΨL (z)† =: eiφL (z) :, ΨR (z̄)† =: eiφR (z̄) : .
How might a human have figured this out? (One human who did was Mandelstam.)
The relation for the charges says that the total fermion charge is the winding number:
1
∂x φ.
2π

(87)

1
∂x φL = jL .
2π

(88)

Ψ†L ΨL + Ψ†R ΨR =Let’s consider the refined statement:Ψ†L ΨL =This means that if we create a left-moving fermion at x = x0 , the value of φL (x) should jumpfor all x < x0 , by a fixed amount so that if we integrate the BHS of (87) we get 1. (Thatamount is 2π.) For example, according to (88) a state with no fermions looks like:+|no fermionsi =Acting on this state with one fermion at x = x0 should do the following:+Ψ†L (x0 )+=.

Such a jump is produced by the conjugate field momentum (recall the translation operator
in QM is e−ipa ψ(x) = ψ(x + a)):R x0

ΨL (x0 ) ∼ ei2π −∞ dxπL (x) .
99

Now who is the canonical momentum πL (x) in terms of φL , φR ? It’s the thing with acanonical commutator with φL . One way to figure this out is to go back to the mode
expansions and show that2i(89)[∂x φL (x), φL (y)] = δ(x − y)Twhich shows that πL (x) = T2 ∂x φL (x) (and similarly πR = − T2 ∂x φR (x)). (Notice that thisequation is roughly the imaginary part of the x derivative of (79).) So (when T = π1 , which
is when this operator is single-valued!) the fermion operator is
R x0

ΨL (x0 ) ∼ ei −∞ dx ∂x φL = eiφL .
We already saw that this kind of object creates winding in (80). So again the fermion is adomain-wall creation operator; notice that because the field φ is periodic, it can be possible
to have a single such domain wall on a circle. We’ll say a bit more about the periodicity of
φL below.
In order for this to be a fermion operator, we require {ΨL (x), ΨL (x0 )} = 0, x 6= x0 . Is it
true? Yes in fact, as we can see by BCH, in the form eA eB = eB eA e−[A,B] if [A, B] is ac-number. A useful intermediate fact comes from the x-integral of (89) (from −∞ to x)
which says[φL (x), φL (y)] = iπsign(x − y)(the additive constant of integration on the RHS is fixed by demanding antisymmetry inx → y which is manifest on the LHS). Therefore:: eiφL (x) :: eiφL (0) := e−iπsign(x) : eiφL (0) :: eiφL (x) := − : eiφL (0) :: eiφL (x) : .
More generally, the same calculation shows that
: einφL (x) :: eimφL (0) := (−1)nm : eimφL (0) :: einφL (x) :So einφL is bosonic for even n and fermionic for odd n.
You can also see the fermi statistics from the OPEs:
1
1=−.
: eiφL (z1 ) :: e−iφL (z2 ) :∼z1 − z2z2 − z1: eiφL (z1 ) :: eiφL (z2 ) :∼ (z1 − z2 ) : eiφL (z1 )+iφL (z2 ) :Notice that this object

1

R x0

R x0

ei π −∞ dx∂x φL = ei −∞ dxjL (x)has a lot in common with the Jordan-Wigner string, which counts the number of particles
to the left.
Periodicity of chiral bosons. Since φ ' φ + 2πm, m ∈ Z and φ = φL + φR it is tempting
to suggest that φL has period π. Such a period would mean that e−iφL ∼ ΨL is not singlevalued. But it is not quite correct. To see the correct statement, notice that the T-dual
100field is φ̃ = φL − φR . With our normalization, this field also has periodicity 2π. The winding
modes of the dual fields are momentum modes of the original one, so   
φL + φR ' φL + φR + 2πm
φL
φL
m+j
=⇒
'
+π
,
φL − φR ' φL + φR + 2πj
φR
φR
m−j
where m and j are independent integers. But then m ± j are not independent integers: whenm + j is odd or even, so is m − j. This is the same fact as my earlier statement about thelabels α, β on Vα,β at Eq. (81).
Note that the T-dual field is often called ϑ ≡ φL − φR in discussions of Luttinger liquids(e.g. in Fradkin’s and Sachdev’s books).
Possibly Depressing comment. So now you are starting to see that this duality businessis actually often a sad story: we thought we could solve two systems (free bosons and freefermions) but since they are really the same system in disguise, it turns out we can only
solve one!
Multiple fermions. Notice that if we had N complex fermions ψ a (for example a could
be a spin index), we could bosonize them by N scalar fields at the free fermion radius:?

aψ a (z) ∼ eiφ (z) , a = 1..N .

(90)

This formula almost works, but it has the problem that the operators on the RHS with different values of a will still commute (since the associated scalars don’t have any singularities intheir OPE, i.e. they commute). This must be fixed by so-called Klein factors (or 2-cocycles)
which introduce the extra necessary signs. The necessary object is constructed by ordering
the species of fermions a = 1..N and defining:ca ≡ (−1)

Pb<a Nbwhere Nb is the total fermion number of species b. The correct version of (90)a

ψ a (z) = ca eiφ (z) , a = 1..N(no sum on a).

(91)

Spin fields in (two copies of ) the Ising model.
Now I will partially fulfil my promise to use bosonization todetermine the scaling dimension of the spin operator at the
Ising critical point. Recall that I claimed without any justification that in the TFIM, the behavior of the magnetizationhZi for g <∼ gc washZi ∼ (gc − g)1/8(as in the figure). This 1/8 is the scaling dimension of the order parameter operator in thecritical Ising CFT, the operator onto which the spin Zj matches.
101

To begin, let’s consider just the holomorphic part of a complex fermion ψ(z) = √12 (χ1 (z) + iχ2 (z)).
This is not exactly the Ising critical theory, which is instead a non-chiral majorana mode,χL (z), χR (z̄), but this is what we know how to bosonize. The bosonization map isψ(z) ∼ eiφ(z) .
Recall that the fermions in our spin chain arose as domain wall operators. They satisfy
Zl cj = −cj Zl ,l > j.

We can interpret this formula to say that the spin Zj creates a branch cut for the fermion
field. The CFT version of this statement is:1

ψ(z)σ(0) ∼ z − 2 µ(0) .

(92)Here σ, µ are spin fields (or more generally, twist fields) whose presence at z = 0 creates abranch cut. The actual location of the branch cut is not physical, but the location of the
branch points (here 0, ∞) is meaningful. Can we find operators which accomplish this?
Yes: they areiiσ(w) = e− 2 φ(w) , µ(w) = e+ 2 φ(w) .

(93)

To check this claim, their OPEs with the fermion operator are:
i

i

1

i

i

1

ψ(z)σ(0) =: eiφ(z) :: e− 2 φ(0) :∼ e+ 2 φ(0) z − 2 + regularψ(z)µ(0) =: eiφ(z) :: e+ 2 φ(0) :∼ e+3 2 φ(0) z + 2 + regular(94)

As a check on my ability to do wick contractions, consider the OPE between (and hence thetwo-point function of) two of the spin fields, which is a special case of (82):eαiφ(z) e−αiφ(0) ∼For α = ± 21 , this gives

3

The dimension of ei 2 φ is 12

1
.
z α2

1
1
∆σ = ∆µ = α2 |α= 1 = .
2
2
8

3 2
= 89 so the dimensional analysis in (94) works.
2

The answer ∆σ = 18 is the dimension of the spin field in the theory of a complex chiralfermion. This is two copies of the left-moving part of the Ising CFT. The spin field in the
Ising model which makes a branch cut for both χL and χR has the same dimension (the one
1that makes a branch cut in just χL has dimension 16). The basic fact I am using here is thatthe spin-field dimension is additive in the number of fermion fields which feel the branch cut,that is, in the number of fermion fields whose boundary conditions are twisted.

102

To see this, let’s make spin fields for N complex fermions ψ a , each bosonized as we justdid by N scalar fieldsaψ a (z) ∼ eiφ (z) , a = 1..N
The operators
1a a
σs (z) ≡ eis φ (z) , {sa = ± }2create branch cuts for all N fermions at once:aψ a (z)σs (0) ∼ z s σs0 (0) + regularwhere (s0 )a = sa + 1, (s0 )b6=a = sb . The dimension of these operators is the sum of thedimensions:
 N
1 2
∆s =
s1 + s22 + · · · + s2N = .
2
8
Now take N = 12 to get the answer for a single majorana fermion. If you are not satisfied
by this argument, there is a more direct way to compute the dimension of the Ising model
spin field, using the defining equation (92) which requires a bit more CFT technology (it’sexplained on p. 74 of Ginsparg’s CFT notes). We’ll come to that soon.
Derivation of spin field.
Notice that the object (93) would not be considered single-valued on the original boson
hilbert space! However, the operator which creates a branch cut for both ψL and ψR at theisame time is e 2 (φL (z)±φR (z̄)) and this is an allowed operator. And in fact, here is a derivationof its relation to the spin operator. Recall the JW formulaσj+ = (−1)

Pi<j nic†j .

(Previously we put the JW string to the right; here I defer to the bosonization convention
we used above.) Here (in the basis we used earlier, where Z is the Ising order parameter),σ + ≡ Z + iY, so σ + + σ − = 2Z. In the continuum, in te rms of the fermions, this is
Rx

0

0

σ + (x) = eiπ −∞ dx j0 (x ) Ψ(x)† .
Now, we must use (65) to expand our fermion operator near the fermi surface. Let’s do justthe R piece first (then we’ll put it all together):
Rx

0

0

eiπ −∞ dx j0 (x ) ΨR (x)†
(87)

=
=
∼

∂y φL +∂y φR

Rx

−∞ dy
2π
e−iφR (x)e−iπφ (x)+φ (x)i L 2 R−iφR (x)e φ (x)−φ (x) e
L
R
2
ei

OK, that was practice; the full expression is:
σ + (x)

=

Rx

0

0

eiπ −∞ dx j0 (x ) Ψ(x)†
103

(95)

(87),(65)Rx=e−iπ −∞ dy∼ei∂y φL +∂y φR
2π

φL (x)−φR (x)2e−ikF x e−iφR (x) + e+ikF x e−iφL (x)e−ikF x + e−iφL (x)−φR (x)2e+ikF x
(96)[End of Lecture 15]
The role of the anomaly in bosonization. Recall that I promised that the chiralanomaly would play a role in bosonization. Here’s one way
Pin which it does. Consider thecanonical commutator of the bose mode operators (∂z φ ∝ n ρn z −n−1 ):[ρn , ρ−n0 ] = δn,n0
in terms of the fermions, ρn =P†
l Ll Ll+n , we have:[ρn , ρ−n0 ] =X[L†l Ll+n , L†l0 Ll0 −n0 ]

l,0 l


X †
?
Ll0 −n Ll0 −n0 − L†l0 Ll0 +n−n0 = 0
=

(97)

l0

This expression has the classic form of an anomaly: it looks like this vanishes since we canredefine the dummy index in the second term by n0 → n0 + n and make the two terms look
the same with opposite sign. However, in the presence of a UV cutoff, i.e. in the regulated
theory, that shift will shift the cutoff and will not be innocuous since the summand is notzero for arbitrarily high energies.32 If, for example, we regulate the short distance behavior
using normal-ordered operators, we have
X †[: Ll Ll+n :, : L†l0 Ll0 −n0 :][ρn , ρ−n0 ] =
l,0 l


X †
†
= δnn0hLl−n Ll−n i − hLn Ln i = δnn0 n.

(98)

l

This calculation is related to the anomaly because it says that the OPE of the chiral current
jR with itself has a singularity with a definite coefficient. Specifically (as you can check on
HW 6),kj(z)j(0) ∼ 2zwith k = 1 is called the level of the (here U(1)) current algebra. This singularity implies
the anomalous nonconservation when coupling to a background gauge field. In terms of thefermions, both of these effects come from the ‘diangle’ diagram:32

To see how the high-energy stuff comes in consider the analogous infinite-volume expression
Z ∞
Z ∞
Z
?
dp(f (p) − f (p)) =dp(f (p + s) − f (p)) = dp(sf 0 (p) + O(s2 )) = s(f (∞) − f (−∞)).
0=
−∞

−∞

104

Twist fields and boundary conditions. Above I used the phrase “fermion fields whoseboundary conditions are twisted” by the spin fields. With this phrase I betrayed that I am
secretly using the following picture, which is called radial quantization (more on this at §4.2).
Imagine that the radial coordinate of the plane is the time coordinate (it is euclidean so thisinterpretation is not problematic). Then equal-time surfaces are circles centered at the origin.
(We’ll see below that in a CFT this is equivalent to the ordinary notion of time evolution by
a conformal transformation.) In this case, the boundary conditions on the spatial coordinate
are determined by going around the origin; if there’s a branch point at the origin, we mustcross the branch cut in doing so, and the field comes back to itself up to the discontinuity
across the branch cut, which here is a minus sign. So the spin fields can be regarded as“boundary-condition-changing operators”: they take the fermions from
P the sector with PBC(and hence a fourier expansion with integer powers of z: ψP BC (z) = n∈Z z n ψn ; this is calledP
1
the Ramond sector) to APBC (half-integer mode expansion ψAP BC (z) = n∈Z z n+ 2 ψn+ 1 ;2this is called the Neveu-Schwarz sector) and vice versa.
Understanding in terms of group theory. I must comment on the very essential grouptheory underlying the previous discussion of spin fields for multiple fermions. The theory
of N complex free chiral fermions has a manifest U(N ) symmetry rotating the complex
fermions into each other: ψ a → Uab ψ b , UU† = 1. The free fermion action actually enjoys alarger symmetry which ignores the grouping into complex fermions: ψ a = √12 (χ2a−1 + iχ2a );R¯ B , these 2N chiral majorana fermions transform in the fundamental ofsince S ∝ δAB χA ∂χT
Ban SO(2N ) symmetry, χA → OA
B χ , O O = 1.
This SO(2N ) acts on the hilbert space. In terms of the fermions, the generators are
just fermion bilinears j AB = χA χB , the currents (recall that the adjoint (the rep that thegenerators transform in ) of SO(2N ) is made from antisymmetrizing two fundamentals).
There are 2N (2N − 1)/2 of these.
In terms of the bosons, the generators come in two types (in the theory of Lie groups
this is called Cartan-Weyl basis). N of them are just the ordinary boson currents ∂φa .
This is called the Cartan subalgebra (the number of them is N = rank(SO(2N )) – they allcommute with each other. The non-commuting (charged) generators are vertex operatorsaeiφ wa where
PN in order for this to be a conserved current, it must have dimension 1, which
requiresa=1 wa wa = 2. (in Lie group terminology, wa is a root vector). Recall thatthe entries of w must be integers. This condition is solved by vectors of the form {wa } ={(0...0, ±1, 0 · · · , 0, ±1, ...)}, of which there are N (N −1)/2×4. Altogether we have recovered
the same number of generators N + 2N (N − 1) = N (2N − 1), X.
The states in the Neveu-Schwarz (APBC) sector are of the formBA
ψ−
1 ψ 3 · · · |0i,
−
2

A
ψr>0
|0i = 0, ∀r

2

and therefore transform in various antisymmetric tensor representations of this SO(2N ).
1

1

1

What is going on with these spin field objects eis φ = ei(± 2 φ ± 2 φ ···± 2 φ ) ? These operators
a a

105

1

2

N

map the Neveu-Schwarz sector to the Ramond sector. For a hint about what they are doing,
notice that the fermion zeromodes in the Ramond (PBC) sector satisfyBAB{ψ0a , ψ0b } = 2δ ab =⇒ {χA
0 , χ0 } = δwhich is the algebra of Dirac gamma matrices for SO(2N ), i.e. Clifford algebra. This algebra
must be represented on the Ramond groundstates, which are therefore not unique, i.e. they
are degenerate. The point in life of Clifford algebra is to construct spinor representations. So
the answer is that these vectors {sa } = {± 12 , ± 12 , · · · , ± 12 } are weight vectors for the spinor
representation(s) of SO(2N ), in which the Ramond groundstates transform. The spin fields
create the Ramond groundstates from the NS groundstate:
σs (0)|0iN S = |si.
A very interesting application of this group theory arises in the demonstration (by Fidkowski and Kitaev) that eight free majorana chains in the topological phase (with the dangling majoranas at the ends) can be adiabatically connected to the trivial phase (only) via
interacting Hamiltonians.

3.1.2

Application 2: Briefly, what is a Luttinger liquid?

[Fradkin p. 152] Let us consider a model which is actually made of fermions (which may
have spin) at finite chemical potential. The actual density of the microscopic fermions Ψ is
XXρactual (x) =Ψ†σ (x)Ψσ (x) = ρ0 + ρL +ρR +e2kF ix ψR† ψL + h.c. ...
|{z}
|{z}
σ

k

=Nspin πF

†
=ψL
ψL

σ

where I am reminding you that Ψ(x) = eikF x ψR (x) + e−ikF x ψL (x). The .... are even-higher
momentum contributions which we will ignore. The quantity 2kF is important in any systemwith a fermi surface as it is the diameter of the fermi surface, and (bosonic) particle-holeexcitations near (just below) this momentum can be cheap despite its large value.
Consider the hamiltonian H = H0 + Hpot where the free Galilean-invariant hamiltonian is


XZ
†
†
H0 =
d̄p pvF ψσR ψσR − ψσL ψσL .
σ

Subject this system to an (at-the-moment) external potential
Z
Z
Z
Hpot = dx eV (x)ρ(x) = −e d̄p d̄q Ṽ (q)Ψ† (p + q)Ψ(p) .
Let us consider the special wavenumber V (x) = V0 cos(2kF x), soZ
Z†
Hpot = −eV0 dpψR (p)ψL (p) + h.c. = −eV0 dxψR† (x)ψL (x) + h.c.
106

In the relativistic language, this is a mass term:
Z
Z


†
3
1
H = dxψ σ vF i∂x + eV0 σ ψ = dx ψ̄vF γ x ∂x ψ + eV0 ψ̄ψ(here I used γ 0 = σ 1 , γ x = σ 1 σ 3 , ψ̄ ≡ ψ † γ 0 ).
Where does this potential come from? It can be generated by (e.g. Coulomb) interactions
between the fermions themselvesZZ ZeV KV .
Hint = dxeV ρ +2
If we treat this interaction in mean field theory, we get back (for some choice of K more on
which more below) the previous background potential, Hint ∼ Hpot . On the other hand, if
we integrate out V , we find a density-density interaction
Z Z
Hint = eρK −1 ρwhich for the right kernel, K, might be the Coulomb interaction. In a metal, the coulombinteraction is screened and therefore short-ranged. Expanding in this range, we find


Z





g4 ρ2R + ρ2L+gu lim ψR† (x)ψR† (y)ψL (x)ψL (y) + h.c. Hint = dx 2g2 ρR (x)ρL (x) +
y→x
{z
}
|
{z
}
|
back-scattering, 1d

forward scattering, ubiquitous

The gs can be written in terms of various (fourier) components of the interaction potential
K −1 . In the last term, you should think of y = x + a, one lattice spacing. The last term
describes a process whereby two L/R fermions turn into two R/L fermions, which clearly
violates momentum conservation by ±2kF ; this can happen in a lattice model if the twosides of the fermi surface differ by a reciprocal lattice vector:
kF ≡ −kF +2π.
aThis is called umklapp scattering. When does it happen? kF is determined by the densityof electrons. The maximum density is when all states are filled and is Ns 2π, and then thereaπ
is no fermi surface. Umklapp happens at half-filling, when kF = a . Otherwise gu can beignored. Notice that this effect explicitly breaks the chiral symmetry ψL/R → eiαL/R ψL/R
which independently rotates the phases of the left-movers and right-movers down to a Z2
symmetry, which is just fermion number (which is hard to violate).
Bosonize. Now
R we bosonize. Using φ = φL + φR , Θ = φL − φR and our dictionary above,
we find H0 ≡ dxh0 (x), withh0 (x) = vF
vF
Π2 + (∂x φ)2 =(∂x Θ)2 + (∂x φ)22
2
107

where in the second step we used the canonical commutators we found previously. Using
1∂x φL etc this isρL = 2π
h0 (x) = πvF (ρL + ρR )2 + (ρL − ρR )2 = 2πvF ρ2L + ρ2R .
The forward scattering hamiltonian density is
hforward (x) = g4 ρ2R + ρ2Lof exactly the same form as the free term (!) and therefore just renormalizes the velocity,
the coefficient of the spatial gradient term. The back-scattering term is less innocuousg2(∂x φ)2 − Π2hback (x) = 2g2 ρR ρL =#Umklapp is a little trickier: to define the bosonized form we must absorb the OPE singularity
1: e−iφR (x) :: e−iφR (y) :∼ (x − y) 2πT : e−2iφR (x) : +...
into the coupling constant:hu = g̃u cos 2φ.
1− 2πT(It is useful to set g̃u = gu awhere a is a short-distance cutoff.) With g̃u 6= 0, theremaining chiral symmetry acts by φ → φ + π.
The full hamiltonian density (without umklapp, so away from half-filling) is


v 1
2
2
(∂x Θ) + K (∂x φ)h=
2 K
where v and K can be written in terms of the gs and vF 33 . Notice that T-duality is clear inthis expression as T : K ↔ K1 , φ ↔ Θ.
A few comments:• The backscattering term g2 affects both the velocity and the ‘stiffness’ K (this is
the parameter we called the tension T earlier). The signs are such that a repulsive
interaction g2 > 0 increases K.
• In the absence of umklapp, the chiral symmetry acts on the boson as φ → φ + α,1α ∈ [0, 2π). Think about how this acts on the fermions: since ρ = ... + 2π
∂φ, thedensity transforms as




α
i2(kF x−α) †ρ(x) → ρ0 + ρR + ρL + e
ψR (x)ψL (x) + h.c. = ρ x −kF– this is a rigid translation of the density profile! It makes sense then that latticeeffects would break this symmetry.
33

The relation is
q
v=

s
2

2

(vF + g4 /π) − (g2 /π) , K =108vF + g4 /π + g2 /π.
vF + g4 /π − g2 /π• We have ignored the fact that the dispersion will not be exactly linear arbitrarily far
from the Fermi surface. Such terms can also be bosonized and become things like(∂x φ)3 which change the story dramatically; in particular they will lead to damping ofthe ‘sound mode’ described by φ even away from half-filling.
• If we study spinful fermions, there are even more 4-fermion terms to consider. A payoffof studying the case with spin is that it provides an example of spin-charge separation:
by bosonizing like φc ≡ φ↓ + φ↑ , φs ≡ φ↓ − φ↑ we can see that the charge sector (φc )and the spin sector (φs ) decouple and can have independent vs and Ks – the spin
excitations of the electron travel at a different speed than its charge! There is a lot
more to say about this and I recommend Fradkin’s treatment in chapter 6.
[End of Lecture 16]109How is a Luttinger liquid different from a Fermi liquid?
Continuing with the model with actual microscopic fermions, we have encountered a gaplessphase in D = 1 + 1 described by free fermions, which we also know how to build in D >1 + 1. The excitations are characterized by long-lived fermionic quasiparticles with the sameU(1) charge as the constituent fermions cj – they are Landau quasiparticles. This has theconsequence that there are sharp peaks in the fermion green’s functions,ZZ
eiωt−ikja hc†j (t)c0 (0)i = G(ω, k) ∼
ω − vF (k − kF )(the RHS is the contribution from the right-mover ΨR (t − x)). There is (therefore) a jumpin the fermion occupation number at k = kF :Z ∞dω†ImG(ω, k) = Zθ(k − kF ).
n(k) ≡ c (k)c(k) =| {z }−∞ π∼Zπδ(ω−vF (k−kF ))[Fradkin §6.8] These last two characteristic features of Fermi liquid theory can be destroyed
by interactions in 1 + 1 dimensions, as we can see using the boson description. As a simple
first step, notice that as we change T away from π1 , the dimension of the operator which at
T = π1 creates the right-moving fermion eiφR changes away from 21 :heiφR (x,t) e−iφR (0,0) i ∼11|t − x| πT.

It is tempting to conclude that the fermion correlation function Ψ(x) = eikF x ΨR (x) + ... is?

hΨ(x, t)† Ψ(0, 0)i ∼sin (kF x)
1|t − x| πT

The fourier transform of this has a broad power-law peak (a branch point) instead of a pole;
the spectral density has no delta function, but rather goes like1?

ImG(ω, k) ∼ |ω − vF (k − kF )| πT +1 θ(ω − vF (k − kF )),and hence there is no jump at k = kF .
This is not quite the correct story however, since the operator eiφR is not quite the operator
which creates an electron, except at T = π1 (for example, it is not a fermionic object!).
Instead what happens is that the interactions mix the R and L operators; the operator
which creates the rightmoving electron is actually the linear(-in-the-exponent) combination:
ΨR = eiα+ φR +iα− φL ,with α± =110(K ± 1)2
4K


1(recall that K = πT); it has total dimension ∆ = 41 K + K1 . The statement about the polebeing replaced by a branch point is correct.
(Pieces of) the fermion Green’s function are (relatively) directly measured in angle-resolved
photoemission spectroscopy (ARPES) experiments. Other combinations of them can bemeasured by scanning-tunneling microscopy (STM).
[Fradkin, 2d ed. p. 170] There is a lot more to say about this subject. In particular, one
can understand in great detail various ordering instabilities of interacting fermions in one
dimension using the technology we have developed.

3.1.3

Application 3: sine-Gordon and Thirring

Let us reconsider the scale invariant theory described by the free scalar again.
L=T(∂µ φ)2 + gn : cos nφ :2

I allow the possibility of a cosine term. For the case of the Luttinger liquid we found thatumklapp produced n = 2. There are several important points to make here:(-1) In the absence of lattice effects or a background gauge field, there is a U(1) × U(1)symmetry acting on φL and φR independently.
(0) The theory is scale invariant for any T : there is a line of fixed points parametrized by
T . That is, T is the coefficient of an exactly marginal operator. This is weird and the only
examples in D > 1 + 1 (that I know) require supersymmetry.
(1) by changing T we can study interacting fermion systems.
(2) the fermion mass corresponds to a cos φ term.
(3) These cos nφ terms have a dimension that varies with T . This means that as we vary
T they can go from irrelevant to relevant. If they are relevant they produce a mass foreverybody (You can see this by Taylor expanding about the minimum of the potential).
This happens when
n2
n2
=K
≤2
∆n =
2πT
2
2
n2
that is, for T > Tc = 4π
(or K < Kc = n2 )Brief comments on phase diagram. [Sachdev Ch. 20 – beware that φSachdev = 2φhere ]The RG equations (for n = 1) aredv= (2 − 4K)v,
d`
111

dK
= −δv 2d`where δ > 0 is some non-universal quantity. This leads to the following RG flow diagram in
the space of these couplings (from Sachdev chapter 20):

The figure at right is the phase diagram we can then infer for the modelXX
H = −J1(Xi Xj + Yi Yj + λZi Zj ) − J2
(Xi Xj + Yi Yj + Zi Zj ) .
hiji

hhijiiwhere hhijii denotes next-nearest neighbors. The simpler model we’ve been discussing moves
along the λ axis.
• When v → −∞, we can find the groundstates bymaking the −|v| cos 2φ happy first. The minima occur
at φ = (n + 12 )π, n ∈ Z. They can be distinguished by
the value of sin φi = (−1)n = h(−1)i~σi · ~σi+1 i. Thisis described by singlets on alternating bonds. Thisis called a valence bond solid state. The low energy
excitations in this state are interesting – they are notjust the gapped φ excitations rolling in the minima of the cosine potential. Rather, the1domain walls between the two groundstates (created in the CFT by ei 2 φ ! ) have lower energy.
We can see this explicitly at the special point K = 14 (which is < Kc ) by ‘refermionizing’. The
1previous original fermions where created by ψL† = eiφR = ei 2 (φ+Θ) . Introduce new fermions
1

Θ

1

Θ

ψ̃L† ≡ ei 2 (2φ+ 2 ) , ψ̃R† ≡ ei 2 (2φ− 2 ) .
At K = 14 these are free fermions, and the v cos 2φ term is just a mass term v ψ̃L† ψ̃R + h.c..
Interestingly, there is a regime where the quantum of the boson can decay into a pair of these(fermionic) domain walls.
• The phase at v → +∞ is similar. Minimizing +|v| cos 2φ gives φ = nπ which are distinguished by cos φ = (−1)n = h(−1)i Zi i. This is ‘Ising Neel’ order, that is: staggered Ising112

spins. (Neel order refers to staggered magnetization in an arbitrary direction, spontaneously
breaking SU(2) spin rotation symmetry.)
The two phases at large |v| are related by the relabelling v → −v, φ → φ + π/2. In termsof the bosons or fermions both of these large |v| states are charge density waves, where thedensity of particles breaks the translation symmetry down to a smaller discrete subgroup.
Most continuous phase transitions occur by tuning the coefficient of a relevant operator
to zero (recall the TFIM example where we had to tune g → gc to get massless fermions).
This is not what happens in the XY model; rather, we are varying a marginal parameterand the dimensions of other operators depend on it and become relevant at some criticalvalue of thatmarginal parameter. This leads to very weird scaling near the transition, of thea− √K−Kc(for example, in the correlation length, the exponential arises from inverting
form e1expressions involving GR (z) = − 4πK
log z) – it is sometimes called an ‘infinite order’ phasetransition, because all derivatives of such a function are continuous. I don’t want to say
more about how this comes about because I know Congjun Wu explained it in the RG class.
In the high-energy literature, this is called duality between the Thirring model and thesine-Gordon model.

3.1.4

XY transition from SF to Mott insulator[This discussion is from Ashvin Vishwanath’s lecture notes.] So far in this chapter, we’ve
applied bosonization to understand systems made of spins, fermions or bosons. It will be
useful to make the application to bosons more direct and explicit. Let’s return for a momentto the Bose-Hubbard model (in any dimension)
HBH = −J˜

 UX
X
X †
ni (ni − 1) − µnibi bj + h.c. +
2 i
i
hiji

with ni ≡ b†i bi , and [bi , b†j ] = δij . This Hilbert space has a useful number-phase representation in terms of[ni , φj ] = −iδij , φi ≡ φi + 2π, ni ∈ Z(where the last statement pertains to the eigenvalues of the operator). The bosons are√√bi = e−iφi ni , b†i = ni e+iφi ;these expressions have the same algebra as the original bs. In terms of these operators, thehamiltonian is
X √
X
 UX
√
HBH = −J˜ni ei(φi −φj ) nj + h.c. +ni (ni − 1) − µ
ni .
2 i
i
hiji

113

√
√
If hni i = n0  1, so that ni = n0 + ∆ni , ∆ni  n0 then bi = e−iφ ni ' e−iφi n0 andXUX˜ 0cos (φi − φj ) +(∆ni )2 ≡ Hrotors
HBH ' − 2Jn| {z }
2
i

hiji

≡J

where we set n0 ≡ µ/U  1. This is a rotor model.
This model has two phases:
U  J : then we must satisfy the U term first and the number is locked, ∆n = 0 in the
groundstate. This is a Mott insulator, with a gap of order U . Since n and φ are conjugate
variables, definite number means wildly fluctuating phase.
U  J : then we must satisfy the J term first and the phase is locked, φ = 0 in the
groundstate, or at least it will try. This is the superfluid (SF). That is, we can try to expand
the cosine potential34

X
X
X
X
1
2
2
2
Hrotors = U
ni − Jcos (φi − φj ) ' U
ni − J
1 − (φi − φj ) + ...
2
i
i
hiji

hiji

which is a bunch of harmonic oscillators and can be solved by Fourier: φi = √1 d
N
so
X
H'
(U πk π−k + J (1 − cos ka) φk φ−k )−ik·xi
φk ,
ie

P

k

This has gapless phonon modes at k = 0. I have written the hamiltonian in 1d notation
but nothing has required it so far. The low energy physics is described by the continuum
lagrangian density!
 2ρs (∂τ φ)2~+ c ∇φLeff =
2cp√with ρs = J/U , c = JU . ρs is called the superfluid stiffness. This is the free scalar
theory above, with T = ρs .
Now 1d comes in: In d > 1, there is long range order – the bosons condense and spontaneously break the phase rotation symmetry φ → φ + α; the variable φ is a Goldstone boson.
In 1d there is no long-range order. The two phases are still distinct however, since one has
a gap and the other does not. The correlators of the boson operator bi ∼ eiφi diagnose thedifference. In the Mott phase they have exponential decay. In the “SF” they haveheiφ(x) e−iφ(y) i =

c0
,
rη

η=

1
1
=
.
2πρs
2K

This is algebraic long range order.
However: what happens as U grows and makes the phase fluctuate more? Our continuum
description in terms of harmonic oscillators hides (but does not ignore) the fact that φ 'φ + 2π. The system admits vortices, aka winding modes.
34

From now on the background density n0 will not play a role and I will write ni for ∆ni .

114

Lattice T-duality. To see their effects let us do T-duality on the lattice. As usual, the
1 3, ....
dual
variables
live
on
the
bonds,
labelled
by
ī
=
2 2
Introducemī ≡φi+1 − φi
,
2π

Θī ≡

X

2πnj

(99)

j<ī

which together imply[mī , Θj̄ ] = −iδīj̄ .
To understand where these expressions come from, notice that the operatoreiΘī = eiPj<ī 2πnj

rotates the phase of the boson on all sites to the left of ī (by 2π). It inserts a vortex in
between the sites i and i + 1. The rotor hamiltonian is
2X
U X Θī+1 − Θī
Hrotors =−J
cos 2πmī
2
2π
ī
ī !
X U  ∆Θ 2 J
(100)
+ (2πmī )2
'
2 2π
2
ī

where the second term should be regarded as a π 2 kinetic energy term. This looks like simple
harmonic oscillators again. BUT: we must not forget that Θ ∈ 2πZ! It’s oscillators withdiscretized positions. We can rewrite it in terms of continuous Θ at the expense of imposing
the condition Θ ∈ 2πZ energetically by adding a term −λ cos Θ35 . The resulting model has
the action1(∂µ Θ)2 − λ cos Θ.
Leff =
2(2π)2 ρsThis Θ = φL − φR is the T-dual variable of our previous discussion; recall the ETCRs[φ(x), Θ(y)] = 2πisign(x − y).

(101)This follows directly from the definition of Θ (99). (101) means that the operator cos Θ(x)jumps the SF phase variable θ by 2π – it inserts a 2π vortex, as we designed it to do.
This system has two regimes, depending on the scaling dimension of the vortex insertion
operator:
35

This step seems scary at first sight, since we’re adding degrees of freedom to our system, albeit gapped
ones. Θī is the number of bosons to the left of ī (times 2π). An analogy that I find useful is to thefact that the number of atoms of air in the room is an integer. This constraint can have some importantconsequences, for example, were they to solidify. But in our coarse-grained description of the fluid phase, we
use variables (the continuum number density) where the number of atoms (implicitly) varies continuously.
The nice thing about this story (both for vortices and for air) is that the system tells us when we can’t
ignore this quantization constraint.

115

• If λ is an irrelevant coupling, we can ignore it in the IR and we get a superfluid, withalgebraic LRO.
• If the vortices are relevant, λ → ∞ in the IR, and we pin the dual phase, Θī = 0, ∀ī.
This is the Mott insulator, since Θī = 0 means ni = 0 – the number fluctuations are frozen.
When is λ relevant? Expanding around the free theory,heiΘ(x) e−iΘ(0) i =

c
x2πρs

this has scaling dimension ∆ = πρs which is relevant if 2 > ∆ = πρs . Since the bose1
, we see that only if η < 41 do we have acorrelators behave as hb† bi ∼ x−η with η = 2πρsstable SF phase. If η > 41 , the SF is unstable to proliferation of vortices and we end up in
the Mott insulator.
[End of Lecture 17]

116

CFT at the edge of Quantum Hall states. I must mention another important situation where the previous analysis is useful, namely at the interface between a quantum Hall
insulator and vacuum. Such an edge hosts a D = 1 + 1 CFT which is determined by the(gapped) bulk state. In simple examples (such as integer QHE and Laughlin states) it is achiral Luttinger liquid.
With hindsight, the nature of this CFT can be determined using very simple reasoning
which can be found, e.g. in the QFT textbook by Tony Zee. The basic idea is first to argue
that the effective field theory of the bulk is a Chern-Simons gauge theory; the edge states ofsuch a theory are then determinable, just as in the final section of Witten’s Jones Polynomial
paper.
The preceding results can then be used to understand, for example, tunneling of electrons
into a quantum Hall sample.
This connection is also the place where the most interesting CFTs may be realized innature.

117

4

CFT

I haven’t been very explicit about it so far, but we’ve been studying CFT for a while now.
I decided that it would be useful to understand the important examples well before saying
anything about the more formal aspects of the structure of CFT and constraints on physics
from conformal invariance.
So far we have been studying examples of CFT which have a free-field description (maybe
not always an obvious one). When this crutch is absent, what do we do? Remarkably, theanswer is not ‘nothing’. In particular we can define what we mean by CFT without recourse
to perturbation theory, and sometimes we can even solve such a thing (I will try to makeclear what ‘solve’ means here).
[D > 2: Rychkov. D = 2: Ginsparg.] Here is the way to think about CFT nonperturbatively: A CFT is a list of operators with definite scaling dimensions {Oa , ∆a } andtheir OPE structure constants, cabc , appearing in (83). From this information you can compute any correlation function of local operators by successive uses of the OPE. So the Isingmodel is36 :¯ ,···}
σ , µ , T ≡ χL ∂χL , T̄ ≡ χR ∂χ{|{z}
1 , χL , χR , ε ≡ iχL χR , |{z}|{z} ||{z} |{z} | {z }
{z
} |
{z R}
(0,0)

( 21 ,0) (0, 12 )( 12 , 12 )1 1
, 16 ) ( 1 , 1 )
( 16

(2,0)

16 16

(0,2)

and cχL χR  = i, cσ,ψ,µ = 1, ..., or more explicitly:σ × σ = 1 + cσσε ε, ε × ε = 1 + cεεε ε, σ × ε = cσσε σ.
|{z}=0Notice that I include the identity operator in this list; it has dimension zero for sure. In thisway of enumerating the operators the list goes on forever: e.g. we can keep appending morederivatives. We’ll learn to how do better soon by listing only the primary operators at the
top of each representation of the conformal group. But in fact, as in effective field theory,
one can often get far just knowing about the few lowest-dimension operators37 A special roleis played by the operator T (z) in this list, the stress tensor.

4.1

The stress tensor and conformal invariance (abstract CFT)

Let’s think a bit about relativistic field theory in the continuum. Any continuum QFT has
a stress tensor38 ; there are two useful perspectives on this operator. The simpler, but worse,
36

Actually, I should leave out the fermions. As this model arises from a spin system, a single fermion
operator is not allowed. From the point of view advocated in §2.2.7, it is not gauge invariant.
37
Naturally, this strategy is generally called ‘effective conformal field theory’. A recent victory in thisdirection can be found in this paper, which solves QCD in D = 2 by diagonalizing a 5 × 5 matrix (!).
38
Maybe this is not true. The following is a very technical comment which you should ignore if you want.
It is certainly true if the field theory has a definition in terms of a lagrangian or a local hamiltonian. Possible

118

one is to regard it as the Nöther current for spacetime translations. The stress tensor Tµν
is conserved if the action is translation-invariant: that is, invariant under the replacement
φ(xµ ) → φ̃(x̃) ≡ φ(xµ + aµ ). This is d + 1 symmetries, so we have d + 1 conserved currents:Tνµ , ∂µ Tνµ = 0. The associated charges are the energy and momentumZZ
H=
T00 , Pi =
T0i .
space

space

These generate translations in time and space by commutators (or classically by Poissonbrackets):δa φ(x) = φ̃(x̃) − φ(x) = aµ [Pµ , φ(x)].
The finite solution of this equation isµµφ(x + a) = eiPµ a φ(x)e−iPµ a .
The current we get from the Nöther method is not symmetric in its indices.
The better way to think about the stress tensor is as the (linear) response of the system toa local, small perturbation of the metric of spacetime:
Z
δS
1
√ µν
.
(102)
gT δgµν , aka T µν (x) = √δS =g δgµν (x)√√Here g ≡ det g is the covariant measure.39 Consider in particular making a scale transformation, which we can accomplish by changing the metric by:Z
δgµν = 2λgµν → δS = Tµµ 2λ(103)where λ is a constant. Therefore we see that if Tµµ = 0 then the theory is scale invariant.
In the other direction, scale invariance (vanishing of (103)) actually only implies that Tµµ =∂ µ Kµ is a total derivative. But the object K is a vector operator whose scaling dimension
must be D − 1 (since that of the stress tensor is D). On general grounds, a vector operatorof dimension D − 1 is a conserved current (at least this is the only way I know to protect its
dimension). And a conserved current has ∂ µ Kµ = 0, so that Tµµ = 0 anyway.
If Tµµ = 0 then nothing we said above depended on the fact that λ is a constant, and weshould also consider the following transformation:Zδgµν = 2δΩ(x)gµν → δS = Tµµ 2δΩ(x)(104)
exceptions come from: (1) field theory in a fixed AdS geometry; the absence of dynamical gravity means nostress tensor. (An attempt to take such a thing seriously as a CFT is in this paper by Heemskerk et al.) (2)In his CFT notes, Rychkov discusses a long-range Ising system with a fixed point with no stress tensor; heclaims without further discussion that this is the same as the first example. (3) Weird things ‘defined’ byscaling limits of string theory. Does a 2d CFT without a stress tensor have a central charge?
39
The T µν from the metric variation can be related to the Nöther current for translations by “improvement,” which means adding boundary terms to the action; this modifies the Nöther current.

119Such a transformation, which rescales the whole metric (and therefore preserves angles between vectors at the same point) in a position-dependent way, is a Weyl transformation,
closely related to a conformal transformation. The difference is that for arbitrary Ω(x),the new metric will be curved (R ∝ ∇2 log Ω). If we are not interested in QFT in curved
spacetime, we should restrict ourselves to choices of Ω which preserve the initial choice of curvature; this means that they can be undone by a coordinate transformation. Infinitesimally,such a transformation is xµ → (x0 )µ = xµ + ξ µ (x), and the metric change is!

δgµν (x) = ∂µ ξν (x) + ∂ν ξµ (x) = 2δΩgµν .

(105)This a set of PDEs for ξµ and δΩ.
For the case of gµν = ηµν , Minkowski spacetime, the stage of special relativity, the constraint(105) is solved by the following formulae40 (and by translations and rotations and boosts,
which don’t change the Minkowski metric at all). The conserved currents and charges of thetransformations above (in flat spacetime) are:ZZνdSµ = x Tµν→ D ≡ S 0 d x = xµ p µ
(106)
Z
2
λ
Cµν = (2xµ xλ − x gµλ )Tν→ Cµ ≡ C0µ dd x
(107)
since both ∂ µ Sµ and ∂ µ Cµν are proportional to Tµµ . In the last equation of the first line, pµ
is the momentum density; the action of the integrand on fields is via xµ ∂µ , a rescaling.
Here is the right way to think about this condition on
finite transformations. We are demanding a change of coordinates which accomplishes the following:!

ηµν (dx0 )µ (dx0 )ν = Ω2 (x)dxµ dxν ηµν .
The jacobian for this change of coordinates must therefore
satisfy∂(x0 )µ= Ω(x)Oµν (x)Jνµ ≡[Fig credit: Rychkov]∂xν
where Oµν (x) is a spacetime-dependent Lorentz transformation (OηO = η). This is aspacetime-dependent rescaling and rotation; you should think of it as a position-dependent
RG transformation.
We conclude from the above discussion that, at least classically, if Tµµ = 0 (and T µν issymmetric) then the theory has both scale invariance and conformal invariance. The precise
logical relation between scale invariance and conformal invariance is a subject of a lot of
discussion. There is no example of a reasonable relativistic fixed point without conformal
invariance.
40

The details are at the beginning of the Ginsparg notes, or on page 3 here.

120

The resulting set of transformations forms an extension of the Poincare group. There are
d + 2 extra generators Cµ and D. In fact the algebra is so(d + 1, 2).
Conformal algebra. Here is how to think about this algebra. Most of the commutators
just say that vectors (like Cµ and Pµ ) transform as vectors and D is a scalar. The importantones are:[D, Pµ ] = iPµ , Pµ is a raising operator for D[D, Cµ ] = −iCµ , Cµ is a lowering operator for D[Cµ , Pν ] = 2i (ηµν D − Mµν )which says you can recover D and the spin from Cµ .
The right way to think about the transformation C generates is:inversion ◦ translation ◦ inversion:

xµ
x0µ
=
+ bµ .
~x.~x
x~0 .x~0(108)µInversion is xµ → − xxν xν ; inversion plus poincaré implies conformal.
So if we want to study representations of this algebra, we can diagonalize D and some of
the rotation generators. (Notice that the familiar Poincaré casimir Pµ Pµ is not a casimir of
the conformal group.) Then we can build representation starting with states that have been
lowered as much as possible, so C|primaryi = 0; the rest of the representation is obtainedby acting with derivatives, i.e. Pµ (and SO(d) rotations).
Fields also form representations of the spacetime symmetry group. A primary field is one
which is local enough that its transformation rule under a conformal transformation is thesame as if it were just a scale transformation:
Φ(x) → Φ0 (x0 ) = Ω(x)−∆ Φ(x)for a scalar field, orΦ(x) → Φ0 (x0 ) = Ω(x)−∆ R[O(x)µν ]Φ(x0 )for a field with spin. R is the representation matrix for the rotation O in the Φ representation;
for a vector field, this is just O(x)µν .
Fields which are not primaries (for example, the derivative of a primary) are called descendants and have more complicated transformation rules. Fields can be organized into
eigenstates of D, of definite scaling; every such field is a primary or a descendant. The
algebra determines the behavior of descendants from that of primaries.

4.1.1

Geometric interpretation of conformal group

Here is a geometric interpretation of the conformal group, called the projective null cone
construction (due to Dirac (!)). The conformal group in Rd,1 is isomorphic to SO(d + 1, 2),121the Lorentz group of a space with two extra dimensions. This space Rd+1,2 has metricηab = diag(− + +... + −)ab

(109)

where the last two dimensions are the ‘extra’ ones. So we can find linear representations of
the conformal group by temporarily adding some extra coordinates. A light ray in this spacecan be parameterized by d + 1 dimensional coordinates xµ in the following way:1
1
ζ a = κ(xµ , (1 − x2 ), (1 + x2 ))a
(110)
2
2
where κ is some arbitrary constant. The group SO(d + 1, 2) moves these light rays around.
We can interpet these transformations as maps on the xµ . In fact these transformations(combined by rescalings to get back to the original slice) are precisely the conformal transformations.
Regard the expression (110) as an embedding of Rd,1 ⊂ Rd+1,2 ; this is an isometric embedding, i.e. the induced metric is again ds2 = −dt2 + d~x2 . To see that Lorentz boosts inthe embedding space are conformal transformations on the null slice, note that on the slice,ξ a ξa = 0 =⇒ ξ a dξa = 0, and so the transformation ξ a → λ(x)ξ a (which preserves a null
subspace) takes dξ a dξa → λ(x)2 dξ a dξa , a conformal transformation.
Invariants in Rd+1,2 should therefore be conformal invariants. Consider the object:
1(111)
ζ1 · ζ2 = ηab ζ1a ζ1b = κ1 κ2 (x1 − x2 )2 .
2
ζ a and λζ a are identified with the same xµ , so κ is a redundant variable. So conformal
invariants actually are cross ratios of invariants in Rd+1,2 , for example
ζ1 · ζ2 ζ3 · ζ4.
ζ1 · ζ3 ζ2 · ζ4

(112)

2
≡ (x1 − x2 )µ (x1 − x2 )µAn extremely useful consequence of this is the statement that r12transforms under a conformal transformation by a rescalingx212 →x212Ω(x1 )Ω(x2 )where Ω(xi ) ≡ 1+2b·xi +b2 x2i is the rescaling of the metric at xi under the special conformal
xµ +bµtransformation xµ → 1+2b·x+b2 x2 .

4.1.2

Infinite conformal algebra in D = 2.

Here there is an important division between D = 1 + 1 and D > 1 + 1. In D = 2, in lightcone
or holomorphic coordinates, tracelessness of T says
Tµµ ∝ Tzz̄ = 0 .
122Conservation of the stress tensor 0 = ∂ µ Tµν then says∂¯z Tzz = 0,∂z Tz̄z̄ = 0.

That is T (z) ≡ Tzz is holomorphic (and T̃ (z̄) ≡ Tz̄z̄ is antiholomorphic). In the quantum
theory, these statements are true as operator equations; that is: they are exactly true awayfrom other operator insertions in the path integral (the lingo for this is ‘up to contact terms’).
This holomorphic factorization has the following dramatic consequence: Given an arbitrary
holomorphic function41 ξ(z), the current
jµ(ξ) = (jz , jz̄ )µ ≡ (ξ(z)T (z), 0)µis also conserved (!):
∂ µ jµ(ξ) = ∂¯z jz − ∂jz̄ = 0.
This is infinitely many conserved currents! Basically, just from scale invariance.
What are these transformations? Recall that the current Tµν generates translations, which
by abuse of notation42 we can write as xµ → xµ +aµ . Accordingly, the current T (z) generates(ξ)‘holomorphic translations’: z → z + ξ with ξ constant. So it’s not too shocking that jµ
generates the local transformation z → z + ξ(z). The finite version of this transformation isjust an arbitrary holomorphic map:z → z 0 (z),z̄ → z̄ 0 (z̄).

(The important thing here is that z 0 does not depend on z̄!) What does this transformation
do to the metric? The flat metric in holomorphic coordinates is
ds2 = dτ 2 + dx2 = dzdz̄ →∂z ∂ z̄dz
dz̄ = f (z, z̄)dzdz̄.
∂z 0 ∂ z̄ 0

The metric has only changed by an overall function. This means that the angle between
any two vectors has not changed. This is the definition of a conformal transformation. The
conformal group is infinite dimensional in D = 2.
ForR the example of the free massless scalar in D = 2 (with curved-worldsheet action S[φ] =√1
dxdt gg µν (x)∂µ φ∂ν φ) the stress tensor is
2πK


1
1
2
Tµν =
∂µ φ∂ν φ − gµν (∂φ) .
2πK
2
41

More precisely, since there is one for z and one for z̄, ξ z , ξ z̄ are components of a holomorphic vector field.
This abuse of notation is both very tempting and very confusing. We are not merely relabeling ourcoordinates; that doesn’t do anything – physics is coordinate-independent. We are transforming our fieldsby
φ → φ0 (x0 ) = φ(x) .
42

If you promise to keep this in mind, then we can use the less cumbersome expressions below.

123Notice that it is traceless: g µν Tµν ≡ Tµµ = 0. In holomorphic coordinates, ds2 = dzdz̄, thenonzero components are
1
1
Tzz ≡ T (z) =∂z φ∂z φ =: ∂z φL ∂z φL :,
2πK
2πK
1 ¯ ¯
1
T̄z̄z̄ ≡ T (z̄) =∂z φ∂z φ =: ∂¯z φR ∂¯z φR :
(113)
2πK
2πK
In the last step I’ve emphasized the factorization into L and R parts, and the fact thatquantumly we must define this composite operator somehow, and we are doing it by normal
ordering. This involves a choice of additive constant, about which there is a bit more to say.
Conformal invariance constrains the operator algebra of a CFT, and since (various momentsof) the stress tensor components are generators of conformal transformations, their OPEs
are highly constrained.

4.2

Radial quantization

Here is an important example of a conformal transformation: Our coordinate on the cylinder
was w = x + iτ ; I call it a cylinder because x ' x + L. Let’s set L = 2π for convenience.
Consider the map
w 7→ z = e−iw = eτ −ix
Since it’s holomorphic (z(w, w̄) doesn’t depend on w̄), this is a conformal transformation.
The image is the complex z-plane. Equal-τ surfaces are circles. The spatial momentum
operator was ∂σ = z∂z − z̄ ∂¯z , which acts by z → eib z – rotations about the origin. The
hamiltonian on the cylinder was ∂τ = z∂z + z̄ ∂¯z . So time translations become rescaling
about the origin of the z-plane: the dilatation operator (the operator S that generates scaletransformations) is the hamiltonian in radial quantization. Time-ordered correlators on thecylinder are radially-ordered. One reason this is nice is that it provides an IR cutoff evenwhen studying infinite-volume physics. (Notice that the operators z∂z ± z̄ ∂¯z commute andso we can label states by their scaling dimension ∆ and spin.
So far, this all goes through for CFT in any D. In D > 2, the equal-radial-time-slices are
D − 1-spheres, and therefore the single quantum number for spin would be replaced by arepresentation of SO(D).
Example. The mode expansion of the free boson (about the origin of the z plane) is now43X∂z φ(z) =ρn z −n−1n∈Z(with ρ0 ≡ p), and we can extract them by a contour integral:Idz nz ∂φ(z)ρn =C0 2πi43

Allow me to use φhere =√1
πT φbefore , so that φhere ' φhere + 2πR and the action is S = − 4π

124

R

¯
d2 z∂φ∂φ.

where C0 is a contour encircling z = 0 (and no other operators!).
Primaries and quasiprimaries.
Def: A primary operator (or field) Φh,h̄ (z, z̄) of weight (h, h̄) transforms under the conformaltransformation44
h̄(z, z̄) → (f (z), f¯(z̄))
byΦh,h̄ (z, z̄) → (∂z f )h ∂¯z f¯ Φ(f (z), f¯(z̄)).
(The way to remember this is that Φ(z)dz h dz̄ h̄ is a scalar.) For example,
scaling:rotations:z → eλ z =⇒ Φ → eλ∆ Φ,
z → eiθ z =⇒ Φ → eisθ Φ,∆ = h + h̄ (scaling dimension)s = h − h̄ (spin)(114)

The infinitesimal transformation f (z) = z + ξ(z) results inδξ Φ(z) = (ξ∂z + h∂z ξ) Φ + antiholomorphic bits.
R dz
This transformation is generated by 2πi
ξ(z)T (z) ≡ L[ξ] (as in 4.1.2) in the sense thatδξ Φ = i[L[ξ], Φ].
Consider for a while a holomorphic operator, with h̄ = 0. It has a mode expansionXΦh (z) =Φn z −n−h .
n∈ZThe shift by h in the moding on the plane comes from the conformal transformation fromthe cylinder:XΦn e−iwn
Φh (w = −i ln z) =nwhere this is just fourier expansion, and the conformal factor isis positive energy.

∂z h= z h . Note that n < 0∂w

The modes of the stress tensor are called Virasoro operatorsIX
dz n+1−n−2
T (z) =Ln z
, Ln =z T (z).
2πiC0n
The definition of primary implies that
[Ln , Φh,h̄ (0)] = 0, ∀n > 0.

[L0 , Φh,h̄ (0)] = hΦh,h̄ (0),44Recall that in general dimension D ≥ 2, the transformation of a primary operator of dimension ∆ isΦ∆ (x) → |∂x0 ∆ 0 0
| Φ∆ (x )∂x0where | ∂x∂x | is the Jacobian of the conformal map.

125

This in turn implies that the state|Φh,h̄ i ≡ Φh,h̄ (0)|0i
is a highest weight state of the Virasoro operators, in the sense that
L0 |Φh,h̄ i = h|Φh,h̄ i,

Ln |Φh,h̄ i = 0, ∀n > 0.

Note that the modes with n > 0 raise the value of L0 . They include the ordinary special
conformal generators Cµ = (L1 , L̄1 )µ .
In D = 1 + 1 it is important to distinguish between Vir primary and ordinary conformal
primary, which is just killed by L1 and not Ln≥2 .
Contours and commutators. You may be bothered by the connection between thealgebra in terms of OPEs on the complex plane(like∂φ(z)∂φ(w) ∼1+ ··· )(z − w)2

and the perhaps-more-familiar algebra of mode operators.
(like[ρn , ρm ] = nδn+m ).

The very direct connection between the two comes from radial quantization. Recall that thepath integral on the plane produces radially ordered correlators:Z[Dφ]
e−S[φ] ...A(z)B(w).....
= hT (....A(z)B(w)...)i|{z}
these are numbers, order doesn’t matter

with(
A(z)B(w),
T (A(z)B(w)) =B(w)A(z),So consider, for example, the commutator
I

I
dz
[L[ξ], Φh (w)] =−ξ(z)T (z)Φh (w)2πi|z|>|w||z|<|w|I=dz
ξ(z)T (z)Φh (w).
Cw 2πiThis is a general rule: the commutator of modes of two
fields is given by the contour integral of one about theother.
126|z| > |w||z| < |w|

The previous expression is of interest since it describes
the conformal transformation of the primary Φh parameterized by the holomorphic vector field ξ:i[L[ξ], Φh (w)] = δξ Φh (w) = (ξ∂z + h∂z ξ)Φ.
Compare this expression to the general OPE
X 1
O(n) (0)T (z)Φh (0) =n+1
zn
and we determine


(1)

hΦh
O =
O(0) =∂Φh .

 (else)O= 0

That is, we have shown that the OPE of the stress tensor with a primary operator of dimension ∆ is
∆O ∂O
.
(115)
T (z)O(0) ∼ 2 +zz
Notice that L0 is scaling, L−1 is translations, and L1 is special conformal.
Some examples to check: (1) T (z)∂φ(0) for the free boson theory. (2) T (z)T (0) for the freeboson theory (this one is a trick question as we’ll see in 4.2.1).
Not all operators of definite scaling dimension behave this way, and this can be taken
as a definition of a primary operator. It implies that under a conformal transformation(z, z̄) → (w(z), w̄(z̄),
 h  h̄∂ w̄
∂w
Oh,h̄ (z, z̄) →Oh,h̄ (w(z), w̄(z̄))∂z∂ z̄State-operator correspondence. (any D)
To get a CFT state |Φi from a local operator Φ(z) just insert that operator at the origin
acting on the CFT vacuum in radial quantization about the origin. It’s easier in equations:
Φ(0)|0i ≡ |Φi.
The state |0i is then the image of the identity operator under this map.
To get an operator from a state, let’s think about the path integral on the ball Σ. Think of
it as a functional of the boundary conditions on the fields (which I’ll call φ; θ is a coordinateon ∂Σ):Z[Dφ]e−S = Ψ0 [φ0 ] = hφ0 |0iφ|∂Σ (θ)=φ0 (θ)

127This is an integral representation of the groundstate wavefunctional. If instead we consider
the path integral with a local operator insertion, we get a wavefunctional for a different state:Z[Dφ]e−S Φ(0) = ΨΦ [φ0 ] = hφ0 |Φi .
φ|∂Σ (θ)φ0 (θ)

In a CFT this is related by a conformal transformation z = e−iw to the path integral on thecylinder with the state |Φi inserted in the far past. So to get the operator corresponding toan arbitrary state just glue this path integral around where you want to put the operator. Toµµmove it to a different place, just act with translation generators: Φ(x) = e−iP xµ Φ(0)eiP xµ .
For more on the beauty of the state-operator correspondence in D-dimensional CFT see
the notes here Rychkov.
Convergence of OPE. With this realization in mind, it is clear that the OPE in CFT
is a convergent expansion: we are simply inserting a resolution of the identity on the hilbert
space in a particular basis of eigenstates of the dilatation operator.

Adjoint in radial quantization. The adjoint is a bit weird: it is just the adjoint
on the cylinder, but two things. First, the state in the far future on the cylinder gets mapped
to z = ∞. Second, we must remember that the map from the cylinder to the sphere produces
a Jacobian. Since the adjoint operation involves an inversion (a particular element of theconformal group), the primary Φ acquires an extra factor from the Jacobian. This produces
an extra factor in the state:|Φi† = lim hΦ|z −h .
z→∞This is sometimes called the BPZ adjoint after its discoverers.

4.2.1

The Virasoro central charge in 2d CFT

The OPE of the stress tensor with itself is
T (z)T (0) ∼∂T
cL /2 2T
+ 2 +
4
z
z
z

(116)

The quantity c appearing here is called the Virasoro central charge. The word ‘central’ isbecause it is a c-number, not an operator. It is a crucial piece of data about the CFT. It
can be extracted fromcL /2
hT (z)T (0)i = 4 .
z
There is also a right-moving central charge which appears in the OPE T̄ T̄ .
1
For free bosons, with T (z) = 4πK(∂φ)2 it is equal to cL = 1 (note that the Ks all cancelout). Notice that it is additive: with N free bosons, the answer is cL = N .

128

For a majorana fermion, with T (z) = χL ∂χL , it is equal to cL = 12 . A basic check ofbosonization is that a complex fermion has central charge 21 + 21 = 1.
It is monotonic under RG flows between CFTs: cU V > cIR . Therefore it is a useful measureof the number of degrees of freedom. In fact, using the definitionc = lim 2z 4 hT (z)T (0)iz→0it can be extended away from RG fixed points; Zamolodchikov proved that this quantity ismonotonic along any flow.
Comparing (115) and (116) we see that nonzero c means that the stress tensor itself is nota primary. It transforms weirdly under conformal transformations, in a very definite way.
This means, for example, that under the map from the cylinder to the plane, the additive
normalization of the stress tensor changes. For an application of this to Casimir energy of
CFT on a sphere S d , see HW 6.
For this reason c can be interpreted as a conformal anomaly or Weyl anomaly.
Tµµ =

cL
R.
4π

The relation between this statement and the stress tensor OPE is the same as that betweenF
the chiral anomaly in the form ∂ µ jµA = k 2π
and the anomalous commutator in (97), or betterthe current-current OPEkj(z)j(0) ∼ 2 .
z
In the presence of a source coupling to the current (A for j, or curvature for Tµν ), thestatement about the two-point function implies a statement about the one-point function.

4.3

Back to general dimensions.

The conformal group in D > 2 is finite dimensional, but it still produces powerful constraints.
So we can make our list of operators specifying the CFT much shorter by simply enumerating the primaries (and their dimensions, spins, and structure constants). All the structure
constants for the descendants are determined from this data by conformal invariance.

4.3.1

Constraints on correlation functions from CFT(any D) We consider correlators of primaries; to get correlators of descendants just takederivatives of those of their primaries. We’ll also focus for simplicity on scalar operators. A
129

correlation function of N primaries transforms under a conformal transformation asYYYh Φ0i (x0i )i =Ω(xi )−∆i h Φi (xi )i.
iii

You should regard the operators Φ0 on the LHS as the images of Φ under a spacetimedependent RG transformation.
Suppose we have in our possession a conformal invariant I which depends on N spacetime
positions.
• Translations imply that I depends only on differences xi − xj (of which there areD(N − 1)).
• Rotations imply that I depends only on distances rij ≡ |xi − xj | of which there are
N (N −1).
2rij• Scale invariance implies that I depends only on ratios of differences rkl.

• Special conformal transformations act on these distances by
0 2
(r12
) =

2
r12
Ω1 Ω2

where
Ωi ≡ 1 + 2b · xi + b2 x2i
is the factor by which the metric at xi rescales under the associated transformation.
(The easiest way to see this is by doing it for an inversion first, and then using theconstruction of special conformal transformations as I ◦ T ◦ I. Alternatively, this iswhere the projective null cone comes into its own.) Only cross-ratios
rij rkl
rik rjl
are invariant under this. There are N (N2−3) of these (and none at all for N < 4).
This discussion of invariants has the following implications.
1. One-point functions of primaries of nonzero dimension vanish if there is a conformallyinvariant vacuum, D|0i = 0. That this is so can be seen as follows. The infinitesimal
scale transformation of a primary isδΦ(x) = Ω−∆ Φ(x0 ) − Φ(x) |Ω=e−λ ,x0 =e−λx = λ (∆ + xµ ∂µ ) Φ(x) = iλ[D, Φ(x)].
The last step is the statement that the dilatation operator D generates scale transformations. Therefore:h0|Φ∆ (0)|0i =ih0|[D, Φ∆ (0)]|0i = 0.
∆130

2. In CFT, two point functions of primaries are only nonzero if the primaries have thesame dimension. This requires conformal invariance; the rest is determined by scale,translation and rotation:
δh1 ,h2h0|φ1 (z)φ2 (w)|0i = 2h
.
r12 1
3.

3
Y
2(−h −h +h ) 2(−h −h +h ) 2(−h −h +h )h φi (zi )i = C123 r12 1 2 3 r23 2 3 1 r31 3 1 2
i=1

C123 is the OPE coefficient between the three operators. Notice that this is more thanscaling, which just says that the sum of the powers should be h1 + h2 + h3 .
4.

4
YY −2(hi +hj +P h /3)k k
h φi (zi )i = f (x1 , x2 )
rij
i=1

i<j

where x12 are the two cross-ratios. In two dimensions, these arex=

z12 z34z13 z24and x̄.
It seems like N ≥ 4-point functions are underconstrained. However, we have see that thesuccessive use of the OPE relates them (via known functions) to sums of 3-point functions.
The unknown functions of cross-ratios are actually determined by this process! The generalname for such objects is conformal blocks.

4.3.2

Thermodynamics of a CFT

Thermodynamics of scale-invariant theories is very constrained. The partition function isdefined as
ZCF T = TrCF T (exp(−H/T )) .
In the thermodynamic limit, ln Z is extensive, i.e. proportional to the volume of the space.
But ln Z is a dimensionless quantity. Hence, we must have ln Z ∼ V T d (d is the number ofspatial dimensions) in the absence of any other energy scales (such as a chemical potential
for some conserved charge). The free energy then will be
F = −T ln Z = cV T d+1 .
where this c should be regarded as a rough measure of the number of degrees of freedomof the CFT. In d = 2 it is proportional to the Virasoro central charge. In other cases, thisstatement is less sharp.

131As with any conservation law, Tµµ = 0 is an operator equation in the full quantum theory.
What happens if we put it inside Tr(e−H/T ), with no other insertions? The operator equation
then translates into the following equation
0 = Tr(Tµµ e−H/T ) = hT00 i − hTii i = E − dP.
This last relation gives the speed of the sound
s
r

∂P
1
cs =
=
∂E S
d

4.3.3

(117)

Too few words about the Conformal Bootstrap

Earlier I said that given a ‘solution’ to a CFT in the form of a list of primaries and theirdimensions and OPE structure constants, you could compute any n-point correlator byrepeated use of the OPE to reduce it to a two-point function, by the following kind ofoperation:[from Rychkov’s CFT notes]

Some interesting questions and previews of the answers:(1) Which sets of scaling dimensions arise from CFTs which actually exist?
The general answer is not known. One set of constraints follows from unitarity in the formof positivity of the inner product on the CFT Hilbert space. For example, the dimension ofa scalar operator (other than 1) must be larger than the free field dimension ( D−2)45 .
2(2) In implementing the reduction to three-point functions and structure constants, youhave a choice about the order in which you group the operators. Do you get the same answer
independent of the order?
45

This follows from the fact that D appears on the RHS of [Cµ , Pν ]. A crucial ingredient comes from thefact that the dagger of a lowering operator is a raising operator (obvious, right?). In radial quantization this
implies the weird-looking formula:P†µ = Cµ .
For more details, I recommend the discussion leading up to eqn (5.57) of Jared Kaplan’s AdS/CFT notes.

132

The condition that you do – the associativity of the OPE – imposes significant constraints
on the structure constants and dimensions.

(3) Are there more constraints from higher-order diagrams?
Claim: no. The analog of the ‘pentagon identity’ is automatic.
To make use of the (many!) associativity constraints, a further ingredient is required, whichis called conformal blocks. For more in this direction, a good starting point is [Rychkov] andthe CFT lectures cited above.
In D = 2 it really works [BPZ]. This approach has led to complete solutions of all unitary
CFTs with c < 1 (this is explained in Ginsparg’s notes), and a number of interesting examples(e.g. Liouville theory) with c > 1. More recently, this program has had success in 2 + 1
dimensions: the 3d Ising model has been cornered.
[End of Lecture 18]

133

5

D =2+1

5.1

Duality for the transverse-field Ising model in D = 2 + 1

On a lattice of any dimension, with qbits at the sites r (H = ⊗r H 1 ), the tranvserse-field2
Ising model Hamiltonian is:X
X
HTFIM = −J 
Zr Zr0 + g
Xr  .
hrr0 i

The Z2 symmetry operator is S =Qrr Xr . J > 0 is a ferromagnetic coupling.

At small g, the groundstates (into which the cat states decohere) spontaneously break the Z2 symmetry; at g = 0, they are simply | ↑↑↑ · · · i and| ↓↓↓ · · · i. There is a gap of order 2J to flipping a spin. Small g mixes instates where some spins are flipped.
So far, everything is as in D = 1 + 1. The first difference is that in
D > 1, a domain wall is not a pointlike object. A domain wall has
codimension one – its location is specified by one real condition. So in
D = 2 + 1, it is a string. More precisely, since it is the boundary ofa region (of down spins in a sea of up spins, say), it is a closed string:
the boundary of a region has no boundary. This is a similar topologicalstatement to the fact that they come in pairs in one dimension.
Q
Clearly at g = ∞, the groundstate is a paramagnet, r | →i. What
happens in between? Can we describe the paramagnet in terms of theproliferation of domain walls as we did in D = 1 + 1?
Wegner duality. The theory can be reformulated in terms of domain walls as follows.
Domain walls cover links of the dual lattice, Γ̂, whose sites correspond to faces of the direct
lattice Γ. Each links of Γ crosses a link of Γ̂.
We can rewrite the ferromagnetic interaction term byτl̄x ≡ Zr Zr0
where ¯l is the link of Γ̂ crossed by the link rr0 . But on a 2d
lattice with N sites, there are 2N links, and hence 2N such τ s– way too many variables! The domain wall variables satisfy
the following local constraints. For each plaquette (little loop)
in the direct lattice, we have:
1 = (Z1 Z2 )(Z2 Z3 )(Z3 Z4 )(Z4 Z1 )134

But in terms of the domain wall operators, this plaquette is associated with a site +, andwe haveY
1 = (Z1 Z2 )(Z2 Z3 )(Z3 Z4 )(Z4 Z1 ) =τl̄x ≡ A+ , ∀ sites + ∈ Γ̂ .
l̄∈+In 1d, there were no loops (except with periodic boundary conditions on the spins).
What about the transverse field term in HTFIM ? Xr flips
a spin at r; on a sea of up spins, this creates a domain wall
surrounding the single resulting down spin. If there was alreadya domain wall there, it annihilates it (X is its own inverse).
More generally it moves the boundaries of the domain walls
adjacent to r. It flips the value of the domain wall operatorson all bonds surrounding r:Yτl̄z ≡ B r
Xr =
¯ r
`∈
where r denotes the set of links of the dual lattice surroundingr.
An important observation: each term in HTFIM commutes with the A+ = 1 constraint:•[B r , A+ ] = 0,∀r, +since each + and 2 share two links (X1 X2 Z1 Z2 = (−1)2 Z1 Z2 X1 X2 ).
• since ZZ = τ x is made from τ x s, just like A+ , the ferromagnetic term is unproblematic.
In the domain wall variables, we have
HTFIM = −J Xτl̄x + g

l̄∈Γ̂

X Y

τlz 

2 of Γ̂ l̄∈2

with the constraint that 1 = A+ = l∈+ τlx , ∀+. This is Ising gauge theory. The gaugetransformation acts asτīzj̄ 7→ sī τīzj̄ sj̄ , τīxj̄ 7→ τīxj̄ , sī = ±1.
Q

The transformation with sī = −1 and all others +1 is generated by A+ (ī), in the sense that
O 7→ A†+ (ī)OA+ (ī) .
And since A+ = 1 on the Hilbert space, this is a (huge!) redundancy of the description. One
more time: it does not act on the physical Hilbert space, which is
Hphys = {|ψi | A+ (ī)|ψi = |ψi, ∀ī} ⊂ ⊗links of Γ̂ H 1 = Hbig .
2135

Global issues. As in the 1+1d case, there are some global subtleties with the quantum
duality transformation: In particular, we must specify the value of a single spin in order forthe locations of domain walls to specify the spin configuration. The gauge theory variables
do not produce the symmetry-breaking degeneracy.
Notice by our counting argument (2N links between N sites) that thisQ means that someof the constraints A+ = 1 must also be redundant. An example is all sites of Γ̂ A+ = 1,
automatically, since each link is included twice.
This mismatch is a crucial thing: the spin system can have symmetry-breaking order andthis gauge theory cannot; as we will see below, the gauge theory can have topological orderand this spin system cannot.
From now on in this discussion, we will take the gauge theory as ‘fundamental’ and ignore
these issues. I will therefore stop writing the bars on the coordinates of Γ̂, the lattice which
the gauge theory inhabits.

5.2

Z2 lattice gauge theory

To make it look more like gauge theory familiar from high energy physics, regard the nongauge-invariant variable τ z as
R j̄

τīzj̄ = ei ī ~a·d~s“

00

the holonomy of some fictitious continuum gauge field integrated along the link. Moreprecisely, letτīzj̄ ≡ eiπaīj̄ , aī,j̄ = 0, 1.
Then the plaquette operator isPYHzi ∂2 ~a·d~liπ ◦ ∂2 a Stokes iπb2
B2 =τl “ = e
”=e
= e
,
l∈2

where b2 is the (discrete)
P magnetic flux through the plaquette 2. In the penultimate expression, the symbol ◦ is intended to emphasize that we are summing the as around theclosed loop.
In the Hamiltonian description of gauge theory, the field momentum for ~a is the electricfield ~e. So, we callτlx ≡ eiπel
The star operator isA+ =Yτlx = eiπ

P

l∈+

136

l∈+ eij≡ eiπ∆·ewhich is a lattice divergence operator. The constraint isY1=
τlx ↔ ∆ · e = 0 mod 2.
l∈+This is binary electrodynamics, electrodynamics mod two. Electric charges are violations of
the Gauss’ Law constraint: if(∆ · e) (i) = 1 mod 2
at some site i, we say there is a Z2 charge at site i. Notice that this is not something we cando in the spin system: such a site is the end of a domain wall.
Emergence of gauge theory from a spin system. We can perform the same toric code
trick as we did in D = 1 + 1 (in §2.2.7) and realize the Gauss’ Law constraint energetically,
on a larger Hilbert space, on which the local gauge transformations actually do act. That
is: The physics of the gauge theory is the same as the low energy physics (E  Egiant ) of
X
X
X
HTC = −Jτl̄x − Jg
B2 − Egiant
A+ .
l̄∈Γ̂

2 of Γ̂+∈Γ̂This maneuver (due to Kitaev and Wen) is conceptually quite important, since it meansthat the gauge theory can emerge from an ordinary spin system, with hilbert space Hbig =⊗links of Γ̂ H 1 .
2Which states satisfy the star condition? In the τ xbasis there is an extremely useful visualization: wesay a link l of Γ̂ is covered with a segment of string(an electric flux line) if el = 1 (so τlx = −1) and is
not covered if el = 0 (so τlx = +1):

These strings are just the domain walls across which
the spins disagree. The star condition is precisely thecondition that the strings don’t end.

5.3

Phase diagram of Z2 gauge theory and topological order

Confinement. g  1 At small g we must minimize −τlx for each l, subject to the constraint
Qxxl∈+ τl = 1, ∀+. This is easy, just set τl = 1 for every link. In the dual spin variables, thisis one of the states with all the spins the same, but we can’t specify which one, so there isa unique state of the gauge theory. In terms of the binary electromagnetic fields, this is a137state where hel i = 0, ∀l: electric flux is costly in energy. (This means that the magnetic field
is wildly fluctuating.)
This is the confined phase. The origin of this termcomes from asking what happens if we insert external
test charges, separated by a distance x, and ask what force one exerts on the other via
the gauge interactions.
Choose two sites i1,2 separated by x lattice sites, and demand thatQx
A+ (i1,2 ) = l∈+(i1,2 ) τl = −1, while A+ (i) = 1 for all other sites i. The latter condition
implies that the links covered with electric flux (el = 1, τlx = −1) form a string, which only
terminates at the external charges. This is the familiar statement that lines of electric flux
terminate and originate at electric charges.
But now to find the potential between
static charges V (x), we need to minimizeXτlx = E0 + 2JL(string) .
H(g = 0) = −Jl

Here E0 = −J2N is the energy of the statewith no electric flux and no external charges, where 2N is the number of sites. L(string) is
the length of the electric flux string: the string can be said to have a nonzero tension (energy
per unit length), 2J. Clearly this minimization is accomplished by a straight line, and thepotential between the charges is
V (x) = +2Jx
which is linearly rising with separation, and implies a constant attractive force
F = −∂x V = 2J.
Deconfinement. g  1 At large g, we minimize the magnetic flux termYτlz
−gJB2 = −gJl∈2

subject to the constraint A+ = 1, ∀+. The groundstate has B2 = 1, ∀2: the magnetic flux
is zero (mod 2). (Note that this satisfies the constraint, since B and A commute and cantherefore be simulataneously diagonalized.)
The excitations arise by flipping B2 = −1 for some 2; this costs energy 2Jg.

This is definitely a different phase from the g  1 phase above, as we will see now, but there
is not local order parameter (note that here the global disagreement in the duality transformation is important) which distinguishes them (this was Wegner’s original motivation forstudying this system).
Definite b means wildly fluctuating e. This is one reason to call this state deconfined:
electric flux is cheap. It will be useful to write the groundstate explicitly in the electric138

variables. The states of the electric variables satisfying the A+ = 1 constraint are describedby configurations of closed unoriented strings C on the lattice, just as before. Explicitly: alink with e = 1 is covered with a string segment.
We can write the groundstate in this basis:
X|gsi =
Ψ(C)|CiC

and we will now determine the groundstate wavefunction Ψ(C). The eigenvalue equation for
B2 is thenX!
|gsi = B2 |gsi =Ψ(C)B2 |Ci .
C

So: how does B2 act on a closed string configuration?
This is just the discussion we had above when we determined (in §5.1) that Xr = B r – the domain walls are theloops of electric flux, and they are moved around by themagnetic field operator:
B2 |Ci = |C + 2i

B2 creates a loop of string around the plaquette 2. If one is there already it annihilates it.
If a string has a segment covering 2, B2 moves the string around 2.

So the eigenvalue equation B2 = 1 says Ψ(C) = Ψ(C 0 ) if C 0 and C can be continuously
deformed into each other by attaching or removing plaquettes. If the lattice is simply connected – if all curves are the boundary of some region contained in the lattice – then thismeans
X|gsi =|CiC

is a uniform superposition of all loops.
String condensation. Notice that the deconfined phase of the gauge theory involves thecondensation of the electric flux strings, in the sense that the operators B2 which create
these strings have a nonzero groundstate expectation value:g=∞hgs|B2 |gsi = 1.
As with an ordinary condensate of bosons, away from the zero-correlation-length limit (g =∞), the condensate will not be exactly 1, since finite g suppresses configurations with electric
flux. But within the deconfined phase it will be nonzero.
Topological order. If the space has non-contractible loops, then the eigenvalue equation
does not determine the relative coefficients of loops of different topology! On a space with
2g independent non-contractible loops, there are 22g independent groundstates.

139

This is (a paradigmatic example of) topological order. No local operator can tell the different groundstates apart! A (nonlocal!) order parameter for the phase transition from small gto large g is the (integer!) number of groundstates on a space with nontrivial topology. Twogapped states where this number is different certainly cannot be adiabatically connected,
since an integer cannot vary continuously.
Consider the cylinder. There is one nontrivial class of loops,
call a representative γ. Let η be a line running along the cylinder. The two groundstates are generated by the action of the‘Wilson loop operator’Y
Wx (η) ≡τlx
l crossed by ηin the sense that 46|gs2 i = Wx (η)|gs1 i .
This is also a groundstate (at g = ∞) since there is no plaquette Q
which violates Bp (moresimply: [Hg=∞ , Wx (η)] = 0). They are distinguished by Wz (γ) ≡ l∈γ τlz in the sense thatWz (γ)|gsα i = (−1)α |gsα i, α = 1, 2.
This follows since {Wx (η), Wz (γ)} = 0 – they share a single link (the one pointed to by theyellow arrow in the figure).
At finite g (and in finite volume), there is tunneling between the topologically degenerate
groundstates, since thenY[H,τlz ] 6= 0.
l∈γThis means
hgs2 |H|gs1 i ≡ Γ 6= 0.
However, the amplitude Γ requires the creation of magnetic flux on some plaquette (i.e. aplaquette P with BP = −1, which costs energy 2Jg), which then must hop (using the −Jτ xterm in H) all the way along the path η, of length L, to cancel the action of Wx (η). The
amplitude for this process goes likeLJhgs2 | (−Jτ1x ) (−Jτ2x ) · · · (−JτLx ) |gs1 i∼= e−L log 2g
Γ∼2gJ · 2gJ · . . . 2gJ
2gJ
46

Notice that in the literal dual of the spin system,Yspinsτlx = (Z1 Z2 ) (Z2 Z3 ) · · · (ZN −1 ZN ) = Z1 ZN
l crossed by η

this operator is trivial. (If we identify the two ends of the cylinder to make a torus, Z1 = ZN , it is just1.) Here is where the global mismatch between the spin system and the ‘dual’ gauge theory is crucial. The
TFIM does not have topological order. However, the toric code construction is one demonstration that it ispossible for an ordinary spin system to exhibit this phenomenon.

140which is extremely tiny in the thermodynamic limit. The way to think about this is thatthe Hamiltonian is itself a local operator, and cannot distinguish the groundstates from eachother. It takes a non-perturbative process, exponentially suppressed in system size, to create
the splitting.

141

5.4

Superfluid and Mott Insulator in D = 2 + 1

Mean field theory gives the famous picture at right,
with lobes of different Mott insulator states with different (integer!) numbers of bosons per site. We haveseen above (in §3.1.4) that the prevention of vortices
is essential to superfluidity, which is the condensation of bosons. In D = 1 + 1, vortices are events inspacetime. In D = 2 + 1, vortices are actual particles,
i.e. localizable objects, around which the superfluid
phase variable winds by 2π (times an integer).
More explicitly, if the boson field which condenses(Fig credit: Roman Lutchyn)is b(x) = veiφ , and we choose polar coordinates inspace x + iy ≡ Reiϕ , then a vortex is a configuration of the order parameter field of theR→∞form b(x) = f (R)eiϕ , where f (R) → v far away: the phase of the order parameter winds
R→0around. Notice that the phase is ill-defined in the core of the vortex where f (R) → 0.
To see the role of vortices in destroying superfluidity more clearly, consider superfluid flow
in a 2d annulus geometry, with the same polar coordinates x + iy = Reiϕ . If the superfluid
phase variable is in the configuration φ(R, ϕ) = nϕ, then the current is~~ = ϕ̌ρs n .
J(R,ϕ) = ρs ∇φ
2πR
The current only changes if the integer n changes. This happens if vortices enter from theoutside; removing the current (changing n to zero) requires n vortices to tunnel all the way
through the sample, which if they are gapped and the sample is macroscopic can take acosmologically long time.
There is a dual statement to the preceding two paragraphs: a state where the bosons
themselves are gapped and localized – that is, a mott insulator – can be described by thecondensation of vortices. To see this, let us introduce dual variables to the 2 + 1d rotor
model
X
X
Hrotors = U
n2i − Jcos (φi − φj ) .
i

hiji

Introduce a dual lattice whose sites are (centered in) the faces of the original (direct) lattice;each link of the dual lattice crosses one link of the direct lattice.
φ −φ• First let eīj̄ ≡ i2π j . Here we define īj̄ by the right hand rule:ij × īj̄ = +ž (ij denotes the unit vector pointing from i to j). This is~ 1 . Defining lattice derivatives ∆x φi ≡a lattice version of ~e = ž × ∇φ2πφi − φi+x̌ , the definition is ex = ∆2πx φ , ey = ∆2πy φ . It is like an electric field
vector.
142• The conjugate variable to the electric field is aīj̄ , which must therefore
be made from the conjugate variable of φi , namely ni : [ni , φj ] = −iδij .
Acting with ni translates φi , which means that it shifts all the eīj̄ fromthe surrounding plaquettes. More precisely:
2πni = a1̄2̄ + a2̄3̄ + a3̄4̄ + a4̄1̄ .
1 ~This is a lattice, integer version of n ∼ 2π
∇ × a · ž. In terms of these variables,

UX
Hrotors =
2 i



∆×a
2π

2
−J

X

cos 2πeīj̄hīj̄i1~with the followingconstraint.
 If it were really true that ~e = 2π ž × ∇φ, with single-valued φ,~ · ~e = ∇~ · ž × ∇φ~then ∇= 0. But there are vortices in the world, where φ is not single

valued. The number of vortices nv (R) in some region R with ∂R = C is determined by thewinding number of the phase around C:ZIStokes~~ · ~e~d2 x ∇d` · ∇φ = 2π2πnv (R) =RC~ · ~e = zij ∂i ∂j φ = [∂x , ∂y ]φ clearly vanishes if φ is single-valued.) Since(More explicitly, 2π ∇this is true for any region R, we have
~ · ~e = 2πδ 2 (vortices).
∇
Actually, the lattice version of the equation has more information (and is true) because it
keeps track of the fact that the number of vortices is an integer:~ · ~e(ī) = 2πnv (ī),∆x ex + ∆y ey ≡ ∆nv (ī) ∈ Z.

It will not escape your notice that this is Gauss’ law.
~ varies continuously, i.e. electric flux is not
Phases of the 2d rotors. Since ~e ∼ ∇φquantized, this is called noncompact electrodynamics. Again we will impose the integer?
constraint a ∈ 2πZ energetically, i.e. let a ∈ R and add (something like) ∆H = −t cos a andsee what happens when we make t finite. The expression in the previous sentence is not quite~ · ~e − 2πnv = 0
right, yet, however: This operator does not commute with our constraint ∆47– it jumps ~e but not nv
47

A set of words which has the same meaning as the above: cos a is not gauge invariant. Understanding~ · ~e − 2πnv as the generator of a transformation,these words requires us to think of the operator G(ī) ≡ ∆XδO =s(ī)[G(ī), O].
īIt can be a useful picture.

143

We can fix this by introducing explicitly the variable which creates vortices:[nv (ī), χ(j̄)] = −iδīj̄ .
Since nv ∈ Z, χ ' χ + 2π lives on a circle. So:
H∼

X
ī

U
2



∆×a
2π

2

!
J
+ (2πe)2 − t cos (∆χ − a)2

~ · ~e = 2πnv .
still subject to the constraint ∆Two regimes:
J  U, t : This suppresses e and its fluctuations, which means a fluctuates. The fluctuating
a is governed by the gaussian hamiltonian

X
2
2
~
H∼
~e + b, which should look familiar. This deconfined phase has a gapless photon; awith b ≡ ∆×a2π
2 + 1d photon has a single polarization state. This is the goldstone mode, and this regime
describes the superfluid phase (note that the parameters work out right in the original variables). The relation between the photon a and the original phase variable, in the continuumisµνρ ∂ν aρ = ∂µ φ.
t  U, J : In this regime we must satisfy the cosine first. Like in D = 1 + 1, this canbe described as the statement that vortices condense. Expanding around its minimum, thecosine term is
h 3 t (a − ∂χ)2
– the photon gets a mass by eating the phase variable χ. There is an energy gap. This is
the Mott phase.
If the vortices carry other quantum numbers, the Mott phase can be more interesting.
Compact electrodynamics in D = 2 + 1. Note that this free photon phase of D = 2 + 1
electrodynamics is not accessible if e is quantized (so-called compact electrodynamics) where
monopole instantons proliferate and gap out the photon. This is the subject of §5.5.

5.5

Compact electrodynamics in D = 2 + 1

Consider a quantum system on a two-dimensional lattice (say, square) with rotors Θl ≡
Θl + 2πm on the links l. (Think of this as the phase of a boson or the direction of an144

easy-plane spin.) The conjugate variable nl is an integer[nl , Θl0 ] = −iδl,l0 .
Here nij = nji , Θij = Θji – we have not oriented our links (yet). We also impose the Gauss’law constraintXGs ≡nl = 0 ∀ sites s,l∈v(s)

where the notation v(s) means the set of links incident upon the site s (‘v’ is for ‘vicinity’).
We’ll demand that the Hamiltonian is ‘gauge invariant’, that is, that
[H, Gs ] = 0∀s. Any terms which depend only on n are OK. The natural
single-valued object made from Θ is eiΘl , but this is not gauge invariant. A combination which is gauge invariant is the plaquette operator,
associated to a face p of the lattice:Yye(−1) iΘl ≡ ei(Θ12 −Θ23 +Θ34 −Θ41 )l∈∂p– we put a minus sign on the horizontal links. ∂p denotes the links running around theboundary of p. So a good hamiltonian is
!
X
X
UX 2
cos
(−1)y Θl .
n −K
H=
2 l l
2
l∈∂2
A useful change of variables gets rid of these annoying signs. Assume
the lattice is bipartite: made of two sublattices A, B each of which only
touches the other. Then draw arrows from A sites to B sites, and let
(
+1,
i∈A
eij ≡ ηi nij.
, ηi ≡aij ≡ ηi Θij−1, i ∈ B
Then the Gauss constraint now reads
0 = eī1̄ + eī2̄ + eī3̄ + eī4̄ ≡ ∆ · e(ī).
This is the lattice divergence operation. The plaquette term readscos (Θ12 − Θ23 + Θ34 − Θ41 ) = cos (a12 + a23 + a34 + a41 ) ≡ cos (∆ × a)
– the lattice curl (more precisely, it is (∆ × a) · ž). In these variables,
H=

X
UX 2
el − Kcos ((∆ × a) · ň2 )2 l
2(in the last term we emphasize that this works in D ≥ 2 + 1 if we remember to take thecomponent of the curl normal to the face in question). This is (compact) lattice U(1) gauge
145theory, with no charges. The word ‘compact’ refers to the fact that the charge is quantized;
the way we would add charge is by modifying the Gauss’ law to∆ · e(ī) = charge at ī
| {z } | {z }
=⇒ ∈Z

∈Z

where the charge must be quantized because the LHS is an integer. (In the noncompact
electrodynamics we found dual to the superfluid, it was the continuous angle variable which
participated in the Gauss’ law, and the discrete variable which was gauge variant.)
What is it that’s compact in compact QED?
The operator appearing in Gauss’ law~ · ~e(x) − 4πn(x)G(x) ≡ ∇(here n(x) is the density of charge) is the generator of gauge transformations, in the sensethat a gauge transformation acts on any operator O by
O 7→ e−iPx α(x)G(x)OeiPx α(x)G(x)(118)

This is a fact we’ve seen repeatedly above, and it is familiar from ordinary QED, where using
the canonical commutation relations[ai (x), ej (y)] = −iδ ij δ(x − y), [φ(x), n(y)] = −iδ(x − y)(φ is the phase of a charged field, Φ = ρeiφ ) in (118) reproduce the familiar gauge transformations~~a → ~a + ∇α,
φ→φ+α .
SO: if all the objects appearing in Gauss’ law are integers (which is the case if charge isquantized and electric flux is quantized), it means that the gauge parameter α itself only
enters mod 2π, which means the gauge transformations live in U(1), as opposed to R. So
it’s the gauge group that’s compact.
This distinction is very important, because (in the absence of matter) this model does nothave a deconfined phase! To see this result (due to Polyakov), first consider strong coupling:U  K : The groundstate has el̄ = 0, ∀¯l. (Notice that this configuration satisfies theconstraint.) There is a gap to excitations where some link has an integer e 6= 0, of orderU . (If e were continuous, there would not be a gap!) In this phase, electric flux is confined,
i.e. costs energy and is generally unwanted.
U  K : The surprising thing is what happens when we make the gauge coupling weak.
146

Then we should first minimize the magnetic flux term: minimizing − cos(∆ × a) means ∆ × a ∈ 2πZ. Near each minimum,
the physics looks like Maxwell, h ∼ e2 + b2 + · · · . BUT: itturns out to be a colossally bad idea to ignore the tunnelling
between the minima. To see this, begin by solving the Gauss law constraint ∆ · e = 0 by
introducing
1(φ2 − φ1 )e1̄2̄ ≡
2π
1(i.e. ~e = ž · ∆φ 2π.) φ is a (discrete!) ‘height variable’. Then the operatorei(∆×a)(ī)
increases the value of eīa for all neighbors a, which means it jumps φ → φ + 2π. So we should
regard(∆ × a) (ī) = Πφ (ī)
in the sense that[Πφ (r), φ(r0 )] = −iδrr0 .
The hamiltonian is now
H=XUX
cos Πφ (r)(∆φ)2 − K
2 lr

with no constraint. Deconfinement limit should be K  U , in which case it looks like we
can Taylor expand the cosine cos Πφ ∼ 1 − 21 Π2φ and get harmonic oscillators. But: tunnelingbetween the neighboring vacua of ∆ × a is accomplished by the flux-insertion operator (ormonopole operator)
eiφ , [eiφ(r) , (∆ × a) (r0 )] = eiφ(r) δrr0
–that is, eiφ is a raising operator for ∆ × a. To analyze whether the Maxwell limit survives
this, let’s go to the continuum and study perturbations of the free hamiltonianZ   2
U ~
K 2
H0 =
∇φ + Πφ2
2
byZ
H1 = −

V0 cos φ .

Notice that again we can think of the addition of this term as energetically imposing thecondition that φ ∈ 2πZ.
So: is V0 irrelevant? Very much no. In facthcos φ(r) cos φ(0)i0 ∼ const(119)

has constant amplitude at large r! That meansthat the operator has dimension zero, andR
the perturbation in the action has [S1 = − V0 cos φd2 xdτ ] ∼ L3 , very relevant. The result147

is that it pins the φ field (the would-be photon mode) to an integer, from which it can’t
escape.
To see (119) begin with the gaussian identity
ss00

heisφ(x) eis φ(0) i = e− 2 hφ(x)φ(0)i ,with s, s0 = ±. The required object isZ
Z ∞ Z 1iei~p·~x2πdpd cos θeipx cos θhφ(x)φ(0)i =d̄p 2 = i3Tp(2π) T 0
| −1
{z
}
px
= 2 sinpx

Z ∞sin px2dp2(2π) T Z0px
2 1 1 ∞ sin p̄= idp̄
2πT x 2 −∞
p̄
{z
}
|
= i

=π

i
=
.
2T x

(120)

(I have set the velocity of propagation to 1, and T = U/K.) So
ss00

heisφ(x) eis φ(0) i = e−i 4xT .
And

1
4T x
which does not decay at long distance, and in fact approaches a constant.
hcos φ(x) cos φ(0)i = cos

Adding matter helps to produce a deconfined phase! In particular, the presence of enoughmassless charged fermions can render the monopole operator irrelevant. I recommend thispaper by Tarun Grover for more on this.
[End of Lecture 19]

Next steps:• phase diagram of Z2 gauge theory.
• entanglement in many body quantum systems.

148