Two-Component Gross-Neveu Model

• The two-component Gross-Neveu model is a family of fermionic quantum field theories defined in 1+1 and 2+1 dimensions using minimal (Dirac or Majorana) spinor representations.
• It elucidates key phenomena such as dynamical mass generation, parity and chiral symmetry breaking, and a rich phase structure through various analytical and numerical methods.
• Modern studies extend the model via large-N analyses, CJT effective action, lattice formulations, and tensor-network simulations to capture bound states, solitons, and universality classes.

Searching arXiv for recent and foundational papers on the two-component Gross–Neveu model. The two-component Gross–Neveu model denotes a family of interacting fermionic quantum field theories in which the elementary fermion field is represented by a minimal spinor in low spacetime dimension, most commonly a two-component Dirac or Majorana field in $1+1$ or $2+1$ dimensions. In its standard form, the model couples $N$ fermion flavors through a local scalar four-fermion interaction $(\bar\psi\psi)^2$, and is studied both as a paradigmatic asymptotically free theory and as an effective description of symmetry breaking, dynamical mass generation, and bound-state formation. In $2+1$ dimensions, the two-component representation is especially consequential because the fermion mass term is parity odd, so dynamical mass generation is directly tied to parity breaking (Lima et al., 2016). In $1+1$ dimensions, the same model and its chiral or flavor-extended variants support exact large-$N$ analyses, lattice worldline and tensor-network formulations, and nontrivial phase structure at finite temperature and density (Maillart et al., 2011, Li et al., 2020, Thies, 19 Jun 2026).

1. Definition, field content, and representations

The standard $1+1$-dimensional Gross–Neveu model with $N$ Dirac flavors is

$$\mathcal{L} = \sum_{i=1}^N \bar\psi_i (\gamma_\mu \partial_\mu + m)\psi_i - \frac{g^2}{2}\left(\sum_{i=1}^N \bar\psi_i\psi_i\right)^2,$$

with each Dirac fermion represented by a two-component spinor (Maillart et al., 2011). A two-flavor massive version used in real-time simulation is

$2+1$0

where $2+1$1, $2+1$2 is a two-component Dirac spinor, and the model has a global $2+1$3 flavor symmetry (Li et al., 2020). In that formulation, the spinor is decomposed as

$2+1$4

so “two-component” refers both to the spinor structure $2+1$5 and, in the two-flavor case, to the presence of two flavor species (Li et al., 2020).

In $2+1$6 dimensions, the relativistic Gross–Neveu model for $2+1$7 flavors of two-component fermions is written as

$2+1$8

with $2+1$9 gamma matrices obeying

$N$0

(Lima et al., 2016). In that setting, a two-component spinor is an irreducible representation of the Lorentz group, and the scalar bilinear $N$1 is parity odd, so a dynamically generated mass implies dynamical parity breaking (Lima et al., 2016).

A recurrent reformulation replaces Dirac fermions by Majorana fields. In $N$2 dimensions, one two-component Dirac fermion equals two two-component Majorana fermions, and the $N$3-flavor model can therefore be rewritten in terms of $N$4 Majorana fields with manifest $N$5 symmetry (Maillart et al., 2011). In $N$6 dimensions, a $N$7-component Dirac field can likewise be decomposed into real Majorana components, exposing an $N$8 structure that organizes mass bilinears into singlet, symmetric-tensor, and adjoint channels (Han et al., 2024, Han et al., 2024).

2. Symmetry structure and mass terms

The defining interaction channel of the canonical model is the scalar four-fermion term $N$9, but the symmetry content depends strongly on dimension and representation. In $(\bar\psi\psi)^2$0 dimensions with two-component fermions, parity may be defined by

$(\bar\psi\psi)^2$1

