Skip to content

Bosonization and Sine-Gordon–Thirring Duality

The previous page ended with a clue: neutral fermion correlators in two dimensions have the same algebraic structure as correlators of exponentials of a free scalar field. This is not an accident, and it is not merely a trick for computing determinants. In two spacetime dimensions, the operator algebra of a massless Dirac fermion can be represented by the operator algebra of a compact scalar. This is bosonization.

Bosonization is one of the places where quantum field theory behaves like a magician who actually explains the trick afterward. A fermion field, which anticommutes and has spin 1/21/2, can be written as an exponential of a bosonic field, provided one treats zero modes, branch cuts, and Klein factors correctly. A four-fermion current-current interaction becomes a change of the scalar radius. A fermion mass term becomes a cosine potential. The massive Thirring model is therefore equivalent, in a precise sense, to the sine-Gordon model.

This page develops the dictionary in the normalization closest to the preceding two-dimensional notes. The goal is not to prove every global subtlety of bosonization on an arbitrary Riemann surface. The goal is to make the local operator identities and the sine-Gordon–Thirring map usable.

Required background. The Schwinger model and gauge-invariant correlators supplies the neutral fermion determinant that becomes a bosonic vertex-operator correlator. Helpful background. Gauge fields in two dimensions supplies the chiral propagators and complex-coordinate conventions.

Two-dimensional bosonization conventions. We use Euclidean complex coordinates

z=x1+ix2,zˉ=x1ix2.z=x^1+i x^2, \qquad \bar z=x^1-i x^2.

After analytic continuation, zz and zˉ\bar z become the two light-cone directions x+x^+ and xx^-. We use two independent chiral bosons and define the relative, sine-Gordon combination by

Φ(z,zˉ)=ϕR(z)ϕL(zˉ),\Phi(z,\bar z)=\phi_R(z)-\phi_L(\bar z),

with two-point functions

ϕR(z)ϕR(w)=logzwa,ϕL(zˉ)ϕL(wˉ)=logzˉwˉa,\langle \phi_R(z)\phi_R(w)\rangle=-\log{z-w\over a}, \qquad \langle \phi_L(\bar z)\phi_L(\bar w)\rangle=-\log{\bar z-\bar w\over a},

and

ϕR(z)ϕL(wˉ)=0.\langle \phi_R(z)\phi_L(\bar w)\rangle=0.

The short-distance length aa is a UV cutoff. With this convention,

Φ(z,zˉ)Φ(w,wˉ)=logzw2a2.\langle \Phi(z,\bar z)\Phi(w,\bar w)\rangle =-\log{|z-w|^2\over a^2}.

The sum Φ~=ϕR+ϕL\widetilde\Phi=\phi_R+\phi_L is the dual combination. Which of Φ\Phi and Φ~\widetilde\Phi is called “the boson” varies across the literature. Keeping the relative combination explicit prevents a sign mismatch between the fermion mass operator and the topological current.

Normal ordering ::: \cdots : is defined with respect to this free scalar. Overall constants in vertex operators depend on the cutoff convention and are absorbed into normalization constants.

Relation to canonical scalar normalization. The chiral convention above corresponds to the nonchiral free-boson action

S=18πd2x(μΦ)2,S={1\over8\pi}\int d^2x\,(\partial_\mu\Phi)^2,

for which :eiϕR::e^{i\phi_R}: has holomorphic weight 1/21/2 and can represent a chiral fermion. If instead one uses the canonical normalization

S=12d2x(μφ)2,S={1\over2}\int d^2x\,(\partial_\mu\varphi)^2,

then

Φ=4πφ.\Phi=\sqrt{4\pi}\,\varphi.

In that convention the free-fermion mass operator corresponds to a sine-Gordon cosine with βSG2=4π\beta_{\mathrm{SG}}^2=4\pi. Many apparent disagreements in bosonization formulas are only this rescaling.

The right-moving free fermion has propagator

ψR(z)ψR(w)=1zw.\langle \psi_R(z)\psi_R^\dagger(w)\rangle={1\over z-w}.

Wick’s theorem gives the nn-particle neutral correlator

ψR(z1)ψR(zn)ψR(w1)ψR(wn)=det1i,jn1ziwj.\left\langle \psi_R(z_1)\cdots\psi_R(z_n) \psi_R^\dagger(w_1)\cdots\psi_R^\dagger(w_n) \right\rangle =\det_{1\le i,j\le n}{1\over z_i-w_j}.

The Cauchy determinant identity rewrites this as

deti,j1ziwj=i<j(zizj)i<j(wjwi)i,j(ziwj).\boxed{ \det_{i,j}{1\over z_i-w_j} ={ \prod_{i<j}(z_i-z_j)\prod_{i<j}(w_j-w_i) \over \prod_{i,j}(z_i-w_j) }. }

