format: text

text:
# 2D MODEL FIELD THEORIES AT FINITE TEMPERATURE AND DENSITY **footnote *Contribution to the Festschrift in honor of Boris Ioffe, edited by M. Shifman.

## 1 Introduction

With the shift in emphasis in field theory from perturbative to non-perturbative phenomena during the last decades, exactly soluble models have acquired a quite respectable status. The limited number of analytical tools available at strong coupling, so painfully felt in quantum chromodynamics (QCD), makes it tempting to look for other sources of inspiration. Except for the celebrated case of supersymmetric $`N=2`$ Yang-Mills theory $`^\mathrm{?}`$, “soluble model” in quantum field theory still means “1+1 dimensional model”. This makes it clear at once that we should not expect such toy models to be in any sense realistic; they can teach us only little (if anything) of direct relevance to the real world. Nevertheless, we find it quite appropriate to include them in a book devoted to analytic approaches to QCD. One can perhaps picture the relationship between QCD and the soluble toy models as follows: If QCD is the sun, then the 1+1 dimensional models are like planets which are gravitating around it, shining in its light. Their distance to the sun may vary and each of them develops a life of its own, although powered by the sun. Only then can one understand why a theory like the Schwinger model $`^\mathrm{?}`$ describing massless electrons which only move along a line, have a linear Coulomb potential and form non-interacting, massive bosons can accumulate close to 1000 citations in the hep-archive up to now. This striking success is shared by the Gross-Neveu (GN) model $`^\mathrm{?}`$ (massless fermions with a point interaction) and the ’t Hooft model $`^\mathrm{?}`$ (QCD2 in the large $`N`$ limit). Whether one likes these models or not — they have become an integral part of the intellectual endeavor to understand the strong interactions and their traces pervade the QCD literature.

There is another aspect which we should like to mention. A soluble model in quantum mechanics is typically solved once and for all, then put into the textbooks where it ceases to play any scientific role. By contrast, a soluble field theory is a laboratory which can be used over and over again for new types of questions — this is the reason why it is so much richer. When browsing through the numerous citations of models like the GN model or the ’t Hooft model, one finds the whole history of QCD mapped out – ranging from questions of asymptotic freedom, quark-hadron duality, spectroscopy, scattering, over issues of confinement, covariance, vacuum structure, chiral symmetry breaking, up to sum rules, supersymmetric extensions, “stringy” aspects, heavy quarks, and finally hot and dense baryonic matter, our main topic here. Also some technical developments like light-cone quantization would be unthinkable without 2d models as firm testing ground. Whether in the end these models will help to solve the real hard problems in QCD is still an open question; they certainly help to prepare us intellectually and widen our horizon about the fascinating world of non-perturbative physics and non-Abelian gauge theories.

As will be seen later on, one of the main themes of this article is chiral symmetry and its spontaneous breakdown (SSB). We are aware that studying SSB in low dimensions is like dancing on thin ice — there is the constant danger of falling into the trap of some no-go theorem and making a fool of oneself. Without claiming mathematical rigor, we believe that the large $`N`$ limit provides a reasonable protection against these perils. Thereby, it significantly enlarges the spectrum of questions which can be addressed within lower dimensional soluble models.

In the spirit of what has been said so far, we have attempted to write this article in a rather self-contained and hopefully pedagogical way. The literature about two-dimensional field theory is not always easy to access for the non-expert, mainly because a number of different, highly specialized techniques are employed by different sub-communities (e.g. bosonization, light-cone quantization, functional methods, conformal and topological field theory). We have refrained from using genuine 2d-technology as much as possible and instead based our presentation upon a simple yet powerful technique familiar to every physicist, the Hartree-Fock (HF) approximation.$`^{\mathrm{?},\mathrm{?}}`$ For the large $`N`$ models which we will consider here, one can actually go quite far with a relativistic extension of this basic tool from many-body theory. This should also make the article more readable for those who want to learn field theory while working on QCD related problems (e.g. students or researchers with a nuclear physics background like ourselves).

Finally we should apologize for the perhaps not always balanced choice of citations: It is strongly biased towards work which fits into our approach and may not do justice to other important activities in the field. Our article is not intended as a comprehensive review article of the usual type, but rather an attempt to convey some of the intellectual pleasure of working on QCD motivated, analytically soluble problems, using as example of topical interest hot and dense matter in two-dimensional large $`N`$ models. This is presumably more adequate for a contribution to a Festschrift like the present one in honor of Boris Ioffe.

## 2 Basics about 1+1 dimensions

The transition from the real world to 1+1 dimensions has drastic consequences. These consequences are responsible for tremendous technical simplifications of quantum field theory, but also for many unrealistic features and even pathologies which have to be kept in mind. In the present section, we have collected the main kinematical and geometrical facts which are independent of the specific dynamics, as well as some characteristic aspects of two-dimensional gauge theories. The selection of topics is guided by the applications which we have in mind in later sections. Therefore, we leave out a large number of important topics which are closely associated with low dimensions but less relevant for us.$`^{\mathrm{?},\mathrm{?}}`$ These include boson-fermion mappings, existence of anyons with arbitrary statistics, conformal invariance, sigma models, integrability questions as well as relations to matrix models, string theories, and gravity. We will also make no use of the technique of light-cone quantization which has become a powerful tool for dealing with two-dimensional model field theories.$`^\mathrm{?}`$ It is not clear to us whether it is of any help in the study of field theory at finite temperature and chemical potential.

### 2.1 Kinematics

We use the metric and antisymmetric tensors

$$g^{\mu \nu }=\left(\begin{array}{cc}\hfill 1& \hfill 0\\ \hfill 0& \hfill 1\end{array}\right),ϵ^{\mu \nu }=\left(\begin{array}{cc}\hfill 0& \hfill 1\\ \hfill 1& \hfill 0\end{array}\right)$$
(1)
in 1+1 dimensions. Clearly, the Lorentz group SO(1,1) consists only of boosts — there are no rotations (and hence there is no spin). The inhomogeneous Lorentz group or Poincaré group ISO(1,1) has three generators (boost $`K`$, momentum $`P`$, Hamiltonian $`H`$) which satisfy the Poincaré algebra

$$[H,P]=0,[K,P]=\mathrm{i}H,[K,H]=\mathrm{i}P.$$
(2)
States are characterized by the Casimir operator

$$P_\mu P^\mu =m^2,$$
(3)
the dispersion relation has the standard form $`E=\sqrt{p^2+m^2}`$. A left-over type of helicity can be defined via Lorentz-boosts

$$\left(\begin{array}{c}x^0\\ x^1\end{array}\right)^{}=\left(\begin{array}{cc}\mathrm{cosh}\xi & \mathrm{sinh}\xi \\ \mathrm{sinh}\xi & \mathrm{cosh}\xi \end{array}\right)\left(\begin{array}{c}x^0\\ x^1\end{array}\right)$$
(4)
with rapidity

$$\xi =\mathrm{artanh}v.$$
(5)
States are said to have helicity $`s`$ if they transform according to

$$|s^{}=\mathrm{e}^{\xi s}|s.$$
(6)
Lorentz scalars, vectors and tensors can be constructed in the usual way. If $`V^\mu `$ is a two-vector, its light-cone components transform irreducibly according to helicity $`\pm 1`$,

$$(V^0\pm V^1)^{}=\mathrm{e}^{\pm \xi }(V^0\pm V^1).$$
(7)
The trace of a tensor $`T_\mu ^\mu `$ is Lorentz scalar, as is an antisymmetric tensor (e.g. the electric field strength). A symmetric, traceless tensor has helicity 2,

$$(T^{00}\pm T^{01})^{}=\mathrm{e}^{\pm 2\xi }(T^{00}\pm T^{01}).$$
(8)
One peculiar feature of massless states in two dimensions is the following: Any number of massless non-interacting particles moving in the same direction with momenta $`p_i`$ behave kinematically like a single massless particle of momentum $`_ip_i`$. This is partly responsible for the proliferation of massless modes in certain models.

### 2.2 Dirac fermions

We start our survey of various aspects of two-dimensional field theory with Dirac fermions which will be the main actors in the following. Canonically, fermions are distinguished from bosons by equal-time anticommutation relations,

$$\{\psi _i(x),\psi _j^{}(y)\}=\delta _{ij}\delta (xy).$$
(9)
Since there is no spin, a Dirac fermion has only two components (particle/antiparticle degrees of freedom). The free Dirac equation reads

$$(\mathrm{i}\gamma ^\mu _\mu m)\psi =0$$
(10)
where the $`\gamma `$-matrices are 2$`\times `$2 matrices satisfying

$$\{\gamma ^\mu ,\gamma ^\nu \}=2g^{\mu \nu }.$$
(11)
Unless stated otherwise, they will be chosen as

$$\gamma ^0=\left(\begin{array}{cc}\hfill 0& \hfill 1\\ \hfill 1& \hfill 0\end{array}\right),\gamma ^1=\left(\begin{array}{cc}\hfill 0& \hfill 1\\ \hfill 1& \hfill 0\end{array}\right).$$
(12)
The free Dirac Hamiltonian (in first quantization) is

$$h=\gamma ^5p+\gamma ^0m$$
(13)
with

$$\gamma ^5=\gamma ^0\gamma ^1=\left(\begin{array}{cc}\hfill 1& \hfill 0\\ \hfill 0& \hfill 1\end{array}\right).$$
(14)
This representation is well suited for the chiral limit $`m0`$. Free spinors (eigenvalues $`\pm E_p`$ with $`E_p=\sqrt{m^2+p^2}`$) are

$`u(p)`$ $`=`$ $`\left(\begin{array}{c}\hfill \mathrm{cos}\theta _p/2\\ \hfill \mathrm{sin}\theta _p/2\end{array}\right)={\displaystyle \frac{1}{\sqrt{2E_p}}}\left(\begin{array}{c}\hfill \sqrt{E_p+p}\\ \hfill \sqrt{E_pp}\end{array}\right),`$ (19)
$`v(p)`$ $`=`$ $`\left(\begin{array}{c}\hfill \mathrm{sin}\theta _p/2\\ \hfill \mathrm{cos}\theta _p/2\end{array}\right)={\displaystyle \frac{1}{\sqrt{2E_p}}}\left(\begin{array}{c}\hfill \sqrt{E_pp}\\ \hfill \sqrt{E_p+p}\end{array}\right).`$ (24)
For later convenience, we have introduced the free Bogoliubov angle as

$$\theta _p=\mathrm{arccot}\left(\frac{p}{m}\right)[0,\pi ].$$
(25)
It appears if one diagonalizes $`h`$ of Eq. (13) by a canonical transformation with the unitary operator

$$U_p=\mathrm{exp}\left\{\frac{\theta _p}{2}(a_p^{}b_pb_p^{}a_p)\right\}$$
(26)
for each mode where $`a_p`$ and $`b_p`$ are second quantized operators for right- and lefthanded fermions with momentum $`p`$. In the massless limit, right- and left-handed fermions are projected out by

$$P_{R,L}=\frac{1\pm \gamma ^5}{2}.$$
(27)
Chiral transformations are defined as

$$\psi \mathrm{e}^{\mathrm{i}\alpha \gamma ^5}\psi ,$$
(28)
and represent a symmetry of the massless theory (note that $`\gamma ^5`$ can be defined in an odd number of space dimensions only). To see the meaning of right/left-handedness in 1+1 dimensions, we specialize the above spinors to $`m=0`$, where

$$\theta _p=\mathrm{\Theta }(p)\pi $$
(29)
and hence

$$u(p)=\left(\begin{array}{c}\mathrm{\Theta }(p)\\ \mathrm{\Theta }(p)\end{array}\right)$$
(30)
($`\mathrm{\Theta }(p)`$ denotes the Heaviside step function). Thus handedness in two dimensions is correlated with the direction of motion (“right- and left movers”). The Lorentz transformation is consistent with the above definition of kinematic helicity $`s=\pm 1/2`$ for right- and left-handed fermions, respectively,

$$\psi (x)\mathrm{e}^{\xi \gamma ^5/2}\psi (x^{}),$$
(31)
with $`x^{}`$ as given in Eq. (4). Another important peculiarity of fermions in two dimensions concerns the vector and axial currents,

$$j_V^\mu =\overline{\psi }\gamma ^\mu \psi ,j_A^\mu =\overline{\psi }\gamma ^\mu \gamma ^5\psi $$
(32)
and their (partial) conservation laws,

$$_\mu j_V^\mu =0,_\mu j_A^\mu =2\mathrm{i}m\overline{\psi }\gamma ^5\psi .$$
(33)
Due to the severe restrictions for $`\gamma `$-matrices in two dimensions, these currents have actually only two independent components,

$$j_V^0=j_A^1,j_V^1=j_A^0,$$
(34)
or

$$j_A^\mu =ϵ^{\mu \nu }j_{V,\nu }.$$
(35)
This intimate relationship between vector and axial currents is the key for understanding the existence of massless mesons and baryons (see Sec. 3.3).

### 2.3 Scalar bosons

Bosons are quantized canonically with equal-time commutation relations

$$[\pi (x),\varphi (y)]=\mathrm{i}\delta (xy).$$
(36)
Whereas the free massive scalar field has no special features as compared to 3+1 dimensions and satisfies the Klein-Gordon equation

$$(_\mu ^\mu +m^2)\varphi =0,$$
(37)
massless bosons are quite delicate. The simplest object in 3+1 dimensional quantum field theory, the massless scalar field, becomes the most subtle one in 1+1 dimensions due to severe infrared divergences.$`^\mathrm{?}`$ As pointed out by Coleman,$`^\mathrm{?}`$ formally, the integral appearing in the free two-point function

$$0|\varphi (x)\varphi (0)|0=\frac{\mathrm{d}^2k}{(2\pi )}\mathrm{e}^{\mathrm{i}kx}\delta (k^2)\mathrm{\Theta }(k_0)$$
(38)
is IR divergent, as can be seen by performing the $`k^0`$-integration,

$$_0^{\mathrm{}}\frac{\mathrm{d}k^1}{2\pi k^1}\mathrm{cos}(k^1x^1)\mathrm{e}^{\mathrm{i}k^1x^0}.$$
(39)
Since this is also the simplest example of a Goldstone boson in 3+1 dimensions (broken symmetry: $`\varphi \varphi +c`$),$`^\mathrm{?}`$ one expects trouble with the Goldstone theorem in two dimensions. This is indeed what happens: There is a rigorous no-go theorem which forbids SSB of a continuous symmetry in two dimensions (Coleman theorem $`^\mathrm{?}`$). A similar phenomenon is known in statistical physics: there is no long range order in two-dimensional systems (at least with short range interactions $`^\mathrm{?}`$). The strong infrared fluctuations destroy the order, for instance in a crystal, a magnet or a superfluid. There is a different type of long range order though which we discuss in the next subsection.

### 2.4 Long range order

The absence of long range order in 1+1 dimensional QFT or two-dimensional statistical physics could be a fatal blow to our investigation — what interesting phase structure can possibly be left? To understand the way out, let us first come back to the essence of two-dimensionality alluded to in the preceding subsection: Long range order cannot be maintained over arbitrarily large distances due to strong fluctuations of the Goldstone mode; there is no SSB. A nice example for this effect are atoms in a plane which want to form a crystal.$`^\mathrm{?}`$ The amplitude of fluctuations of each atom around its equilibrium position grows logarithmically with $`L`$, if $`L`$ is the extension of the crystal. The calculation involves a similar infrared divergence as in the 2-point function (38),

$$\stackrel{}{u}^2\frac{\mathrm{d}^2k}{(2\pi )^2}\frac{1}{\omega ^2(k)}_{2\pi /L}^\mathrm{\Lambda }\frac{\mathrm{d}k}{k},$$
(40)
($`\mathrm{\Lambda }=1/a`$, inverse lattice spacing). In order to circumvent this no-go theorem one can go to the large $`N`$ limit, invoking a diverging number of degrees of freedom at each point. This was demonstrated by Witten $`^{\mathrm{?},\mathrm{?}}`$ for the chiral GN model, a case very pertinent to the present study. After integrating out the fermions with the help of a complex scalar auxiliary field $`\mathrm{\Phi }=\mathrm{\Phi }_0\mathrm{e}^{\mathrm{i}\theta }`$, one gets the following effective low-energy action for the Goldstone mode $`\theta `$,

$$_{\mathrm{eff}}=\frac{N}{4\pi }(_\mu \theta )^2.$$
(41)
Chiral symmetry breaking can be probed by the correlator

$$\overline{\psi }\psi (x)\overline{\psi }\psi (y)\mathrm{\Phi }^{}(x)\mathrm{\Phi }(y).$$
(42)
It approaches a constant $`|\mathrm{\Phi }_0|^2`$ for $`|xy|\mathrm{}`$ if the symmetry is broken. Without SSB, one expects an exponential decay

$$\mathrm{\Phi }^{}(x)\mathrm{\Phi }(y)\mathrm{e}^{M|xy|}.$$
(43)
Here, a different behavior is found,

$$\mathrm{\Phi }^{}(x)\mathrm{\Phi }(y)\mathrm{e}^{\mathrm{i}\theta (x)}\mathrm{e}^{\mathrm{i}\theta (y)}\mathrm{e}^{2\pi G(xy)/N}\frac{1}{|xy|^{1/N}}.$$
(44)
$`G`$ is the free massless boson propagator, logarithmic in two dimensions. Exponentiation of this propagator yields a power law behavior for the propagator. Evidently, there is no SSB at any finite $`N`$, in agreement with the no-go theorem, but there is a loophole at $`N=\mathrm{}`$ since

$$\mathrm{\Phi }^{}(x)\mathrm{\Phi }(y)1+\frac{1}{N}\mathrm{ln}|xy|+O(1/N^2).$$
(45)
The correlator becomes more and more long-range until it cannot be distinguished from a constant, in the limit $`N\mathrm{}`$. A similar phenomenon is known from two-dimensional $`xy`$ spin-model $`^{\mathrm{?},\mathrm{?},\mathrm{?}}`$ under the name of topological order.

### 2.5 Gauge fields

Gauge symmetry implies redundant variables; this is nowhere as clear as in a two dimensional world where, due to the absence of transverse directions, one might guess that all degrees of freedom of a gauge field are “pure gauge”. This is not quite true — there are left-over quantum mechanical degrees of freedom stemming from gauge invariant zero-mode fields. They can only be discussed reasonably well on a finite interval of length $`L`$. Evidently, pure gauge theories are a totally different story in 1+1 and 3+1 dimensions, as are the corresponding interactions between static charges.

#### 2.5.1 Abelian gauge fields

For the U(1) gauge theory, Lagrangian and field strength tensor look as usual,

$$=\frac{1}{4}F_{\mu \nu }F^{\mu \nu },F_{\mu \nu }=_\mu A_\nu _\nu A_\mu .$$
(46)
In two dimensions however, an antisymmetric tensor has only one independent component, a Lorentz scalar (the electric field); there is no magnetic field. Canonical quantization in the Weyl gauge ($`A_0=0)`$ reveals that the electric field $`E=_0A^1`$ is conjugate to (minus) the vector potential $`A^1`$,

$$[E(x),A^1(y)]=\mathrm{i}\delta (xy).$$
(47)
The Gauss law in this gauge has to be implemented as a constraint on the physical states (this is possible since $`G(x)`$ is conserved),

$$G(x)|=_1E(x)|=0.$$
(48)
This leaves only the 0-mode of $`A_1`$, the spatially constant electric field, as physical variable. Switching to finite $`L`$, we have

$$a=\frac{1}{L}_0^LdxA^1(x),e=_0^LdxE(x)$$
(49)
with

$$H=\frac{1}{2L}e^2,[e,a]=\mathrm{i}.$$
(50)
Since the integer part of $`\frac{gLa}{2\pi }`$, $`g`$ being the electric charge of static sources, can also be gauged away by a periodic gauge transformation

$$U=\mathrm{e}^{\mathrm{i2}\pi nx/L},$$
(51)
pure QED reduces to quantum mechanics of a particle on a circle.$`^\mathrm{?}`$ If we couple static charges to the gauge field, this changes Gauss’s law into $`G=_xEe\rho `$. The “longitudinal” electric field (i.e., the $`x`$-dependent part) is the Coulomb field of the static charge, which is linear in one space dimension. For comparison with Yang-Mills theory below, assume two charges $`\pm g`$ on a circle of length $`L`$ at points $`x,y`$; then the change in energy is

$$\mathrm{\Delta }E=g^2K(xy)$$
(52)
with the periodic Coulomb potential

$`K(xy)`$ $`=`$ $`x|{\displaystyle \frac{1}{_x^2}}|y={\displaystyle \frac{L}{4\pi ^2}}{\displaystyle \underset{n0}{}}{\displaystyle \frac{1}{n^2}}\mathrm{e}^{\mathrm{i2}\pi n(xy)/L}`$ (53)
$`=`$ $`{\displaystyle \frac{L}{12}}+{\displaystyle \frac{1}{2}}|xy|{\displaystyle \frac{1}{2L}}(xy)^2.`$
Pure QED is confining in 1+1 dimensions for rather trivial reasons. In the limit $`L\mathrm{}`$, the fact that the potential to order $`g^2`$ is linear also follows on purely dimensional grounds, since $`[g]=L^1`$. Finally, we note that one has to impose the residual Gauss law ($`Q|=0`$); only the charge 0 sector survives on the circle.

#### 2.5.2 Non-Abelian gauge fields

Consider the Lagrangian of Yang-Mills theory with a SU(N) gauge group on a circle,$`^{\mathrm{?},\mathrm{?},\mathrm{?}}`$

$$=\frac{1}{4}F_{\mu \nu }^aF^{a\mu \nu },$$
(54)
and

$$F_{\mu \nu }^a=_\mu A_\nu ^a_\nu A_\mu ^agf^{abc}A_\mu ^bA_\nu ^c.$$
(55)
Whereas canonical quantization in the Weyl gauge is as straightforward as in the Abelian case, the resolution of Gauss’s law becomes more intricate due to its non-linearity

$$G^a(x)=(DE)^a(x)=_xE^agf^{abc}A^{1,b}E^c.$$
(56)
Generically, the covariant derivative has $`N1`$ zero modes; the projections of $`E`$ onto these are the physical (quantum mechanical) variables $`e^p`$.$`^\mathrm{?}`$ They are non-hermitean due to the projection onto $`A^1`$-dependent basis vectors. The corresponding gauge invariant coordinates are the eigenphases of the path ordered integrals (untraced Polyakov loops) around the compact space direction,

$$𝒫=P\mathrm{e}^{\mathrm{i}g_0^LdxA^1(x)}=V\mathrm{e}^{\mathrm{i}gaL}V^{},a=a^pt^p.$$
(57)
One finds

$$[e^p,a^q]=\mathrm{i}\delta _{pq}$$
(58)
and the Hamiltonian

$$H_a=\frac{1}{2L}e^pe^p=\frac{1}{2L}\frac{1}{𝒥(a)}\frac{}{a^p}𝒥(a)\frac{}{a^p}.$$
(59)
This form makes explicit the fact that the variables are curvilinear coordinates — $`H_a`$ is just the Laplacian on the SU(N) group manifold, the Jacobian $`𝒥(a)`$ being the reduced Haar measure (Vandermonde determinant for unitary matrices)

$$𝒥(a)=\underset{i>j}{}\mathrm{sin}^2\left(\frac{1}{2}gL(a_{ii}a_{jj})\right).$$
(60)
Convenient angular variables are

$$\phi _i=gLa_{ii}$$
(61)
with the constraint $`_i\phi _i=0`$ for SU(N). If one introduces the analogue of a “radial wavefunction”

$$\mathrm{\Phi }(\phi )=\frac{u(\phi )}{\sqrt{𝒥(\phi )}},$$
(62)
the Laplacian is reduced to cartesian form (except for a center of mass correction in the SU(N) case, as opposed to U(N)), and a constant effective potential appears,

$$H=\frac{g^2L}{4}\left(\underset{i=1}{\overset{N}{}}\frac{^2}{\phi _i^2}\frac{1}{N}\left(\underset{i=1}{\overset{N}{}}\frac{}{\phi _i}\right)^2\right)\frac{Ng^2L(N^21)}{48}.$$
(63)
The reduced wave function $`u(\phi )`$ now has to satisfy the boundary conditions

$$u(\phi )=0\mathrm{if}\phi _i=\phi _j.$$
(64)
The ground state is simply $`\mathrm{\Phi }_0=`$ const., normalizable owing to the compactness of the variables, with $`E_0=0`$. To get the excited states for any $`N`$, it is advantageous to reinterpret this problem in terms of free, non-relativistic fermions on a circle $`^\mathrm{?}`$ (the Pauli principle is enforced by Eq. (64)) with periodic (antiperiodic) boundary conditions for odd (even) $`N`$. Then the wavefunctions are Slater determinants of single particle wave functions

$$\psi _m(\phi )=\frac{1}{\sqrt{2\pi }}\mathrm{e}^{\mathrm{i}m\phi }$$
(65)
with $`m`$ integer for odd, half integer for even $`N`$. The energy is

$$E=\frac{g^2L}{4}\left(\underset{m}{\overset{\mathrm{occ}}{}}m^2\frac{1}{N}\left(\underset{m}{\overset{\mathrm{occ}}{}}m\right)^2\right)\frac{Ng^2L(N^21)}{48},$$
(66)
where the sums run over all occupied orbits. By filling the lowest $`N`$ single particle states, one can easily ascertain that the ground state wavefunction is constant and that the ground state energy vanishes. The configuration space is determined by the fundamental domain, the smallest region bounded by zeros of the Jacobian. We can choose

$$\phi _1<\phi _2<\mathrm{}<\phi _N$$
(67)
and the center of mass constraint $`_i\phi _i=0`$ for $`SU(N)`$ (this represents a more complete gauge fixing). Since the ordering is on a circle rather than a line, it still leaves a residual gauge freedom, namely cyclic permutations of color labels. This residual $`Z_N`$ symmetry is expected on topological grounds due to

$$\mathrm{\Pi }_1(SU(N)/Z_N)=Z_N.$$
(68)
This is the well-known center symmetry of pure Yang Mills theory which plays an important role at finite temperature.$`^\mathrm{?}`$ By way of illustration, take SU(2): There is only one independent angular variable, $`\phi =(\phi _1\phi _2)/2[0,\pi ]`$. Then, YM2 reduces to quantum mechanics of a particle in an infinite square well with

$$H=\frac{g^2L}{8}\frac{^2}{\phi ^2}\frac{g^2L}{8}$$
(69)
and

$$u_n(\phi )=\sqrt{\frac{2}{\pi }}\mathrm{sin}n\phi ,𝒥(\phi )=\mathrm{sin}^2\phi ,E_n=\frac{g^2L}{8}(n^21).$$
(70)
The Z2 center symmetry is just parity. A detailed discussion of the SU(3) case with two independent variables can be found in Ref. 27. At very large $`N`$, since the $`\phi _i`$ repel each other (cf. Eq. 64)) and live on a circle, they will tend to distribute themselves along the circle like a pearl necklace, with fluctuations which are $`1/N`$ suppressed in amplitude. This observation will become important later on for the ’t Hooft model at finite temperature.

