{"id":"41a81f73-3ebc-4e75-b7b4-b6be9ef6a20c","arxiv_id":"2505.02408","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"At next-to-leading order, the real scalar singlet dark matter model is viable only for masses above about 20 TeV, and the complex scalar version is excluded in its perturbative regime.","lead":"This paper calculates higher-order quantum corrections for the simplest dark matter model, a scalar particle interacting only with the Higgs boson. Combined with the latest LZ direct-detection limits, it finds the real version survives only above about 20 TeV, while the complex version is essentially or fully excluded.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'entire perturbative region excluded' claim for the complex scalar rests on the β-function 25% criterion λφH<15; with the tree-level unitarity bound λφH<8π the paper's own Table I leaves a viable window, so the abstract overstates the model verdict.","rationale":"In good faith, the paper accomplishes a careful NLO calculation: the real-scalar lower bound is supported by three renormalization schemes spanning only 22–24 TeV, the scheme-independent combination Δ_FO − Δ_DD explains this stability, and the analytic formulae in Appendix A are explicit enough to be checked. The paper is also transparent, quoting both the β-function criterion (λ ≈ 15) and the unitarity bound (λ ≈ 8π) and flagging the scheme sensitivity in the complex case. The load-bearing weakness is the definition of the perturbative region itself. The central claims are not internally inconsistent; they are convention-dependent. Since the reader identified precisely this premise and already returned CONDITIONAL, my stress test confirms the reader's weakest-assumption analysis rather than moving the verdict. The recommended revision is to re-state the complex-scalar exclusion and the real-scalar upper window as functions of the perturbativity cutoff, e.g. 'excluded for λφH < 15; a narrow window reappears if the unitarity limit λφH < 8π is adopted.'","tokens_in":9579,"tokens_out":11864,"duration_ms":155589,"concrete_test":"Recompute the Section IV/Table I perturbativity ceilings for the complex and real scalars with the cutoff varied over λφH ∈ {15, 20, 25}, using the same NLO relic-density and direct-detection rates. For each renormalization scheme (MSbar, OS-DD, OS-FO), compare the maximum mφ consistent with the relic density at the chosen ceiling with the NLO direct-detection lower bound in that scheme. If any scheme shows mφ(DD) < mφ(pert) for λ = 25, the abstract's 'entire perturbative region excluded' claim is falsified by the authors' own rates; if no scheme shows a window even at the unitarity limit, the exclusion survives but should be re-stated as valid up to λ = 8π. No new diagram computation is needed: only the relic-density curves already used in Figure 1 need to be re-evaluated at the alternative ceilings.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The complex-scalar exclusion and the upper edge of the real-scalar window are both controlled by the perturbativity cutoff defined in Section IV. Equation (9), λφH < 32π²/21 ≈ 15, is obtained by demanding that the two-loop NLL term in the β function be no more than 25% of the one-loop term. This is a defensible literature-based prescription, but it is a convention, not a rigorous bound, and the same section quotes the tree-level unitarity bound |λφH| < 8π ≈ 25. The abstract's claim that the entire perturbative region is excluded for the complex scalar is directly sensitive to this choice: Table I already shows narrow allowed windows in the MSbar (26 TeV direct-detection bound vs. 29 TeV perturbativity ceiling) and OS-FO (32 vs. 33 TeV) schemes at λ = 15, with the text describing these cases as 'above or very close'. If the perturbativity ceiling were taken at the unitarity value λ ≈ 25, the relic-density ceiling would rise well above the direct-detection bound in every scheme, leaving a viable complex-scalar window. The real-scalar statement 'possibly up to 40-50 TeV' similarly uses λ = 15 and would shift upward. Because the conclusion is phrased as a model verdict rather than as a function of the adopted cutoff, the strongest claim is not robust to this convention.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript revisits the simplest scalar singlet dark matter models, real and complex, coupled to the Higgs boson through a portal interaction, and confronts them with the latest LZ direct-detection limits and the Planck relic-density measurement. The authors compute next-to-leading-order corrections to both the freeze-out annihilation cross section and the direct-detection cross section, working in three renormalization schemes (MSbar, OS-DD, and OS-FO) and providing analytic expressions in the m_h/m_phi expansion. Their main numerical results are that the real scalar is excluded below roughly 20-24 TeV at NLO (compared with about 30 TeV at LO), with a possible perturbative window up to approximately 40-50 TeV, while the complex scalar is claimed to be excluded in its entire perturbative regime, with only a narrow window near the Higgs resonance m_phi ~ m_h/2 remaining. The paper also discusses vacuum stability, Landau-pole scales, and the effect of a potentially large quartic self-coupling.","tokens_in":9864,"tokens_out":4378,"duration_ms":56774,"significance":"If the results hold, the paper is a useful and timely quantitative update for two benchmark WIMP models in the multi-TeV regime. Its main strengths are the explicit analytic NLO expressions, the transparent three-scheme comparison, and the robustness of the real-scalar lower bound, which varies only between 22 and 24 TeV across schemes; this is a solid, scheme-stable result. The paper also identifies a concrete experimental target: a factor-three improvement in spin-independent sensitivity would fully probe the real-scalar perturbative window. The complex-scalar verdict, however, is less robust than the abstract suggests, because it depends on the adopted perturbativity cutoff rather than on a scheme-independent calculation; this should be fixed before the paper is accepted.","major_comments":[{"comment":"The abstract claim that \"the entire perturbative region is excluded\" for the complex scalar is stronger than Table I supports. In the MSbar scheme the direct-detection bound is m_phi > 26 TeV while the perturbativity ceiling for lambda = 15 is 29 TeV, and in the OS-FO scheme the corresponding numbers are 32 TeV and 33 TeV. Thus, at the paper's own adopted cutoff lambda = 15, there are formally allowed windows of several TeV in two of the three schemes. The conclusion is only strictly true in the OS-DD scheme, where both numbers coincide at 41 TeV. The abstract and conclusions should be qualified, for example by saying that the complex scalar is excluded at the edge of perturbativity or that any allowed window is narrow and scheme-dependent.","section":"Abstract and Section IV, Table I"},{"comment":"The central exclusion claim for the complex scalar and the upper edge of the real-scalar allowed window depend directly on the perturbativity criterion lambda_phiH < 15 obtained by requiring the two-loop NLL contribution to the beta function to be at most 25% of the one-loop term. This is a defensible convention, but it is not a rigorous bound, and the same section quotes the tree-level unitarity bound |lambda_phiH| < 8*pi ~ 25. If the cutoff were taken at the unitarity value, the relic-density ceiling in Table I would rise well above the direct-detection bound in every scheme, leaving a viable complex-scalar window. The paper should present the allowed mass windows as a function of the adopted maximum coupling, or at least state explicitly that the model verdict is contingent on this 25% perturbativity convention rather than being a scheme-independent conclusion.","section":"Section IV, Eq. (9)"},{"comment":"There is a numerical inconsistency and a conceptual tension in the treatment of the complex scalar. Footnote 2 states that in the OS-FO and MSbar schemes the direct-detection constraint turns back at about 34 TeV and 26 TeV, but Table I lists the OS-FO direct-detection bound as 32 TeV, not 34 TeV. More importantly, if the fact that the NLO direct-detection bound crosses the perturbativity ceiling is taken as evidence that non-perturbativity is close in those schemes, then the OS-DD case, where the two numbers coincide at 41 TeV, should be interpreted with the same standard. The authors should state precisely what criterion they use to declare the complex scalar excluded, and apply it uniformly across schemes.","section":"Section IV, footnote 2 and Table I"}],"minor_comments":[{"comment":"The abstract and the Figure 1 caption should specify that the displayed exclusion for the complex scalar corresponds to the OS-DD scheme; as written, the claim appears scheme-independent, which is not the case in Table I.","section":"Abstract and Figure 1 caption"},{"comment":"The table would be easier to read if the masses were explicitly labeled in TeV and if the heading made clear that the \"Pert.\" columns are the maximum masses consistent with the relic density at the perturbativity cutoff lambda_phiH = 15.","section":"Table I"},{"comment":"The notation x_h is used before it is defined; x_h = m_h/m_phi should be introduced immediately before Eq. (6), and the phrase \"s-wave contribution\" should be clarified, since the Boltzmann calculation includes p-wave contributions as well.","section":"Section III, Eqs. (6) and (7)"},{"comment":"The number 34 TeV quoted for the OS-FO scheme disagrees with the value 32 TeV in Table I; please correct this inconsistency.","section":"Section IV, footnote 2"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the technical calculation appears to be careful and reproducible. The real-scalar lower bound is a solid and useful result, and the scheme-dependence analysis is a genuine contribution. However, the headline claim about the complex scalar being excluded in its entire perturbative region is stronger than the paper's own Table I supports and is sensitive to the adopted perturbativity convention. This is a load-bearing issue for the abstract and conclusions, but it is fixable by rephrasing the claim and presenting the results as a function of the perturbativity cutoff. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real headline is that NLO corrections lower the direct detection bound for the real scalar singlet from ~30 TeV to ~22-24 TeV, and that bound is stable across their three renormalization schemes. The second headline—that the complex scalar's entire perturbative region is excluded—is weaker than it sounds. It sits on a perturbativity cutoff λ<15, derived from a 25% next-to-leading-log beta-function rule. If you instead take the tree-level unitarity bound λ<8π, their own Table I leaves a viable window for the complex scalar in the MSbar and OS-FO schemes. So the abstract overstates the model verdict.\n\nWhat's genuinely new is the systematic NLO treatment of both freeze-out and direct detection in three schemes, applied to the 2024 LZ limits. The scheme comparison is the paper's best feature: the real-scalar lower bound only shifts by a couple of TeV, and the scheme-invariance of ΔFO−ΔDD explains why. The analytic approximations in Eqs. (6)-(7) and the appendix formulas are useful. The resonant mφ≈mh/2 region is handled carefully with DRAKE, and the Landau-pole estimates are a nice extra. I also appreciate the explicit statement that the λφ quartic correction cannot reduce the lower bound because it enters with the opposite sign.\n\nThe soft spots are real but not fatal. The perturbativity criterion is a convention, and the paper treats it as a sharp threshold. That affects both the complex-scalar exclusion and the claim that the real scalar is perturbative up to 40-50 TeV. The text does include footnote 1, but the abstract and conclusions don't carry the caveat. I would want the revised version to say \"excluded under the λ<15 criterion\" rather than categorically excluded. Secondarily, the NLO calculation is presented through formulas with no code or independent cross-check; for a calculation-heavy paper a supplementary notebook would help, but this is minor. The fN-dependence check is fine.\n\nWho gets value: anyone setting WIMP direct detection targets, and phenomenologists working on scalar singlet extensions. The real-scalar bound will be cited. It deserves a serious referee. I would recommend acceptance with revision, mainly to qualify the complex-scalar claim and make the perturbativity dependence explicit.","headline":"Solid NLO update on scalar singlet dark matter; the real-scalar lower bound is robust across schemes, but the 'entire perturbative region excluded' claim for the complex scalar depends on a perturbativity convention and overstates the model verdict.","tokens_in":10397,"tokens_out":2834,"would_cite":true,"duration_ms":33019,"reading_group":"yes","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 argues that NLO corrections push the simplest scalar dark matter models to the edge of perturbativity, leaving only a narrow viable window for the real scalar singlet.","keywords":["singlet scalar dark matter","Higgs portal","thermal freeze-out","direct detection","next-to-leading order","perturbativity","LZ experiment","WIMP"],"falsifier":"A spin-independent direct-detection search with three times LZ's current sensitivity that finds no recoil events for WIMP masses between 20 and 50 TeV would contradict the paper's conclusion that the real scalar singlet has a surviving perturbative window.","tokens_in":9358,"feed_emoji":"🌌","tokens_out":8317,"duration_ms":95007,"temperature":0.7,"pith_summary":"This paper asks how far the simplest possible thermal dark matter model, a single real or complex scalar field coupled only to the Higgs boson, can survive once next-to-leading-order radiative corrections are included. Using the latest LZ direct-detection limits, it finds the complex scalar version is excluded across the whole range where perturbative calculations can be trusted. For the real scalar version, NLO corrections lower the required dark matter mass to roughly 20 TeV, leaving a narrow viable window from about 20 TeV up to roughly 40-50 TeV with a Higgs-portal coupling of order 7-15. A three-fold improvement in direct-detection sensitivity would fully test this remaining window, and only a narrow resonant region near half the Higgs mass survives at low masses.","feed_headline":"Complex scalar WIMP ruled out; real survives above 20 TeV","feed_subtitle":"After next-to-leading-order corrections, only the real scalar survives, in a 20-50 TeV window.","key_machinery":"The load-bearing object is the Higgs-portal interaction $\\lambda_{\\phi H} H^\\dagger H \\phi^\\dagger\\phi$, the single coupling that controls both thermal freeze-out of the dark scalar and its spin-independent scattering off nucleons. The argument proceeds by computing the NLO corrections to both rates in three renormalization schemes (MSbar, OS-DD, OS-FO), chosen so that scheme dependence reveals the size of missing higher-order terms. Perturbativity is quantified by requiring the two-loop next-to-leading-log contribution to the $\\beta$ function to be at most 25% of the one-loop contribution, which gives $\\lambda_{\\phi H}\\lesssim 15$, and Landau-pole locations are estimated from the same running.","core_discovery":"The central claim is that, at next-to-leading order, the complex scalar singlet model is excluded by direct detection over its entire perturbative parameter space, while the real scalar singlet model survives only in a narrow high-mass window. For the real scalar the lower mass bound moves from about 31 TeV at leading order to roughly 22-24 TeV at NLO depending on renormalization scheme, with the upper edge of the perturbative regime around 40-50 TeV; below about 20 TeV the model is firmly excluded. For the complex scalar the NLO lower bound is 26-40 TeV depending on scheme, which lies at or above the mass reachable with a perturbative coupling, so no perturbative solution remains. The paper also identifies a scheme-independent combination: the difference between the NLO correction to freeze-out and the NLO correction to direct detection is independent of the renormalization scheme.","pith_inferences":["Editorial: The exclusion of the complex scalar stands or falls with the $\\lambda_{\\phi H}\\lesssim 15$ perturbativity threshold; if the boundary were instead the tree-level unitarity limit $|\\lambda_{\\phi H}|<8\\pi\\approx 25$, the complex scalar would keep a viable window around 26-50 TeV in the MSbar scheme.","Editorial: The paper's Landau-pole estimates imply that, if the real scalar is indeed the dark matter, new physics must enter below roughly 300-500 TeV (or below about 7400 TeV in the most conservative scheme) to tame the running coupling, making the scenario testable in principle through signatures of that new sector.","Editorial: The scheme-independence of $\\Delta^{\\rm FO}_{\\rm NLO}-\\Delta^{\\rm DD}_{\\rm NLO}$ suggests a robust relation between relic-density and direct-detection rates that could be carried over to other Higgs-portal dark matter models."],"forward_implications":["The real scalar singlet model can be a complete thermal dark matter candidate only with $m_\\phi$ between about 20 TeV and roughly 40-50 TeV and Higgs-portal coupling $\\lambda_{\\phi H}$ between about 7 and 15.","The complex scalar singlet version is ruled out in its perturbative regime, so in this model dark matter cannot be a complex singlet unless non-perturbative dynamics or additional fields intervene.","A three-fold improvement in spin-independent direct-detection sensitivity will fully probe, and if null exclude, the real scalar's remaining perturbative window.","Both models still allow a narrow resonant region near $m_\\phi\\simeq 59-63$ GeV, which would need roughly a factor 12-25 improvement in direct-detection sensitivity to be fully tested.","Because the annihilation that sets the relic density proceeds through the same $H^\\dagger H \\phi^\\dagger\\phi$ operator, these conclusions carry over to other dark matter models whose freeze-out is driven by the Higgs portal."],"supporting_citations":[{"why":"Provides the 90% C.L. LZ spin-independent direct-detection limit that sets the lower mass bounds in both scalar cases.","marker":"[1]"},{"why":"Underlies the relic-abundance calculation procedure used to translate the observed dark matter density into a coupling-versus-mass relation.","marker":"[15]"},{"why":"Supplies the Higgs-nucleon coupling value $f_N=0.308$ used to convert direct-detection limits into constraints on the portal coupling.","marker":"[20]"},{"why":"Statistical treatment used to turn the LZ upper limit into a 90% C.L. bound on the portal coupling.","marker":"[23]"},{"why":"Basis for the perturbativity criterion that the next-to-leading-log beta-function contribution be no more than 25% of the leading-log one.","marker":"[24]"},{"why":"Lattice results used to calibrate where the perturbative description breaks down in quartic-scalar theories.","marker":"[25]"},{"why":"Provides the application of this criterion to the SM quartic coupling and the Landau-pole prescription used here.","marker":"[26]"},{"why":"Used to compute the one- and two-loop beta functions of the Higgs-portal coupling.","marker":"[27, 28]"},{"why":"Underpins the resonant-region relic-abundance calculation including kinetic-equilibrium effects.","marker":"[37, 38]"}],"fun_headline_variants":["NLO excludes complex WIMP; real WIMP survives only 20-50 TeV","Complex scalar WIMP dead after NLO; real left at 20-50 TeV","NLO shrinks scalar WIMP window: real only, 20-50 TeV","Real scalar WIMP: NLO mass window 20-50 TeV, complex excluded","LZ limits plus NLO: complex scalar ruled out, real at 20-50 TeV"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the definition of where perturbativity ends: the paper takes the theory to be perturbative only while $\\lambda_{\\phi H}\\lesssim 15$, a threshold obtained by requiring the two-loop contribution to the $\\beta$ function to be at most 25% of the one-loop contribution; if the true boundary is higher, the complex scalar could still have an allowed window.","fun_headline_variants_meta":{"raw":{"variants":["NLO excludes complex WIMP; real WIMP survives only 20-50 TeV","Complex scalar WIMP dead after NLO; real left at 20-50 TeV","NLO shrinks scalar WIMP window: real only, 20-50 TeV","Real scalar WIMP: NLO mass window 20-50 TeV, complex excluded","LZ limits plus NLO: complex scalar ruled out, real at 20-50 TeV"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001033,"raw_usage":{"total_tokens":4348,"prompt_tokens":938,"completion_tokens":3410,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":554,"completion_tokens_details":{"reasoning_tokens":3292}},"tokens_in":554,"tokens_out":3410,"duration_ms":25001,"temperature":1.0,"reasoning_tokens":3292,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:52:25.706540+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A spin-independent direct-detection search with three times LZ's current sensitivity that finds no recoil events for WIMP masses between 20 and 50 TeV would contradict the paper's conclusion that the real scalar singlet has a surviving perturbative window.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the 90% C.L. LZ spin-independent direct-detection limit that sets the lower mass bounds in both scalar cases."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the application of this criterion to the SM quartic coupling and the Landau-pole prescription used here."}],"review_version":1}