TeV Scale Research Articles

Higgs physics at ILC

International Linear Collider (ILC) is a viable option for realization of a future Higgs factory. With proven accelerator technology and technically mature detector design, ILC offers well-understood physics program in the top, EW and Higgs sectors. Clean environment of electron–positron collisions, together with beam polarization of both beams and center-of-mass reach at a TeV scale, will make ILC a precision tool to probe realizations of new physics. The above will be illustrated on examples of fully simulated measurements in the Higgs sector, including Higgs couplings, Higgs self-coupling, CP violation measurements and Higgs exotic decay studies at ILC.

Di-Higgs production via axion-like particles

Due to the pseudo-scalar nature of the axion-like particle (ALP), the CP-conserving production of two Higgs bosons via the ALP necessarily involves an additional Z or γ boson. We examine the existing constraints from di-Higgs searches at Run 2 of the LHC and find that, despite the presence of extra objects in the final state, these searches are sensitive to a combination of ALP couplings to gluons and three-bosons in the TeV scale range. Additionally, we propose a specialized search strategy incorporating an energetic leptonic Z boson. This refined ALP-induced production process would allow for the identification of the h h → 4 b-jet final state and could potentially probe the TeV scale using data from Run 2 of the LHC. This production process can also occur through a coupling between the top quark and the ALP. We translate the current constraints on di-Higgs production into new limits on the ALP-top coupling.

Spontaneous symmetry breaking, gauge hierarchy, and electroweak vacuum metastability

The so-called metastability bound asserts that an unnaturally small Higgs mass is a necessary condition for electroweak vacuum metastability, offering a new approach toward solving the hierarchy problem. So far, this result relies on the assumption of a negative Higgs mass parameter, or equivalently, on electroweak spontaneous symmetry breaking. We derive a new, corresponding bound for the case of a positive mass parameter. When the new bound is significantly more restrictive than or comparable to its established counterpart, it may offer an explanation for the sign of the Higgs mass parameter, and thus, spontaneous symmetry breaking itself. New physics at scales O(1–10) TeV can lower these bounds as far as the TeV scale. As an illustration, we consider vacuum stability in the presence of additional TeV-scale fermions with Yukawa couplings to the Higgs, as well as a dimension-six term parametrizing new physics in the UV. This scenario requires new physics that couples strongly to the Higgs, and can potentially be probed at future colliders. Finally, to allow for comparison with concrete mechanisms predicting metastability, we provide the mass-dependent lifetime of the electroweak vacuum for this model. Published by the American Physical Society 2024

Interplay between vector-like lepton and seesaw mechanism: oblique corrections

The non-vanishing neutrino mass strongly hints the existence of right-handed neutrinos (RHNs), singlets of the standard model (SM). However, they are highly decoupled from the SM and difficult to probe. In this work, we consider the Majorana RHNs from the type-I seesaw mechanism may well mix with the heavy neutral lepton dwelling in certain vector-like lepton (VLL), thus acquiring a sizable electroweak charge. Such a simple scenario yields many interesting consequences, and the imprint on oblique corrections, well expected from the mass splitting between components of VLL by virtue of VLL-RHN mixing, is our focus here. We analytically calculate the Peskin-Takeuchi parameters S, T and U with full details, carefully treating the Majorana loop to obtain the self consistent expressions free of divergence. Then, we constrain on the VLL-RHN system which only gives a sizable T parameter using the PDG-2021 data and CDF-II data, separately, by imposing T ≲ O\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{O} $$\\end{document}(0.1). It is found that for the RHN and VLL below the TeV scale, with a properly large mixing, stands in the frontier of the electroweak precision test such as W-boson mass.

Common origin of dark matter and leptogenesis in U(1)B−L

In this paper, we investigate the common parameter space of dark matter and leptogenesis in the U(1)B−L symmetry. This model involves a complex scalar ϕ, sterile neutrinos N, and Majorana dark matter χ, where only dark matter χ is charged under the Z2 symmetry. Masses of N and χ are generated via the Yukawa interactions to ϕ after breaking of the U(1)B−L symmetry. TeV scale sterile neutrinos N are responsible for the generation of baryon asymmetry through the resonance leptogenesis mechanism. The new particles in the U(1)B−L have a significant impact on the dilution of N, thus on leptogenesis. Meanwhile, the annihilation processes of dark matter χ are almost identical to that of N, which indicates that both leptogenesis and dark matter are closely related to satisfying the observed results simultaneously. Under various theoretical and experimental constraints, the viable common parameter space of dark matter and leptogenesis is obtained for both global and local U(1)B−L symmetry.