The existence of zero mode gluons also complicates the interaction between static charges, at least for finite $`L`$. If one puts a static $`q\overline{q}`$-pair on a circle, the Coulomb potential evidently involves $`D^2`$ (covariant derivatives) rather than $`^2`$, so that the Coulomb interaction gets modified by the gauge field remnants. The Coulomb potential can be transformed to coordinate space,$`^{\mathrm{?},\mathrm{?}}`$

$`K_{ij}(z)`$ $`=`$ $`L{\displaystyle \underset{n}{}}{\displaystyle \frac{1}{(2\pi n\phi _j+\phi _i)^2}}\mathrm{e}^{\mathrm{i2}\pi nz/L}`$
$`=`$ $`{\displaystyle \frac{L}{4}}\mathrm{e}^{\mathrm{i}(\phi _j\phi _i)z/L}\left({\displaystyle \frac{1}{\mathrm{sin}^2(\phi _i\phi _j)/2}}{\displaystyle \frac{2}{L}}|z|+\mathrm{i}{\displaystyle \frac{2}{L}}z\mathrm{cot}(\phi _i\phi _j)/2\right)`$
for $`ij`$, whereas it is identical to the Abelian one, Eq. (53), for $`i=j`$. Since the $`\phi _i`$ are dynamical variables, one is not yet through but still has to solve a fairly complicated quantum mechanical $`(N1)`$-body problem with a two-body potential defined by the Coulomb interaction. This reflects the fact that due to color spin dynamics, charges are never static in the non-Abelian case, even if they are nailed down in space. A full analytic solution is available for the SU(2) case,$`^\mathrm{?}`$ where a dynamical quantum mechanical supersymmetry is at work.$`^\mathrm{?}`$ It explains why the potential is strictly linear (at all $`L`$, not just in the limit $`L\mathrm{}`$), and why – for vanishing separation of the charges – all excited states are doubly degenerate whereas the ground state is non-degenerate and has strictly zero energy, see Fig. 1.

In the limit $`L\mathrm{}`$, all excited states move up to infinite energy, the gauge field dynamics get frozen, and the difference between QED2 and YM2 disappears.

The discussion about the zero mode gauge fields may seem a bit pedantic: Why bother about finite $`L`$ at all? The reason is of course that we are interested among other things in the finite temperature problem. The physics relevance of the Polyakov loop is beyond doubt there.

Finally, the fact that the left-over gluon variables are eigenphases of path-ordered integrals has been used in the literature to emphasize the string picture: Closed strings that wrap around the compact direction in pure Yang-Mills theory, open strings with quarks at the ends if one includes matter fields in the fundamental representation.$`^{\mathrm{?},\mathrm{?},\mathrm{?}}`$ We shall not exploit this picture here any further.

### 2.6 Renormalization

From the kinetic term in the Lagrangian in two dimensions, one infers that Dirac fermion, scalar boson and vector fields have dimensions

$$[\psi _i]=L^{1/2},[\varphi ]=[A^\mu ]=L^0.$$
(72)
In self-interacting fermion theories, a quartic self-interaction has dimensionless coupling constant and is therefore perturbatively renormalizable (GN model) unlike in four dimensions (NJL model). Self-interacting bosons can have interactions of arbitrary high order with couplings of positive mass dimension and are all super renormalizable (e.g. sine-Gordon model). The gauge coupling has dimension of mass, therefore gauge theories like the Schwinger model or QCD2 are also super-renormalizable.

In the U(1) theory, an axial anomaly occurs, as is well known from the Schwinger model. In one-flavor SU(N) QCD2, there is no axial anomaly.

Pure fermionic theories of GN type are asymptotically free; they share this property with QCD and are rather exceptional in this respect. The $`\beta `$-function is negative (see Sec. 3.1.1). If the bare mass is put equal to zero (as in the original GN model), dimensional transmutation occurs like in pure Yang-Mills theory in four dimensions: A massless theory generates mass dynamically. All dimensionless ratios can be predicted without free parameter, only the physical fermion mass is needed as an input.$`^\mathrm{?}`$

In the massive GN model, there are two distinct parameters of dimension mass (matching the two bare parameters $`Ng^2,m_0`$): The physical fermion mass $`m`$ and the ratio $`m_0/(Ng^2)`$. The first parameter is directly observable via the spectrum. The second parameter is also observable: In the chiral GN model, it determines the mass of the would-be Goldstone boson, cf. Sec. 3.3. In the non-chiral GN model, it can be measured via the $`\overline{q}q`$ scattering length which is proportional to $`Ng^2/m_0`$. In the limit $`m_00`$, the pion mass goes to zero and the scattering length in the non-chiral model diverges. In this way, the second parameter disappears and all observables can be expressed in terms of $`m`$ only.

QCD2 is trivially asymptotically free, being super-renormalizable. The perturbative (short distance) interaction is determined by dimensional considerations.

### 2.7 Soluble models of interest for QCD

The first generation of soluble models like the Thirring model $`^\mathrm{?}`$ or the massless Schwinger model $`^\mathrm{?}`$ from $`\pm `$1960 were not yet rich enough for the questions we are interested in. Although they can be solved exactly and have certain interesting non-perturbative aspects, they essentially describe free theories with a trivial $`S`$-matrix. If one tries to make them non-trivial by adding e.g. a mass term or flavor, one loses solvability. The Schwinger model has confinement, but provides no model for hadrons.

From this point of view, ’t Hooft’s idea of the $`1/N`$ expansion $`^\mathrm{?}`$ is very nice indeed. By invoking a large number of fermion species and using $`1/N`$ as a small parameter, one finds highly non-trivial results already to leading order. Semiclassical methods become applicable and yield systematic, reliable results. All symmetries like Poincaré invariance, chiral symmetry, gauge invariance are preserved. The phenomenological arguments in favor of the large $`N`$ limit for QCD (Zweig rule etc.$`^{\mathrm{?},\mathrm{?},\mathrm{?}}`$) naturally play no role in two dimensions, where it is a purely theoretical tool.

Models based on the large $`N`$ expansion can be considered as second generation field theory models, also historically. The two models we are particularly interested in came up in 1974, the GN $`^\mathrm{?}`$ and the ’t Hooft model.$`^\mathrm{?}`$ The GN model generalizes the Thirring model to $`N`$ flavors, while the ’t Hooft model generalizes the Schwinger model to the non-Abelian SU(N) gauge group with a large number of colors. They acquire a rich physical content, as we hope to convince the reader in the remainder of this article. Here, we only mention the following basic facts: The GN models describe self-interacting massless fermions with a point-like four-fermion interaction. Depending on the particular variant, they can have either a discrete or a continuous chiral symmetry; in the latter case, the GN model may be thought of as the NJL model $`^\mathrm{?}`$ in two dimensions. Mass terms which violate chiral symmetry explictly can be included as well without loss of solvability (at large $`N`$). The ’t Hooft model is QCD2 with SU(N) color in the limit $`N\mathrm{}`$, with massless or massive quarks in the fundamental representation. As mentioned above, the model on the infinite line has no gluon degrees of freedom, but is also a self-interacting fermion theory, here with a linear Coulomb potential characteristic for the two dimensional world (in an appropriate gauge). Unlike the GN models, this model shows confinement, which makes it particularly interesting from the point of view of QCD. The fact that also in real life (in four dimensions), the NJL model and QCD are used in parallel nowadays for certain questions of strong interaction physics makes this pair of models a good match and immediately suggests many questions to be asked. Throughout this work, we do not consider any embellishments or generalizations (e.g. flavor, more complicated interactions, different matter fields) of these models which have occasionally been invoked. The nice thing about the original models is just the contrast between the simplicity of the Lagrangian and the complexity of the phenomena it produces, and we do not want to spoil this property.

Finally, we should point out that the main goal of the present study are questions of hot and dense matter in these toy models. For this purpose, the large $`N`$ limit is not only a technical trick leading to solvability of the models, but is in fact instrumental for enabling chiral symmetry breakdown in two dimensions. This is the reason why another aspect of the GN model will play no role here, its integrability.$`^{\mathrm{?},\mathrm{?},\mathrm{?}}`$ A lot of progress has been made towards solving the finite $`N`$ model exactly. Here, we shall restrict ourselves to the large $`N`$ limit so that the corresponding difference between the non-integrable ’t Hooft model and the integrable GN models will not play any role, but a common approach for both models will be used throughout.

## 3 Two dimensional models at zero temperature and chemical potential

### 3.1 Gross-Neveu models

We cannot discuss finite temperature or density without first reviewing the particle content of the theory. Originally, the GN model has been solved with the help of semi-classical functional methods,$`^\mathrm{?}`$ supplemented by inverse scattering techniques from soliton theory for the baryons.$`^\mathrm{?}`$ This approach can be rigorously justified in the large $`N`$ limit. Typically, one introduces auxiliary bosonic fields and integrates out the fermions (Gaussian integral over Grassmann variables), then applies a saddle point approximation to the remaining bosonic functional integral. More recently, an alternative to the inverse scattering method has been proposed, based on the Gel’fand-Dikii equation.$`^\mathrm{?}`$ Here, in order to have a coherent scheme for all questions that interest us as well as for pedagogical reasons, we shall make use of time honored concepts well known from atomic and nuclear physics, such as the HF and random phase approximation (RPA), although in a fully relativistic setting.$`^{\mathrm{?},\mathrm{?},\mathrm{?}}`$ In contrast to the functional integral methods, we thereby handle fermions more directly and, since we work canonically, can also more readily address questions of the vacuum structure (Dirac sea). Needless to say, the equations one eventually solves are always the same (HF or Schwinger-Dyson, RPA or Bethe-Salpeter equations). More than anything else the choice is a question of language.

#### 3.1.1 Vacuum, physical fermion

The Lagrangian of the GN model family, including a bare fermion mass term and suppressing color labels, reads $`^\mathrm{?}`$

$$=\overline{q}\mathrm{i}\gamma ^\mu _\mu q+\frac{1}{2}g^2\left[\left(\overline{q}q\right)^2\lambda \left(\overline{q}\gamma ^5q\right)^2\right]m_0\overline{q}q.$$
(73)
The original model with discrete chiral symmetry $`q\gamma ^5q`$ is recovered for $`\lambda =0,m_0=0`$; the choice $`\lambda =1,m_0=0`$ corresponds to the NJL-type model with continuous chiral symmetry; the $`m_0`$ term breaks both chiral symmetries explicitly. To leading order in $`1/N`$, the Hartree approximation can be used, replacing $`\overline{q}q`$ in the Euler-Lagrange equation

$$\mathrm{i}\gamma ^\mu _\mu q+g^2\left(\overline{q}q\lambda (\overline{q}\gamma ^5q)\gamma ^5\right)qm_0q=0$$
(74)
by its vacuum expectation value (the Fock term is $`1/N`$ suppressed). The non-zero value of the order parameter for chiral symmetry signals SSB. We first assume $`\overline{q}\gamma ^5q=0`$,

$$\left\{\mathrm{i}\gamma ^\mu _\mu m_0+Ng^2\rho _s\right\}q(x)=0.$$
(75)
A physical fermion mass appears,

$$m=m_0Ng^2\rho _s,\rho _s=\frac{1}{N}\overline{q}q,$$
(76)
as dictated by covariance. The essence of the Hartree approximation is to impose a self-consistency condition, here on the scalar density $`\rho _s`$ of the vacuum. Using the free fermion spinors of Eq. (24) above we get

$`\rho _s={\displaystyle \frac{(mm_0)}{Ng^2}}`$ $`\stackrel{!}{=}`$ $`{\displaystyle _{\mathrm{\Lambda }/2}^{+\mathrm{\Lambda }/2}}{\displaystyle \frac{\mathrm{d}k}{2\pi }}\overline{v}(k)v(k)`$ (77)
$`=`$ $`{\displaystyle _{\mathrm{\Lambda }/2}^{+\mathrm{\Lambda }/2}}{\displaystyle \frac{\mathrm{d}k}{2\pi }}{\displaystyle \frac{m}{\sqrt{k^2+m^2}}}{\displaystyle \frac{m}{2\pi }}\mathrm{log}{\displaystyle \frac{\mathrm{\Lambda }^2}{m^2}}.`$
As usual, this is also a variational solution obtained by minimizing the vacuum energy in the space of Slater determinants. The Hartree- or gap equation,

$$m=m_0+m\frac{Ng^2}{2\pi }\mathrm{ln}\frac{\mathrm{\Lambda }^2}{m^2},$$
(78)
contains much information. For $`m_0=0`$ (discrete or continuous model), it has two solutions, either $`m=0`$ or

$$1=\frac{Ng^2}{2\pi }\mathrm{ln}\frac{\mathrm{\Lambda }^2}{m^2}.$$
(79)
The broken symmetry solution is always lower in energy density by

$$_{\mathrm{vac}}_{\mathrm{vac}}^{\mathrm{free}}=\frac{Nm^2}{4\pi }.$$
(80)
The gap equation for the broken phase does not determine $`m`$, as is clear from the fact that the theory contains no dimensionful parameter. If $`m`$ is adjusted (to the “observed” fermion mass, whatever that means in two dimensions), then the gap equation is a renormalization condition which tells us how $`Ng^2`$ (the bare coupling constant) depends on the cutoff $`\mathrm{\Lambda }`$. Correspondingly, the physical mass has the dependence

$$m=\mathrm{\Lambda }\mathrm{exp}\left\{\frac{\pi }{Ng^2}\right\}$$
(81)
familiar from similar dependences in real QCD and in line with the renormalization group. All physical quantities with dimension of mass have to depend in the same way on $`g,\mathrm{\Lambda }`$. The essential singularity in $`g`$ at $`g=0`$ underlines the power of the $`1/N`$ expansion as compared to ordinary perturbation theory in $`g`$. Dimensionless quantities can be predicted without any parameter, at $`m_0=0`$. Apart from this phenomenon of dimensional transmutation, the other remarkable feature which we find here is asymptotic freedom (AF).$`^{\mathrm{?},\mathrm{?}}`$ One can determine the $`\beta `$-function of the GN model in an elementary way as follows: The physical fermion mass $`m`$ should be cut-off independent, therefore

$$\mathrm{\Lambda }\frac{\mathrm{d}}{\mathrm{d}\mathrm{\Lambda }}m(\mathrm{\Lambda },g(\mathrm{\Lambda }))=\left(\mathrm{\Lambda }\frac{}{\mathrm{\Lambda }}+\mathrm{\Lambda }\frac{\mathrm{d}g}{\mathrm{d}\mathrm{\Lambda }}\frac{}{g}\right)m=0.$$
(82)
Set

$$\mathrm{\Lambda }\frac{\mathrm{d}g}{\mathrm{d}\mathrm{\Lambda }}=\beta (g)$$
(83)
and find, by inserting $`m`$ from Eq. (81),

$$\beta (g)=\frac{Ng^3}{2\pi }<0$$
(84)
or AF (to this order, the result is independent of the renormalization scheme). More elaborate computations of the higher loop beta-functions $`^{\mathrm{?},\mathrm{?}}`$ do not play a role in the large $`N`$ limit. We shall return to the issue of AF in connection with $`q\overline{q}`$ scattering, where a running coupling constant will be defined.

For $`m_0=0`$ and the case of continuous chiral symmetry, the assumption $`\overline{q}\gamma ^5q=0`$ is too restrictive; one gets a continuum of vacua related by U(1) transformations and labeled by a chiral angle $`\phi `$,

$$\overline{q}q=\frac{m}{g^2}\mathrm{cos}\phi ,\overline{q}\mathrm{i}\gamma ^5q=\frac{m}{g^2}\mathrm{sin}\phi .$$
(85)
For non-vanishing bare mass $`m_0`$, the gap equation (78) signals the appearance of a second physical parameter $`\frac{m_0}{Ng^2}`$, in addition to the overall scale set by $`m`$. Its relevance for scattering and bound state observables will become clear later on.

#### 3.1.2 Baryons

It would be premature to conclude from the preceding section that the GN model at $`N\mathrm{}`$ reduces to a free, massive fermion theory. As observed by Witten,$`^\mathrm{?}`$ baryons have to be considered at large $`N`$ and can be treated in HF approximation, like the vacuum. They correspond to a different kind of HF solution which breaks translational invariance. Their mass diverges like $`N`$, the baryon-baryon interaction also scales like $`N`$, and these effects have to be taken into account in leading order. Baryons were found originally with functional techniques and inverse scattering methods.$`^\mathrm{?}`$ Since we now know the scalar potential, we can use it to verify Witten’s picture of baryons as relativistic HF solutions — this is what we propose to do here.$`^\mathrm{?}`$ We consider the discrete chiral model first, since the continuous model has in addition totally different, light baryons, see Sec. 3.3 below.

If we restrict ourselves to baryons made up of $`nN`$ quarks, it is sufficient to take into account one positive energy “valence” level filled with $`n`$ fermions, together with the completely filled Dirac sea. The Hartree equation without invoking translational invariance becomes

$$\left\{\mathrm{i}\gamma ^\mu _\mu +Ng^2\rho _s(x)\right\}q(x)=0$$
(86)
where the scalar density now refers to the baryon state $`|B`$ ,

$$\rho _s(x)=\frac{1}{N}\underset{i=1}{\overset{N}{}}B|\overline{q}_i(x)q_i(x)|B.$$
(87)
We introduce the single particle energies and eigenfunctions of the Dirac Hamiltonian corresponding to Eq. (86),

$$\left\{\gamma ^5\frac{1}{\mathrm{i}}\frac{}{x}\gamma ^0Ng^2\rho _s(x)\right\}\psi _\alpha ^{(\pm )}(x)=\pm ϵ_\alpha \psi _\alpha ^{(\pm )}(x).$$
(88)
Then $`\rho _s(x)`$ is given by

$$\rho _s(x)=\frac{n}{N}\overline{\psi }_0^{(+)}(x)\psi _0^{(+)}(x)+\underset{\alpha }{}\overline{\psi }_\alpha ^{()}(x)\psi _\alpha ^{()}(x)$$
(89)
where $`\psi _0^{(+)}`$ is the valence state. Dashen et al. $`^\mathrm{?}`$ have found two different types of solutions for the scalar density, the “kink” and the “double kink” (or rather kink-antikink). Let us start with the less exotic double kink, characterized by the following scalar potential:

$$S(x)Ng^2\rho _s(x)=m\left\{1+y\left(\mathrm{tanh}\xi _{}\mathrm{tanh}\xi _+\right)\right\}$$
(90)
with the definitions

$$y=\mathrm{sin}\theta ,\theta =\frac{\pi }{2}\frac{n}{N},$$
(91)
$$\xi _\pm =ymx\pm \frac{1}{2}\text{artanh}y.$$
(92)
$`S(x)`$ has the form of a potential well and approaches for $`x\pm \mathrm{}`$ the asymptotic value $`m`$, i.e., the physical fermion mass. Varying $`y`$ or equivalently the degree of occupation $`n/N`$, $`S(x)`$ changes from a $`1/\mathrm{cosh}^2(ymx)`$ shape at small $`y`$ to a Woods-Saxon like shape at $`y1`$,

$$S(x)\left(1+e^{(|x|R)/c}\right)^1,$$
(93)
with skin thickness $`c`$ and radius $`R`$ given by

$$c=\frac{1}{2m},R=\frac{1}{2m}\mathrm{log}\frac{\mathrm{cos}\theta }{2}(y1),$$
(94)
respectively. We proceed to verify the self-consistency of this scalar potential. After a straightforward solution of the Dirac equation (86) with potential (90), one finds the following results: The normalized wavefunctions of the discrete (valence) states are given by (in a basis where $`\gamma ^0=\sigma _1,\gamma ^1=\mathrm{i}\sigma _3`$)

$$\psi _0^{(\pm )}(x)=\frac{\sqrt{ym}}{2}\left(\begin{array}{c}\hfill \frac{1}{\mathrm{cosh}\xi _{}}\\ \hfill \frac{1}{\mathrm{cosh}\xi _+}\end{array}\right)$$
(95)
with eigenvalues

$$E_0^{(\pm )}=\pm m\sqrt{1y^2}.$$
(96)
The continuum states are

$$\psi _k^{(\pm )}(x)=\frac{1}{\sqrt{2}E(k)(\mathrm{i}k+ym)}\left(\begin{array}{c}\hfill (\mathrm{i}km)(\mathrm{i}kym\mathrm{tanh}\xi _{})\\ \hfill \pm E(k)(\mathrm{i}kym\mathrm{tanh}\xi _+)\end{array}\right)\mathrm{e}^{\mathrm{i}kx},$$
(97)
with

