pith. machine review for the scientific record. sign in

arxiv: 2512.24215 · v2 · submitted 2025-12-30 · ✦ hep-ph

Recognition: no theorem link

Neutrino Mass, Vacuum Stability and Higgs Inflation with Vector-Like Quarks and a Single Right-Handed Neutrino

Authors on Pith no claims yet

Pith reviewed 2026-05-16 18:55 UTC · model grok-4.3

classification ✦ hep-ph
keywords neutrino massvacuum stabilityHiggs inflationvector-like quarksright-handed neutrinorenormalization group evolutionType-I seesawelectroweak vacuum
0
0 comments X

The pith

A Standard Model extension with n degenerate vector-like quarks and one right-handed neutrino stabilizes the electroweak vacuum to the Planck scale while generating viable neutrino masses and matching Higgs inflation observables to data.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper investigates an extension of the Standard Model that adds n degenerate down-type vector-like quarks together with a single right-handed neutrino. These additions modify the renormalization group running of the Higgs quartic coupling so that it stays positive up to the Planck scale. The right-handed neutrino simultaneously accounts for the observed light neutrino masses through the Type-I seesaw and improves the high-scale behavior of the potential. Using two-loop Standard Model running plus one-loop new particle effects with threshold matching, the authors locate parameter regions that achieve vacuum stability, produce acceptable neutrino masses, and yield inflationary predictions for the spectral index and tensor ratio that agree with the latest Planck and ACT data.

Core claim

The SM+(n)VLQ+RHN model produces RG trajectories that keep the Higgs potential stable to the Planck scale, generate phenomenologically viable neutrino masses via Type-I seesaw, and deliver n_s and r values consistent with Planck-LB-BK18 and ACT-LB-BK18 observations in the metric formulation of non-minimal Higgs inflation.

What carries the argument

Renormalization group evolution of the Higgs quartic coupling including one-loop contributions from degenerate vector-like quarks and the right-handed neutrino, with threshold matching, to obtain the RG-improved potential for non-minimal inflation.

If this is right

  • Specific regions of the parameter space (n, y_D, M_D) ensure the Higgs quartic remains positive to the Planck scale.
  • The Type-I seesaw mechanism from the single right-handed neutrino produces light neutrino masses consistent with experiment.
  • The computed inflationary observables n_s and r lie within the allowed contours of recent CMB data.
  • The VLQ contributions provide the primary stabilization effect, with the RHN playing a supporting role at high scales.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Vector-like quarks in the stabilizing mass range could be searched for at high-energy colliders.
  • Removing the degeneracy assumption on the VLQs could enlarge the viable parameter space.
  • The same RG stabilization technique might address related issues such as coupling unification at high scales.
  • Future precision measurements of the Higgs self-coupling could test the required threshold corrections.

Load-bearing premise

The vector-like quarks must be degenerate with identical masses and Yukawa couplings, and the one-loop threshold matching plus the specific non-minimal coupling must suffice for accurate high-scale running.

What would settle it

Detection of a Higgs quartic coupling that turns negative at scales below the Planck scale in improved calculations, or inflationary observables outside the model's predicted range, or no evidence for vector-like quarks in the relevant mass window.

Figures

Figures reproduced from arXiv: 2512.24215 by Canan Karahan.