Probing bottom-associated production of a TeV scale scalar decaying to a top quark and dark matter at the LHC

A minimal non-thermal dark matter model that can explain both the existence of dark matter and the baryon asymmetry in the universe is studied. It requires two color-triplet, iso-singlet scalars with OTeV\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{O}\\left(\ extrm{TeV}\\right) $$\\end{document} masses and a singlet Majorana fermion with a mass of OGeV\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{O}\\left(\ extrm{GeV}\\right) $$\\end{document}. The fermion becomes stable and can play the role of the dark matter candidate. We consider the fermion to interact with a top quark via the exchange of QCD-charged scalar fields coupled dominantly to third generation fermions. The signature of a single top quark production associated with a bottom quark and large missing transverse momentum opens up the possibility to search for this type of model at the LHC in a way complementary to existing monotop searches.

Ladder top-quark condensation imprints in supercooled electroweak phase transition

The electroweak (EW) phase transition in the early Universe might be supercooled due to the presence of the classical scale invariance involving Beyond the Standard Model (BSM) sectors and the supercooling could persist down till a later epoch around which the QCD chiral phase transition is supposed to take place. Since this supercooling period keeps masslessness for all the six SM quarks, it has simply been argued that the QCD phase transition is the first order, and so is the EW one. However, not only the QCD coupling but also the top Yukawa and the Higgs quartic couplings get strong at around the QCD scale due to the renormalization group running, hence this scenario is potentially subject to a rigorous nonperturbative analysis. In this work, we employ the ladder Schwinger-Dyson (LSD) analysis based on the Cornwall-Jackiw-Tomboulis formalism at the two-loop level in such a gauge-Higgs-Yukawa system. We show that the chiral broken QCD vacuum emerges with the nonperturbative top condensate and the lightness of all six quarks is guaranteed due to the accidental U(1) axial symmetry presented in the top-Higgs sector. We employ a quark-meson model-like description in the mean field approximation to address the impact on the EW phase transition arising due to the top quark condensation at the QCD phase transition epoch. In the model, the LSD results are encoded to constrain the model parameter space. We then observe the cosmological phase transition of the first-order type and discuss the induced gravitational wave (GW) productions. We find that in addition to the conventional GW signals sourced from an expected BSM at around or over the TeV scale, the dynamical topponium-Higgs system can yield another power spectrum sensitive to the BBO, LISA, and DECIGO, etc.

Gamma-ray signal from ZN≥3 dark matter-companion models

In Ref. [10], we proposed to replace the final dark matter (DM) particle in the semi-annihilation mode DM+DM→antiDM+Higgs boson with its ZN≥3 companion, thus reducing DM number density without DM-nucleon scattering. In this work, we study the indirect detection signals from DM annihilation, the Higgs boson pair with one of them from the companion decay being on- or off- shell, depending on the DM-companion mass splitting. We generate the photon spectrum by using PYTHIA8 and study the properties of the spectrum, to find that the hard part of the spectrum in our model is mainly shaped by the direct Higgs boson and thus does not differ much from that of the conventional semi-annihilation mode. Using the Fermi-LAT data of dwarf spheroidal galaxies, we derive the current limit of the DM annihilation cross section for DM+DM→companion⁎+Higgs boson, and for the relatively light DM, it reaches the typical thermal cross section. However, for the TeV scale DM, we have to rely on the Cherenkov Telescope Array, which is able to rule out the whole parameter space except for the coannihilation region.

Dark branes for dark matter