$$E(k)=\sqrt{k^2+m^2}.$$
(98)
They are reflectionless, a characteristic feature of time independent solutions.$`^\mathrm{?}`$ $`^{\mathrm{?},\mathrm{?}}`$ The spectrum of the Dirac Hamiltonian therefore consists of the usual positive and negative energy continua starting from $`\pm m`$ and a pair of discrete states inside the mass gap at energies given by Eq. (96), see Fig. 2.

These latter states are of course a prerequisite for obtaining bound baryons. The discrete negative energy state has to be filled completely ($`N`$ fermions), its positive energy partner with $`nN`$ fermions (remember that, unlike in QCD, there is no restriction to color singlets). Reinserting Eqs. (95,97) into (90), one can then verify the self-consistency of the solution. This requires the use of the vacuum gap equation (78) to eliminate infinities.$`^\mathrm{?}`$

Finally, one has to compute the mass of the baryon. This is somewhat delicate because of the vacuum subtraction. The single particle energies of two systems must be subtracted: One of the systems has only continuum states, whereas the other system has one extra pair of discrete states. The easiest way of doing this is to enclose the system in a finite box, thereby discretizing all states, and to use the well-known relation between phase shifts and density of states. In this way, one recovers the result of Dashen et al. $`^\mathrm{?}`$ for the baryon mass,

$$M_B=\frac{2Nym}{\pi }=nm\left(\frac{\mathrm{sin}\theta }{\theta }\right).$$
(99)
Let us briefly compare this exact result to the non-relativistic approach (in the present model, the number of valence quarks governs the degree of relativity). In the non-relativistic Hartree approximation, we start from the Schrödinger equation for a particle of mass $`m`$, interacting with the other particles through a $`\delta `$ -function potential. The Hartree equation then reduces to the non-linear Schrödinger equation,$`^\mathrm{?}`$

$$\left\{\frac{1}{2m}\frac{^2}{x^2}ng^2|\varphi _0|^2\right\}\varphi _0(x)=E_0\varphi _0(x),$$
(100)
with the following bound state solution,

$$\varphi _0(x)=\sqrt{\frac{\kappa }{2}}\frac{1}{\mathrm{cosh}\kappa x},$$
(101)
$$\kappa =\frac{1}{2}ng^2m,E_0=\frac{\kappa ^2}{2m}.$$
(102)
The baryon mass turns out to have the value

$$M_B=nm\left(1\frac{1}{6}\left(\frac{ng^2}{2}\right)^2\right).$$
(103)
If we choose

$$g^2=\frac{\pi }{N},$$
(104)
both $`M_B`$ and the scalar density in the full calculation and in the non-relativistic limit agree up to O($`y^2`$). The reason why the GN model corresponds to a specific choice of the coupling constant is that it involves only one parameter, whereas the non-relativistic model has two independent parameters, the fermion mass and the coupling constant.

Now turn to the opposite, extreme relativistic limit. As $`\theta `$ approaches $`\pi /2`$, the effects of the Dirac sea become overwhelming. Eventually, the baryon goes over into a well separated kink anti-kink pair. Although the scalar density becomes flat, the fermion density for the discrete state is concentrated near the surface, since

$$\rho _f^0=\psi _0^{}(x)\psi _0(x)=\frac{ym}{4}\left(\frac{1}{\mathrm{cosh}^2\xi _{}}+\frac{1}{\mathrm{cosh}^2\xi _+}\right).$$
(105)
In the limit $`y1`$, the two kinks are completely decoupled, see Fig. 3.

If one looks at a single kink now, one finds another type of baryon of the GN model, referred to as Callan-Coleman-Gross-Zee kink in the literature.$`^\mathrm{?}`$ It has the scalar potential

$$S(x)=m\mathrm{tanh}\left(xm\right).$$
(106)
This solution has many unusual features: $`S(x)`$ interpolates between the two different vacua related by the discrete $`\gamma ^5`$-transformation. Unlike the double kink baryon which is only stabilized by the binding of fermions, the single kink is topologically stable. The fact that a kink type potential can bind particles is a well-known, purely relativistic effect which has no analogue in non-relativistic quantum mechanics. The discrete state corresponding to Eq. (95) becomes

$$\psi _0(x)=\sqrt{\frac{m}{2}}\left(\begin{array}{c}\frac{1}{\mathrm{cosh}mx}\\ 0\end{array}\right),$$
(107)
with eigenvalue $`E_0=0`$, while the continuum wavefunctions replacing Eq. (97) are

$$\psi _k^{(\pm )}(x)=\frac{1}{\sqrt{2}E(k)}\left(\begin{array}{c}\mathrm{i}km\mathrm{tanh}mx\\ \pm E(k)\end{array}\right)\mathrm{e}^{\mathrm{i}kx}.$$
(108)
The valence state has evidently vanishing scalar density so that there is no feedback at all to the Hartree potential. This is consistent with the fact that the mass of the kink baryon is independent of the number $`n`$ of valence fermions. The self-consistency for the scalar potential can again be verified using the vacuum Hartree condition (78). As expected intuitively, the mass is $`1/2`$ of the double kink mass in the limit $`y1`$

$$M_B|_{\mathrm{kink}}=\frac{Nm}{\pi },$$
(109)
in agreement with Dashen et al. $`^\mathrm{?}`$

#### 3.1.3 Mesons

In this section, we shall rederive the known fermion-antifermion bound state (meson) of the GN model, using the relativistic RPA. This meson has zero binding energy in the large $`N`$ limit, i.e., its mass is just twice the physical fermion mass.$`^\mathrm{?}`$ We discuss only the discrete chiral symmetry case here, without bare mass; the NJL type model with its massless Goldstone mode will be treated in Sec. 3.3. We derive the RPA equations using an equations of motion approach $`^{\mathrm{?},\mathrm{?}}`$ and solve them analytically for mesons with arbitrary momentum, thereby demonstrating explicitly the covariance of the spectrum.$`^\mathrm{?}`$

In RPA, the meson is described in terms of particle-hole excitations on top of the HF vacuum. We therefore start from the equation of motion of the following operator bilinear in the fermion fields,

$$Q_{\alpha \beta }(x,y)=\frac{1}{N}\underset{i}{}q_{i\beta }^{}(y)q_{i\alpha }(x).$$
(110)
Using the equation of motion for $`q`$, Eq. (74), and neglecting $`1/N`$ suppressed terms, we find

$`\mathrm{i}{\displaystyle \frac{}{t}}Q(x,y)`$ $`=`$ $`\mathrm{i}\left\{{\displaystyle \frac{}{y}}Q(x,y)\gamma ^5+{\displaystyle \frac{}{x}}\gamma ^5Q(x,y)\right\}`$
$`Ng^2\left\{\text{tr}\left(\gamma ^0Q(x,x)\right)\gamma ^0Q(x,y)Q(x,y)\gamma ^0\text{tr}\left(\gamma ^0Q(y,y)\right)\right\}`$
Furthermore, we expand the operator $`Q(x,y)`$ around its vacuum expectation value

$$Q(x,y)=0|Q(x,y)|0+\frac{1}{\sqrt{N}}\stackrel{~}{Q}(x,y)+\mathrm{}$$
(112)
where the $`c`$-number part corresponds to the density matrix of the vacuum,

$$0|Q(x,y)|0=\rho (xy).$$
(113)
Inserting (112) into (3.1.3) and equating terms with the same power in $`1/\sqrt{N}`$, we get to zeroth order the following alternative formulation of the Hartree equation,

$$\mathrm{i}\frac{}{t}\rho (x)=\mathrm{i}\frac{}{x}[\gamma ^5,\rho (x)]Ng^2\text{tr}\left(\gamma ^0\rho (0)\right)[\gamma ^0,\rho (x)]=0$$
(114)
and to order $`1/\sqrt{N}`$ the linearized equation of motion for the fluctuation $`\stackrel{~}{Q}`$,

$`\mathrm{i}{\displaystyle \frac{}{t}}\stackrel{~}{Q}(x,y)`$ $`=`$ $`\mathrm{i}\left\{{\displaystyle \frac{}{y}}\stackrel{~}{Q}(x,y)\gamma ^5+{\displaystyle \frac{}{x}}\gamma ^5\stackrel{~}{Q}(x,y)\right\}`$
$`Ng^2\{\text{tr}\left(\gamma ^0\rho (0)\right)\gamma ^0\stackrel{~}{Q}(x,y)\stackrel{~}{Q}(x,y)\gamma ^0\text{tr}\left(\gamma ^0\rho (0)\right)`$
$`+\text{tr}\left(\gamma ^0\stackrel{~}{Q}(x,x)\right)\gamma ^0\rho (xy)\rho (xy)\gamma ^0\text{tr}\left(\gamma ^0\stackrel{~}{Q}(y,y)\right)\}.`$
Switching to momentum variables,

$`\rho (k)`$ $`=`$ $`{\displaystyle dx\mathrm{e}^{\mathrm{i}kx}\rho (x)}`$
$`\stackrel{~}{Q}(k^{},k)`$ $`=`$ $`{\displaystyle dxdy\mathrm{e}^{\mathrm{i}k^{}x+\mathrm{i}ky}\stackrel{~}{Q}(x,y)},`$ (116)
the Hartree equation simplifies to

$$[h(k),\rho (k)]=[k\gamma ^5+m\gamma ^0,\rho (k)]=0$$
(117)
where $`h(k)`$ denotes the single particle Dirac Hamiltonian and where we have again introduced the physical fermion mass, Eq. (76). The Hartree solution of Sec. 3.1.1 corresponds to

$$\rho (k)=v(k)v^{}(k)=\frac{1}{2}\frac{m}{2E(k)}\gamma ^0\frac{k}{2E(k)}\gamma ^5$$
(118)
and obviously satisfies Eq. (117). In order to convert Eq. (3.1.3) into the canonical form of the RPA, we first have to project the operator $`\stackrel{~}{Q}`$ onto self-consistent Hartree spinors. In the present case this simply amounts to a change of basis from massless to massive free spinors. We define the “particle-hole” and “hole-particle” components of the operator $`\stackrel{~}{Q}`$ as

$`\stackrel{~}{Q}_{12}(k^{},k)`$ $`=`$ $`u^{}(k^{})\stackrel{~}{Q}(k^{},k)v(k)`$
$`\stackrel{~}{Q}_{21}(k^{},k)`$ $`=`$ $`v^{}(k^{})\stackrel{~}{Q}(k^{},k)u(k).`$ (119)
Using the commutation relations, the corresponding “particle-particle” and “hole-hole” operators $`\stackrel{~}{Q}_{11}`$ and $`\stackrel{~}{Q}_{22}`$ can be shown to be of higher order in $`1/\sqrt{N}`$ and can therefore be neglected. The resulting system of coupled equations is sandwiched between the vacuum $`|0`$ and a one meson state $`n,P|`$, where $`P`$ denotes the total momentum and $`n`$ the type of meson. We define the transition matrix-elements

$`n,P|\stackrel{~}{Q}_{21}(k^{},k)|0`$ $`=`$ $`2\pi \delta (Pk+k^{})X_n(P,k),`$
$`n,P|\stackrel{~}{Q}_{12}(k^{},k)|0`$ $`=`$ $`2\pi \delta (Pk+k^{})Y_n(P,k).`$ (120)
Then, using the Heisenberg equation of motion in order to replace time derivatives by energies, we obtain

$`_n(P)X_n(P,k)`$ $`=`$ $`E(kP,k)X_n(P,k)`$
$`Ng^2\overline{v}(kP)u(k){\displaystyle }{\displaystyle \frac{\mathrm{d}k^{}}{2\pi }}\{\overline{v}(k^{})u(k^{}P)Y_n(P,k^{})`$
$`+\overline{u}(k^{})v(k^{}P)X_n(P,k^{})\}`$
$`_n(P)Y_n(P,k)`$ $`=`$ $`E(kP,k)Y_n(P,k)`$
$`+Ng^2\overline{u}(kP)v(k){\displaystyle }{\displaystyle \frac{\mathrm{d}k^{}}{2\pi }}\{\overline{v}(k^{})u(k^{}P)Y_n(P,k^{})`$
$`+\overline{u}(k^{})v(k^{}P)X_n(P,k^{})\}`$
with

$$E(k^{},k)=E(k^{})+E(k).$$
(123)
These are the final RPA equations. Taking into account the symmetry relation

$$\overline{u}(k)v(k^{})=\overline{v}(k)u(k^{})$$
(124)
which holds for the spinors (24), they can be cast into the familiar matrix form of the RPA,

$$\left(\begin{array}{cc}\hfill A& \hfill B\\ \hfill B& \hfill A\end{array}\right)\left(\begin{array}{c}X\\ Y\end{array}\right)=\left(\begin{array}{c}X\\ Y\end{array}\right).$$
(125)
It is illustrated graphically in Fig. 4.

Characteristic of the RPA (as opposed to the Tamm-Dancoff approximation (TDA)) are the backward going bubbles (corresponding to $`Y`$-components) which account for ground state correlations. Since the integral kernels are separable, one can solve these equations algebraically. After a few straightforward steps $`^\mathrm{?}`$ one arrives at the eigenvalue equation

$$1=\frac{Ng^2}{2}\frac{\mathrm{d}k}{2\pi }\left(\frac{1}{E(kP)}+\frac{1}{E(k)}\right)\frac{4m^2+P^2E^2(kP,k)}{_n^2(P)E^2(kP,k)}.$$
(126)
Clearly, the choice

$$_0^2(P)=4m^2+P^2$$
(127)
converts this into the vacuum Hartree equation. This shows that there is a marginally bound meson with mass $`2m`$, and that the spectrum is covariant.

It is instructive to compute also the $`q\overline{q}`$ scattering matrix in this approach $`^\mathrm{?}`$. The scattering matrix can easily be found, since we are dealing with a separable potential. With the Mandelstam variable $`s=2P`$, find

$$\frac{1}{\tau (s)}=\frac{\sqrt{2}}{\pi m^2}\sqrt{1\eta }\left(\mathrm{log}\frac{1\sqrt{1\eta }}{1+\sqrt{1\eta }}+\mathrm{i}\pi \right)$$
(128)
in the scattering regime with $`\eta =4m^2/s<1`$. ($`\tau `$ is normalized such that $`\tau =Ng^2m^2/\sqrt{8}`$ in tree approximation). $`\tau (s)`$ has the expected pole at $`s=4m^2`$ corresponding to the threshold bound state. By going to large $`s`$, we can check asymptotic freedom in a very physical manner,

$$\tau (s)\frac{\pi m^2}{\sqrt{2}}\frac{1}{\mathrm{log}(m^2/s)}.$$
(129)
The interaction shows the logarithmic decrease in the UV characteristic for asymptotically free theories. A running coupling constant can be defined by comparison with the tree approximation ($`s=\mu ^2`$),

$$Ng_{\mathrm{eff}}^2(\mu )=\frac{2\pi }{\mathrm{log}\mu ^2/m^2}.$$
(130)
This illustrates once again the usefulness of the GN model as toy model for QCD.

### 3.2 ’t Hooft model

The ’t Hooft model $`^\mathrm{?}`$ is defined as the large $`N`$ limit of 1+1 dimensional SU($`N`$) gauge theory with quarks in the fundamental representation,

$$=\overline{q}\mathrm{i}D/q\frac{1}{2}\mathrm{tr}F_{\mu \nu }F^{\mu \nu },D/=\gamma ^\mu (_\mu +\mathrm{i}gA_\mu ).$$
(131)
In the axial gauge, the gluons are gauged away, leaving behind a theory of fermions interacting via a Coulomb potential which is linearly rising in one space dimension. ’t Hooft originally solved this model for the meson spectrum, using light-cone coordinates and the light-cone gauge. This was the first practical demonstration of the power of the $`1/N`$ expansion for gauge theories.$`^\mathrm{?}`$ Working diagrammatically, he identified planar diagrams without quark loops as being of leading order in $`1/N`$ (nested rainbow graphs and ladders) and summed them up by solving (light-cone) Schwinger-Dyson and Bethe-Salpeter equations. He thus found an infinite tower of mesons lying on a kind of Regge trajectory. Early applications of this model dealt with issues of confinement, asymptotic freedom, form factors, scattering etc.$`^{\mathrm{?},\mathrm{?},\mathrm{?},\mathrm{?},\mathrm{?}}`$ For some time, the precise character of the ’t Hooft limit was somewhat disputed (strong vs. weak coupling, order of limits, ’t Hooft boundary condition, IR regularization of the quark self-energy, Wu’s alternative meson equation $`^\mathrm{?}`$). Studies in ordinary coordinates followed where the vacuum and baryons can be dealt with in a more direct way; $`^{\mathrm{?},\mathrm{?},\mathrm{?},\mathrm{?},\mathrm{?}}`$ they fully confirmed the original light-cone results for the meson spectrum. During the last few years, various extensions of the ’t Hooft model were found useful as laboratory for QCD related questions, such as adjoint quarks $`^{\mathrm{?},\mathrm{?}}`$ (a toy model for gluons and strings), heavy quarks, operator product expansion, duality etc.$`^{\mathrm{?},\mathrm{?},\mathrm{?}}`$ In line with our main subject, we propose to go through the same three topics as in Sec. 3.1 for the GN model — the vacuum, baryons and mesons. The light sector (mesons and baryons which become massless in the chiral limit) will be deferred to Sec. 3.3 where we will cover this aspect in a unified way for both models.

#### 3.2.1 Vacuum, confinement

We start with a brief reminder of the vacuum, chiral symmetry breaking and the role of confinement in the ’t Hooft model.$`^{\mathrm{?},\mathrm{?}}`$ Since this latter aspect is intimately related to infrared singularities (absent in the GN model), we work in a finite interval of length $`L`$, taking the limit $`L\mathrm{}`$ at the end whenever possible. Fermions are assumed to obey antiperiodic boundary conditions. We introduce a bilinear operator which describes color singlet dynamics to leading order in large $`N`$ (cf. GN model in Sec. 3.1.3),

$$Q(p^{},p)=\frac{1}{N}\underset{i}{}\left(\begin{array}{cc}a_i^{}(p)a_i(p^{})& b_i^{}(p)a_i(p^{})\\ a_i^{}(p)b_i(p^{})& b_i^{}(p)b_i(p^{})\end{array}\right).$$
(132)
The $`a_i,b_i`$ are annihilation operators for right- and left-handed quarks. The Hamiltonian for the ’t Hooft model in the axial gauge has the form

$`H`$ $`=`$ $`{\displaystyle \underset{p,i}{}}{\displaystyle \frac{2\pi }{L}}(p+1/2)\left(a_i^{}(p)a_i(p)b_i^{}(p)b_i(p)\right)`$
$`+{\displaystyle \frac{g^2L}{16\pi ^2}}{\displaystyle \underset{ij,n0}{}}{\displaystyle \frac{j_{ij}(n)j_{ji}(n)}{n^2}}+m{\displaystyle \underset{p,i}{}}\left(a_i^{}(p)b_i(p)+b_i^{}(p)a_i(p)\right).`$
The currents $`j_{ij}(n)`$ can be taken in the U($`N`$) form at large $`N`$,

$$j_{ij}(n)=\underset{p}{}\left(a_j^{}(p)a_i(p+n)+b_j^{}(p)b_i(p+n)\right).$$
(134)
$`H`$ is expressed in terms of $`Q`$ as

$`H`$ $`=`$ $`{\displaystyle \frac{2\pi }{L}}N{\displaystyle \underset{p}{}}\left(p+{\displaystyle \frac{1}{2}}\right)\mathrm{tr}\gamma ^5Q(p,p)+mN{\displaystyle \underset{p}{}}\mathrm{tr}\gamma ^0Q(p,p)`$
$`+{\displaystyle \frac{N^2g^2L}{16\pi ^2}}{\displaystyle \underset{n0,rs}{}}{\displaystyle \frac{1}{n^2}}\mathrm{tr}\{\delta _{rs}Q(r,r)Q(s,r)Q(rn,sn)\},`$
whereas the equation of motion for $`Q`$ becomes

$`\mathrm{i}_tQ(p^{},p)`$ $`=`$ $`{\displaystyle \frac{2\pi }{L}}\left((p^{}+1/2)\gamma ^5Q(p^{},p)(p+1/2)Q(p^{},p)\gamma ^5\right)`$
$`+m[\gamma ^0,Q(p^{},p)]+{\displaystyle \frac{Ng^2L}{8\pi ^2}}{\displaystyle \underset{n0,p^{\prime \prime }}{}}{\displaystyle \frac{1}{n^2}}\{Q(p^{},p^{\prime \prime })Q(p^{\prime \prime }n,pn)`$
$`Q(p^{}n,p^{\prime \prime }n)Q(p^{\prime \prime },p)\}.`$
Expansion of $`Q`$ around its $`c`$-number part (vacuum expectation value)

$`Q(p^{},p)`$ $`=`$ $`0|Q(p^{},p)|0+{\displaystyle \frac{1}{\sqrt{N}}}\stackrel{~}{Q}(p^{},p)+\mathrm{}`$ (137)
$`=`$ $`\delta _{p^{}p}\rho (p)+{\displaystyle \frac{1}{\sqrt{N}}}\stackrel{~}{Q}(p^{},p)+\mathrm{}`$
yields, to leading order, the HF equation

$$[h_{HF},\rho (p)]=0$$
(138)
with the single particle HF Hamiltonian

$$h_{HF}=\frac{2\pi }{L}(p+1/2)\gamma ^5+m\gamma ^0\frac{Ng^2L}{8\pi ^2}\underset{n0}{}\frac{1}{n^2}(\rho (pn)1/2).$$
(139)
In contrast to the GN model where the Hartree term dominates, here only the Fock term survives at large $`N`$, see Fig. 5

(the Coulomb interaction acts like a color exchange force). Solving the HF equation self-consistently is equivalent to summing up nested rainbow graphs in a Schwinger-Dyson equation.$`^{\mathrm{?},\mathrm{?}}`$ On general grounds, the vacuum density matrix $`\rho `$ can be parametrized as

$$\rho (p)=\frac{1}{2}+\gamma ^0\rho _0(p)\mathrm{i}\gamma ^1\rho _1(p)+\gamma ^5\rho _5(p)$$
(140)
where $`\rho _1`$ is only allowed if $`m=0`$. Charge conjugation requires $`\rho _{0,1}`$ to be even, $`\rho _5`$ to be odd functions of momenta. The Slater determinant condition $`\rho ^2(p)=\rho (p)`$ characteristic for HF then translates into

$$\rho _0^2(p)+\rho _1^2(p)+\rho _5^2(p)=\frac{1}{4}.$$
(141)
Consider first the case of massive quarks. For $`m0`$, $`\rho _1`$ vanishes and we can set

$$\rho _0=\frac{1}{2}\mathrm{sin}\theta (p),\rho _5=\frac{1}{2}\mathrm{cos}\theta (p).$$
(142)
Comparing

$$\rho (p)=v(p)v^{}(p)$$
(143)
with the HF spinors of Sec. 2.2 (Eq. 24), $`\theta (p)`$ is recognized as the Bogoliubov angle. The HF equation now becomes

$$\frac{2\pi }{L}(p+1/2)\mathrm{sin}\theta (p)m\mathrm{cos}\theta (p)+\frac{Ng^2L}{16\pi ^2}\underset{n0}{}\frac{1}{n^2}\mathrm{sin}\left(\theta (p)\theta (pn)\right)=0,$$
(144)
i.e., a nonlinear integral equation in the limit $`L\mathrm{}`$. For $`m=0`$, $`\rho _1`$ can be non-zero. As in the chiral GN model, there is now a continuum of vacua parametrized by a chiral angle $`\phi `$,

$$\left(\begin{array}{c}\rho _0(p)\\ \rho _1(p)\\ \rho _5(p)\end{array}\right)=\frac{1}{2}\left(\begin{array}{c}\mathrm{sin}\theta (p)\mathrm{cos}\phi \\ \mathrm{sin}\theta (p)\mathrm{sin}\phi \\ \mathrm{cos}\theta (p)\end{array}\right).$$
(145)
The gap equation (144) remains valid provided we set $`m=0`$. The chiral condensates are

