{"id":"8c3ef55d-f3ab-4acb-a9ca-f11df03b7bd9","arxiv_id":"2509.10851","paper_version":2,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Derives exact non-Gaussian PDF for curvaton curvature perturbations and computes PBH abundances plus induced GW spectra in a self-consistent fluctuation model without per-scale variance matching.","lead":"The paper derives the exact probability density function for curvature perturbations in the curvaton scenario from the sudden-decay relation using a branchwise change of variables on Gaussian field fluctuations, then computes the non-perturbative tail for primordial black hole formation fractions. It replaces scale-by-scale fitting with a self-consistent model to generate PBH mass functions and induced gravitational wave spectra for comparison with current constraints and LSA","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Sudden-decay relation defines the ζ(σ) map before branchwise change of variables; gradual decay would alter the far tail","rationale":"The reader's weakest assumption is precisely the load-bearing step. The full manuscript supplies the analytic steps and self-consistent spectrum, but the physical mapping remains an idealization whose violation directly impacts the tail used for both PBH abundance and induced GW spectra. The proposed numerical check is decisive and independent of the paper's internal derivations.","tokens_in":1827,"tokens_out":386,"duration_ms":21711,"concrete_test":"Implement a numerical integration of the curvaton background and perturbation equations with finite decay width Γ/H ≈ 0.1–1 around decay; sample 10^6 Gaussian initial σ fluctuations, evolve each to the end of decay, extract the resulting ζ distribution, and compare its cumulative probability above the PBH threshold δ_c ≈ 0.67 to the analytic sudden-decay tail; a discrepancy exceeding a factor of 3 in β would falsify the exact non-perturbative claim.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central construction obtains the full non-Gaussian PDF of ζ by an explicit change-of-variables from Gaussian σ fluctuations using the exact sudden-decay formula ζ = (2/3) ln(1 + (3/2)σ/σ_d) (or equivalent). This mapping is inserted directly prior to computing the Jacobian and support of the distribution. In any realistic model the curvaton decay rate Γ is finite, so the relation between the pre-decay field value and the final ζ acquires additional smoothing from the integrated decay dynamics; this smoothing is known to suppress the extreme positive tail that controls the PBH collapse fraction. The paper replaces scale-by-scale variance matching with a self-consistent spectrum but retains the instantaneous-decay premise throughout the analytic and numerical pipeline.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript studies the curvaton scenario for primordial black hole formation by deriving the curvature perturbation from the exact sudden-decay relation. It obtains the full probability density function of the curvature perturbation through a branchwise change of variables from Gaussian field fluctuations and evaluates the PBH formation fraction from the exact non-perturbative tail. The results are compared to Gaussian and local quadratic non-Gaussian cases. A self-consistent model is used where the dimensionless field fluctuation spectrum is specified once, allowing direct computation of smoothed variances and collapse fractions without additional per-scale fitting. The induced gravitational wave background is generated from the linear spectrum, and the mass functions are compared to constraints with GW spectra shown against observational windows.","tokens_in":2005,"tokens_out":609,"duration_ms":103808,"significance":"If the central derivations hold, the work provides a non-perturbative analytic treatment of non-Gaussianity in the curvaton model for PBH formation, valuable due to the exponential sensitivity of PBH abundance to the tail. The explicit step-by-step analytic forms for the Jacobian, normalization, support, and small-fluctuation expansion, together with the self-consistent fluctuation model that avoids per-scale variance matching, are strengths. Confrontation with current constraints and display of induced GW spectra against PTA/LISA/DECIGO windows add observational relevance.","major_comments":[{"comment":"The central construction obtains the PDF of ζ by an explicit change of variables using the sudden-decay mapping ζ = (2/3) ln(1 + (3/2)σ/σ_d) (or equivalent) inserted directly before the Jacobian and support calculation. This assumes instantaneous decay. Realistic models with finite decay rate Γ introduce smoothing from integrated decay dynamics that suppresses the extreme positive tail controlling the PBH collapse fraction. This premise is load-bearing for the claim of an 'exact non-perturbative tail' and the resulting PBH mass functions.","section":"Derivation of the curvature perturbation from the exact sudden-decay relation"}],"minor_comments":[{"comment":"The self-consistent model reduces per-scale fitting but still treats the dimensionless field fluctuation spectrum as a free parameter whose value is chosen once; clarify how this choice is motivated or varied to test robustness of the final mass functions and GW spectra.","section":"Self-consistent curvaton fluctuation model"},{"comment":"Ensure that the support of the distribution and the branchwise handling in the change-of-variables step are illustrated with an explicit equation or figure for the reader to follow the normalization.","section":"Probability density function derivation"},{"comment":"The figures displaying mass functions against the conservative constraint envelope and GW spectra against sensitivity windows are useful; verify that all curves are labeled with the corresponding model (exact curvaton, Gaussian benchmark, quadratic) for clarity.","section":"Results and figures"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the careful reading and constructive comments. We address the major comment below and indicate the revisions we will implement.","responses":[{"response":"We agree that the sudden-decay approximation is central to our analytic derivation and that a finite decay width would introduce smoothing that could suppress the far positive tail. Our manuscript explicitly adopts the sudden-decay relation as a standard analytic limit in the curvaton literature, allowing the exact branchwise change-of-variables PDF to be derived without truncation. We will revise the text to add a dedicated paragraph in Section 2 (or the discussion) that (i) states the approximation explicitly, (ii) cites literature on finite-Γ effects, and (iii) notes that the reported tail represents the limiting case of maximal non-Gaussian enhancement. This clarifies the scope without altering the central derivations.","revision_made":"yes","referee_comment":"The central construction obtains the PDF of ζ by an explicit change of variables using the sudden-decay mapping ζ = (2/3) ln(1 + (3/2)σ/σ_d) (or equivalent) inserted directly before the Jacobian and support calculation. This assumes instantaneous decay. Realistic models with finite decay rate Γ introduce smoothing from integrated decay dynamics that suppresses the extreme positive tail controlling the PBH collapse fraction. This premise is load-bearing for the claim of an 'exact non-perturbative tail' and the resulting PBH mass functions."}],"tokens_in":1484,"tokens_out":353,"duration_ms":32519,"standing_objections":["Quantitative evaluation of the tail suppression under a finite decay rate Γ, which would require a separate integration over the decay dynamics outside the present sudden-decay framework."]},"desk_editor":{"model":"grok-4.3","letter":"The main takeaway is that they derive the exact non-Gaussian tail for ζ by mapping Gaussian σ fluctuations through the sudden-decay formula ζ = (2/3) ln(1 + (3/2)σ/σ_d) and handling the branches, Jacobian, support, and normalization analytically. They then drop the usual scale-by-scale variance matching in favor of a single dimensionless field fluctuation spectrum that determines the smoothed variance, the collapse fraction, and the linear two-point function for the gravitational waves all at once. That setup lets them compare the exact result directly to the Gaussian benchmark and to an exact local quadratic case without Edgeworth truncation, and they plot the resulting mass functions against current PBH bounds plus the induced spectra against PTA, LISA, and DECIGO windows from one consistent pipeline. The analytic display of the small-fluctuation limit and the normalization step is useful for anyone who wants to check the tail behavior by hand. The self-consistent spectrum choice removes one layer of ad-hoc adjustment, which is a modest but real improvement over the usual approach. The soft spot is the sudden-decay premise. The mapping is inserted before the change of variables, and the paper carries that assumption through the entire calculation. In models with finite decay rate the integrated dynamics smooth the relation and suppress the extreme positive tail that sets the PBH fraction, so the abundances here are likely higher than they would be under gradual decay. The initial form chosen for the dimensionless spectrum remains a model parameter even if it is not retuned per scale. This is for people already working on curvaton non-Gaussianity or small-scale PBH constraints who need a concrete non-perturbative PDF rather than another Edgeworth expansion. A reader who cares about the technical step of the branchwise derivation will get value from the explicit steps. It deserves a serious referee because the derivation is laid out clearly enough to check and the self-consistent pipeline is a step forward, even if the decay approximation requires discussion in review.","headline":"The paper works out the full non-Gaussian PDF for the curvaton curvature perturbation via explicit branchwise change of variables from the sudden-decay map, then runs PBH and induced-GW calculations from one self-consistent spectrum instead of per-scale fitting.","tokens_in":2529,"tokens_out":492,"would_cite":false,"duration_ms":26169,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":{"model":"grok-4.3","evidence":[{"relation":"unclear","rs_module":"IndisputableMonolith/Cost/FunctionalEquation.lean","rs_theorem":"washburn_uniqueness_aczel","paper_passage":"Solving Eq. (15) gives e^{3ζ_χ} = Y(ζ) ≡ 1/Ω (e^{3ζ} + (Ω-1)e^{-ζ}); combining with e^{3ζ_χ} = (1 + δχ/χ̄)^2 yields the mapping δχ(ζ) = χ̄ (±√(Y(ζ)−1)) and the change-of-variables PDF P_ζ(ζ) = Σ_s P_χ(δχ_s(ζ)) |dδχ_s/dζ|."}],"headline":"Curvaton sudden-decay PDF via branchwise Jacobian change-of-variables is standard multi-field inflation machinery","alignment":"orthogonal","rationale":"The paper's core construction derives the full non-Gaussian P(ζ) from the exact sudden-decay algebraic relation (Eq. 15-21) by explicit branchwise inversion and Jacobian transformation of a Gaussian δχ field. This is conventional curvaton phenomenology with no invocation of recognition cost J(x), ratio symmetry, φ-ladder spacings, or parameter-free forcing. The exponential forms appearing in Y(ζ) are incidental to the model and do not parallel RS cosh-cost or J-cost identities.","tokens_in":47187,"confidence":"moderate","tokens_out":321,"duration_ms":11529,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"The exact sudden-decay relation in the curvaton scenario yields the full non-Gaussian probability distribution of curvature perturbations from which primordial black hole abundances are computed directly.","keywords":["primordial black holes","curvaton scenario","non-Gaussianity","curvature perturbation","sudden decay","probability density function","induced gravitational waves"],"falsifier":"A numerical simulation or observation showing that the actual mapping from curvaton field to curvature perturbation deviates substantially from the sudden-decay formula on the relevant scales would invalidate the exact tail used for the black-hole fraction.","tokens_in":2711,"feed_emoji":"","tokens_out":703,"duration_ms":33362,"temperature":0.7,"pith_summary":"The paper establishes that primordial black hole formation depends exponentially on the far tail of the curvature perturbation distribution, making it a sensitive test of non-Gaussianity. The authors derive the curvature perturbation exactly from the sudden-decay relation in the curvaton model and obtain the complete probability density function by performing an explicit branchwise change of variables on the underlying Gaussian curvaton field fluctuations. This non-perturbative method is then used to evaluate the black hole formation fraction without Edgeworth truncations or other approximations. The same consistent fluctuation spectrum also determines the induced gravitational wave background, allowing direct comparison with current constraints and future detector sensitivities.","feed_headline":"Exact curvaton decay fixes non-Gaussian tails for black hole formation","feed_subtitle":"Branchwise change of variables yields the full curvature distribution and consistent gravitational wave signals without scale-by-scale fits.","key_machinery":"The exact sudden-decay relation between the curvaton field value and the curvature perturbation, which enables an explicit branchwise change of variables that maps the Gaussian field distribution onto the complete non-Gaussian curvature perturbation probability density function.","core_discovery":"In the curvaton scenario the curvature perturbation is derived from the exact sudden-decay relation, the full probability density function is obtained through an explicit branchwise change of variables from the Gaussian curvaton-field fluctuation, and the primordial black-hole formation fraction is evaluated from the exact non-perturbative tail; the same model supplies a self-consistent fluctuation spectrum that replaces scale-by-scale variance matching and generates the induced gravitational-wave background.","pith_inferences":["The branchwise method could be extended to other sources of non-Gaussianity such as multi-field inflation to obtain similarly exact tails.","Precision measurements of the stochastic gravitational wave background by future detectors could distinguish this curvaton signature from other early-universe scenarios.","The exact tail calculation may tighten or relax existing limits on curvaton parameters once applied to updated observational envelopes."],"forward_implications":["The primordial black hole mass function follows directly from the exact non-perturbative tail without perturbative truncation or additional fitting.","Induced gravitational wave spectra are generated from the linear curvaton two-point function within the same self-consistent model.","Predictions differ measurably from both the pure Gaussian case and the exact local quadratic benchmark.","Current conservative constraints on primordial black holes can be confronted with mass functions computed once for the entire spectrum."],"fun_headline_variants":["Exact curvaton decay maps non-Gaussian PBH tails","Branchwise variables yield full curvaton PBH PDF","Self-consistent curvaton spectrum predicts PBH and GW signals","Exact sudden decay computes precise PBH formation fractions"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The curvaton decays suddenly so that the curvature perturbation is given directly by the field value before decay without any gradual transition or smoothing corrections.","fun_headline_variants_meta":{"raw":{"variants":["Exact curvaton decay maps non-Gaussian PBH tails","Branchwise variables yield full curvaton PBH PDF","Self-consistent curvaton spectrum predicts PBH and GW signals","Exact sudden decay computes precise PBH formation fractions"]},"model":"grok-4.3","cost_usd":0.008817,"raw_usage":{"total_tokens":4000,"prompt_tokens":731,"num_sources_used":0,"completion_tokens":62,"cost_in_usd_ticks":88174500,"prompt_tokens_details":{"text_tokens":731,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":3207,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":731,"tokens_out":62,"duration_ms":30767,"temperature":1.0,"reasoning_tokens":3207,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-05-18T16:58:07.904334+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A numerical simulation or observation showing that the actual mapping from curvaton field to curvature perturbation deviates substantially from the sudden-decay formula on the relevant scales would invalidate the exact tail used for the black-hole fraction.","supporting_citations":[],"review_version":1}