{"id":"91c49051-a114-49dd-8dc7-75a98427cb97","arxiv_id":"2502.07988","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In a two-heavy-neutral-lepton seesaw with large imaginary Casas-Ibarra angle, the cLFV matrices S and R become rank-one related, and the resulting baryon asymmetry vanishes at leading order.","lead":"Two heavy sterile neutrinos can make the low-energy matrices that describe charged lepton flavour violation satisfy a simple rank-one relation, even when the sterile neutrino masses differ. The paper also ties the baryon asymmetry of the Universe to these same matrices, though the resulting bounds turn out to be far too weak to constrain experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Improved bounds (9) follow from positivity alone; the new saturation equalities (8) are not load-bearing for the paper's headline numerical claim.","rationale":"The paper's central new algebraic result is Eq. (8), the saturation equalities for S and R in the large-Imω limit with non-degenerate heavy masses, and the advertised consequence is the stronger bounds in Eq. (9). I checked whether Eq. (9) actually requires Eq. (8). Table 1 gives Ŝ_ee + Ŝ_μμ ≤ 0.53·10^{-3} and Ŝ_ττ ≤ 0.64·10^{-3}. Since S is a sum of Hermitian positive semidefinite rank-one terms, it is positive semidefinite, and Cauchy-Schwarz gives |Ŝ_αβ|^2 ≤ Ŝ_αα Ŝ_ββ. Combining these immediately yields |Ŝ_eτ| ≤ sqrt(0.53·0.64)·10^{-3} ≈ 0.58·10^{-3}, and likewise for μτ. Thus the headline bound (9) is an elementary consequence of the input constraints and does not depend on assumption (4), non-degeneracy, or saturation. If the bounds in (9) are the paper's claim to 'strengthen the existing cLFV constraints', that strengthening is not a validation of the new relations. The reader's formal weakest assumption, the accuracy of (4), is therefore not the most load-bearing issue: even if (4) is violated the numerical bounds stand, while the conceptual claim about saturation loses its advertised payoff. The paper should either reframe (9) as a trivial positivity bound or demonstrate a consequence of (8) that is not already implied by positive semidefiniteness. Equation (10) is also a black box from Ref. [11] with undefined φ and M0 and should be checked, but the overclaim about (9) is the more directly load-bearing concern. Since this is a framing and overclaim issue rather than a demonstrated algebraic error, the reader's conditional verdict remains appropriate.","tokens_in":6144,"tokens_out":13649,"duration_ms":119016,"concrete_test":"Recompute the right-hand sides of Eq. (9) using only the Schwarz inequality |S_αβ|^2 ≤ S_αα S_ββ for the positive semidefinite matrix S = F M^{-2} F† and the two diagonal entries in Table 1 (Ŝ_ee + Ŝ_μμ ≤ 0.53·10^{-3}, Ŝ_ττ ≤ 0.64·10^{-3}), omitting Eq. (8) and assumption (4) entirely. If the resulting bounds are exactly |Ŝ_eτ|, |Ŝ_μτ| ≤ 0.58·10^{-3}, then the saturation equalities are not load-bearing for the paper's headline strengthening. If the result is weaker, the paper's presentation is justified and the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper presents Eq. (9) as a consequence of the new saturation equalities (8) under assumption (4). This is the advertised payoff of the paper. However, Table 1 already contains the diagonal bounds Ŝ_ee + Ŝ_μμ ≤ 0.53·10^{-3} and Ŝ_ττ ≤ 0.64·10^{-3}, and S = F M^{-2} F† is positive semidefinite. Cauchy-Schwarz therefore gives |Ŝ_eτ|^2 ≤ Ŝ_ee Ŝ_ττ ≤ (Ŝ_ee + Ŝ_μμ) Ŝ_ττ ≤ 0.53·10^{-3} · 0.64·10^{-3}, so |Ŝ_eτ| ≤ 0.58·10^{-3}. The identical argument yields |Ŝ_μτ| ≤ 0.58·10^{-3}. No use of Eq. (8), of non-degeneracy, or of assumption (4) is required. The 'strengthened constraints' (9) are therefore a direct corollary of the same input bounds and do not test the new saturation relations. This does not show that Eq. (8) is false, but it removes the principal motivation: the numerical improvement that Eq. (8) is said to enable is already implied by the table the paper starts from.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the νMSM with two heavy neutral leptons and derives approximate algebraic forms for the effective cLFV matrices Sαβ and Rαβ under the large-Imω condition (4). The authors obtain the proportionality relation (7), the saturation equalities (8), and then claim improved off-diagonal bounds (9). In the second part, they import a baryon asymmetry expression from Ref. [11], translate it into an upper bound (10) on nB/s in terms of Ŝ entries, evaluate it with Table 1, and conclude that baryogenesis and accelerator constraints differ by many orders of magnitude. The paper contains no free parameters; the inputs are external experimental limits and standard seesaw relations.","tokens_in":6353,"tokens_out":19250,"duration_ms":172572,"significance":"The algebraic observation that non-degenerate heavy neutrino masses still produce a rank-one form for S and R in the large-Imω limit is potentially useful and would extend earlier degenerate/massless results. The derivation is not shown, however, and the advertised numerical improvement (9) is not actually powered by the new relations. The baryon-asymmetry part does not yield the claimed restrictions and contains internal logical inconsistencies. With the missing derivation supplied and the overclaims corrected, the paper could be a modest but acceptable contribution; in its present form its significance is limited.","major_comments":[{"comment":"The bounds |Ŝeτ| ≤ 0.58×10^{-3} and |Ŝμτ| ≤ 0.58×10^{-3} do not test or use the new saturation equalities (8). They follow directly from the diagonal bounds already quoted in Table 1 together with positive semidefiniteness of Ŝ: |Ŝeτ|² ≤ ŜeeŜττ ≤ (Ŝee+Ŝμμ)Ŝττ ≤ 0.53×10^{-3} × 0.64×10^{-3}, and likewise |Ŝμτ|² ≤ ŜμμŜττ ≤ (Ŝee+Ŝμμ)Ŝττ. Thus Eq. (9) is a corollary of the input bounds and of the ordinary Cauchy–Schwarz inequality, not of the rank-one limit (4). The paper should either demonstrate a quantitative consequence of Eq. (8) that goes beyond positivity or explicitly present Eq. (9) as the standard positivity bound and reframe the novelty claim.","section":"§2, Eq. (9) and Table 1"},{"comment":"The central derivation is not shown. The paper does not give the explicit Casas-Ibarra R matrix used, nor the asymptotic expansion that leads to the factorized forms (5)–(6), nor an estimate of the error made when assumption (4) is imposed. Equation (8) is an exact equality only in the limit e^{2 Imω} → ∞; for finite Imω the ratio |Sαβ|²/(SααSββ) differs from unity by terms of relative order e^{-2 Imω}, with coefficients that depend on the other parameters and on the mass ratio. Since the experimental region 'above the seesaw line' may not satisfy (4) with the required accuracy, the practical validity of Eq. (8) is not established. A short appendix giving the R-matrix parameterization, the expansion, and a numerical example of the deviation would be necessary to support the claim that Eq. (8) 'holds true with sufficient accuracy'.","section":"§2, Eqs. (5)–(8)"},{"comment":"Equation (10) is quoted from Ref. [11] with no derivation and with undefined symbols (φ, M0, TW, and ΔM21 in Eqs. (11)–(12)). The accompanying statement that the baryon asymmetry is proportional to Im[S∗αβRαβ] needs clarification: since S and R defined in Eqs. (2)–(3) are Hermitian matrices, the sum over α≠β of Im[S∗αβRαβ] vanishes identically; if the relevant CP-violating combination is a fixed off-diagonal term or a different contraction, this must be stated and derived. The text also says that Eq. (10) is obtained 'without assumption (4)', although the preceding derivation of Eq. (7), which is the only visible link between S and R, relies on (4); this tension needs to be resolved.","section":"§3, Eq. (10)"},{"comment":"The estimates (11)–(12) are upper bounds on nB/s whose right-hand sides are much larger than the observed value for M ≳ 1 GeV and M/ΔM21 ≫ 1. An upper bound that is far above the observed asymmetry imposes no restriction on Ŝ or R̂, and it does not imply that 'the actual values of the observed elements of the Ŝ and R matrices are much lower than the experimental limits'. The Conclusions similarly refer to 'lower limits (baryon asymmetry)', but no lower limit on the seesaw parameters is derived in the paper. This part should be rewritten to state only what the inequality actually shows, and the abstract's claim of 'new restrictions' should be withdrawn or replaced by a correct statement.","section":"§3 and Conclusions"}],"minor_comments":[{"comment":"The approximation in Eq. (4) implicitly requires Im ω > 0; for Im ω < 0, sinh 2 Im ω is negative and the stated chain of approximations fails. Please state the sign condition explicitly.","section":"§2, Eq. (4)"},{"comment":"The 'Future experiments' column is empty for Ŝee+Ŝμμ and Ŝττ; please clarify whether no future projections exist or whether the entries are simply omitted.","section":"Table 1"},{"comment":"The notation '(M/1GeV)' should be typeset as '(M/GeV)' or 'M/(1 GeV)', and 'ΔM21' should be defined explicitly, presumably as M2−M1.","section":"§3, Eqs. (11)–(12)"},{"comment":"The phrase 'as one can effortlessly see' is informal; the reality and positivity of the diagonal elements follow from the displayed factorized form and should simply be stated.","section":"§2, before Eq. (5)"},{"comment":"The statement that Eq. (8) is 'independent of the mass difference between sterile neutrinos' should be qualified as holding only to leading order in the asymptotic expansion (4), not as an exact statement for arbitrary heavy neutrino masses.","section":"Conclusions"}],"recommendation":"major_revision","confidential_remarks":"The referee's main concern is that the paper's headline quantitative claims are either already implied by existing bounds or not derived. Requiring an appendix with the R-matrix derivation and a rewrite of the baryon section should be a condition of acceptance. The paper is short and could be revised without changing its scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper has a genuinely new algebraic observation—the S-R linear relation (7) and the saturation equalities (8) extended to non-degenerate HNL masses with massive active neutrinos. But the advertised payoff, the improved bounds (9), does not actually test that result: it follows from trivial positivity of S and the diagonal limits already in Table 1. The BAU part is a black box imported from [11] and, as the paper admits, is not even valid under the same assumption (4) that makes (8) hold.\n\nThe good parts: extending the equality |S_αβ|^2 = S_αα S_ββ beyond the degenerate-mass/massless-active-neutrino cases is a useful contribution. If (4) holds, (7)-(8) give a compact description of the cLFV parameter space of the minimal seesaw. The paper is also honest in noting that the baryon asymmetry bound has limited practical value (right side >> 1) and that under (4) the asymmetry vanishes; the nonzero asymmetry requires going beyond the approximation. That is a welcome admission.\n\nThe soft spots are real. First, the derivation of (5)-(8) is not shown despite the 'effortlessly see' comment, and assumption (4) is left unquantified. Second, the central numerical claim (9) is an overreach: with Ŝ_ee + Ŝ_μμ ≤ 0.53·10^-3 and Ŝ_ττ ≤ 0.64·10^-3, positivity alone gives |Ŝ_eτ|^2 ≤ (Ŝ_ee + Ŝ_μμ) Ŝ_ττ, so |Ŝ_eτ| ≤ 0.58·10^-3; no saturation equalities needed. Same for μτ. Thus (9) is a corollary of the table, not a test of (8). That should be stated clearly. Third, the baryon asymmetry expression (10) uses symbols (sin φ, M0, TW?) that are never defined and is simply quoted from [11]. Since the whole BAU link is built on it, the paper needs at least a derivation or a precise reference with definitions.\n\nFourth, and maybe most interesting: the regime where (8) saturates (Im ω large) is exactly where baryon asymmetry vanishes, whereas the BAU formula requires the deviation from (4). The paper notes this but does not discuss what it means for combining the two constraints. That is a real physics point a referee should push on.\n\nVerdict: worth sending to peer review. The algebraic generalization is plausible and deserves a careful check. But the referee should insist on a step-by-step derivation, quantification of the approximation, and a rewriting that separates the positivity bound from the saturation bound.","headline":"The saturation equalities are a real generalization, but the claimed numerical improvement is just positivity and the BAU link is a black box.","tokens_in":6958,"tokens_out":5350,"would_cite":false,"duration_ms":35537,"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":"In a two-heavy-neutrino seesaw with unequal masses, the flavour-violating couplings obey $|S_{\\alpha\\beta}|^2 = S_{\\alpha\\alpha} S_{\\beta\\beta}$ for large complex mixing angle, tightening $\\tau$–$e$ and $\\tau$–$\\mu$ bounds tenfold and…","keywords":["heavy neutral leptons","charged lepton flavour violation","seesaw mechanism","Casas-Ibarra parameterization","neutrino Minimal Standard Model","baryon asymmetry","leptogenesis","cLFV operators"],"falsifier":"Evaluate the exact ratio $|S_{\\alpha\\beta}|^2/(S_{\\alpha\\alpha}S_{\\beta\\beta})$ from equations (5)--(6) without imposing (4), scanning over $\\mathrm{Im}\\,\\omega$ and the heavy-mass splitting $\\Delta M/M$ with realistic active-neutrino masses: if the ratio stays well below 1 for the values of $\\mathrm{Im}\\,\\omega$ that still allow the observed baryon asymmetry, the improved bounds (9) do not hold in the baryogenesis-viable region, and a future measurement of $\\tau\\to\\mu\\gamma$ finding $|\\hat{S}_{\\mu\\tau}| > 0.58\\times10^{-3}$ would directly contradict the saturation bound.","tokens_in":5919,"feed_emoji":"⚛️","tokens_out":13377,"duration_ms":109993,"temperature":0.7,"pith_summary":"The paper asks how tightly experiments can bound the flavour-violating couplings that heavy sterile neutrinos induce, and what those couplings imply for the origin of matter. It claims that in a seesaw extension with two heavy neutral leptons of different masses and massive active neutrinos, the effective couplings $S_{\\alpha\\beta}$ and $R_{\\alpha\\beta}$ obey $|S_{\\alpha\\beta}|^2 = S_{\\alpha\\alpha} S_{\\beta\\beta}$ (and likewise for $R$) whenever the Casas-Ibarra mixing angle has a large imaginary part. That saturation — previously known only for massless active neutrinos or degenerate heavy masses — strengthens the experimental upper limits on $|\\hat{S}_{e\\tau}|$ and $|\\hat{S}_{\\mu\\tau}|$ to $0.58\\times10^{-3}$. The same operator language turns the $\\nu$MSM baryon asymmetry into a bound on the $S$ and $R$ entries, and shows that the saturation regime makes the asymmetry vanish. The result matters because it tells experimenters where to look and exposes a tension between the parameter region that gives the stronger collider bounds and the region that can generate the observed matter–antimatter asymmetry.","feed_headline":"Two sterile neutrinos tighten lepton-flavour bounds tenfold","feed_subtitle":"The tau–electron and tau–muon couplings drop below 0.58×10⁻³, and baryon asymmetry is set by the same operators.","key_machinery":"The load-bearing object is the Casas-Ibarra parameterization of the Yukawa matrix, packaged here as $X = (i/v)\\,U_\\nu\\sqrt{m^{\\rm diag}_\\nu}$. Under the large-$\\mathrm{Im}\\,\\omega$ assumption (4), the sums over the two heavy neutrinos in (2)--(3) collapse into a single rank-one factor $(X_{\\alpha 2}-iX_{\\alpha 3})(X^*_{\\beta 2}+iX^*_{\\beta 3})$ times a mass-dependent coefficient. This factorization is what turns the Schwarz inequality into the saturation equality for both $S_{\\alpha\\beta}$ and $R_{\\alpha\\beta}$, and it is also why $\\mathrm{Im}[S^*_{\\alpha\\beta}R_{\\alpha\\beta}]$, the combination controlling baryon asymmetry, vanishes identically in the same limit.","core_discovery":"The central claim is that the seesaw observables $S_{\\alpha\\beta}$ and $R_{\\alpha\\beta}$ — the effective couplings generating charged lepton flavour violation — saturate the Schwarz inequality, $$|S_{\\$\\alpha$\\$\\beta$}|^2 = S_{\\$\\alpha$\\$\\alpha$} S_{\\$\\beta$\\$\\beta$}, \\qquad |R_{\\$\\alpha$\\$\\beta$}|^2 = R_{\\$\\alpha$\\$\\alpha$} R_{\\$\\beta$\\$\\beta$},$$ not only in the previously studied limits (massless active neutrinos and degenerate heavy masses) but also when the active neutrinos are massive and the two heavy neutral leptons have different masses, provided the complex Casas-Ibarra angle $\\omega$ obeys (4), i.e. $\\cosh 2\\,\\mathrm{Im}\\,\\omega \\simeq \\sinh 2\\,\\mathrm{Im}\\,\\omega \\simeq e^{2\\,\\mathrm{Im}\\,\\omega}/2 \\gg 1$. From these equalities the paper derives the improved bounds $|\\hat{S}_{e\\tau}| \\le 0.58\\times10^{-3}$ and $|\\hat{S}_{\\mu\\tau}| \\le 0.58\\times10^{-3}$. Re-expressing the $\\nu$MSM baryon asymmetry in the same language yields Eq. (10), and shows that under (4) the asymmetry, being proportional to $\\mathrm{Im}[S^*_{\\alpha\\beta}R_{\\alpha\\beta}]$, vanishes, so generating the observed asymmetry requires leaving the saturation regime.","pith_inferences":["If the same large-parameter limit is realized with more than two heavy neutrinos, an analogous rank-one factorization should appear, so the saturation relations and the resulting bound improvement are likely not special to the two-flavour case.","Since the saturation regime forces $\\mathrm{Im}[S^*_{\\alpha\\beta}R_{\\alpha\\beta}]$ to vanish, the goals of maximal cLFV sensitivity and successful baryogenesis pull toward opposite corners of parameter space; a positive cLFV signal near the improved bounds would count against the $\\nu$MSM leptogenesis explanation of the baryon asymmetry.","A numerical scan of $|S_{\\alpha\\beta}|^2/(S_{\\alpha\\alpha}S_{\\beta\\beta})$ as a function of $\\mathrm{Im}\\,\\omega$ and $\\Delta M/M$, with finite active-neutrino masses, would turn the qualitative condition (4) into a quantitative statement of how large the large-parameter limit must be."],"forward_implications":["The current upper limits on the tau–electron and tau–muon seesaw couplings tighten from the $10^{-3}$–$10^{-2}$ level in Table 1 to $|\\hat{S}_{e\\tau}|, |\\hat{S}_{\\mu\\tau}| \\le 0.58\\times10^{-3}$.","A future observation of charged lepton flavour violation in any one channel would, within this model, fix the values in the other channels through the saturation equalities.","The baryon asymmetry formula (10) expresses a cosmological observable in terms of the same $S$ and $R$ operators that control lepton flavour violation, making collider and cosmology constraints comparable.","In the large-$\\mathrm{Im}\\,\\omega$ regime that produces the improved bounds, the baryon asymmetry vanishes, so the parameter region able to explain the matter–antimatter asymmetry is disjoint from the region with the strongest collider constraints.","Because the baryogenesis lower bound and the experimental upper bounds differ by many orders of magnitude, successful $\\nu$MSM baryogenesis implies the actual $S$ and $R$ entries lie far below the limits in Table 1."],"supporting_citations":[{"why":"Provides the Casas-Ibarra parameterization with a complex angle used to express all S and R elements in terms of neutrino masses and mixings.","marker":"[24]"},{"why":"Gives the effective-operator framework, the previous saturation result in limiting cases, and the experimental upper bounds in Table 1 that the paper improves.","marker":"[25]"},{"why":"Supplies the existing bounds on lepton non-unitarity and heavy-neutrino mixing that the new saturation relations generalize.","marker":"[26]"},{"why":"Provides the nuMSM baryogenesis calculation from which Eq. (10) for the baryon asymmetry in terms of S and R is taken.","marker":"[11]"},{"why":"Shows that leptogenesis can work without degenerate heavy-neutrino masses, motivating the non-degenerate case treated here.","marker":"[23]"}],"fun_headline_variants":["Seesaw saturation yields tenfold tighter lepton-flavour bounds","Lepton-flavour bounds and baryon asymmetry tied to seesaw saturation","Non-degenerate sterile neutrinos tighten lepton bounds tenfold","Baryon asymmetry vanishes as lepton flavour operators saturate"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the imaginary part of the complex neutrino mixing angle is very large, so that $\\cosh 2\\,\\mathrm{Im}\\,\\omega$, $\\sinh 2\\,\\mathrm{Im}\\,\\omega$ and $e^{2\\,\\mathrm{Im}\\,\\omega}/2$ are all interchangeable and much larger than one; if that limit is not satisfied, the saturation equalities break, and the parameter region that yields a non-zero baryon asymmetry is exactly where it breaks.","fun_headline_variants_meta":{"raw":{"variants":["Seesaw saturation yields tenfold tighter lepton-flavour bounds","Lepton-flavour bounds and baryon asymmetry tied to seesaw saturation","Non-degenerate sterile neutrinos tighten lepton bounds tenfold","Baryon asymmetry vanishes as lepton flavour operators saturate"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000833,"raw_usage":{"total_tokens":3637,"prompt_tokens":948,"completion_tokens":2689,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":2615}},"tokens_in":564,"tokens_out":2689,"duration_ms":17409,"temperature":1.0,"reasoning_tokens":2615,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T11:14:53.210588+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the exact ratio $|S_{\\alpha\\beta}|^2/(S_{\\alpha\\alpha}S_{\\beta\\beta})$ from equations (5)--(6) without imposing (4), scanning over $\\mathrm{Im}\\,\\omega$ and the heavy-mass splitting $\\Delta M/M$ with realistic active-neutrino masses: if the ratio stays well below 1 for the values of $\\mathrm{Im}\\,\\omega$ that still allow the observed baryon asymmetry, the improved bounds (9) do not hold in the baryogenesis-viable region, and a future measurement of $\\tau\\to\\mu\\gamma$ finding $|\\hat{S}_{\\mu\\tau}| > 0.58\\times10^{-3}$ would directly contradict the saturation bound.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the effective-operator framework, the previous saturation result in limiting cases, and the experimental upper bounds in Table 1 that the paper improves."},{"cited_title":"Drewes, B","cited_arxiv_id":null,"evidence_quote":"Shows that leptogenesis can work without degenerate heavy-neutrino masses, motivating the non-degenerate case treated here."}],"review_version":1}