$$\overline{q}q=\frac{N}{L}\underset{p}{}\mathrm{sin}\theta (p)\mathrm{cos}\phi ,\overline{q}\mathrm{i}\gamma ^5q=\frac{N}{L}\underset{p}{}\mathrm{sin}\theta (p)\mathrm{sin}\phi .$$
(146)
Although the gap equation (144) for the ’t Hooft model has to be solved numerically,$`^{\mathrm{?},\mathrm{?},\mathrm{?}}`$ the value of the quark condensate is known analytically, owing to an indirect determination via sum rules and the ’t Hooft equation for mesons; $`^\mathrm{?}`$ it is (for $`\phi =0`$)

$$\overline{q}q_v=\frac{N}{\sqrt{12}}\left(\frac{Ng^2}{2\pi }\right)^{1/2}.$$
(147)
(This formula has been generalized to the case of finite bare quark mass.$`^\mathrm{?}`$) In Fig. 6, we show a few typical results for the Bogoliubov angle for several values of the quark mass and compare them to the free theory. For larger quark masses, the interaction effects are very small and can be treated perturbatively.

So far, there seems to be little difference between the GN and ’t Hooft models — superficially, the vacuum looks in both cases as if the fermions had acquired a dynamical mass, although momentum dependent in the ’t Hooft case. How does confinement of quarks manifest itself? In QCD2, there are no massive quarks in the spectrum, and we have to understand how their appearance can be avoided in an independent particle picture like the HF approximation. As already suspected by ’t Hooft, the answer lies in the IR divergence of the quark self-energy.

We first remark that the gap equation (144) has no IR problem in the limit $`L\mathrm{}`$: The Coulomb interaction term has only a simple pole which can be regulated with a principal value prescription (or symmetric cut-off). The IR divergence discussed in the literature concerns the quark self-energies or, in our language, the HF single particle energies which have to be interpreted physically as removal energies (Koopmans’ theorem $`^\mathrm{?}`$). We can deduce them by varying the vacuum energy functional

$`E_{HF}/N`$ $`=`$ $`{\displaystyle \underset{p}{}}{\displaystyle \frac{2\pi }{L}}(p+1/2)\mathrm{cos}\theta (p)m{\displaystyle \underset{p}{}}\mathrm{sin}\theta (p)`$ (148)
$`+{\displaystyle \frac{Ng^2L}{32\pi ^2}}{\displaystyle \underset{n0,p}{}}{\displaystyle \frac{1}{n^2}}\left(1\mathrm{cos}(\theta (p)\theta (pn))\right)`$
with respect to the occupation number of level $`p`$ and find

$$\omega (p)=\frac{2\pi }{L}(p+1/2)\mathrm{cos}\theta (p)+m\mathrm{sin}\theta (p)+\frac{Ng^2L}{16\pi ^2}\underset{n0}{}\frac{1}{n^2}\mathrm{cos}(\theta (p)\theta (pn))$$
(149)
or, in the thermodynamic limit $`L\mathrm{}`$ (setting $`q=2\pi (p+1/2)/L`$),

$$\omega (q)=q\mathrm{cos}\theta (q)+m\mathrm{sin}\theta (q)+\frac{Ng^2}{4}\frac{\mathrm{d}q^{}}{2\pi }\frac{\mathrm{cos}\left(\theta (q)\theta (q^{})\right)}{(qq^{})^2}.$$
(150)
Due to the double pole, the integral is badly IR divergent. To exhibit the divergence, let us isolate the divergent part of the sum in Eq. (149) before taking the limit $`L\mathrm{}`$,

$$\omega (q)=q\mathrm{cos}\theta (q)+m\mathrm{sin}\theta (q)+\frac{Ng^2}{4}\frac{\mathrm{d}q^{}}{2\pi }\frac{\mathrm{cos}\left(\theta (q)\theta (q^{})\right)1}{(qq^{})^2}+\frac{Ng^2L}{48}.$$
(151)
The principal value prescription for the quadratic pole advocated by Bars and Green $`^\mathrm{?}`$ and Li et al. $`^\mathrm{?}`$ amounts to throwing away the divergent constant $`Ng^2L/48`$. As a consequence, one gets an awkward sign change (tachyonic behavior) in the single particle energies at small $`p`$ and $`m`$. The usual excuse is to say that these self-energies have no physical meaning, being gauge dependent and non-covariant. Indeed, this difficulty has no direct consequence for the vacuum properties — as pointed out above, the gap equation is free of IR divergences. Moreover, the calculation of color singlet mesons is also IR safe, as pointed out by ’t Hooft (his IR cutoff parameter $`\lambda `$ is related to our box size $`L`$ by $`\lambda =12/\pi L`$) $`^\mathrm{?}`$: The infinite self-energy term is cancelled by an equally infinite piece in the Coulomb interaction.

Since the emergence of the constant $`Ng^2L/48`$ in the single particle energy (151) but not in the vacuum energy (148) is a subtle but rather important point, let us briefly recall the argument given in Ref. 45 (see also Ref. 74). In many-body language, the Coulomb interaction produces both one- and two-body terms in $`H`$. If the single particle energies are decomposed accordingly as

$$\omega (p)=\omega ^{(1)}(p)+\omega ^{(2)}(p),$$
(152)
then the HF ground state energy is

$$0|H|0=N\underset{p}{}\left(\omega ^{(1)}(p)+\frac{1}{2}\omega ^{(2)}(p)\right).$$
(153)
As is well known, the factor 1/2 is necessary to avoid double counting of the 2-body interaction term. In the ’t Hooft model, one finds

$`\omega ^{(1)}(p)`$ $`=`$ $`{\displaystyle \frac{2\pi }{L}}\left(p+{\displaystyle \frac{1}{2}}\right)\mathrm{cos}\theta (p)+m\mathrm{sin}\theta (p){\displaystyle \frac{Ng^2L}{48}},`$
$`\omega ^{(2)}(p)`$ $`=`$ $`{\displaystyle \frac{Ng^2L}{16\pi ^2}}{\displaystyle \underset{n0}{}}{\displaystyle \frac{1}{n^2}}\left(1+\mathrm{cos}(\theta (p)\theta (p+n))\right).`$ (154)
Remembering that $`_{n0}1/n^2=\pi ^2/3`$, we see that in the ground state energy (153) the IR tamed combination $`(\mathrm{cos}(\theta (p)\theta (p+n))1)/n^2`$ appears whereas in the quark self-energy (152), the 1-term in $`\omega ^{(2)}`$ is cancelled and an IR divergence (quadratic pole for $`L\mathrm{}`$) survives.

Why bother about this IR divergence at all? As one proceeds to other applications of the ’t Hooft model than the meson spectrum, it is no longer true that the quark removal energies do not show up in any physical quantity. The prime example is finite temperature field theory in the large $`N`$ expansion by means of the standard generalization of HF to finite $`T`$.$`^\mathrm{?}`$ Here, the thermal occupation factors depend crucially on the single particle energies.$`^\mathrm{?}`$ If these are made finite by some ad-hoc prescription, quarks will manifest themselves through contributions to the pressure of order $`N`$, in conflict with confinement. We shall return to this issue in more detail in Sec. 4.3. Summarizing, we believe that the divergence of the quark single particle energies simply reflects confinement. In fact this may be the only way how confinement can be realized in the independent particle picture characteristic for the leading order in a $`1/N`$ expansion.

#### 3.2.2 Baryons

Since their structure is very different, we distinguish between those baryons which become massless in the chiral limit and those which stay massive. The first ones are closely related to Goldstone bosons and will be deferred to Sec. 3.3, where they will be discussed for the chiral GN model as well. Here, we briefly sketch what is known about the “normal” type of baryons in the ’t Hooft model.

Several papers deal with QCD2 baryons for a small number of colors, like $`N=2`$ (where they are degenerate with mesons $`^\mathrm{?}`$) or $`N=3`$.$`^\mathrm{?}`$ We concentrate on the large $`N`$ limit instead where baryons are described in relativistic HF approximation. An analytic solution like in the GN model seems to be out of the question. The work closest in spirit to our description of the GN model baryons (cf. Sec. 3.1.2) is Ref. 44. In this work, a numerical approach based on a lattice discretization of space and the use of staggered fermions has been employed. The HF equation is solved in position space, both at strong ($`Ng^2m`$) and weak ($`Ng^2m`$) coupling. One valence level is taken into account in addition to the Dirac sea. Unlike in the GN model where the filling fraction $`n/N`$ can be varied, here one has to fill the valence level completely with $`N`$ quarks in order to form a color singlet baryon. For heavy quarks, the HF equation reduces to the non-relativistic one with a linearly rising potential and the sea becomes irrelevant; the HF equation for the valence level then becomes

$$\left(\frac{1}{2m}\frac{^2}{x^2}+\frac{Ng^2}{4}dy|xy||\phi (y)|^2\right)\phi (x)=ϵ\phi (x).$$
(155)
For moderately heavy quarks, significant relativistic effects were found. The $`B=2`$ sector was also explored, with the result that two $`B=1`$ baryons were observed interacting via an exponentially decreasing repulsive potential.

#### 3.2.3 Mesons

Originally, ’t Hooft determined the meson spectrum by solving the Bethe-Salpeter equation in light-cone gauge, summing up all planar diagrams without quark loops. He found an infinite tower of mesons with masses which behave asymptotically like

$$\mu ^2(n)\pi ^2n+\mathrm{O}(\mathrm{ln}n).$$
(156)
This should be contrasted to the GN models with one or two mesons only; the difference obviously reflects confinement (Regge behavior). In the chiral limit, a massless meson appears (see Sec. 3.3). These results have subsequently been confirmed in normal coordinates, using the axial gauge. The necessary equations were first derived by Bars and Green $`^\mathrm{?}`$ and later solved numerically by several groups.$`^{\mathrm{?},\mathrm{?}}`$ The equal-time approach turns out to be significantly more involved than the light-cone approach. Nevertheless, since we cannot use light-cone quantization for all the questions which interest us and wish to stay within one common framework, we have to outline this approach here.

As before, we consider small oscillations around the HF solution in the space of Slater determinants. The quantized modes of these oscillations are just the RPA modes. The RPA equation is identical to the Bethe-Salpeter equation in a diagrammatic large $`N`$ approximation.$`^\mathrm{?}`$ To derive it, we return to the equation of motion for the bilinear quark operator $`Q(p,p^{})`$, Eq. (3.2.1), insert the expansion (137) and focus onto the next-to-leading order in $`1/\sqrt{N}`$,

$`\mathrm{i}_t\stackrel{~}{Q}(p^{},p)`$ $`=`$ $`h_{HF}(p^{})\stackrel{~}{Q}(p^{},p)\stackrel{~}{Q}(p^{},p)h_{HF}(p)`$
$`+{\displaystyle \frac{Ng^2L}{8\pi ^2}}{\displaystyle \underset{n0}{}}{\displaystyle \frac{1}{n^2}}\{\rho (p^{})\stackrel{~}{Q}(p^{}n,pn)\stackrel{~}{Q}(p^{}n,pn)\rho (p)\}.`$
Expanding $`\stackrel{~}{Q}`$ in the HF spinors, as in Sec. 3.1.3 Eq. (119) and (120), and setting

$`n,K|\stackrel{~}{Q}_{21}(p^{},p)|0`$ $`=`$ $`\delta _{K,pp^{}}\mathrm{\Phi }_+^n(K,p),`$
$`n,K|\stackrel{~}{Q}_{12}(p^{},p)|0`$ $`=`$ $`\delta _{K,pp^{}}\mathrm{\Phi }_{}^n(K,p)`$ (158)
then yields the two coupled equations

$`\pm E(K)\mathrm{\Phi }_\pm (K,p)`$ $`=`$ $`(\omega _0(pK)+\omega _0(p))\mathrm{\Phi }_\pm (K,p)`$ (159)
$`{\displaystyle \frac{Ng^2L}{8\pi ^2}}{\displaystyle \underset{p^{}(p)}{}}{\displaystyle \frac{1}{(pp^{})^2}}\{s(p,p^{},K)\mathrm{\Phi }_{}(K,p^{})`$
$`+c(p,p^{},K)\mathrm{\Phi }_\pm (K,p^{})\}`$
with the definitions

$`s(p,p^{},K)`$ $`=`$ $`\mathrm{sin}\left({\displaystyle \frac{\theta (p^{})\theta (p)}{2}}\right)\mathrm{sin}\left({\displaystyle \frac{\theta (pK)\theta (p^{}K)}{2}}\right),`$
$`c(p,p^{},K)`$ $`=`$ $`\mathrm{cos}\left({\displaystyle \frac{\theta (p^{})\theta (p)}{2}}\right)\mathrm{cos}\left({\displaystyle \frac{\theta (pK)\theta (p^{}K)}{2}}\right).`$
Eqs. (159) have the standard RPA form of Eq. (125), with

$`A_{pp^{}}(K)`$ $`=`$ $`\delta _{pp^{}}(\omega _0(pK)+\omega _0(p)){\displaystyle \frac{Ng^2L}{8\pi ^2}}{\displaystyle \frac{(1\delta _{pp^{}})}{(pp^{})^2}}c(p,p^{},K),`$
$`B_{pp^{}}(K)`$ $`=`$ $`{\displaystyle \frac{Ng^2L}{8\pi ^2}}{\displaystyle \frac{(1\delta _{pp^{}})}{(pp^{})^2}}s(p,p^{},K).`$ (161)
The diagrammatic content of these equations is sketched in Fig. 5 (see Sec. 3.2.1). Using a slightly different language, Bars and Green $`^\mathrm{?}`$ interpret them in terms of forward and backward going strings. The covariance of the spectrum is a non-trivial issue if one works canonically in normal coordinates, but can be checked either numerically or via the $`1/N`$ expansion of the Poincaré algebra. Meson wavefunctions (or rather RPA amplitudes) can be found in the literature for a few low-lying states.$`^\mathrm{?}`$ In the limit of heavy quarks, the RPA reduces to the non-relativistic Schrödinger equation with a linear potential. The light-cone TDA (’t Hooft equation) on the other hand can be recovered by going to the infinite momentum frame $`^\mathrm{?}`$ (backward going bubbles are suppressed, see also Sec. 3.3).

### 3.3 Chiral limit and Goldstone modes

Since QCD2 or the chiral GN model in the massless limit break chiral symmetry spontaneously (at least at infinite $`N`$), the appearance of a boson whose mass behaves like the square root of the bare quark mass has to be expected. This has a lot in common with four dimensional theories and makes these models good toy models for QCD. However, the chiral limit is somewhat obscured by a certain pathology of two dimensional theories which gives rise to a much greater variety of (decoupled) massless states. In particular, non-interacting massless mesons and baryons already appear at finite $`N`$ where the Coleman theorem certainly forbids talking about Goldstone bosons. Our main interest here is in the large $`N`$ limit where we will derive both massless mesons and baryons. However, before doing so, we also briefly explain the more pathological aspects of two dimensional theories leading to massless states at finite $`N`$ as well.

The identification of massless excitations in two dimensional field theories is one of the tasks where light-cone quantization is particularly elegant. Due to the light-cone dispersion relation $`p_+=\frac{m^2}{2p_{}}`$, a massless state in two dimensions has zero energy at all momenta. Therefore its presence can be associated with a local symmetry of the Hamiltonian. In the chiral limit the right-handed fermions drop out of the light-cone Hamiltonian $`P^+=HP`$.$`^{\mathrm{?},\mathrm{?},\mathrm{?}}`$ Since the Coulomb interaction admits local chiral transformations, $`P^+`$ is invariant under local phase transformations of right-handed quarks. Free, massless bosons then must appear and it is sufficient to identify the winding number of their field with baryon number to infer the existence of free, massless baryons as well. The IR trouble of Coleman’s theorem is avoided by the fact that the massless modes decouple.

The essence of the argument can of course also be phrased without recourse to light-cone quantization.$`^\mathrm{?}`$ Let us return to the fermion currents discussed in Sec. 2.2. Vector current conservation allows us to set

$$j_V^\mu =ϵ^{\mu \nu }_\nu \varphi $$
(162)
with a scalar field $`\varphi `$. If the axial current is non-anomalous, we get in the limit $`m0`$ (see $`j_A^\mu `$ above, Eq. (33))

$$_\mu j_A^\mu =\mathrm{}\varphi =0,$$
(163)
i.e., the announced massless boson. This purely classical argument can easily be promoted to the quantum level within the canonical framework.$`^\mathrm{?}`$ Define the right-handed current as

$$j_R(q,t)=\frac{1}{2}(j_V^0(q,t)+j_V^1(q,t))$$
(164)
where Heisenberg operators, Fourier transformed with respect to the space coordinate $`x`$ are used. In the chiral limit, vector and axial vector current conservation imply

$$[j_R(q,t),H]=qj_R(q,t).$$
(165)
Besides, if $`P`$ denotes the momentum operator, we have trivially

$$[j_R(q,t),P]=qj_R(q,t).$$
(166)
For any energy-momentum eigenstate $`|i`$ with eigenvalues ($`E_i,P_i`$), we thus get

$`Hj_R(q,t)|i`$ $`=`$ $`(E_iq)j_R(q,t)|i,`$
$`Pj_R(q,t)|i`$ $`=`$ $`(P_iq)j_R(q,t)|i.`$ (167)
Thereby decoupled massless bosons can be added to or removed from any stationary state. Similarly, left moving massless bosons can be related to the left-handed current.

Let us now turn to the large $`N`$ limit and discuss the emergence of massless particles there. Again we benefit from the experience in many-body theory. The Goldstone theorem has a close correspondence which we may invoke here (remember that it also holds in non-relativistic theories): If the HF ground state breaks a continuous symmetry spontaneously, the RPA develops a gapless mode. The basic argument can be given in two lines, using the language of linear response theory (without external force).$`^\mathrm{?}`$ Take the equation of motion for the one-body density matrix,

$$\mathrm{i}\dot{\rho }=[h(\rho ),\rho ].$$
(168)
Expand $`\rho `$ around its ground state expectation value, $`\rho =\rho ^{(0)}+\delta \rho `$, and linearize

$$\mathrm{i}\delta \dot{\rho }=[h(\rho ^{(0)}),\delta \rho ]+[\frac{\delta h}{\delta \rho }\delta \rho ,\rho ^{(0)}].$$
(169)
This equation has the form of the RPA, cf. the above derivation in the equations of motion approach in Sec. 3.1.3. If the ground state density matrix $`\rho `$ breaks a continuous symmetry, we find that the deformation $`\delta \rho `$ corresponding to the symmetry transformation (e.g. a chiral rotation) solves the RPA equation with zero energy. The backward going ph-bubbles are crucial for getting the zero modes; as is well known, the TDA does not have this property.

These insights enable us to construct the Goldstone boson type solution of the RPA equation explicitly. In the case of QCD2, the RPA amplitudes of the “pion” can easily be obtained by projecting the right-handed current onto HF spinors. Since these depend only on the Bogoliubov angle, the pion amplitudes are intimately related to the broken symmetry vacuum. One finds $`^{\mathrm{?},\mathrm{?}}`$

$`\mathrm{\Phi }_+^\pi (K,p)`$ $`=`$ $`\sqrt{{\displaystyle \frac{2\pi }{K}}}\mathrm{cos}\left(\theta (p)/2\right)\mathrm{sin}\left(\theta (pK)/2\right),`$
$`\mathrm{\Phi }_{}^\pi (K,p)`$ $`=`$ $`\sqrt{{\displaystyle \frac{2\pi }{K}}}\mathrm{sin}\left(\theta (p)/2\right)\mathrm{cos}\left(\theta (pK)/2\right).`$ (170)
The chiral GN model can be analyzed in a similar manner. Generalizing the derivation of the RPA equation, Sec. 3.1.3, to the chirally invariant model, one finds an integral equation with a two-term separable kernel (instead of one-term in the model with discrete chiral symmetry) and correspondingly two bound states, the threshold bound state at 2$`m`$ (the “$`\sigma `$”) and a massless one (the “$`\pi `$”). The pion has the same structure as in Eq. (170), except that the Bogoliubov angles now refer to the GN vacuum. They are just the Bogoliubov angles of the free, massive fermions given in Eq. (25). In this particular case, we are in the rare position of having complete analytical control over a meson (the $`\pi `$) in an arbitrary Lorentz frame.

In Fig. 7, we illustrate the disappearance of the backward components as one boosts the pion to higher momenta (in the rest frame, both are exactly equal). At very high momenta, the constant wave function familiar from the light-cone meson equation emerges.$`^\mathrm{?}`$ In the infinite momentum limit, the ’t Hooft and GN model pions become identical, although their “wave functions” (170) differ in any other frame due to different Bogoliubov angles. This illustrates nicely how the vacuum decouples in the infinite momentum frame, one of the celebrated advantages of the light-cone approach (triviality of the vacuum). Connected to this is the disappearance of the backward amplitudes — the RPA reduces to the TDA in the infinite momentum frame, the other characteristic benefit on the light-cone.

We now turn to the issue of baryons which become massless in the chiral limit, a phenomenon which does not happen in 3+1 dimensions. The best way of thinking about them is in terms of coherent, topological excitations of the pion field. Like the massive baryons discussed in Secs. 3.1.2 and 3.2.2, they are accessible through the relativistic HF approximation. In order to find these solutions in the first place, another point of view has turned out to be very useful. It is closely related to bosonization,$`^\mathrm{?}`$ but since only the chiral field is bosonized, it is conceptually and technically much simpler. It is also similar in spirit to the Skyrme picture of baryons in 3+1 dimensions.$`^\mathrm{?}`$

Following Salcedo et al.,$`^\mathrm{?}`$ we start from a variational ansatz for the baryon near the chiral limit. The interaction part of the ground state energy functional

$`E[\rho ]`$ $`=`$ $`N{\displaystyle dx\mathrm{tr}\left[\left(\mathrm{i}\gamma ^5\stackrel{}{}_x+m\gamma ^0\right)\rho (x,x)\right]}`$ (171)
$`+{\displaystyle \frac{N^2g^2}{8}}{\displaystyle dxdy|xy|\mathrm{tr}\left(\rho (x,y)\rho (y,x)\right)}`$
(where $`\stackrel{}{}_x`$ acts only on the first variable of $`\rho (x,x)`$) is invariant under any local chiral transformation. Let us perform a local chiral rotation of the vacuum density matrix,

$$\rho ^{}(x,y)=\mathrm{e}^{\mathrm{i}\chi (x)\gamma ^5}\rho ^v(xy)\mathrm{e}^{\mathrm{i}\chi (y)\gamma ^5}.$$
(172)
We can use the following general decomposition of the baryon density matrices $`\rho ^v`$ and $`\rho ^{}`$

$$\rho =\frac{1}{2}\rho _e+\gamma ^0\rho _0\mathrm{i}\gamma ^1\rho _1+\gamma ^5\rho _5$$
(173)
with corresponding superscripts. In the vacuum, $`\rho _e=\delta (xy)`$, due to translational invariance all components depend only on $`xy`$ and $`\rho _1`$ vanishes unless $`m=0`$. If the vacuum breaks chiral symmetry, the ansatz (172) yields a new density matrix which can be inserted into Eq. (171); $`\chi (x)`$ is then determined by minimizing $`E`$. In order to understand the emergence of baryon number out of winding number, we need the short distance singularities of $`\rho ^v`$ which are the same as in the free theory,

$`\rho _0(x)`$ $`=`$ $`{\displaystyle \frac{m}{2\pi }}K_0(m|x|){\displaystyle \frac{m}{2\pi }}\mathrm{ln}(m|x|),`$
$`\rho _5(x)`$ $`=`$ $`{\displaystyle \frac{m}{2\pi \mathrm{i}}}\mathrm{sgn}(x)K_1(m|x|){\displaystyle \frac{1}{2\pi \mathrm{i}x}}.`$ (174)
Using

