--- title: 'Gross-Neveu Model: Theory and Applications' url: https://www.emergentmind.com/topics/gross-neveu-model type: topic --- # Gross-Neveu Model: Theory and Applications The Gross–Neveu model is a relativistic quantum field theory of interacting fermions with a local four-fermion interaction. Its canonical form in \(1+1\) dimensions contains \(N\) Dirac or Majorana fermions and a dimensionless coupling, for example \[ \mathcal L =i\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2, \] with repeated flavor indices summed. The model is asymptotically free, exhibits dimensional transmutation, dynamically generates a fermion mass, and possesses a discrete chiral symmetry that is spontaneously broken in the massless theory. It is also integrable in several important \(1+1\)-dimensional formulations, with exact Hartree–Fock, soliton, and scattering solutions. Generalizations include chiral, isospin, tensorial, anisotropic, lattice, higher-dimensional, and fermion-pairing variants. ## 1. Canonical formulation and fundamental properties ### Lagrangian and symmetries For \(N\) Majorana fermions transforming in the vector representation of \(O(N)\), a conventional two-dimensional Lagrangian is \[ \mathcal L =\frac{i}{2}\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2. \] For Dirac fermions, the corresponding internal symmetry is \(U(N)\), whereas the Majorana theory has \(O(N)\) symmetry. The coupling is dimensionless in two dimensions. A fermion mass term is excluded by the discrete chiral transformation \[ \psi_i\longrightarrow\gamma^5\psi_i, \qquad \bar\psi_i\psi_i\longrightarrow-\bar\psi_i\psi_i. \] The model is asymptotically free: the interaction decreases at short distances and becomes stronger toward the infrared. Dimensional transmutation generates a mass scale despite the absence of a bare mass. The dynamically generated mass is accompanied by spontaneous breaking of the discrete chiral symmetry, producing two degenerate vacua related by the chiral transformation. In the large-\(N\) formulation, the scalar order parameter is \[ S(x,t)=-g^2\langle\bar\psi\psi\rangle, \] and the homogeneous vacua are \(S=\pm m\) [1011.3835]. The ordinary model has infinitely many conserved quantities in its integrable \(1+1\)-dimensional realization. Its large-\(N\) dynamics can be formulated through relativistic Hartree–Fock or time-dependent Hartree–Fock equations. These equations admit exact kink, baryon, breather, multikink, and scattering configurations [1408.5506]. ### Auxiliary-field representation The four-fermion interaction can be linearized with an auxiliary scalar field \(\sigma\): \[ \mathcal L = \bar\psi i\!\not\!\partial\psi +\sigma\bar\psi\psi -\frac{N}{2g}\sigma^2. \] Integrating out \(\sigma\) reproduces the four-fermion interaction. At a critical coupling in dimensions \(23)\), analogous cubic currents are not classically conserved because additional flavor indices produce terms of the form \[ \psi_{lL}\psi_{lR}r^\pm_{ijk}. \] Anomaly-improved currents and contracted charge algebras can be constructed, but current conservation, closure on physical translation generators, and genuine supersymmetry are not established [1003.5040]. ### Higher-spin currents and conformal criticality At the critical \(U(N)\)-symmetric Gross–Neveu point in \(20. \] The ground-state entropy inherits the nonperturbative dependence of the dynamically generated mass. Although the Gross–Neveu coupling is asymptotically free, the entropy has its own renormalization-group flow and is not asymptotically free in the same sense [1512.00023]. ## 4. Finite temperature, density, and magnetic backgrounds ### Finite-temperature transitions In \(2+1\)-dimensional Gross–Neveu models with massless fermions, finite temperature produces a transition from a massive chirally ordered phase to a symmetric phase. A functional renormalization-group treatment retaining both fermionic and composite bosonic fluctuations finds a second-order transition in the two-dimensional Ising universality class for every finite \(N_f\) [1212.4624]. At mean-field level, the critical temperature is \[ T_c^{\rm MF} =\frac{m}{2\ln2} \simeq0.7213\,m. \] For finite \(N_f\), bosonic zero-mode fluctuations dominate sufficiently close to the transition. The Ginzburg region narrows approximately as \[ w_\beta(N_f)\sim \frac1{N_f}, \] so mean-field behavior is recovered outside an increasingly narrow critical region at large \(N_f\). The same FRG analysis identifies a local-order regime above \(T_c\), in which the flow is broken at intermediate renormalization scales but symmetric in the infrared. It also identifies a classical regime below \(T_c\), where the bosonic zero Matsubara mode dominates. The finite-temperature critical theory dimensionally reduces to a two-dimensional real scalar theory with \(\mathbb Z_2\) symmetry. ### Massive three-dimensional model A \(2+1\)-dimensional, one-component massive Gross–Neveu model has been studied as an effective theory using Matsubara sums, dimensional regularization, and generalized Epstein functions. Its thermal mass satisfies \[ m(T)=m_0+\Sigma(T), \] and the one-loop thermal self-energy gives \[ \Sigma_R(3,T,1) =\frac{\lambda_0m_0^2}{2\pi}f(T), \] with \[ f(T) = -\frac{T}{m_0\ln(1+e^{-m_0/T})}. \] A resummed four-point function develops a pole when \[ \lambda_0\Pi_R(3,T_c)=1. \] The authors interpret the vanishing thermal mass, the four-point singularity, and a nontrivial infrared fixed point as evidence for a second-order transition. In the strong-coupling limit, the critical temperature approaches \[ T_c\simeq0.87\,m_0. \] Because the model contains an explicit mass and \(N=1\), its transition is a thermal transition in an effective massive theory rather than the standard exact chiral-restoration transition of the massless large-\(N\) model [1009.4459]. ### Chemical potential and crystalline matter At finite density, the \(O(2N)\)-symmetric Gross–Neveu model develops distinct neutral and charged nonperturbative scales: \[ \Lambda_{\rm n} \propto h\exp\left[-\frac{\pi}{(N-a-1)g^2(h)}\right], \] and \[ \Lambda_{\rm c} \propto h\exp\left[-\frac{2\pi}{a\,g^2(h)}\right]. \] The first is associated with neutral low-energy scattering; the second is associated with charged particle-hole interactions near the charged Fermi points. Above the threshold \[ h_*^{\rm finite\,N}=\frac{m_a}{a}, \] where \(m_a\) is the mass of the rank-\(a\) bound state, charged bound states condense. At infinite \(N\), the ground state becomes an inhomogeneous crystal. At high density, \[ \langle\sigma(x)\rangle \simeq \Lambda_{\rm n} +2\Lambda_{\rm c}\sin(2hx). \] Thus the spatial average and oscillation amplitude are controlled by different dynamically generated scales. At finite \(N\), translational long-range order is replaced by quasi-long-range order, while the gapless phonon excitation remains [2505.23388]. ### Magnetic field and finite chemical potential In \(2+1\) dimensions, a homogeneous magnetic field quantizes the fermion spectrum into Landau levels, \[ E_n=\sqrt{\sigma^2+2|eB|n}. \] At large \(N_f\), the lowest Landau level can enhance chiral symmetry breaking at low chemical potential, while inverse magnetic catalysis and multiple transitions may appear at larger chemical potential. One-flavor overlap-fermion lattice simulations instead find that the magnetic field enhances chiral symmetry breaking throughout the broken phase, with a single transition whose critical chemical potential increases with \(B\). The mean-field multiple-transition structure and inverse magnetic catalysis are not observed at \(N_f=1\). Apparent first-order behavior at \(B=0\) remains subject to finite-volume and finite-flavor effects [2310.13455]. ## 5. Variants of the Gross–Neveu model ### Chiral Gross–Neveu and NJL\(_2\) The chiral Gross–Neveu or NJL\(_2\) model adds the pseudoscalar channel: \[ \mathcal L_{\rm NJL} = i\bar\psi\slashed{\partial}\psi +\frac{g^2}{2} \left[ (\bar\psi\psi)^2 +(\bar\psi i\gamma^5\psi)^2 \right]. \] Its order parameter is complex, \[ \Delta=S-iP, \] and the discrete chiral symmetry is enlarged to continuous chiral \(U(1)\). The vacuum manifold is a circle rather than two isolated points. At finite density, the preferred massless solution is the chiral spiral \[ \Delta(x)=Ae^{i\mu x}. \] The complex condensate is governed by the AKNS hierarchy, whereas the real Gross–Neveu condensate is governed by mKdV [1011.3835]. ### Isospin and chiral Heisenberg Gross–Neveu models The isoGN model contains an \(SU(2)\) isospin triplet: \[ \mathcal L_{\rm isoGN} = \bar\psi i\!\not\!\partial\psi +\frac{g^2}{2}(\bar\psi\vec\tau\psi)^2. \] Its order parameter is \[ \vec S=-g^2\langle\bar\psi\vec\tau\psi\rangle. \] At large \(N_c\), the vacuum manifold is \[ \mathcal M_{\rm vac}=S^2, \] with spontaneous breaking \[ SU(2)\longrightarrow U(1). \] The spectrum contains two massless scalar bosons and one massive scalar boson of mass \(2M\). A duality maps the model to a restricted isoNJL theory in which the isovector pseudoscalar field is constrained to a plane. This maps ordinary Gross–Neveu condensates and crystals into isoGN configurations with fixed isospin direction [2008.13119]. The \(O(2)_L\times O(2)_R\) chiral Gross–Neveu model has vacuum manifold \[ O(2)\simeq S^1\sqcup S^1. \] One connected component represents particle–hole or chiral condensation, while the other represents Cooper pairing under duality to the perfect Gross–Neveu model. Kinks interpolate between the two disconnected components, with mass \[ M_{\rm kink}=\frac{Nm}{2\pi}. \] The model is dual to a quantum version of the Zakharov–Mikhailov integrable spinor model and provides an exactly solvable setting for competing chiral and Cooper-pair order [2302.07660]. ### Pairing, self-duality, and the perfect Gross–Neveu model A one-chirality particle–hole transformation maps fermion–antifermion condensates into fermion–fermion Cooper-pair condensates. Applying this transformation to the ordinary GN and NJL models yields dual integrable pairing theories. Adding the GN interaction to its dual gives the self-dual GN model, containing both a real scalar condensate \(\mathcal S\) and a real Cooper condensate \(\mathcal B\). Its Hartree–Fock–Bogoliubov Hamiltonian can be block-diagonalized into two ordinary Gross–Neveu Hamiltonians with mean fields \[ S_I=\mathcal S-\mathcal B, \qquad S_{II}=\mathcal S+\mathcal B. \] Consequently, the self-dual GN model is exactly two independent \(O(N)\) Gross–Neveu models and is integrable. Its four vacua include chirally broken and superconducting configurations. Kinks can interpolate between these vacua, and some configurations are analogous to Josephson junctions [1408.5506]. The perfect Gross–Neveu model is the self-dual version of the NJL model. It contains complex particle–hole and Cooper-pair condensates and preserves the full Pauli–Gürsey symmetry of free massless fermions. Its general integrability is not proven in the earlier construction, although its classical interaction coincides with the Zakharov–Mikhailov model. The \(O(2)\) chiral Gross–Neveu analysis provides large-\(N\) evidence that the perfect model is quantum mechanically solvable [1408.5506; 2302.07660]. ### Tensorial Gross–Neveu models Tensorial generalizations replace the vector flavor index by a rank-three tensor, \[ \psi_{abc}, \qquad a,b,c=1,\ldots,N, \] with \(U(N)^3\) or \(O(N)^3\) symmetry. The theories retain two-dimensional four-fermion interactions and a discrete chiral symmetry while admitting melonic large-\(N\) expansions. For complex fermions, the interactions include double-trace and pillow invariants. For Majorana fermions, a non-bipartite tetrahedral invariant appears. In the Dirac theory, the large-\(N\) Schwinger–Dyson equation yields a mass gap \[ m^2 = \Lambda^2 \exp\left[ -\frac{\pi}{\lambda_0^S+3\lambda_1^S} \right]. \] A formally massless subspace exists when \[ \lambda_0^S+3\lambda_1^S=0, \] but it is infrared unstable and does not define an interacting conformal fixed point. The Majorana tetrahedral coupling produces a momentum-dependent self-energy and has beta function \[ \beta_2=\frac{3}{\pi^2}\lambda_2^3. \] It is marginally irrelevant in two dimensions. The large-\(N\) stable vacuum breaks \(U(N)^3\) at the saddle-point level, but the Coleman–Mermin–Wagner theorem restores continuous symmetry at finite \(N\); the apparent breaking is interpreted as a large-\(N\) artifact associated with a fluctuating Grassmannian sigma model [1710.10253]. ## 6. Higher dimensions, anisotropy, lattices, and numerical simulation ### \(2+1\)-dimensional and Lifshitz-like models The conventional \(2+1\)-dimensional Gross–Neveu model has an interacting ultraviolet fixed point and a finite-temperature transition in the two-dimensional Ising universality class. Functional renormalization-group calculations retain both fermionic and composite bosonic fluctuations and show that the Ginzburg region shrinks as \(1/N_f\), while the Ising critical regime persists for every finite \(N_f\) [1212.4624]. A Hořava–Lifshitz-like extension replaces the spatial Dirac operator by \[ (i\gamma^i\partial_i)^z \] under anisotropic scaling \[ t\to b^zt, \qquad x^i\to bx^i. \] For even \(z\), the gap equation becomes scaleless and produces no dynamical mass. For odd \(z\), a nonzero mass is generated for the appropriate sign of the coupling, dynamically breaking parity. At finite temperature, parity breaking disappears at a critical temperature. The auxiliary-field theory is renormalizable in the \(1/N\) expansion for every \(z\). For odd \(z\), coupling to a vector field produces an anisotropic Chern–Simons-like term; the Lorentz-invariant result is recovered at \(z=1\) [1612.05900]. A recent analysis of the \(2+1\)-dimensional canonical model rewrites the interaction in terms of singlet, symmetric-tensor, and adjoint nematic channels of \(SO(2N)\). The conventional Gross–Neveu transition at \(g_{c1}\sim1/N\) is a genuine continuous transition that preserves \(SO(2N)\) and breaks a discrete Ising symmetry. Mean-field theory suggests a second transition at \(g_{c2}\sim -Ng_{c1}\), with \[ SO(2N)\longrightarrow SO(N)\times SO(N). \] A continuation through an auxiliary \(L>1\) theory shows that the symmetric-tensor critical point loses its divergent susceptibility as \(L\to1\), becoming equivalent to the Gaussian fixed point. The resulting first-order interpretation is physically inferred rather than rigorously established [2403.09627]. ### Gross–Neveu–Wilson models and correlated topological phases A Wilson discretization in \(1+1\) dimensions produces the Gross–Neveu–Wilson model. Its phase diagram contains: - a trivial symmetry-preserving insulator; - a symmetry-preserving topological insulator in the \(\mathsf{BDI}\) class; - a parity-breaking Aoki phase. The Wilson masses at the two lattice Dirac points are \[ \widetilde m_0=m, \qquad \widetilde m_\pi=m+\frac2a. \] The noninteracting Zak phase is \[ \varphi_{\rm Zak} = \frac{N\pi}{2} \left[ \operatorname{sgn}(\widetilde m_\pi) -\operatorname{sgn}(\widetilde m_0) \right]. \] Interactions renormalize these masses and generate a correlated symmetry-protected topological phase. At strong coupling, a pseudoscalar condensate \[ \Pi_0=\langle\bar\Psi\gamma^5\Psi\rangle \] defines an Aoki phase. The trivial, topological, and Aoki regions meet at a tricritical point. For \(N=1\), matrix-product-state calculations find Ising critical behavior at the Aoki boundaries and Dirac criticality at the weak-coupling trivial–topological transition [1807.03202]. Borici–Creutz fermions provide another lattice realization. Their preferred matrix \(\Gamma\) requires both scalar and \(\Gamma\)-direction condensates. Strong-coupling analysis and hybrid Monte Carlo simulations find a continuous chiral transition, with the scalar condensate vanishing continuously at the critical boundary [1409.7999]. ### Quantum and tensor-network simulation The Gross–Neveu model maps naturally to qubit Hamiltonians after spatial discretization and a Jordan–Wigner transformation. In a staggered formulation, \[ H^{(N)} = \sum_n \left[ i\sum_a\chi^{a\dagger}(n) \bigl(\chi^a(n+1)-\chi^a(n-1)\bigr) +m(-1)^n\chi^{a\dagger}(n)\chi^a(n) +G^2\left(\sum_a\chi^{a\dagger}(n)\chi^a(n)\right)^2 \right]. \] The interaction becomes local same-site \(ZZ\) coupling between different flavors. For \(N=2\), the model is related to the one-dimensional Hubbard model. Small systems have been studied with exact diagonalization, DMRG, TEBD, VQE, and hardware experiments on IBM and Quantinuum platforms. Current demonstrations focus on short-time evolution, local densities, wave-packet propagation, and scattering protocols rather than precision continuum observables [2208.05906]. A separate MPS study of a two-flavor massive lattice model uses a Wilson term to remove fermion doubling and examines one-, two-, and four-particle real-time evolution. Strong coupling and sufficiently small mass produce synchronized constituent dynamics that are interpreted as bound-state-like, although no bound-state mass or scattering pole is extracted [2011.07538]. ### Conceptual scope and limitations Several distinctions are essential when comparing Gross–Neveu results: - Large-\(N\) spontaneous breaking of continuous symmetries in two dimensions is not a finite-\(N\) statement because of the Coleman–Mermin–Wagner theorem. - A formal supersymmetry-like charge algebra does not establish genuine supersymmetry without conservation and closure on physical generators. - Mean-field phase diagrams can differ qualitatively from finite-\(N\) lattice results, particularly in magnetic fields and at finite chemical potential. - Crystalline condensates are exact large-\(N\) solutions, but translational long-range order is altered by infrared fluctuations at finite \(N\). - Quantum-simulation results generally use finite lattices, short times, Trotter approximations, and classical tensor-network benchmarks; they do not by themselves establish continuum mass gaps or scattering amplitudes. - The rigorous Polchinski-flow construction applies to the massive two-dimensional model at sufficiently small positive coupling and does not prove dynamical mass generation in the massless theory [2403.18562]. The Gross–Neveu model therefore occupies several interconnected roles: an asymptotically free and dynamically massive quantum field theory; an exactly solvable integrable system; a laboratory for discrete and continuous symmetry breaking; a prototype for finite-density crystals and solitons; a conformal theory with higher-spin data; a field theory of correlated topological phases; and a benchmark for tensor-network and quantum computation. Source: https://www.emergentmind.com/topics/gross-neveu-model