{"id":"bd3766e0-b14e-43d4-9b5c-b93a3d5c32a2","arxiv_id":"2505.00324","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For scalaron dark matter in f(R) gravity, the two-photon decay rate is confirmed at one loop, the resulting extragalactic background spectrum is derived, and the thermal scalaron abundance is shown to be negligible.","lead":"The paper recalculates how a scalar particle from modified gravity, the scalaron, interacts with photons and other matter, confirming a previously disputed decay rate. It also estimates the faint cosmic background this decay would create and shows that hot early-universe scalarons are a negligible dark matter fraction.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Decay rate (A.7) depends on the vanishing of scalaron-derivative couplings in the fermion loop; Appendix B's proof is heuristic and the paper itself notes the Pauli–Villars argument fails beyond one subtraction, leaving a possible finite correction to Γ.","rationale":"The reader's weakest-assumption analysis correctly identifies the field-redefinition/Jacobian issue as the most load-bearing point. The decay rate (A.7) is the quantitative core of the paper, driving both the astrophysical bounds and the comparison with earlier work. Appendix B makes a plausible case using dimensional regularization and a one-subtraction Pauli–Villars scheme, but the authors explicitly concede that the Pauli–Villars argument is limited to a single subtraction, and the two different Jacobian results (B.7)–(B.8) demonstrate sensitivity to the choice of spinor variables. This does not mean the result is wrong; the direct diagrammatic calculation may well be correct and matches the independent result of Cembranos. It does mean that the central claim is not yet settled beyond reasonable doubt. My read does not change the reader's CONDITIONAL verdict: the concern is real, but it is a technical caveat that can be resolved by an explicit two-subtraction check or by downgrading the claim of having eliminated all ambiguities. I agree with the reader's identification of the same weakest assumption.","tokens_in":15568,"tokens_out":17121,"duration_ms":183635,"concrete_test":"Re-evaluate the triangle diagrams from the first term in Eq. (B.1) without relying on momentum-shift cancellation: compute the integral (B.2) using an explicit ultraviolet cutoff with a gauge-invariant Pauli–Villars regulator implemented with two subtraction fields, and verify whether the residual finite part is exactly zero. If a non-zero surface term remains, recompute F in Eq. (A.6) including this contribution and compare the resulting Γ with the stated 5.2×10^-30 (m/MeV)^3 s^-1 at m = 2m_e.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim, Eq. (A.7), follows from the effective interaction (A.6), which is obtained by removing the scalaron from the fermion kinetic term through the field redefinition ψ = e^{3φ/4M} ψ̃ and retaining only the mass couplings (3.7)–(3.8). The load-bearing step, argued in Appendix B, is that the derivative-coupling contributions to the φγγ loop vanish: the integral (B.2) cancels in dimensional regularization, and a Pauli–Villars construction with one subtraction is used as additional support. This is exactly the point where the previous version of the paper, following ref. [6], found a non-zero Jacobian/anomaly contribution; the authors now attribute that to a regularization artifact. Their own text concedes that the Pauli–Villars argument breaks down when more than one subtraction is needed, while the direct diagrammatic argument relies on momentum-shift invariance in a quadratically divergent integral. If a finite surface term survives, L_φγγ acquires an extra term of the form (B.3) with an undetermined F_an. Near m = 2m_e, |F| ≈ 4.2, so even an O(1) F_an changes Γ by an order-unity factor, shifting all cosmological rates and bounds derived in Sec. 4. The claim that Eq. (A.7) unambiguously resolves the literature discrepancy therefore rests on a regularization-sensitive step that is asserted rather than rigorously established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the scalaron of f(R) gravity as a dark-matter candidate. It computes the one-loop scalaron decay into two photons, obtaining Γ_{φ→γγ} = α² m³ |F|²/(2^10 π³ M²) (Eq. A.7), argues that a direct loop calculation with dimensional and Pauli–Villars regularization removes the Jacobian ambiguities associated with spinor field redefinitions, derives the resulting cosmological photon background spectrum, and estimates the thermal scalaron abundance from the hot early universe, finding n_b/n_s ≲ 10^-12 (Eq. 5.18). The intended message is that the model is predictive with essentially one parameter, the scalaron mass, and that the coherent-condensate picture is self-consistent.","tokens_in":15925,"tokens_out":11065,"duration_ms":127124,"significance":"If Eq. (A.7) is accepted, the paper resolves a genuine literature discrepancy and provides a closed-form, falsifiable prediction for scalaron decay into photons, with concrete observational targets in Sec. 4. The cosmological radiation calculation is careful and transparent, and the thermal bound (5.18) is robust to many orders of magnitude of uncertainty in the cross-section estimates. The authors are also honest about the heuristic character of the regularization argument in Appendix B, which is the main risk to the central claim.","major_comments":[{"comment":"The central claim that the derivative-coupling and Jacobian contributions to the φγγ amplitude vanish is not established rigorously. The integral (B.2) is quadratically divergent in d=4, and the cancellation in dimensional regularization depends on the validity of the momentum shift, which is exactly the regulator-sensitive step. The Pauli–Villars construction uses a single subtraction, and the paper itself concedes in the last paragraph of Appendix B that the argument breaks down when more than one subtraction is needed. Since |F(2m_e)|≈4.2, an omitted O(1) F_an in Eq. (B.3) would change Γ in Eq. (A.7) by an order-unity factor and would shift all Sec. 4 background predictions, including Eq. (4.1) and the intensity spectrum (4.13). Please provide a regulator-independent proof, or an explicit calculation with at least two independent regulators demonstrating that the surface term vanishes, or alternatively quantify the residual uncertainty and propagate it into the mass bounds and background estimates.","section":"Appendix B, Eqs. (B.2)–(B.8)"},{"comment":"The cross-sections used for the thermal scalaron abundance are asserted as order-of-magnitude estimates with no derivation or error estimate, and the prefactors 10^-3 and 10^-1 are not explained. The final bound (5.18) is so small that an O(1) or even O(100) uncertainty does not affect the qualitative conclusion, so this is not a fatal issue. However, the text should state more explicitly that these prefactors are not derived from the Lagrangian and that the conclusion assumes no exponentially larger production mechanism beyond the processes considered.","section":"Section 5.1, Eqs. (5.1)–(5.3)"}],"minor_comments":[{"comment":"The statement that the model has 'only one essential free parameter' is overstated: the thermal production estimates in Sec. 5 depend on the reheating temperature T_i and on the undetermined prefactors in Eqs. (5.1)–(5.3). The sentence should be qualified to refer to the scalaron's tree-level couplings and decay rates.","section":"Section 1 and Section 6"},{"comment":"The ratio n^(b)/n^(f) is said to be 'typically much larger than unity'; it would be helpful to give the explicit numerical value for the fiducial parameters T_i=10^15 GeV, α_*=0.1, and m_f=m_e, so the reader can see the margin directly.","section":"Section 5.4, Eq. (5.14)"},{"comment":"The ratio R_Γ is plotted over a range that includes the W and top thresholds; a brief mention of where these thresholds occur would make the flatness of the ratio easier to interpret.","section":"Figure 5"},{"comment":"The denominator in the final expression is written without parentheses; adding parentheses would improve readability, e.g., 1/√(Ω_m (E_max/E)^3 + 1-Ω_m).","section":"Eq. (4.13)"}],"recommendation":"major_revision","confidential_remarks":"The main risk is the Appendix B regularization claim. If the authors cannot supply a fully rigorous proof, I would consider an explicit acknowledgment of the residual uncertainty and a sensitivity analysis of the Sec. 4 predictions to an O(1) anomaly coefficient sufficient for publication. The self-citation to the earlier preprint [35] is handled transparently and is not a concern."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid, honest paper. The decay rate itself is not new — it confirms Cembranos — but the analytic extragalactic background spectrum (Eq. 4.13) and the estimate that thermal scalarons are negligible (Eq. 5.18) are genuinely useful additions. The loop calculation is careful, and the authors deserve credit for explicitly reversing their earlier preprint position and for flagging the limitations of their own argument.\n\nWhere I would push: the whole edifice rests on Appendix B, where the authors argue that derivative-coupling contributions from the spinor kinetic term vanish and no Jacobian anomaly survives. The dimensional-regularization argument involves a shift in a quadratically divergent integral; the Pauli–Villars construction uses one subtraction and the authors themselves note it breaks down beyond that. If a surface term survives, the effective interaction acquires an undetermined piece of the form (B.3). Near m = 2m_e, |F| ≈ 4.2, so an O(1) F_an changes the decay rate by an order-unity factor and all the Section 4 rates and bounds shift correspondingly. So the paper is not wrong in an obvious way, but its central claim to have resolved the literature discrepancy is asserted more strongly than the proof warrants.\n\nThe thermal section is order-of-magnitude, as the authors admit; that is fine because the conclusion — thermal scalarons are a negligible fraction of the dark matter — is robust to the large uncertainties. The self-citations [15,16] supply initial-condition premises and do not prop up the loop calculation, so I do not see a citation problem.\n\nWho this is for: anyone working on scalaron or modified-gravity dark matter, and on line or continuum photon signals from f(R) models. A referee should spend time on Appendix B and ask for a more rigorous treatment of the Jacobian — or an explicit statement that the result is regularization-dependent within the allowed renormalization freedom. If that step is fixed or honestly downgraded, the paper is publishable; the cosmological spectrum and thermal bound are useful additions.","headline":"Careful re-derivation of the scalaron-to-two-photon rate with useful new cosmological applications, but the central claim to have eliminated the Jacobian ambiguity is not airtight.","tokens_in":16427,"tokens_out":9687,"would_cite":false,"duration_ms":114186,"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":"In f(R) gravity the scalaron can be all of dark matter; its one-loop decay to two photons has a definite rate and thermal production is negligible.","keywords":["f(R) gravity","scalaron","dark matter","one-loop decay","two-photon decay","photon background","thermal production","modified gravity"],"falsifier":"Compute the coefficient of $\\phi F_{\\mu\\nu}F^{\\mu\\nu}$ in the one-loop effective action directly from the fermion determinant in a background $\\phi$ and electromagnetic field, without any field redefinition; if the coefficient differs from $(\\alpha/16\\pi)F(m)/M$, the rate formula is wrong. Observationally, a dark-matter halo should show a two-photon line at $E_\\gamma=m/2$ with the intensity and spectrum derived here, and a meaningful upper limit below that prediction would exclude the scenario once halo modelling is included.","tokens_in":15360,"feed_emoji":"🌌","tokens_out":13827,"duration_ms":129984,"temperature":0.7,"pith_summary":"The scalaron of $f(R)$ gravity—the scalar degree of freedom hidden in the modified metric—can account for all of the dark matter if its mass lies between a few meV and about 1 MeV. This paper re-derives the scalaron's interaction with the Standard Model and settles a disputed one-loop result: its decay to two photons has the definite rate $\\Gamma_{\\phi\\to\\gamma\\gamma}=\\alpha^2 m^3 |F|^2/(2^{10}\\pi^3 M^2)$, with $\\Gamma\\approx 5.2\\times 10^{-30}(m/\\mathrm{MeV})^3\\,\\mathrm{s}^{-1}$ near $m=2m_e$. The paper computes the cosmological photon background this decay would create and shows that scalarons produced thermally in the early universe contribute at most one part in $10^{12}$ of the dark-matter density. These results support the original picture in which scalaron dark matter is a coherently oscillating classical field rather than a thermal relic.","feed_headline":"Scalaron dark matter's photon-decay rate is now unambiguous","feed_subtitle":"The two-photon decay rate sets the background spectrum; thermal scalarons are negligible.","key_machinery":"The load-bearing object is the one-loop effective interaction $$\\mathcal{L}_{\\phi\\gamma\\gamma}=\\frac{\\$\\alpha$}{16\\pi}F(m)\\frac{\\phi}{M}F_{\\mu\\nu}$F^{{\\mu\\nu}}$,$$ where $F(m)$ is the same Standard Model loop form factor that governs $H\\to\\gamma\\gamma$ ($W$-boson and charged-fermion loops). It carries the argument because the scalaron enters each diagram at a single vertex with the replacement $\\chi/v\\to -\\phi/2M$, so the entire disputed question reduces to whether the spinor field redefinition used to diagonalise the kinetic term generates extra contributions. The paper argues it does not, at one loop: redefinition-dependent kinetic terms are proportional to the free spinor equation of motion and vanish inside the triangle diagram, while the Jacobian is a renormalisation-scheme artefact; a ghost-field regularisation shows the Jacobians cancel between fields of opposite statistics. Only the mass couplings survive, making the amplitude unambiguous.","core_discovery":"The paper's central claim is that the ambiguity in the scalaron–photon vertex is an artifact of regularisation-dependent Jacobians and disappears when the one-loop amplitude is computed directly. Because the scalaron couplings to fermion mass terms and to $W^\\pm$ and $Z^0$ boson masses have the same form as the Higgs couplings with $\\chi/v\\to -\\phi/2M$, the Standard Model $H\\to\\gamma\\gamma$ form factor carries over. The spinor field redefinition needed to remove $\\phi$ from kinetic terms produces no net contribution: the kinetic part of the redefined Lagrangian is proportional to the free spinor equation of motion, so it cannot contribute to the triangle diagram, and the Jacobian is a gauge-choice artifact rather than a physical term. The resulting decay width, $\\Gamma_{\\phi\\to\\gamma\\gamma}=\\alpha^2 m^3 |F|^2/(2^{10}\\pi^3 M^2)$, is real for $m<2m_e$ and gives $\\Gamma\\approx 5.2\\times10^{-30}(m/\\mathrm{MeV})^3\\,\\mathrm{s}^{-1}$ near threshold; with this rate, the diffuse background from cosmological decays has a calculable spectrum and the thermal scalaron fraction is $n_b/n_s\\lesssim 10^{-12}$.","pith_inferences":["The paper does not work out the corresponding gluon final state, but the same one-loop matching applies directly: scalaron decay into two gluons should have an equally definite rate built from quark-loop contributions, which would extend the model's observability into cosmic-ray and gamma-ray channels.","The Jacobian-cancellation claim is a statement about the fermion determinant in the background $\\phi F_{\\mu\\nu}F^{\\mu\\nu}$; a regulator-independent evaluation of that determinant would separate the physical amplitude from scheme-dependent artifacts and is a natural independent check.","Because the decay rate is fixed once $m$ is known and the mass window is narrow, the model is directly falsifiable by a line search at $E_\\gamma=m/2$; a detection would measure $m$ and test the minimal-coupling assumption at the same time."],"forward_implications":["If the scalaron is all the dark matter, its decays produce a diffuse photon background with maximum energy $E_{\\rm max}=m/2$ and a spectrum that rises toward that maximum; near the upper allowed mass bound the line intensity scales as $(m/\\mathrm{MeV})^3$, making that regime the most visible.","The 511-keV positron-annihilation bound on scalaron decay to $e^+e^-$ confines a dark-matter scalaron to roughly $1.04\\,\\mathrm{MeV}\\lesssim m\\lesssim 1.15\\,\\mathrm{MeV}$ if it makes all of the dark matter.","Thermal scalarons from the hot plasma contribute at most $10^{-12}$ of the dark-matter number density, so the dark matter is a coherent condensate rather than a thermal relic in this model.","The rate is real and unambiguous for $m<2m_e$, so the formula stays valid across the whole allowed mass window and can be used as a prediction for line searches."],"supporting_citations":[{"why":"The original proposal that the $R^2$-gravity scalaron is all dark matter; its decay estimate is the one this paper confirms and refines.","marker":"[3]"},{"why":"Supplies the 511-keV positron bound and mass-window argument used to constrain a dark-matter scalaron.","marker":"[4]"},{"why":"An earlier Jacobian-based calculation of scalaron interactions whose ambiguity the direct loop calculation is meant to remove.","marker":"[6]"},{"why":"Sets the initial-conditions framework for the scalaron condensate assumed by the dark-matter scenario.","marker":"[15]"},{"why":"Provides the scalaron thermal history and the perturbative-decay assumption that avoids parametric resonance.","marker":"[16]"},{"why":"Supplies current cosmological parameters used in the numerical estimates of the photon background and thermal fraction.","marker":"[34]"},{"why":"The earlier version of this paper whose Jacobian-based result is superseded by the direct loop calculation.","marker":"[35]"},{"why":"Provides the Standard Model $H\\to\\gamma\\gamma$ amplitude whose form factor is transferred to the scalaron.","marker":"[37]"},{"why":"Justifies replacing the Higgs vertex by the scalaron vertex in the one-loop diagrams.","marker":"[38]"},{"why":"Supplies the path-integral and anomaly formalism whose regularisation dependence Appendix B examines.","marker":"[39]"}],"fun_headline_variants":["Scalaron dark matter's photon decay rate: ambiguity resolved","Two-photon decay of scalaron dark matter: no more ambiguity","Scalaron dark matter: photon decay rate calculation unambiguous","Scalaron's two-photon decay rate: unambiguous at last"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes the spinor field redefinition $\\psi=e^{3\\phi/4M}\\tilde{\\psi}$ leaves the S-matrix unchanged at one loop, so the Jacobian and kinetic-term contributions to $\\phi\\gamma\\gamma$ vanish and only the mass couplings matter; if a Jacobian term survived, the decay rate and background predictions would shift.","fun_headline_variants_meta":{"raw":{"variants":["Scalaron dark matter's photon decay rate: ambiguity resolved","Two-photon decay of scalaron dark matter: no more ambiguity","Scalaron dark matter: photon decay rate calculation unambiguous","Scalaron's two-photon decay rate: unambiguous at last"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000483,"raw_usage":{"total_tokens":2412,"prompt_tokens":996,"completion_tokens":1416,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":612,"completion_tokens_details":{"reasoning_tokens":1347}},"tokens_in":612,"tokens_out":1416,"duration_ms":9803,"temperature":1.0,"reasoning_tokens":1347,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:46:21.175069+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the coefficient of $\\phi F_{\\mu\\nu}F^{\\mu\\nu}$ in the one-loop effective action directly from the fermion determinant in a background $\\phi$ and electromagnetic field, without any field redefinition; if the coefficient differs from $(\\alpha/16\\pi)F(m)/M$, the rate formula is wrong. Observationally, a dark-matter halo should show a two-photon line at $E_\\gamma=m/2$ with the intensity and spectrum derived here, and a meaningful upper limit below that prediction would exclude the scenario once halo modelling is included.","supporting_citations":[],"review_version":1}