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, characterized by its asymptotic freedom, dimensional transmutation, dynamic fermion mass generation, and discrete chiral symmetry.
• Applications and significance include its role as an integrable system with exact solutions, a prototype for studying symmetry breaking, and a benchmark for quantum simulations and tensor-network computations.
• Variants such as the chiral and isospin Gross-Neuveu models, as well as higher-dimensional and tensorial generalizations, extend its utility in exploring various quantum field theory phenomena, and interacting field theory diagnostics.

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$0 formulation, the scalar order parameter is

$N$1

and the homogeneous vacua are $N$2 (Basar et al., 2010).

The ordinary model has infinitely many conserved quantities in its integrable $N$3-dimensional realization. Its large-$N$4 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 (Thies, 2014).

Auxiliary-field representation

The four-fermion interaction can be linearized with an auxiliary scalar field $N$5: $N$6 Integrating out $N$7 reproduces the four-fermion interaction. At a critical coupling in dimensions $N$8, this formulation defines a conformal Gross–Neveu fixed point with a nontrivial Wilson–Fisher-type structure (Manashov et al., 2016).

In $N$9 dimensions the four-fermion coupling has negative canonical mass dimension, but the model possesses an interacting non-Gaussian ultraviolet fixed point in the large-$$\mathcal L =i\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2,$$0 formulation. A renormalized trajectory can be specified by fixing one physical scale, such as the zero-temperature fermion mass (Scherer et al., 2012).

Dynamical mass generation

At leading large $$\mathcal L =i\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2,$$1, the fermions propagate in a self-consistent scalar background: $$\mathcal L =i\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2,$$2 For a homogeneous condensate, the gap equation produces a dynamically generated fermion mass. In $$\mathcal L =i\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2,$$3 dimensions, the renormalized effective potential has two nonzero minima, and the mass depends nonperturbatively on the coupling through an exponential of the inverse coupling squared (Cotler et al., 2015). The same mechanism appears in lattice Gross–Neveu–Wilson models, where the dynamically generated mass behaves as

$$\mathcal L =i\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2,$$4

in the weak-coupling large-$$\mathcal L =i\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2,$$5 regime (Bermudez et al., 2018).

A rigorous construction of the massive Euclidean Gross–Neveu model in infinite volume has been obtained using the Polchinski flow equation and a Banach fixed-point theorem. This construction yields nonzero truncated four-point Schwinger functions and proves that the limiting theory is interacting rather than Gaussian (Duch, 2024).

2. Large-$$\mathcal L =i\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2,$$6 dynamics, condensates, and integrability

Hartree–Fock and nonlinear Dirac equations

In the ’t Hooft limit,

$$\mathcal L =i\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2,$$7

the fermions propagate in a classical self-consistent background. The Hartree–Fock equation is

$$\mathcal L =i\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2,$$8

with

$$\mathcal L =i\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2,$$9

A particularly useful class consists of “Type I” solutions, for which each occupied mode separately satisfies a nonlinear Dirac equation,

$1+1$0

The Hartree–Fock potential is reconstructed by summing the mode bilinears. The static and time-dependent solutions discussed in the integrability literature belong to this class (Basar et al., 2010).

Static crystalline condensates

Restricting the gap equation to homogeneous condensates gives an incomplete finite-temperature and finite-density phase diagram. The massless $1+1$1-dimensional model possesses a self-consistent inhomogeneous crystal,

$1+1$2

where $1+1$3 sets the scale and $1+1$4 controls the period and shape.

The principal phases are:

• the symmetric phase, $1+1$5;
• the homogeneous massive phase, $1+1$6;
• the inhomogeneous crystalline phase, $1+1$7 and $1+1$8.

The parameters are determined by minimizing the grand potential at temperature $1+1$9 and chemical potential $N$0. The crystalline phase is a thermodynamic phase rather than merely a formal solution of the self-consistency equations (Basar et al., 2010).

For a real static condensate, the functional gap equation reduces to

$N$1

with first integral

$N$2

The elliptic crystal is a finite-gap solution of this equation. The conserved functionals appearing in the Ginzburg–Landau expansion are the conserved quantities of the modified Korteweg–de Vries hierarchy. For the chiral Gross–Neveu or NJL$N$3 model, whose order parameter is complex,

$N$4

the corresponding hierarchy is AKNS rather than mKdV (Basar et al., 2010).

Time-dependent solutions and geometric integrability