We propose a setup for the origin of dark matter based on spacetime with a warped extra dimension and three branes: the Planck brane, the TeV brane, at a (few) TeV scale ρT, and a dark brane, at a (sub-)GeV scale ρ1≲100 GeV≪ρT. The Standard Model (SM) is localized in the TeV brane, thus solving the Higgs hierarchy problem, while the dark matter χ, a Dirac fermion with mass mχ<ρ1, is localized in the dark brane. The radion, with mass mr<mχ, interacts strongly [∼mχ/ρ1∼O(1)] with dark matter and very weakly (∼mfρ1/ρT2≪1) with the Standard Model matter f. The generic conflict between the bounds on its detection signatures and its proper relic abundance is avoided as dark matter annihilation is p-wave suppressed. The former is determined by its very weak interactions with the SM and the latter by its much stronger annihilation into radions. Therefore, there is a vast range in the dark matter’s parameter space where the correct relic abundance is achieved consistently with the existing bounds. Moreover, for the dark brane with ρ1≲3 GeV, a confinement/deconfinement first order phase transition, where the radion condensates, produces a stochastic gravitational wave background at the nanohertz frequencies, which can be identified with the signal detected by the Pulsar Timing Array (PTA) experiments. In the PTA window, for 0.15≲mχ≲2 GeV the relic abundance is reproduced and all constraints are satisfied. Published by the American Physical Society 2024

Minimal U(1) two-Higgs-doublet models for quark and lepton flavor

In the context of the 2HDM, and assuming that neutrinos acquire masses via the Weinberg operator, we perform a systematic analysis to determine the minimal quark and lepton flavor patterns, compatible with masses, mixing and CP violation data, realizable by Abelian symmetries. We determine four minimal models for quarks, where the number of independent parameters matches the number of observables. For the lepton sector, three minimal predictive models are identified. Namely, we find scenarios with a preference for the upper/lower octant of the θ23 atmospheric mixing angle, that exhibit lower bounds on the lightest neutrino masses currently probed by cosmology and testable at future neutrinoless double beta decay experiments, even for a normally ordered neutrino masses. We investigate the phenomenology of each model taking into account all relevant theoretical, electroweak precision observables, scalar sector constraints, as well as stringent quark flavor processes such as B¯→Xsγ, Bs→μ−μ+ and meson oscillations, and the charged lepton flavor-violating decays eα−→eβ−eγ+eδ− and eα→eβγ. We show that, in some cases, Abelian flavor symmetries provide a natural framework to suppress flavor-changing neutral couplings and lead to scenarios featuring heavy neutral/charged scalar masses below the TeV scale within the reach of current experiments. Published by the American Physical Society 2024

Interplay of freeze-in and freeze-out: Lepton-flavored dark matter and muon colliders

We study a lepton-flavored dark matter model and its signatures at a future muon collider. We focus on the less-explored regime of feeble dark matter interactions, which suppresses the dangerous lepton-flavor-violating processes, gives rise to dark matter freeze-in production, and leads to long-lived particle signatures at colliders. We find that the interplay of dark matter freeze-in and its mediator freeze-out gives rise to an upper bound of around TeV scales on the dark matter mass. The signatures of this model depend on the lifetime of the mediator and can range from generic prompt decays to more exotic long-lived particle signals. In the prompt region, we calculate the signal yield, study useful kinematics cuts, and report tolerable systematics that would allow for a 5σ discovery. In the long-lived region, we calculate the number of charged tracks and displaced lepton signals of our model in different parts of the detector and uncover kinematic features that can be used for background rejection. We show that, unlike in hadron colliders, multiple production channels contribute significantly, which leads to sharply distinct kinematics for electroweakly charged long-lived particle signals. Ultimately, the collider signatures of this lepton-flavored dark matter model are common among models of electroweak-charged new physics, rendering this model a useful and broadly applicable benchmark model for future muon collider studies that can help inform work on detector design and studies of systematics. Published by the American Physical Society 2024

Scaling laws of elastic proton-proton scattering differential cross sections

We show that elastic scattering pp→pp differential cross sections as function of the four-momentum transfer square |t| have a universal property, such that the ratio of bump-to-dip positions is constant from the energies of the ISR to the LHC, from tens of GeV and up to the TeV scale. We explore this property to compare the geometrical scaling observed at the ISR with the recently proposed scaling law at the LHC. We argue that, at the LHC, within present experimental uncertainties, there is in fact a family of scaling laws. We discuss the constraints that the scaling laws impose on the parametrization of the elastic pp→pp scattering amplitude.

Larmor power limit for cyclotron radiation of relativistic particles in a waveguide

