{"id":"898c1707-6f00-4c69-b5dd-57b4ac716c2b","arxiv_id":"2505.13074","paper_version":4,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Future GW experiments retain sensitivity to several SMEFT operators despite renormalization scale uncertainties, provided the leading (H†H)^3 operator is measured by colliders.","lead":"Future gravitational wave observatories can still measure new physics in the Higgs potential even when renormalization scale uncertainties are included. The catch is that the dominant (H†H)^3 operator must first be pinned down by future collider experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Scale uncertainty is compared at two fixed scales rather than marginalised over; the advertised sensitivity 'in the presence of scale uncertainties' is therefore not actually demonstrated.","rationale":"The reader's weakest assumption concerns the validity of the dimension-six SMEFT truncation for the SFO-EWPT, which the paper itself acknowledges as a benchmark limitation. That is a real concern but it is external to the paper's explicitly conditional claim: the abstract says sensitivity survives 'provided that' CH is precisely measured, and the analysis adopts CH as the driver. The more load-bearing internal issue is whether the renormalization-scale uncertainty is actually incorporated into the statistical statement. The paper evaluates Fisher contours at two fixed scales and reports them side by side; it never constructs a likelihood that includes μ as an unknown. Consequently, the claim that GW observations 'remain sensitive' under scale uncertainty is not directly supported by the figures. This is not an accusation of error in the fixed-scale Fisher computations; it is a gap between the stated conclusion and the method. The proposed check is deliberately narrow: adding μ as a nuisance parameter for a single representative operator would settle whether the two-scale comparison is representative of the true sensitivity. If the check passes, the paper's qualitative conclusion is strengthened; if it fails, the central claim would need substantial qualification. I therefore keep the reader's CONDITIONAL verdict unchanged rather than escalating to rejection, because the paper does provide enough information to define the missing test and the underlying numerical framework is standard.","tokens_in":32956,"tokens_out":5257,"duration_ms":63032,"concrete_test":"Recompute the 95% C.L. interval for one representative operator (e.g., CH□ at Λ/√|cH| = 600 GeV, vb = 0.5, DECIGO, Tobs = 1 yr) with ln μ treated as a nuisance parameter: either profile over μ ∈ [Tn/2, 2πTn] with a flat prior, or add a theory-error term (∂ln Sh/∂ln μ)^2 σ_μ^2 to the covariance, with σ_μ = ln 4. If the marginalized interval remains within roughly a factor of two of the fixed-scale widths in Fig. 5, the claim survives; if it inflates by an order of magnitude or becomes unbounded, the advertised sensitivity is an artifact of conditioning on the unknown scale.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is conditional: once CH is fixed by colliders, GW sensitivity to other dimension-six operators survives renormalization-scale uncertainty. The analysis tests this by computing Fisher contours at μPT = 2πTn and Tn/2 (Figs. 5–8, Section 3.3) and showing that the interval widths are similar. But the Fisher likelihood in Eqs. (3.11)–(3.12) does not include μ as a parameter. Sh(f, μ, {p}) is evaluated at one fixed μ for each contour; the two contours are never combined. The scale uncertainty is thus not integrated out or profiled; it is only sampled at two endpoints. If the true μ is unknown, an observed spectrum can be fit by shifting μ as well as the Wilson coefficients. Because the peak amplitude changes by orders of magnitude between the two adopted scales (Fig. 4), the derivative ∂Sh/∂ln μ is large and may be strongly degenerate with, e.g., CH□ or CuH. Showing that fixed-μ intervals at two endpoints are similar does not establish that the marginalized interval is finite or comparable. The abstract's wording 'even in the presence of renormalization scale uncertainties' therefore overstates what Eqs. (3.11)–(3.12) and Figs. 5–8 demonstrate.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits the question of whether gravitational-wave (GW) observations can perform precision measurements of SMEFT Wilson coefficients when the strong theoretical scale uncertainty of the conventional daisy-resummed effective potential is taken into account. Using the SMEFT with the (H^†H)^3 operator as the source of a strongly first-order electroweak phase transition, the authors compute the one-loop improved Higgs potential, include one-loop RGE running of dimension-six operators between M_Z and the phase-transition scale, and evaluate phase-transition parameters and GW spectra at the two renormalization scales μ_PT = 2πT_n and μ_PT = T_n/2. They then perform Fisher-matrix forecasts for LISA, DECIGO, and BBO, assuming that future colliders fix the central value of the (H^†H)^3 Wilson coefficient. Their central claim is that, under this collider precondition, future GW observations remain sensitive to many other dimension-six operators even in the presence of renormalization-scale uncertainties.","tokens_in":33173,"tokens_out":3779,"duration_ms":45936,"significance":"If the main claim is established, the paper resolves a genuine tension in the literature: Ref. [47] showed that daisy-resummed GW predictions have large scale dependence, while earlier Fisher studies [48,82] ignored that uncertainty. The paper makes a useful step toward a more honest forecast by combining one-loop SMEFT RGE running, one-loop matching, and a two-scale comparison in the statistical analysis. It is also commendably transparent about its assumptions, including the conditional reliance on a collider measurement of C_H, the variation of the bubble wall velocity, and the known limitations of the dimension-six SMEFT benchmark. However, as discussed below, the statistical treatment of the scale uncertainty is not yet a marginalization or profiling over μ, and the pairwise Fisher analysis leaves open multi-operator degeneracies; these need to be addressed before the abstract's claim is fully supported.","major_comments":[{"comment":"The Fisher likelihood in Eqs. (3.11)–(3.12) does not include the renormalization scale μ as a parameter. The analysis evaluates Sh(f, μ, {p}) at two fixed values μ_PT = 2πT_n and μ_PT = T_n/2 and compares the resulting contours, but never marginalizes or profiles over μ. Since Fig. 4 shows that the peak amplitude changes by orders of magnitude between these two scales, an unknown μ is a potentially large systematic direction in parameter space that can be partially absorbed by Wilson coefficients. The abstract's wording, 'even in the presence of renormalization scale uncertainties', is therefore stronger than what Figs. 5–8 demonstrate. I recommend adding a profile-likelihood or nuisance-marginalized Fisher analysis in which μ is varied (e.g., over the interval [T_n/2, 2πT_n]), and reporting whether the marginalized widths remain comparable to the fixed-μ widths.","section":"§3.3, Eqs. (3.11)–(3.12)"},{"comment":"The Fisher forecasts are performed pairwise: only the GW source C_H and one other Wilson coefficient are taken as active, with all other operators set to zero. The real parameter space has many active operators simultaneously, and degeneracies among them (for example, between C_H□, C_HD, and operator effects entering through the top sector) could broaden the inferred intervals. Since the paper claims sensitivity to 'various' dimension-six operators, at least one multi-operator Fisher computation, or an explicit argument that the pairwise treatment is conservative, is needed to support the conclusion.","section":"§3.3, Figs. 5–8"},{"comment":"The scale uncertainty is sampled at only two endpoints, μ_PT = 2πT_n and T_n/2, rather than explored continuously. The phase-transition parameters and GW spectrum need not vary monotonically between these two points, and the two endpoints are also not accompanied by a clear criterion for why they bound the systematic error. A scan over μ_PT, or a statement of why the endpoints are representative, would make the robustness claim much stronger.","section":"§2.5 and Fig. 1"}],"minor_comments":[{"comment":"The first panel of each figure labels the operator 'CW', while Table 1 and the text consistently use 'OW' for the triple-gauge-boson operator. This should be corrected for consistency.","section":"Figs. 6–8"},{"comment":"The notation H*R* and H(T_n)R* is used in the same equation without defining H* at that point; it should be stated that H* is the Hubble parameter at the nucleation (or percolation) temperature, and the two notations should be made consistent.","section":"§3.1, Eq. (3.6)"},{"comment":"Fig. 2 plots v_c/T_c as a continuous function of μ_PT, while the rest of the paper evaluates only μ_PT = 2πT_n and T_n/2. A sentence clarifying that Fig. 2 is illustrative and not used in the Fisher analysis would be helpful.","section":"§2.5 and Fig. 2"},{"comment":"The use of Λ/√|c_H| alongside C_H = c_H/Λ² is sometimes confusing because Λ/√|c_H| has mass dimension one while C_H has mass dimension −2. Defining both notations in one place near Eq. (2.11) would improve readability.","section":"§3.3"},{"comment":"No code or numerical inputs are released. Given that the paper's quantitative conclusions depend on the bounce solver, RGE running, and detector noise curves, making the analysis scripts available, or at least specifying all numerical input values and software versions, would aid reproducibility.","section":"General"},{"comment":"The paper correctly emphasizes that the dimension-six SMEFT benchmark is limited by the results of Refs. [24,65]. However, the abstract and conclusions should state more explicitly that the quoted sensitivities apply to this benchmark scenario, not to generic new-physics models, to avoid over-generalization by readers.","section":"§1 and §4"}],"recommendation":"major_revision","confidential_remarks":"The paper is a competent and transparent extension of previous Fisher-matrix GW studies, and the central conditional claim is plausible. The main technical shortcoming is that the scale uncertainty is sampled, not statistically treated; this is fixable within the scope of the manuscript. The pairwise-only Fisher analysis is a second issue that should be addressed at least with a representative multi-operator example. No concerns about novelty or scope: the paper fits hep-ph and would be of interest to the phase-transition and GW community once the robustness analysis is completed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the paper does exactly what it says: it adds one-loop SMEFT RGE running and renormalization-scale variation to the earlier Fisher-matrix analysis of Ref. [48], and finds that if the dominant (H†H)^3 operator (CH) is pinned down by colliders, the projected GW sensitivities to the other dimension-six operators remain roughly the same at μPT = 2πTn and Tn/2. Second, that conclusion is conditional and honestly labeled: the paper does not claim to measure CH from GWs, and it flags the known EFT-truncation problems for the SFO-EWPT by citing Refs. [24,65].\n\nWhat is genuinely new: the one-loop RGE running of the SMEFT operators into the thermal potential, and the explicit check that the Fisher ellipses for each subleading operator are similar at the two scales. The computations look standard and mostly transparent; the long appendix tables of self-energies and RGE corrections are useful. The paper is careful with the central-value dependence of CH and varies vb too.\n\nSoft spots, in order of seriousness. First, the renormalization scale is sampled at two endpoints, not marginalized or profiled in the Fisher likelihood. The two contours in Fig. 5 are separate forecasts conditional on μ; they are never combined into a single error budget. The abstract's 'in the presence of renormalization scale uncertainties' overstates things. The stress-test note has a point: if μ were treated as unknown in a joint fit, degeneracy with the overall amplitude could inflate some of the quoted intervals. That said, this does not sink the paper; the similar widths at the two scales are meaningful evidence against a drastic degradation, and because CH is fixed externally, the overall amplitude anchors μ. But the authors should add a caveat, or better, do a simple profile over μ. Second, the Fisher analysis is pairwise, one operator at a time; simultaneous fits could change the picture. The authors note this implicitly but don't address it. Third, no code or data release; reproducibility suffers a bit. Fourth, the EFT truncation issue is openly acknowledged and is a shared problem of the whole SMEFT-GW literature; it shouldn't count against this paper in itself.\n\nThe paper deserves a serious referee. It is a solid, incremental contribution to a specific open question raised in Ref. [47], and the negative result about scale-robustness—if it holds up—is useful for planning forecasts. Recommend sending to peer review, with the scale-marginalization issue as the main request for revision.","headline":"A solid, honest extension of the authors' earlier Fisher analysis: once the (H†H)^3 coefficient is fixed externally, GW sensitivities to subleading SMEFT operators survive the two scale choices tested, though the scale uncertainty is sampled, not yet marginalized.","tokens_in":33717,"tokens_out":4302,"would_cite":true,"duration_ms":45767,"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":"Future GW observatories can still measure new-physics Higgs operators despite renormalization-scale uncertainty—if colliders first pin down the $(H^\\dagger H)^3$ operator.","keywords":["gravitational waves","electroweak phase transition","SMEFT","Fisher matrix","renormalization scale uncertainty","daisy resummation","Higgs potential","dimension-six operators"],"falsifier":"Measure the coefficient of $(H^\\dagger H)^3$ at a future collider with the projected precision; if its uncertainty is comparable to or larger than the spread between the $\\mu_{\\rm PT}=2\\pi T_n$ and $T_n/2$ Fisher contours, the claimed robustness fails. A cleaner test is to repeat the same Fisher forecast in a dimensionally reduced or lattice effective theory: if the renormalization-scale dependence of the peak amplitude is not reduced there, then externally fixing $(H^\\dagger H)^3$ is not enough to make the GW precision measurements reliable.","tokens_in":32718,"feed_emoji":"🌊","tokens_out":5111,"duration_ms":54746,"temperature":0.7,"pith_summary":"This paper asks whether gravitational-wave (GW) observations can still perform precision new-physics measurements once the recognized renormalization-scale ambiguity of the daisy-resummed effective potential is taken into account. Using SMEFT as a benchmark, the paper includes one-loop RGE running of dimension-six operators in the Higgs potential, computes GW spectra at the two scales $\\mu_{\\rm PT}=2\\pi T_n$ and $T_n/2$, and runs a Fisher matrix forecast. It finds that while the absolute GW peak amplitude shifts with the renormalization scale, the relative sensitivities to subleading Wilson coefficients survive, provided the $(H^\\dagger H)^3$ operator is fixed by other experiments. This restores the GW precision program, but only conditionally: the transition-driving operator must be known externally.","feed_headline":"GW detectors stay sensitive to new physics despite scale uncertainty","feed_subtitle":"Once colliders measure the Higgs six-point operator, gravitational-wave forecasts keep their precision.","key_machinery":"The carrying mechanism is the full one-loop effective potential $V_{\\rm full}=V+V_{\\rm CW}+V_T+V_{\\rm daisy}$, evaluated at two renormalization scales that differ by a factor of $4\\pi$ ($\\mu_{\\rm PT}=T_n/2$ and $2\\pi T_n$), with one-loop SMEFT RGE running of $m^2$, $\\lambda$, the top Yukawa, the gauge couplings, and the Wilson coefficients from $\\mu=M_Z$ to the transition scale. This scale variation quantifies the theoretical uncertainty. On top of it, the Fisher information matrix $F_{ab}=2T_{\\rm obs}\\int df\\,\\partial_{p_a}S_h\\,\\partial_{p_b}S_h/[S_{\\rm eff}+S_h]^2$ converts derivatives of the GW power spectrum with respect to Wilson coefficients into 95% confidence contours, giving the projected measurement precision.","core_discovery":"The central claim is that future GW observations, such as DECIGO and BBO, can remain sensitive to a wide set of dimension-six SMEFT operators even when the renormalization-scale uncertainty of the daisy-resummed approach is explicitly included. The paper shows that the scale dependence makes it impossible to determine the Wilson coefficient $C_H$ of the $(H^\\dagger H)^3$ operator from GW data alone; the peak amplitude varies by one to two orders of magnitude between $\\mu_{\\rm PT}=2\\pi T_n$ and $T_n/2$. However, once $C_H$ is fixed by future collider measurements, the 95% confidence contours for the other operators stay narrow, with sensitivity to new physics scales above roughly 10 TeV, and this robustness holds for bubble wall velocities $v_b=0.2$, $0.5$, and $1$. The paper therefore concludes that precision new-physics searches via GW observations remain viable under renormalization-scale uncertainties, conditional on external determination of the phase-transition-driving operator.","pith_inferences":["The paper does not say this, but its results imply a practical division of labor: colliders fix the strength of the transition, while GW detectors constrain the shape of the potential. If the collider uncertainty on $C_H$ is comparable to the spread between the two renormalization-scale contours, that division breaks down and the GW constraints degrade accordingly.","The cited limitations of the dimension-six SMEFT benchmark suggest that the sensitivities may not transfer to realistic ultraviolet completions; a direct extension would repeat the Fisher forecast for a singlet-extended model or with dimension-eight operators to test whether the robustness survives.","The two-scale comparison is a conservative envelope rather than a statistical error; one could combine the two contours into a single band, but that would still not capture systematic errors from the choice of phase-transition parameters such as the percolation temperature.","Because the forecasts rely on the acoustic GW peak, foreground subtraction of compact white dwarf binaries in the millihertz band is likely to set the practical floor for the achievable precision, which the paper includes but does not vary as a nuisance parameter."],"forward_implications":["If the central claim is correct, GW observations alone cannot determine the coefficient of the $(H^\\dagger H)^3$ operator; a collider measurement of that coefficient becomes a prerequisite for GW-based new-physics searches.","The precision for subleading dimension-six operators is largely preserved across the two renormalization scales and across bubble wall velocities from 0.2 to 1.","The projected sensitivities reach new physics scales above roughly 10 TeV for DECIGO and BBO with one year of observation.","The analysis suggests that a dimensionally reduced effective theory, applied to the same Fisher forecast, should shrink the renormalization-scale uncertainty substantially.","Varying the assumed central value of $C_H$ between the collider-reachable points (600 to 700 GeV for $\\Lambda/\\sqrt{|c_H|}$) does not erase the GW sensitivity to the other operators."],"supporting_citations":[{"why":"Establishes the daisy-resummed renormalization-scale uncertainty in the GW peak amplitude that this paper explicitly incorporates and works around.","marker":"[47]"},{"why":"Previous Fisher matrix analysis of SMEFT effects on the GW spectrum without renormalization-scale uncertainties; the method extended here.","marker":"[48]"},{"why":"Provides future collider projections for measuring $C_H$, which the paper assumes as the external input that fixes the GW source.","marker":"[71]"},{"why":"Shows the limited applicability of the SMEFT dimension-six description for the singlet-extended Standard Model; a load-bearing caveat on the benchmark.","marker":"[24]"},{"why":"Argues that the dimension-six SMEFT truncation fails to capture the phase transition in a wide range of theories; the main limitation on the benchmark scenario.","marker":"[65]"},{"why":"Supplies the Fisher information matrix formalism for stochastic GW background correlation analysis used in the statistical forecast.","marker":"[63]"},{"why":"Earlier Fisher matrix analysis selecting models through the synergy between collider and GW experiments, providing the methodological baseline.","marker":"[82]"},{"why":"Defines the dimension-six SMEFT operator basis and the normalization of the Wilson coefficients used throughout the paper.","marker":"[12]"}],"fun_headline_variants":["GW forecasts survive scale uncertainty once colliders pin down Higgs operator","Scale uncertainty can't kill GW new-physics reach if colliders fix Higgs term","GW precision survives scale drift when colliders measure Higgs six-point operator","Colliders hold the key: GW new-physics reach robust to scale uncertainty"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole forecast rests on the assumption that the dimension-six SMEFT, with the $(H^\\dagger H)^3$ operator driving the transition, is a valid stand-in for real new physics; the paper itself cites studies showing this fails for many complete theories of new physics at high energy.","fun_headline_variants_meta":{"raw":{"variants":["GW forecasts survive scale uncertainty once colliders pin down Higgs operator","Scale uncertainty can't kill GW new-physics reach if colliders fix Higgs term","GW precision survives scale drift when colliders measure Higgs six-point operator","Colliders hold the key: GW new-physics reach robust to scale uncertainty"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000707,"raw_usage":{"total_tokens":3203,"prompt_tokens":977,"completion_tokens":2226,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":593,"completion_tokens_details":{"reasoning_tokens":2147}},"tokens_in":593,"tokens_out":2226,"duration_ms":17482,"temperature":1.0,"reasoning_tokens":2147,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:20:47.509305+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the coefficient of $(H^\\dagger H)^3$ at a future collider with the projected precision; if its uncertainty is comparable to or larger than the spread between the $\\mu_{\\rm PT}=2\\pi T_n$ and $T_n/2$ Fisher contours, the claimed robustness fails. A cleaner test is to repeat the same Fisher forecast in a dimensionally reduced or lattice effective theory: if the renormalization-scale dependence of the peak amplitude is not reduced there, then externally fixing $(H^\\dagger H)^3$ is not enough to make the GW precision measurements reliable.","supporting_citations":[{"cited_title":"Correlation analysis of stochastic gravitational wave background around 0.1-1Hz","cited_arxiv_id":"gr-qc/0510067","evidence_quote":"Supplies the Fisher information matrix formalism for stochastic GW background correlation analysis used in the statistical forecast."}],"review_version":1}