2026-02-10 | | Total: 37
We employ a simple model for nuclear modification of ordinary parton densities in a proton to evaluate the Transverse Momentum Dependent gluon and quark distributions in nuclei (nTMDs) within the popular Kimber-Martin-Ryskin/Watt-Martin-Ryskin approach. The model is based on a global analysis of available deep inelastic scattering data for different nuclear targets within the rescaling model, incorporating Fermi motion effects. The derived nTMDs are tested with latest CMS data on inclusive $b$-jet and $B^+$ meson production in proton-lead collisions collected at $\sqrt s = 5.02$ and $8.16$~TeV using the High Energy Factorization framework. We predict the corresponding nuclear medium modification factors to be about of $0.8 - 1.2$ in the probed kinematical region, which is consistent with other estimations. Specially we highlight a possibility to investigate the nuclear modification of parton densities by applying different cuts on the final states in such processes.
It is demonstrated that, within the eikonal approach, the Coulomb corrections to the elastic electromagnetic non-flip and spin-flip proton--nucleus amplitudes are identical when the two amplitudes share the same exponential form factors. This result allows Coulomb corrections to be computed numerically, and with high precision, for both electromagnetic and hadronic elastic $p^{\uparrow}A$ amplitudes in the massless-photon limit, including the effects of soft magnetic photon exchange. The method relies on analytical expressions and numerical integrations over a finite impact-parameter range with nonsingular integrands, providing a practical and systematically controlled framework for phenomenological applications.
The Born-Oppenheimer approximation is widely used to investigate the properties of hydrogen-like systems and doubly heavy hadrons. However, the extent to which this approximation reliably captures the true features of such systems remains an open question. In this work, we adopt the results obtained with the Gaussian expansion method as a benchmark to assess the validity of the Born-Oppenheimer approximation in hadronic systems. We also investigate the dependence of the Born-Oppenheimer approximation results on the choice of trial wave functions. A comprehensive study of the Born-Oppenheimer approximation is carried out by performing calculations using Slater-type functions and Gaussian-type functions as trial wave functions, and by comparing the resulting predictions with those obtained from the Gaussian expansion method. We find that the calculations performed within the Born-Oppenheimer approximation are close to those obtained with the Gaussian expansion method when the heavy-quark mass is relatively small. However, as the heavy-quark mass increases, calculations employing Slater-type functions yield larger values than those from the Gaussian expansion method, whereas those using Gaussian-type functions lead to smaller ones. The use of Slater-type functions generally leads to an enhanced binding energy. The underestimation observed in Born-Oppenheimer approximation calculations with Gaussian-type functions primarily stems from the neglect of non-adiabatic corrections. This comparative study provides deeper insight into the structure of doubly heavy hadrons and clarifies the applicability and limitations of the Born-Oppenheimer treatment in these systems.
The Affleck-Dine mechanism is a leading baryogenesis scenario in which scalar condensates form coherently during inflation along supersymmetric flat directions that are lifted by supersymmetry-breaking effects. We update the viable parameter space for baryogenesis using recent Cosmic Microwave Background constraints on baryon-density isocurvature perturbations, taking the quantum fluctuations of the scalar condensate generated during inflation as initial conditions. We then show that primordial features arising from the inflaton sector can serve as a unique probe of baryogenesis models, whose mechanisms are otherwise difficult to access directly due to their high energy scales. These primordial features leave correlated imprints, such as sharp feature signals and clock signals, on both the curvature and baryon-density isocurvature perturbations, providing direct evidence for the existence of both light and heavy modes involved in the Affleck-Dine mechanism.
The photoleptonic decay $B^-\to γ\ell^-\barν$ is the simplest low-energy process that probes the substructure of the $B$ meson, making it an excellent candidate to determine the parameters of $B$-meson light-cone distribution amplitudes from experimental data. More recently, the decay $B^-\to γ^*(\to \ell'^-\ell'^+)\ell^-\barν$ has received attention as an alternative probe of these parameters. Both decays are described through a common set of hadronic form factors, which if computed within the framework of QCD factorization give rise to the sensitivity to the $B$-meson light-cone distribution amplitudes. Nevertheless, in this case the form factors still receive so-called soft contributions that can only be estimated but not rigorously computed. In this work, we provide results for the QCD factorization expressions for the form factors, including terms at next-to-leading power in the $b$ quark mass and in the photon energy. We further use these results to estimate the soft contributions within a light-cone sum rule setup. Finally, using numerical results obtained from a benchmark model of the light-cone distribution amplitudes, we show that the soft contributions are under significantly better theoretical control at a mildly spacelike photon virtuality.
This study investigates photon entanglement generated from para-positronium decay by analyzing azimuthal correlations after the double Compton scattering with stationary electrons. We introduce a normalized correlation observable $\mathcal{O}_1 = \cos(2φ_1 - 2φ_2)/C_1$ to witness entanglement. In the absence of decoherence, $\langle\mathcal{O}_1\rangle = -1$, corresponding to a maximally entangled Bell state. With decoherence parameterized by $ρ$, the expectation becomes $-(1-ρ)$, allowing direct experimental quantification of coherence loss. A prior symmetry analysis of the Compton scattering process within the quantum field theory (QFT) is provided, which establishes the mirror-symmetric nature of the single-photon angular distribution. We further examine a local hidden-variable theory (LHVT) under the angular-momentum conservation. Imposing the mirror symmetry with respect to the plane defined by the photon spin and momentum leads to a non-negative LHVT prediction for $\langle \sin^2θ_1 \sin^2θ_2 \cos(2φ_1-2φ_2)\rangle$, contradicting the negative QFT prediction value for any $ρ< 1$. Thus, mirror symmetry serves as a novel criterion to exclude LHVT descriptions of the entangled state, whereas without preserving this symmetry, LHVTs can reproduce the correlations.
The observation of collectivity in small and large collision systems challenges our understanding of thermalization and particle production. EPOS4 models this via a dynamical core--corona separation, where high-density regions form a collectively expanding core while low-density regions hadronize via string fragmentation. Its microcanonical core hadronization improves the description of transverse momentum and multiplicity-dependent observables. We present EPOS4 predictions for pp, O-O and Pb-Pb collisions, with and without UrQMD, showing non-universal $\langle p_T\rangle$ scaling, significant hadronic-phase effects, and system-size-dependent $R_{AA}$ suppression. Charged-particle and transverse-energy densities show participant scaling; the transverse energy per charged particle is systematically larger in O--O than in Pb--Pb at comparable participant fraction, indicating a harder effective production in the lighter system. Identified-hadron spectra harden with event multiplicity with mass ordering and increasing core fractions. The mean transverse momentum exhibits a strong system dependence, with the steepest multiplicity evolution in pp, demonstrating that $\langle p_T\rangle$ does not follow universal multiplicity scaling. The $p/π$ ratio shows an enhanced intermediate-$p_T$ region; the suppression of the integrated $p/π$ at the highest Pb--Pb multiplicities is reproduced only with UrQMD, highlighting hadronic-phase effects. The nuclear modification factor shows sizeable suppression in Pb--Pb and substantial suppression in central O--O collisions. Blast-wave fits exhibit the anti-correlation between $T_{\rm kin}$ and $\langleβ_T\rangle$, with UrQMD shifting the parameters towards lower $T_{\rm kin}$ and higher $\langleβ_T\rangle$. These results provide a timely baseline for Run~3 measurements and for constraining the onset of medium-like effects across system size.
We present a light dark matter model with sizable invisible rare meson decays to accommodate possible larger than standard model (SM) predictions for both $B^+\to K^+ν\barν$ by Belle II and $K^+\toπ^+ν\bar ν$ by NA62. Since the emitted neutrinos are unobserved, the excesses may be due to light dark matter pairs in the final states. We show that it is possible to realize this by a UV-complete model based on a two Higgs doublet model. The neutral spin-zero Higgs bosons (scalar or pseudoscalar) mediating dark matter interaction, can also produce some effects at low energies which modify the SM predictions significantly. In particular we find large testable consequences for $B_s - \bar B_s$ mixing due to splitting of scalar and pseudoscalar masses, and also neutron EDM. We find that there is a cancellation due to exchange of neutral spin-zero particle for neutron EDM in general, but QCD running will lift this cancellation.This is generally true for any neutral Higgs contributions. However, we find that such a cancellation does not happen for charged scalar contribution. The allowed CP violating phase in the Yukawa sector can produce a neutron EDM experimental bound.
The effects of isotopic invariance violation in the processes $e^{+}e^{-}\to B\bar{B}$, $e^{+}e^{-}\to B^{*}\bar{B}$, and $e^{+}e^{-}\to B^{*}\bar{B}^{*}$ are considered in the energy range between the thresholds of $B^{+}B^{-}$ and $B_{s}^{0}\bar{B}_{s}^{0}$ production. The analysis is based on taking into account the final-state interaction in a six-channel problem. Our approach allowed us to obtain good agreement with recent Belle-II results for the ratio of the $B^{0}\bar{B}^{0}$ and $B^{+}B^{-}$ production cross sections in $e^{+}e^{-}$ annihilation in the vicinity of $Υ(4S)$. It is shown that at higher energies the ratios of the cross sections for pairs of neutral and charged $B^{(*)}$ mesons can differ significantly from unity. A detailed measurement of this effect will provide evidence that the nontrivial energy dependence of the $B^{(*)}\bar{B}^{(*)}$ production cross sections is a consequence of interference of the particle production amplitudes in the multichannel problem.
Neutrinoless double beta ($0νββ$) decays of spin-1/2 hyperons are investigated in a covariant baryon chiral perturbation theory framework, extended by a $ΔL=2$ operator proportional to the Majorana neutrino mass, where $L$ denotes the lepton number. Within the light Majorana neutrino exchange mechanism, the decay amplitudes are found to emerge at the one-loop level, representing the long-range contribution. The extended-on-mass-shell scheme is employed to renormalize the one-loop amplitudes and restore consistent chiral power counting. Consequently, the differential decay rates for all accessible hyperon $0νββ$ channels are predicted and the corresponding branching ratios are more than 20 orders of magnitude smaller than the current experimental upper bounds. Interestingly, it is found that the leading contribution to hyperon $0νββ$ decay is actually from short-range counterterm operators, as required by the renormalization argument. Neutrinoless transition form factors are proposed to determine this leading contribution through future lattice QCD simulations.
In this study, we derive the longitudinal and transverse gluon self-energies and the corresponding dielectric functions for a nonextensive QGP, based on nonextensive statistical mechanics and a kinetic theory framework. The nonextensive parameter $q$ enters these quantities primarily through the modification of the Debye mass. Utilizing the derived dielectric functions, we then calculate the collisional energy loss for a heavy quark using two established formalisms: the plasma physics-based Thoma-Gyulassy formula and the thermal field theory-originated Kirzhnits-Thoma formula. Our results show that for both formalisms, the collisional energy loss increases with the nonextensive parameter $q$ with this enhancement being more significant at higher incident quark momenta and suppressed for a heavier quark mass. The energy loss predicted from the Kirzhnits-Thoma formula is substantially larger than that from the Thoma-Gyulassy formula, and the nonextensive effect on the energy loss is more pronounced in the former. Furthermore, the mass suppression of the nonextensive effect on the energy loss is weaker in the Kirzhnits-Thoma approach. These calculations demonstrate that nonextensive statistics can significantly alter the energy loss in the QGP.
We consider the FIMPs scenario in a two-component dark matter model with $Z_2 \times Z_4$ symmetry, where a singlet scalar $S$ and a Majorana fermion $χ$ are introduced as dark matter candidates. We also introduce another singlet scalar $S_0$ with a non-zero vacuum expectation value to the SM so that the fermion dark matter can obtain mass after spontaneous symmetry breaking. The model admits six free parameters in the decoupling limit: three masses and three dimensionless parameters. Depending on the mass hierarchies between dark matter particles with half of the new Higgs mass, the DM relic density will be determined by different channels, where $χ$ and $S$ production can be generated individually. We numerically study the relic density as a function of the model's free parameters and determine the regions consistent with the dark matter constraint for four possible cases. Our results show that this scenario is viable over a wide range of couplings and dark matter masses, where the coupling $λ_{ds}$ can be as tiny as $\mathcal{O}(10^{-20})$ level. We stress that even for such tiny couplings, the new Higgs can still play a dominant role in determining dark matter production.
In this work we investigate the production of a dark fermionic particle $χ$ in the neutrino scattering experiments. In the framework of effective field theory, such process can be induced by the effective four-fermion interactions involving neutrinos, the dark particle $χ$ and standard model particles. In particular, we examine the constraints on the effective couplings from the neutrinos scattering off nuclei in the COHERENT and CONUS+ experiments as well as the prospects at the DUNE near detector from neutrino-electron scattering. It turns out the current COHERENT and CONUS+ constraints on the cutoff scales are less stringent than those from the existing Large Hadron Collider data. However, the DUNE near detector could probe the cutoff scales beyond the existing CHARM II and LEP limits up to roughly 1 TeV, for the dark particle mass up to roughly 50 MeV.
The analytical expressions of the shear viscosity of both one and two particle species with Boltzmann statistics and $2 \rightarrow 2$ elastic scatterings are known from the Chapman-Enskog method and have been shown to be quite accurate. The expression for a multi-species hadronic gas under $2 \rightarrow 2$ elastic scatterings is also known. Here we use the Chapman-Enskog method to derive the shear viscosity of a massless quark-gluon gas of $N_f$ quark flavors in chemical equilibrium subjected to all $2 \rightarrow 2$ parton scatterings including for the first time inelastic scatterings. We then verify the relation in a general single-species limit, where the shear viscosity of the quark-gluon gas should reduce to the result for a single particle species. In addition, we show the explicit analytical result in terms of the seven independent cross sections for the special case of isotropic and energy-independent cross sections. The analytical relations derived here can be useful for determining the shear viscosity of parton transport models with any $2 \rightarrow 2$ scattering cross sections. They can also be coupled with finite temperature QCD cross sections to help study the shear viscosity of the quark gluon plasma.
Searches for Lorentz invariance violation (LIV) in the neutrino sector have traditionally focused on non-standard neutrino oscillations induced by LIV in vacuum. In this work, however, we study anisotropic LIV in matter. First, we review vacuum LIV phenomenology, explaining the energy and direction dependence of sidereal modulations for anisotropic coefficients in the Standard Model Extension. We then demonstrate that for high-energy neutrinos, the interplay between anisotropic LIV operators and the Earth's matter potential produces, distinct, observable signatures absent in the vacuum case. We identify a crossover regime where the energy-dependent LIV Hamiltonian becomes comparable to the matter potential, leading to strong interference effects. By analyzing the propagation of neutrinos through a realistic Earth model, we establish three key phenomenological consequences: (1) direction-dependent resonant enhancements of oscillation probabilities, (2) a macroscopic breakdown of neutrino-antineutrino symmetry for CPT-even operators, and (3) a significant increase of the $ν_τ$ flux due to LIV-driven injection of high-energy neutrinos into the $τ$ regeneration cycle. These results highlight that accounting for the interplay between LIV and matter is essential for future LIV searches at large-scale neutrino telescopes.
Quantum estimation theory provides a fundamental framework for analyzing how precisely physical parameters can be estimated from measurements. Neutrino oscillations are characterized by a set of parameters inferred from experiments conducted in different production and detection environments. The two solar oscillation parameters, $Δm^2_{21}$ and $θ_{12}$, can be estimated using both solar neutrino experiments and reactor neutrino experiments. In reactor experiments, neutrinos are detected after coherent vacuum evolution, while solar neutrinos arrive at the detector as incoherent mixtures. In this work, we use Quantum Fisher Information (QFI) to quantify and compare the information content accessible in these two experimental setups. We find that for reactor neutrinos, flavor measurements saturate the QFI bound for both parameters over specific energy ranges, demonstrating their optimality and explaining the high precision achieved by these experiments. In contrast, for solar neutrinos the phase-based contribution to the QFI, originating from the quantum coherence, is absent, rendering the estimation of $Δm_{21}^2$ purely population-based and effectively classical, while the QFI for $θ_{12}$ is dominated by basis rotation at high energies and is nearly saturated by flavor measurements. Consequently, solar neutrino experiments are intrinsically more sensitive to $θ_{12}$ than to $Δm_{21}^2$. This analysis highlights a fundamental distinction between the two estimation problems and accounts for their differing achievable precisions.
We study the structure of the $X(3915)$ meson in a quark model and explore how its production in heavy ion collisions depends on its internal structure. We first analyze the $X(3915)$ as a $c\bar{c}s\bar{s}$ state and solve the Hamiltonian with color-spin interactions within the quark model. We find that the ground state of the $c\bar{c}s\bar{s}$ with total spin 0 obtained from the quark model analysis favors a separated $D_s \bar{D}_s$ state. To probe its structure further, we study its production in relativistic heavy ion collisions for various proposed configurations. We calculate the transverse momentum distributions and yields for the $X(3915)$ assuming its structure to be either a charmonium, a tetraquark, or a hadronic molecular state. We argue that by measuring the transverse momentum distributions and yields of the $X(3915)$ produced in heavy ion collisions, one can identify the structure of the $X(3915)$.
We present a systematic analysis of Coulomb corrections for leptonic ($B^0_{d,s}\to \ell^+\ell^-$), semileptonic ($B^0_{d,s}\to h^0\,\ell^+\ell^-$, $B^0_{d,s}\to V^0\ell^+\ell^-$) and radiative leptonic ($B^0_{d,s}\to γ\ell^+\ell^-$) decays of neutral $B$-mesons. The relativization of the Coulomb factor was performed by comparing the Gamow-Sommerfeld-Sakharov factor, the exact relativistic approach of Crater-Alstine-Sazdjian applied by us to scalar systems, and well-known one-loop QED calculations. Coulomb corrections are calculated for differential, angular, and double-differential distributions, as well as for partial decay widths. For the $B_s^0 \to μ^+μ^-$ channel, Coulomb corrections improve the prediction of the partial width to $δ= |\mathcal{B}^{(exp)} - \mathcal{B}^{(theory)}|/\mathcal{B}^{(exp)} = 2\%$. This improvement brings the prediction closer to the LHCb/CMS experimental results within the current experimental (11\%) and theoretical (5\% lattice QCD) errors. In the decays $B^0\to K^0μ^+μ^-$ and $B^0 \to K^{0*}μ^+μ^-$, Coulomb effects also reduce the discrepancies between theoretical predictions and experimental data (to less than $δ= 1\%$ and from $δ= 11\%$ to $δ= 4\%$ respectively). Finally, for the decays involving $τ$-leptons, the Coulomb correction $\mathcal{K} = \mathcal{B}^{(Coulomb)}/ \mathcal{B}^{(free)}$ reaches 4\%. While currently smaller than the dominant form-factor uncertainties and experimental errors, the Coulomb correction represents a non-negligible systematic effect. It should be accounted for in the high-precision era of $B$-physics, where such effects may become significant for the interpretation of potential New Physics signals.
We look at our since long studied model for dark matter as being pearls of a speculated new vacuum containing highly compressed ordinary matter, with so much ordinary in it that the content of ordinary matter in the dark matter pearls dominate. Most dark matter models have the dark matter consisting mainly of new-physics-matter such as WIMPs being supersymmetric partners of possibly known particles or, as in Maxim Khlopovs model, a doubly negatively charged new-physics-particle with a helium nucleus attached. But usually the new-physics matter makes up weight-wise the major content. It is only in our model that the ordinary matter content in the dark matter dominates. We here expose some weak phenomenological evidence that, in truth, dark matter should be of the type with a dominant component of ordinary matter (weight-wise), thus favoring as the typical example our previously so much studied vacuum type 2 model. The main such evidence is that we manage a fit to data in which the 3.5 keV X-rays, presumed to result from dark matter, come both from collisions of dark matter with dark matter and from dark matter with ordinary matter! Both mechanisms are of so similar an order of magnitude that they are both seen, indicating that their similarity is due to a significant similarity between dark with ordinary matter. The fact that the amounts of ordinary and dark matter only deviate by a factor 6 points in the same direction. Using the information obtained from this fitting, we develop our speculation that the main content weight-wise of dark matter is ordinary matter to the very DAMA experiment. Actually we found three spots on the sky in which we fit the observed production of 3.5 keV X-rays with ordinary plus dark scattering.
We study an exotic type of single spin asymmetry in unpolarized electron-proton scattering, in which the outgoing electron momentum exhibits a left-right asymmetry relative to the transverse spin of the leading baryon $B$ in the target fragmentation region. We lay out two theoretical frameworks for describing this effect: The twist-three fracture function at high-$Q^2$ and the spin-dependent odderon in the high energy limit.
This PhD Thesis is devoted to the study of the emergence of attractors, universality and collectivity in ultra-relativistic collisions by means of relativistic kinetic theory. After an introduction about Quantum Chromodynamics (QCD), Quark-Gluon Plasma (QGP) and the importance of heavy-ion collisions to investigate both, we give an overview about the two main models able to describe the hot QCD matter collective behaviour, namely kinetic theory and hydrodynamics. Afterwards, the Relativistic Boltzmann Transport (RBT) model, which has been employed to obtain most part of the results of this thesis, is carefully described, from the numerical and physical perspectives. The study of attractors and universality proceeds then by starting from a simple one-dimensional massless model, moving to increasingly more complex scenarios, involving the full 3+1D setup, non-conformal systems and realistic event-by-event fluctuations. Particular attention is paid to the physical scales which govern the system collectivity and their interplay. We show that a very good description of collective behaviour can be carried out by means of a few variables which characterise the systems under study.
We study a strongly interacting complex scalar dark matter candidate $(χ)$, subject to an attractive potential mediated by a vector boson $A^μ$. Such interactions allow $χ$ to form bound states: SIMPonium. In this work, we systematically investigate the bound state dynamics of our dark sector, including the formation and decay of SIMPonium and it's influence on the thermal history. Our analysis shows in absence of any bound state formation $χ$ freezes out at $x\approx 16$ and the presence of SIMPonium in the thermal bath slightly modifies the freeze out behaviour of the free $χ$ particles, which freezes out at $x \approx 20$. While the bound state itself remains in chemical equilibrium for a longer duration and freezes out at a significantly later time, $x \approx 250$. We compute the indirect energy spectra arising from free dark matter annihilation and SIMPonium decay, where the resulting Standard Model particles subsequently produce both final state radiation and radiative decay spectra. We find that the total differential photon flux from dark matter with mass $m_χ= 150\,\mathrm{MeV}$ lies in the range $E_γ^2 \frac{dφ}{dΩ\, dE_γ} \in [10^{-27},\,10^{-17}] \,\mathrm{MeV\,cm^{-2}\,s^{-1}\,sr^{-1}}$ , for photon energies in the interval $E_γ\in [10^{-3},\,10^{2}] \,\mathrm{MeV}$. The predicted signal is therefore exceedingly feeble and remains well below the sensitivity of current experimental facilities.
We reinterpret classic radiochemical solar-neutrino measurements as ``rate meters'' for additional, non-negative capture-like contributions induced by fermionic dark matter absorption. Using the chlorine and gallium production-rate data, we build a Bayesian likelihood that accounts for the dominant uncertainties in the solar-neutrino capture-rate prediction (solar fluxes, oscillation parameters, and capture cross sections). Solar-model metallicity systematics are made explicit by presenting results for both the B16--GS98 and B16--AGSS09met solar-model realizations. From the 1D marginalized posteriors of the joint $(R_{χ,\mathrm{Cl}},R_{χ,\mathrm{Ga}})$ analysis, we obtain 90\% upper limits on additional capture-like rate contributions, dominated by chlorine: $R_{χ,\mathrm{Cl},90}\simeq 0.388~\mathrm{SNU}$ (B16--GS98) and $0.588~\mathrm{SNU}$ (B16--AGSS09met). In the charged-current V--A benchmark, we map these constraints onto upper bounds on $y\equiv m_χ^2/(4πΛ^4)$ for $m_χ$ above the ${}^{71}$Ga and ${}^{37}$Cl capture thresholds, using a pep-normalized operator mapping anchored to solar-neutrino capture inputs, where $m_χ$ is the dark matter mass and $Λ$ is the effective scale suppressing the charged-current operator. At $m_χ\simeq 1~\mathrm{MeV}$, we find $y_{90}\simeq 4.88\times 10^{-49}~\mathrm{cm}^2$ (B16--GS98) and $7.08\times 10^{-49}~\mathrm{cm}^2$ (B16--AGSS09met). These radiochemical bounds are complementary to xenon-based absorption searches and collider interpretations by probing distinct nuclear targets with minimal reliance on spectral reconstruction.
In particle physics model building, a discrete $\mathcal{Z}_2$ symmetry is often spontaneously broken for phenomenological reasons. When this breaking occurs dynamically in the early Universe, stable domain wall networks are formed, which can eventually dominate the cosmic energy density. To avoid this problem, explicit $\mathcal{Z}_2$-breaking terms in the scalar potential are usually introduced in an ad hoc manner. In this Letter, we show that if the same $\mathcal{Z}_2$ symmetry is also responsible for generating light Dirac neutrino masses, such explicit breaking terms can instead arise radiatively from the particles involved in the Dirac mass generation. We find that the resulting bias term scales inversely with the cube of the Dirac neutrino mass, leading to a gravitational wave spectrum proportional to the sixth power of the Dirac neutrino mass. This establishes a nontrivial connection between the Dirac seesaw scale, the domain wall annihilation epoch, and the resulting stochastic gravitational wave signal. We further demonstrate that a wide range of Dirac seesaw scales can be probed by upcoming gravitational wave and cosmic microwave background experiments, while part of the parameter space simultaneously explains the observed baryon asymmetry via Dirac leptogenesis.
The minimal scotogenic model where small neutrino masses arise via radiative seesaw, is known to provide a unified framework for neutrino mass and origin of matter via leptogenesis. However if the inflation reheat temperature of the universe is below the sphaleron reheating temperature, then leptogenesis fails and one way to understand the origin of matter would be to add an effective interaction involving the right handed neutrino (RHN) $N$ of the form $\frac{1}{Λ^2}N u_Rd_Rd_R$. This model can lead to observable neutron-anti-neutron ($n-\bar{n}$) oscillation. We show that if RHNs are produced non-thermally, we can get a cosmological upper limit on the transition time $τ_{n-\bar{n}}$, which is within the reach of the planned ESS HIBEAM/NNBAR experiment. The proton stability is guaranteed by the scotogenic $Z_2$ invariance, which prevents the appearance of the Dirac mass term for the neutrino