{"id":"9f5d4e77-572f-4fb9-b356-9c4b4cadb5c9","arxiv_id":"2506.15772","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Graviton bremsstrahlung from right-handed neutrino decays during leptogenesis yields a gravitational wave spectrum with peak frequency inversely proportional to the Yukawa coupling and amplitude proportional to the neutrino mass squared.","lead":"This paper calculates gravitational waves emitted as graviton bremsstrahlung when heavy right-handed neutrinos decay during leptogenesis, including a phase of early matter domination. It predicts a spectrum whose peak frequency depends on the Yukawa coupling and whose amplitude grows with the square of the neutrino mass, offering a possible future observational test of leptogenesis.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The early matter-dominated phase requires N1 to decouple before Boltzmann suppression; the paper assumes thermalization via gauge interactions but never computes the N1 yield, so the M^2-scaled spectrum may rest on an unestablished abundance.","rationale":"The reader's weakest assumption already pointed to thermalization and the survival of N1 until it dominates, which is the right area. I sharpen that concern: the same gauge interactions invoked to populate N1 can keep it in chemical equilibrium through T ~ M, causing Boltzmann suppression and preventing the required early matter domination. The paper's Eq. (3.5) is only a lifetime bound; it does not ensure a sufficient abundance. This is a genuine gap because the M^2 scaling and the Y^{-1} peak frequency are derived from the matter-dominated phase, not from the decay process itself. I do not flag the matrix-element disagreement with Refs. [19,20] as the primary concern because the paper provides internal consistency checks (soft graviton theorem, angular momentum conservation, gauge invariance), whereas the abundance issue has no corresponding check in the text. Nevertheless, the concern is addressable: a dedicated Boltzmann calculation for a specified gauge sector would settle it. Since the paper is already CONDITIONAL and this concern is a condition rather than a demonstrated contradiction, I leave the reader's verdict unchanged.","tokens_in":21333,"tokens_out":34592,"duration_ms":335625,"concrete_test":"For a minimal U(1)_B-L realization, solve the Boltzmann equation for the N1 yield with parameters varied over g in [0.01, 1] and M_Z' in [10^13, 10^17] GeV, using T_rh = 10^16 GeV, M = 10^15 GeV, and Y = 10^-4, and include N1 scattering, decays, and inverse decays. Compute rho_N(T_D)/rho_rad(T_D) and the time of matter-radiation equality. If no parameter point that satisfies thermalization at T_rh gives rho_N(T_D) >= rho_rad(T_D), then the early matter-dominated phase, and with it the M^2-scaled graviton bremsstrahlung spectrum, is not established for the claimed benchmark.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central observable, h^2 Omega_gw ~ 2.8e-16 (M/10^15 GeV)^2 in Eq. (4.6) and the peak frequency scaling in Eq. (4.5), depends on the universe passing through an early matter-dominated phase driven by right-handed neutrinos. The only condition stated for this phase, Eq. (3.5), Y^2 < 0.23 M/m_Pl, is a lifetime constraint: it ensures N1 decays after matter-radiation equality, but it says nothing about whether the N1 abundance is large enough to dominate. In Sec. 3 the authors assume that gauge interactions such as U(1)_B-L or SU(2)_R keep N1 in thermal equilibrium, but they never compute N1 decoupling or freeze-out. If those interactions maintain chemical equilibrium down to T ~ M or below, the N1 number density is Boltzmann-suppressed and rho_N never overtakes rho_SM; in that case there is no early matter domination, no M^2 scaling, and no CGMB degeneracy breaking. The danger is concrete: for a light Z' (M_Z' < T), the relativistic interaction rate is Gamma ~ g^4 T, so Gamma/H ~ g^4 m_Pl/T grows as the universe cools; a species thermalized at T_rh stays in equilibrium through T ~ M and is exponentially depleted. Early matter domination instead requires decoupling while T >~ M, or a residual yield large enough to dominate, and this depends on the gauge coupling and Z' mass, neither of which is specified or varied. Thus the cosmological premise of the headline signal is an assumption, not a derived property of the benchmark.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies stochastic gravitational wave backgrounds generated during seesaw leptogenesis. It computes graviton bremsstrahlung in the decays of the lightest right-handed neutrino, presents a detailed gauge-invariant matrix element calculation in Appendix A, and derives the resulting GW spectrum for a cosmological history that includes an early matter-dominated phase. The headline results are the scalings h^2 Omega_gw ~ 2.8 x 10^-16 (M/10^15 GeV)^2 and f_peak ~ 7.9 x 10^12 Hz (M/10^15 GeV)^{1/2} (10^-4/Y). The paper also re-evaluates the Cosmic Gravitational Microwave Background (CGMB) in the presence of early matter domination and argues that a joint observation of both spectra could break degeneracies among M, Y, T_rh, and g*. The central claim is that a high-frequency GW background from right-handed neutrino decays would constitute strong evidence for leptogenesis.","tokens_in":21645,"tokens_out":43997,"duration_ms":384745,"significance":"If the early matter-dominated phase is real, the paper provides a careful and testable prediction, with a matrix element that is checked against the soft graviton theorem and angular momentum conservation. The CGMB analysis raises an important degeneracy that is worth pointing out, and the paper is transparent about the benchmark assumptions. However, the observability is far below current and near-future sensitivities, and, more importantly, the existence of the early matter-dominated phase is assumed rather than derived from a concrete particle physics model. The significance is therefore moderate: the paper is a useful phenomenological target for high-frequency GW efforts, but its headline signal is contingent on an unestablished cosmological premise.","major_comments":[{"comment":"The condition Y^2 < 0.23 M/m_Pl is necessary for an early matter-dominated phase only if the right-handed neutrino number density at T = M is close to the thermal relativistic value and then freezes in as matter. This is not established. If the U(1)_B-L or SU(2)_R gauge interactions invoked in Sec. 2 keep N1 in chemical equilibrium below T ~ M, the equilibrium abundance is Boltzmann-suppressed and rho_N never overtakes rho_SM; if N1 decouples earlier, the yield depends on the gauge coupling and gauge boson mass, neither of which is specified or varied. No Boltzmann equation for n_N, including annihilations and inverse decays, is presented. Because the M^2 scaling in Eq. (4.6) and the peak-frequency scaling in Eq. (4.5) are contingent on the early matter-dominated phase, the paper should either compute the N1 yield/decoupling for a concrete gauge model or state the thermal abundance as an explicit assumption and estimate the gauge couplings for which it holds.","section":"Sec. 3, Eq. (3.5)"},{"comment":"The alpha parameter does not scale as stated. Combining Eq. (4.10) with Eq. (3.8) gives T_D proportional to M^{1/2} Y^{1/4}, whereas Eq. (2.10) gives T_D proportional to M^{1/2} Y; for M = 10^15 GeV and Y = 10^-4 the two expressions differ by about a factor 1.5 in T_D, and the discrepancy grows as Y is lowered. A direct estimate from t_D = 8 pi/(Y^2 M) and t_MD = t_M (g*_SM/g*_N)^2 gives alpha proportional to (M/m_Pl)^{2/3} Y^{-4/3}, not Y^{-1/3}. Consequently the CGMB peak scalings in Eqs. (4.12)-(4.13) and the alpha values used in Fig. 6 are not correct as written and should be revised before the degeneracy-breaking claim is assessed.","section":"Sec. 4.2, Eq. (4.10)"}],"minor_comments":[{"comment":"The figures use T_rh = 10^16 GeV, which exceeds the instantaneous-reheating upper bound T_rh,max = 5.7 x 10^15 GeV derived in Eq. (2.7). A more optimistic inflationary model should be cited, or the plots should use the stated bound.","section":"Figs. 6 and 7"},{"comment":"The discussion of dropping the n_B(r) and n_F(q) factors is clear, but it would be helpful to state explicitly that the suppression for soft lepton or Higgs momenta is controlled by the matrix element, as the phase space alone would not justify neglecting these terms.","section":"Sec. 4.1, Eq. (4.3)"},{"comment":"There are several typographical and readability issues, including missing parentheses in Eq. (2.13) and inconsistent formatting of g* and h^2 Omega_gw; a careful proofread is recommended.","section":"General presentation"}],"recommendation":"major_revision","confidential_remarks":"The matrix-element calculation in Appendix A is a genuine strength: the gauge-invariance check, the soft-graviton theorem, and the angular-momentum arguments give independent support for the claimed disagreement with Refs. [19,20]. The main risk is not the graviton calculation but the cosmological premise: the early matter-dominated phase is assumed without a freeze-out calculation. The alpha-scaling error in Sec. 4.2 is concrete and fixable, but it currently affects the quantitative CGMB predictions. In my view this is a major revision rather than a rejection, provided the authors either derive the N1 abundance in a specified gauge model or clearly frame the early-MD benchmark as a conditional assumption."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing you should know first: this is a solid correction-and-extension paper, not a detection paper. The new pieces that matter are the re-derived graviton bremsstrahlung matrix element for N1 -> l H h, with consistency checks the earlier work on this process lacked, and the CGMB degeneracy argument: an early matter-dominated era shifts the thermal graviton background in a way that mimics extra relativistic species and a higher reheating temperature, and the bremsstrahlung spectrum can break that degeneracy. Both are real contributions.\n\nThe calculation deserves credit. The appendix computes the spin-averaged amplitude with explicit gauge invariance, recovers Weinberg's soft graviton theorem in the x_G -> 0 limit, and explains the (1 - x_G) vs (1 - x_G)^2 disagreement with Refs [19,20] via angular momentum conservation in the collinear limit. I could not independently verify the overall phase-space normalization in Eq. (A.45) on a quick read, but the checks are the right ones and the argument against the older results is specific rather than hand-wavy. That is worth a referee's time.\n\nThe soft spots are in proportion. First, the signal sits four to six orders of magnitude below any planned sensitivity, and the authors say this plainly; the paper is a roadmap for future high-frequency detectors, not a current observable. Second, and more concerning, the early matter domination that gives the headline M^2 scaling is assumed rather than derived. Section 3 assumes N1 stays in thermal equilibrium through gauge interactions like U(1)_B-L or SU(2)_R, but never computes the N1 yield or decoupling temperature. If those interactions keep N1 chemically coupled down to T ~ M, the abundance is Boltzmann-suppressed and there is no matter-dominated era, no M^2 spectrum, and no degeneracy breaking. The stress-test note is correct on this point. The fix is straightforward in principle — pick a gauge sector, compute decoupling — but it is load-bearing. Third, 'observing leptogenesis in action' overstates the uniqueness: any heavy decaying particle with similar mass and Yukawa coupling would produce a similar spectrum. The authors explicitly acknowledge that, so it is a framing issue, not a hidden flaw.\n\nThis is a paper for the high-frequency GW and early-universe cosmology crowd; experimentalists should read it as a target-setting exercise. Bottom line: it deserves a serious referee. The matrix-element correction and the CGMB degeneracy result are worth having, and the missing abundance calculation is fixable rather than fatal. I would cite the CGMB part, and I would tell the authors to add the N1 yield calculation before publication.","headline":"A careful matrix-element correction and a useful CGMB degeneracy warning, but the headline spectrum rests on an uncomputed N1 abundance and is far from detection.","tokens_in":22218,"tokens_out":4559,"would_cite":true,"duration_ms":47865,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Heavy right-handed neutrino decays during leptogenesis are argued to emit a stochastic gravitational-wave background whose amplitude scales with neutrino mass squared and peak frequency inversely with Yukawa coupling, offering a direct…","keywords":["leptogenesis","right-handed neutrinos","seesaw mechanism","graviton bremsstrahlung","stochastic gravitational-wave background","early matter domination","cosmic gravitational microwave background","high-frequency gravitational waves"],"falsifier":"An independent recomputation of the graviton-bremsstrahlung phase-space integral would settle a decisive point: if the squared amplitude vanishes as $(1-x_G)^2$ rather than $(1-x_G)$ near the kinematic endpoint, the predicted peak amplitude and frequency would change. Observationally, a future CMB spectral-distortion survey reaching a few nK sensitivity around the predicted $10^{13}$ Hz peak for $M=10^{15}$ GeV and $Y=10^{-4}$ would either find the peak or place a direct exclusion on this leptogenesis benchmark.","tokens_in":21096,"feed_emoji":"📡","tokens_out":8501,"duration_ms":84321,"temperature":0.7,"pith_summary":"This paper argues that the decays of the heavy right-handed neutrinos responsible for leptogenesis also emit gravitons, producing a stochastic gravitational-wave background that is a direct fossil of the decay process. If those neutrinos drove an early matter-dominated phase, the spectrum's peak amplitude scales as the square of the neutrino mass and its peak frequency scales inversely with the neutrino Yukawa coupling. The paper also shows that the same early matter era reshapes the gravitational-wave background from the thermal plasma, and that combining the two spectra can disentangle the reheating temperature, the effective number of relativistic degrees of freedom, the right-handed neutrino mass, and its Yukawa coupling. For a mass $M = 10^{15}$ GeV and coupling $Y = 10^{-4}$, the predicted peak sits near $10^{13}$ Hz with $h^2\\Omega_{\\rm gw}\\sim 10^{-16}$, far above current detectors but within the conceptual reach of proposed high-frequency and CMB spectral-distortion techniques.","feed_headline":"Gravitational waves could catch leptogenesis in the act","feed_subtitle":"A predicted high-frequency background from right-handed neutrino decays would expose the particles that seeded the baryon asymmetry.","key_machinery":"The load-bearing object is the spin- and polarization-averaged squared amplitude for graviton bremsstrahlung in $N_1\\to \\ell H h$, computed from four Feynman diagrams and written as $|{\\cal M}|^2_{\\rm av} = \\tfrac12 |Y|^2(\\kappa/8)^2\\,16M^2(1-x_G)(2-x_G-x_Gx_L)/x_G^2$, where $x_G=2k/M$ and $x_L=2q/M$ are the graviton and lepton energy fractions. The paper checks this amplitude in two limits: as $x_G\\to 0$ it reproduces the soft-graviton theorem, and as $x_G\\to 1$ the amplitude vanishes only linearly with $(1-x_G)$, as angular-momentum conservation demands, in contrast to earlier results that vanish quadratically. The phase-space integral then gives the spectral shape $x(x-2)^2(1-x)$, and an instantaneous-decay approximation yields the peak frequency and amplitude. For the thermal plasma contribution, the key object is the gain term $\\hat\\eta(T,k/T)$ from the standard CGMB calculation, modified by the scale-factor ratio $\\alpha=a_D/a_{MD}$ that encodes how long the early matter-dominated phase lasted.","core_discovery":"The paper's central claim is that graviton bremsstrahlung in the decay $N_1\\to \\ell H h$ (where $h$ is the graviton) during leptogenesis gives a stochastic gravitational-wave background whose amplitude and peak frequency carry the mass $M$ and Yukawa coupling $Y$ of the decaying right-handed neutrino. Including the early matter-dominated era that a long-lived $N_1$ induces, the spectrum peaks at $f_{\\rm peak}\\simeq 7.9\\times 10^{12}\\,{\\rm Hz}\\,(M/10^{15}\\,{\\rm GeV})^{1/2}(10^{-4}/Y)$ with $h^2\\Omega_{\\rm gw}(f_{\\rm peak})\\simeq 2.8\\times 10^{-16}(M/10^{15}\\,{\\rm GeV})^2$. The calculation uses a newly derived spin-and-polarization averaged matrix element checked against the soft-graviton theorem and angular-momentum conservation, and it disagrees with earlier attempts at this rate. The paper further computes how an early matter-dominated era modifies the Cosmic Gravitational Microwave Background, finding that it can mimic extra relativistic degrees of freedom and a higher reheating temperature; the bremsstrahlung spectrum breaks that degeneracy. The authors do not claim this alone proves leptogenesis, but they argue that detection would make the case far more convincing when combined with other circumstantial evidence.","pith_inferences":["Beyond the paper, if the bremsstrahlung background were observed together with the cosmic-string signal expected from $U(1)_{B-L}$ breaking, the two backgrounds would pin down both the symmetry-breaking scale and the decay parameters, effectively reconstructing the full early-universe history around leptogenesis.","Beyond the paper, an independent recomputation of the phase-space integral could settle the discrepancy with earlier decay-rate calculations; the two candidate scalings predict different peak amplitudes, so the check is a closed-form calculation that does not require new data.","Beyond the paper, a natural extension is to relax the single-generation assumption and include washout and flavor effects; the gravitational-wave signal would then be correlated with the surviving lepton asymmetry, turning the spectrum into a quantitative consistency check on the leptogenesis mechanism itself.","Beyond the paper, the spectrum's shape distinguishes production scenarios: if the right-handed neutrinos are not thermally populated, the early matter-dominated phase is absent and both the peak amplitude and frequency change, so an observed or absent peak can discriminate among reheating and production models."],"forward_implications":["If the predicted spectrum exists, a single detection would be direct evidence that a heavy particle decayed in the early universe at the right time and with the right couplings to be the source of the observed baryon asymmetry.","From the observed peak frequency and amplitude one can read off the right-handed neutrino mass $M$ and Yukawa coupling $Y$; combined with the CGMB spectrum, that information fixes the reheating temperature and effective degrees of freedom and removes the degeneracy the CGMB alone has.","The early matter-dominated phase makes the bremsstrahlung peak amplitude scale as $M^2$ and its frequency shift with $1/Y$, a distinctive pattern that separates this source from cosmic strings, inflation, and other gravitational-wave backgrounds.","Smaller $Y$ moves the bremsstrahlung peak to higher frequencies while moving the CGMB peak to lower frequencies, making the two signals easier to separate observationally.","The predicted frequencies, in the GHz-to-THz range, lie beyond current broadband sensitivity, so the result gives a concrete target for future high-frequency gravitational-wave detectors and CMB spectral-distortion experiments."],"supporting_citations":[{"why":"Supplies the gain-term calculation for gravitational-wave production from the thermal plasma that the CGMB section builds on.","marker":"[16]"},{"why":"Completes the leading-order Standard Model CGMB calculation used as the baseline.","marker":"[17]"},{"why":"Introduces the $\\hat\\eta$ parameterization and the dark-radiation bound used to evaluate the CGMB spectrum.","marker":"[18]"},{"why":"The earlier graviton-bremsstrahlung rate calculation whose endpoint behavior this paper argues is incorrect.","marker":"[19]"},{"why":"The other prior rate calculation, which this paper argues copied a scalar-decay expression inappropriate for neutrinos.","marker":"[20]"},{"why":"The soft-graviton theorem used to validate the new matrix element at low graviton energies.","marker":"[45]"}],"fun_headline_variants":["Gravitational waves could expose the neutrino that seeded matter","Right-handed neutrino decays leave a gravitational wave signature","Gravitational wave background could fingerprint leptogenesis's agent","Heavy neutrino decays leave a specific graviton bremsstrahlung signal"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument requires that the right-handed neutrinos actually reach thermal equilibrium in the early universe and that the lightest one, with a small enough Yukawa coupling (about $10^{-2}$ or below) and satisfying $Y^2 < 0.23\\,M/m_{\\rm Pl}$, lives long enough to dominate the energy density before decaying; if thermalization fails or another particle takes over the universe, the early matter-dominated phase and the predicted $M^2$ scaling do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Gravitational waves could expose the neutrino that seeded matter","Right-handed neutrino decays leave a gravitational wave signature","Gravitational wave background could fingerprint leptogenesis's agent","Heavy neutrino decays leave a specific graviton bremsstrahlung signal"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.002102,"raw_usage":{"total_tokens":8190,"prompt_tokens":992,"completion_tokens":7198,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":608,"completion_tokens_details":{"reasoning_tokens":7129}},"tokens_in":608,"tokens_out":7198,"duration_ms":48345,"temperature":1.0,"reasoning_tokens":7129,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:32:42.923091+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An independent recomputation of the graviton-bremsstrahlung phase-space integral would settle a decisive point: if the squared amplitude vanishes as $(1-x_G)^2$ rather than $(1-x_G)$ near the kinematic endpoint, the predicted peak amplitude and frequency would change. Observationally, a future CMB spectral-distortion survey reaching a few nK sensitivity around the predicted $10^{13}$ Hz peak for $M=10^{15}$ GeV and $Y=10^{-4}$ would either find the peak or place a direct exclusion on this leptogenesis benchmark.","supporting_citations":[],"review_version":2}