Time-dependent Hartree–Fock solutions include boosted kinks, kink–antikink scattering, breathers, and multiparticle baryon configurations. Bäcklund transformations generate more complicated scattering solutions. The time-dependent condensate is related to the Sinh–Gordon equation. For a Minkowski surface,

$N$5

The generalized Weierstrass representation identifies the Gross–Neveu nonlinear Dirac equation with the spinor representation of constant-mean-curvature surfaces in $N$6. In this correspondence,

$N$7

and the Gauss–Codazzi compatibility condition becomes the Sinh–Gordon equation. Constant-mean-curvature surfaces in $N$8 are related through a curvature tradeoff to minimal string worldsheets in $N$9, which obey the same integrable Sinh–Gordon system (Basar et al., 2010).

Explicit soliton constructions

Classical nonzero-background solitons can be generated through a Cauchy-matrix construction based on a rank-one Sylvester equation. Under the Gross–Neveu reduction, the scalar condensate is

$O(N)$0

where

$O(N)$1

The associated tau functions satisfy a Sinh–Gordon equation. For a one-soliton solution,

$O(N)$2

and the fermion-density-like bilinear contains a localized depletion on a nonzero background. These configurations are therefore dark solitons in the integrable-systems sense and kinks in field-theoretic terminology (Vekslerchik, 2021).

3. Quantum currents, higher-spin structure, and entanglement

Anomaly-improved currents

In the $O(N)$3 Gross–Neveu model, the classical chiral cubic currents

$O(N)$4

are conserved because repeated chiral fermions vanish classically. At the quantum level, anomalies modify the current. The conserved current is taken to be

$O(N)$5

where $O(N)$6 is an anomaly contribution containing mixed left- and right-moving cubic operators.

The corresponding charges satisfy

$O(N)$7

The algebra therefore has the form expected of a supersymmetry algebra, with formally vanishing central charge. The interpretation is algebraic rather than a proof that the original Gross–Neveu model possesses ordinary supersymmetry. For $O(N)$8, analogous cubic currents are not classically conserved because additional flavor indices produce terms of the form

$O(N)$9

Anomaly-improved currents and contracted charge algebras can be constructed, but current conservation, closure on physical translation generators, and genuine supersymmetry are not established (Morales, 2010).

Higher-spin currents and conformal criticality

At the critical $$\mathcal L =\frac{i}{2}\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2.$$0-symmetric Gross–Neveu point in $$\mathcal L =\frac{i}{2}\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2.$$1, symmetric traceless fermion bilinears acquire anomalous dimensions. The leading non-singlet anomalous dimension is

$$\mathcal L =\frac{i}{2}\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2.$$2

where $$\mathcal L =\frac{i}{2}\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2.$$3. Singlet currents differ for even spin because singlet operators can couple to fermion loops. For odd spin, singlet and non-singlet anomalous dimensions coincide.

The spin-one flavor current and spin-two stress tensor remain protected: $$\mathcal L =\frac{i}{2}\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2.$$4 At order $$\mathcal L =\frac{i}{2}\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2.$$5, the conservation of these currents persists. In $$\mathcal L =\frac{i}{2}\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2.$$6, the anomalous dimensions provide precision data for the conjectured duality between the critical Gross–Neveu model and Type-B higher-spin theory in $$\mathcal L =\frac{i}{2}\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2.$$7. The bulk higher-spin mass shifts are determined by the boundary anomalous dimensions, while the graviton remains massless because the stress tensor is conserved (Manashov et al., 2016).

Entanglement entropy

A Grassmann Gaussian variational state gives the half-line entropy of $$\mathcal L =\frac{i}{2}\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2.$$8 massive fermions: $$\mathcal L =\frac{i}{2}\bar\psi_i\gamma^\mu\partial_\mu\psi_i +g(\bar\psi_i\psi_i)^2.$$9 The mass is determined self-consistently by the variational gap equation. At fixed renormalized auxiliary field, the entropy decreases monotonically and is convex as the coupling increases: $U(N)$0 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 (Cotler et al., 2015).

4. Finite temperature, density, and magnetic backgrounds

Finite-temperature transitions

In $U(N)$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 $U(N)$2 (Scherer et al., 2012).

At mean-field level, the critical temperature is

$U(N)$3

For finite $U(N)$4, bosonic zero-mode fluctuations dominate sufficiently close to the transition. The Ginzburg region narrows approximately as

