Quartic Yukawa Model Overview

• The quartic Yukawa model is defined by coupling scalar and fermionic fields through Yukawa and quartic self-interactions, forming a cornerstone of quantum and statistical field theories.
• Renormalization group analyses reveal a Gaussian fixed point and rich flow dynamics that bridge perturbative and nonperturbative regimes.
• Extensions including tensor, supersymmetric, and lattice variants enable studies of conformal dynamics, operator mixing, and cosmological implications.

The quartic Yukawa model denotes a class of quantum field theories in which scalar and (often) fermionic matter fields interact via both a standard Yukawa term and a quartic (or higher-degree) scalar self-interaction. In the broadest sense, this encompasses renormalizable models with Lagrangians of the form $\mathcal{L} \sim -g\,\phi\,\bar\psi\psi$ and $-(\lambda/4!)\,\phi^4$, as well as their extensions to multiple fields, higher symmetries, nontrivial topology (e.g., tensor structure), and inclusion of other couplings or sectors. Such models form a backbone of particle theory, statistical field theory, and the study of strongly coupled conformal behavior, with central importance in both perturbative and nonperturbative regimes.

1. Model Definition and Core Lagrangian Structures

In canonical form, the quartic Yukawa model consists of a real scalar $\phi(x)$ and a Dirac fermion $\psi(x)$ described by the Lagrangian

$$\mathcal{L} = \tfrac{1}{2}\,\partial_\mu\phi\,\partial^\mu\phi - \tfrac{1}{2}\,m_\phi^2\,\phi^2 - \tfrac{\lambda}{4!}\,\phi^4 + i\,\bar\psi\gamma^\mu\partial_\mu\psi - m_\psi\,\bar\psi\psi - g\,\phi\,\bar\psi\psi,$$

where $\lambda$ is the quartic scalar self-coupling and $g$ the Yukawa coupling (Toms, 2018). Generalizations include the addition of further quartic or higher monomials, multi-field extensions, and non-minimal kinetic or curvature couplings. In four-dimensional flat spacetime, both $g$ and $\lambda$ are dimensionless and the theory is renormalizable.

Specialized extensions involve:

• Tensor models with $O(N)^3$ symmetry, coupling rank-3 real scalar and Dirac fields via $-(\lambda/4!)\,\phi^4$0-invariant operators such as $-(\lambda/4!)\,\phi^4$1, and including sextic interactions (e.g., $-(\lambda/4!)\,\phi^4$2) (Fraser-Taliente et al., 2024).
• Higgs-like mediator models with Lagrangians featuring additional cubic and quartic self-interactions for the mediator field, such as $-(\lambda/4!)\,\phi^4$3 (Chigodaev et al., 2011).
• Scalar-doublet extensions, e.g., real singlet $-(\lambda/4!)\,\phi^4$4 and complex $-(\lambda/4!)\,\phi^4$5 doublet $-(\lambda/4!)\,\phi^4$6 fields with all quartic and Yukawa couplings (Saad et al., 29 Jun 2025).
• Supersymmetric variants, such as $-(\lambda/4!)\,\phi^4$7 SQCD, with both gluino-quark-squark Yukawa interactions and four-squark quartic terms (Costa et al., 2022).

2. Renormalization, β-Functions, and Fixed Points

The renormalization group (RG) behavior of the quartic Yukawa model is controlled by the coupled beta functions for the Yukawa and quartic couplings: $-(\lambda/4!)\,\phi^4$8 with, at one loop in the minimal model (Toms, 2018, Davies et al., 2021): $-(\lambda/4!)\,\phi^4$9

$\phi(x)$0

Higher-loop results have been computed in full generality up to three loops for Yukawa and two loops for quartic couplings (Davies et al., 2021): $\phi(x)$1 The only perturbative fixed point is the Gaussian $\phi(x)$2 in four dimensions. In large-$\phi(x)$3 tensorial versions, the RG flow is much richer, exhibiting multiple interacting fixed points: melonic (fermion sector), prismatic (bosonic higher-order), and mixed fixed points, connected via RG flows in $\phi(x)$4 (Fraser-Taliente et al., 2024). Supersymmetric points also emerge under appropriate dof balancing.

3. Melonic Tensor Extensions and Conformal Dynamics

