{"id":"2d8de857-1f74-4fa2-bac3-d742259fae7b","arxiv_id":"2411.15849","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Closed-form probability distributions for the curvature perturbation are derived for four exactly solvable curvaton potentials, and the mass-decay-rate parameter space is split into regions by the probability of producing 10% to 100% of the observed signal.","lead":"This paper derives analytical probability distributions for the curvature perturbation created by a curvaton field in four exactly solvable models of stochastic inflation. It then maps the curvaton mass and decay rate into regions where the curvaton is likely to produce a non-negligible part of the observed perturbations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Non-quadratic curvaton distributions rely on the quadratic conversion (3.3), which the paper itself flags as approximate; the high-probability parameter regions may be an artifact.","rationale":"The central claim is that exact stochastic distributions can be converted into closed-form curvature distributions and that parameter-space regions with high probability exist. This claim is strongest for the non-quadratic potentials, which are the main novelty compared to earlier work on the quadratic case. The conversion step (3.3) is the only place where the potential shape enters the map from φ to ζ; the stochastic distributions f_∞(φ) themselves are exact for each potential. If the quadratic conversion is corrected, the distributions (3.15), (3.17), (3.24) change, and the parameter-space topology (Fig. 5) and probability maps (Figs. 4, 6) are not reliable. The paper's own caveat in Sec. 3.1 and the Conclusion confirms the approximation is uncontrolled for the anharmonic cases. The proposed test directly quantifies the error by recomputing the observable probability with the full potential. The harmonic case (Eq. 3.12, Fig. 3) does not suffer from this issue and reproduces earlier results, so the paper retains value; hence the verdict should remain conditional pending this check.","tokens_in":39152,"tokens_out":16994,"duration_ms":147846,"concrete_test":"Compute the exact ζ(φ) for the Scarf potential (s=0.5, t=0) using the full energy density V(φ) = (3/4π²) H⁴ [−s ln cos(α√(8π²/H²) φ)] in the sudden-decay formalism, including the non-zero equation of state during oscillations, and re-derive P(ζ) by the change of variables. Then recompute P(0.1 ζ_obs < ζ < ζ_obs) over the (Γ, m) grid of Fig. 6 (left) and compare the locus where the probability exceeds 0.8 with the black region shown. If this locus shifts by more than a factor of 2 in Γ or m (or vanishes), the quadratic conversion is load-bearing and the claimed parameter-space result is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing assumption is that the ζ–φ conversion, Eq. (3.3) (and its inverse Eq. (3.4)), derived for a purely quadratic curvaton potential under the sudden-decay approximation, remains valid for the radial harmonic oscillator and the trigonometric Scarf/Rosen–Morse potentials. The paper itself states in Sec. 3.1 that this relation 'ceases to be an exact relation' for non-quadratic potentials and that corrections 'can be relevant' (citing Refs. [18,26,27,29,33,34,37,38,90,91]); the Conclusion repeats this caveat. Yet the closed-form distributions (3.15), (3.17), (3.24) and the parameter-space maps (Figs. 4–6), including the claim of >80% probability in certain (Γ, m) regions, are derived entirely from this approximate conversion. For the Scarf potential with s=0.5, the anharmonicity is not small: expanding V(φ) about the minimum gives a quartic term with coefficient (2π²/3)(m/H)², and the typical field displacement ⟨φ²⟩^{1/2} ~ H/m in the high-probability region makes the leading correction of order (φ/H)², which is O(1) for m/H ~ 0.1. Thus the derived P(ζ) is not the physical distribution for those cases, and the claimed high-probability regions may be an artifact of the quadratic approximation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper applies four exactly solvable stationary solutions of stochastic inflation (harmonic oscillator, radial harmonic oscillator, trigonometric Scarf, and trigonometric Rosen–Morse) to the curvaton scenario. It uses the standard quadratic-curvaton relation ζ(φ), Eq. (3.3), to change variables and obtain closed-form distributions P(ζ) for each case (Eqs. 3.12, 3.15, 3.17, and 3.24), imposes normalization on each distribution, and computes the probability P(cζ_obs < ζ < ζ_obs) with c = 0.1 in the (Γ, m) parameter space. The paper finds regions where the probability of a non-negligible curvaton contribution exceeds 80%. For the trigonometric potentials, a finite-field-domain condition x < 1 is introduced, leading to a six-region parameter-space classification (Fig. 5 and Table 2).","tokens_in":39528,"tokens_out":10499,"duration_ms":101743,"significance":"If the results hold, the paper provides a useful analytical demonstration of how exact non-Gaussian spectator distributions can be propagated into curvature-perturbation statistics, extending the earlier quadratic-potential analysis of Refs. [64,65]. Strengths include the explicit closed-form expressions, the consistently enforced normalizations, the transparent order-of-magnitude checks, and the detailed appendix on the trigonometric-Scarf moments. The main limitation is that the central novelty for the three non-quadratic potentials rests on an approximation that the manuscript itself flags, so the quantitative parameter-space statements for those cases need additional support or qualification.","major_comments":[{"comment":"The conversion relation ζ(φ) in Eq. (3.3) is acknowledged by the authors to be exact only for a purely quadratic potential, yet it is used to derive the distributions (3.15), (3.17), and (3.24) for the radial harmonic oscillator and the two trigonometric potentials. The issue is quantitative, not merely formal: for the Scarf potential with s = 0.5, the potential about the minimum is V(φ) = (1/2)m²φ² + (2π²/3)(m/H)²φ⁴ + O(φ⁶), and the typical displacement from Eq. (3.23) gives ⟨φ²⟩^{1/2}/H ~ 0.24 for the field values dominating the distribution. The quartic term is therefore an O(1) correction to the quadratic term in this regime, so the derived P(ζ) and the >80% probability regions in Figs. 4–6 cannot be considered the physical distributions for the non-quadratic cases without further justification. I request either that the main claims be restricted to the quadratic harmonic case, or that the paper provide a quantitative estimate or validation of the correction terms over the plotted parameter region.","section":"§3.1, Eq. (3.3)"},{"comment":"The parameter-space maps assume the curvaton has reached the stationary distribution. For the harmonic potential, Eq. (2.13b) implies a relaxation timescale of order 3H²/m² e-folds. The first high-probability branch in Fig. 3, and the corresponding peak in Fig. 2, is located at m/H ~ 10⁻⁵, which requires roughly 10¹⁰ e-folds to relax. That is orders of magnitude larger than the 'much longer than N★ = 60' stated below Eq. (2.14). Since Ref. [64] shows that finite-duration effects can be important in stochastic spectator dynamics, the paper should either state the required total inflation duration for the plotted regions or demonstrate that the >80% probability claim is robust under finite-duration initial conditions.","section":"§3.3.1 and Fig. 3"}],"minor_comments":[{"comment":"The distribution P(ζ) is defined and normalized on the positive interval [0, ζM], while the original stationary φ-distribution is supported on the full real line (for the harmonic cases) and is symmetric. The quantity plotted is therefore effectively the distribution of |ζ|. This should be stated explicitly in the probability interpretation.","section":"§3.2"},{"comment":"The statement after Fig. 3 that 'implementing a different scale merely shifts the high-probability regions' is not exact for the trigonometric cases, because the allowed-region condition x < 1 in Eq. (3.19) depends on H through α/ζM.","section":"§3.3.2, Eq. (3.19)"},{"comment":"Eq. (B.2) is introduced as a conjecture, but it is used to derive the variance formula (3.23) and the fourth-moment formula (B.15), which are quoted as analytical results in the main text. Please provide a proof or clearly mark these moments as numerically verified rather than exact.","section":"Appendix B, Eq. (B.2)"},{"comment":"Minor editorial issues: in §3.3.1, 'excess the observational amplitude' should be 'exceed'; in the same section the abbreviation 'r_dec' is used inconsistently with 'r_decay' elsewhere; and the color/grayscale conventions in Figs. 3, 4, and 6 are described only in the text, not in the captions.","section":"§3.3.1 and §3.3.2"},{"comment":"The analysis of the trigonometric potentials is restricted to the symmetric t = 0 limit, and the asymmetric cases of Eqs. (3.16) and (3.24) are not explored. The conclusions should be phrased as applying to the symmetric slice, not to the full trigonometric Scarf or Rosen–Morse parameter space.","section":"§3.3.2, Fig. 5"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest and well structured, and the authors explicitly flag the two main limitations I raise: the non-quadratic conversion approximation and the stationarity assumption. These are therefore fixable within the scope of a revision, by either narrowing the claims or adding quantitative control. I would not reject the paper, but the central quantitative claims for the non-quadratic potentials are not yet established."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a mostly careful, transparent application of exact stochastic-inflation distributions to the curvaton, and the closed-form P(zeta) for the three non-quadratic potentials are genuinely new. The harmonic case reproduces Lerner and Melville qualitatively, so that piece is not new; the new content is the radial harmonic oscillator, trigonometric Scarf, and trigonometric Rosen–Morse distributions, plus the six-region topology of the (Gamma, m) plane that arises from the bounded field range. The derivation is straightforward: change variables in the stationary PDF, enforce normalization, and integrate. No parameter is fit to zeta_obs; the inputs are fixed. That is good, and so is the willingness to state in Sec. 3.1 and the Conclusion that the conversion 'ceases to be exact' for non-quadratic potentials and that corrections 'can be relevant.'\n\nSoft spots, in order of severity. (1) The load-bearing conversion is the main issue. The stress-test note quantifies what 'can be relevant' means: for the Scarf potential with s=0.5 and m/H~0.1, the quartic term in the potential gives an O(1) correction to the field-displacement relation, so the derived P(zeta) and the >80% probability regions in Figs 4–6 are not the physical distributions. The authors are honest about this limitation, but the abstract presents the distributions without that caveat. A serious referee should ask for a quantitative estimate of the error, or a restricted claim. (2) Appendix B's conjectured integral (B.2) is unproved; the variance and fourth-moment results rest on it. The limits match expectations, so it is probably right, but 'conjectured' is a weak foundation for a printed result. (3) The paper also assumes the stationary regime and ignores f_NL; those are stated limitations, not flaws.\n\nThe paper does not resolve any long-standing problem, but it supplies useful closed forms for model builders and for PBH/non-Gaussianity studies. It deserves a serious referee. I would send it to peer review with a request to quantify the approximation error and to prove or numerically check the conjectured integral. The harmonic case alone is solid enough to justify the effort.","headline":"Genuinely new closed-form curvaton distributions from stochastic inflation, but the trigonometric cases lean on the quadratic zeta–phi conversion the paper itself flags as approximate, so the high-probability maps are provisional.","tokens_in":39995,"tokens_out":3868,"would_cite":true,"duration_ms":36509,"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":"The paper derives closed-form probability distributions for the curvature perturbation produced by a curvaton in stochastic inflation and shows that, in a broad region of the decay-rate–mass plane, a single curvaton can account for the…","keywords":["curvaton","stochastic inflation","curvature perturbation distribution","exact stationary solutions","Fokker-Planck equation","Starobinsky-Yokoyama equilibrium","parameter space","primordial perturbations"],"falsifier":"Run a lattice or $\\delta N$ simulation of a spectator field in the trigonometric Scarf potential $v(y) = -s \\ln\\cos y$ with the same $(\\Gamma,m)$ values, converting to $\\zeta$ through the actual energy-density transfer rather than Eq. (3.3), and check whether $P(0.1\\,\\zeta_{\\rm obs} < \\zeta < \\zeta_{\\rm obs}) > 0.8$ in the claimed strips; if not, the central parameter-space claim fails for the non-quadratic cases.","tokens_in":38927,"feed_emoji":"🌌","tokens_out":9251,"duration_ms":74129,"temperature":0.7,"pith_summary":"The paper aims to show that the exact stationary probability distributions of a spectator scalar field in stochastic inflation can be turned, by a change of variables, into closed-form distributions for the curvature perturbation produced by a curvaton. It derives $P(\\zeta)$ for four exactly solvable potentials and uses the resulting expressions to map the $(\\Gamma, m)$ parameter space, asking where a curvaton has high probability of contributing between 10% and 100% of the observed $\\zeta_{\\rm obs} = 10^{-5}$. The central finding is that such high-probability regions exist, including a branch where the curvaton dominates the universe before decay and a branch where it does not. A sympathetic reader would care because the paper turns a model-dependent early-universe question into an analytically tractable probability statement.","feed_headline":"Curvaton can explain all CMB ripples with >80% probability","feed_subtitle":"A change of variables turns stochastic-inflation distributions into probability maps for the curvaton's share.","key_machinery":"The master formula is Eq. (3.9), obtained from $P(\\zeta) = \\sum_{\\lambda=\\pm} f_\\infty(\\phi_\\lambda)\\,|\\mathrm{d}\\phi_\\lambda/\\mathrm{d}\\zeta|$ using the inverted relation $\\phi_\\pm = (H/6\\pi\\zeta_{\\rm M})\\sqrt{(1\\pm Y)/(1\\mp Y)}$, where $Y = \\sqrt{1-(\\zeta/\\zeta_{\\rm M})^2}$ and $\\zeta_{\\rm M}^2 = (H/M_{\\rm P})^2/[(12\\pi)^2\\sqrt{\\Gamma/m}]$. The Jacobian has a square-root singularity at $\\zeta = \\zeta_{\\rm M}$, and the two branches correspond to curvaton domination ($r_{\\rm decay}\\simeq 1$) and subdominance. The input distributions $f_\\infty(\\phi)$ are the late-time equilibria of the Fokker–Planck equation, the Starobinsky–Yokoyama formula, and the model dependence enters through the mass parameter $\\alpha$ and the shape parameters $\\ell$, $s$, and $t$.","core_discovery":"The paper claims that, for each of the four stationary spectator distributions — harmonic oscillator, radial harmonic oscillator, trigonometric Scarf, and trigonometric Rosen–Morse — substituting the inverted field–curvature relation (3.4) into probability conservation yields a closed-form $P(\\zeta)$: Eqs. (3.12), (3.15), (3.17), and (3.24). It further claims that in the $(\\Gamma,m)$ plane, with $H/M_{\\rm P} = 10^{-5}$ and $\\zeta_{\\rm obs} = 10^{-5}$, there are parameter regions where $P(0.1\\,\\zeta_{\\rm obs} < \\zeta < \\zeta_{\\rm obs})$ exceeds 80%, meaning the curvaton alone can account for the observed curvature perturbation with high probability. For the three non-quadratic potentials the conversion relation (3.3), exact only for a quadratic potential, is used as a stated first approximation.","pith_inferences":["Applying the same change-of-variables step to the time-dependent transition probabilities would reveal how the >80% regions depend on the initial curvaton value; the paper explicitly leaves this for future work.","The Jacobian singularity at $\\zeta=\\zeta_{\\rm M}$ concentrates probability near the maximum allowed curvature perturbation; if that tail feeds primordial-black-hole formation thresholds, the stochastic distributions could also predict an enhanced PBH abundance even when the mean $\\zeta$ is small.","The exact distributions could be used to derive the full probability density of $f_{\\rm NL}$ or of the curvaton's energy fraction, rather than just the integrated probability interval, giving sharper observational tests.","A numerical computation of $\\zeta$ from the exact nonlinear field equations for the cosine/logarithmic potentials would show whether the 80% contours persist; until then the trigonometric-case maps should be read as illustrative of the method."],"forward_implications":["If the central claim is correct, a curvaton with $m/M_{\\rm P}\\sim 10^{-10}$ at $H/M_{\\rm P}=10^{-5}$ can generate the full observed curvature perturbation with probability above 80% on the curvaton-dominated branch, essentially independent of the decay rate.","The probability maps give a quantitative naturalness criterion: parameter sets in the low-probability regions are disfavored because, with high probability, they produce less than 10% of $\\zeta_{\\rm obs}$ or exceed it.","For the radial harmonic oscillator, $\\ell>0$ broadens the high-probability region, while for the trigonometric Scarf potential increasing $s$ narrows it, with $s\\to\\infty$ recovering the quadratic-potential result.","For the trigonometric potentials the finite field range splits the parameter space into six regions; parameters in region F necessarily overproduce the total curvature perturbation and are excluded.","Because the stationary state is assumed, the results apply only when inflation lasts much longer than the relaxation time of the curvaton; otherwise initial conditions would matter."],"supporting_citations":[{"why":"Establishes the curvaton paradigm in which a spectator field can generate the curvature perturbation.","marker":"[11]"},{"why":"Provides the stochastic-inflation formalism, the Langevin equation, and the Fokker–Planck equation used throughout.","marker":"[59]"},{"why":"Supplies the closed-form stationary distributions of the four exactly solvable spectator potentials that the paper converts into P(ζ).","marker":"[60]"},{"why":"Introduced the method of deriving the curvature perturbation distribution from a stochastic curvaton distribution, which this paper follows.","marker":"[64]"},{"why":"Performed the analogous (Γ,m) parameter-space probability study for a quadratic-potential curvaton, the result this paper extends.","marker":"[65]"},{"why":"Gives the equilibrium (stationary) distribution that serves as the input f∞(φ) in the conversion.","marker":"[85]"},{"why":"Provides the leading-order quadratic-potential relation ζ = r_decay H/(3π φ) used in the conversion.","marker":"[89]"}],"fun_headline_variants":["Curvaton odds: >80% for all CMB ripples","Closed-form curvaton distributions pinpoint >80% regions","Stochastic inflation yields closed-form curvaton odds","Curvaton alone: >80% odds for all CMB ripples"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conversion formula (3.3), exact only for a purely quadratic potential with sudden decay, is applied to the three non-quadratic potentials; if the neglected corrections are sizable, the derived $P(\\zeta)$ and probability maps for those cases are not the physical distributions.","fun_headline_variants_meta":{"raw":{"variants":["Curvaton odds: >80% for all CMB ripples","Closed-form curvaton distributions pinpoint >80% regions","Stochastic inflation yields closed-form curvaton odds","Curvaton alone: >80% odds for all CMB ripples"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001645,"raw_usage":{"total_tokens":6462,"prompt_tokens":796,"completion_tokens":5666,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":412,"completion_tokens_details":{"reasoning_tokens":5595}},"tokens_in":412,"tokens_out":5666,"duration_ms":38216,"temperature":1.0,"reasoning_tokens":5595,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:50:08.194838+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a lattice or $\\delta N$ simulation of a spectator field in the trigonometric Scarf potential $v(y) = -s \\ln\\cos y$ with the same $(\\Gamma,m)$ values, converting to $\\zeta$ through the actual energy-density transfer rather than Eq. (3.3), and check whether $P(0.1\\,\\zeta_{\\rm obs} < \\zeta < \\zeta_{\\rm obs}) > 0.8$ in the claimed strips; if not, the central parameter-space claim fails for the non-quadratic cases.","supporting_citations":[{"cited_title":"Quantifying the 'naturalness' of the curvaton model","cited_arxiv_id":"1402.3176","evidence_quote":"Performed the analogous (Γ,m) parameter-space probability study for a quadratic-potential curvaton, the result this paper extends."}],"review_version":1}