{"id":"b4e7d7b4-8a76-4a09-a605-e004c3fe2bef","arxiv_id":"2509.03905","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"In weak-washout heavy-particle decay leptogenesis, the lepton asymmetry is approximately independent of the decaying particle mass and the couplings of thermalized decay products, because the nonthermal distribution scales inversely with both.","lead":"This paper finds that in a common type of weak-washout leptogenesis, the final lepton asymmetry barely changes when the decaying particle mass or the couplings of thermalized decay products are varied. The result matters because it tells model builders that tuning these parameters to boost the baryon asymmetry generally fails, while a working benchmark opens a broad parameter plateau.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No final-asymmetry integral is computed; cancellation is inferred from large-z δf scaling, but the z^4-weighted integrand peaks near z=O(1), where Eq. (26) is not established for low-momentum modes.","rationale":"The paper's argument is structurally sound: after factoring m^4 from the CP source, the only m,y dependence in Eq. (1) is in Im[y'^2y^2], MPl/m, and δf; if δf ∝ m/|y|^2 over the integration region, cancellation is exact. The steady-state balance behind Eq. (26) is transparent, and the two-mode numerical solutions support the scaling at large z. The weakest point is that the claim is about the final asymmetry, not about δf or Sbar_CP separately; no z-integrated Yℓ is presented. Because the z^4 weight biases the integral to z=O(1), one must verify the scaling at exactly the z where the integrand peaks. The paper explicitly notes a possible breakdown at z≈1 for rϕ≪1 and large m, and the low-momentum contribution is dismissed by an expectation rather than a computation. A direct numerical integration over z and rϕ would settle this. This does not invalidate the mechanism; it means the strong 'general phenomenon' statement is not yet established, consistent with the reader's CONDITIONAL verdict.","tokens_in":10990,"tokens_out":11111,"duration_ms":103036,"concrete_test":"Evaluate Yℓ directly from Eq. (1) with Sbar_CP from Eqs. (29)–(30), solving Eq. (22) on a dense rϕ grid (or better, solving the full momentum-dependent Boltzmann equation without the linearization Eq. (11)) for mϕ = 10^2...10^5 GeV and |y| = 10^-5...10^-2. Plot Yℓ/[Im(y'^2y^2)] versus mϕ and |y|. If the spread exceeds a factor of 2, the claimed 'basically independent' behavior fails. Also weight the integrand to identify which z and rϕ dominate, to test whether the low-momentum breakdown of Eq. (26) is kinematically suppressed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires the scaling δfϕ ∝ (1/|y|^2)(mϕ/Mbar_Pl), Eq. (26), to hold over the z-range that dominates the integral in Eq. (1). The paper demonstrates the scaling only through ratios Ri/j and R~i/j in Figs. 2–3, which approach the asymptotic values at z≫1. The z^4 weight in Eq. (1), together with Boltzmann suppression, makes the integrand peak near z=O(1)–few; at z=1 the scaling is not fully established, especially for rϕ=0.1 with mϕ=10^5 GeV, where D in Eq. (23) is smaller than the z coth(ξ/2)/ξ term and the paper itself concedes the simple scaling may break down. The authors show Sbar_CP in Fig. 4 but never compute Yℓ = cmd Im[y'^2y^2] (MPl/m) ∫ z^4 Sbar_CP dz; a partially compensating source can still change the final asymmetry by O(1). The linearization Eq. (11) is also uncontrolled where δfϕ∼O(10), but the more direct gap is the missing integrated test of the claimed insensitivity.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the mass and Yukawa-coupling dependence of the final lepton asymmetry in a class of weak-washout leptogenesis scenarios from scalar decay. The setup has a coupling hierarchy: a small Yukawa coupling keeps one decay product out of equilibrium, while a larger Yukawa coupling |y| keeps other flavors thermalized. The authors write the final asymmetry as Y_l = c_md Im[y'^2 y^2](M_Pl/m)∫ dz z^4 Sbar_CP(z) (Eq. (1)). They derive a linearized Boltzmann equation for the chemical-potential perturbation δf_φ, Eq. (22), solve it numerically for two momentum modes r_φ=0.1,1 and for m_φ=10^2–10^5 GeV, |y|=10^−5–10^−2, and observe from ratios that δf_φ ∝ |y|^−2(m_φ/Mbar_Pl) at large z. Since this scaling cancels the explicit |y|^2 and 1/m factors in Eq. (1), they conclude that varying the decaying-particle mass and the thermalized Yukawa coupling does not noticeably change the final asymmetry. They further plot the explicit CP-violating source Sbar_CP from Ref. [21] and claim that this demonstrates the cancellation. Exceptions are noted for early termination of leptogenesis and for flavor effects.","tokens_in":11214,"tokens_out":6300,"duration_ms":68287,"significance":"If the result is established, it is a genuinely useful conceptual correction: in this specific weak-washout class, parameter directions that are commonly expected to enhance leptogenesis are approximately flat, so a single benchmark opens a broad mass/coupling region while tuning to boost the asymmetry is ineffective. The paper is largely parameter-free up to the model-dependent constant c_md, and the ratios in Figs. 2–3 give direct numerical evidence for an asymptotic scaling of δf_φ. The main weakness is quantitative: the headline claim is about the integrated asymmetry in Eq. (1), but the paper never computes that integral and instead infers the cancellation from large-z scaling and from unweighted plots of |Sbar_CP|. The phase-space implementation of the explicit source is also not fully specified. These gaps are fixable and do not, by themselves, invalidate the idea, but the present version does not yet prove the claimed insensitivity.","major_comments":[{"comment":"The central conclusion requires that the weighted integral ∫ z^4 Sbar_CP dz scales as m_φ/(|y|^2 Mbar_Pl) times a constant, so that the prefactor in Eq. (1) is flat. Fig. 4 shows |Sbar_CP| as a function of z, not the z^4-weighted integral. Because z^4 peaks at z=O(1)–few, the large-z scaling in Figs. 2–3 is not sufficient. Please compute and report the actual Y_l ratios from Eq. (1), e.g. Y_l(m=10^5)/Y_l(m=10^2) at fixed |y| and Y_l(|y|=10^-2)/Y_l(|y|=10^-5) at fixed m, including the Im[y'^2 y^2] and 1/m prefactors.","section":"§V, Eq. (1)"},{"comment":"The explicit source in Eq. (29) is a three-dimensional momentum integral over r_1,r_2,r_3, and F in Eq. (30) contains δf_φ evaluated at some momentum, but §IV solves Eq. (22) only at the two representative momenta r_φ=0.1 and 1. The text does not state how δf_φ is extended to the full momentum range required by Eq. (29), or whether Fig. 4 uses a single representative mode. Without this information Fig. 4 is not reproducible and cannot quantitatively support the cancellation claim. Please specify the interpolation/procedure or provide the full momentum-dependent solution.","section":"§V, Eq. (29)"},{"comment":"The ansatz (11) is a leading-order perturbation with a small chemical potential. Fig. 2 shows δf_φ values that reach O(10) or larger for m_φ=10^5 GeV and r_φ=0.1 near z~1. For such values the linearized expression for δf_φ and hence Eq. (22) are no longer controlled, precisely in a parameter region used to infer the scaling (26). Either solve the full momentum-dependent distribution without the linearization, or restrict the numerical evidence to the region where Eq. (11) is valid and show that the conclusion is unchanged.","section":"§III–IV, Eqs. (11), (22), Fig. 2"},{"comment":"The paper concedes that the simple scaling may break down near z=1 for r_φ≪1. This is exactly the region that the z^4 weight in Eq. (1) can make relevant. The statement that r_φ≪1 modes contribute little because the dominant phase-space integration occurs at E(p)≃T is plausible but not quantified. Since Eq. (26) is demonstrated only at z≫1, the gap between the asymptotic scaling and the z=O(1) peak of the integrand must be closed before the integrated asymmetry claim is established.","section":"§IV, final paragraph; §I"}],"minor_comments":[{"comment":"The notation |y| is used both as the summed squared matrix elements in Eq. (19) and as a single dimensionless number |y|=10^{-i} in Fig. 3 and Eq. (28). Please clarify how the flavor-dependent Im[y'^2 y^2] in Eq. (1) is related to the |y| used in the Boltzmann equation, since the cancellation argument treats them as the same parameter.","section":"§II–IV, Eqs. (19), (28)"},{"comment":"In the definitions of R_{i/j} and Rtilde_{i/j}, the reference point is not stated in the equation itself; e.g., R_{2/3} is δf(m=10^2)/δf(m=10^3). Adding this explicitly would improve readability.","section":"§IV, Eqs. (27)–(28)"},{"comment":"The statement that the first term in D is smaller than the second for 'the typical parameter space' is not quantified. Since the balance between these two terms controls where Eq. (26) holds, giving the boundary in (m_φ,|y|,r_φ) would be helpful.","section":"§III, Eq. (23)"},{"comment":"Minor typographical issues: 'y′2y2' is written without parentheses, and Eq. (1) uses Im[y'^2y^2] while the text speaks of quartic Yukawa couplings; standardizing the notation would reduce ambiguity.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is short, clearly written, and the idea is interesting and likely correct in some form. The most serious technical gap is the missing evaluation of the integrated asymmetry in Eq. (1), together with the unclear momentum-extension in Eq. (29). Both seem straightforward for the authors to address, given that they already have the solver and the explicit source. If the integrated ratios confirm the cancellation in the controlled region, the paper would be a solid contribution. I do not see grounds for rejection, but the present form does not establish the headline insensitivity claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper identifies a real cancellation mechanism: in weak-washout decay leptogenesis, the nonthermal distribution of the decaying scalar scales as roughly 1/y^2 and m, cancelling the explicit y^2 and 1/m factors in the CP-asymmetry integral. That observation is genuinely useful, and it runs against the usual boost-by-large-Yukawa / low-mass intuition. The paper is also honest: it lists the escape hatches (early termination, flavor effects, strong washout) rather than hiding them.\n\nWhat's not there is the direct computation of the final asymmetry. The authors show scaling of δf and of the CP source Sbar_CP in the large-z limit, but the integral in Eq. (1) has a z^4 weight and the integrand peaks at z around a few. The stress-test note names exactly this gap, and it holds up on reading: Fig. 2 shows that at z ~ 1 the simple scaling is not yet established, especially for r_phi = 0.1 with large masses. So the central claim—that varying m and y leaves Y_ell basically unchanged—remains an inference, not a demonstrated result.\n\nThe second soft spot is the linearized chemical-potential ansatz, Eq. (11), used in regimes where δf reaches O(10) or more (Fig. 2, top left). The small-μ expansion is not controlled there, so the numerical solution itself could be an artifact. The authors need either a full Boltzmann treatment or a demonstration that the asymmetry integral is dominated by the small-δf region.\n\nThird, the abstract's 'general phenomenon' overstates the case. The authors themselves list exceptions, so the strong wording should be tempered.\n\nOn citations: the CP source is taken from their earlier papers, which is acceptable here because the scaling argument is independent of the explicit source form. I don't see circularity.\n\nBottom line: this is a clear, useful paper for model builders in leptogenesis, and it deserves a serious referee. But the published version should include a direct integration of Eq. (1) over a mass-coupling grid and a controlled check of the linearization. Without that, it's an intriguing observation, not a proof.","headline":"Useful counterpoint to the usual leptogenesis enhancement lore, but the claimed insensitivity of the final asymmetry is inferred from large-z scaling rather than demonstrated by integrating Eq. (1).","tokens_in":11753,"tokens_out":4686,"would_cite":true,"duration_ms":46514,"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 broad class of heavy-particle-decay leptogenesis models, the final lepton asymmetry is nearly insensitive to the thermalized Yukawa coupling and the decaying particle mass, because the out-of-equilibrium distribution scales as y^-2 m a","keywords":["leptogenesis","baryon asymmetry","CP violation","heavy particle decay","nonthermal distribution","weak washout","Yukawa coupling","mass-coupling cancellation"],"falsifier":"Solve the full, nonlinear Boltzmann equation without the linearized ansatz of Eq. (11) over m_phi ∈ [10^2, 10^5] GeV and |y| ∈ [10^-5, 10^-2], and compute the final asymmetry Y_L. The cancellation predicts a flat plateau; observing a systematic slope, for example Y_L ∝ |y|^2/m at z ≈ 1 where δf reaches O(10) for m_phi = 10^5 GeV, would falsify the central claim.","tokens_in":1801,"feed_emoji":"⚛️","tokens_out":1888,"duration_ms":73097,"temperature":0.7,"pith_summary":"This paper tries to establish that a common assumption in leptogenesis model building is wrong: in the class of scenarios where a decaying heavy particle provides the out-of-equilibrium condition, the final lepton asymmetry does not respond to the size of the thermalized Yukawa couplings or the particle mass. The reason is a cancellation: a larger coupling (or a smaller mass) drags the decaying particle back toward thermal equilibrium, so the nonthermal distribution at the epoch that matters scales roughly as 1/|y|^2 times m, which neutralizes the explicit |y|^2 and 1/m factors in the standard asymmetry formula. If true, hitting one working benchmark point automatically opens a broad plateau of viable masses and couplings, but trying to boost the asymmetry by tuning these parameters will not help. The claim is deliberately general, built on a simplified but nontrivial Boltzmann evolution rather than a scan of one specific model.","feed_headline":"Mass and coupling tweaks fail to boost leptogenesis","feed_subtitle":"The final asymmetry stays put because the nonthermal distribution cancels the coupling and mass factors.","key_machinery":"The central object is the damping (friction) force D in the Boltzmann equation for the decaying particle's departure from equilibrium, which scales as D ∝ |y|^2 (Mbar_Pl/m). This force controls how δf_phi approaches the large-z scaling δf ∝ (1/|y|^2)(m/Mbar_Pl). That scaling, combined with the explicit prefactor Im(y'^2 y^2) M_Pl/m in the asymmetry formula, produces the mass-coupling cancellation on which the whole paper rests.","core_discovery":"The paper argues that in weak-washout leptogenesis from heavy scalar (or fermion) decay, combined with a coupling hierarchy where some decay products are thermalized and others stay out of equilibrium, the final lepton asymmetry is approximately invariant under rescaling of the thermalized Yukawa coupling y and the decaying particle mass m. The asymmetry formula contains an explicit factor Im(y'^2 y^2) M_Pl/m times an integral of the CP-violating source, which is proportional to the departure from thermal equilibrium δf_phi. The evolution of δf_phi is governed by a damping force proportional to |y|^2 Mbar_Pl/m, and numerically this drives δf_phi into the simple scaling δf ∝ (1/|y|^2)(m/Mbar_","pith_inferences":["If this flatness holds beyond the studied class, the baryon asymmetry essentially constrains a combination of the CP phases and the quartic Yukawa product, but not the individual mass or the thermalized coupling; model builders should look for constraints from collider production, gauge interactions, or flavor structure rather than from the abundance itself.","The same scaling argument suggests a testable nonthermal-level prediction: a momentum-resolved measurement or simulation of the decaying particle's deviation from equilibrium should trace the m/|y|^2 trajectory across parameter space at z ≈ 1, which could confirm or break the plateau.","Because the structure of the damping force is stated to be common to fermion singlet decay, the cancellation should survive in right-handed-neutrino-like scenarios, making low-scale leptogenesis less tunable but more forgiving in parameter space."],"forward_implications":["Realizing leptogenesis at one benchmark point in this class automatically opens a much broader region of viable masses and Yukawa couplings, because the final asymmetry stays on a plateau.","Varying the decaying particle mass or the thermalized Yukawa couplings will not, by itself, enhance the asymmetry; boosting leptogenesis by tuning these parameters will be challenging.","The usual mass effect and large-Yukawa-coupling enhancement mechanisms fail to operate in this class, contrary to naive expectations.","The cancellation is not exact: scenarios where leptogenesis terminates early at z < 1, or scenarios with significant flavor effects in the Yukawa matrix, can evade the flat behavior.","The conclusion is derived under the weak-washout approximation; the strong-washout regime requires the full integro-differential Boltzmann equation, where the dependence on mass and couplings may be more complicated."],"supporting_citations":[{"why":"Review that supplies the standard Boltzmann-equation treatment of leptogenesis used to write the general parameterization of Eq. (1).","marker":"[7]"},{"why":"Review that supplies the kinetic-equation framework and model-specific formulas for the lepton asymmetry that Eq. (1) is meant to capture.","marker":"[8]"},{"why":"Scalar decay leptogenesis scenario with Dirac neutrinos, one of the concrete realizations of the class analyzed here.","marker":"[11]"},{"why":"Introduces the coupling hierarchy between small and large Yukawa elements, the enhancement mechanism this paper shows is canceled.","marker":"[15]"},{"why":"Hierarchical Dirac neutrino model realizing the parameter enhancement mechanism that the present paper's flatness claim directly targets.","marker":"[16]"},{"why":"Provides the explicit CP-violating source Sbar_CP used in Section V to numerically confirm the δf scaling and cancellation.","marker":"[21]"},{"why":"Complete thermal leptogenesis calculation giving the model-dependent formulas that motivate the factorized form of Eq. (1).","marker":"[25]"},{"why":"Higgs-doublet decay leptogenesis with one-loop self-energy CP asymmetry, the prototype scalar-decay setup whose Boltzmann evolution is solved here.","marker":"[30]"}],"fun_headline_variants":["Leptogenesis ignores mass and coupling tweaks","Mass and coupling changes don't move lepton asymmetry","Leptogenesis immune to mass and coupling tuning","Weak washout leptogenesis resists parameter tuning","Heavy decay asymmetry cancels mass and coupling effects"],"cache_read_input_tokens":13440,"weakest_assumption_plain":"The whole flatness result rests on the claim that the lepton asymmetry is generated at the epoch where the approximate scaling δf ∝ m/(|y|^2 Mbar_Pl) holds while the linearized chemical-potential perturbation of Eq. (11) is still accurate; for the largest masses considered, the plotted δf reaches order 10, which strains the small-perturbation assumption exactly where the cancellation is being invoked.","fun_headline_variants_meta":{"raw":{"variants":["Leptogenesis ignores mass and coupling tweaks","Mass and coupling changes don't move lepton asymmetry","Leptogenesis immune to mass and coupling tuning","Weak washout leptogenesis resists parameter tuning","Heavy decay asymmetry cancels mass and coupling effects"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000285,"raw_usage":{"total_tokens":1484,"prompt_tokens":680,"completion_tokens":804,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":424,"completion_tokens_details":{"reasoning_tokens":731}},"tokens_in":424,"tokens_out":804,"duration_ms":5581,"temperature":1.0,"reasoning_tokens":731,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T10:32:57.485925+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the full, nonlinear Boltzmann equation without the linearized ansatz of Eq. (11) over m_phi ∈ [10^2, 10^5] GeV and |y| ∈ [10^-5, 10^-2], and compute the final asymmetry Y_L. The cancellation predicts a flat plateau; observing a systematic slope, for example Y_L ∝ |y|^2/m at z ≈ 1 where δf reaches O(10) for m_phi = 10^5 GeV, would falsify the central claim.","supporting_citations":[],"review_version":1}