{"id":"a5e3fd6b-7f2d-4fd7-8da9-180e9d2f86db","arxiv_id":"2501.04774","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":7,"one_line_summary":"The paper derives analytic dark matter freeze-in yields for arbitrary power-law reheating histories and maps the gravitational production parameter space.","lead":"This paper computes how much dark matter is created during the poorly understood reheating phase after inflation, using a general two-parameter model of that phase. It provides formulas that let physicists scan many reheating histories and pin down the dark matter abundance for gravitational and inflaton-decay production.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reheating background defined by Eqs. (2.2)–(2.7) does not solve the stated Boltzmann pair (2.4)–(2.5): nonzero Γ depletes ρφ, so the decay yields (3.22)/(3.24) can exceed the available inflaton number and the relic contours inherit an uncontrolled O(1) error.","rationale":"The reader's weakest assumption was that a single power-law reheating is physically representative; my concern is more specific and arguably more damaging. Even granting the power-law ansatz, the paper's background is not a solution of the Boltzmann equations it writes down. The derivation of Γ from Eq. (2.5) is fine as a definition, but Eq. (2.4) is then violated for every Γ≠0, because a pure power-law ρφ has zero decay-loss term. The numerical consequence is visible in the α=3/8 limit: computed decay yields exceed the total initial inflaton particle number. This is not a matter of external consensus; it is an internal consistency failure. The scaling exponents k_c and the log/power-law classification may survive a correct treatment, but the prefactors that set the relic-density contours in Figs. 4–6 would change by O(1) or worse, and Eq. (3.22) as printed can even produce a negative yield. Therefore I would not reject the paper outright, but I would make ACCEPT conditional on a numerical two-fluid check and a corrected derivation of the decay formulas.","tokens_in":18636,"tokens_out":29237,"duration_ms":285844,"concrete_test":"Solve Eqs. (2.4)–(2.5) numerically with Γ(a) from Eq. (2.6) for ω=0, α=3/8 (and, for robustness, α=0 and 3/4), starting at the same T_I and H_I, and compute the DM yield from γ_d=C Γ ρφ at the epoch where ρφ=ρR. Compare with Eq. (3.22). If the full-system yield differs from Eq. (3.22) by more than ~10%, or if Eq. (3.22) is negative for α=3/8, the decoupled power-law background is inconsistent and the decay-based contours in Figs. 4–6 must be revised.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Section 2 assumes ρφ(a) ∝ a^{-3(1+ω)} and H(a)=H_rh(a_rh/a)^{3(1+ω)/2} throughout reheating, while Eq. (2.6) gives nonzero Γ. Substituting these into Eq. (2.4) leaves the left-hand side equal to zero, not -Γρφ; energy is not conserved unless the power law already includes decay, which it does not. The paper fixes Γ from the radiation equation only and never enforces Eq. (2.4). The error is not confined to the boundary: for the nominal constant-width case ω=0, α=3/8, Eq. (2.6) gives Γ/H_rh=5/2, so ∫Γ dt from a_I to a_rh equals 5/3. The decay-generated comoving yield from Eq. (3.8) is then (5/3)×(Br x Nφ,initial), exceeding the maximum number of DM particles that the initial inflaton population can produce when Br=1. The same inconsistency feeds the normalization of the scattering yields near a_rh, where ρR is no longer negligible, and the printed Eq. (3.22) already becomes negative for ω=0, α=3/8. Because Figs. 4–6 convert these prefactors into relic-density contours, the quantitative maps are not secured by the stated assumptions.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a general power-law parametrization of the reheating epoch, with inflaton energy density scaling as a^{-3(1+ω)} and SM temperature scaling as a^{-α}, and derives analytic expressions for UV freeze-in dark matter production from both thermal scatterings and direct inflaton decays. The authors distinguish light and heavy DM, recover several known limits, and apply the framework to gravitational DM production from SM and inflaton scatterings, presenting parameter-space maps for the observed relic density subject to CMB constraints.","tokens_in":19008,"tokens_out":19504,"duration_ms":160193,"significance":"If correct, the paper provides a useful unified treatment linking reheating dynamics parametrized by (ω, α) to UV freeze-in yields, with explicit classification of order-one, logarithmic, and power-law enhancements controlled by the critical exponent k_c = 3(ω+3)/(2α). The recovery of earlier results, the systematic mapping of viable parameter regions, and the concrete gravitational DM examples are valuable for future phenomenological studies. The analytic derivations are transparent enough for independent checking, and the inclusion of CMB bounds on T_I is a strength.","major_comments":[{"comment":"The assumed power-law background does not solve the coupled Boltzmann system (2.4)–(2.5). Substituting ρ_φ ∝ a^{-3(1+ω)} into Eq. (2.4) gives a zero left-hand side, not -Γρ_φ, so the inflaton energy density is never depleted by the decay that sources radiation in Eq. (2.5). This is not merely a formal issue: for the standard case ω=0, α=3/8, Eq. (2.6) yields Γ/H_rh = 5/2, and integrating Γ dt from a_I to a_rh gives 5/3, implying an average decay probability greater than one. Consequently, the decay yields in Eqs. (3.22) and (3.24) can exceed the total number of DM particles available from the initial inflaton population, and the normalization of all yields near a_rh inherits an uncontrolled O(1) error.","section":"Sec. 2, Eqs. (2.2)–(2.7)"},{"comment":"For the physically central case ω=0, α=3/8, the first branch of Eq. (3.22) evaluates to (4α-1)/(4α-3(1+ω)) = (1/2)/(-3/2) = -1/3, giving a negative DM yield. Since a yield cannot be negative, this signals a sign error in the derivation; the denominator should presumably be 3(1+ω)-4α on this branch. Because Section 3.2.2 and the decay-based contour plots in Figs. 4 and 6 rely on Eqs. (3.22) and (3.24), the decay phenomenology as presented is not correct and needs to be rederived.","section":"Eq. (3.22), Sec. 3.2.2"},{"comment":"The Hubble rate during reheating in Eq. (2.2) omits ρ_R, despite the definition of T_rh requiring ρ_R(a_rh) = ρ_φ(a_rh). Thus for a approaching a_rh, H is underestimated by a factor up to sqrt(2), which introduces order-one errors in the scattering yields of Eqs. (3.15) and (3.17) in the regime k < k_c where production peaks near T_rh. The paper should either solve the background equations consistently or explicitly restrict the power-law parametrization to the era with Γ ≪ H and match onto radiation domination.","section":"Sec. 2, Eq. (2.2)"}],"minor_comments":[{"comment":"The condition 'α ≤ 1' and the possible sign of the exponent (8α-3(ω+1))/2 deserve a brief discussion, since Γ(a) should remain nonnegative and finite across the parameter space shown in Fig. 1.","section":"Eq. (2.6) and surrounding text"},{"comment":"The red region labeled 'non-viable reheating' is defined by α > 3(1+ω)/4, but the caption does not state this condition; please add it for clarity.","section":"Fig. 1 caption"},{"comment":"The statement that 'the dependence in ω cancels out' in Eq. (4.2) is not immediately obvious from the displayed expression and deserves a short derivation or a clarifying comment, particularly because ρ_φ and m_φ both depend on a through ω.","section":"Sec. 4.2, Eq. (4.2)"}],"recommendation":"major_revision","confidential_remarks":"The reader's report appears to have missed the sign error in Eq. (3.22) and the background inconsistency in Section 2. The scattering analysis is likely salvageable, but the decay section as written is not reliable, and the O(1) background errors affect the quantitative maps. I recommend major revision with a request to rederive the decay yields and either solve or clearly restrict the background system."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know that this paper gives a clean analytic treatment of UV freeze-in with arbitrary power-law reheating, but the central background parametrization has a real inconsistency, and Eq. (3.22) produces a negative yield for the standard case ω=0, α=3/8. I think the reader's ACCEPT is too generous.\n\nThe parametrization of H and T with independent ω and α is genuinely useful. It recovers known limits (Refs. [104, 105]) and extends to α=0, negative α, and general decays. The separation into k below, equal to, or above k_c is neat, and the gravitational DM examples are well chosen. The derivations are transparent and the paper is clearly written.\n\nThe problem is that Eqs. (2.2)–(2.7) do not solve the Boltzmann pair (2.4)–(2.5). If ρϕ evolves as a^{-3(1+ω)}, then the left-hand side of Eq. (2.4) vanishes, but Eq. (2.6) gives nonzero Γ. Energy is not conserved unless Γ=0 or the power law already includes the decay backreaction. For ω=0, α=3/8, Γ/H_rh = 5/2, not a small correction. This is not a boundary artifact. Plugging the standard case into Eq. (3.22) gives (4α−1)/(4α−3) = (1.5−1)/(1.5−3) = −1/3, a negative abundance. That is a red flag that the integral has been mis-evaluated or the normalization mishandled. The stress-test's point is well taken: the comoving decay yield can exceed the number of inflaton particles initially available, which is unphysical. The scattering formulas near the end of reheating inherit O(1) errors when ρ_R is no longer negligible, especially for k<k_c where production peaks near T∼T_rh. So the quantitative maps in Figs. 4–6 are not secured by the stated assumptions.\n\nThat said, the paper is a useful starting point. The general framework is worth having, and the scattering formulas for k>k_c, where production happens early, may survive with modest corrections. But as written, the decay results are wrong in the standard case and the background is not self-consistent. I would recommend major revision: either solve the coupled equations or state the approximation and show it is under control in the relevant regimes. This paper is for phenomenologists computing gravitational DM or inflaton-decay freeze-in during reheating; they should wait for a corrected version.\n\nI would not desk-reject it. Send it to referees with the expectation of major revision.","headline":"Useful general framework for UV freeze-in during reheating, but the background is internally inconsistent and the decay-yield formula has a negative coefficient in the canonical case; needs major revision before the parameter-space plots can be trusted.","tokens_in":19502,"tokens_out":6112,"would_cite":false,"duration_ms":57714,"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":"Dark matter from ultraviolet freeze-in is set by two numbers that describe how the universe reheated: the equation-of-state $\\omega$ and the temperature-scaling exponent $\\alpha$.","keywords":["dark matter","ultraviolet freeze-in","FIMP","reheating","inflaton decays","gravitational production","relic abundance","reheating temperature"],"falsifier":"Take any explicit inflaton potential and decay or annihilation operator, solve the coupled background and DM Boltzmann equations with time-dependent $\\omega(a)$ and $\\alpha(a)$, and compare the final yield to the corresponding branch of Eqs. (3.17)--(3.24). If the ratio differs from one by more than an order-one factor for the same nominal $(\\omega,\\alpha)$, the constant-power-law description of reheating fails and the paper's central classification is not realized in that model.","tokens_in":18444,"feed_emoji":"🌌","tokens_out":9477,"duration_ms":81340,"temperature":0.7,"pith_summary":"The paper asks how dark matter produced by ultraviolet freeze-in—feebly interacting particles generated from the hottest part of the plasma without ever reaching equilibrium—depends on the unknown way the universe reheated after inflation. It parametrizes the reheating epoch by two power-law exponents, the inflaton equation-of-state $\\omega$ and the temperature-scaling exponent $\\alpha$, so that $H(a)\\propto a^{-3(1+\\omega)/2}$ and $T(a)=T_{\\mathrm{rh}}(a_{\\mathrm{rh}}/a)^{\\alpha}$. Within that parametrization it derives closed-form relic yields for scattering rates $\\gamma_a=T^k/\\Lambda^{k-4}$ and for direct inflaton decays. The central result is a critical exponent $k_c=\\frac{3}{2}\\frac{\\omega+3}{\\alpha}$: production during reheating is boosted by an order-one factor, a logarithm, or a power of $T_I/T_{\\mathrm{rh}}$ according to whether $k<k_c$, $k=k_c$, or $k>k_c$. The framework recovers earlier special cases and maps the viable $T_{\\mathrm{rh}}$--$T_I$ regions for gravitational dark matter production.","feed_headline":"How the universe reheated decides dark matter's abundance","feed_subtitle":"Two power-law exponents control whether freeze-in relic density is boosted logarithmically or by a power of the temperature ratio.","key_machinery":"The load-bearing object is a two-parameter model of reheating: $H(a)=H_{\\mathrm{rh}}(a_{\\mathrm{rh}}/a)^{3(1+\\omega)/2}$ and $T(a)=T_{\\mathrm{rh}}(a_{\\mathrm{rh}}/a)^{\\alpha}$ during $a_I\\le a\\le a_{\\mathrm{rh}}$, together with an effective inflaton decay width $\\Gamma(a)=4(1-\\alpha)H_{\\mathrm{rh}}(a_{\\mathrm{rh}}/a)^{(8\\alpha-3(\\omega+1))/2}$ chosen to realize that temperature law. This ansatz turns the dark-matter Boltzmann equation into a single integral that can be evaluated in closed form. The comparison point is the critical exponent $k_c=\\frac{3}{2}\\frac{\\omega+3}{\\alpha}$; it determines whether the production is dominated near $T_{\\mathrm{rh}}$ ($k<k_c$), spread logarithmically over the whole reheating interval ($k=k_c$), or dominated near $T_I$ with a power-law enhancement ($k>k_c$). For decay production the analogous dividing line is $\\alpha=\\frac{3}{4}(1+\\omega)$.","core_discovery":"The authors establish that the relic abundance from UV freeze-in during reheating is controlled by the reheating dynamics through two parameters, and they give analytic expressions for the yield at the end of reheating. For a scattering rate $\\gamma_a=T^k/\\Lambda^{k-4}$, the contribution produced during reheating is the post-reheating yield prefactor multiplied by $1/[\\alpha(k_c-k)]$, $(1/\\alpha)\\ln(T_I/T_{\\mathrm{rh}})$, or $(1/\\alpha)(T_I/T_{\\mathrm{rh}})^{k-k_c}/(k-k_c)$ depending on whether $k<k_c$, $k=k_c$, or $k>k_c$. Heavy dark matter with mass between $T_{\\mathrm{rh}}$ and $T_I$ is handled by cutting off the integral at $T=m$, yielding analogous formulas in terms of $m/T_{\\mathrm{rh}}$ or $T_I/m$. For inflaton decays the same distinction is governed by $\\alpha=\\frac{3}{4}(1+\\omega)$: the yield gains a logarithm at that value, while the would-be power-law branch lies in the non-viable reheating region. Applied to gravitational production from SM scatterings ($k=8$) and from inflaton scatterings (effective $k=6/\\alpha$), the formulas produce $(T_{\\mathrm{rh}},T_I)$ parameter-space maps for the observed relic density, with CMB constraints on $T_I$ included.","pith_inferences":["If future measurements pin down the reheating exponents $\\omega$ and $\\alpha$—for example through a gravitational-wave background—the same formulas convert that information directly into a prediction for the UV freeze-in relic density; conversely, fixing the DM abundance turns the $T_{\\mathrm{rh}}$–$T_I$ maps into a reconstruction of the reheating history.","The single-power-law assumption is the main limitation: a realistic reheating with preheating, backreaction, or time-varying exponents would require numerical integration, and the sharp $k_c$ taxonomy could become a crossover.","The same $k_c$ classification should apply to other UV-sensitive relics produced by non-renormalizable operators, such as gravitinos or axions, whenever their rates can be written in the $T^k/\\Lambda^{k-4}$ form.","A natural extension is to include momentum-dependent or threshold-suppressed production rates, where the clean power-law/log/order-one split would be smoothed out by additional scales."],"forward_implications":["For fixed coupling scale and dark-matter mass, the reheating temperature required to match $\\Omega h^2\\simeq 0.12$ drops by powers of $T_I/T_{\\mathrm{rh}}$ when $k>k_c$ and $\\alpha>0$, opening regions of much lower $T_{\\mathrm{rh}}$ than instantaneous-reheating estimates allow.","When $k=k_c$ the boost is only logarithmic, $\\ln(T_I/T_{\\mathrm{rh}})$, so $T_{\\mathrm{rh}}$ moves weakly as $T_I$ changes; when $k<k_c$ the yield is essentially independent of $T_I$, so the required $T_{\\mathrm{rh}}$ is nearly a vertical line in the ($T_{\\mathrm{rh}},T_I$) plane.","For gravitational SM production ($k=8$), the standard case $\\omega=0$, $\\alpha=3/8$ is recovered, while $\\alpha=3/4$ gives $k>k_c$ and allows $T_{\\mathrm{rh}}$ values orders of magnitude lower.","For gravitational inflaton scattering the production exponent is $k=6/\\alpha$, so the critical condition becomes a critical equation of state $\\omega_c=1$; power-law enhancement appears for both $\\omega<1$ and $\\omega>1$.","CMB constraints on the inflationary Hubble scale restrict the high-$T_I$ end of the allowed parameter space, so tensor-mode measurements can directly shrink the viable dark-matter regions."],"supporting_citations":[{"why":"Supplies the $T(a)\\propto a^{-\\alpha}$ scaling and the fact that the bath can reach $T_I>T_{\\mathrm{rh}}$ during reheating, motivating production at the hottest epoch.","marker":"[13]"},{"why":"Defines the viable reheating region $\\alpha\\le 3(1+\\omega)/4$ and the constant- or increasing-temperature scenarios that the parametrization is designed to cover.","marker":"[27]"},{"why":"Supplies the power-law forms for $H(a)$ and $T(a)$ and the effective decay width used throughout the paper.","marker":"[77]"},{"why":"Relates the equation-of-state $\\omega$ to the oscillation of the inflaton in a monomial potential, linking the abstract parameters to microphysical models.","marker":"[80]"},{"why":"Provides the field-dependent inflaton mass and decay behaviors during reheating that justify Eq. (3.7) and the decay production formulas.","marker":"[89]"},{"why":"Gives the $\\alpha$ values for bosonic and fermionic decay and annihilation channels, the basis for the scenario map in Fig. 1.","marker":"[91]"},{"why":"Earlier evaluation of DM production before reheating for $\\omega=0$, $\\alpha=3/8$, which the paper recovers as a special case of its general yield.","marker":"[104]"},{"why":"Earlier ultraviolet freeze-in analysis in non-standard cosmologies for general $\\omega$ and $\\alpha=3(1+\\omega)/8$, recovered and generalized here.","marker":"[105]"},{"why":"Source of the gravitational SM-annihilation rate coefficients $C_g$ used in Section 4.1.","marker":"[49]"},{"why":"Source of the gravitational inflaton-scattering rate density used in Section 4.2.","marker":"[47]"}],"fun_headline_variants":["Reheating dynamics decide dark matter's freeze-in abundance","Two reheating exponents govern dark matter relic density","UV freeze-in dark matter tied to reheating's power laws","Reheating era sets dark matter's ultraviolet freeze-in yield","Dark matter's freeze-in hinges on reheating expansion rate"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole argument rests on reheating being a single power-law epoch with constant $\\omega$ and $\\alpha$; if the inflaton's equation of state or the bath-temperature exponent changes with scale factor, or if reheating is non-perturbative, the analytic yields no longer apply.","fun_headline_variants_meta":{"raw":{"variants":["Reheating dynamics decide dark matter's freeze-in abundance","Two reheating exponents govern dark matter relic density","UV freeze-in dark matter tied to reheating's power laws","Reheating era sets dark matter's ultraviolet freeze-in yield","Dark matter's freeze-in hinges on reheating expansion rate"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000222,"raw_usage":{"total_tokens":1477,"prompt_tokens":990,"completion_tokens":487,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":606,"completion_tokens_details":{"reasoning_tokens":406}},"tokens_in":606,"tokens_out":487,"duration_ms":4829,"temperature":1.0,"reasoning_tokens":406,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:28:26.566219+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take any explicit inflaton potential and decay or annihilation operator, solve the coupled background and DM Boltzmann equations with time-dependent $\\omega(a)$ and $\\alpha(a)$, and compare the final yield to the corresponding branch of Eqs. (3.17)--(3.24). If the ratio differs from one by more than an order-one factor for the same nominal $(\\omega,\\alpha)$, the constant-power-law description of reheating fails and the paper's central classification is not realized in that model.","supporting_citations":[],"review_version":1}