$`\rho _e^{}(x,y)`$ $`=`$ $`\mathrm{cos}(\chi (x)\chi (y))\rho _e^v(xy)+2\mathrm{i}\mathrm{sin}(\chi (x)\chi (y))\rho _5^v(xy),`$
$`\rho _5^{}(x,y)`$ $`=`$ $`\mathrm{cos}(\chi (x)\chi (y))\rho _5^v(xy)+{\displaystyle \frac{\mathrm{i}}{2}}\mathrm{sin}(\chi (x)\chi (y))\rho _e^v(xy),`$
one then finds

$`\underset{yx}{lim}[\rho _e^{}(x,y)\rho _e^v(xy)]`$ $`=`$ $`{\displaystyle \frac{1}{\pi }}_x\chi (x),`$
$`\underset{yx}{lim}_x[\rho _5^{}(y,x)\rho _5^v(xy)]`$ $`=`$ $`{\displaystyle \frac{\mathrm{i}}{4\pi }}(_x\chi (x))^2.`$ (176)
Similarly,

$$\rho _0^{}(x,x)\mathrm{i}\rho _1^{}(x,x)=\frac{\overline{q}q_v}{2N}\mathrm{e}^{\mathrm{i2}\chi (x)}$$
(177)
where $`\overline{q}q_v`$ refers to the vacuum. If we invoke a finite spatial box of length $`L`$ for these topological considerations, the baryon number gets identified with the winding number of the chiral field,

$$B=_0^Ldx[\rho _e^{}(x,x)\rho _e^v(0)]=\frac{1}{\pi }(\chi (L)\chi (0)).$$
(178)
This is an integer if the baryon goes over into the vacuum at $`0`$ and $`L`$. For small quark mass, find

$$E[\rho ^{}]=E[\rho ^v]+Ndx\left\{\frac{1}{2}(_x\varphi )^2+m\frac{\overline{q}q_v}{N}\left[\mathrm{cos}(\sqrt{4\pi }\varphi (x))1\right]\right\}$$
(179)
with the rescaled field

$$\varphi (x)=\frac{1}{\sqrt{\pi }}\chi (x)$$
(180)
and $`\overline{q}q_v`$ referring to the vacuum in the chiral limit. $`E[\rho ^{}]`$ is minimized if $`\varphi `$ satisfies the sine-Gordon equation $`^\mathrm{?}`$

$$_x^2\varphi +\sqrt{4\pi }m\frac{\overline{q}q_v}{N}\mathrm{sin}(\sqrt{4\pi }\varphi (x))=0.$$
(181)
The solution with baryon number one is just the famous sine-Gordon kink,

$$\varphi (x)=\frac{2}{\sqrt{\pi }}\mathrm{arctan}\mathrm{exp}\left(\frac{x}{x_B}\right).$$
(182)
It is localized with spatial extent $`x_B`$,

$$x_B=\left(4\pi m\frac{\overline{q}q_v}{N}\right)^{1/2}$$
(183)
and has mass

$$M_B=8N\left(\frac{m}{4\pi }\frac{\overline{q}q_v}{N}\right)^{1/2}.$$
(184)
If one linearizes the sine-Gordon equation (in the limit $`m0`$), one can immediately read off the Gell-Mann, Oakes, Renner relation (GOR) $`^\mathrm{?}`$

$$m_\pi ^2=4\pi m\frac{\overline{q}q_v}{N}$$
(185)
for the Goldstone boson. We have thus derived an effective low energy theory near the chiral limit. Since the potential energy does not enter Eq. (179) but only the kinetic energy and mass term, the resulting low energy effective theory is universal for a whole class of such chiral models.

In the limit $`m0`$, both the pion and the baryon become massless; the baryon gets completely delocalized. This could actually have been anticipated from current conservation: $`_1j_A^1=0`$ holds for any stationary state and in two dimensions, $`j_A^1`$ is just the baryon density $`j_V^0`$ (cf. Sec. 2.2).

An attempt to illustrate the baryon structure in the chiral limit has been made in Fig. 8. For $`m=0`$, $`\chi (x)`$ becomes simply a linear function and the scalar and pseudoscalar condensates

$$\overline{q}q=\overline{q}q_v\mathrm{cos}(2\pi (xx_0)/L),\overline{q}\mathrm{i}\gamma ^5q=\overline{q}q_v\mathrm{sin}(2\pi (xx_0)/L),$$
(186)
can be regarded as projections of a spiral of radius $`\overline{q}q_v`$ winding around once. This picture will be generalized later on and turn out to be useful for characterizing baryonic matter in two dimensional chiral models.

In Ref. 44 the Skyrme type picture of the baryon was compared to the exact numerical solution for light quarks and good agreement was found also for small, finite quark masses. It was pointed out that these results support the Skyrme picture of the baryon more than the bag picture in two dimensions.

Finally, let us come back to the question of validity of the variational calculation based on the ansatz (172). In the massless limit, it is easy to convince oneself that Eq. (172) with a linear function $`\chi (x)`$ solves the HF equation exactly. Thus for instance in the ’t Hooft model, the massless HF equation reads $`^\mathrm{?}`$

$$\mathrm{i}(\gamma _5)_{\alpha \beta }\frac{}{x}\phi _\beta ^{(n)}(x)+\frac{Ng^2}{4}dy|xy|\rho _{\alpha \beta }^{}(x,y)\phi _\beta ^{(n)}(y)=\omega _n\phi _\alpha ^{(n)}(x).$$
(187)
Upon substituting

$$\phi _\alpha ^{(n)}(x)=\left(\mathrm{e}^{\mathrm{i}\pi x\gamma _5/L}\right)_{\alpha \beta }\stackrel{~}{\phi }_\beta ^{(n)}(x)$$
(188)
as we are instructed to do by Eq. (172), we discover that $`\stackrel{~}{\phi }^{(n)}`$ does indeed solve the HF equation, the only change being that the single particle energy $`\omega _n`$ gets replaced by $`\omega _n+\pi /L`$. The same argument goes through in the chiral GN model, or in any field theory where the interaction term has a local chiral invariance for that matter. This shows that the result becomes exact in the chiral limit (to leading order in the $`1/N`$ expansion, of course).

## 4 Finite temperature

### 4.1 Finite temperature versus finite extension

We briefly discuss the relation between field theories at finite extension and finite temperature.$`^\mathrm{?}`$ Equivalence between relativistic theories at finite extension and finite temperature is formally almost trivial but gives rise to some of the most intriguing consequences of covariance. By rotational invariance in the Euclidean, the value of the partition function of a system with finite extension $`L`$ in one space direction and $`\beta `$ in the time direction is invariant under the exchange of these two parameters,

$$Z(\beta ,L)=Z(L,\beta ),$$
(189)
provided bosonic (fermionic) fields satisfy periodic (antiperiodic) boundary conditions in both time and compact space coordinate. Sending one of the parameters ($`\beta ,L`$) to infinity, relativistic covariance connects the thermodynamic properties of a canonical ensemble with the properties of the pure state of the vacuum corresponding to the same physical system but at finite extension. In particular, as a consequence of (189), energy density and pressure are related by

$$ϵ(\beta ,L)=p(L,\beta ).$$
(190)
For a system of non-interacting particles this equation can be used to relate quantitatively the Stefan-Boltzmann law with the Casimir effect — phenomena which at first glance would not seem to be related at all.$`^\mathrm{?}`$ Another application of this discrete relativistic symmetry worth mentioning are quark propagators and the interpretation of lattice data in “funny space”.$`^\mathrm{?}`$ In four dimensional QCD, the confinement-deconfinement phase transition and the chiral phase transition, when quarks are present, appear as “quantum phase transitions” if one swaps space and Euclidean time. They are driven by changes in quantum rather than thermal fluctuations which in turn are induced by changes of a parameter of the system ($`L`$).$`^{\mathrm{?},\mathrm{?}}`$ We have found it useful to apply this technique to the two dimensional models at hand to see how it works in detail in the presence of interaction effects. Consistency between finite extension and thermal field theory calculations provides a non-trivial test of the approximations used and will help us to clarify the role of confinement in the ’t Hooft model. For the case of two dimensions with one extension taken to be infinite, the symmetry is illustrated in Fig. 9.

To leading order in the large $`N`$ limit, it is easy to include finite temperature and chemical potential since the required generalization of our main tool, the HF approximation, is well known.$`^\mathrm{?}`$ Let us focus on finite temperature here since finite density and chemical potential will be the subject of Secs. 5 and 6. The main difference to the temperature zero case is that for sytems in thermal equilibrium, one is no longer dealing with a pure state but with a mixed state. The concepts of independent particles and single particle states remain valid. Whereas at $`T=0`$ each orbit is either empty or occupied (sharp Fermi surface), at finite $`T`$ the occupation number is given by the thermal occupation probability (Dirac-Fermi distribution $`^\mathrm{?}`$). The one-body density matrix no longer satisfies the projector property $`\rho ^2=\rho `$ characteristic for a single Slater determinant but assumes the form

$$\rho =N_+uu^{}+N_{}vv^{},$$
(191)
with occupation numbers

$$N_\pm =\frac{1}{\mathrm{e}^{\pm \beta ϵ}+1}.$$
(192)
Using

$$N_++N_{}=1,N_+N_{}=\mathrm{tanh}\frac{\beta ϵ}{2},$$
(193)
one can write

$$\rho =\frac{1}{2}+\left(\gamma ^0\rho _0+\gamma ^5\rho _5\right)\mathrm{tanh}\frac{\beta ϵ}{2}.$$
(194)
where $`\rho _0,\rho _5`$ are given in terms of Bogoliubov angles in the same way as at zero temperature (in the chiral limit, there could be a $`\rho _1`$-term in addition.) In contrast to the $`T=0`$ case now the single particle energies $`ϵ`$ become physically relevant and will appear explicitly in the self-consistency equation, simply because excited states do play a role in the thermal equilibrium. The alternative scenario, namely HF at $`T=0`$ but finite extension, does not require any new tools, at least for the GN model. In the case of the ’t Hooft model, the situation will turn out to be more involved due to the different gluon dynamics in a finite box (cf. Sec. 4.3).

### 4.2 GN model on a cylinder

Before embarking on the details of the finite temperature GN model, we have to come back once more to the issue of SSB in lower dimensional systems.$`^\mathrm{?}`$ Since the thermodynamic limit can be taken only with respect to one dimension now, there are further restrictions as compared to those of Sec. 2.4. According to the Mermin-Wagner theorem, even a discrete symmetry should be immediately restored at $`T=0^+`$. The restoration of chiral symmetry at finite temperature is expected to be driven by the presence of kinks and antikinks $`^{\mathrm{?},\mathrm{?}}`$ for any large but finite $`N`$. If however the limit $`N\mathrm{}`$ is taken before the thermodynamic limit, these configurations are suppressed and one is left with the mean field theory where a second order phase transition at a critical temperature $`T_c0`$ becomes possible. Similarly, for the continuous model, for large but finite $`N`$ at $`T0`$, the almost long range order disappears. Only if $`N\mathrm{}`$ is taken first, mean field theory applies and one finds the same critical temperature as in the discrete model. The non-analytic structure revealed by the phase diagram of the infinite $`N`$ GN model depends on the diverging number of fermion components like in reduced models.$`^{\mathrm{?},\mathrm{?}}`$

Barducci et al. $`^\mathrm{?}`$ advocate to include a small bare fermion mass, thereby eliminating kink-antikink configurations in the thermodynamic limit. Since their results go over rather smoothly into those obtained in mean field theory at $`m=0`$, we see nothing wrong in considering the chiral limit directly for some questions.

First consider finite extension which requires only minimal modifications of the HF approach for the vacuum discussed in Sec. 3.1.1. Since the fermion momenta are discretized due to the antiperiodic boundary conditions, the gap equation (77) (self-consistency requirement for the scalar density) has to be replaced by

$$m=\frac{2Ng^2}{L}\underset{n=0}{\overset{N_\mathrm{\Lambda }}{}}\frac{m}{\sqrt{m^2+k_n^2}}$$
(195)
with

$$N_\mathrm{\Lambda }=\frac{L\mathrm{\Lambda }}{4\pi },k_n=\frac{\pi }{L}(2n+1).$$
(196)
As the finite extension does not affect any UV properties, we can renormalize the theory in the limit $`L\mathrm{}`$. Denoting the physical mass on the infinite line by $`m_0`$ and using Eqs. (77) and (195), we obtain

$$\underset{\mathrm{\Lambda }\mathrm{}}{lim}\left(_0^{\mathrm{\Lambda }/2}\frac{\mathrm{d}k}{2\pi }\frac{1}{\sqrt{m_0^2+k^2}}\frac{1}{L}\underset{n=0}{\overset{N_\mathrm{\Lambda }}{}}\frac{1}{\sqrt{m^2+k_n^2}}\right)=0.$$
(197)
This equation is now free of divergences but cannot be solved analytically in general. To find the critical length below which chiral symmetry gets restored it is sufficient to set $`m=0`$ and solve for $`L`$,

$$\frac{1}{2}\mathrm{ln}\frac{\mathrm{\Lambda }}{m_0}=\underset{n=0}{\overset{N_\mathrm{\Lambda }}{}}\frac{1}{2n+1}.$$
(198)
The sum can be performed,

$$\underset{n=0}{\overset{N_\mathrm{\Lambda }}{}}\frac{1}{2n+1}=\frac{1}{2}(\mathrm{C}+\mathrm{ln}4N_\mathrm{\Lambda })+O(1/N_\mathrm{\Lambda }^2)$$
(199)
($`\mathrm{C}0.5772`$ Euler constant) yielding the critical length

$$L_c=\frac{\pi }{m_0}e^\mathrm{C}.$$
(200)
A bare fermion mass $`m_b`$ can easily be taken into account if one replaces the gap equation (197) by

$$\left(\frac{m_b}{2Ng^2}\right)\left(\frac{1}{m_0}\frac{1}{m}\right)+\underset{\mathrm{\Lambda }\mathrm{}}{lim}(_0^{\mathrm{\Lambda }/2}\frac{\mathrm{d}k}{2\pi }\frac{1}{\sqrt{m_0^2+k^2}}\frac{1}{L}\underset{n=0}{\overset{N_\mathrm{\Lambda }}{}}\frac{1}{\sqrt{m^2+k_n^2}})=0.$$
(201)
Note the appearance of the new parameter $`\frac{m_b}{Ng^2}`$ in addition to $`m_0`$. Numerical solutions of this gap equation for different quark masses are shown in Fig. 10, exhibiting a second order phase transition in the chiral limit and a cross-over phenomenon for non-zero bare masses. Since it explicitly violates chiral symmetry, the bare quark mass plays a role analogous to that of an external magnetic field in the ferromagnetic phase transition. Finally, we reemphasize that to this order in the $`1/N`$ expansion the calculation is identical for the GN models with discrete or continuous chiral symmetry.

We now switch to the equivalent finite temperature calculation to familiarize ourselves with the corresponding generalization of the HF approximation. As indicated in Sec. 4.1, the gap equation of thermal HF has to be modified by $`T`$-dependent occupation numbers which “smear out” the Fermi surface,

$$mm_b=Ng^2\frac{\mathrm{d}k}{2\pi }\frac{m}{ϵ(k)}\mathrm{tanh}\frac{\beta ϵ(k)}{2},ϵ(k)=\sqrt{m^2+k^2}.$$
(202)
The $`\mathrm{tanh}`$ was derived in Eq. (194). Renormalizing at $`T=0`$ then yields the following equation for the $`T`$-dependence of the physical fermion mass,

$$0=m_b\left(\frac{1}{m_0}\frac{1}{m}\right)+Ng^2\frac{\mathrm{d}k}{2\pi }\left(\frac{\mathrm{tanh}\frac{1}{2}\beta \sqrt{m^2+k^2}}{\sqrt{m^2+k^2}}\frac{1}{\sqrt{m_0^2+k^2}}\right).$$
(203)
The integral in (203) converges and we recover the well known results obtained by functional integration methods,$`^{\mathrm{?},\mathrm{?}}`$ in particular a second order phase transition at a critical temperature $`T_c`$ which coincides with $`L_c^1`$ from Eq. (200). Bosonic excitations are $`1/N`$ suppressed and play no role at this stage. Hence just like in the finite extension case, there is no distinction yet between GN models with continuous and discrete chiral symmetry.

It is interesting to compare the two approaches leading to the gap equations (201) and (203) respectively at an intermediate stage to understand how covariance “works”. The crucial identity which ensures the equivalence of the finite $`L`$ and finite $`T`$ calculations is

$$_{\mathrm{\Lambda }/2}^{\mathrm{\Lambda }/2}\frac{\mathrm{d}k}{\sqrt{m^2+k^2}}\mathrm{tanh}\left(\frac{1}{2}\beta \sqrt{m^2+k^2}\right)\stackrel{!}{=}\frac{4\pi }{L}\underset{n=0}{\overset{N_\mathrm{\Lambda }}{}}\frac{1}{\sqrt{m^2+k_n^2}}$$
(204)
for $`\beta =L`$. Eq. (204) can be verified with the help of Cauchy’s theorem. The integrand has simple poles in the upper half plane at

$$k=\mathrm{i}\sqrt{m^2+((2n+1)\pi /\beta )^2},(n=0\mathrm{}N_\mathrm{\Lambda })$$
(205)
with residues $`2\mathrm{i}/(\beta \sqrt{m^2+((2n+1)\pi /\beta )^2})`$. Thereby, a continuous average involving thermal occupation numbers can be converted into a sum characteristic for a finite interval. This is of course reminiscent of the imaginary time approach with Matsubara frequencies although different in detail (in the canonical approach) due to the interchange of frequencies and momenta. In the Euclidean path integral formalism, we would expect no difference whatsoever.

So far we have concentrated on the dynamical mass or equivalently the chiral condensate. Other bulk thermodynamic observables can easily be calculated from the thermodynamic potential. Let us therefore consider the two functions which are minimized with respect to $`m`$ in the two approaches. At finite extension in the chiral limit, it is the HF vacuum energy density

$$\frac{}{N}=\frac{m^2}{2Ng^2}\frac{1}{L}\underset{n}{}\sqrt{m^2+k_n^2}$$
(206)
(cf. Sec. 3.1.1). The gap equation (195) follows from the condition $`/m=0`$. In the finite $`T`$ case on the other hand, we have to minimize the thermodynamic potential $`\varphi `$ or free energy density $``$ related to the partition function via

$$\frac{\varphi }{L}==\frac{1}{\beta L}\mathrm{ln}\mathrm{tr}\mathrm{e}^{\beta H}$$
(207)
and equal to the negative of the pressure. In the GN model at the HF level, the Hamiltonian is given by that of the free, massive theory plus a $`c`$-number term correcting for double counting of the two-body interaction,

$$H=H_m^{(0)}+\frac{Lm^2}{2g^2}.$$
(208)
Therefore the partition function can be evaluated following the textbook approach for a free Fermi gas,$`^\mathrm{?}`$

$$\frac{}{N}=\frac{m^2}{2Ng^2}\frac{\mathrm{d}k}{2\pi }\left[ϵ(k)+\frac{2}{\beta }\mathrm{ln}\left(1+\mathrm{e}^{\beta ϵ(k)}\right)\right].$$
(209)
Differentiating with respect to $`m`$ reproduces Eq. (203) (here for $`m_b=0`$). Since we have already verified that the two gap equations are identical, it is sufficient to check that $``$ and $``$ coincide for one value of $`m`$ (upon identifying $`L`$ and $`\beta `$) to establish that the two functions agree everywhere as expected on general grounds (cf. Eq. (190)). We choose $`m=0`$ where all calculations can be done in closed form. For the energy density in the finite interval we find, using heat kernel regularization (as usually done in treating the Casimir effect),

$$\frac{}{N}=\frac{\pi }{6L^2}.$$
(210)
For finite temperature Eq. (209) yields

$$\frac{}{N}=\frac{\pi }{6\beta ^2}.$$
(211)
i.e., the Stefan-Boltzmann law. In both cases, the same $`L`$\- resp. $`\beta `$-independent quadratically divergent term has been dropped. This illustrates the above mentioned correspondence between the Casimir effect and the Stefan-Boltzmann law and completes the demonstration that compressed and hot relativistic systems are fully related by covariance for the case of the GN model.

In order to derive other bulk thermodynamic observables it is useful to renormalize the thermodynamic potential.$`^{\mathrm{?},\mathrm{?}}`$ Using the gap equation at zero temperature one can write

$$\frac{}{N}|_{\mathrm{ren}}=\frac{m^2}{4\pi }\left(\mathrm{ln}\left(\frac{m^2}{m_0^2}\right)1\right)+\frac{m_0^2}{4\pi }\frac{2}{\beta }\frac{\mathrm{d}k}{2\pi }\mathrm{ln}\left(1+\mathrm{e}^{\beta ϵ(k)}\right).$$
(212)
Following Barducci et al. $`^\mathrm{?}`$ we have added a constant “bag term” $`m_0^2/4\pi `$ in order to to normalize the (unobservable) vacuum pressure to zero. The derivative $`_{\mathrm{ren}}/m`$ reproduces the renormalized gap equation (203). From this, the bulk thermodynamical observables can be obtained. For temperatures below the phase transition, pressure and internal energy density are

$`{\displaystyle \frac{p}{N}}`$ $`=`$ $`{\displaystyle \frac{m^2}{4\pi }}\left(\mathrm{ln}\left({\displaystyle \frac{m^2}{m_0^2}}\right)1\right){\displaystyle \frac{m_0^2}{4\pi }}+{\displaystyle \frac{2}{\beta \pi }}{\displaystyle _0^{\mathrm{}}}dk\mathrm{ln}\left(1+\mathrm{e}^{\beta ϵ(k)}\right)`$ (213)
$`{\displaystyle \frac{\epsilon }{N}}`$ $`=`$ $`{\displaystyle \frac{m^2}{4\pi }}\left(\mathrm{ln}\left({\displaystyle \frac{m^2}{m_0^2}}\right)1\right)+{\displaystyle \frac{m_0^2}{4\pi }}+{\displaystyle \frac{2}{\pi }}{\displaystyle _0^{\mathrm{}}}dk{\displaystyle \frac{ϵ(k)}{1+\mathrm{e}^{\beta ϵ(k)}}}`$ (214)
where $`m`$ has to be evaluated at the minimum of $``$ (solution of the gap equation). In the region above $`T_c`$, the results are simply

$`{\displaystyle \frac{p}{N}}`$ $`=`$ $`{\displaystyle \frac{\pi T^2}{6}}{\displaystyle \frac{m_0^2}{4\pi }},`$ (215)
$`{\displaystyle \frac{\epsilon }{N}}`$ $`=`$ $`{\displaystyle \frac{\pi T^2}{6}}+{\displaystyle \frac{m_0^2}{4\pi }}.`$ (216)
In Fig. 11, pressure and internal energy density have been plotted after dividing out the $`T^2`$ dependence of the Stefan-Boltzmann law, so as to enhance interaction effects.

In conclusion, both finite extension and finite temperature calculations agree perfectly for the GN model. Covariance of the single particle energies is certainly the key feature behind this result.

### 4.3 ’t Hooft model on a cylinder, decompactification

The behavior of the ’t Hooft model at finite temperature is rather controversial. McLerran and Sen $`^\mathrm{?}`$ argue that there is no deconfining phase transition, except possibly at infinite temperature. Ming Li,$`^\mathrm{?}`$ using standard finite temperature field theory methods, concludes that chiral symmetry may get restored in the limit $`T\mathrm{}`$. In both of these studies, severe infrared problems were encountered, either in the form of divergent diagrams or ambiguous quark self energies. Hansson and Zahed $`^\mathrm{?}`$ address the question of the $`N`$-dependence of thermodynamic quantities. Several studies have considered QCD2 on a spatial circle at zero temperature. As explained above, this can also be reinterpreted as finite temperature calculations for a spatially extended system. Lenz et al. $`^\mathrm{?}`$ observe a chiral phase transition in the massless ’t Hooft model at some critical length, strongly reminiscent of the GN model. Dhar et al. $`^{\mathrm{?},\mathrm{?},\mathrm{?}}`$ treat the zero mode gluons in a more ambitious way than Ref. 45 using technology developed in the framework of matrix models and string theory, but are not able to fully solve the resulting complicated equations. They propose that the gauge variables get decompactified by the fermions in complete analogy with the Schwinger model,$`^\mathrm{?}`$ a claim which has recently been disputed by Engelhardt.$`^\mathrm{?}`$