$U(N)$5

so mean-field behavior is recovered outside an increasingly narrow critical region at large $U(N)$6.

The same FRG analysis identifies a local-order regime above $U(N)$7, in which the flow is broken at intermediate renormalization scales but symmetric in the infrared. It also identifies a classical regime below $U(N)$8, where the bosonic zero Matsubara mode dominates. The finite-temperature critical theory dimensionally reduces to a two-dimensional real scalar theory with $U(N)$9 symmetry.

Massive three-dimensional model

A $O(N)$0-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

$O(N)$1

and the one-loop thermal self-energy gives

$O(N)$2

with

$O(N)$3

A resummed four-point function develops a pole when

$O(N)$4

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

$O(N)$5

Because the model contains an explicit mass and $O(N)$6, its transition is a thermal transition in an effective massive theory rather than the standard exact chiral-restoration transition of the massless large-$O(N)$7 model (Khanna et al., 2010).

Chemical potential and crystalline matter

At finite density, the $O(N)$8-symmetric Gross–Neveu model develops distinct neutral and charged nonperturbative scales: $O(N)$9 and

$$\psi_i\longrightarrow\gamma^5\psi_i, \qquad \bar\psi_i\psi_i\longrightarrow-\bar\psi_i\psi_i.$$0

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

$$\psi_i\longrightarrow\gamma^5\psi_i, \qquad \bar\psi_i\psi_i\longrightarrow-\bar\psi_i\psi_i.$$1

where $$\psi_i\longrightarrow\gamma^5\psi_i, \qquad \bar\psi_i\psi_i\longrightarrow-\bar\psi_i\psi_i.$$2 is the mass of the rank-$$\psi_i\longrightarrow\gamma^5\psi_i, \qquad \bar\psi_i\psi_i\longrightarrow-\bar\psi_i\psi_i.$$3 bound state, charged bound states condense. At infinite $$\psi_i\longrightarrow\gamma^5\psi_i, \qquad \bar\psi_i\psi_i\longrightarrow-\bar\psi_i\psi_i.$$4, the ground state becomes an inhomogeneous crystal. At high density,

$$\psi_i\longrightarrow\gamma^5\psi_i, \qquad \bar\psi_i\psi_i\longrightarrow-\bar\psi_i\psi_i.$$5

Thus the spatial average and oscillation amplitude are controlled by different dynamically generated scales. At finite $$\psi_i\longrightarrow\gamma^5\psi_i, \qquad \bar\psi_i\psi_i\longrightarrow-\bar\psi_i\psi_i.$$6, translational long-range order is replaced by quasi-long-range order, while the gapless phonon excitation remains (Benini et al., 29 May 2025).

Magnetic field and finite chemical potential

In $$\psi_i\longrightarrow\gamma^5\psi_i, \qquad \bar\psi_i\psi_i\longrightarrow-\bar\psi_i\psi_i.$$7 dimensions, a homogeneous magnetic field quantizes the fermion spectrum into Landau levels,

$$\psi_i\longrightarrow\gamma^5\psi_i, \qquad \bar\psi_i\psi_i\longrightarrow-\bar\psi_i\psi_i.$$8

At large $$\psi_i\longrightarrow\gamma^5\psi_i, \qquad \bar\psi_i\psi_i\longrightarrow-\bar\psi_i\psi_i.$$9, 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 $N$00. The mean-field multiple-transition structure and inverse magnetic catalysis are not observed at $N$01. Apparent first-order behavior at $N$02 remains subject to finite-volume and finite-flavor effects (Mandl et al., 2023).

5. Variants of the Gross–Neveu model

Chiral Gross–Neveu and NJL$N$03

The chiral Gross–Neveu or NJL$N$04 model adds the pseudoscalar channel: $N$05 Its order parameter is complex,

$N$06

and the discrete chiral symmetry is enlarged to continuous chiral $N$07. The vacuum manifold is a circle rather than two isolated points. At finite density, the preferred massless solution is the chiral spiral

$N$08

The complex condensate is governed by the AKNS hierarchy, whereas the real Gross–Neveu condensate is governed by mKdV (Basar et al., 2010).

Isospin and chiral Heisenberg Gross–Neveu models

The isoGN model contains an $N$09 isospin triplet: $N$10 Its order parameter is

$N$11

At large $N$12, the vacuum manifold is

$N$13

with spontaneous breaking

$N$14

