Gross-Neveu Equation: Integrability & Phase Transitions

• Gross-Neveu equation is a nonlinear Dirac system that models fermions with local four-fermion interactions in 1+1 dimensions and spontaneous chiral symmetry breaking.
• It provides exact solutions for solitonic excitations, crystalline condensates, and finite-density phase transitions using techniques like the thermodynamic Bethe Ansatz and finite-gap methods.
• The model’s applications span quantum field theory and condensed matter, offering insights through mean-field approximations, integrable hierarchies, and inverse scattering techniques.

The Gross-Neveu equation, in the context of quantum field theory and integrable systems, refers to the nonlinear Dirac equation governing the dynamics of fermions with local four-fermion interactions in 1+1 dimensions (“Gross–Neveu model”). This integrable model plays a central role in many-body theory, providing exact results on symmetry breaking, solitonic excitations, finite-density phase transitions, and the emergence of nonlinear structures such as crystalline condensates and kink solutions.

1. Definition and Fundamental Equations

The Gross-Neveu model is defined in 1+1-dimensional spacetime for an $N$-component Dirac spinor $\psi_i(x)$ coupled via a local four-fermion interaction. The Lagrangian density is conventionally written as

$$\mathcal{L}[\psi, \sigma] = \bar{\psi}_i\, (i\gamma^\mu\partial_\mu + \mu\gamma^0 - \sigma)\,\psi_i - \frac{N}{2\lambda}\,\sigma^2$$

where $i = 1, \dots, N$ labels fermion “flavors,” $\mu$ is the chemical potential for the $U(1)$ baryon number, $\lambda$ is the ’t Hooft coupling, and $\sigma(x)$ is a Hubbard–Stratonovich scalar field (Melin et al., 2024).

The Euler-Lagrange equations yield:

• The nonlinear Dirac (Gross-Neveu) equation:

$$(i\gamma^\mu\partial_\mu + \mu\gamma^0 - \sigma(x))\,\psi_i(x) = 0$$

• The self-consistency (gap) condition (“Hartree-Fock equation”):

$$\sigma(x) = -\frac{\lambda}{N}\sum_i\,\bar\psi_i(x)\psi_i(x)$$

This system describes interacting Dirac fields with dynamical mass generation, asymptotic freedom, and spontaneous breaking of discrete chiral symmetry—$\psi \rightarrow \gamma^5\psi$, $\sigma \rightarrow -\sigma$—at vanishing chemical potential. In the chiral version, the condensate is complex and the gap equation generalizes accordingly (Basar et al., 2010, 0803.1501).

2. Integrability, Bethe Ansatz, and Thermodynamics

The Gross-Neveu equation is integrable in 1+1 dimensions. Quantization proceeds via the thermodynamic Bethe Ansatz (TBA), which yields the exact spectrum:

• Elementary “vector” fermions with mass $m_a$,
• Topological “kinks” (spinors) with mass $m_s$ and baryon charge $B$ ($B = N/2$ per kink).

At finite density, the ground state is a filled Fermi sea of kinks. The Bethe Ansatz equations for the rapidity density $\rho(\theta)$ and dressed energy $\epsilon(\theta)$ are (Melin et al., 2024): $$\begin{aligned} \rho(\theta) - \int_{-B}^{B} K(\theta - \eta)\,\rho(\eta)\,d\eta &= \cosh\theta \ \epsilon(\theta) - \int_{-B}^{B} K(\theta - \eta)\,\epsilon(\eta)\,d\eta &= \cosh\theta - (\mu/m_s) \end{aligned}$$ Here $B$ defines the Fermi boundary and $K(\theta)$ derives from the exact kink–kink $S$-matrix. The large-$N$ limit simplifies these to singular integral equations involving elliptic integrals, resulting in an explicit spectral curve and ground-state potential. Excitations include phonon-like holes and particle branches corresponding to anti-kinks and fermions, with an underlying genus-one algebraic curve reflecting the band structure of the spectrum.

3. Mean-Field, Finite-Gap Solutions, and Integrable Hierarchies