Tensorial quartic Yukawa models introduce real rank-$\phi(x)$5 tensors (typically $\phi(x)$6) for both $\phi(x)$7 and $\phi(x)$8, with $\phi(x)$9 invariance, and admit large-$\psi(x)$0 limits dominated by “melonic” diagrams: $\psi(x)$1 In the scaling $\psi(x)$2, $\psi(x)$3, the Schwinger–Dyson equations for $\psi(x)$4 and $\psi(x)$5 admit conformal solutions in the deep IR,

$\psi(x)$6

fixing anomalous dimensions and enabling a full resummation of melonic diagrams (Fraser-Taliente et al., 2024). Nonperturbative analysis recovers the perturbative $\psi(x)$7-expansions for scaling dimensions, and establishes the existence of IR “windows of stability” for real scaling dimensions in non-integer $\psi(x)$8.

A conformal operator spectrum with Regge trajectories emerges: bilinear singlets have analytic scaling dimensions with shadow symmetry $\psi(x)$9. Characteristic features are unitarity wedges, supersymmetric special points, and an exact mapping between perturbative and melonic SDE results.

4. Nonlinear Interactions, Few-Body Bound States, and Potentials

The quartic Yukawa model admits extensions in which the scalar mediator is treated as an independent field with self-interactions, yielding nontrivial few-body dynamics. For instance, with complex scalar matter $$\mathcal{L} = \tfrac{1}{2}\,\partial_\mu\phi\,\partial^\mu\phi - \tfrac{1}{2}\,m_\phi^2\,\phi^2 - \tfrac{\lambda}{4!}\,\phi^4 + i\,\bar\psi\gamma^\mu\partial_\mu\psi - m_\psi\,\bar\psi\psi - g\,\phi\,\bar\psi\psi,$$0 and massive mediator $$\mathcal{L} = \tfrac{1}{2}\,\partial_\mu\phi\,\partial^\mu\phi - \tfrac{1}{2}\,m_\phi^2\,\phi^2 - \tfrac{\lambda}{4!}\,\phi^4 + i\,\bar\psi\gamma^\mu\partial_\mu\psi - m_\psi\,\bar\psi\psi - g\,\phi\,\bar\psi\psi,$$1: $$\mathcal{L} = \tfrac{1}{2}\,\partial_\mu\phi\,\partial^\mu\phi - \tfrac{1}{2}\,m_\phi^2\,\phi^2 - \tfrac{\lambda}{4!}\,\phi^4 + i\,\bar\psi\gamma^\mu\partial_\mu\psi - m_\psi\,\bar\psi\psi - g\,\phi\,\bar\psi\psi,$$2 (Chigodaev et al., 2011). Integrating out $$\mathcal{L} = \tfrac{1}{2}\,\partial_\mu\phi\,\partial^\mu\phi - \tfrac{1}{2}\,m_\phi^2\,\phi^2 - \tfrac{\lambda}{4!}\,\phi^4 + i\,\bar\psi\gamma^\mu\partial_\mu\psi - m_\psi\,\bar\psi\psi - g\,\phi\,\bar\psi\psi,$$3 produces nonlocal potentials:

• Two-body: Yukawa $$\mathcal{L} = \tfrac{1}{2}\,\partial_\mu\phi\,\partial^\mu\phi - \tfrac{1}{2}\,m_\phi^2\,\phi^2 - \tfrac{\lambda}{4!}\,\phi^4 + i\,\bar\psi\gamma^\mu\partial_\mu\psi - m_\psi\,\bar\psi\psi - g\,\phi\,\bar\psi\psi,$$4, $$\mathcal{L} = \tfrac{1}{2}\,\partial_\mu\phi\,\partial^\mu\phi - \tfrac{1}{2}\,m_\phi^2\,\phi^2 - \tfrac{\lambda}{4!}\,\phi^4 + i\,\bar\psi\gamma^\mu\partial_\mu\psi - m_\psi\,\bar\psi\psi - g\,\phi\,\bar\psi\psi,$$5;
• Three- and four-body: genuine multipoint corrections from $$\mathcal{L} = \tfrac{1}{2}\,\partial_\mu\phi\,\partial^\mu\phi - \tfrac{1}{2}\,m_\phi^2\,\phi^2 - \tfrac{\lambda}{4!}\,\phi^4 + i\,\bar\psi\gamma^\mu\partial_\mu\psi - m_\psi\,\bar\psi\psi - g\,\phi\,\bar\psi\psi,$$6 and $$\mathcal{L} = \tfrac{1}{2}\,\partial_\mu\phi\,\partial^\mu\phi - \tfrac{1}{2}\,m_\phi^2\,\phi^2 - \tfrac{\lambda}{4!}\,\phi^4 + i\,\bar\psi\gamma^\mu\partial_\mu\psi - m_\psi\,\bar\psi\psi - g\,\phi\,\bar\psi\psi,$$7 terms, with analytic representations for special kinematics, always preserving net attraction at large $$\mathcal{L} = \tfrac{1}{2}\,\partial_\mu\phi\,\partial^\mu\phi - \tfrac{1}{2}\,m_\phi^2\,\phi^2 - \tfrac{\lambda}{4!}\,\phi^4 + i\,\bar\psi\gamma^\mu\partial_\mu\psi - m_\psi\,\bar\psi\psi - g\,\phi\,\bar\psi\psi,$$8.