under which $(\bar\psi\psi)^2$2 changes sign (Lima et al., 2016). A nonzero fermion mass therefore breaks parity dynamically. The same parity sensitivity persists in the Hořava–Lifshitz-like generalization, where the spatial kinetic operator is $(\bar\psi\psi)^2$3: for odd $(\bar\psi\psi)^2$4 the kinetic term is parity even and the mass term parity odd, whereas for even $(\bar\psi\psi)^2$5 the parity properties differ and no dynamical mass is generated (Lima et al., 2016).

In $(\bar\psi\psi)^2$6 dimensions, the conventional model has a global $(\bar\psi\psi)^2$7 charge symmetry and $(\bar\psi\psi)^2$8 flavor symmetry, and, in chiral variants, additional axial structure. The generalized Gross–Neveu model with both scalar and pseudoscalar channels,

$(\bar\psi\psi)^2$9

has a continuous axial symmetry on the line $2+1$0, where the model becomes the chiral Gross–Neveu model (Roose et al., 2021). In that case, the bilinears $2+1$1 and $2+1$2 form a doublet under axial rotation, and bosonization shows that the infrared theory contains a compact boson whose shift symmetry realizes the continuous chiral symmetry (Roose et al., 2021). The same work emphasizes a mixed ’t Hooft anomaly between the continuous chiral symmetry and charge conservation, implying the existence of a massless mode in the continuum theory (Roose et al., 2021).

The O$2+1$3 and $2+1$4 reformulations refine this symmetry perspective. In the $2+1$5-dimensional Wilson-fermion loop formulation, the continuum theory

$2+1$6

makes the $2+1$7 invariance manifest (Maillart et al., 2011). In $2+1$8 dimensions, the $2+1$9 tensor program identifies the canonical Gross–Neveu model as the $1+1$0 slice of a broader theory where a real symmetric traceless tensor order parameter unifies all Lorentz-invariant mass-gap orders for $1+1$1 two-component Dirac fermions except the $1+1$2-singlet anomalous quantum Hall state (Han et al., 2024). This suggests that the “two-component Gross–Neveu model” is not a single universality class once tensor and adjoint channels are admitted, but rather the singlet member of a larger family (Han et al., 2024, Han et al., 2024).

3. Dynamical mass generation and phase structure

The standard nonperturbative phenomenon associated with the model is dynamical mass generation. In the Hořava–Lifshitz-like $1+1$3-dimensional theory

$1+1$4

the gap equation displays a sharp parity-dependent dichotomy: for even $1+1$5, dimensional regularization yields only the trivial solution $1+1$6, while for odd $1+1$7 one obtains

$1+1$8

so a nonzero dynamical mass exists for attractive coupling $1+1$9 and parity is broken dynamically (Lima et al., 2016). At finite temperature, the same theory admits a critical temperature $N$0 at which the parity-breaking solution disappears; for $N$1 the paper gives an explicit closed form for $N$2 in terms of $N$3, $N$4, and $N$5 (Lima et al., 2016).

A distinct nonperturbative treatment based on the Cornwall–Jackiw–Tomboulis effective action finds three different possible mass-generating phases in the massless $N$6-dimensional model, depending on the cutoff dependence chosen for the bare coupling $N$7 (Khunjua et al., 2021). At first order in $N$8, the CJT stationarity equations admit: a Dirac mass $N$9 phase with $1+1$0, a Haldane mass $1+1$1 phase with $1+1$2, and a mixed $1+1$3 phase with

$1+1$4

leading to $1+1$5 for $1+1$6 (Khunjua et al., 2021). The paper interprets these as three different nontrivial phases distinguished by which discrete symmetries are broken and by which dynamical mass appears (Khunjua et al., 2021).

In $1+1$7 dimensions, finite-density large-$1+1$8 analysis has recently introduced a different meaning of “two-component”: only a fraction $1+1$9 of flavors couples to the chemical potential, while the remaining fraction is neutral (Thies, 19 Jun 2026). The grand potential in the homogeneous phase becomes

$N$0

and the full phase diagram in $N$1 exhibits a qualitative change at the critical filling fraction

$N$2