The mean-field (Hartree–Fock) approximation for static configurations—the self-consistent Dirac equation with a spatially periodic condensate—exhibits a direct link to classical soliton theory. The gap equation admits solutions in the form of finite-gap (cnoidal) potentials: $$\langle\sigma(x)\rangle = m \sqrt{k} \,\mathrm{sn}(m x \sqrt{k}\,|\,k)$$ where $\mathrm{sn}$ is the Jacobi elliptic sine and $k$ the elliptic modulus (Melin et al., 2024, Basar et al., 2010). The mean-field spectral problem for the Dirac Hamiltonian maps precisely to the KdV/mKdV integrable hierarchies. The resulting Dirac spectrum exhibits band structures with two gaps, described by a genus-one (elliptic) spectral curve: $$y^2 = (\epsilon^2 - \epsilon_+^2)(\epsilon^2 - \epsilon_-^2)$$ with band edges $\epsilon_\pm = m\,(\cosh B,\, \sinh B)$. The full hierarchy of integrals of motion matches those in the mKdV theory, and the stationary gap equation reduces to an ordinary nonlinear Schrödinger equation in imaginary time (Basar et al., 2010).

4. Nonlinear PDE Structure and Global Solutions

The Gross-Neveu equation, in component form, is a nonlinear Dirac system with cubic nonlinearity: $$i\,\gamma^\mu\partial_\mu \psi - m\psi + g^2(\bar\psi\psi)\psi = 0$$ In two-component notation, this reads: $$\begin{aligned} &-i\,(u_t + u_x) = -\,m\,v + N_1(u, v) \ & i\,(v_t - v_x) = -\,m\,u + N_2(u, v) \end{aligned}$$ with $N_1(u, v), N_2(u, v)$ constructed from the quartic potential $W(u, v) = \frac{1}{4}(\bar u v + u \bar v)^2$ (Zhang et al., 2014). For initial data in $L^2(\mathbb{R}) \times L^2(\mathbb{R})$, there exists a unique strong global solution in $$C([0,\infty); L^2(\mathbb{R}) \times L^2(\mathbb{R}))$$. This result holds without any small-data assumption, demonstrating that the nonlinear Dirac (Gross-Neveu) equation has robust global well-posedness in one spatial dimension (Zhang et al., 2014).

5. Crystalline Phases, Peierls Instability, and Phase Structure

The Gross-Neveu equation underlies a true second-order crystalline (Peierls-type) quantum phase transition at large $N$ and high baryon density. For chemical potential $\mu < \mu_c = 2m/\pi$, the vacuum is uniform ($\sigma = \pm m$). When $\mu$ exceeds $\mu_c$, it becomes energetically favorable to populate kink states, and the system transitions to a spatially inhomogeneous, periodic (finite-gap) condensate—a “chiral crystal” (Melin et al., 2024). The nature of the transition is governed by the closing of spectral gaps and the opening of a Peierls gap at the Fermi surface. Mean-field and Bethe Ansatz calculations show complete quantitative agreement in the large-$N$ limit, both yielding identical mass–charge relations, ground-state energies, and spectral data.

6. Multi-Flavor, Twisted-Kink, and Time-Dependent Solutions

The chiral and multiflavor extensions of the Gross-Neveu equation introduce a complex or matrix-valued condensate and greatly enrich the solution space. In the multiflavor chiral Gross-Neveu model (continuous chiral symmetry $U_L(N_f) \times U_R(N_f)$), the time-dependent Hartree-Fock (TDHF) equations admit soliton, twisted-kink, and breather solutions using inverse scattering. Analytical multi-kink and breather solutions can be constructed via a Zakharov–Shabat formalism and reduced to compact expressions where all $(x, t)$-dependence is explicit, and constant flavor-twist matrices encode scattering and internal rotations (Thies, 2021).

Crystalline condensates, chiral spirals, twisted kinks, and kink-antikink scattering all occur as special limits of the general nonlinear Dirac system. Time-dependent solutions are often governed by additional integrable structures, such as the Sinh-Gordon equation for certain sectors (Basar et al., 2010).

7. Connections to Integrable Geometry and Physical Applications

The deep integrability of the Gross-Neveu equation connects it to a wide range of mathematical structures—KdV/mKdV hierarchies (stationary gap equations), finite-gap (algebro-geometric) theory, and the embedding of nonlinear Dirac systems into geometric representation via Weierstrass spinor techniques. In some cases, explicit solutions correspond to minimal surfaces or classical string worldsheets in $AdS_3$ (Basar et al., 2010).

Physically, the Gross-Neveu equation and model provide canonical solvable examples for the study of dynamical mass generation, asymptotic freedom, and quantum phase transitions in lower dimensions. Applications extend to condensed matter physics (e.g., Peierls instability, LOFF phases in superconductivity, and one-dimensional polymers), and provide models for baryon-like soliton excitations and the study of exact many-body spectra in quantum field theory (0803.1501, Thies, 2021).