Let us first consider finite temperature HF. Here, the single particle energies enter through the equilibrium thermal occupation numbers. If we used the principal value prescription with tachyonic behavior at small momenta $`^\mathrm{?}`$ (Sec. 3.2.1), we would necessarily get a non-trivial temperature dependence in leading order of the $`1/N`$ expansion ($`T`$-dependent condensate and pressure of order $`N`$). Physically, this would imply that the quarks get deconfined. Color singlet mesons which might contribute to thermodynamic observables can show up only in next-to-leading order in $`1/N`$. We have argued before that the HF single particle energies (interpreted as removal energies) should include the constant $`Ng^2L/48`$ diverging in the thermodynamic limit (cf. Sec. 3.2.1). If this is kept, $`\mathrm{e}^{\beta ϵ}0`$ and all thermal factors become independent of $`T`$. As a consequence,

$$\frac{}{T}\underset{N\mathrm{}}{lim}\frac{\overline{q}q}{N}=0,\underset{N\mathrm{}}{lim}\frac{p}{N}=0.$$
(217)
A system of independent quarks cannot be heated up due to confinement. This leaves no room for a chiral phase transition, not even in the limit $`T\mathrm{}`$.

Let us now turn to the finite extension alternative. The ’t Hooft model on a line reduces to a purely fermionic theory with linear Coulomb interaction, see Sec. 3.2. If we simply take this theory and put it on a finite interval with antiperiodic boundary conditions for the quarks, we can repeat the HF calculation of the vacuum (Sec. 3.2.1) numerically, in analogy to what we did analytically in the GN model (Sec. 4.2). This calculation has in fact been done long ago, although for reasons totally unrelated to finite temperature.$`^{\mathrm{?},\mathrm{?}}`$ A second order phase transition to a chirally restored phase was found at a critical interval length of $`L_c19.4/\sqrt{Ng^2}`$. Here, the quark self-energies cannot be responsible since they do not enter the gap equation which is perfectly well behaved in the IR, see Sec. 3.2.1. So something must have gone wrong with covariance — the equivalence of finite extension and finite temperature which we have verified in the GN model seems to have been lost in the gauge theory.

The solution to this puzzle was found only recently.$`^\mathrm{?}`$ It turns out that as soon as one considers finite $`L`$, it is not true that the gluons in QCD2 can be completely eliminated in favor of a static Coulomb potential between quarks. The large $`N`$ limit of the zero mode gluons (discussed for pure Yang-Mills in Sec. 2.5) needs to be refined. Eventually, the finite interval calculation yields the same null-result as the thermal HF approach. In effect, the zero mode gluons make it impossible to compress the quark system, in the same way as the IR-divergent single particle energies make it impossible to heat the system up (to leading order in $`1/N`$). We are familiar with the fact that gauge theories can hide their microscopic degrees of freedom; it is perhaps less well known that they can even hide the structure of space-time, making a cylinder look like an infinite plane (“decompactification” $`^\mathrm{?}`$). Since this is a novel and unexpected phenomenon, it may be worth going briefly through the analysis.

It is easy to include dynamical fermions into the pure Yang-Mills theory of Sec. 2.5. We again work canonically on a spatial circle of length $`L`$ in the gauge $`_1A_1=0`$, $`(A_1)_{ij}=\delta _{ij}\frac{\phi _i}{gL}`$ diagonal in color. The Hamiltonian reads $`^{\mathrm{?},\mathrm{?},\mathrm{?}}`$

$$H=H_\mathrm{g}+H_\mathrm{f}+H_\mathrm{C},$$
(218)
with the gauge field kinetic energy

$$H_\mathrm{g}=\frac{g^2L}{4}\underset{i}{}\frac{^2}{\phi _i^2},$$
(219)
the quark kinetic energy

$`H_\mathrm{f}`$ $`=`$ $`{\displaystyle \underset{n,i}{}}{\displaystyle \frac{2\pi }{L}}\left(n+{\displaystyle \frac{\phi _i}{2\pi }}\right)\left(a_i^{}(n)a_i(n)b_i^{}(n)b_i(n)\right)`$ (220)
$`+m{\displaystyle \underset{n,i}{}}\left(a_i^{}(n)b_i(n)+b_i^{}(n)a_i(n)\right),`$
and the Coulomb interaction

$$H_\mathrm{C}=\frac{g^2L}{4}\underset{n,i,j}{}\frac{j_{ij}(n)j_{ji}(n)}{\left(2\pi n\phi _j+\phi _i\right)^2}.$$
(221)
The notation for the quark operators and currents is the same as in Sec. 3.2. As discussed in Sec. 2.5 due to the curved configuration space of the $`\phi _i`$ and the SU($`N`$) Haar measure originally appearing in $`H_\mathrm{g}`$, this Hamiltonian has to be supplemented by the following boundary condition for the wavefunctionals,

$$\mathrm{\Psi }(\phi _1,\mathrm{},\phi _N;\mathrm{fermions})=0\mathrm{if}\phi _i=\phi _j\mathrm{mod}2\pi .$$
(222)
In Ref. 45 quantization was performed after complete classical gauge fixing. Thereby, the fact that the $`\phi _i`$ are curvilinear coordinates was ignored. This led to the assumption that all the $`\phi _i`$ are frozen at the value $`\pi `$ in the large $`N`$ limit. In the resulting purely fermionic theory, the only remnant of the gluons are antiperiodic boundary conditions for the quarks in the compact space direction. In the meantime, this whole approach has been put on a more rigorous basis by first quantizing in the Weyl gauge and then resolving the Gauss law quantum mechanically.$`^\mathrm{?}`$ This made it clear that the $`\phi _i`$ are parameters on the group manifold with corresponding Jacobian, the SU($`N`$) reduced Haar measure. When solving the theory, it is then possible to restrict oneself to the smallest region in field space bounded by zeros of the Jacobian, see Ref. 27 where the consequences for SU(2), SU(3) have been explored in the case of adjoint fermions. How can this be generalized to the large $`N`$ limit? A definite choice of “fundamental domain” obviously means that the $`\phi _i`$ always remain ordered, say $`0\phi _1\phi _2\mathrm{}\phi _N2\pi `$. If we think of the gluons as particles on a circle, they cannot cross each other and become closely packed in the limit $`N\mathrm{}`$. Their fluctuations are suppressed by $`1/N`$, simply due to lack of space. The only degree of freedom left, the collective rotation of this “pearl necklace”, is a U(1) factor which anyway is not present in the SU($`N`$) theory. This suggests that the correct choice for the gluon background field as seen by the fermions is not $`\phi _i=\pi `$, but rather the continuum limit of the lattice points

$$\phi _i=2\pi \frac{i}{N},(i=1\mathrm{}N)$$
(223)
(see Fig. 12).

Instead of antiperiodic boundary conditions, the fermions then acquire color dependent boundary conditions which interpolate smoothly between the phases 0 and $`2\pi `$,

$$\psi _k(L)=\mathrm{e}^{\mathrm{i2}\pi k/N}\psi _k(0),(k=1\mathrm{}N).$$
(224)
In the thermodynamic limit $`L\mathrm{}`$, both of these choices of the gluon field configuration, $`\phi _i=\pi `$ or Eq. (223), will become indistinguishable and yield the well-known results. We now show that at finite $`L`$ the effects of the gluons on the quarks is radically different in these two cases, and that it is the spread out distribution (223) which is in fact the correct one.

In such a gluonic background field, the fermions can be treated in a relativistic HF approximation along the lines explained in Sec. 3.2.1 except that the Bogoliubov angle acquires a color index. The gap equation then becomes

$$\frac{2\pi }{L}(n+\alpha _i)\mathrm{sin}\theta _i(n)m\mathrm{cos}\theta _i(n)+\frac{g^2L}{16\pi ^2}\underset{n^{},j}{}\frac{\mathrm{sin}\left(\theta _i(n)\theta _j(nn^{})\right)}{(n^{}\alpha _j+\alpha _i)^2}=0$$
(225)
where we have switched to the slightly more convenient variable $`\alpha _i=\frac{\phi _i}{2\pi }[0,1]`$ for the gluons. If $`\alpha _i=1/2`$ as chosen in Ref. 45, $`\theta _i(n)`$ becomes $`i`$-independent and we recover the old gap equation considered in that work. Now, we assume $`\alpha _i=i/N`$ and perform the large $`N`$ limit before solving the gap equation. Since $`\alpha _i`$ becomes a continuous variable, we replace $`\theta _i(n)\theta _\alpha (n)`$ and $`_jN_0^1d\alpha ^{}`$, with the result

$`{\displaystyle \frac{2\pi }{L}}(n+\alpha )\mathrm{sin}\theta _\alpha (n)m\mathrm{cos}\theta _\alpha (n)`$
$`+{\displaystyle \frac{Ng^2L}{16\pi ^2}}{\displaystyle \underset{n^{}}{}}{\displaystyle _0^1}d\alpha ^{}{\displaystyle \frac{\mathrm{sin}\left(\theta _\alpha (n)\theta _\alpha ^{}(nn^{})\right)}{(n^{}\alpha ^{}+\alpha )^2}}=0.`$ (226)
This infinite set of coupled integral equations collapses into a single, one-dimensional integral equation, if we set

$$\theta _\alpha (n)=\theta (n+\alpha ).$$
(227)
Since $`n`$ is integer and $`\alpha [0,1]`$, this step in effect decompactifies the original spatial circle. With this ansatz, the notation $`n+\alpha =\nu `$, $`nn^{}+\alpha ^{}=\nu ^{}`$ (where $`\nu ,\nu ^{}`$ are dimensionless, continuous variables) and the substitution $`_n^{}_0^1d\alpha ^{}_{\mathrm{}}^{\mathrm{}}d\nu ^{}`$, we obtain

$$\frac{2\pi }{L}\nu \mathrm{sin}\theta (\nu )m\mathrm{cos}\theta (\nu )+\frac{Ng^2L}{16\pi ^2}\mathrm{d}\nu ^{}\frac{\mathrm{sin}(\theta (\nu )\theta (\nu ^{}))}{(\nu \nu ^{})^2}=0.$$
(228)
The simple pole in the integral has been regularized by the principal value prescription. After rescaling the variables via

$$\frac{2\pi }{L}\nu :=p,\frac{2\pi }{L}\nu ^{}:=p^{},$$
(229)
where $`p,p^{}`$ have the dimension of momenta, we recover exactly the continuum version of the HF equation, namely

$$p\mathrm{sin}\theta \left(\frac{Lp}{2\pi }\right)m\mathrm{cos}\theta \left(\frac{Lp}{2\pi }\right)+\frac{Ng^2}{4}\frac{\mathrm{d}p^{}}{2\pi }\frac{\mathrm{sin}\left(\theta \left(\frac{Lp}{2\pi }\right)\theta \left(\frac{Lp^{}}{2\pi }\right)\right)}{(pp^{})^2}=0.$$
(230)
Denoting the Bogoliubov angle of the continuum ’t Hooft model by $`\theta _{\mathrm{cont}}(p)`$, we conclude that

$$\theta (\nu )=\theta _{\mathrm{cont}}\left(\frac{2\pi }{L}\nu \right),$$
(231)
or, in terms of the original, color dependent Bogoliubov angle,

$$\theta _i(n)\theta _{\mathrm{cont}}\left(\frac{2\pi }{L}\left(n+\frac{i}{N}\right)\right),(N\mathrm{}).$$
(232)
This last relation becomes exact in the limit $`N\mathrm{}`$ only. In this limit, the color- and $`L`$-dependent Bogoliubov angles for the ’t Hooft model on the circle of length $`L`$ are all given by one universal function, namely the momentum dependent Bogoliubov angle of the ’t Hooft model on the infinite line. We emphasize that this universality only holds for the “pearl necklace” type distribution of gauge variables, Eq. (223). If the $`\phi _i`$ are all set equal to $`\pi `$, there is no analytic way known how to relate $`\theta (n)`$ for different $`L`$ values, but one has to solve the gap equation numerically for each $`L`$.$`^\mathrm{?}`$

The upshot of this simple exercise is the following: In the large $`N`$ limit, the gluon variables influence the fermion boundary conditions in such a way that the circle gets replaced by a line — they decompactify space-time (see Fig. 13; this only works for gauge invariant variables, where one sums over all colors). The length $`L`$ of the spatial circle becomes an irrelevant parameter. To confirm this last point, let us evaluate the condensate as an example of an observable,

$$\overline{q}q=\frac{1}{L}\underset{n,i}{}\mathrm{sin}\theta _i(n)=\frac{N}{L}\underset{n}{}_0^1d\alpha \mathrm{sin}\theta (n+\alpha )=N\frac{\mathrm{d}p}{2\pi }\mathrm{sin}\theta _{\mathrm{cont}}(p).$$
(233)
The sum over the discrete momenta and the color sum in the large $`N`$ limit conspire to produce the continuum result, independently of the starting $`L`$ value. In the alternative thermodynamic interpretation, the condensate becomes independent of temperature to leading order in $`1/N`$, now in perfect agreement with the finite temperature HF result.

Notice that so far, we have discussed the influence of the gauge fields on the quarks which is indeed dramatic. Vice versa, we do not expect the quarks to influence significantly the gluon zero point motion, again due to the effects of the Jacobian. The kinetic energy $`H_\mathrm{g}`$ will give the same result as in pure Yang Mills theory. Since this contribution to the energy density is $`L`$-independent, it again yields zero pressure, reflecting the absence of physical gluonic excitations in 1+1 dimensions.

### 4.4 $`1/N`$ corrections at finite $`T`$, role of pions

In the present section, we wish to discuss the role of the pions at low temperature. Since this has to do with the next order in the $`1/N`$ expansion and cannot be covered systematically here, we will content ourselves with briefly reviewing the state of the art and indicating the issues which have attracted most attention recently.

The importance of pions for low temperature QCD was pointed out by Leutwyler and collaborators $`^{\mathrm{?},\mathrm{?}}`$ some time ago. In the vicinity of the chiral limit, pions are naturally the lightest bosons and as such preferred degrees of freedom at low temperature. Using chiral perturbation theory, interaction effects in the pion gas can be computed reliably to rather high order. An extrapolation of the temperature dependence of the chiral condensate even gives a surprisingly good prediction for the critical temperature of the QCD chiral phase transition. Glueballs (if they exist in nature) are much too heavy to matter at low temperatures, quarks are anyway ruled out due to confinement.

The situation is somewhat different in the case of the NJL model which does not exhibit confinement. However also here thermal excitation of a $`q\overline{q}`$ pair is strongly suppressed as compared to pions at low $`T`$. It was soon recognized that the popular $`1/N`$ expansion is somewhat misleading physically: To leading order, the pions do not contribute at all to thermodynamic observables, they are dominated by the quarks due to their color weight factor (entropy). For any finite $`N`$, one can estimate a temperature below which pions dominate, but for $`N\mathrm{}`$ this temperature moves to $`0`$. Zhuang, Klevansky and Hüfner $`^{\mathrm{?},\mathrm{?}}`$ and others have developed techniques to go to next-to-leading order in the $`1/N`$ expansion, e.g. in the thermodynamic potential. They take into account fluctuations about the mean field which correspond diagrammatically to a ring sum and physically to mesonic contributions. Both $`q\overline{q}`$ bound states and scattering phase shifts enter. Since a meson in the NJL model can decay into quarks, a “mixed phase” arises at intermediate temperatures where mesons as well as quarks play a role. This problem was also investigated by Barducci et al. $`^{\mathrm{?},\mathrm{?}}`$ including the temperature dependence of meson masses.

The same issues are relevant for the lower dimensional models which are the subject of the present work as well. Since these are toy models anyway, one usually takes the limit $`N\mathrm{}`$ seriously (not as an approximation of $`N=3`$), at least this is the point of view which we have adopted here. In this case, one might argue that the pions have no influence on the phase diagram of the chiral GN model. In the strict chiral limit, this is a rather extreme picture: $`N`$ has to be so large that massive fermions win over massless pions in spite of their exponential suppression through the Boltzmann factors, but as a mathematical limit it is perhaps acceptable.

It is instructive to think about the $`1/N`$ corrections (meson effects) in the “rotated” picture of finite extension field theory. There, in next-to-leading order one would have to evaluate the RPA vacuum in a finite interval, in clear correspondence to the bubble sum of Zhuang et al. $`^\mathrm{?}`$ Unlike the thermodynamic calculation, this does not require any new formalism since the RPA technique is well established. Fermionic and mesonic degrees of freedom are then treated consistently and it is also clear that scattering states as well as bound states enter. Such a calculation has not yet been carried out for the GN model. The thermodynamic treatment of the $`1/N`$ corrections has been presented by Barducci et al. $`^\mathrm{?}`$ showing low temperature dominance of pion-like excitations in the massive GN model at finite $`N`$. Introducing a bare quark mass enables the authors to vary the degree of explicit chiral symmetry breaking.

Let us finally come back to the ’t Hooft model. Here, due to confinement, the quarks do not contribute at all to the pressure at order $`N`$ so that the mesonic contributions are in fact the leading ones. Since in the large $`N`$ limit one has non-interacting mesons, one can use the thermodynamics of the ideal Bose gas to analyze this situation. Near the chiral limit, the pion will naturally stick out at low temperatures.

A theoretical question of some interest is the corresponding finite interval RPA calculation in the ’t Hooft model. Here, we do not expect that the mechanism which leads to “decompactification” in leading order remains intact, otherwise one would also not get any meson contribution to the pressure. This must mean that the fluctuations of the zero mode gluons come into the picture. It would be interesting to study this in more detail since it would teach us something about the gluon dynamics at large but finite $`N`$ and perhaps help to settle some of the difficult issues raised by Dhar et al. $`^{\mathrm{?},\mathrm{?},\mathrm{?}}`$ in connection with the string picture of QCD2.

## 5 Finite density

In the present section, we address the problem of finite baryon density or chemical potential in the GN and ’t Hooft models. Since the physics issues are quite different from those at finite temperature (Sec. 4), we restrict ourselves to the $`T=0`$ case here and try to clarify the structure of cold baryonic matter first. Ultimately, we are of course interested in the phase diagram of our soluble models in the whole ($`T,\mu `$)-plane; this will be the subject of Sec. 6. Whereas the GN model at finite density has received a lot of attention in the literature,$`^{\mathrm{?},\mathrm{?},\mathrm{?},\mathrm{?},\mathrm{?}}`$ the corresponding extension of the ’t Hooft model has remained essentially unexplored until very recently.$`^{\mathrm{?},\mathrm{?}}`$ Throughout this section, our main focus will be on the relationship between the structure of the single baryon and that of baryonic matter.

### 5.1 GN model as a Fermi gas

To determine the ground state of the GN model at finite density, we can again use the HF method expected to become exact in the large $`N`$ limit. As we shall see, this does not mean that one should trust the results of a HF calculation blindly: If the ground state breaks some symmetry, the HF equations will develop different solutions and one may easily miss the solution with the lowest energy due to erroneous assumptions. Anyway, to start our investigation let us first assume translational invariance and only offer the system the chance to break chiral symmetry. We thus simply fill a number of positive energy plane wave states on top of the Dirac sea until we reach the desired baryon density. This is the standard approach to the finite density GN model, phrased in many-body language. Although we have in mind primarily the chiral GN model, the corresponding calculation would in fact be identical for the model with discrete symmetry only.

We denote the fermion density per color (or baryon density) by $`\rho _B=p_f/\pi `$ ($`p_f`$: Fermi momentum). At the mean field level, the fermions acquire a physical mass $`m`$ which has to be determined self-consistently. The ground state energy density per color is then given by

$$\frac{}{N}=2_{p_f}^{\mathrm{\Lambda }/2}\frac{\mathrm{d}k}{2\pi }\sqrt{m^2+k^2}+\frac{m^2}{2Ng^2}$$
(234)
where $`\mathrm{\Lambda }`$ is the UV cutoff. The first term sums up the single particle energies for all occupied states (the Dirac sea plus all positive energy states with $`|p|<p_f`$), the second term corrects for double counting of the interaction. We renormalize the theory at $`p_f=0`$, cf. Sec. 3.1.1, and denote the physical fermion mass in the vacuum by $`m_0`$. Using the vacuum gap equation (78) to renormalize the matter ground state energy density, Eq. (234), we find (dropping an irrelevant infinity)

$$\frac{}{N}=\frac{m^2}{4\pi }+\frac{p_f}{2\pi }\sqrt{p_f^2+m^2}+\frac{m^2}{2\pi }\mathrm{ln}\left(\frac{p_f+\sqrt{m^2+p_f^2}}{m_0}\right).$$
(235)
The energy is minimal provided $`m`$ satisfies

$$m\mathrm{ln}\left(\frac{p_f+\sqrt{m^2+p_f^2}}{m_0}\right)=0,$$
(236)
i.e., for

$$m=0\text{or}m=m_0\sqrt{1\frac{2p_f}{m_0}}\left(p_f<\frac{m_0}{2}\right).$$
(237)
The corresponding energy densities are

$`{\displaystyle \frac{}{N}}|_{m=0}`$ $`=`$ $`{\displaystyle \frac{p_f^2}{2\pi }},`$
$`{\displaystyle \frac{}{N}}|_{m0}`$ $`=`$ $`{\displaystyle \frac{m_0^2}{4\pi }}+{\displaystyle \frac{p_fm_0}{\pi }}{\displaystyle \frac{p_f^2}{2\pi }}\left(p_f<{\displaystyle \frac{m_0}{2}}\right).`$ (238)
The physical quark masses (237) and the energy densities (238) are plotted in Figs. 15 and 15.

From these figures one might be tempted to conclude that chiral symmetry is broken at low densities and gets restored in a second order phase transition at $`p_f=m_0/2`$. As is well known, this does not occur, rather there is a first order chiral phase transition at $`p_f=m_0/\sqrt{2}`$. This can easily be inferred by inspection of the thermodynamic potential of the GN model.$`^{\mathrm{?},\mathrm{?}}`$ At zero temperature, it is possible to understand it also more directly in terms of the following instability. Let us compare the energy densities (238) with the energy density for a system of size $`L`$ divided into two homogeneous regions I (size $`\mathrm{}`$) and II (size $`L\mathrm{}`$). In region I chiral symmetry is restored; it contains the extra fermions needed to get the prescribed average density (the “MIT bag” $`^\mathrm{?}`$). Region II consists of the physical vacuum with broken chiral symmetry, void of excess fermions. The mean energy density obtained in this way is

$$\frac{}{N}=\left(\frac{L\mathrm{}}{L}\right)\frac{m_0^2}{4\pi }+\frac{Lp_f^2}{2\pi \mathrm{}}.$$
(239)
Minimization with respect to $`\mathrm{}`$ yields

$$\mathrm{}=\frac{\sqrt{2}p_fL}{m_0}$$
(240)
valid for $`p_f<m_0/\sqrt{2}`$, and hence the optimal energy density

$$\frac{}{N}=\frac{m_0^2}{4\pi }+\frac{p_fm_0}{\sqrt{2}\pi }\left(p_f<\frac{m_0}{\sqrt{2}}\right).$$
(241)
As shown in Fig. 15, this solution is lower in energy than the homogeneous one; moreover, it yields the convex hull of $``$. It ends exactly at the first order phase transition point $`p_f=m_0/\sqrt{2}`$ where all space is filled with one big bag. This should be contrasted to the scenario underlying Fig. 15 where the fermion mass decreases continuously. We thus recover the well known mixed phase interpretation of the GN model at finite density, in the zero temperature limit. Notice also that only the total size of regions I and II matters, not how they are subdivided; there could be baryon “droplets” as well.$`^\mathrm{?}`$