(Thies, 19 Jun 2026). For $N$3, the homogeneous theory has a tricritical point; for $N$4, that tricritical point is replaced by a critical end point inside the massive phase (Thies, 19 Jun 2026). Once inhomogeneous condensates are included, the phase diagram contains crystal phases bounded by second-order instability lines obtained from stability analyses rather than a full thermal Hartree–Fock solution (Thies, 19 Jun 2026).

The broadest recent extension of the mass-generation problem is the $N$5 tensor framework in $N$6 dimensions. Mean-field theory suggests a first transition at $N$7, where a singlet mass condenses and a discrete symmetry is broken, and a second transition at $N$8, where $N$9 (Han et al., 2024). Renormalization-group analysis of an $$\mathcal{L} = \sum_{i=1}^N \bar\psi_i (\gamma_\mu \partial_\mu + m)\psi_i - \frac{g^2}{2}\left(\sum_{i=1}^N \bar\psi_i\psi_i\right)^2,$$0-copy Majorana extension shows that while symmetric-tensor and adjoint-nematic fixed points exist at large $$\mathcal{L} = \sum_{i=1}^N \bar\psi_i (\gamma_\mu \partial_\mu + m)\psi_i - \frac{g^2}{2}\left(\sum_{i=1}^N \bar\psi_i\psi_i\right)^2,$$1, they lose criticality as $$\mathcal{L} = \sum_{i=1}^N \bar\psi_i (\gamma_\mu \partial_\mu + m)\psi_i - \frac{g^2}{2}\left(\sum_{i=1}^N \bar\psi_i\psi_i\right)^2,$$2, and at $$\mathcal{L} = \sum_{i=1}^N \bar\psi_i (\gamma_\mu \partial_\mu + m)\psi_i - \frac{g^2}{2}\left(\sum_{i=1}^N \bar\psi_i\psi_i\right)^2,$$3 become equivalent to the Gaussian fixed point, leaving only the standard Gross–Neveu singlet fixed point as genuinely critical (Han et al., 2024). A complementary Gross–Neveu–Yukawa analysis of the $$\mathcal{L} = \sum_{i=1}^N \bar\psi_i (\gamma_\mu \partial_\mu + m)\psi_i - \frac{g^2}{2}\left(\sum_{i=1}^N \bar\psi_i\psi_i\right)^2,$$4 symmetric-tensor order parameter finds that for $$\mathcal{L} = \sum_{i=1}^N \bar\psi_i (\gamma_\mu \partial_\mu + m)\psi_i - \frac{g^2}{2}\left(\sum_{i=1}^N \bar\psi_i\psi_i\right)^2,$$5, corresponding to the canonical model, the tensor transition is fluctuation-induced first order, and only for $$\mathcal{L} = \sum_{i=1}^N \bar\psi_i (\gamma_\mu \partial_\mu + m)\psi_i - \frac{g^2}{2}\left(\sum_{i=1}^N \bar\psi_i\psi_i\right)^2,$$6 does a stable continuous tensor fixed point appear (Han et al., 2024).

4. Bound states, solitons, and real-time dynamics

The model’s nonperturbative content is not limited to mass generation; it also includes bound states and solitons. In the $$\mathcal{L} = \sum_{i=1}^N \bar\psi_i (\gamma_\mu \partial_\mu + m)\psi_i - \frac{g^2}{2}\left(\sum_{i=1}^N \bar\psi_i\psi_i\right)^2,$$7-dimensional O$$\mathcal{L} = \sum_{i=1}^N \bar\psi_i (\gamma_\mu \partial_\mu + m)\psi_i - \frac{g^2}{2}\left(\sum_{i=1}^N \bar\psi_i\psi_i\right)^2,$$8 Gross–Neveu model with Wilson fermions, the fermion-loop formulation rewrites the partition function as a sum over closed, non-backtracking Majorana loops and monomers with positive weights (Maillart et al., 2011). In this representation, open worms sample fermion two-point functions, while composite worms sample bound-state correlators such as

