{"id":"0f0abf1a-e3d1-44a9-a7ed-58d211010e74","arxiv_id":"2412.06626","paper_version":3,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The study connects the reheating temperature to the dark matter mass in two gravitational production scenarios and derives narrow viable mass ranges for each.","lead":"This paper derives ranges for dark matter particle masses from two ways the early universe could have reheated after inflation: by decay of heavy particles created by gravity, or by the inflaton itself decaying. It finds that the allowed dark matter mass depends strongly on which reheating path is assumed, from below the TeV scale in one scenario to about 10^11 GeV in the other.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (24), the imported production efficiency Theta_A ≈ 2.26 (m_A/Mpl)^(5/2), fixes H_END = 10^-6 Mpl and is the sole input behind every Section III mass window; its WKB derivation is not reproduced, so the quoted ranges need an independent check.","rationale":"The reader's weakest-assumption pick is the right one: Eq. (24) is the load-bearing input for Section III, and it is imported rather than derived or stress-tested in this note. I checked the delayed-decay algebra and found it internally coherent: Eq. (74) follows from Eq. (66) with z = exp(-Gamma_phi t_reh/2), and the 25/9 factor in Eq. (78) is correct, so the new delayed-decay correction is not the main risk. I also note a wording inconsistency that is not the central technical issue: the conclusion says the quintessential-inflation maximum mass is 'below the TeV scale,' but Eq. (34)'s upper end 10^-13 Mpl is about 2.4 x 10^5 GeV, several hundred TeV; this should be fixed in revision. The primary concern is that the mass windows inherit an unvalidated production formula and a fixed H_END without sensitivity quantification. That is precisely the kind of issue that conditional acceptance should require the authors to address, but it is not a demonstrated internal contradiction, so I would not move the verdict to rejection.","tokens_in":13925,"tokens_out":19964,"duration_ms":200580,"concrete_test":"Independently compute Theta_A for a conformally coupled scalar in a flat FLRW background with H_END = 10^-6 Mpl and m_A/H_END in [10^-4, 10^-1], by numerically integrating the mode equation through the end of inflation and extracting the Bogoliubov coefficient; compare the resulting Theta_A with Eq. (24) over that mass range. If the numerical values differ by more than about 30% anywhere in the range used for Eqs. (25)-(43), the Section III mass windows should be re-derived. In the same run, repeat for H_END = 3 x 10^-7 Mpl and 3 x 10^-6 Mpl to quantify the H_END sensitivity that the paper currently suppresses.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central mass windows of Section III — Eqs. (25), (34), (40), and (43) — are all rescalings of Eq. (24), which is imported from the companion paper [26]. The displayed numerical form hides an H_END dependence: the coefficient is (1/(12 pi^3)) (m_A/Mpl)^(5/2) sqrt(Mpl/(sqrt(2) H_END)), and the value 2.26 is obtained only after fixing H_END ≈ 10^-6 Mpl. The paper does not derive this formula, does not bound the WKB-in-the-complex-plane approximation, and does not state how the mass ranges would shift if H_END differs from the chosen representative value. A factor-of-10 change in H_END shifts the prefactor by about 3.2 and the derived m_Y windows by a fractional power, up to about 10^(1/5) ≈ 1.6 for the n = infinity quintessential case; the Section IV window (58) is likewise anchored to m_Y << H_END ≈ 10^-6 Mpl. The same import problem affects Eq. (55), the n = 1 production density used for the 10^11 GeV window. Because the paper's headline is a set of fairly narrow mass ranges, not a qualitative existence claim, the unvalidated status of Eq. (24) is the most load-bearing point: if that formula is off by an order of magnitude, all the Section III bounds and the claimed low-mass window shift correspondingly.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This note studies gravitational production of a conformally coupled scalar dark-matter candidate in two reheating scenarios: (i) gravitational reheating through the production and decay of heavy X-particles, and (ii) reheating via direct inflaton decay. In both scenarios the paper relates the dark-matter mass to the reheating temperature and, using the observed abundance Omega_Y h^2 = 0.12 together with BBN and gravitino bounds on the reheating temperature, quotes allowed mass windows: roughly sub-TeV for quintessential inflation (Eq. 34) and around 10^11 GeV for inflaton decay with a quadratic minimum (Eq. 58). A delayed-decay analysis is also presented, leading to a corrected reheating temperature (Eq. 81) that is the same order as the instantaneous-decay result.","tokens_in":14285,"tokens_out":17888,"duration_ms":161823,"significance":"If the imported production efficiencies (Eq. 24 and Eq. 55) are correct, the note provides a compact and useful mapping between reheating dynamics and viable dark-matter masses, with explicit falsifiable mass windows. The delayed-decay calculation in Section IV A is a genuine improvement: it treats the coupled Boltzmann system with an exponential decay factor and is internally consistent, yielding an analytic result that differs by less than an order of magnitude from earlier treatments. The main limitation is that the central mass bounds are rescalings of the imported Eq. (24), so the new results are constraints contingent on that formula rather than independent derivations.","major_comments":[{"comment":"The central production efficiency Theta_A = (1/12 pi^3)(m_A/Mpl)^{5/2} sqrt(Mpl/(sqrt(2) H_END)) ~ 2.26 (m_A/Mpl)^{5/2} is imported from the companion paper [26] without derivation. Every mass window in Section III (Eqs. 25, 27, 29, 34, 40, 43, 50, 51) is a rescaling of this formula, and the displayed numerical coefficient 2.26 fixes H_END = 10^{-6} Mpl. The paper does not quantify how the quoted ranges shift if H_END differs from this representative value or if the WKB-in-the-complex-plane approximation loses accuracy. Because this formula is load-bearing, the authors should either reproduce its essential derivation in an appendix or explicitly state the general H_END dependence and the regime of validity.","section":"Section III, Eq. (24)"},{"comment":"There is an internal numerical inconsistency in the quintessential-inflation bounds. Setting n = infinity in Eq. (27) gives m_X/Mpl >= 0.7 x 10^{-72/5} ~ 2.8 x 10^{-15}, whereas Eq. (32) states m_X/Mpl >= 5.51 x 10^{-15}, and Eq. (34) inherits the latter through Eq. (30). The intermediate algebra is not shown, so the lower end of the headline sub-TeV dark-matter window is uncertain by nearly a factor of two. Please reconcile the two expressions and display the calculation.","section":"Section III, Eqs. (27), (32), and (34)"},{"comment":"The production density in Eq. (55) is quoted for m_Y << H_END, but the allowed window in Eq. (58) extends up to m_Y/Mpl ~ 10^{-6}, which is comparable to H_END ~ 10^{-6} Mpl, and the representative value m_Y ~ 10^{-7} Mpl is only an order of magnitude below H_END. The paper should state whether Eq. (55) remains accurate at m_Y/H_END ~ 0.1 and, if not, how the upper part of the quoted mass range is modified.","section":"Section IV, Eqs. (55) and (58)"}],"minor_comments":[{"comment":"There are several typographical slips, including 'witha(t)' in the opening paragraph and the duplicated 'the' in the sentence defining rho_r(t) after Eq. (1).","section":"Throughout"},{"comment":"References [8] and [9] appear to be the same paper (identical title and arXiv number) listed twice; one duplicate should be removed or replaced with a distinct citation.","section":"References"},{"comment":"The text says 'n is a natural number' but later treats n = infinity for quintessential inflation and uses inequalities for general n; it would help to state explicitly that n may formally be taken to infinity in the kination limit.","section":"Section II"},{"comment":"The figures would benefit from axis labels and a caption explaining the fixed value of H_END used for the numerical curves.","section":"Figures 1 and 2"},{"comment":"The notation for the constraint on Theta_X^{n/(2(n-1))} is easy to misread; writing the exponent as a separate factor or using parentheses would improve clarity.","section":"Section III, Eq. (26)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a short note that relies heavily on the authors' companion papers [25,26] for the production efficiencies. The referee could not access those papers; the editor may wish to confirm that Eq. (24) and Eq. (55) appear there with the same normalization and assumptions, since the present note does not reproduce them. The internal inconsistency in the quintessential mass bounds should also be verified before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this note. First, the genuinely new part is Section IV A, the delayed inflaton decay calculation. The z_reh = 3/5 solution and the 25/9 correction to the dark matter density parameter are not in the earlier literature, and they fix a real inconsistency in Mambrini-Olive's treatment. The rest of the paper is mostly a repackaging of the authors' own companion paper [26] plus the classic Chung et al. result. Second, every mass window in Sections III and IV is a rescaling of Eq. (24), the production efficiency Θ_A ≈ 2.26 (m_A/Mpl)^(5/2), which is imported from [26] with no derivation in this note and with H_END fixed to 10^-6 Mpl. If that formula is off by a factor, all the bounds shift by roughly a fifth power. That is the load-bearing element, and the paper does not quantify the H_END sensitivity.\n\nWhat the paper does well: the derivations that are shown are coherent, the delayed-decay calculation is internally consistent, and the connection between reheating temperature and dark matter mass is clearly laid out. The figures and the explicit formulas make it easy for a model builder to plug in their own n and masses. I also found the discussion around Eq. (64) honest about why the earlier delayed-decay approach in [12] was flawed.\n\nThe soft spots are real but addressable. The main one is the unvalidated import of Eq. (24). The paper should either reproduce its derivation in an appendix or at least state explicitly how the mass ranges shift if H_END differs by a factor of 10. Second, several numerical constraints (Eqs. 23, 27, 32, 58) appear without intermediate algebra; that makes verification tedious, not impossible. Third, the authors use the observed dark matter abundance to normalize the mass-temperature relations, so the results are constraints, not parameter-free predictions. That is fine, but the abstract overstates it slightly.\n\nWho is this for? People working on gravitational production of scalar dark matter, especially in quintessential inflation or inflaton-decay reheating. It is not a new mechanism, but it sharpens the viable mass windows and corrects a small error in earlier work. I would send it to peer review—it deserves a serious referee—and I would ask the referee to check Eq. (31) and the H_END dependence carefully. I would probably cite it if I worked in this area, mostly for the delayed-decay result.","headline":"A useful little note: mostly a re-derivation of the authors' own earlier work, with one genuinely new delayed-decay correction that shifts the reheating temperature by a small factor; the mass windows all hinge on an imported production-efficiency formula that deserves independent checking.","tokens_in":14825,"tokens_out":1524,"would_cite":true,"duration_ms":18152,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.20.-q","98.80.Jk","98.80.Bp"],"model":"deepseek-v4-flash","headline":"A conformally coupled scalar produced only by gravity can be all of the dark matter only in two narrow mass windows, one light and one near $10^{11}$ GeV.","keywords":["gravitational dark matter production","gravitational reheating","reheating temperature","quintessential inflation","conformally coupled scalar","inflaton decay","WKB Bogoliubov coefficients","dark matter mass bounds"],"falsifier":"Compute the Bogoliubov coefficient numerically for a conformally coupled scalar with $m_A\\simeq10^{-12}M_{\\rm pl}$ and $H_{\\rm END}\\simeq10^{-6}M_{\\rm pl}$: if the resulting efficiency differs from $2.26(m_A/M_{\\rm pl})^{5/2}$, the mass windows shift; or find a reheating history consistent with the quoted $T_{\\rm reh}$ bounds that yields $\\Omega_Yh^2=0.12$ for a scalar mass between the two windows, which would refute the dichotomy.","tokens_in":13710,"feed_emoji":"🌌","tokens_out":17902,"duration_ms":152961,"temperature":0.7,"pith_summary":"This note sets out to show that a conformally coupled scalar produced purely by gravity—with no couplings beyond the curvature coupling—cannot have an arbitrary mass. Combining the analytic efficiency $\\Theta_A\\simeq2.26(m_A/M_{\\rm pl})^{5/2}$ with the reheating temperature bounds from Big Bang nucleosynthesis and the gravitino constraint, it derives two narrow viable windows. In gravitational reheating with a stiff post-inflationary phase the dark matter is comparatively light (for quintessential inflation, about $10^1$–$10^2$ GeV in the general case, with an upper bound of $10^{-13}M_{\\rm pl}$ in the maximum-temperature case), while in inflaton-decay reheating from a quadratic potential it must weigh around $10^{11}$ GeV. The paper also re-derives the delayed-decay reheating temperature as $T_{\\rm reh}\\simeq1.22\\times10^{-25}(M_{\\rm pl}/m_Y)^2M_{\\rm pl}$, within an order of magnitude of the instantaneous result. If these bounds hold, measuring the dark matter mass would discriminate between the two reheating histories.","feed_headline":"Gravity alone fixes dark matter mass to two narrow ranges","feed_subtitle":"The lighter window comes from gravitational reheating; the heavier one, from inflaton decay.","key_machinery":"The central object is the heating efficiency $\\Theta_A = \\rho_{A,\\rm END}/\\rho_{B,\\rm END}$ of a conformally coupled scalar, evaluated analytically by WKB Bogoliubov coefficients in the complex plane as $\\Theta_A \\simeq 2.26 (m_A/M_{\\rm pl})^{5/2}$ when $m_A\\ll H_{\\rm END}$ and $H_{\\rm END}\\simeq10^{-6}M_{\\rm pl}$. This one formula converts a scalar mass into the fraction of background energy that becomes dark matter, and all mass bounds in Sections III and IV are rescalings of it. The argument's second engine is the exact Boltzmann integral for the radiation produced by decaying X-particles or by the inflaton, Eq. (65); imposing late-time radiation conservation produces the corrected delayed-decay reheating temperature and the factor $25/9$ in Eq. (79).","core_discovery":"The paper claims that the observed dark-matter abundance $\\Omega_Y h^2=0.12$ can be reproduced by a massive scalar field that is conformally coupled to curvature and otherwise interacts only gravitationally, provided its mass lies in one of two disjoint intervals fixed by the reheating mechanism. In the gravitational-reheating branch, where the inflaton potential behaves like $\\varphi^{2n}$ near its minimum with $n\\ge3$, the maximum-temperature case gives $1.48\\times10^{-17}\\le m_Y/M_{\\rm pl}\\ll10^{-13}$, and the general decay-before-radiation case restricts quintessential inflation to roughly $10^{-17}$–$10^{-16}M_{\\rm pl}$, about $10^1$–$10^2$ GeV. In the inflaton-decay branch, where the potential is nearly quadratic, the paper obtains $2.59\\times10^{-8}\\le m_Y/M_{\\rm pl}\\ll10^{-6}$, i.e., around $10^{11}$ GeV, both for instantaneous decay and for the delayed decay it re-derives with a corrected reheating temperature $T_{\\rm reh}\\simeq1.22\\times10^{-25}(M_{\\rm pl}/m_Y)^2M_{\\rm pl}$.","pith_inferences":["A detection of dark matter with mass between the two windows would falsify the minimal conformally coupled scalar picture, independent of the reheating assumption.","Repeating the efficiency calculation for fermions, vectors, or non-conformally coupled scalars should shift the windows, so the inferred mass would become a probe of the dark sector's spin and curvature coupling.","A full numerical solution of the two-fluid Boltzmann system could check whether the corrected delayed-decay factor $25/9$ changes the predicted abundance enough to matter for CMB or BBN observables."],"forward_implications":["In gravitational reheating with a stiff post-inflationary phase (quintessential inflation), the dark matter mass is confined below $10^{-13}M_{\\rm pl}$, with the delayed-domination subcase in the $10^1$–$10^2$ GeV range.","In reheating by inflaton decay from a near-quadratic potential, the dark matter mass must be around $10^{11}$ GeV; the instantaneous and delayed decay treatments give the same order of magnitude.","The reheating temperature and the dark matter mass determine each other through inverse-square relations ($T_{\\rm reh}\\propto m_Y^{-2}$), so measuring one fixes the other within the model.","The allowed mass windows are selected entirely by the BBN lower bound and the gravitino upper bound on reheating temperature, not by any particle physics coupling.","The corrected delayed-decay reheating temperature ($T_{\\rm reh}\\simeq1.22\\times10^{-25}(M_{\\rm pl}/m_Y)^2M_{\\rm pl}$) differs by less than an order of magnitude from earlier values, so previous constraints on the model remain roughly unchanged."],"supporting_citations":[{"why":"Supplies the central production efficiency $\\Theta_A\\simeq2.26(m_A/M_{\\rm pl})^{5/2}$ and the reheating constraints used for the mass bounds.","marker":"[26]"},{"why":"Provides the gravitational reheating formulas and bounds in oscillating backgrounds that Section II reviews and extends.","marker":"[25]"},{"why":"Gives the delayed-decay integration approach that Section IV corrects to obtain the revised reheating temperature.","marker":"[12]"},{"why":"Supplies the instantaneous inflaton-decay reheating temperature $T_{\\rm reh}\\simeq5.4\\times10^{-1}\\sqrt{\\Gamma_\\varphi M_{\\rm pl}}$ used in Eq. (52).","marker":"[58]"},{"why":"Supplies Eq. (55) for the end-of-inflation energy density of gravitationally produced Y-particles in the $n=1$ case.","marker":"[4]"},{"why":"Supplies the observed values $\\Omega_Yh^2=0.12$, $\\Omega_rh^2\\simeq2.47\\times10^{-5}$, $h$, and $T_0$ used to convert ratios into mass bounds.","marker":"[24]"},{"why":"Establishes gravitational reheating through conformally coupled superheavy scalar particles, the mechanism the two-field setup builds on.","marker":"[22]"},{"why":"Gives the relation $\\Omega_Yh^2=\\Omega_rh^2\\Theta_Y T_{\\rm reh}/T_0$ used for instantaneous inflaton decay.","marker":"[2]"}],"fun_headline_variants":["Two mass windows for dark matter from gravity","Gravity alone sets dark matter mass in two ranges","Dark matter mass pinned to two ranges by reheating","Reheating branches select dark matter mass ranges","Gravitational production fixes dark matter mass intervals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The numerical bounds all trace back to the production efficiency $\\Theta_A \\simeq 2.26 (m_A/M_{\\rm pl})^{5/2}$, which presumes a conformally coupled scalar much lighter than $H_{\\rm END}$ and fixes $H_{\\rm END}$ at $10^{-6}M_{\\rm pl}$.","fun_headline_variants_meta":{"raw":{"variants":["Two mass windows for dark matter from gravity","Gravity alone sets dark matter mass in two ranges","Dark matter mass pinned to two ranges by reheating","Reheating branches select dark matter mass ranges","Gravitational production fixes dark matter mass intervals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000411,"raw_usage":{"total_tokens":2104,"prompt_tokens":899,"completion_tokens":1205,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":515,"completion_tokens_details":{"reasoning_tokens":1133}},"tokens_in":515,"tokens_out":1205,"duration_ms":11734,"temperature":1.0,"reasoning_tokens":1133,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T19:27:50.907293+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the Bogoliubov coefficient numerically for a conformally coupled scalar with $m_A\\simeq10^{-12}M_{\\rm pl}$ and $H_{\\rm END}\\simeq10^{-6}M_{\\rm pl}$: if the resulting efficiency differs from $2.26(m_A/M_{\\rm pl})^{5/2}$, the mass windows shift; or find a reheating history consistent with the quoted $T_{\\rm reh}$ bounds that yields $\\Omega_Yh^2=0.12$ for a scalar mass between the two windows, which would refute the dichotomy.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the delayed-decay integration approach that Section IV corrects to obtain the revised reheating temperature."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies Eq. (55) for the end-of-inflation energy density of gravitationally produced Y-particles in the $n=1$ case."},{"cited_title":"Gravitational production of dark matter in the Peebles-Vilenkin model","cited_arxiv_id":"1904.02393","evidence_quote":"Establishes gravitational reheating through conformally coupled superheavy scalar particles, the mechanism the two-field setup builds on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the relation $\\Omega_Yh^2=\\Omega_rh^2\\Theta_Y T_{\\rm reh}/T_0$ used for instantaneous inflaton decay."}],"review_version":1}