One important point to which we would like to draw the attention of the reader is the behavior of $``$ near $`\rho _B=p_f/\pi =0`$. Since ultimately, at very low density, the fermionic matter problem must reduce to the problem of a single baryon, one would expect

$$\frac{}{\rho _B}|_{\rho _B=0}=M_B$$
(242)
where $`M_B`$ is the baryon mass. In the present calculation, $`M_B`$ is not the physical baryon mass, but the mass of an alleged “delocalized” baryon. This is inherent in the translationally invariant HF approach, i.e., the assumption that the single particle orbitals are momentum eigenstates. Using Eq. (242) we obtain in the homogeneous, single phase calculation, Eq. (238), $`M_B=Nm_0`$, consistent with a short range force and a delocalized baryon. The (physically more viable) mixed phase approach, Eq. (241), predicts a baryon mass lower by a factor of $`1/\sqrt{2}`$ which can be understood in terms of the bag model for the baryon.$`^\mathrm{?}`$

However, the GN model possesses bound baryons with lowest mass $`Nm_0/\pi `$ (kink solution for the model with discrete chiral symmetry, cf. Sec. 3.1.2), or even massless baryons (model with continuous chiral symmetry, cf. Sec. 3.3). These binding effects are not $`1/N`$ suppressed and should be correctly reproduced in a HF approach, in the low density limit. They have obviously been missed here due to our tacit assumption of translational invariance. There is no good reason why such effects should not play a role at higher densities as well. Moreover, differences between the continuous and discrete chirally symmetric GN models based on their different baryon structure and spectra are not at all captured by the Fermi gas approach. In Sec. 5.3 we shall present a cure for this disease. Before that however, let us first repeat the naive calculation for the ’t Hooft model where we can also investigate the role of confinement.

### 5.2 ’t Hooft model as a Fermi gas

We have discussed the vacuum and baryons in the ’t Hooft model in Sec. 3.2. Assuming translational invariance, it is again straightforward to include a finite baryon density into the HF calculation. If $`p_f`$ denotes the Fermi momentum we only need to replace the density matrix (143) by

$`\rho (p)`$ $`=`$ $`\mathrm{\Theta }(p_f|p|)u(p)u^{}(p)+v(p)v^{}(p)`$ (243)
$`=`$ $`\mathrm{\Theta }(p_f|p|)+\mathrm{\Theta }(|p|p_f)v(p)v^{}(p)`$
where we have used the completeness relation for the spinors in the second step. In the expression for the HF ground state energy density at $`m=0`$ (148), according to the second line of Eq. (243), we must exclude the region $`[p_f,p_f]`$ from the momentum integrations and pick up an additional term due to the change in the baryon density $`\mathrm{tr}\rho `$,

$`{\displaystyle \frac{}{N}}`$ $`=`$ $`{\displaystyle \frac{\mathrm{d}p}{2\pi }\mathrm{\Theta }(|p|p_f)p\mathrm{cos}\theta (p)}`$ (244)
$`{\displaystyle \frac{Ng^2}{8}}{\displaystyle \frac{\mathrm{d}p}{2\pi }\frac{\mathrm{d}p^{}}{2\pi }\mathrm{\Theta }(|p|p_f)\mathrm{\Theta }(|p^{}|p_f)\frac{\mathrm{cos}(\theta (p)\theta (p^{}))1}{(pp^{})^2}}`$
$`+{\displaystyle \frac{Ng^2}{4}}{\displaystyle \frac{\mathrm{d}p}{2\pi }\frac{\mathrm{d}p^{}}{2\pi }\mathrm{\Theta }(p_f|p|)\mathrm{\Theta }(|p^{}|p_f)\frac{1}{(pp^{})^2}}.`$
This yields at once the following finite density generalization of the gap equation,

$$p\mathrm{sin}\theta (p)+\frac{Ng^2}{4}\frac{\mathrm{d}p^{}}{2\pi }\mathrm{\Theta }(|p^{}|p_f)\frac{\mathrm{sin}\left(\theta (p)\theta (p^{})\right)}{(pp^{})^2}=0,(|p|>p_f),$$
(245)
whereas the condensate now becomes

$$\overline{q}q=N\frac{\mathrm{d}p}{2\pi }\theta (|p|p_f)\mathrm{sin}\theta (p).$$
(246)
The gap equation (245) can easily be solved numerically for various $`p_f`$.

The resulting condensate is shown in Fig. 16. We find that it decreases monotonically with increasing density, disappearing at a critical Fermi momentum

$$p_f^c0.117\left(\frac{Ng^2}{2\pi }\right)^{1/2}.$$
(247)
This behavior is strikingly similar to the corresponding result for the GN model depicted in Fig. 15, again suggesting some phase transition with restoration of chiral symmetry at high density. Unfortunately, here we are not able to assess whether we are dealing with a first or second order phase transition. The reason lies in the following problem: If we compute the energy density (244) for the ’t Hooft model we discover that subtraction of the value at $`p_f=0`$ is not sufficient to give a finite result. Unlike in the GN model, the difference is still IR divergent. To be able to proceed, we enclose the system in a box of length $`L`$. We then find that the divergence is due to the last term in Eq. (244) (the one which does not involve the Bogoliubov angles) which now contributes the following double sum to the energy per color,

$$\frac{E}{N}|_{\mathrm{div}}=\frac{Ng^2L}{16\pi ^2}\underset{pI}{}\underset{n0,(pn)I}{}\frac{1}{n^2}.$$
(248)
Here antiperiodic boundary conditions for fermions have been employed in the box regularization, and correspondingly the interval $`I`$ is defined in the following way,

$$I=[n_f,n_f]\text{for}B=2n_f+1\text{odd},I=[n_f1,n_f]\text{for}B=2n_f+2\text{even}.$$
(249)
The result (248) is even more alarming than the non-convex behavior of $``$ in the GN model, Fig. 15, due to its $`L`$-dependence. Adding quarks to the vacuum causes the energy to increase by an infinite amount in the limit $`L\mathrm{}`$. Evaluating the double sums in Eq. (248) for low values of $`B`$, we obtain information on the origin of this divergent behavior. For $`B=1`$ ($`I=[0,0]`$) in particular, the calculated baryon mass (to leading order in $`L`$) is

$$M_B=N\left(\frac{Ng^2L}{48}\right).$$
(250)
This is the same relation as $`M_B=Nm_0`$ in the GN model except that the physical fermion mass is replaced by the infinite constant $`Ng^2L/48`$ characteristic of confinement, cf. Eq. (151).

Summarizing, the problems encountered in the GN model with translationally invariant baryonic matter again show up in the ’t Hooft model, although in aggravated form. The physics reason is clear: In the GN model the cost of distributing $`N`$ fermions over the whole space is governed by their physical mass; in the ’t Hooft model, due to confinement of quarks, the corresponding quark effective mass diverges with the volume. On the other hand, it is known that both models do possess massless, delocalized baryons in the chiral limit. Evidently, this has to be accounted for, and we conclude that the naive, translationally invariant HF approximation cannot be trusted.

### 5.3 Breakdown of translational invariance

#### 5.3.1 Strict chiral limit

We now return to the massless baryons introduced in Sec. 3.3 in the chiral limit of two-dimensional field theories. By using literally the same techniques (except for a different value of the baryon number) we can easily find the ground state of the system for any baryon density. As discussed in Sec. 3.3 and illustrated in Fig. 8, the single baryon can be visualized in terms of a “chiral spiral” parametrized by $`\chi (x)`$ with one single turn in the whole volume. A finite density $`\rho _B=B/L=p_f/\pi `$ on the other hand implies that

$$\chi (x)=p_f(xx_0),$$
(251)
i.e., one full rotation over a physical distance which now has a well defined limit for $`L\mathrm{}`$, namely 2/$`\rho _B`$. The baryon density remains constant in space for the reasons discussed in Sec. 3.3, but the condensates are modulated as

$`\overline{q}q`$ $`=`$ $`\overline{q}q_v\mathrm{cos}2p_f(xx_0),`$
$`\overline{q}\mathrm{i}\gamma _5q`$ $`=`$ $`\overline{q}q_v\mathrm{sin}2p_f(xx_0).`$ (252)
They can again be viewed as projections of a chiral spiral of radius $`|\overline{q}q_v|`$ onto two orthogonal planes, see Fig. 17.

This state breaks translational symmetry; it is a crystal. In fact, it may be viewed as the simplest possible realization of the old idea of a Skyrme crystal,$`^\mathrm{?}`$ here in the context of large $`N`$ two-dimensional field theories. It is meaningless to ask where one baryon begins and ends; each full turn of the spiral adds one unit to the baryon number. Only the condensates reveal that translational symmetry has been broken down to a discrete subgroup. The energy density of this unusual kind of “nuclear matter” is simply (after subtracting the vacuum energy density)

$$\frac{}{N}=\frac{p_f^2}{2\pi }.$$
(253)
Surprisingly, this is exactly what one would get for a free Fermi gas of massless quarks although Eq. (253) holds for interacting theories where the vacuum has lower energy due to chiral symmetry breaking. In Fig. 18 we compare the energy density for this state to the ones discussed above for the GN model, where translational symmetry had been assumed. The crystal is always energetically favored, the dependence on $`p_f`$ is now convex, and there is no trace of a phase transition, neither first nor second order, at any density. The horizontal slope at $`p_f=0`$ correctly signals the presence of massless baryons and eliminates the above-mentioned problems with the spurious massive, delocalized baryons. We cannot even draw the corresponding picture for the ’t Hooft model, simply because in this case the quark Fermi gas is infinitely higher in energy than the Skyrme crystal for $`L\mathrm{}`$. Nevertheless, all the results for baryonic matter discussed in this section apply to the ’t Hooft model as well.

In the high density limit the oscillations of the condensates become more and more rapid. If we are interested only in length scales large as compared to $`1/p_f`$, the condensates average to zero. In this sense, one might argue that chiral symmetry gets restored at high density, although not in the naive way suggested by Figs. 15 or 15.

Finally, we remark that the “chiral spiral” ground state for fixed baryon density still preserves one continuous, unbroken symmetry, namely the combination of translation and chiral rotation generated by $`P+p_fQ_5`$ ($`P`$: momentum operator, $`Q_5`$: axial charge). One would therefore expect that RPA excitations on this ground state $`^{\mathrm{?},\mathrm{?}}`$ (or mesons in nuclear matter) will have only one collective, gapless mode, a hybrid of a “phonon” and a “pion”.

#### 5.3.2 Non-vanishing bare quark masses

As discussed in Sec. 3.3 following Ref. 44, a finite bare quark mass can be included in the Skyrme type approach to the baryon; it changes the effective action for the chiral phase $`\chi (x)`$ from that of a free massless field into a sine-Gordon model. The baryon was identified with the kink solution. It is then clear which type of solution is likely to be a good candidate for baryonic matter: The sine-Gordon kink crystal. It goes over into the results of the previous section in the limit $`m0`$ and is an approximate solution to the HF equations (presumably a very accurate one, judging from the results of Ref. 44 for the single baryon). Luckily, the sine-Gordon kink crystal has already been studied thoroughly in the literature, first in solid state physics $`^{\mathrm{?},\mathrm{?}}`$ and more recently as a toy model for the Skyrme crystal,$`^\mathrm{?}`$ in terms of Jacobi elliptic functions and elliptic integrals.$`^\mathrm{?}`$ We take over the results from Ref. 116 which is close in spirit to the present study although the authors apparently did not have in mind two-dimensional large $`N`$ field theories. Adapting the formulae of their work to our notation, the following steps allow us to generalize the Skyrme crystal of the previous section to small, finite bare quark masses: Let $`m_\pi `$ denote the mass of the Goldstone boson, Eq. (185), and $`\overline{\rho }_B=p_f/\pi `$ the average baryon density (this is our definition of $`p_f`$ for the case of broken translational symmetry). We then first have to solve the transcendental equation

$$\frac{\pi m_\pi }{p_f}=2k𝐊(k)$$
(254)
for $`k`$ where $`𝐊(k)`$ is the complete elliptic integral of the first kind. The sine-Gordon kink crystal is then given by the following solution of Eq. (181),

$$\chi (x)=\frac{\pi }{2}+\mathrm{am}(\xi ,k),\xi =\frac{m_\pi }{k}(xx_0),$$
(255)
($`\mathrm{am}(\xi ,k)`$ is the Jacobian elliptic amplitude function). From this, we can express the baryon density and the various condensates in terms of further Jacobian elliptic functions ($`\mathrm{dn},\mathrm{sn},\mathrm{cn}`$) as follows,

$`\rho _B(x)`$ $`=`$ $`{\displaystyle \frac{1}{\pi }}_x\chi (x)={\displaystyle \frac{m_\pi }{\pi k}}\mathrm{dn}(\xi ,k),`$
$`\overline{q}q`$ $`=`$ $`\overline{q}q_v\mathrm{cos}2\chi (x)=\overline{q}q_v\left(\mathrm{cn}^2(\xi ,k)\mathrm{sn}^2(\xi ,k)\right),`$
$`\overline{q}\mathrm{i}\gamma _5q`$ $`=`$ $`\overline{q}q_v\mathrm{sin}2\chi (x)=+\overline{q}q_v2\mathrm{s}\mathrm{n}(\xi ,k)\mathrm{cn}(\xi ,k).`$ (256)
Here $`\xi `$ is as defined in Eq. (255). Finally, the energy per volume of this kind of matter is given by

$$\frac{}{N}=\frac{m_\pi p_f}{4\pi ^2}\left\{\frac{8}{k}𝐄(k)+4k\left(1\frac{1}{k^2}\right)𝐊(k)\right\},$$
(257)
$`𝐄(k)`$ denoting the complete elliptic integral of the second kind.

Let us now illustrate these results in two regimes of interest, namely at low and high density. At low density ($`p_fm_\pi `$), $`k`$ in Eq. (254) approaches 1 exponentially, and the baryon density features a chain of well resolved lumps whose shape is determined by the single kink solution (Fig. 19). Likewise, the condensates behave like those of a single baryon: $`\overline{q}q`$ changes from the vacuum value outside the baryons to its negative in their center whereas $`\overline{q}\mathrm{i}\gamma _5q`$ is peaked in the surface region of each baryon (Figs. 21, 21). These condensates are projections of the distorted “chiral spiral” shown in Fig. 22. The energy (257) for low densities behaves as

$$N\frac{2m_\pi p_f}{\pi ^2}=M_B\rho _B,$$
(258)
showing the expected connection to the baryon mass. At high densities ($`p_fm_\pi `$), $`k`$ approaches 0 like

$$k\frac{m_\pi }{p_f}.$$
(259)
Thus $`\xi `$ in Eq. (255) becomes $`p_f(xx_0)`$. Moreover, for $`k0`$, the Jacobian elliptic functions $`\mathrm{am}(\xi ,k),\mathrm{sn}(\xi ,k),\mathrm{cn}(\xi ,k)`$ are known to reduce to the argument $`\xi `$ and the ordinary trigonometric functions $`\mathrm{sin}\xi `$ and $`\mathrm{cos}\xi `$, respectively. We thus recover the results for the simple chiral spiral in Sec. 3.3 (the parameter $`x_0`$ has to be readjusted to take care of the shift by $`\pi /2`$ in Eq. (255)). The energy in this case is approximately

$$\frac{}{N}\frac{p_f^2}{2\pi }+\frac{m_\pi ^2}{8\pi }.$$
(260)
The condensates look very much like the sin- and cos-functions of the massless case and need not be plotted. The baryon density oscillates around a constant value, reflecting the strong overlap of the baryons, and are well approximated at high density by

$$\rho _B(x)\frac{p_f}{\pi }\left(1\frac{1}{2}\left(\frac{m_\pi }{p_f}\right)^2\mathrm{sin}^2p_f(xx_0)\right).$$
(261)
The behavior of the baryon density $`\rho _B(x)`$ as one increases $`p_f`$ (i.e., the mean density) is illustrated in Fig. 23. In the chiral or high-density limit ($`m_\pi /p_f0`$) $`\rho _B(x)`$ eventually becomes $`x`$-independent. This provides us with another way of understanding the structure of matter described in the previous section, namely as arising from a chain of very extended, strongly overlapping lumps.

#### 5.3.3 Comparison to other works

It is noteworthy that a similar chiral structure of fermionic matter has been reported previously in a variety of models different from the ones discussed here. This indicates that the basic results are more generally valid than our derivation might suggest. Let us briefly go through these works to see what is common and what is different from our case.

The first mention of nonuniform structures we are aware of is Ref. 118 on the massive Schwinger model. Since in the standard Schwinger model baryon charge is confined, one either has to add an inert, uniform background charge to neutralize the system or introduce flavor. In the first case, Fischler et al. find that a spatially varying, periodic charge density is induced which they interpret as 1+1 dimensional Wigner crystal.$`^{\mathrm{?},\mathrm{?}}`$ More recently, the massless Schwinger model at finite chemical potential has also been examined and oscillating chiral condensates with wavenumber $`k=2\mu `$ were found, by bosonization $`^\mathrm{?}`$ as well as by functional methods.$`^\mathrm{?}`$ In the massless case, the charge density becomes uniform. This phenomenon has obviously a lot in common with our findings, although it is different in detail. In particular, the way to avoid conflict with the no-go theorem for spontaneous symmetry breaking in two dimensions must be different since one cannot invoke large $`N`$ arguments; presumably, in the U(1) case, the anomaly and the long range Coulomb interaction play a crucial role. In the generalized Schwinger model with two flavors, one can have baryonic charge for neutral systems. Here, Fischler et al. $`^\mathrm{?}`$ also arrive at a sine-Gordon equation for the light meson but point out that a possible crystal solution would be destroyed immediately by quantum fluctuations. This is of course avoided in our case by the large $`N`$ limit.

QCD2 with flavor at finite chemical potential has been discussed by Christiansen et al. $`^\mathrm{?}`$ in the chiral limit. Oscillatory condensates where found within the path-integral approach. Since these authors work at finite $`N`$, it is not clear how the no-go theorem is avoided. On the other hand, since we know that QCD2 does have massless mesons and baryons at finite $`N`$, one cannot rule out that there is yet another way of bypassing the Coleman theorem.

Even more surprising are perhaps quite a number of reports about spatially inhomogeneous chiral condensates with exactly the same wave number as in our case, but in 3+1 dimensions. Kutschera et al. $`^\mathrm{?}`$ have studied quark matter with pion condensates within the $`\sigma `$-model, using a mean-field approximation. Whereas their quark matter is closely analogous to our Fermi gas, the pion condensed state resembles the chiral spiral, especially as far as the condensates are concerned. Another important work is the one by Deryagin et al. $`^\mathrm{?}`$ on large $`N`$ QCD in 3+1 dimensions. With the help of a variational calculation, these authors find an instability of cold and dense matter with respect to formation of an inhomogenous, anisotropic condensate, a standing wave with wave number $`2\mu `$. They point out that unlike color superconductivity, this condensation survives in the $`N\mathrm{}`$ limit. Physically, they interpret their results as condensation of pairs of particles and holes on opposite points of the Fermi surface $`|\stackrel{}{p}|=\mu `$ (“Overhauser effect” $`^\mathrm{?}`$). This scenario was subsequently taken up by Shuster et al. $`^\mathrm{?}`$ and Park et al. $`^\mathrm{?}`$ who elaborated on the competition between BCS (pp or hh) pairing and Overhauser (ph) pairing. They confirm that the large $`N`$ limit favors the chiral wave, although it seems very unlikely to be relevant for $`N`$ values as low as 3.

Let us see how this pairing picture fits into the 1+1 dimensional models considered here. Since relativistic HF is very similar to BCS theory, we can interpret the vacuum result (as encoded in the Bogoliubov angles, cf. Fig. 6, Sec. 3.2.1) as a ph-pairing in a momentum region with a width of order $`m`$ or $`\sqrt{Ng^2}`$ (for the GN and ’t Hooft models, respectively) and centered around $`p=0`$. The crystal solution is obtained by substituting

$$q(x)\mathrm{e}^{\mathrm{i}p_fx\gamma ^5}q(x),$$
(262)
thereby splitting momenta between right-handed and left-handed fermion components by $`\pm p_f`$. Then, provided that $`p_f(m,\sqrt{Ng^2})`$, we can indeed describe the result as ph pairing on opposite sides of the Fermi “sphere” (here an interval). However our results also hold for small $`p_f`$ where this picture is not really applicable.

## 6 Phase diagram in the ($`T,\mu `$)-plane

Most of the work done so far on the phase diagram for the two dimensional large $`N`$ models has been devoted to the original, non-chiral GN model. Its amazingly rich phase diagram has been discussed comprehensively by Wolff $`^\mathrm{?}`$ (for $`m_b=0`$) and by Barducci et al. $`^\mathrm{?}`$ (for finite bare quark masses). It has also been invoked for understanding some questions about the statistical physics of polyacetylene,$`^\mathrm{?}`$ a linear polymer. We are primarily interested in the role played by chiral symmetry and will therefore concentrate on the chiral GN model. Actually, as far as the conventional Fermi gas approach is concerned, both variants of the GN model yield identical results to leading order in $`1/N`$.$`^\mathrm{?}`$ Apart from briefly reviewing the common lore about the GN model, we will also investigate the novel “crystalline” phase (chiral spiral, cf. Sec. 5.3) and determine which phase is thermodynamically stable. In order to keep things technically as simple as possible, we consider only the $`m_b=0`$ case and work to leading order in the $`1/N`$ expansion throughout this section. Incidentally, the ’t Hooft model will be no issue here: Since in the limit $`N\mathrm{}`$ confinement suppresses any $`T`$-dependence (cf. Sec. 4.3), there is nothing we could add to what has already been said about $`T=0`$ and finite density in Sec. 5.

### 6.1 Chiral GN model with translational invariance

Let us assume unbroken translational invariance and, without loss of generality, $`\overline{q}\mathrm{i}\gamma ^5q=0`$. If we are interested in equilibrium thermodynamics at finite temperature and density, it is convenient to work with the grand canonical ensemble where the baryon density is adjusted via a Lagrange multiplier, the chemical potential $`\mu `$. The HF method can be generalized without any difficulty to finite $`\mu `$.$`^\mathrm{?}`$ Thereby, the HF equation remains unchanged and equal to the free Dirac equation,

$$\left(\gamma ^5\frac{1}{\mathrm{i}}\frac{}{x}+m\gamma ^0\right)\phi _n(x)=ϵ_n\phi _n(x),$$
(263)
whereas the chemical potential only enters the gap equation through the standard Fermion occupation numbers

$$\frac{m}{Ng^2}=\frac{\overline{q}q}{N}=\underset{n}{}\overline{\phi }_n(x)\phi _n(x)\frac{1}{\mathrm{e}^{\beta (ϵ_n\mu )}+1}.$$
(264)
The solution of the HF equation is once more given by the free, massive Dirac theory with $`ϵ_n=\pm ϵ(k)=\sqrt{k^2+m^2}`$. The relevant thermodynamic function here is the grand canonical potential often referred to as “effective potential” in this context,

$$V_{\mathrm{eff}}=\frac{1}{\beta }\mathrm{ln}Z=\frac{1}{\beta }\mathrm{ln}\mathrm{tr}\mathrm{e}^{\beta (H\mu N)}.$$
(265)
Recalling that the Hamiltonian of the GN model differs from that of the free massive Dirac theory only by a $`c`$-number (in leading order in $`1/N`$) and using the partition function for a free Fermi gas as follows,