$$\mathcal{L} = \sum_{i=1}^N \bar\psi_i (\gamma_\mu \partial_\mu + m)\psi_i - \frac{g^2}{2}\left(\sum_{i=1}^N \bar\psi_i\psi_i\right)^2,$$9

and the corresponding correlators show clean exponential decay, allowing extraction of both fermion and bosonic bound-state masses near the continuum limit (Maillart et al., 2011). Because the worm algorithm updates topological sectors and can work directly at the massless point, it provides determinant-free access to the spectrum with essentially no critical slowing down (Maillart et al., 2011).

A different route to bound-state physics appears in real-time simulations of the two-flavor massive model using matrix product states. On an 11-site lattice with 4 fermionic modes per site, the authors Trotterize the Hamiltonian and simulate up to four particles initially localized on one site (Li et al., 2020). They observe that when $2+1$00 is small and $2+1$01 is large, two- and four-fermion initial states remain localized and their constituent probability densities evolve “synchronistically,” which the paper interprets as dynamical evidence for bound-state behavior (Li et al., 2020). Flavor symmetry implies that particles of different flavor but the same spin evolve identically, and the Wilson term is required to suppress doubling artifacts (Li et al., 2020). A plausible implication is that the two-component model furnishes a useful benchmark for nonequilibrium algorithms because the same parameters that favor bound-state formation in continuum studies also leave a visible signature in real-space dynamics.

The richest exact soliton structure in the data arises in the $2+1$02-dimensional two-flavor O$2+1$03 chiral Gross–Neveu model,

$2+1$04

whose vacuum manifold consists of two disjoint circles (Thies, 2023). In large $2+1$05, the model supports a massless Goldstone mode and three massive mesons of mass $2+1$06, as well as explicit kink, baryon, and breather solutions (Thies, 2023). Via a Majorana decomposition, this model is dual to the “perfect Gross–Neveu model,” a one-flavor theory with both chiral and Cooper-pairing channels and full Pauli–Gürsey symmetry, itself dual to the Zakharov–Mikhailov model (Thies, 2023). This demonstrates that two-component Gross–Neveu variants can interpolate between conventional particle–hole condensation and particle–particle pairing while remaining integrable in large $2+1$07 (Thies, 2023).

5. Lattice formulations and numerical methods

Lattice formulations of the two-component Gross–Neveu model differ substantially in their fermion discretization and in the observables they make most accessible. Wilson fermions provide a direct path to loop/worldline algorithms. In the O$2+1$08 Majorana formulation, the lattice action

$2+1$09

with $2+1$10 and $2+1$11 yields positive loop weights and a worm algorithm that samples both fluctuating topological sectors and bound-state correlators (Maillart et al., 2011). Partition-function combinations such as $2+1$12 vanish at the critical point, allowing the critical Wilson mass to be located numerically by a zero crossing (Maillart et al., 2011).

Borici–Creutz fermions provide a minimally doubled, chirally invariant alternative. In $2+1$13 dimensions, the Gross–Neveu action with Borici–Creutz fermions is

$2+1$14

with $2+1$15 in 2D (Goswami et al., 2014). Strong-coupling analysis and hybrid Monte Carlo both find a second-order chiral phase transition, with a critical value $2+1$16 in the massless limit (Goswami et al., 2014). This formulation reproduces the qualitative chiral phase structure of the continuum model while introducing an additional $2+1$17 channel tied to the special Borici–Creutz direction (Goswami et al., 2014).

Staggered-fermion formulations make it possible to study competition between ordinary bilinear mass generation and more exotic symmetric mass generation. A $2+1$18-dimensional lattice model with two flavors of massless staggered fermions, a nearest-neighbor current–current interaction $2+1$19, and an on-site four-fermion interaction $2+1$20,

$2+1$21

