{"id":"3b8dff15-74b9-4ab6-9169-1cd13e35edfe","arxiv_id":"2608.09814","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"With one or two stiff kination epochs before radiation domination, vanilla leptogenesis can produce the observed baryon asymmetry while imprinting a blue-tilted gravitational-wave spectrum, leaving a small DECIGO-visible window.","lead":"This paper studies whether rapid 'kination' epochs before radiation domination, motivated by string-theory scalar fields, can rescue high-energy leptogenesis and leave a detectable gravitational-wave signature. It maps the parameter space where successful leptogenesis, CMB bounds on extra radiation, and future detector sensitivities overlap.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The assumption of initial thermal RHN abundance is unjustified in the stiff era; for several benchmark points the effective decay parameter is below unity, so a zero initial abundance would suppress or remove the claimed successful-leptogenesis regions.","rationale":"The reader identified the unchanged sphaleron conversion factor and thermal initial abundance as the weakest assumption. I agree that the initial-abundance part is the more serious issue, but the paper's specific failing is sharper: the quoted strong-washout regime (K = 600) does not apply at the epoch where the asymmetry is generated because the modified Hubble rate reduces the effective decay parameter below unity for several benchmarks. In weak effective washout, the final yield is not insensitive to initial conditions, so the assumed Y_N1^ini = Y_N1^eq is load-bearing rather than a harmless standard choice. The concrete zero-initial-abundance integration would settle whether the successful-leptogenesis regions survive. The PGW detectability claims also lack a stated H_e or tensor-to-scalar ratio, but that is a missing input that could be fixed by specification; the thermalization assumption is a physics assumption that, if wrong, invalidates the central leptogenesis result. I therefore recommend keeping the reader's CONDITIONAL verdict, with the condition being a zero-initial-abundance check and an explicit statement of the reheating temperature scale. This is a technical concern about the argument, not about the authors' intent; the paper is otherwise internally consistent and the qualitative mechanism is plausible.","tokens_in":20663,"tokens_out":11384,"duration_ms":105249,"concrete_test":"Re-solve Eqs. (2.19) and (2.21) with Y_N1(z_in) = 0 at z_in = M1/T_rh, taking T_rh from the upper bound in Eq. (2.11) (or from any stated reheating temperature above max(M1, few x T_R)), for the benchmark points in Tables 1-4 while keeping K and epsilon fixed. As a diagnostic, compute K_eff = K/J_1(z=1) for each point; for BM1 of Table 1 this is about 0.2. If the zero-initial-abundance Y_L falls below the required Y_L ~ (8-10) x 10^-11 x 79/28 ~ 2.3-2.8 x 10^-10 for any claimed point, that successful-leptogenesis region must be removed or relabeled as freeze-in production. This check also determines whether the equilibrium initial condition is conservative or significantly overproduces the asymmetry.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The most load-bearing assumption is that right-handed neutrinos begin with thermal equilibrium abundance, Y_N1^ini = Y_N1^eq, in the modified expansion described by Eqs. (2.19)-(2.20) and used throughout Sec. 3.1. In the stiff era, the effective decay parameter is K_eff = Gamma_1/(H_1 J_1) = K/J_1, not the quoted K = 600. For benchmark BM1 of Table 1 (w = 0.6, M1 = 1e9 GeV, T_R/M1 = 2.4e-9), J_1(z=1) ~ sqrt((M1/T_R)^0.8) ~ 2800, so K_eff ~ 0.2 < 1. Similar or even smaller values hold for most Table 1-4 benchmarks. With K_eff < 1, washout is weak and the final lepton asymmetry depends on the initial RHN abundance; there is no strong-washout memory erasure. The paper does not establish that the reheating temperature is high enough, or that the relevant decay and scattering rates are fast enough relative to the enhanced H_NS, to populate the RHNs up to equilibrium during or before the stiff epoch. If the RHN abundance starts near zero (freeze-in regime), the lepton yield is suppressed, and the successful-leptogenesis curves in Figs. 2, 5, and 7 would shift or partially disappear. This concern targets the leptogenesis half of the central claim and is independent of the unstated H_e used in the PGW amplitude.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies an early-universe cosmology in which one or two stiff epochs (equations of state w > 1/3, motivated by D-brane moduli) intervene between inflation and radiation domination. It solves the standard unflavored type-I seesaw leptogenesis Boltzmann equations with the modified Hubble rate, finding that the enhanced expansion suppresses washout and can produce the observed baryon asymmetry for lightest right-handed neutrino masses M1 in the range 10^9-10^11 GeV. It then computes the spectral energy density of primary gravitational waves for the same histories, showing that the stiff epochs imprint a blue tilt, and analyzes detectability by LISA, ET, DECIGO and other proposed detectors together with the Delta N_eff bound. The central combined result is a small region in the two-field case where successful leptogenesis, compatibility with Delta N_eff, and DECIGO detectability overlap.","tokens_in":21050,"tokens_out":23671,"duration_ms":204384,"significance":"If the assumptions are valid, the paper provides a useful joint analysis of two previously known effects: kination-assisted leptogenesis and stiff-era enhancement of the primary gravitational-wave spectrum. The analytic framework for the modified Hubble rates and the piecewise PGW spectra is transparent and reduces to standard limits, and the paper explicitly identifies benchmark points where the two observables overlap. The main novelty is the combination of the leptogenesis requirement with PGW detectability and N_eff constraints, which sharpens the parameter-space statements. However, the leptogenesis half of the analysis relies on an initial-condition assumption that is not established, and the PGW amplitude depends on an unstated inflationary Hubble scale; these issues must be resolved before the central claim can be regarded as robust.","major_comments":[{"comment":"The paper assumes Y_N1^ini = Y_N1^eq throughout, but this is not justified in the weak-washout regime actually realized by the benchmark points. The effective decay parameter is K_eff = K/J1 = 600/J1, not K = 600. For BM1 in Table 1 (w = 0.6, M1 = 10^9 GeV, TR/M1 = 2.4e-9), J1(z=1) ~ (M1/TR)^0.4 ~ 2.8e3, so K_eff ~ 0.2, and the other benchmarks also have K_eff < 1. With K_eff < 1, the final lepton asymmetry retains memory of the initial RHN abundance; if the RHN population starts near zero, the yield is set by freeze-in and is strongly suppressed. Note that the source and washout terms in Eq. (2.19) are both suppressed by 1/J1, so the faster expansion does not protect an existing asymmetry; it also suppresses production. The paper should either demonstrate that RHNs are thermalized before or during the stiff epoch (for example by specifying the reheating temperature and comparing the relevant decay and scattering rates with H_NS), or repeat the analysis with Y_N1^ini = 0. This is a decisive check for the successful-leptogenesis regions in Figs. 2, 5 and 7.","section":"Sec. 3.1, Eqs. (2.19)-(2.20), Tables 1-4"},{"comment":"The absolute amplitude of the PGW spectrum is never fixed because the inflationary Hubble scale H_e (or equivalently the tensor-to-scalar ratio r) is not stated. Eq. (2.28) gives Omega_GW^flat proportional to H_e^2/M_P^2, and H_e also enters the cutoff scale k_e in Eqs. (2.29), (2.33)-(2.34). Since the SNR curves in Figs. 4 and 7 and the Delta N_eff exclusion lines in Figs. 3, 4, 6 and 7 scale with this amplitude, all the observability statements are conditional on an unstated input. The authors should specify the value of H_e (or r) used for every plot, and should ideally show how the claimed detection regions change when H_e is varied over the allowed range up to the current constraint H_e < 5e13 GeV.","section":"Sec. 2.2.1, Eq. (2.28), Figs. 3-7"}],"minor_comments":[{"comment":"The mass scale M appearing in the upper bounds on T_rh is never defined; presumably it is the Planck mass or an O(M_P) scale, but this should be stated explicitly for the bounds to be checkable.","section":"Eqs. (2.11)-(2.12)"},{"comment":"The notation for the radiation energy density is inconsistent: Eq. (2.4) writes rho_NS in terms of rho_RD(T_R) with exponent 3w+3, while Eq. (2.8) appears to use rho_RD(T) with exponent 3w1-1. Please define rho_RD(T) explicitly and consistently to avoid confusion.","section":"Eqs. (2.4) and (2.8)"},{"comment":"The numerical values of the effective degrees of freedom g_* and g_s used in Eqs. (2.28)-(2.29) are not listed; the authors should specify them, especially because T_R can be as low as a few GeV where g_* changes significantly from its high-temperature value.","section":"Sec. 2.2, numerical inputs"},{"comment":"The text says for w = 0.6 that BM1 'touches' the N_eff-excluded region, while the summary in the same paragraph says all w = 0.6 benchmark points are safe from the N_eff bound; this apparent tension should be clarified.","section":"Fig. 4 and Sec. 3.1"},{"comment":"The tables contain only four benchmark points per scenario, but Sec. 4 repeatedly uses global language such as 'the complete parameter space' or 'the entire range'; a denser scan over (M1, TR, w) would better support these statements.","section":"Tables 1-4 and Sec. 4"},{"comment":"The PGW spectral shape depends only on the kination equation-of-state history, not directly on the leptogenesis parameters, so phrases such as 'indirect probe of high scale leptogenesis' should be phrased more carefully as a compatibility/model-selection statement rather than a direct probe of leptogenesis.","section":"Abstract and Sec. 4"},{"comment":"There are several typographical errors, including 'sigle-field' in the Fig. 4 caption and 'Shakharov' in the Introduction; these should be corrected during revision.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of JCAP and the PGW formalism is largely standard. The main technical risk is the leptogenesis initial-condition issue, which is load-bearing because all benchmark points operate in the weak-washout regime once the kination factor is included. The missing specification of H_e is also a serious reproducibility problem. Both issues are fixable in revision: the authors can add a freeze-in comparison and state the inflationary Hubble scale used in the plots. The benchmark-fitting approach is acceptable as long as it is clearly labeled, but the global language should be softened."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. The paper does what it says: it takes a modified expansion with one or two stiff eras, solves the unflavored leptogenesis Boltzmann equations, and computes the PGW spectrum from inflation using standard transfer functions. The two-step kination spectrum with a break between w1 and w2 is the genuinely new element, and it's a small but real addition to the kination-PGW literature. The qualitative physics—faster expansion suppresses washout, stiff epochs tilt the spectrum blue—is correct and well explained. The authors are also honest about the relation to refs. [48] and [56]; there is no attempt to oversell the novelty.\n\nThe soft spot is load-bearing. The paper assumes thermal equilibrium abundance for the lightest RHNs throughout (Sec. 3.1). In standard strong-washout leptogenesis with K=600 that's harmless, because final asymmetry is independent of initial conditions. But in the modified expansion the effective decay parameter is K_eff = K/J1, and with J1 ~ thousands for the low-T_R benchmarks, K_eff drops below one. For BM1 in Table 1 (w=0.6, M1=1e9 GeV, T_R/M1=2.4e-9), K_eff ~ 0.2. That's weak washout: the final asymmetry is proportional to the initial abundance, and if the RHNs start near zero (freeze-in), the yield is suppressed. The paper does not show that scattering or decay rates can populate the RHNs up to equilibrium against the enhanced Hubble rate. This directly threatens the leptogenesis half of the central claim. The PGW spectra are unaffected—they depend only on the expansion history—but the overlap regions that are said to be detectable rely on the same benchmarks.\n\nTwo smaller issues. The absolute PGW amplitude is never fixed: H_e (or r) doesn't appear in the text, so the SNR curves and the DECIGO-visible windows are not reproducible as stated. And the sphaleron conversion factor 28/79 is applied unchanged even for T_R as low as a few GeV; that's probably a minor correction but it should be checked.\n\nWho gets value: people mapping out early-universe thermal histories and GW backgrounds. The paper deserves a serious referee; it's clearly written, internally consistent, and the issue I raised is addressable. If I were the editor I'd send it out, with the referee instructed to press on the initial abundance question and the amplitude normalization.","headline":"A clean, honest parameter-space study linking kination-era leptogenesis to PGW spectra, but the thermal RHN abundance assumption makes the successful-leptogenesis regions shaky.","tokens_in":21612,"tokens_out":3324,"would_cite":false,"duration_ms":30480,"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":"A stiff pre-radiation epoch rescues high-scale leptogenesis and imprints a blue-tilted gravitational-wave background.","keywords":["leptogenesis","kination","primordial gravitational waves","spectral energy density","stiff equation of state","brane cosmology","right-handed neutrino","baryon asymmetry"],"falsifier":"A future DECIGO-class measurement that sees no blue-tilted gravitational-wave excess above the predicted spectral energy density for the $w=0.6$ benchmarks in the frequency band where the paper predicts the rise would falsify the single-field scenario; a tightened $N_{\\rm eff}$ bound that excludes the $T_R$ values in Table 1 would do the same from the other side.","tokens_in":20405,"feed_emoji":"🌊","tokens_out":9153,"duration_ms":76404,"temperature":0.7,"pith_summary":"In the standard cosmology, unflavored leptogenesis with a lightest right-handed neutrino mass between $10^9$ and $10^{11}$ GeV produces too little baryon asymmetry because the large decay parameter $K\\simeq 600$ washes the lepton number out. The paper argues that inserting one or two stiff epochs ($w>1/3$, kination-like, driven by scalar fields motivated by brane cosmology) before radiation domination fixes this: the faster Hubble expansion suppresses the washout, so the observed baryon asymmetry $Y_B\\simeq 8\\times10^{-11}$ is reproduced. The same epochs tilt the otherwise flat spectral energy density of primordial gravitational waves to the blue, with a slope fixed by $w$. The identifiable payoff is a concrete target: for a single stiff epoch with $w=0.6$, the whole mass window satisfies the $N_{\\rm eff}$ bound and is detectable by DECIGO, turning high-scale leptogenesis into a gravitational-wave astronomy question.","feed_headline":"Kination rescues high-scale leptogenesis and gives GWs a blue tilt","feed_subtitle":"A pre-radiation stiff epoch produces the observed baryon asymmetry and a spectrum DECIGO could detect.","key_machinery":"The load-bearing object is the modified Hubble rate, hence the factors $J_1$ (single field) and $J_2$ (multiple fields) that enter the Boltzmann equations for the right-handed neutrino and lepton abundances. For one field, $J_1=[1+(M_1/(T_R z))^{3w-1}]^{1/2}$ with $z=M_1/T$; for two fields, $J_2$ additionally carries the ratio $x=T_2/T_R$, so it encodes the length of the second stiff epoch. These factors reduce the effective washout parameter $\\Gamma_1/(H_1J_i)$, allowing the produced lepton asymmetry to survive. The companion mechanism is the piecewise power-law spectral energy density of primordial gravitational waves, Eq. (2.27) and Eq. (2.32), whose slope $2(3w-1)/(1+3w)$ between $k_R$ and $k_e$ is the observable fingerprint of the same stiff epochs. The cutoff at $k_e$ uses a regularized ultraviolet tail, and the integrated gravitational-wave energy density feeds the $N_{\\rm eff}$ constraint that closes the parameter space.","core_discovery":"On the paper's own terms, the discovery is that a modified pre-BBN expansion history with stiff equation-of-state epochs makes vanilla (unflavored) leptogenesis viable in a mass range where standard radiation-dominated leptogenesis fails, and simultaneously makes the primordial gravitational-wave spectrum carry a direct record of that history. Concretely, replacing $H_{\\rm RD}$ with $H_{\\rm NS}=H_{\\rm RD}[1+(T/T_R)^{3w-1}]^{1/2}$ in the Boltzmann equations inserts a factor $J_1$ that lowers the effective washout from $\\Gamma_1/H_1\\simeq600$ to $\\Gamma_1/(H_1J_1)$; for $w=0.6$ this yields the correct $Y_B$ for $M_1=10^9$--$10^{11}$ GeV with $T_R/M_1$ between $2.4\\times10^{-9}$ and $1.8\\times10^{-6}$. The gravitational-wave spectral energy density becomes $\\Omega_{\\rm GW}(k)=\\Omega_{\\rm flat}(k/k_R)^{2(3w-1)/(1+3w)}$ between the kination-to-radiation transition scale $k_R$ and the end-of-inflation scale $k_e$, so a stiff epoch is a blue tilt whose slope is set by $w$. Combining both, the paper identifies benchmark points for $w=0.6$ that are simultaneously safe under the $N_{\\rm eff}$ bound and above the SNR=10 threshold for DECIGO, and a two-epoch case ($w_2=1$, $w_1=2/3$) in which only a small triangular region near $M_1\\simeq10^{10}$ GeV survives both cuts. The claimed novelty is indirect: observing the shape of the primordial gravitational-wave spectrum could probe the same high-scale leptogenesis that cannot be tested in the laboratory.","pith_inferences":["The paper assumes equilibrium initial abundance for the right-handed neutrinos, but several successful benchmarks have $T_R$ as low as a few GeV, well below $M_1$; a non-thermal production channel would change the required $T_R/M_1$ and could shift the plotted windows.","Because the fast expansion delays charged-lepton Yukawa equilibration, the same mechanism should extend into flavored leptogenesis and relax the bound on $M_1$ even further; the paper explicitly restricts to unflavored leptogenesis.","The clean two-slope spectrum predicted for two stiff epochs could be used as a template search in future gravitational-wave data: fitting the frequency of the slope break would directly measure $T_2/T_R$, a parameter that is otherwise invisible in the baryon asymmetry alone."],"forward_implications":["A single stiff epoch with $w=0.6$ makes the whole $M_1=10^9$--$10^{11}$ GeV range of unflavored leptogenesis consistent with the observed baryon asymmetry and safe under the $N_{\\rm eff}$ bound, so no flavor model or additional CP source is needed.","If DECIGO or a comparable detector measures a blue-tilted stochastic gravitational-wave background with slope $2(3w-1)/(1+3w)$, the value of $w$ read off the slope would identify the equation of state of the pre-radiation universe.","For $w=1$ with a single field, no benchmark is simultaneously safe from the $N_{\\rm eff}$ bound and observable, so the allowed region is driven toward smaller $w$.","In the two-field case, the spectrum has two distinct slopes; observing the break frequency $k_2$ would distinguish one stiff epoch from two, which the baryon abundance alone cannot do.","Future detectors with better sensitivity could push the probe to lightest right-handed neutrino masses above $10^{11}$ GeV, making even higher-scale leptogenesis accessible."],"supporting_citations":[{"why":"Introduces the leptogenesis mechanism whose heavy-neutrino decay produces the lepton asymmetry.","marker":"[38]"},{"why":"Provides the sphaleron conversion relation $Y_B=(28/79)Y_L$ and the standard leptogenesis framework the paper modifies.","marker":"[41]"},{"why":"Derives leptogenesis in a fast-expanding universe, the basis for the modified Boltzmann equations with the $J$ factors.","marker":"[48]"},{"why":"Supplies the Casas-Ibarra parametrization used to bound the CP asymmetry $\\varepsilon$ in terms of $M_1$ and neutrino masses.","marker":"[53]"},{"why":"Shows minimal leptogenesis in brane-inspired cosmology, the direct predecessor whose scalar-field setup and parameter space are extended here.","marker":"[56]"},{"why":"Supplies the blue-tilted primordial gravitational-wave spectrum formula whose slope is set by the equation of state.","marker":"[35]"},{"why":"Provides the regularization of the short-wavelength gravitational-wave spectrum used to cut off modes above $k_e$.","marker":"[13]"},{"why":"Gives the $N_{\\rm eff}$ bound used to exclude parameter space.","marker":"[2]"},{"why":"Gives the observed baryon abundance $Y_B\\simeq(8$--$10)\\times10^{-11}$ used as the target.","marker":"[55]"}],"fun_headline_variants":["Stiff epoch leptogenesis blueprints detectable GW spectrum","Brane-modified cosmology couples leptogenesis to GW tilt","Early stiff epoch unlocks leptogenesis, signs in gravitational waves","Leptogenesis via kination, GWs turn blue for DECIGO","Modified expansion history connects baryons and gravitational waves"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that right-handed neutrinos start in thermal equilibrium and that the usual sphaleron conversion of lepton asymmetry into baryon asymmetry is unchanged, even when radiation domination begins at temperatures as low as a few GeV.","fun_headline_variants_meta":{"raw":{"variants":["Stiff epoch leptogenesis blueprints detectable GW spectrum","Brane-modified cosmology couples leptogenesis to GW tilt","Early stiff epoch unlocks leptogenesis, signs in gravitational waves","Leptogenesis via kination, GWs turn blue for DECIGO","Modified expansion history connects baryons and gravitational waves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000255,"raw_usage":{"total_tokens":1653,"prompt_tokens":1107,"completion_tokens":546,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":723,"completion_tokens_details":{"reasoning_tokens":475}},"tokens_in":723,"tokens_out":546,"duration_ms":5154,"temperature":1.0,"reasoning_tokens":475,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T10:15:56.608105+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A future DECIGO-class measurement that sees no blue-tilted gravitational-wave excess above the predicted spectral energy density for the $w=0.6$ benchmarks in the frequency band where the paper predicts the rise would falsify the single-field scenario; a tightened $N_{\\rm eff}$ bound that excludes the $T_R$ values in Table 1 would do the same from the other side.","supporting_citations":[],"review_version":1}