Figure 1
Figure 1. Figure 1: FIG. 1. Vacuum stability in the SM+( [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Running of the Higgs quartic coupling [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Vacuum stability in the SM+( [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Running of the Higgs quartic coupling [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Predictions for the spectral index [PITH_FULL_IMAGE:figures/full_fig_p013_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Predictions for the spectral index [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
read the original abstract

We investigate a Standard Model extension containing $n$ degenerate down-type isosinglet vector-like quarks (VLQs) with masses $M_{\mathcal D}$ and Yukawa couplings $y_{\mathcal D}$, supplemented by a single right-handed neutrino (RHN), aiming to simultaneously address neutrino mass generation, electroweak vacuum stability, and Higgs inflation. The VLQs play the dominant role in stabilizing the Higgs potential through their impact on the renormalization-group evolution, while the RHN generates light neutrino masses via a Type-I seesaw mechanism and smooths the high-scale running of the Higgs quartic coupling in the inflationary regime. We perform a two-loop Standard Model renormalization-group equation analysis supplemented by the one-loop contributions of the VLQs and the RHN, with proper matching across their mass thresholds. Using these RG trajectories, we identify the regions in $(n,\, y_{\mathcal D},\, M_{\mathcal D})$ that stabilize the Higgs potential up to the Planck scale while satisfying experimental constraints. Employing the RG-improved Higgs potential in the metric formulation of non-minimal Higgs inflation, we compute the inflationary observables $n_s$ and $r$. The SM+$(n)$VLQ+RHN framework yields predictions consistent with the latest Planck-LB-BK18 and ACT-LB-BK18 data, while simultaneously ensuring electroweak vacuum stability and phenomenologically viable neutrino masses within well-defined regions of parameter space. For comparison, we also investigate the SM+$(n)$VLQ limit and present its vacuum stability and Higgs inflation predictions as a reference to quantify the stabilizing role of the VLQ sector.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript investigates a Standard Model extension with n degenerate down-type vector-like quarks (VLQs) of mass M_D and Yukawa y_D, plus a single right-handed neutrino (RHN). It performs a two-loop SM RGE analysis supplemented by one-loop VLQ and RHN contributions with threshold matching to identify parameter regions (n, y_D, M_D) that stabilize the Higgs potential to the Planck scale, generate viable neutrino masses via Type-I seesaw, and yield n_s and r from the RG-improved non-minimal Higgs potential (metric formulation) consistent with Planck-LB-BK18 and ACT-LB-BK18 data. A comparison to the SM+nVLQ case is included to isolate the VLQ stabilizing role.

Significance. If the results hold, this provides a minimal unified framework simultaneously addressing neutrino mass, electroweak vacuum stability, and Higgs inflation, with explicit viable parameter regions and a direct comparison quantifying each sector's contribution. The RG-improved potential approach for inflation and the use of latest cosmological data contours are standard strengths that make the predictions falsifiable.

major comments (2)
  1. [RG evolution section (beta functions for λ)] The renormalization-group analysis (abstract and RG evolution section) employs two-loop SM beta functions for λ but only one-loop contributions from the VLQs and RHN, with one-loop threshold matching at M_D and the RHN mass. Near the would-be instability scale, two-loop corrections from the new fields can be comparable in size to the SM two-loop pieces and can alter the sign of β_λ or shift the location of the minimum; this directly impacts the central claim of stability up to the Planck scale and the subsequent n_s, r predictions.
  2. [Parameter scan and degeneracy assumption] The VLQs are assumed degenerate with identical masses and Yukawas (abstract and parameter scan). This collapses the beta-function coefficients; any realistic mass or Yukawa splitting would modify the running of λ and the high-scale behavior, potentially moving the stability regions and the (n_s, r) points relative to the Planck-LB-BK18 and ACT-LB-BK18 contours.
minor comments (2)
  1. [Abstract and results section] The abstract states that the RHN 'smooths the high-scale running' but the explicit numerical impact on β_λ from the RHN one-loop term versus the VLQ terms is not quantified in a dedicated plot or table.
  2. [Inflation section] Notation for the non-minimal coupling ξ and the precise matching procedure across thresholds could be clarified with an explicit equation for the threshold correction to λ.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment point by point below, indicating planned revisions where appropriate.

read point-by-point responses
  1. Referee: The renormalization-group analysis (abstract and RG evolution section) employs two-loop SM beta functions for λ but only one-loop contributions from the VLQs and RHN, with one-loop threshold matching at M_D and the RHN mass. Near the would-be instability scale, two-loop corrections from the new fields can be comparable in size to the SM two-loop pieces and can alter the sign of β_λ or shift the location of the minimum; this directly impacts the central claim of stability up to the Planck scale and the subsequent n_s, r predictions.

    Authors: We agree that a complete two-loop treatment of the BSM contributions would be desirable for maximal precision near the would-be instability scale. Our analysis follows the common approach in the literature of retaining two-loop SM beta functions while adding one-loop BSM terms, which capture the dominant positive contribution to β_λ from the VLQ Yukawa y_D. For the parameter regions we identify, the stabilization occurs sufficiently below the Planck scale that the leading one-loop BSM effects control the sign of β_λ. In the revised manuscript we will add a dedicated paragraph in the RG section estimating the expected size of the omitted two-loop BSM corrections (based on the magnitude of the one-loop terms) and explicitly noting that a full two-loop BSM calculation is left for future work. revision: partial

  2. Referee: The VLQs are assumed degenerate with identical masses and Yukawas (abstract and parameter scan). This collapses the beta-function coefficients; any realistic mass or Yukawa splitting would modify the running of λ and the high-scale behavior, potentially moving the stability regions and the (n_s, r) points relative to the Planck-LB-BK18 and ACT-LB-BK18 contours.

    Authors: The degeneracy assumption is introduced to keep the parameter space minimal while isolating the stabilizing role of the VLQ sector. For small relative splittings the beta-function coefficients receive averaged contributions, so the qualitative location of the stability bands and the resulting (n_s, r) predictions are expected to remain similar. We will insert a short clarifying remark in the parameter-scan section stating that the degenerate case is adopted as a benchmark and that non-degenerate spectra would require a broader numerical scan, which we leave for future study. revision: partial

Circularity Check

1 steps flagged

Parameter scan over n, y_D, M_D selects regions that enforce stability and data consistency, rendering ns/r agreement fitted rather than predicted

specific steps
  1. fitted input called prediction [Abstract (final paragraph) and sections on RG trajectories + inflationary observables]
    "The SM+(n)VLQ+RHN framework yields predictions consistent with the latest Planck-LB-BK18 and ACT-LB-BK18 data, while simultaneously ensuring electroweak vacuum stability and phenomenologically viable neutrino masses within well-defined regions of parameter space."

    The 'well-defined regions' are obtained by scanning n, y_D, M_D until both vacuum stability and the ns/r values fall inside the experimental contours; the quoted consistency is therefore the result of that selection, not an a-priori prediction of the RGEs.

full rationale

The paper scans the discrete parameter space (n, y_D, M_D) to locate intervals that keep the Higgs quartic positive to the Planck scale and produce ns, r inside the cited Planck/ACT contours after RG improvement. Once those intervals are chosen, the inflationary observables are read off from the same RG trajectories; the agreement is therefore a direct consequence of the selection criterion rather than an independent output of the equations. The underlying two-loop SM + one-loop new-physics RGEs themselves contain no circularity, but the load-bearing claim that the framework 'yields predictions consistent with data' reduces to the fitting step. No self-citation chains, self-definitional equations, or smuggled ansatze are present.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 2 invented entities

The central claim rests on standard QFT renormalization-group methods plus two new particle sectors whose effects are added at one-loop; no machine-checked proofs or external data anchors are supplied.

free parameters (3)
  • n
    Integer number of degenerate VLQs chosen to achieve stability up to Planck scale.
  • y_D
    Yukawa coupling of the VLQs scanned over ranges that keep the quartic positive.
  • M_D
    Common mass of the VLQs used for threshold matching in the RG flow.
axioms (2)
  • standard math Two-loop Standard Model RGEs supplemented by one-loop VLQ and RHN contributions with proper threshold matching
    Invoked in the abstract as the computational backbone.
  • domain assumption Type-I seesaw mechanism generates light neutrino masses from the single RHN
    Standard assumption for neutrino mass generation in the model.
invented entities (2)
  • n degenerate down-type vector-like quarks no independent evidence
    purpose: Dominant stabilization of the Higgs quartic via RG evolution
    New particles introduced to alter the running of the Higgs self-coupling.
  • single right-handed neutrino no independent evidence
    purpose: Generates neutrino masses and smooths high-scale Higgs running
    Postulated to implement Type-I seesaw and assist inflation regime.

pith-pipeline@v0.9.0 · 5598 in / 1584 out tokens · 69405 ms · 2026-05-16T18:55:42.378813+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

64 extracted references · 64 canonical work pages · 23 internal anchors

  1. [1]

    computations of the Higgs effective potential. In literature, numerous extensions of the Standard Model have been proposed to address this issue, such as introducing additional scalar fields to raise the Higgs quartic coupling at high energies [17–20], or new fermions/gauge sectors that alter the renormalization-group (RG) running ofλ H [21–23]. In partic...

  2. [2]

    For anSU(2) L singlet with hyperchargeY=− 1 3, the covariant derivative isD µ =∂ µ −ig 3T aGa µ −ig ′Y Bµ

    Down-Type Isosinglet Vector-Like Quark Sector The kinetic and mass terms of down-type isosinglet VLQ sector are LD = Di /DD −M D DD,(3) whereM D is a gauge-invariant Dirac mass. For anSU(2) L singlet with hyperchargeY=− 1 3, the covariant derivative isD µ =∂ µ −ig 3T aGa µ −ig ′Y Bµ. Gauge invariance permits for a single renormalizable Yukawa operator: LY...

  3. [3]

    Here, ˜H≡iσ 2H ∗ denotes theSU(2) L- conjugate Higgs doublet

    Right-Handed Neutrino Sector The RHN kinetic and Majorana mass terms are LN = 1 2 N i /∂ N− 1 2 MN N cN.(6) The neutrino Yukawa interaction is LYuk N =−y N L eH N+ h.c.,(7) producing a Dirac massm D =y N v/ √ 2 after electroweak symmetry breaking. Here, ˜H≡iσ 2H ∗ denotes theSU(2) L- conjugate Higgs doublet. A light neutrino mass arises from the Type-I se...

  4. [4]

    This relation follows from the standard diagonalization of the VLQ–SM mass matrix and is consistent with the formalism developed in vector-like quark studies [36, 37]

    Experimental constraints on the Yukawa couplingy D and massM D In the SM extended by down-type isosinglet VLQs with the massM D and Yukawa interactiony D, the mass terms after electroweak symmetry breaking lead to a left-handed mixing angle defined by sinθ L ≃ yD v√ 2M D ,(11) wherev≃246 GeV is the Higgs vacuum expectation value. This relation follows fro...

  5. [5]

    These constraints ensure that the theory remains perturbative, unitary, and free of Landau poles up to the Planck scale

    Theoretical constraints on the Yukawa couplingy D In addition to the experimental limits discussed above, the down-type VLQ Yukawa couplingy D is subject to several theoretical consistency conditions. These constraints ensure that the theory remains perturbative, unitary, and free of Landau poles up to the Planck scale. a. Perturbativity and Landau–pole c...

  6. [6]

    The point at which a given trajectory intersects the lineλ min = 0 signals the onset of vacuum instability for the corresponding (n, M D)

    For fixed values of (n, M D), the minimum value of the Higgs quartic coupling,λ min, exhibits a monotonic decrease as the VLQ Yukawa couplingy D is increased. The point at which a given trajectory intersects the lineλ min = 0 signals the onset of vacuum instability for the corresponding (n, M D). Parameter regions with λmin >0 correspond to an absolutely ...

  7. [7]

    Comparing the three panels, the range of Yukawa couplingsy D that keepsλ min >0 becomes more restrictive as VLQ massM D increase

  8. [8]

    For sufficiently small Yukawa couplings, increasingnshiftsλ min upward and improves vacuum stability

    At all three benchmark masses, the dependence on the number of VLQsnshows a characteristic two–regime behaviour. For sufficiently small Yukawa couplings, increasingnshiftsλ min upward and improves vacuum stability. However, beyond a certainy D value in each panel, the destabilizing∼ −n y 4 D contribution dominates the running, and the trend reverses: larg...

  9. [9]

    Forn≥3, vacuum stability can be achieved within a finite interval of smally D, but this interval becomes progressively narrower asM D is increased

    In particular, the casen= 1 andn= 2 never yields a stable electroweak vacuum for any of the benchmark 8 masses. Forn≥3, vacuum stability can be achieved within a finite interval of smally D, but this interval becomes progressively narrower asM D is increased. Taken together, these results show that, in the presence of a single RHN withy N = 0.42, restorin...

  10. [10]

    Consequently, for eachn,λ min stays higher than in the full SM+(n)VLQ+RHN model

    Removing the RHN Yukawa eliminates the negative contribution that previously pushedλ min downward at high scales. Consequently, for eachn,λ min stays higher than in the full SM+(n)VLQ+RHN model. This confirms that VLQs alone have a net stabilizing tendency at small Yukawa couplings

  11. [11]

    However, without the RHN, this stabilizing effect is stronger and persists to larger values ofy D

    As in the full model, increasing the number of VLQs shifts the running upward, delaying (or preventing) the 10 point at whichλ min turns negative. However, without the RHN, this stabilizing effect is stronger and persists to larger values ofy D

  12. [12]

    In other words, the destabilizing regime sets in more softly because only the VLQ sector drives it, not the RHN

    The characteristic transition seen previously — where large values ofneventually accelerate the fall ofλ min still appears, but with a reduced rate and at later scales. In other words, the destabilizing regime sets in more softly because only the VLQ sector drives it, not the RHN

  13. [13]

    Here, for SM+(n)VLQ, even smallnvalues lead to a visibly less negativeλ min, and several trajectories remain close to or above zero

    In the SM+(n)VLQ+RHN case,n= 1 andn= 2 never achieved stability in the allowed parameter space. Here, for SM+(n)VLQ, even smallnvalues lead to a visibly less negativeλ min, and several trajectories remain close to or above zero. Forn= 2, a vacuum–stable region exists only for the lightest VLQ massM D = 1.5 TeV, while forM D = 3.0 and 5.0 TeV then= 2 curve...

  14. [14]

    In this case, vacuum stability is still governed primarily by the interplay betweennandy D

    Overall, the SM+(n)VLQ scenario leads to improved electroweak vacuum stability compared to the SM+(n)VLQ +RHN model for all considered values ofnandM D, as reflected by larger values ofλ min and a broader region withλ min >0. In this case, vacuum stability is still governed primarily by the interplay betweennandy D. The absence of the RHN Yukawa contribut...

  15. [15]

    Standard Model (1-loop) At one loop the SMβ-functions for the gauge couplings, top Yukawa and Higgs quartic coupling read [42] βSM,1L g1 = 41 10 g3 1,(A2) βSM,1L g2 =− 19 6 g3 2,(A3) βSM,1L g3 =−7g 3 3,(A4) βSM,1L yt =y t 9 2 y2 t −8g 2 3 − 9 4 g2 2 − 17 20 g2 1 ,(A5) βSM,1L λ = 24λ2 −6y 4 t + 9 8 g4 2 + 27 200 g4 1 + 9 20 g2 2g2 1 +λ 12y2 t −9g 2 2 − 9 5...

  16. [16]

    Standard Model (2-loop top-sector terms) We also include the dominant two-loop terms involving the top Yukawa, matching exactly the structure used in the numerical implementation [43]: βSM,2L g1 =g 3 1 199 50 g2 1 + 27 10 g2 2 + 44 5 g2 3 − 17 10 y2 t ,(A7) βSM,2L g2 =g 3 2 9 10 g2 1 + 35 6 g2 2 + 12g2 3 − 3 2 y2 t ,(A8) βSM,2L g3 =g 3 3 11 10 g2 1 + 9 2 ...

  17. [17]

    Down-type singlet VLQ contributions (1 loop) In the following, we present the RGE fory D together with the corresponding one-loop contributions relevant for the present analysis [24, 25]. βVLQ yD =y D 3 2 y2 t + 3 2 + 3n y2 D −8g 2 3 − 9 4 g2 2 − 1 4 g2 1 ,(A12) ∆βVLQ g1 = 4 15 ng3 1,(A13) ∆βVLQ g2 = 0,(A14) ∆βVLQ g3 = 2 3 ng3 3,(A15) ∆βVLQ yt = 3 2 n y2 ...

  18. [18]

    [44, 45] βRHN yN =y N 5 2 y2 N + 3n y2 D + 3y2 t − 9 20 g2 1 − 9 4 g2 2 ,(A18) ∆βRHN λ =−2y 4 N + 4y2 N λ.(A19) 17

    Majorana RHN contributions (1 loop) In the following, we present the RGE for the RHN Yukawa coupling and retain only its dominant one-loop con- tributions to the Higgs quartic coupling, which are sufficient to capture the impact of the RHN sector on vacuum stability. [44, 45] βRHN yN =y N 5 2 y2 N + 3n y2 D + 3y2 t − 9 20 g2 1 − 9 4 g2 2 ,(A18) ∆βRHN λ =−...

  19. [19]

    Full RG system in theSM+ (n)V LQ+RHNFramework In this part, we summarize the RGEs for the couplings in theSM+ (n)V LQ+RHNframework. Collecting all contributions, the fullβ-functions used in the numerical analysis read βg1 = 1 16π2 βSM,1L g1 + 1 (16π2)2 βSM,2L g1 + 1 16π2 ∆βVLQ g1 ,(A20) βg2 = 1 16π2 βSM,1L g2 + 1 (16π2)2 βSM,2L g2 ,(A21) βg3 = 1 16π2 βSM,...

  20. [20]

    Measurement of the rate of nu_e + d --> p + p + e^- interactions produced by 8B solar neutrinos at the Sudbury Neutrino Observatory

    Q. R. Ahmadet al.[SNO Collaboration], “Measurement of the Rate ofν e +d→p+p+e − Interactions Produced by 8B Solar Neutrinos at the Sudbury Neutrino Observatory,” Phys. Rev. Lett.87, 071301 (2001), [arXiv:nucl-ex/0106015]

  21. [21]

    Evidence for oscillation of atmospheric neutrinos

    Y. Fukudaet al.[Super-Kamiokande Collaboration], “Evidence for Oscillation of Atmospheric Neutrinos,” Phys. Rev. Lett. 81, 1562–1567 (1998), [arXiv:hep-ex/9807003]. 18

  22. [22]

    Higgs mass and vacuum stability in the Standard Model at NNLO

    G. Degrassi et al., “Higgs mass and vacuum stability in the Standard Model at NNLO,” JHEP08, 098 (2012), [arXiv:1205.6497 [hep-ph]]

  23. [23]

    Investigating the near-criticality of the Higgs boson

    D. Buttazzo et al., “Investigating the near-criticality of the Higgs boson,” JHEP12, 089 (2013), [arXiv:1307.3536 [hep-ph]]

  24. [24]

    A New Type of Isotropic Cosmological Models Without Singularity,

    A. A. Starobinsky, “A New Type of Isotropic Cosmological Models Without Singularity,” Phys. Lett. B91, 99 (1980). doi:10.1016/0370-2693(80)90670-X

  25. [25]

    The Standard Model Higgs boson as the inflaton

    F. L. Bezrukov and M. Shaposhnikov, “The Standard Model Higgs boson as the inflaton,” Phys. Lett. B659, 703 (2008), [arXiv:0710.3755 [hep-th]]

  26. [27]

    Complex Spinors and Unified Theories

    M. Gell-Mann, P. Ramond, and R. Slansky, “Complex Spinors and Unified Theories,” inSupergravity, edited by P. van Nieuwenhuizen and D. Z. Freedman (North Holland, Amsterdam, 1979), p. 315 [arXiv:1306.4669 [hep-th]]

  27. [28]

    Yanagida, inProceedings of the Workshop on Unified Theory and the Baryon Number of the Universe, edited by O

    T. Yanagida, inProceedings of the Workshop on Unified Theory and the Baryon Number of the Universe, edited by O. Sawada and A. Sugamoto (KEK, Tsukuba, 1979)

  28. [29]

    Neutrino Mass and Spontaneous Parity Nonconservation,

    R. N. Mohapatra and G. Senjanovi´ c, “Neutrino Mass and Spontaneous Parity Nonconservation,” Phys. Rev. Lett.44, 912 (1980), doi:10.1103/PhysRevLett.44.912

  29. [30]

    Neutrino Masses, Mixings and Oscillations in SU(2)×U(1) Models of Electroweak Interactions,

    T. P. Cheng and L. F. Li, “Neutrino Masses, Mixings and Oscillations in SU(2)×U(1) Models of Electroweak Interactions,” Phys. Rev. D22, 2860 (1980), doi:10.1103/PhysRevD.22.2860

  30. [31]

    Neutrino masses and mixings in gauge models with spontaneous parity violation,

    R. N. Mohapatra and G. Senjanovi´ c, “Neutrino masses and mixings in gauge models with spontaneous parity violation,” Phys. Rev. D23, 165 (1981), doi:10.1103/PhysRevD.23.165

  31. [32]

    Seesaw Neutrino Masses Induced by a Triplet of Leptons,

    R. Foot, H. Lew, X. G. He and G. C. Joshi, “Seesaw Neutrino Masses Induced by a Triplet of Leptons,” Z. Phys. C44, 441–444 (1989), doi:10.1007/BF01415558

  32. [33]

    Type-3/2 seesaw mechanism,

    D. Demir, C. Karahan and O. Sargın, “Type-3/2 seesaw mechanism,” Phys. Rev. D104, no.7, 075038 (2021), [arXiv:2105.06539 [hep-ph]]

  33. [34]

    Electroweak Higgs Potentials and Vacuum Stability,

    M. Sher, “Electroweak Higgs Potentials and Vacuum Stability,” Phys. Rept.179, 273–418 (1989), doi:10.1016/0370- 1573(89)90061-6

  34. [35]

    Three-loop \beta-functions for top-Yukawa and the Higgs self-interaction in the Standard Model

    K. G. Chetyrkin and M. F. Zoller, “Three-loop beta functions for the Higgs self-interaction and the top-Yukawa coupling,” JHEP06, 033 (2012), [arXiv:1205.2892 [hep-ph]]

  35. [36]

    Vacuum Stability, Perturbativity, and Scalar Singlet Dark Matter

    M. Gonderinger, Y. Li, H. Patel and M. J. Ramsey-Musolf, “Vacuum Stability, Perturbativity, and Scalar Singlet Dark Matter,” JHEP01, 053 (2010), [arXiv:0910.3167 [hep-ph]]

  36. [37]

    Vacuum stability vs. positivity in real singlet scalar extension of the standard model,

    P. Ghorbani, “Vacuum stability vs. positivity in real singlet scalar extension of the standard model,” Nucl. Phys. B971, 115533 (2021), [arXiv:2104.09542 [hep-ph]]

  37. [38]

    Vacuum stability and scalar masses in the superweak extension of the Standard Model,

    Z. P´ eli and Z. Tr´ ocs´ anyi, “Vacuum stability and scalar masses in the superweak extension of the Standard Model,” Phys. Rev. D106, 055045 (2022), [arXiv:2204.07100 [hep-ph]]

  38. [39]

    Vacuum Stability in the Standard Model and Beyond,

    G. Hiller, T. H¨ ohne, D. F. Litim and T. Steudtner, “Vacuum Stability in the Standard Model and Beyond,” Phys. Rev. D 110, 115017 (2024), [arXiv:2401.08811 [hep-ph]]

  39. [40]

    Impact of massive neutrinos on the Higgs self-coupling and electroweak vacuum stability

    W. Rodejohann and H. Zhang, “Impact of massive neutrinos on the Higgs self-coupling and electroweak vacuum stability,” JHEP06, 022 (2012), [arXiv:1203.3825 [hep-ph]]

  40. [41]

    Higgs Vacuum Stability in $B-L$ extended Standard Model

    A. Datta, A. Elsayed, S. Khalil and A. M. Moretti, “Higgs vacuum stability in theB−Lextended Standard Model,” Phys. Rev. D88, 053011 (2013), [arXiv:1308.0816 [hep-ph]]

  41. [42]

    S. Oda, N. Okada and D. s. Takahashi, ‘Classically conformal U(1)’ extended standard model and Higgs vacuum stability,” Phys. Rev. D92, no.1, 015026 (2015), [arXiv:1504.06291 [hep-ph]]

  42. [43]

    Higgs Vacuum Stability with Vector-like Fermions

    S. Gopalakrishna and A. Velusamy, “Higgs vacuum stability with vectorlike fermions,” Phys. Rev. D99, 115020 (2019), [arXiv:1812.11303 [hep-ph]]

  43. [44]

    Vacuum stability in the Standard Model with vector-like fermions,

    A. Arsenault, K. Y. Cingiloglu and M. Frank, “Vacuum stability in the Standard Model with vector-like fermions,” Phys. Rev. D107, 036018 (2023), [arXiv:2207.10332 [hep-ph]]

  44. [45]

    Gauge singlet scalar as inflaton and thermal relic dark matter

    R. N. Lerner and J. McDonald, “Gauge singlet scalar as inflaton and thermal relic dark matter,” Phys. Rev. D80, 123507 (2009), [arXiv:0909.0520 [hep-ph]]

  45. [46]

    Running Kinetic Inflation

    K. Nakayama and F. Takahashi, “Running kinetic inflation,” JCAP1011, 009 (2010), [arXiv:1008.2956 [hep-ph]]

  46. [47]

    Inflation scenario via the Standard Model Higgs boson and LHC

    A. O. Barvinsky, A. Y. Kamenshchik and A. A. Starobinsky, “Inflation scenario via the Standard Model Higgs boson and LHC,” JCAP11, 021 (2008), [arXiv:0809.2104 [hep-ph]]

  47. [48]

    Inflation with Non-Minimal Coupling: Metric vs. Palatini Formulations

    F. Bauer and D. A. Demir, “Inflation with Non-Minimal Coupling: Metric versus Palatini Formulations,” Phys. Lett. B 665, 222-226 (2008), [arXiv:0803.2664 [hep-ph]]

  48. [49]

    Higgs inflation

    J. Rubio, “Higgs inflation,” Front. Astron. Space Sci.5, 50 (2019), [arXiv:1807.02376 [hep-ph]]

  49. [50]

    Review of particle physics,

    R. L. Workmanet al.(Particle Data Group), “Review of particle physics,” Prog. Theor. Exp. Phys. 2022, 083C01 (2022), doi:10.1093/ptep/ptac097

  50. [51]

    Neutrino Mass Ordering from Oscillations and Beyond: 2018 Status and Future Prospects

    P. F. De Salas, S. Gariazzo, O. Mena, C. A. Ternes and M. T´ ortola, “Neutrino Mass Ordering from Oscillations and Beyond: 2018 Status and Future Prospects,” Front. Astron. Space Sci.5, 36 (2018), [arXiv:1806.11051 [hep-ph]]

  51. [52]

    M. Aaboudet al.[ATLAS], “Search for pair production of heavy vector-like quarks decaying into high-p T Wbosons and top quarks in the lepton-plus-jets final state inppcollisions at √s= 13 TeV with the ATLAS detector,” JHEP08, 048 (2018), [arXiv:1806.01762 [hep-ex]]

  52. [53]

    Search for pair production of vector-like quarks in leptonic final states in proton-proton collisions at √s= 13 TeV,

    A. Tumasyanet al.[CMS], “Search for pair production of vector-like quarks in leptonic final states in proton-proton collisions at √s= 13 TeV,” JHEP07, 020 (2023), [arXiv:2209.07327 [hep-ex]]

  53. [54]

    Vector-like Fermions and Higgs Effective Field Theory Revisited

    C. Y. Chen, S. Dawson and E. Furlan, “Vectorlike fermions and Higgs effective field theory revisited,” Phys. Rev. D96, no.1, 015006 (2017), [arXiv:1703.06134 [hep-ph]]. 19

  54. [55]

    A handbook of vector-like quarks: mixing and single production

    J. A. Aguilar-Saavedra, R. Benbrik, S. Heinemeyer and M. P´ erez-Victoria, “Handbook of vectorlike quarks: Mixing and single production,” Phys. Rev. D88, no.9, 094010 (2013), [arXiv:1306.0572 [hep-ph]]

  55. [56]

    Bounds on the mixing of the down-type quarks with vectorlike singlet quarks,

    L. Lavoura and J. P. Silva, “Bounds on the mixing of the down-type quarks with vectorlike singlet quarks,” Phys. Rev. D 47, 1117 (1993), doi:10.1103/PhysRevD.47.1117

  56. [57]

    Theoretical constraints on models with vectorlike fermions,

    A. Adhikary, M. Olechowski, J. Rosiek and M. Ryczkowski, “Theoretical constraints on models with vectorlike fermions,” Phys. Rev. D110, 075029 (2024), [arXiv:2406.16050 [hep-ph]]

  57. [58]

    Upper bounds on the values of masses in unified gauge theories,

    D. A. Dicus and V. S. Mathur, “Upper bounds on the values of masses in unified gauge theories,” Phys. Rev. D7, 3111 (1973), doi:10.1103/PhysRevD.7.3111

  58. [59]

    Weak Interactions of Ultraheavy Fermions. 2.,

    M. S. Chanowitz, M. A. Furman and I. Hinchliffe, “Weak Interactions of Ultraheavy Fermions. 2.,” Nucl. Phys. B153, 402-430 (1979), doi:10.1016/0550-3213(79)90606-0

  59. [60]

    Improved Con- straints on Primordial Gravitational Waves using Planck, WMAP, and BICEP/Keck Observations through the 2018 Observing Season,

    P. A. R. Adeet al.[BICEP and Keck], “Improved Constraints on Primordial Gravitational Waves using Planck, WMAP, and BICEP/Keck Observations through the 2018 Observing Season,” Phys. Rev. Lett.127, no.15, 151301 (2021), [arXiv:2110.00483 [astro-ph.CO]]

  60. [61]

    T. P. Cheng and L. F. Li,Gauge Theory of Elementary Particle Physics, Oxford University Press (1984)

  61. [62]

    Machacek and M.T

    M.E. Machacek and M.T. Vaughn, ”Two-loop Renormalization Group Equations in a General Quantum Field Theory (I). Gauge Fields”, Nucl. Phys. B222 (1983) 83, doi:10.1016/0550-3213(83)90610-7

  62. [63]

    Neutrino Mass Operator Renormalization Revisited

    S. Antusch, M. Drees, J. Kersten, M. Lindner and M. Ratz, “Neutrino mass operator renormalization revisited,” Phys. Lett. B519, 238-242 (2001), [arXiv:hep-ph/0108005 [hep-ph]]

  63. [64]

    Neutrino Mass Operator Renormalization in Two Higgs Doublet Models and the MSSM

    S. Antusch, M. Drees, J. Kersten, M. Lindner and M. Ratz, “Neutrino mass operator renormalization in two Higgs doublet models and the MSSM,” Phys. Lett. B525, 130-134 (2002), [arXiv:hep-ph/0110366 [hep-ph]]

  64. [65]

    Review of particle physics

    S. Navaset al.(Particle Data Group), “Review of Particle Physics,” Phys. Rev. D110, 030001 (2024), doi:10.1103/PhysRevD.110.030001