is free of sign problems in the fermion bag representation (Maiti et al., 2021). Based on earlier limits, the authors expect three phases: a PMW phase with massless fermions, an FM phase where fermions become massive through spontaneous symmetry breaking and a bilinear condensate, and a PMS phase where fermions are massive without any bilinear condensate (Maiti et al., 2021). This suggests that even in a two-component Gross–Neveu–Thirring setting, the bilinear condensate is not the only possible mass-generating mechanism.

Tensor-network methods furnish a complementary Hamiltonian approach. The two-flavor lattice Hamiltonian used in real-time MPS simulation contains kinetic, Wilson, mass, and on-site interaction terms and acts on $2+1$22 fermionic modes, encoded as qubits via Jordan–Wigner (Li et al., 2020). For $2+1$23, $2+1$24, and $2+1$25, the full state would naively live in a $2+1$26-dimensional Hilbert space, but MPS compression keeps the simulation tractable, with reported MPS errors below $2+1$27 for the few-particle states considered (Li et al., 2020). The same work quotes a Trotter error of order $2+1$28, highlighting the balance between circuit fidelity and tensor-network compression in explicit time evolution (Li et al., 2020).

6. Modern extensions and related universality classes

Several recent developments situate the two-component Gross–Neveu model within a broader web of fixed points and universality classes. In $2+1$29 dimensions, the bosonized Gross–Neveu–Yukawa theory has been pushed to three loops in $2+1$30 dimensions for general flavor number $2+1$31, yielding explicit $2+1$32-expansions for $2+1$33, $2+1$34, and $2+1$35 at the chiral Ising fixed point (Mihaila et al., 2017). For $2+1$36, the paper quotes

$2+1$37

and for $2+1$38,

$2+1$39

(Mihaila et al., 2017). In the special $2+1$40 case, corresponding to a single-component fermion in the paper’s conventions, the exponents satisfy emergent super-scaling relations order by order up to three loops, consistent with $2+1$41 supersymmetry (Mihaila et al., 2017).

At the lower-critical end, the fate of the non-supersymmetric Gross–Neveu–Yukawa fixed point with a two-component Majorana fermion continued to two dimensions has been analyzed through fermionic minimal models (Nakayama et al., 2022). Under the assumptions that the fixed point is a fermionic minimal model with a chiral $2+1$42 symmetry and just two relevant singlet operators, only four candidates survive; matching topological defect line spin content under the assumed flow to the supersymmetric fermionic tricritical Ising model rules out two of them, leaving the fermionic $2+1$43 and fermionic $2+1$44 models (Nakayama et al., 2022). An additional double-braiding constraint favors the non-unitary fermionic $2+1$45 model (Nakayama et al., 2022). This suggests that the two-dimensional continuation of the non-supersymmetric two-component Majorana GNY fixed point may be non-unitary, even though the higher-dimensional parent theory is treated perturbatively.

The $2+1$46 generalization provides a final unifying viewpoint. Rewriting the canonical $2+1$47-dimensional Gross–Neveu interaction in terms of a single $2+1$48-component Majorana field exposes three quartic channels: a singlet scalar, a symmetric tensor, and an adjoint nematic (Han et al., 2024). Extending the theory to $2+1$49 copies reveals three corresponding fixed points at large $2+1$50, but only the singlet Gross–Neveu fixed point remains physically critical at $2+1$51 (Han et al., 2024). A distinct $2+1$52 analysis of an $2+1$53-symmetric, Fierz-complete Majorana theory reaches a closely related conclusion: the Gross–Neveu–Ising fixed point remains critical for all $2+1$54, whereas the symmetric-tensor fixed point loses criticality below

$2+1$55

so for the physical $2+1$56 case only the singlet Gross–Neveu–Ising transition survives as a genuine continuous transition (Hawashin et al., 27 Oct 2025). This reinforces a common modern reading: the canonical two-component Gross–Neveu model is robustly associated with singlet mass generation, while more elaborate tensor or nematic channels are generically fluctuation-destabilized unless the fermion flavor content is enlarged (Han et al., 2024, Han et al., 2024, Hawashin et al., 27 Oct 2025).