These corrections enhance binding at short distances, potentially develop local wells for clusters, and never induce repulsive long-range potentials when $$\mathcal{L} = \tfrac{1}{2}\,\partial_\mu\phi\,\partial^\mu\phi - \tfrac{1}{2}\,m_\phi^2\,\phi^2 - \tfrac{\lambda}{4!}\,\phi^4 + i\,\bar\psi\gamma^\mu\partial_\mu\psi - m_\psi\,\bar\psi\psi - g\,\phi\,\bar\psi\psi,$$9. Variational Fock-space analyses yield relativistic Bethe–Salpeter-like equations for few-body bound state wavefunctions.

5. Quantum Corrections in Cosmology and Critical Phenomena

Quartic Yukawa models have direct implications for the inflationary power spectrum and reheating dynamics: $\lambda$0 for a massless, minimally coupled inflaton with quartic and Yukawa interactions (Onemli, 2015).

One-loop corrections induce opposite spectral tilts:

• Yukawa loop: blue-tilted corrections, $\lambda$1.
• Quartic self-interaction: red-tilted corrections, $\lambda$2.

In reheating with a quartic inflaton potential $\lambda$3 and Yukawa decay to Dirac fermions, nonperturbative effects such as parametric resonance of $\lambda$4, kinematic blocking, and Pauli suppression dominate for $\lambda$5, substantially suppressing the reheating efficiency unless $\lambda$6 (Bhusal et al., 18 Dec 2025).

6. Lattice Formulations, Operator Mixing, and Triviality

The quartic Yukawa model on the lattice reveals critical phenomena and operator mixing structures. In lattice studies with a real scalar singlet and an $\lambda$7 doublet, the action incorporates all quartic vertices and a Yukawa-type cubic (Saad et al., 29 Jun 2025). Key features include:

• Non-trivial scalar propagator dressing functions, strongly enhanced at low $\lambda$8 (“non-triviality”).
• Yukawa vertex exhibits weak momentum dependence, no Landau pole.
• Quartic self-coupling effects are mitigated by the presence of other cubic and quartic terms.
• The spectrum of $\lambda$9 states uncovers an ultra-light band with an observed scarcity around hundreds of GeV, reflecting significant operator mixing.
• No evidence for first- or second-order phase transitions within the explored parameter space; the model remains in a single analytically connected phase.

Supersymmetric lattice formulations require precise fine-tuning of the Yukawa and quartic couplings to restore continuum SUSY Ward identities. Operator mixing under renormalization is controlled via discrete symmetries and explicit one-loop calculation of $g$0-factors and mixing matrices (Costa et al., 2022).

7. Outlook: Universal Properties and Research Frontiers

The quartic Yukawa model, both in minimal and extended forms, provides a laboratory for the analytic study of strongly coupled field theories:

• Large-$g$1 melonic CFTs universally exhibit dominance of solvable melon diagrams, nontrivial infrared conformal spectra, IR/UV stability “wedges,” and analytic matching between perturbative ε-expansion and nonperturbative SDE machinery (Fraser-Taliente et al., 2024).
• Triviality in the pure scalar sector is typically resolved in coupled models due to enhanced interaction structure.
• Multiloop RG functions now available for wide classes of models, facilitating controlled studies of fixed points, stability, and operator mixing (Davies et al., 2021).
• Cosmological and particle-physics applications depend sensitively on the interplay between quartic and Yukawa couplings, particularly in nonequilibrium dynamics.
• Supersymmetric cases require and enable nonperturbative restoration of full SUSY at the lattice level, with direct implications for lattice gauge theory and beyond (Costa et al., 2022).

Ongoing directions include further exploration of operator mixing, nontrivial infrared dynamics in large-$g$2 limits, and the phenomenology of extended and supersymmetric quartic Yukawa systems.