Cyclotron radiation emission spectroscopy (CRES) is a modern technique for high-precision energy spectroscopy, in which the energy of a charged particle in a magnetic field is measured via the frequency of the emitted cyclotron radiation. The He6-CRES collaboration aims to use CRES to probe beyond the standard model physics at the TeV scale by performing high-resolution and low-background beta-decay spectroscopy of 6He and 19Ne . Having demonstrated the first observation of individual, high-energy (0.1–2.5 MeV) positrons and electrons via their cyclotron radiation, the experiment provides a novel window into the radiation of relativistic charged particles in a waveguide via the time-derivative (slope) of the cyclotron radiation frequency, dfc/dt . We show that analytic predictions for the total cyclotron radiation power emitted by a charged particle in circular and rectangular waveguides are approximately consistent with the Larmor formula, each scaling with the Lorentz factor of the underlying e± as γ 4. This hypothesis is corroborated with experimental CRES slope data.

Experimental overview of recent study on anomalous chiral effects at TeV scale

Experimental overview of recent study on anomalous chiral effects at TeV scale

SMEFT matching to Z′ models at dimension eight

Heavy neutral gauge bosons arise in many motivated models of beyond the Standard Model Physics. Experimental searches require that such gauge bosons are above the TeV scale in most models which means that the tools of effective field theories, in particular the Standard Model effective field theory (SMEFT), are useful. We match the SMEFT to models with heavy Z′ bosons, including effects of dimension-8 operators, and consider the restrictions on model parameters from electroweak precision measurements and from Drell Yan invariant mass distributions and forward-backward asymmetry, AFB, measurements at the LHC. The results demonstrate the model dependence of the resulting limits on SMEFT coefficients and the relatively small impact of including dimension-8 matching. In all cases, the limits from invariant mass distributions are stronger than from AFB measurements in the Z′ models we consider. Published by the American Physical Society 2024

The flavor of a minimal composite S1 leptoquark and the b→cτν¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ b\ o c\ au \\overline{\ u} $$\\end{document} anomaly

Several experiments have measured a deviation in B → D(∗) semileptonic decays, that point to new physics at the TeV scale violating lepton flavor universality. A scalar leptoquark S1, with a suitable structure of couplings in flavor space, is known to be able to solve this anomaly modifying b→cτν¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ b\ o c\ au \\overline{\ u} $$\\end{document}. In the context of composite Higgs models, we consider a theory containing H and S1 as Nambu-Goldstone bosons (NGBs) of a new strongly interacting sector, with ordinary resonances at a scale O10\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{O}(10) $$\\end{document} TeV. Assuming anarchic partial compositeness of the Standard Model (SM) fermions we calculate the potential of the NGBs that is dominated by the fermions of the third generation, we compute RD∗\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {R}_{D^{\\left(\\ast \\right)}} $$\\end{document} and estimate the corrections to flavor observables by the presence of S1. We find that the SM spectrum and mS1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {m}_{S_1} $$\\end{document} ∼ TeV can be obtained with a NGB decay constant of order ∼ 5 TeV. We obtain a robust correlation between the main corrections to RD∗\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {R}_{D^{\\left(\\ast \\right)}} $$\\end{document}, BK∗νν\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {B}_{K^{\\left(\\ast \\right)}\ u \ u} $$\\end{document} and gτ/gμ, that leads to a sever bound on RD∗\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {R}_{D^{\\left(\\ast \\right)}} $$\\end{document}, roughly 2σ below the experimental value. Besides the bounds on the flavor observables gτW\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {g}_{\ au}^W $$\\end{document}, BR(τ → μγ) and ∆mBs\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\Delta {m}_{B_s} $$\\end{document} are saturated, with the first one requiring a coupling between resonances g* ≲ 2, whereas the second one demands mS1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {m}_{S_1} $$\\end{document} ≳ 1.7 TeV, up to corrections of O1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{O}(1) $$\\end{document}.

Gravitational Waves from First-Order Phase Transition in an Electroweakly Interacting Vector Dark Matter Model

