{"id":"ea08dddb-b44a-4530-8735-267f5b529907","arxiv_id":"2505.05206","paper_version":2,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A full thermal-resummation computation shows that plasmon and plasmino corrections leave the TeV-scale singlet scalar annihilation cross section essentially unchanged, reconfirming existing relic density constraints.","lead":"This paper computes how hot plasma effects, extra polarization states of gauge bosons and fermions, alter the annihilation rate of a heavy singlet scalar dark matter candidate. It finds these effects are tiny for TeV scale masses, so standard freeze-out calculations and the resulting dark matter constraints stand unchanged.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unquantified pole truncation of HTL spectral functions could hide an O(m_D^2/m_phi^2 log) continuum; a direct spectral integral would settle whether the temperature-independence claim survives.","rationale":"The reader identified the pole truncation of HTL spectral functions as the weakest assumption; I agree. My work here is to sharpen it: the continuum spectral weight is parametrically m_E^2/q^2 log, not suppressed by a Boltzmann factor, and the asymptotic expansions used for the poles do not control the continuum at q0-q ~ q. This is a genuine soft spot in the proof. However, the numerical estimate suggests the effect is at the 10^-3 level for the TeV benchmark, and the paper's final relic-density curve agrees with earlier unresummed computations, which would be difficult to maintain if the continuum were O(1). The missing piece is a direct numerical check, not a new conceptual ingredient. I therefore do not change the reader's ACCEPT verdict: the concern is worth testing but not, on the current evidence, enough to reject or make acceptance conditional.","tokens_in":30978,"tokens_out":18788,"duration_ms":212359,"concrete_test":"Evaluate the WW contribution to Eq. (4.18) for the benchmark m_phi,phys = 2 TeV, kappa = 1, at T = 100 and 160 GeV, using the exact HTL spectral functions rho_T,E = -Im G_T,E from Eqs. (3.10)-(3.11) and (3.16)-(3.18), with the q0 and r0 integrals performed numerically, and compare with the same integral using the delta-function approximations (4.30)-(4.31). The fractional difference measures the pole-truncation error directly. Repeat for the ZZ and tt channels (Eqs. (3.56)-(3.57)). If the difference is below 1%, the central claim is unchanged; if it is larger, the temperature-independence statement requires revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim (Sec. 4.4 and Conclusions) rests on replacing the full HTL spectral functions of final-state W/Z/t by on-shell delta functions at their asymptotic masses, Eqs. (4.30)-(4.31), (C.8), (D.6)-(D.8). The dropped part is the Landau-damping continuum |q0|<q. At hard q ~ m_phi its integrated weight is not exponentially small: for q0 = q - delta, Im Π_T ~ m_E^2 delta/(2q), while |G_T^{-1}|^2 ~ (2q delta + m_WT^2)^2, so ∫ rho_cont dq0 ~ (m_E^2/q^2) log(q^2/m_WT^2). For the benchmark m_phi,phys = 2 TeV, T ~ 160 GeV this is about 10^-3, so the effect is likely harmless; but it is never computed. The asymptotic expansions (4.26)-(4.27) used to justify pole dominance are valid only for q0 - q ≪ q, so the continuum at q0 ~ q/2 is outside the stated domain of validity and is not covered by the paper's power counting. Because the continuum enters the same 1->2 phase-space integral (4.18) as the pole terms, it could be as large as the retained (b)-(d) corrections. A direct numerical evaluation of the spectral integral is needed to show it is not.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper computes the thermally averaged annihilation cross section for a TeV-scale singlet scalar dark matter candidate with Hard Thermal Loop (HTL) resummed Standard Model propagators. It collects the HTL-resummed propagators for the W, Z, and top quark in the electroweak crossover region, including plasmon and plasmino branches, and uses them to evaluate the imaginary part of the singlet scalar self-energy. The central derivation is carried out in Sections 3 and 4, where the authors show gauge-parameter cancellation, derive power-counting estimates for the various channels, and present numerical results for the WW, ZZ, and t-tbar contributions. The resulting relic density is then computed in Section 5 and found to agree with earlier unresummed calculations. The main claim is that for m_phi,phys in the few-TeV range, the inclusive annihilation cross section is dominated by a temperature-independent longitudinal gauge channel, with thermal corrections suppressed as powers of the Debye mass over the dark matter mass, so cosmological constraints are unchanged.","tokens_in":31268,"tokens_out":16650,"duration_ms":174980,"significance":"If the result holds, the paper is a valuable technical reference. It provides a complete set of HTL-resummed Standard Model Feynman rules and propagators in the crossover region, demonstrates explicitly that the gauge dependence cancels in the annihilation rate, and gives transparent power-counting estimates that are numerically confirmed in Fig. 2. The relic-density curve is an output of the calculation, not an input, and the comparison with earlier unresummed results is a useful crosscheck. The main weakness is that the replacement of full HTL spectral functions by on-shell poles is not quantitatively controlled in the paper, as discussed in the major comments.","major_comments":[{"comment":"The central quantitative claim is established only after replacing the full HTL spectral functions by delta functions at their asymptotic masses (Eqs. (4.30)-(4.31), and similarly for ZZ and t-tbar in Appendices C-D). The Landau-damping continuum |q0|<q is dropped without an estimate. The asymptotic expansions (4.26)-(4.27) are derived under the condition q0-q << q and do not control the continuum, which enters the same 1->2 phase-space integral (4.18) through the spectral densities. A crude estimate for the benchmark m_phi,phys = 2 TeV, T = 160 GeV gives a continuum weight of order (m_E2^2/q^2) log(qT/m_W^2) ~ 10^-2, which is small but orders of magnitude larger than the retained (b)-(d) channels in Fig. 2. Since the manuscript claims in Sec. 6 that the computation takes into account the \"full structures predicted by HTL effective theories,\" the authors should either evaluate the continuum contribution numerically with the full spectral densities or explicitly state and justify the kinematic domain in which it is negligible.","section":"Sec. 4.4, Eqs. (4.23)-(4.31) and Fig. 2"},{"comment":"The power-counting estimates in Eqs. (4.41)-(4.44) assume q0 ≈ q for both final-state particles. However, in the lab frame the decay products carry momenta that can differ from m_phi by up to O(sqrt(m_phi T)) (compare Eq. (4.22)), so contributions with q0 - q of order T or larger are kinematically allowed and are not exponentially suppressed by the Boltzmann factor e^{-beta p0}. This means the statement in Eq. (4.23) that q0 - q ~ m_infty^2/m_phi is too restrictive for the purpose of bounding the spectral continuum. The conclusion that thermal corrections are power-suppressed needs an explicit estimate over the full integration region, not only the near-pole region.","section":"Sec. 4.2, Eq. (4.18) and Sec. 4.4, power counting"}],"minor_comments":[{"comment":"The symbol m_h is used for different quantities in the two phases: the physical Higgs-like excitation below T_c and the resummed scalar mass parameter above T_c. Renaming the high-temperature scalar mass, for example to m_phi, would avoid confusion.","section":"Sec. 3.1, Eqs. (3.5)-(3.6)"},{"comment":"The curve labelled m_ZT- is the thermal asymptotic mass of the hypercharge/photon-like excitation; the caption could state this more explicitly, since it is not obvious from the figure alone.","section":"Fig. 1"},{"comment":"The agreement with refs. [21, 23] is stated as being \"within plot resolution\" and the quoted benchmark for m_phi,phys = 9.79 TeV lies mostly in the T > T_c regime, where the new HTL machinery is not exercised. A direct quantitative comparison for a benchmark with freeze-out in the Higgs phase, such as m_phi,phys = 2 TeV, would strengthen the claim.","section":"Sec. 5, Fig. 3"},{"comment":"The replacement of the averaged squared matrix element Phi_hh by the expression in Eq. (B.4) is non-trivial and the derivation is only sketched. A short explanation or a more explicit reference to the corresponding step in ref. [20] would improve readability.","section":"Appendix B, Eq. (B.4)"}],"recommendation":"major_revision","confidential_remarks":"The concern about the unquantified spectral continuum is the main obstacle to acceptance. I do not think it is a fatal flaw: the final relic-density constraints are likely unchanged, and the authors already have the full spectral functions in hand, so a direct numerical evaluation of the continuum contribution appears feasible within the scope of the manuscript. If that numerical estimate confirms the smallness of the continuum, the paper would be acceptable after a revision. The comparison with the literature is somewhat weakened by the fact that the quoted high-mass benchmark freezes out mostly above T_c; a Higgs-phase benchmark comparison would be more compelling."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe thing to know: this is a serious, well-executed computation that actually answers a question that had been begged in the singlet scalar literature — whether Hard Thermal Loop resummation changes the freeze-out of a TeV-scale scalar near the electroweak crossover. It does not. The authors show this carefully, with real crosschecks, and the final relic-density constraints are unchanged from the unresummed results of GAMBIT and Cline et al.\n\nWhat is genuinely new is the full set of HTL-resummed Standard Model propagators in the Higgs phase, including the plasmon/plasmino branches and mixings, and the explicit cut computations for the WW, ZZ and top-antitop channels. The gauge-parameter cancellation in eq. (3.54) is a strong consistency check, and the power-counting estimates in eqs. (4.41)–(4.44) are confirmed by the numerics. The appendices are detailed enough that a patient reader can redo the algebra. The agreement with the earlier relic-density results is an independent check, not a restatement, because the technical path here is different.\n\nThe main soft spot is the treatment of the final-state HTL spectral functions at hard momenta. The paper replaces them by on-shell poles at the asymptotic masses and drops the Landau-damping continuum with |q0|<q. The asymptotic expansions used to justify pole dominance are only valid for q0−q much less than q, so the continuum at q0 about q/2 sits outside the stated domain of validity. The stress-test estimate puts the continuum weight at roughly 10^-3 of the leading term for the 2 TeV benchmark, so the central conclusion very likely survives. But because the authors already integrate numerically, it would have been straightforward to keep the full spectral functions and remove the doubt. That omission is a minor flaw, not a load-bearing one. Also minor: no code or data release, although the appendices partially compensate.\n\nThe paper is honest about its limits: it explicitly says the approximation breaks down below about 0.5 TeV, where thermal effects become O(1), and it does not oversell the phenomenological implications. The citation pattern is fine; the self-citations point to directly relevant earlier work.\n\nWho should read it: practitioners computing dark matter freeze-out around the electroweak crossover, and anyone who wants the HTL-resummed SM propagator catalog for other problems. I would send it to peer review and accept after the authors add a direct estimate of the continuum contribution. The paper is worth referee time.","headline":"A solid, carefully crosschecked HTL computation that reconfirms the standard singlet scalar limits; the dropped spectral-function continuum is a real but likely harmless gap.","tokens_in":31826,"tokens_out":3453,"would_cite":true,"duration_ms":36194,"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":"For a TeV-scale singlet scalar dark matter candidate, fully resummed hot-plasma final states leave the inclusive annihilation cross section essentially at its vacuum value, so the relic-density constraint on the model is unchanged.","keywords":["singlet scalar dark matter","Hard Thermal Loop resummation","plasmon","plasmino","dark matter freeze-out","electroweak crossover","relic density","thermal field theory"],"falsifier":"Evaluate the full HTL spectral integrals in eq. (4.18) without replacing the spectral functions by on-shell poles, at $m_\\phi\\simeq2$ TeV and $T\\simeq160$ GeV, and compare the result with the pole approximation; if the continuum weight contributes at order $(m_E/m_\\phi)^2$ rather than being negligible at leading power, the paper's power counting would need revision.","tokens_in":30751,"feed_emoji":"🌌","tokens_out":10267,"duration_ms":97737,"temperature":0.7,"pith_summary":"This paper asks whether the hot plasma at the electroweak crossover can change the annihilation cross section that sets the relic density of a heavy singlet-scalar dark matter particle. Final-state Standard Model particles in a plasma are not vacuum particles: Hard Thermal Loop resummation gives gauge bosons a third polarization state (plasmon) and fermions an extra branch (plasmino). The paper assembles the full resummed spectrum for the Standard Model near $T_c\\approx160$ GeV and computes the inclusive $2\\to2$ annihilation cross section for a TeV-scale singlet scalar. It finds that the dominant longitudinal gauge channel is essentially temperature-independent, that thermal corrections enter only as powers of the Debye mass over the dark matter mass, and that the resulting relic density curve agrees with previous unresummed computations. The upshot is that cosmological constraints on TeV-scale singlet scalars are reconfirmed.","feed_headline":"Dark matter freeze-out unchanged by hot plasma quasiparticles","feed_subtitle":"Resummed annihilation stays near its vacuum value, so relic-density constraints on TeV singlet scalars stand.","key_machinery":"The central machinery is the HTL-resummed propagators together with the cut of the retarded singlet-scalar self-energy. The cross section is obtained from $\\Gamma^{\\mathrm{max}}_{K;\\phi}=\\mathrm{Im}\\,\\Pi_{K;\\phi}/\\omega$ and a two-particle phase-space representation in which a heavy off-shell Higgs of momentum $P$ decays into two Standard Model quasiparticles. At hard momenta $q\\sim m_\\phi\\gg\\pi T$, each HTL propagator reduces to pole-like spectral functions with asymptotic masses: the transverse $W$ acquires $m_{WT}^2=m_W^2+m_{E2}^2/2$, while the electric (longitudinal) $W$ keeps $m_W^2$ with no thermal shift; the $Z$ splits into three channels, and the top quark into two chiral channels. Power counting with $q\\sim m_\\phi$, $p\\sim\\sqrt{m_\\phi T}$, and $p_0\\sim2m_\\phi$ then shows that the longitudinal gauge channel dominates, with thermal corrections suppressed by powers of $m_E^2/q^2$.","core_discovery":"The paper's central claim is that in the freeze-out of a TeV-scale singlet scalar with mass $m_{\\phi,\\mathrm{phys}}\\sim$ few TeV and coupling $\\kappa\\sim1$, the thermally averaged inclusive annihilation cross section $\\langle\\sigma v_{\\mathrm{rel}}\\rangle$ computed with the full Hard-Thermal-Loop-resummed Standard Model is dominated by the longitudinal gauge channel, whose contribution is independent of both the temperature and the Higgs mechanism. The largest correction to that leading term is the ordinary vacuum correction induced by the Higgs mechanism, and genuine thermal effects appear only as power corrections: they are proportional to the second or fourth power of the Debye mass, i.e. suppressed by $(m_E/m_\\phi)^2$ or $(m_E/m_\\phi)^4$, with $m_E\\sim gT$. Because these corrections are numerically negligible for $m_\\phi\\sim2$ TeV, the relic-density contours in the $(m_{\\phi,\\mathrm{phys}},\\kappa)$ plane agree with earlier computations that did not implement HTL resummation, reproducing the benchmark $\\Omega_\\phi h^2\\approx0.114$ for $m_{\\phi,\\mathrm{phys}}=9.79$ TeV and $\\kappa=3.1$.","pith_inferences":["We infer that the dominance of the temperature-independent longitudinal channel follows from kinematics rather than from the scalar's quantum numbers, so other heavy Higgs-portal dark matter candidates should show the same insensitivity of the inclusive annihilation rate to HTL corrections.","We infer that the paper's pole approximation predicts the continuum part of the HTL spectral functions at $q\\sim m_\\phi$ to be negligible at leading power; a direct numerical check of that weight would bound the neglected corrections.","We infer that once the relic-density constraint is confirmed to be resummation-insensitive, the dominant theoretical uncertainty in the constraint shifts to the inputs $v(T)$ near the crossover and the Standard Model equation of state, rather than to the thermal-resummation procedure."],"forward_implications":["The relic-density contour in the $(m_{\\phi,\\mathrm{phys}},\\kappa)$ plane is the same one obtained without HTL resummation, including the benchmark $\\Omega_\\phi h^2\\approx0.114$ at $m_{\\phi,\\mathrm{phys}}=9.79$ TeV and $\\kappa=3.1$.","Genuine thermal corrections to the inclusive cross section scale as $(m_E/m_\\phi)^2$ or $(m_E/m_\\phi)^4$, so they are numerically negligible for TeV-scale masses and do not shift the cosmological constraint.","The exclusive channels individually acquire new plasmon and plasmino contributions, but their sum remains close to the vacuum-like result, so an exclusive view is viable even though the intermediate steps are substantially modified.","Below $m_{\\phi,\\mathrm{phys}}\\sim0.5$ TeV the hierarchy $m_\\phi\\gg\\pi T$ breaks down, thermal masses can close the $2\\to2$ channel, and a computation including initial equilibration and $1\\leftrightarrow2$ processes would be required."],"supporting_citations":[{"why":"Supplies the Hard Thermal Loop effective theory and resummed propagators that define the plasmon and plasmino states; the starting point for the computation.","marker":"[13–16]"},{"why":"Provides the Lorentz-invariant averaging of initial-state momenta used in appendix B for the Higgs channels.","marker":"[20]"},{"why":"Provides an earlier scalar-singlet relic-density computation that the paper's figure 3 is compared with.","marker":"[21]"},{"why":"Provides the comparison benchmark with $\\Omega_\\phi h^2\\approx0.114$ at $m_{\\phi,\\mathrm{phys}}=9.79$ TeV and $\\kappa=3.1$.","marker":"[23]"},{"why":"Defines the maximal interaction rate used to connect the self-energy cut to the annihilation cross section.","marker":"[36]"},{"why":"Provides the lattice determination of the electroweak crossover temperature used to parametrize $v(T)$ and locate the freeze-out regime.","marker":"[39]"},{"why":"Provides the HTL-resummed fermion propagator formalism that yields the plasmino branch and the top and bottom spectral functions.","marker":"[46]"},{"why":"Provides the Standard Model equation-of-state tabulations used in the freeze-out integration.","marker":"[48]"},{"why":"Provides multiloop running couplings and masses used for the Standard Model inputs.","marker":"[49]"}],"fun_headline_variants":["Thermal corrections vanish in TeV singlet freeze-out","Plasmon and plasmino states leave freeze-out unchanged","Hot plasma quasiparticles invisible in dark matter relic density","Freeze-out survives Hard Thermal Loop resummation intact","TeV-scale singlet relic density confirmed despite thermal effects"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The computation assumes that, at the high momenta typical of the annihilation products, each hot-plasma particle behaves like an ordinary particle with a shifted mass, so the continuous part of its spectral function can be neglected; if that continuous part carried non-negligible weight at momenta of order $m_\\phi$, the suppression of thermal corrections could fail.","fun_headline_variants_meta":{"raw":{"variants":["Thermal corrections vanish in TeV singlet freeze-out","Plasmon and plasmino states leave freeze-out unchanged","Hot plasma quasiparticles invisible in dark matter relic density","Freeze-out survives Hard Thermal Loop resummation intact","TeV-scale singlet relic density confirmed despite thermal effects"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000201,"raw_usage":{"total_tokens":1391,"prompt_tokens":969,"completion_tokens":422,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":585,"completion_tokens_details":{"reasoning_tokens":341}},"tokens_in":585,"tokens_out":422,"duration_ms":4593,"temperature":1.0,"reasoning_tokens":341,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:10:48.219104+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the full HTL spectral integrals in eq. (4.18) without replacing the spectral functions by on-shell poles, at $m_\\phi\\simeq2$ TeV and $T\\simeq160$ GeV, and compare the result with the pole approximation; if the continuum weight contributes at order $(m_E/m_\\phi)^2$ rather than being negligible at leading power, the paper's power counting would need revision.","supporting_citations":[{"cited_title":"The real singlet scalar dark matter model","cited_arxiv_id":"1006.2518","evidence_quote":"Provides the Lorentz-invariant averaging of initial-state momenta used in appendix B for the Higgs channels."},{"cited_title":"Resonant $s$-channel dark matter annihilation at NLO","cited_arxiv_id":"2211.06008","evidence_quote":"Defines the maximal interaction rate used to connect the self-energy cut to the annihilation cross section."},{"cited_title":"Weldon, Effective fermion masses of order gT in high-temperature gauge theories with exact chiral invariance, Phys","cited_arxiv_id":null,"evidence_quote":"Provides the HTL-resummed fermion propagator formalism that yields the plasmino branch and the top and bottom spectral functions."},{"cited_title":"Laine, M","cited_arxiv_id":null,"evidence_quote":"Provides the Standard Model equation-of-state tabulations used in the freeze-out integration."}],"review_version":1}