The signs in the numerator depend on the ordering convention for the ww operators, but the determinant expression is unambiguous.

Now compare this with a Gaussian scalar. For vertex operators

Vα(z)=:eiαϕR(z):,V_\alpha(z)=:e^{i\alpha\phi_R(z)}:,

Wick’s theorem gives

kVαk(zk)={i<j(zizja)αiαj,kαk=0,0,kαk0,\left\langle \prod_k V_{\alpha_k}(z_k)\right\rangle =\begin{cases} \displaystyle \prod_{i<j}\left({z_i-z_j\over a}\right)^{\alpha_i\alpha_j}, & \sum_k\alpha_k=0,\\[1.3em] 0, & \sum_k\alpha_k\ne 0, \end{cases}

where the second line is the zero-mode selection rule on the plane. Taking charges +1+1 at the ziz_i and charges 1-1 at the wjw_j gives

i=1n:eiϕR(zi):j=1n:eiϕR(wj):=ani<j(zizj)i<j(wiwj)i,j(ziwj)\left\langle \prod_{i=1}^n :e^{i\phi_R(z_i)}: \prod_{j=1}^n :e^{-i\phi_R(w_j)}: \right\rangle =a^n{ \prod_{i<j}(z_i-z_j) \prod_{i<j}(w_i-w_j) \over \prod_{i,j}(z_i-w_j) }

up to a sign fixed by the ordering of fermions. The factor ana^n is required by dimensional analysis and is canceled by the normalization of the 2n2n fermion operators. After that cancellation, this is exactly the same rational function as the fermion determinant.

Cauchy determinant fermion correlator represented as a neutral gas of scalar vertex operators

A neutral right-moving fermion correlator is a Cauchy determinant. The same function is produced by a free chiral scalar with vertex charges +1+1 at fermion insertions and 1-1 at conjugate insertions.

Thus the local identification

ψR(z)FRa:eiϕR(z):,ψR(z)FRa:eiϕR(z):\boxed{ \psi_R(z)\sim {F_R\over\sqrt a}:e^{i\phi_R(z)}:, \qquad \psi_R^\dagger(z)\sim {F_R^\dagger\over\sqrt a}:e^{-i\phi_R(z)}: }

reproduces all neutral right-moving correlators. Similarly,

ψL(zˉ)FLa:eiϕL(zˉ):,ψL(zˉ)FLa:eiϕL(zˉ):.\boxed{ \psi_L(\bar z)\sim {F_L\over\sqrt a}:e^{i\phi_L(\bar z)}:, \qquad \psi_L^\dagger(\bar z)\sim {F_L^\dagger\over\sqrt a}:e^{-i\phi_L(\bar z)}:. }

These prefactors match the preceding convention ψR(z)ψR(w)=1/(zw)\langle\psi_R(z)\psi_R^\dagger(w)\rangle=1/(z-w). If the fermion propagator is normalized instead as 1/[2π(zw)]1/[2\pi(z-w)], each prefactor acquires the familiar additional factor 1/2π1/\sqrt{2\pi}.

The factors FR,FLF_R,F_L are Klein factors. They commute with the bosonic oscillators but anticommute among themselves, ensuring

{ψR,ψL}=0.\{\psi_R,\psi_L\}=0.

For many local correlation functions with equal numbers of each species, the Klein factors only provide the overall fermionic sign. For operator algebra, Hilbert-space construction, boundary conditions, or finite-size systems, they are not optional.

The determinant formula assumes a complex chiral fermion. For a real Majorana fermion, Wick’s theorem instead gives a Pfaffian of the antisymmetric matrix of two-point functions. Two Majorana fields combine into one Dirac field, which is why the charged Cauchy-determinant form is the most direct doorway to bosonization.

Scaling dimensions and the neutrality condition

Section titled “Scaling dimensions and the neutrality condition”

The same Gaussian calculation immediately gives the scaling dimension of a vertex operator. For a chiral operator,

Vα(z)Vα(w)=(azw)α2.\langle V_\alpha(z)V_{-\alpha}(w)\rangle =\left({a\over z-w}\right)^{\alpha^2}.

A holomorphic primary field of weight hh has two-point function proportional to (zw)2h(z-w)^{-2h}, so

h[Vα]=α22.\boxed{h[V_\alpha]={\alpha^2\over2}.}

The fermion operator corresponds to α=1\alpha=1 and therefore has h=1/2h=1/2, as it should.

For the nonchiral vertex operator

Vα,αˉ(z,zˉ)=:eiαϕR(z)+iαˉϕL(zˉ):,V_{\alpha,\bar\alpha}(z,\bar z) =:e^{i\alpha\phi_R(z)+i\bar\alpha\phi_L(\bar z)}:,

we have

h=α22,hˉ=αˉ22,h={\alpha^2\over2}, \qquad \bar h={\bar\alpha^2\over2},

so the scaling dimension and Euclidean spin are

Δ=h+hˉ=α2+αˉ22,s=hhˉ=α2αˉ22.\Delta=h+\bar h={\alpha^2+\bar\alpha^2\over2}, \qquad s=h-\bar h={\alpha^2-\bar\alpha^2\over2}.

A right-moving spinor has spin +1/2+1/2, so any bosonic representation of an interacting right-moving fermion must satisfy

α2αˉ2=1.\boxed{\alpha^2-\bar\alpha^2=1.}

This simple equation is the seed of anomalous fermion dimensions in the massless Thirring model.

The neutrality condition is just as important. The nonchiral correlator

k:eiαkΦ(xk):\left\langle\prod_k :e^{i\alpha_k\Phi(x_k)}:\right\rangle

contains the zero mode of Φ\Phi. Since the free boson action depends only on derivatives, the constant mode is integrated over. The integral

dΦ0exp(iΦ0kαk)\int d\Phi_0\,\exp\left(i\Phi_0\sum_k\alpha_k\right)

vanishes unless

kαk=0.\boxed{\sum_k\alpha_k=0.}

Equivalently, if one regulates the theory in a box of size RR, nonneutral correlators carry powers of RR and disappear in the infinite-volume limit after normalizing the vacuum. This is the infrared counterpart of charge conservation.

The operator product expansion of vertex operators follows from separating the contraction between nearby points from the normal-ordered product:

:eiαϕR(z)::eiβϕR(w):=exp[αβϕR(z)ϕR(w)]:eiαϕR(z)+iβϕR(w):.:e^{i\alpha\phi_R(z)}::e^{i\beta\phi_R(w)}: =\exp\left[-\alpha\beta\langle\phi_R(z)\phi_R(w)\rangle\right] :e^{i\alpha\phi_R(z)+i\beta\phi_R(w)}:.

Using

ϕR(z)ϕR(w)=logzwa,\langle\phi_R(z)\phi_R(w)\rangle=-\log{z-w\over a},

and expanding the slowly varying field ϕR(z)\phi_R(z) around ww, we get

:eiαϕR(z)::eiβϕR(w):(zwa)αβ:ei(α+β)ϕR(w):[1+iα(zw)ϕR(w)+].\boxed{ :e^{i\alpha\phi_R(z)}::e^{i\beta\phi_R(w)}: \sim \left({z-w\over a}\right)^{\alpha\beta} :e^{i(\alpha+\beta)\phi_R(w)}: \left[1+i\alpha(z-w)\partial\phi_R(w)+\cdots\right]. }

For α=1\alpha=1 and β=1\beta=-1,

:eiϕR(z)::eiϕR(w):azw[1+i(zw)ϕR(w)+],:e^{i\phi_R(z)}::e^{-i\phi_R(w)}: \sim {a\over z-w}\left[1+i(z-w)\partial\phi_R(w)+\cdots\right],

which reproduces the fermion short-distance singularity. For α=β=1\alpha=\beta=1,

:eiϕR(z)::eiϕR(w):zwa:e2iϕR(w):+.:e^{i\phi_R(z)}::e^{i\phi_R(w)}: \sim {z-w\over a}:e^{2i\phi_R(w)}: +\cdots.

The zero as zwz\to w is the bosonic representation of Pauli exclusion for two identical chiral fermions. The anticommutation sign is encoded by the branch choice and Klein factors; the vanishing of the local product is already visible in the OPE.

OPE fusion rule for two scalar vertex operators

When two vertex operators approach each other, their charges add. The singular or vanishing prefactor is determined by the product of charges, (zw)αβ(z-w)^{\alpha\beta}.

This OPE is the local engine behind bosonization. The global theory contains more data: compactification radius, allowed charges, spin structures, and Klein factors. But the short-distance algebra already knows why exponentials of a free scalar can behave like fermions.

Current dictionary and charge normalization

Section titled “Current dictionary and charge normalization”

The vector current is the derivative of the relative boson. In the chiral normalization used above, and up to the overall sign fixed by the orientation convention for ϵμν\epsilon^{\mu\nu},

jμ=12πϵμννΦ.j^\mu={1\over2\pi}\epsilon^{\mu\nu}\partial_\nu\Phi.

Equivalently, in the canonical free-fermion normalization Φ=4πφ\Phi=\sqrt{4\pi}\,\varphi,

jμ=1πϵμννφ.j^\mu={1\over\sqrt\pi}\epsilon^{\mu\nu}\partial_\nu\varphi.

The charge is therefore a winding number:

Q=dxj0=12π[Φ(+)Φ()]=1π[φ(+)φ()].Q=\int dx\,j^0={1\over2\pi}\left[\Phi(+\infty)-\Phi(-\infty)\right] ={1\over\sqrt\pi}\left[\varphi(+\infty)-\varphi(-\infty)\right].

At nonzero Thirring coupling the canonically normalized sine-Gordon field has a coupling-dependent radius, and the corresponding dictionary becomes jμ=(βSG/2π)ϵμννφj^\mu=(\beta_{\mathrm{SG}}/2\pi)\epsilon^{\mu\nu}\partial_\nu\varphi. The free point βSG2=4π\beta_{\mathrm{SG}}^2=4\pi reproduces the coefficient 1/π1/\sqrt\pi.

This formula is the bridge between the local vertex-operator dictionary and the soliton interpretation below. A fermion is local in fermion variables but becomes a kink in the bosonic field.

The massless Thirring model as a free boson with a shifted radius

Section titled “The massless Thirring model as a free boson with a shifted radius”

The massless Thirring model is a Dirac fermion in two dimensions with a current-current interaction:

LTh,0=ψˉiγμμψgTh2jμjμ,jμ=ψˉγμψ.\mathcal L_{\mathrm{Th},0} =\bar\psi i\gamma^\mu\partial_\mu\psi -{g_{\mathrm{Th}}\over2}j_\mu j^\mu, \qquad j^\mu=\bar\psi\gamma^\mu\psi.

In two dimensions this interaction is special. It is classically marginal: since [ψ]=1/2[\psi]=1/2, the current has dimension 11, so jμjμj_\mu j^\mu has dimension 22, the same as the Lagrangian density. Quantum mechanically, the model remains exactly solvable. The interaction does not create an arbitrary complicated non-Gaussian bosonic theory; it changes the normalization, or radius, of the boson.

A convenient way to encode this is to represent the interacting right- and left-moving fermions by

ψRFR:eiaϕR+ibϕL:,ψLFL:eibϕR+iaϕL:,\psi_R\sim F_R:e^{i a\phi_R+i b\phi_L}:, \qquad \psi_L\sim F_L:e^{i b\phi_R+i a\phi_L}:,

with

a2b2=1.a^2-b^2=1.

Then

ψR(z,zˉ)ψR(0,0)1za2zˉb2=1z1z2b2.\langle\psi_R(z,\bar z)\psi_R^\dagger(0,0)\rangle \propto {1\over z^{a^2}\bar z^{b^2}} ={1\over z}\,{1\over |z|^{2b^2}}.

The first factor gives the spinor transformation law; the second is an anomalous power. In Lorentzian notation this has the form

TψR(x)ψR(0)1x+f(x+x),\langle T\psi_R(x)\psi_R^\dagger(0)\rangle \propto {1\over x^+}\,f(x^+x^-),

where the function ff is a power law in the conformal massless theory. Thus the Thirring interaction changes the scaling dimension of the fermion without changing its spin.

There are several equivalent normalizations in the literature. In the canonical Coleman normalization one writes a scalar φ\varphi with kinetic term

12(μφ)2,{1\over2}(\partial_\mu\varphi)^2,

and the Thirring coupling is related to the sine-Gordon coupling by

4πβSG2=1+gThπ.\boxed{ {4\pi\over \beta_{\mathrm{SG}}^2}=1+{g_{\mathrm{Th}}\over\pi}. }

In the chiral normalization above, the corresponding parameter is

βSG2=4πλ2,λ=ab.\beta_{\mathrm{SG}}^2=4\pi\lambda^2, \qquad \lambda=a-b.

At gTh=0g_{\mathrm{Th}}=0, one has b=0b=0, a=1a=1, and βSG2=4π\beta_{\mathrm{SG}}^2=4\pi, the free-fermion point.

The massive Thirring model adds

ΔL=mψˉψ.\Delta\mathcal L=m\bar\psi\psi.

In chiral components,

ψˉψ=ψLψR+ψRψL.\bar\psi\psi=\psi_L^\dagger\psi_R+\psi_R^\dagger\psi_L.

Using the interacting-fermion representation,

ψLψR:ei(ab)(ϕRϕL):.\psi_L^\dagger\psi_R \sim :e^{i(a-b)(\phi_R-\phi_L)}:.

Using the relative field already introduced above, define only the remaining exponent

λ=ab.\lambda=a-b.

Then the mass operator becomes

ψˉψ  Cmcos(λΦ),\boxed{ \bar\psi\psi\ \longleftrightarrow\ C_m\cos(\lambda\Phi), }

where CmC_m is a cutoff-dependent normalization constant. Therefore the massive Thirring model maps to a sine-Gordon theory,

SSG,E=d2x[12(μφ)2+μβSG2(1cos(βSGφ))]\boxed{ S_{\mathrm{SG},E} =\int d^2x\left[ {1\over2}(\partial_\mu\varphi)^2 +{\mu\over\beta_{\mathrm{SG}}^2}\left(1-\cos(\beta_{\mathrm{SG}}\varphi)\right) \right] }

in canonical normalization, or equivalently to

SE=18πd2x(μΦ)2+μ~d2x[1cos(λΦ)]S_E={1\over8\pi}\int d^2x\,(\partial_\mu\Phi)^2 +\widetilde\mu\int d^2x\,[1-\cos(\lambda\Phi)]

in the chiral normalization of the discussion, with βSG=4πλ\beta_{\mathrm{SG}}=\sqrt{4\pi}\lambda. The precise relation between mm, μ\mu, μ~\widetilde\mu, and the cutoff is nonuniversal; the relation between the operator algebras and the dimensionless coupling βSG\beta_{\mathrm{SG}} is universal once the renormalization convention is fixed.

Massless Thirring model becomes a free boson with shifted radius, and the fermion mass becomes a sine-Gordon cosine

The current-current interaction of the massless Thirring model changes the compact-boson radius, represented by KK in the diagram and by the exponents a,ba,b in the operator dictionary. The fermion mass operator then becomes the sine-Gordon cosine.

This is the massive Thirring/sine-Gordon equivalence. It is strong–weak in an important sense. From

4πβSG2=1+gThπ,{4\pi\over \beta_{\mathrm{SG}}^2}=1+{g_{\mathrm{Th}}\over\pi},

large positive Thirring coupling corresponds to small sine-Gordon coupling. In that regime the sine-Gordon semiclassical soliton is a useful description of the fermionic excitation.

The mass perturbation also provides a direct consistency check. Because

:eiλΦ(x)::eiλΦ(0):1x2λ2,\left\langle :e^{i\lambda\Phi(x)}: :e^{-i\lambda\Phi(0)}: \right\rangle \propto {1\over |x|^{2\lambda^2}},

the cosine has scaling dimension

Δcos=λ2=βSG24π=11+gTh/π.\boxed{ \Delta_{\cos}=\lambda^2 ={\beta_{\mathrm{SG}}^2\over4\pi} ={1\over1+g_{\mathrm{Th}}/\pi}. }

At the free-fermion point, Δcos=1\Delta_{\cos}=1, exactly the dimension of ψˉψ\bar\psi\psi. The cosine is relevant when βSG2<8π\beta_{\mathrm{SG}}^2<8\pi, marginal at 8π8\pi, and irrelevant on the other side of that boundary. This is the renormalization-group content hidden inside the apparently simple replacement of a fermion mass by a periodic potential.

The sine-Gordon potential is periodic. Its classical vacua are

βSGφ=2πn,nZ.\beta_{\mathrm{SG}}\varphi=2\pi n, \qquad n\in\mathbb Z.

A finite-energy static configuration must approach vacua at spatial infinity:

φ(x)=2πnβSG,φ(x+)=2πn+βSG.\varphi(x\to-\infty)={2\pi n_-\over\beta_{\mathrm{SG}}}, \qquad \varphi(x\to+\infty)={2\pi n_+\over\beta_{\mathrm{SG}}}.

The integer

Qtop=n+n=βSG2π[φ(+)φ()]Q_{\mathrm{top}}=n_+-n_- ={\beta_{\mathrm{SG}}\over2\pi}\left[\varphi(+\infty)-\varphi(-\infty)\right]

is the soliton number. The one-soliton solution for the potential

V(φ)=mSG2βSG2(1cosβSGφ)V(\varphi)={m_{\mathrm{SG}}^2\over\beta_{\mathrm{SG}}^2} \left(1-\cos\beta_{\mathrm{SG}}\varphi\right)

is

φsol(xX)=4βSGarctanemSG(xX).\boxed{ \varphi_{\mathrm{sol}}(x-X) ={4\over\beta_{\mathrm{SG}}}\arctan e^{m_{\mathrm{SG}}(x-X)}. }

It interpolates from 00 to 2π/βSG2\pi/\beta_{\mathrm{SG}}. For a static finite-energy solution, the first integral of the field equation is

12(dφdx)2=V(φ),{1\over2}\left({d\varphi\over dx}\right)^2 =V(\varphi),

where the integration constant vanishes because both terms go to zero in a vacuum. The energy can therefore be evaluated without inserting the detailed profile:

Msol=02π/βSGdφ2V(φ).M_{\mathrm{sol}} =\int_0^{2\pi/\beta_{\mathrm{SG}}} d\varphi\,\sqrt{2V(\varphi)}.

Evaluating the elementary integral gives

Msol=8mSGβSG2.\boxed{ M_{\mathrm{sol}}={8m_{\mathrm{SG}}\over\beta_{\mathrm{SG}}^2}. }

The interacting bosonized current,

jμ=βSG2πϵμννφ,j^\mu={\beta_{\mathrm{SG}}\over2\pi} \epsilon^{\mu\nu}\partial_\nu\varphi,

identifies fermion number with topological charge:

QF=dxj0=βSG2π[φ(+)φ()]=Qtop.Q_F=\int dx\,j^0 ={\beta_{\mathrm{SG}}\over2\pi}\left[\varphi(+\infty)-\varphi(-\infty)\right] =Q_{\mathrm{top}}.

Thus the elementary fermion of the Thirring description is represented, in the sine-Gordon description, by a soliton. The antifermion is the antisoliton.

Topological charge alone does not prove fermionic statistics. The nonlocal soliton operator and its branch prescription supply that extra information; in the full Coleman–Mandelstam dictionary, those operators obey the massive Thirring anticommutation relations.

Sine-Gordon periodic vacua, soliton interpolation, and tilted potential

The sine-Gordon field has infinitely many equivalent vacua. A soliton interpolates between neighboring vacua and carries one unit of fermion number. A background electric field or theta-like source tilts the periodic potential, making adjacent vacua energetically inequivalent and foreshadowing pair creation.

This identification is one of the sharpest lessons of two-dimensional field theory: what is elementary in one set of variables may be extended and topological in another. The word “duality” should be taken literally here. The two descriptions organize the same Hilbert space in different coordinates.

Tilted vacua and the electric-field interpretation

Section titled “Tilted vacua and the electric-field interpretation”

The sine-Gordon potential also gives a compact way to visualize electric backgrounds in two-dimensional gauge theory. A background electric field, or equivalently a theta-like source, can add a term linear in the boson,

V(φ)mSG2βSG2(1cosβSGφ)hEφ.V(\varphi)\to {m_{\mathrm{SG}}^2\over\beta_{\mathrm{SG}}^2} \left(1-\cos\beta_{\mathrm{SG}}\varphi\right) -h_E\varphi.

The coefficient hEh_E is proportional to the applied electric field; its exact normalization depends on the gauge-field and bosonization conventions.

The shift

φφ+2πβSG\varphi\to\varphi+{2\pi\over\beta_{\mathrm{SG}}}

no longer relates exactly degenerate vacua. Neighboring minima are tilted relative to each other. In the fermionic language, producing a soliton–antisoliton pair changes the electric flux between them. If the electric field is strong enough, pair production is energetically favored.

This is not yet the full confinement/screening story of QED₂, but it is already the right picture. The bosonic field keeps track of charge through its winding. A region between two solitons can sit in a different branch of the electric field. The next page uses this idea to discuss linear potentials, screening by massless fermions, and confinement in one spatial dimension.

Bosonization begins from a concrete identity: the Cauchy determinant of free chiral fermion correlators equals a neutral correlator of scalar vertex operators. With the chiral scalar normalized by

ϕR(z)ϕR(w)=log(zw),\langle\phi_R(z)\phi_R(w)\rangle=-\log(z-w),

the operator :eiϕR::e^{i\phi_R}: has weight 1/21/2 and reproduces a right-moving fermion. The zero mode enforces charge neutrality, while the vertex-operator OPE encodes both the fermion pole and the short-distance exclusion of identical chiral fermions.

The massless Thirring interaction changes the boson radius and gives the fermion an anomalous dimension while preserving its spin. The fermion mass term becomes a cosine perturbation. This turns the massive Thirring model into the sine-Gordon model, with the Coleman relation

4πβSG2=1+gThπ{4\pi\over\beta_{\mathrm{SG}}^2}=1+{g_{\mathrm{Th}}\over\pi}

in a standard normalization. In the sine-Gordon variables, fermion number is the topological winding number of the bosonic field, and the elementary fermion is a soliton.

Confusing the two beta symbols. Here βSG\beta_{\mathrm{SG}} is the sine-Gordon coupling, not an RG beta function. The chiral exponents were called a,ba,b to avoid the worst collision.

Dropping Klein factors from the full operator algebra. They often disappear from neutral local correlators, but they are needed so that right- and left-moving fermions anticommute.

Identifying the sum and relative bosons. With the fermion phases used here, the mass operator and vector current involve Φ=ϕRϕL\Phi=\phi_R-\phi_L; the sum Φ~=ϕR+ϕL\widetilde\Phi=\phi_R+\phi_L is the dual field. Changing the sign convention for the left-moving vertex exchanges these names, but the dictionary must change everywhere at once.

Treating the mass coefficient as universal. The relation between the Thirring mass parameter mm and the sine-Gordon cosine coefficient μ\mu depends on the UV normalization of the composite operator ψˉψ\bar\psi\psi. The relation between scaling dimensions and the dimensionless coupling is the universal part.

Reading bosonization as a classical field identity. The massive Thirring/sine-Gordon relation is a quantum operator equivalence. A classical Dirac field does not literally equal a classical exponential of a scalar.

Using

ϕR(z)ϕR(w)=logzwa,\langle\phi_R(z)\phi_R(w)\rangle=-\log{z-w\over a},

show that

:eiαϕR(z)::eiαϕR(w):=(azw)α2.\left\langle :e^{i\alpha\phi_R(z)}::e^{-i\alpha\phi_R(w)}:\right\rangle =\left({a\over z-w}\right)^{\alpha^2}.

What is the conformal weight of :eiαϕR::e^{i\alpha\phi_R}:?

Solution

For normal-ordered exponentials of a Gaussian field,

:eA::eB:=eAB.\langle :e^A::e^B:\rangle=e^{\langle AB\rangle}.

Here

A=iαϕR(z),B=iαϕR(w),A=i\alpha\phi_R(z), \qquad B=-i\alpha\phi_R(w),

so

AB=α2ϕR(z)ϕR(w)=α2logzwa.\langle AB\rangle =\alpha^2\langle\phi_R(z)\phi_R(w)\rangle =-\alpha^2\log{z-w\over a}.

Therefore

:eiαϕR(z)::eiαϕR(w):=exp[α2logzwa]=(azw)α2.\left\langle :e^{i\alpha\phi_R(z)}::e^{-i\alpha\phi_R(w)}:\right\rangle =\exp\left[-\alpha^2\log{z-w\over a}\right] =\left({a\over z-w}\right)^{\alpha^2}.

A holomorphic primary of weight hh has two-point function proportional to (zw)2h(z-w)^{-2h}, so

h=α22.h={\alpha^2\over2}.

Exercise 2: The two-by-two Cauchy identity

Section titled “Exercise 2: The two-by-two Cauchy identity”

Prove the n=2n=2 Cauchy identity

det(1z1w11z1w21z2w11z2w2)=(z1z2)(w2w1)(z1w1)(z1w2)(z2w1)(z2w2).\det \begin{pmatrix} {1\over z_1-w_1} & {1\over z_1-w_2}\\[0.5em] {1\over z_2-w_1} & {1\over z_2-w_2} \end{pmatrix} ={ (z_1-z_2)(w_2-w_1) \over (z_1-w_1)(z_1-w_2)(z_2-w_1)(z_2-w_2)}.
Solution

Compute the determinant directly:

D=1(z1w1)(z2w2)1(z1w2)(z2w1).D={1\over(z_1-w_1)(z_2-w_2)} -{1\over(z_1-w_2)(z_2-w_1)}.

Putting the two terms over a common denominator gives

D=(z1w2)(z2w1)(z1w1)(z2w2)(z1w1)(z1w2)(z2w1)(z2w2).D={ (z_1-w_2)(z_2-w_1)-(z_1-w_1)(z_2-w_2) \over (z_1-w_1)(z_1-w_2)(z_2-w_1)(z_2-w_2)}.

The numerator is

(z1z2z1w1w2z2+w2w1)(z1z2z1w2w1z2+w1w2)(z_1z_2-z_1w_1-w_2z_2+w_2w_1) -(z_1z_2-z_1w_2-w_1z_2+w_1w_2)

or

z1(w2w1)+z2(w1w2)=(z1z2)(w2w1).z_1(w_2-w_1)+z_2(w_1-w_2)=(z_1-z_2)(w_2-w_1).

This proves the formula.

Let

Va,b(z,zˉ)=:eiaϕR(z)+ibϕL(zˉ):.V_{a,b}(z,\bar z)=:e^{ia\phi_R(z)+ib\phi_L(\bar z)}:.

Show that its Euclidean spin is

s=a2b22.s={a^2-b^2\over2}.

Use this to derive the condition for Va,bV_{a,b} to represent a right-moving spinor.

Solution

The right-moving and left-moving weights are

h=a22,hˉ=b22.h={a^2\over2}, \qquad \bar h={b^2\over2}.

The Euclidean spin is the difference

s=hhˉ=a2b22.s=h-\bar h={a^2-b^2\over2}.

A right-moving spinor has spin +1/2+1/2, so

a2b22=12.{a^2-b^2\over2}={1\over2}.

Therefore

a2b2=1.a^2-b^2=1.

Using

ψRFR:eiaϕR+ibϕL:,ψLFL:eibϕR+iaϕL:,\psi_R\sim F_R:e^{ia\phi_R+ib\phi_L}:, \qquad \psi_L\sim F_L:e^{ib\phi_R+ia\phi_L}:,

show that the mass operator ψLψR+ψRψL\psi_L^\dagger\psi_R+\psi_R^\dagger\psi_L is proportional to cos[(ab)(ϕRϕL)]\cos[(a-b)(\phi_R-\phi_L)].

Solution

Ignoring Klein factors and normalization constants, the first term is

ψLψR:eibϕRiaϕL::eiaϕR+ibϕL:.\psi_L^\dagger\psi_R \sim :e^{-ib\phi_R-ia\phi_L}::e^{ia\phi_R+ib\phi_L}:.

At coincident points the product must be renormalized, but the exponent is simply

i(ab)ϕR+i(ba)ϕL=i(ab)(ϕRϕL).i(a-b)\phi_R+i(b-a)\phi_L =i(a-b)(\phi_R-\phi_L).

Thus

ψLψR:ei(ab)(ϕRϕL):.\psi_L^\dagger\psi_R\sim :e^{i(a-b)(\phi_R-\phi_L)}:.

Similarly,

ψRψL:ei(ab)(ϕRϕL):.\psi_R^\dagger\psi_L\sim :e^{-i(a-b)(\phi_R-\phi_L)}:.

Adding them gives

ψLψR+ψRψLcos[(ab)(ϕRϕL)].\psi_L^\dagger\psi_R+\psi_R^\dagger\psi_L \propto \cos[(a-b)(\phi_R-\phi_L)].

Exercise 5: The classical sine-Gordon soliton

Section titled “Exercise 5: The classical sine-Gordon soliton”

For the sine-Gordon energy functional

E=dx[12(dφdx)2+m2β2(1cosβφ)],E=\int_{-\infty}^{\infty}dx\left[ {1\over2}\left({d\varphi\over dx}\right)^2 +{m^2\over\beta^2}(1-\cos\beta\varphi) \right],

show that the one-soliton solution

φsol(x)=4βarctanemx\varphi_{\mathrm{sol}}(x)={4\over\beta}\arctan e^{mx}

interpolates from 00 to 2π/β2\pi/\beta. Then verify that it satisfies the first-order equation

dφdx=2mβsinβφ2.{d\varphi\over dx}={2m\over\beta}\sin{\beta\varphi\over2}.
Solution

As xx\to-\infty, emx0e^{mx}\to0, so

φsol()=0.\varphi_{\mathrm{sol}}(-\infty)=0.

As x+x\to+\infty, arctanemxπ/2\arctan e^{mx}\to\pi/2, so

φsol(+)=4βπ2=2πβ.\varphi_{\mathrm{sol}}(+\infty)={4\over\beta}{\pi\over2}={2\pi\over\beta}.

Now differentiate. Let u=emxu=e^{mx}. Then

dφdx=4β11+u2mu=4mβu1+u2.{d\varphi\over dx} ={4\over\beta}{1\over1+u^2}mu ={4m\over\beta}{u\over1+u^2}.

Since

βφ2=2arctanu,{\beta\varphi\over2}=2\arctan u,

we use

sin(2arctanu)=2u1+u2.\sin(2\arctan u)={2u\over1+u^2}.

Therefore

2mβsinβφ2=2mβ2u1+u2=4mβu1+u2=dφdx.{2m\over\beta}\sin{\beta\varphi\over2} ={2m\over\beta}{2u\over1+u^2} ={4m\over\beta}{u\over1+u^2} ={d\varphi\over dx}.

This proves the first-order equation.

  • S. Coleman, “Quantum sine-Gordon equation as the massive Thirring model,” Physical Review D 11 (1975), 2088–2097.
  • S. Mandelstam, “Soliton operators for the quantized sine-Gordon equation,” Physical Review D 11 (1975), 3026–3030.
  • D. C. Mattis and E. H. Lieb, “Exact solution of a many-fermion system and its associated boson field,” Journal of Mathematical Physics 6 (1965), 304–312.
  • S. Coleman, Aspects of Symmetry: Selected Erice Lectures (Cambridge University Press, 1985).
  • T. Giamarchi, Quantum Physics in One Dimension (Oxford University Press, 2004).
  • A. O. Gogolin, A. A. Nersesyan, and A. M. Tsvelik, Bosonization and Strongly Correlated Systems (Cambridge University Press, 1998).
  • A. M. Polyakov, Gauge Fields and Strings (Harwood Academic Publishers, 1987).
  • A. Zee, Quantum Field Theory in a Nutshell, 2nd ed. (Princeton University Press, 2010).