Abstract We discuss gravitational waves (GWs) in an electroweakly interacting vector dark matter (DM) model. In the model, the electroweak gauge symmetry is extended to SU(2)$_0 \times$SU(2)$_1 \times$SU(2)$_2 \times$U(1)$_Y$ and spontaneously broken into SU(2)$_L \times$U(1)$_Y$ at TeV scale. The model has an exchange symmetry between SU(2)$_0$ and SU(2)$_2$. This symmetry stabilizes some massive vector bosons associated with the spontaneous symmetry breaking described above, and an electrically neutral one is a DM candidate. In a previous study, it was found that the gauge couplings of SU(2)$_0$ and SU(2)$_1$ are relatively large to explain the measured value of the DM energy density via the freeze-out mechanism. With the large gauge couplings, the gauge bosons potentially have a sizable effect on the scalar potential. In this paper, we focus on the phase transition of SU(2)$_0 \times$SU(2)$_1 \times$SU(2)$_2 \rightarrow$ SU(2)$_L$. We calculate the effective potential at finite temperature and find that the phase transition is first-order and strong in a wide range of the parameter space. The strong first-order phase transition generates GWs. We calculate the GW spectrum and find that it will be possible to detect the GWs predicted in the model by future space-based GW interferometers. We explore the regions of the parameter space probed by the GW detection. We find that the GW detection can probe the region where the mass of $h^{\prime }$, a CP-even scalar in the model, is a few TeV.

A closer look in the mirror: reflections on the matter/dark matter coincidence

We argue that the striking similarity between the cosmic abundances of baryons and dark matter, despite their very different astrophysical behavior, strongly motivates the scenario in which dark matter resides within a rich dark sector parallel in structure to that of the standard model. The near cosmic coincidence is then explained by an approximate ℤ2 exchange symmetry between the two sectors, where dark matter consists of stable dark neutrons, with matter and dark matter asymmetries arising via parallel WIMP baryogenesis mechanisms. Taking a top-down perspective, we point out that an adequate ℤ2 symmetry necessitates solving the electroweak hierarchy problem in each sector, without our committing to a specific implementation. A higher-dimensional realization in the far UV is presented, in which the hierarchical couplings of the two sectors and the requisite ℤ2-breaking structure arise naturally from extra-dimensional localization and gauge symmetries. We trace the cosmic history, paying attention to potential pitfalls not fully considered in previous literature. Residual ℤ2-breaking can very plausibly give rise to the asymmetric reheating of the two sectors, needed to keep the cosmological abundance of relativistic dark particles below tight bounds. We show that, despite the need to keep inter-sector couplings highly suppressed after asymmetric reheating, there can naturally be order-one couplings mediated by TeV scale particles which can allow experimental probes of the dark sector at high energy colliders. Massive mediators can also induce dark matter direct detection signals, but likely at or below the neutrino floor.

Right-handed neutrino as a common progenitor of baryon number asymmetry and dark matter

Cosmological and astrophysical observations suggest that both energy densities of baryon and dark matter take the same order values in the present Universe. We propose a scenario to give an answer for this problem in a scotogenic model and its extension. The model naturally provides dark matter candidates as its crucial ingredients to explain the small neutrino mass. If a neutral component of an inert doublet scalar plays a role of dark matter and leptogenesis occurs at TeV scales, both dark matter abundance and baryon number asymmetry could be explained with a same mother particle. Coincidence between the order of their energy densities might be understood through such a background. Published by the American Physical Society 2024

Leptoquark-mediated two-loop neutrino mass in unified theory

Scalar leptoquarks naturally arise within unified theories, offering a promising avenue for addressing one of the most significant challenges of the Standard Model–the existence of non-zero neutrino masses. In this work, we present a unified theory based on the SU(5) gauge group, where neutrino mass appears at the two-loop level via the propagation of scalar leptoquarks. Due to the unified framework, the charged fermion and neutrino masses and mixings are entangled and determined by a common set of Yukawa couplings. These exotic particles not only shed light on the neutrino mass generation mechanism but also help to achieve the unification of gauge couplings and are expected to lead to substantial lepton flavor violating rates, offering tangible opportunities for experimental verification. Reproducing the observed neutrino mass scale necessitates that a set of leptoquarks reside a few orders of magnitude below the unification scale–a specific feature of the proposed scenario. Moreover, maximizing the unification scale implies TeV scale new physics states, making them accessible at colliders. The diverse roles that leptoquarks play highlight the elegance and predictive ability of the proposed unified model.