The spectrum contains two massless scalar bosons and one massive scalar boson of mass $N$15. 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 (Thies, 2020).

The $N$16 chiral Gross–Neveu model has vacuum manifold

$N$17

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

$N$18

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 (Thies, 2023).

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 $N$19 and a real Cooper condensate $N$20.

Its Hartree–Fock–Bogoliubov Hamiltonian can be block-diagonalized into two ordinary Gross–Neveu Hamiltonians with mean fields

$N$21

Consequently, the self-dual GN model is exactly two independent $N$22 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 (Thies, 2014).

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 $N$23 chiral Gross–Neveu analysis provides large-$N$24 evidence that the perfect model is quantum mechanically solvable (Thies, 2014, Thies, 2023).

Tensorial Gross–Neveu models

Tensorial generalizations replace the vector flavor index by a rank-three tensor,

$N$25

with $N$26 or $N$27 symmetry. The theories retain two-dimensional four-fermion interactions and a discrete chiral symmetry while admitting melonic large-$N$28 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$29 Schwinger–Dyson equation yields a mass gap

$N$30

A formally massless subspace exists when

$N$31

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

$N$32

It is marginally irrelevant in two dimensions. The large-$N$33 stable vacuum breaks $N$34 at the saddle-point level, but the Coleman–Mermin–Wagner theorem restores continuous symmetry at finite $N$35; the apparent breaking is interpreted as a large-$N$36 artifact associated with a fluctuating Grassmannian sigma model (Benedetti et al., 2017).

6. Higher dimensions, anisotropy, lattices, and numerical simulation

$N$37-dimensional and Lifshitz-like models

The conventional $N$38-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 $N$39, while the Ising critical regime persists for every finite $N$40 (Scherer et al., 2012).

A Hořava–Lifshitz-like extension replaces the spatial Dirac operator by

$N$41

under anisotropic scaling

$N$42

For even $N$43, the gap equation becomes scaleless and produces no dynamical mass. For odd $N$44, 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 $N$45 expansion for every $N$46. For odd $N$47, coupling to a vector field produces an anisotropic Chern–Simons-like term; the Lorentz-invariant result is recovered at $N$48 (Lima et al., 2016).

A recent analysis of the $N$49-dimensional canonical model rewrites the interaction in terms of singlet, symmetric-tensor, and adjoint nematic channels of $N$50. The conventional Gross–Neveu transition at $N$51 is a genuine continuous transition that preserves $N$52 and breaks a discrete Ising symmetry. Mean-field theory suggests a second transition at $N$53, with

$N$54

A continuation through an auxiliary $N$55 theory shows that the symmetric-tensor critical point loses its divergent susceptibility as $N$56, becoming equivalent to the Gaussian fixed point. The resulting first-order interpretation is physically inferred rather than rigorously established (Han et al., 2024).

Gross–Neveu–Wilson models and correlated topological phases

A Wilson discretization in $N$57 dimensions produces the Gross–Neveu–Wilson model. Its phase diagram contains:

• a trivial symmetry-preserving insulator;
• a symmetry-preserving topological insulator in the $N$58 class;
• a parity-breaking Aoki phase.

The Wilson masses at the two lattice Dirac points are

$N$59

The noninteracting Zak phase is

$N$60

Interactions renormalize these masses and generate a correlated symmetry-protected topological phase. At strong coupling, a pseudoscalar condensate

$N$61

defines an Aoki phase. The trivial, topological, and Aoki regions meet at a tricritical point. For $N$62, matrix-product-state calculations find Ising critical behavior at the Aoki boundaries and Dirac criticality at the weak-coupling trivial–topological transition (Bermudez et al., 2018).

Borici–Creutz fermions provide another lattice realization. Their preferred matrix $N$63 requires both scalar and $N$64-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 (Goswami et al., 2014).

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,

$N$65

The interaction becomes local same-site $N$66 coupling between different flavors. For $N$67, 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 (Asaduzzaman et al., 2022).

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 (Li et al., 2020).

Conceptual scope and limitations

Several distinctions are essential when comparing Gross–Neveu results:

• Large-$N$68 spontaneous breaking of continuous symmetries in two dimensions is not a finite-$N$69 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$70 lattice results, particularly in magnetic fields and at finite chemical potential.
• Crystalline condensates are exact large-$N$71 solutions, but translational long-range order is altered by infrared fluctuations at finite $N$72.
• 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 (Duch, 2024).

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.