$`\mathrm{ln}Z`$ $`=`$ $`\mathrm{ln}{\displaystyle \underset{n}{}}\left(1+\mathrm{e}^{\beta (ϵ_n\mu )}\right)`$
$`=`$ $`L{\displaystyle \frac{\mathrm{d}k}{2\pi }\mathrm{ln}\left[\mathrm{e}^{\beta (ϵ(k)+\mu )}\left(1+\mathrm{e}^{\beta (ϵ(k)\mu )}\right)\left(1+\mathrm{e}^{\beta (ϵ(k)+\mu )}\right)\right]},`$
the effective potential density ($`𝒱_{\mathrm{eff}}=V_{\mathrm{eff}}/L`$) per color becomes

$$\frac{𝒱_{\mathrm{eff}}}{N}=\frac{m^2}{2Ng^2}\frac{\mathrm{d}k}{2\pi }\left\{ϵ(k)+\mu +\frac{1}{\beta }\mathrm{ln}\left[(1+\mathrm{e}^{\beta (ϵ(k)\mu )})\left(1+\mathrm{e}^{\beta (ϵ(k)+\mu )}\right)\right]\right\}.$$
(267)
In this way, the “vacuum part” ($`T=\mu =0`$) can easily be identified and separated from the “matter part” which in turn is made up of particle and anti-particle contributions. Incidentally, the term

$$\frac{\mathrm{d}k}{2\pi }\mu $$
(268)
in eq. (267) is a pure vacuum term which is usually omitted in the literature; it arises from the infinite fermion density of the Dirac sea. The gap equation (264) can be further worked out to yield

$$m=\frac{Ng^2}{2\pi }dk\frac{m}{ϵ(k)}\left(1\frac{1}{\mathrm{e}^{\beta (ϵ(k)\mu )}+1}\frac{1}{\mathrm{e}^{\beta (ϵ(k)+\mu )}+1}\right),$$
(269)
where we have used

$$\overline{u}(k)u(k)=\overline{v}(k)v(k)=\frac{m}{ϵ(k)}.$$
(270)
Eq. (269) is equivalent to the condition that the effective potential (267) is stationary with respect to $`m`$,

$$\frac{𝒱_{\mathrm{eff}}}{m}=0.$$
(271)
Renormalization can be performed at $`T=0,\mu =0`$ (physical fermion mass $`m_0`$) resulting in the renormalized potential

$`{\displaystyle \frac{𝒱_{\mathrm{eff}}}{N}}|_{\mathrm{ren}}`$ $`=`$ $`{\displaystyle \frac{m^2}{4\pi }}\mathrm{ln}\left({\displaystyle \frac{m^2}{m_0^2}}\right)+{\displaystyle \frac{m_0^2m^2}{4\pi }}`$
$`{\displaystyle \frac{1}{\beta }}{\displaystyle \frac{\mathrm{d}k}{2\pi }\mathrm{ln}\left[\left(1+\mathrm{e}^{\beta (ϵ(k)\mu )}\right)\left(1+\mathrm{e}^{\beta (ϵ(k)+\mu )}\right)\right]}.`$
By adding the “bag pressure” $`m_0^2/(4\pi )`$ and dropping the term in Eq. (268), it has been adjusted to vanish in the vacuum. The effective potential (6.1) contains all the information about the conventional phase diagram of the GN model. (For the path integral derivation and the representation of $`𝒱_{\mathrm{eff}}`$ in terms of sums over Matsubara frequencies instead of integrals over thermal occupation numbers, we refer to Refs. 95, 110.) For the sake of completeness, we record the expressions for standard thermodynamical bulk observables: The free energy density $``$ (or, equivalently, minus the pressure $`p`$) is equal to $`𝒱_{\mathrm{eff}}`$, Eq. (6.1). The quark density $`\rho `$ and the entropy density $`s`$ are obtained by differentiating $`p`$ with respect to $`\mu `$ and $`T`$,

$`{\displaystyle \frac{\rho }{N}}`$ $`=`$ $`{\displaystyle \frac{1}{\pi }}{\displaystyle _0^{\mathrm{}}}dk\left({\displaystyle \frac{1}{\mathrm{e}^{\beta (ϵ(k)\mu )}+1}}{\displaystyle \frac{1}{\mathrm{e}^{\beta (ϵ(k)+\mu )}+1}}\right),`$ (273)
$`{\displaystyle \frac{s}{N}}`$ $`=`$ $`{\displaystyle \frac{1}{\pi }}{\displaystyle _0^{\mathrm{}}}dk\mathrm{ln}\left[\left(1+\mathrm{e}^{\beta (ϵ(k)\mu )}\right)\left(1+\mathrm{e}^{\beta (ϵ(k)+\mu )}\right)\right]`$ (274)
$`+{\displaystyle \frac{\beta }{\pi }}{\displaystyle _0^{\mathrm{}}}dk\left({\displaystyle \frac{ϵ(k)\mu }{\mathrm{e}^{\beta (ϵ(k)\mu )}+1}}+{\displaystyle \frac{ϵ(k)+\mu }{\mathrm{e}^{\beta (ϵ(k)+\mu )}+1}}\right),`$
whereas the internal energy density can now be deduced via

$$\epsilon =Tsp+\mu \rho .$$
(275)
Needless to say, all of these observables have to be evaluated in the absolute minimum of the effective potential with respect to $`m`$, i.e., by solving the gap equation.

The phase diagram which one obtains in this way is depicted in Fig. 24, adapted from Ref. 95. The line AB is a critical line of second order transitions ($`𝒱_{\mathrm{eff}}`$ changes from one minimum to a maximum and a minimum). The point B is a tricritical point located at $`1/\beta =0.3183,\mu =0.6082`$ (all numbers here are in units of $`m_0`$); it separates the second order line from a first order line BD along which $`𝒱_{\mathrm{eff}}`$ has two degenerate minima and a maximum. The endpoint D lies at $`\mu =1/\sqrt{2}`$ where the $`T=0`$ phase transition occurs, cf. Sec. 5.1. Lines BC and BE are boundaries of metastability; when crossed, the potential acquires or loses a second minimum; point C is at $`\mu =1/2`$. In region OABD, chiral symmetry is broken and the quarks are massive; the outside region has unbroken chiral symmetry and massless fermions.

In order to highlight the different behavior of the system when going through a first and second order phase transition, respectively, we have illustrated in Figs. 26 and 26 the evolution of $`𝒱_{\mathrm{eff}}`$ along the $`\mu =0`$ and $`T=0`$ axes (this corresponds to the special cases discussed in Secs. 4 and 5, respectively). Also shown is the behavior of the physical fermion mass. In this way, the unstable HF solution which we have found at finite density (with continuously decreasing mass, see Sec. 5.1) can be identified as corrsponding to a local maximum of $`𝒱_{\mathrm{eff}}`$. In preparing Fig. 26, we have taken advantage of the fact that the integration in Eq. (6.1) can be carried out at $`T=0`$ with the result

$`{\displaystyle \frac{𝒱_{\mathrm{eff}}}{N}}`$ $`=`$ $`{\displaystyle \frac{m^2}{4\pi }}\mathrm{ln}\left({\displaystyle \frac{m^2}{m_0^2}}\right)+{\displaystyle \frac{m_0^2m^2}{4\pi }}`$
$`+{\displaystyle \frac{1}{2\pi }}\theta (\mu m)\left[m^2\mathrm{ln}\left({\displaystyle \frac{\mu +\sqrt{\mu ^2m^2}}{m}}\right)\mu \sqrt{\mu ^2m^2}\right].`$
### 6.2 Chiral spiral at finite temperature

We have seen before that the translationally non-invariant finite density solution (constant baryon density, but oscillating condensates in the chiral limit) of the HF equation at $`T=0`$ can be obtained from the vacuum solution by a space dependent chiral rotation of the quarks, cf. Sec. 3.3. As already noted by Kao et al. $`^\mathrm{?}`$ for the massless Schwinger model and by Sadzikowski et al. $`^\mathrm{?}`$ for the $`\sigma `$\- and NJL models, the same procedure works at finite temperature as well. We show that the Fermi gas solution at zero chemical potential and finite temperature can be mapped onto the chiral spiral type solution at finite $`\mu `$ and the same $`T`$. This will allow us to follow the fate of the crystalline solution as we heat up baryonic matter.

We start from the HF equation at finite $`T,\mu `$, but neither invoking translational invariance nor assuming that $`\overline{q}\mathrm{i}\gamma ^5q`$ vanishes,

$$\left\{\gamma ^5\frac{1}{\mathrm{i}}\frac{}{x}g^2\left(\overline{q}(x)q(x)\gamma ^0+\overline{q}(x)\mathrm{i}\gamma ^5q(x)\mathrm{i}\gamma ^1\right)\right\}\phi _n(x)=\omega _n\phi _n(x).$$
(277)
As before, temperature and chemical potential only show up in the self-consistency relations which now read

$`\overline{q}(x)q(x)`$ $`=`$ $`N{\displaystyle \underset{n}{}}\overline{\phi }_n(x)\phi _n(x){\displaystyle \frac{1}{\mathrm{e}^{\beta (\omega _n\mu )}+1}},`$
$`\overline{q}(x)\mathrm{i}\gamma ^5q(x)`$ $`=`$ $`N{\displaystyle \underset{n}{}}\overline{\phi }_n(x)\mathrm{i}\gamma ^5\phi _n(x){\displaystyle \frac{1}{\mathrm{e}^{\beta (\omega _n\mu )}+1}}.`$ (278)
Guided by the $`T=0`$ results (cf. Sec. 5.3.1), we make the following ansatz for the $`x`$-dependence of the scalar and pseudo-scalar condensates,

$$\overline{q}(x)q(x)=\frac{m}{g^2}\mathrm{cos}(2\mu x),\overline{q}(x)\mathrm{i}\gamma ^5q(x)=\frac{m}{g^2}\mathrm{sin}(2\mu x).$$
(279)
The HF equation for this particular potential can be turned into

$$\left\{\gamma ^5\frac{1}{\mathrm{i}}\frac{}{x}+m\mathrm{e}^{\mathrm{i}\mu x\gamma ^5}\gamma ^0\mathrm{e}^{\mathrm{i}\mu x\gamma ^5}\right\}\phi _n(x)=\omega _n\phi _n(x)$$
(280)
where we have used the fact that

$$\mathrm{e}^{\mathrm{i}\mu x\gamma ^5}\left(\begin{array}{c}\gamma ^0\\ \mathrm{i}\gamma ^1\end{array}\right)\mathrm{e}^{\mathrm{i}\mu x\gamma ^5}=\left(\begin{array}{cc}\hfill \mathrm{cos}2\mu x& \hfill \mathrm{sin}2\mu x\\ \hfill \mathrm{sin}2\mu x& \hfill \mathrm{cos}2\mu x\end{array}\right)\left(\begin{array}{c}\gamma ^0\\ \mathrm{i}\gamma ^1\end{array}\right).$$
(281)
A local chiral rotation of the form

$$\phi _n(x)=\mathrm{e}^{\mathrm{i}\mu x\gamma ^5}\chi _n(x)$$
(282)
then eliminates the $`x`$-dependence from the HF potential,

$$\{\gamma ^5\frac{1}{\mathrm{i}}\frac{}{x}+\mu +m\gamma ^0)\chi _n(x)=\omega _n\chi _n(x),$$
(283)
and enables us to solve the HF equation trivially,

$`\chi _n(x)`$ $`=`$ $`\mathrm{e}^{\mathrm{i}kx}u(k)(\omega _n=ϵ(k)+\mu ),`$
$`\chi _n(x)`$ $`=`$ $`\mathrm{e}^{\mathrm{i}kx}v(k)(\omega _n=ϵ(k)+\mu ),`$ (284)
with $`ϵ(k)=\sqrt{k^2+m^2}`$ and $`u(k),v(k)`$ free, massive spinors. Since

$$\chi _n^{}(x)\mathrm{i}\gamma ^1\chi _n(x)=0,$$
(285)
the two self-consistency conditions (278) collapse into the single condition

$$\frac{m}{Ng^2}=\underset{n}{}\chi _n^{}(x)\chi _n(x)\frac{1}{\mathrm{e}^{\beta (\omega _n\mu )}+1}.$$
(286)
The additive term $`\mu `$ in $`\omega _n`$, Eq. (284), just cancels the chemical potential and we have mapped our problem exactly onto the $`\mu =0`$, translational invariant HF problem at finite $`T`$,

$$\frac{m}{Ng^2}=\underset{n}{}\chi _n^{}(x)\chi _n(x)\frac{1}{\mathrm{e}^{\beta ϵ_n}+1},ϵ_n=\omega _n\mu =\pm \sqrt{k^2+m^2}.$$
(287)
The parameter $`m`$ therefore has to be identified with the physical fermion mass of the translationally invariant vacuum at $`\mu =0`$ and the same finite $`T`$. We can thereby construct a translationally broken solution of the HF problem for any $`(\beta ,\mu )`$; it can be pictured as a chiral spiral whose radius shrinks with increasing temperature until it disappears at $`T_c=m_0\mathrm{e}^\mathrm{C}/\pi `$, cf. Sec. 4.2.

Let us determine the free energy of this new phase. The density of the effective potential per color is still given by eq. (267), or equivalently by

$$\frac{𝒱_{\mathrm{eff}}}{N}=\frac{m^2}{2Ng^2}\frac{1}{L\beta }\underset{n}{}\mathrm{ln}\left(1+\mathrm{e}^{\beta (\omega _n\mu )}\right).$$
(288)
The double counting term depends only on the radius of the chiral circle and is therefore unchanged (cf. Eq. (279)),

$$\overline{q}q^2+\overline{q}\mathrm{i}\gamma ^5q^2=\frac{m^2}{g^4}.$$
(289)
In view of the relation $`\omega _n\mu =ϵ_n`$, it appears superficially as if nothing depended on $`\mu `$ and we got the same result as for the translationally invariant system at zero chemical potential. As we know from the $`T=0`$ case, this argument is too rough — the $`\mu `$-dependence is intimately related to the UV regularization. Let us impose an UV-cutoff on the single particle energies at the bottom of the Dirac sea,

$$\omega _n>\mathrm{\Lambda }/2|k|<\mathrm{\Lambda }/2+\mu .$$
(290)
The entire $`\mu `$-dependence of $`𝒱_{\mathrm{eff}}`$ arises from this shift in the cutoff. We can then easily evaluate the difference between the free energy density of the “chiral spiral” at $`(\beta ,\mu )`$ and the standard solution at $`(\beta ,0)`$,

$`{\displaystyle \frac{𝒱_{\mathrm{eff}}(\beta ,\mu )}{N}}|_{\mathrm{spir}}{\displaystyle \frac{𝒱_{\mathrm{eff}}(\beta ,0)}{N}}|_{\mathrm{stand}}`$ $`=`$ $`{\displaystyle _{(\mathrm{\Lambda }/2+\mu )}^{\mathrm{\Lambda }/2+\mu }}{\displaystyle \frac{\mathrm{d}k}{2\pi }}\sqrt{k^2+m^2}`$ (291)
$`+{\displaystyle _{\mathrm{\Lambda }/2}^{\mathrm{\Lambda }/2}}{\displaystyle \frac{\mathrm{d}k}{2\pi }}\sqrt{k^2+m^2}`$
$``$ $`{\displaystyle \frac{(\mathrm{\Lambda }+\mu )\mu }{2\pi }}.`$
This result does not depend on temperature, a feature characteristic for an UV effect. The fermion density per color can be obtained by differentiation with respect to $`\mu `$; we find the same result as at zero temperature,

$$\frac{\rho }{N}=\frac{}{\mu }\frac{𝒱_{\mathrm{eff}}}{N}=\frac{\mathrm{\Lambda }}{2\pi }+\frac{\mu }{\pi },$$
(292)
hence, with the definition of the Fermi momentum adopted in Sec. 5.3.1,

$$\mu =p_f$$
(293)
(the term $`\mathrm{\Lambda }`$ in Eq. (292) is the fermion density of the Dirac sea). Summarizing, we can now compare the effective potentials for two solutions of the HF equation at finite $`(\beta ,\mu )`$. The Fermi gas solution with unbroken translational invariance has already been given in Eq. (6.1). The crystalline phase on the other hand yields

$$\frac{𝒱_{\mathrm{eff}}}{N}|_{\mathrm{spir}}=\frac{m^2}{4\pi }\mathrm{ln}\left(\frac{m^2}{m_0^2}\right)+\frac{m_0^2m^2}{4\pi }\frac{2}{\beta }\frac{\mathrm{d}k}{2\pi }\mathrm{ln}\left(1+\mathrm{e}^{\beta ϵ(k)}\right)\frac{\mu ^2}{2\pi },$$
(294)
where we have again dropped the term $`\mathrm{\Lambda }\mu /(2\pi )`$. In both cases, one still has to minimize the effective potential with respect to $`m`$. The depth of the respective minima (or physically minus the pressure) allows us to decide which phase is the stable one. The chiral spiral is favored everywhere over the homogeneous phase, below the critical temperature. This can be shown by a straightforward numerical evaluation, but also deduced analytically from the simple relation

$$𝒱_{\mathrm{eff}}(\beta ,\mu )|_{\mathrm{spir}}=𝒱_{\mathrm{eff}}(\beta ,0)N\frac{\mu ^2}{2\pi }.$$
(295)
As the additional term in Eq. (295) does not depend on $`m`$, both $`𝒱_{\mathrm{eff}}(\beta ,\mu )|_{\mathrm{spir}}`$ and $`𝒱_{\mathrm{eff}}(\beta ,0)`$ are minimized by the same $`m`$

$$m_{\mathrm{spir}}(\beta ,\mu )=m(\beta ,0).$$
(296)
The free energy density of the spiral at $`\beta `$ and $`\mu `$ evaluated at its minimum $`m_{\mathrm{spir}}(\beta ,\mu )`$ is thus equal to

$$𝒱_{\mathrm{eff}}(\beta ,\mu )|_{\mathrm{spir},m=m_{\mathrm{spir}}}=𝒱_{\mathrm{eff}}(\beta ,0)|_{m=m(\beta ,0)}N\frac{\mu ^2}{2\pi }.$$
(297)
$`𝒱_{\mathrm{eff}}(\beta ,\mu )|_{m=m(\beta ,\mu )}`$ is a monotonously decreasing function of $`\mu `$ at fixed $`\beta `$ (its derivative with respect to $`\mu `$ at fixed $`\beta `$ is just $`\rho `$).

Furthermore from the phase diagram (Fig. 24) it is clear that for any fixed $`\beta `$ there is a critical chemical potential $`\mu _c`$, where chiral symmetry gets restored. Above this critical chemical potential, $`𝒱_{\mathrm{eff}}(\beta ,\mu )`$ evaluated at its minimum just behaves like the energy density of free fermions in 1+1 dimensions

$$𝒱_{\mathrm{eff}}(\beta ,\mu )|_{m=m(\beta ,\mu )}=N\frac{\mu ^2}{2\pi }\mathrm{for}\mu \mu _c.$$
(298)
Comparing Eq. (297) to $`𝒱_{\mathrm{eff}}(\beta ,\mu )|_{m=m(\beta ,\mu )}`$ it is clear that the chiral spiral is always lower in free energy than the homogeneous phase.

Eq. (295) is also useful if one wants to evaluate the thermodynamic observables for the crystal phase

$`\rho (\beta ,\mu )|_{\mathrm{spir}}`$ $`=`$ $`N{\displaystyle \frac{\mu }{\pi }},`$
$`p(\beta ,\mu )|_{\mathrm{spir}}`$ $`=`$ $`p(\beta ,0)+N{\displaystyle \frac{\mu ^2}{2\pi }},`$
$`s(\beta ,\mu )|_{\mathrm{spir}}`$ $`=`$ $`s(\beta ,0),`$
$`\epsilon (\beta ,\mu )|_{\mathrm{spir}}`$ $`=`$ $`\epsilon (\beta ,0)+N{\displaystyle \frac{\mu ^2}{2\pi }}.`$ (299)
(At $`\mu =0`$, there is no distinction between the Fermi gas and the other phase.) Since the translationally broken phase is the thermodynamically stable phase wherever it exists, the phase diagram of the chiral GN model has only one horizontal second order line at $`T=T_c`$ where chiral symmetry is restored and the radius of the chiral spiral shrinks to zero (Fig. 27); no phase transition occurs if we increase the chemical potential. Note that our results are rather different from a recent study of a non-uniform chiral phase in 3+1 dimensional effective field theories ($`\sigma `$-model, NJL model $`^\mathrm{?}`$). In our language, these authors use the radius and the helix angle of the chiral spiral as independent variational parameters. Since in 1+1 dimensions baryon number is topologically connected to winding number, we do not seem to have this option but can vary only with respect to the radius (or $`m`$).

As a last remark, we wish to comment on the original GN model with only discrete chiral symmetry. The chiral spiral does not exist here and one might think that our analysis does not have anything to say about it. However, the criticism of Sec. 5.1 also applies here. As the non-chiral GN model has (massive) bound baryons, the low density behavior of the energy obtained in standard HF approximation cannot be correct and the phase diagram has to be reconsidered. This is a harder task than for the chiral model since here dynamics is not completely dominated by chiral symmetry. It is quite possible that the first order phase transition survives, but this still has to be worked out in detail.

## 7 Closing remarks

In this paper, we have illustrated the kind of insights one gets from studies of exactly soluble field theories by means of two classic examples, the GN and the ’t Hooft model. The undeniable success of these models stems from a fortunate combination of two very bold substitutions: Replacing the number of spatial dimensions by 1 and the number of colors by $`\mathrm{}`$ (ironically, both are 3 in real QCD). Without the first choice, there is no solvability; without the 2nd, there are hardly enough interesting phenomena due to rigorous no-go theorems. If one is willing to swallow both of these bitter pills, one is awarded with a striking number of analogies between these toy models and the real world which have led to a long list of QCD-motivated applications. These analogies include such central features of QCD as asymptotic freedom, dimensional transmutation, chiral symmetry breaking, and confinement. In view of the fact that these models have the simplest Lagrangians of certain types which one can write down and are not contrived in any sense, this is quite amazing and underlines their pedagogical value for QCD related questions. Nowadays, it may well be that the study of lower dimensional field theories is a better training ground for students interested in QCD than a standard field theory course based on Feynman diagrams and perturbative QED.

Throughout this work we have tried to emphasize the fact that the solution of these models can be found by rather elementary means — we have deliberately chosen the relativistic HF approach as our main tool. This does not mean though that there are no subtleties: We only mention the emergence of baryon number from the UV regularization, or how zero mode gluons affect the thermodynamics of quarks. Such features make the use of analytical methods particularly important and set the scale for the even more subtle properties to be expected from 4-dimensional QCD.

As far as our more specific subject — finite temperature and density behavior — is concerned, we have found that even after more than 25 years these models are still good for some surprises. The strongest impression for us is the way in which chiral symmetry governs not only the hadron spectrum, but also the phase structure at finite density and temperature. Confinement is most clearly seen in the temperature dependence of observables, whereas the dependence on the chemical potential is dominated by chiral symmetry. Other non-trivial insights are related to the realization of the Skyrme picture in a precise fashion, a triple point in the phase diagram of the non-chiral GN model and inhomogeneous chiral condensates. As a by-product, we have resolved some puzzles in the literature and identified issues where further work is needed, for instance the phase diagram of the non-chiral GN model.

In conclusion, we see no indication that models such as the GN or ’t Hooft model are anywhere close to having exhausted their potential as a theoretical laboratory for QCD.

## Acknowledgements

We are grateful to Misha Shifman for the invitation to contribute to this Festschrift in honor of Boris Ioffe, thereby forcing us to clarify some issues before a certain deadline. M.T. thanks his collaborators, colleagues and students for most enjoyable and useful interactions on 2d field theory, in particular M. Burkardt, V. Dolman, M. Engelhardt, O. Jahn, M. Kurz, F. Lenz, S. Levit, R. Oderkerk, K. Ohta, R. Pausch, C. Satzinger, B. Schreiber, M. Seeger, A. Smilga, D. Stoll, R. Weingärtner, S. Wörlen, K. Yazaki, and I. Zahed.

## References