{"id":"01dbb401-502c-40e1-84f4-d7e79b639cf1","arxiv_id":"2502.00460","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The paper reports preliminary gradient-flow lattice results for the Weinberg and mixed susceptibilities, but statistical errors are too large to extract a useful Peccei-Quinn theta contribution.","lead":"Lattice QCD calculations of the topological, Weinberg, and mixed CP-violating susceptibilities are presented using the gradient flow scheme with continuum and chiral extrapolations. The work is a step toward determining the theta contribution induced by the Weinberg operator, which is relevant for electric dipole moment searches beyond the Standard Model.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The load-bearing step is the chiral-continuum extrapolation of chi_W and chi_M; it rests on AIC-selected fits to 7 of 11 ensembles with unquantified M_K/finite-volume terms, and the displayed chi_Q fit is not the AIC-best form.","rationale":"The reader's weakest_assumption identifies the same Achilles heel. I considered the strong flow-time dependence of chi_W reported in Sec. 8, but that is a symptom; the direct threat to the central claim is that the extrapolated chi_W and chi_M can change with fit form and ensemble set. The paper's own Table 2 shows the displayed chi_Q fit is not the AIC-best, and Section 7.2 admits M_K/FV ignorability. These choices are not supported by the data alone. A single inclusive refit with M_K^2 and finite-volume terms would settle the matter. Independent support: chi_Q ~ (71 MeV)^4 matches prior determinations, and the gradient-flow scheme is theoretically sound; the concern is about the extrapolation, not the formalism. Hence the verdict should remain CONDITIONAL: the extrapolated Weinberg and mixed numbers should not be used until the inclusive refit is done.","tokens_in":9447,"tokens_out":17601,"duration_ms":190024,"concrete_test":"Refit the full 11-ensemble data set with Eqs. (7)-(8) extended by an M_K^2 term and a finite-volume term, without excluding C13, D220, D5L, and D6; then check whether the continuum physical-mass values of chi_W and chi_M, and the ratio chi_M/chi_Q, shift by more than the quoted 1-sigma uncertainties. If they do, the displayed AIC-selected 7-ensemble extrapolations are not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the continuum and physical-mass extrapolations of the mixed and Weinberg susceptibilities are controlled. Section 7.2 models the a and M_pi^2 dependence by the low-order polynomials in Eqs. (7)-(8), with coefficients as cubic splines in flow time, and states that M_K and finite-volume effects are 'assumed negligible'. Only 7 of the 11 ensembles in Table 1 are used; C13, D220, D5L, and D6 are excluded, although D6 (a=0.091 fm, M_pi=175 MeV, L~4.4 fm) does not obviously fit the stated criteria of coarse spacing, heavy pion, or small volume. The displayed chi_Q result is the Eq. (7) a^2 form, which has AIC 142.9 and chi^2/dof=2.19, while the AIC-best chi_Q fit in Table 2 is the Eq. (8) a form with AIC 106.2; the choice is justified by the theoretical expectation that chi_Q vanishes in the chiral-continuum limit, not by data. Since chi_W and chi_M are new and their extrapolated values are the main products, the fit-form choice, the ensemble exclusions, and the neglected M_K/FV terms could each shift the central extrapolations by more than the quoted statistical errors. This is the load-bearing weak point of the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a preliminary lattice-QCD study of three CP-violating gluonic susceptibilities: the topological susceptibility chi_Q, the Weinberg-operator susceptibility chi_W, and the mixed susceptibility chi_M. The calculation uses the gradient-flow scheme on 11 ensembles of 2+1-flavor clover fermions, models the flow-time dependence with cubic splines, and extrapolates to the continuum and physical-pion limits using low-order polynomials in the lattice spacing and M_pi^2 (Eqs. 7 and 8). The main reported results are that chi_Q is approximately (71 MeV)^4 at the physical point and consistent with prior determinations, that chi_W has a strong flow-time dependence, and that the ratio chi_M/chi_Q entering the Peccei-Quinn induced theta (Eq. 9) has statistical uncertainties too large for a useful estimate. The paper explicitly labels the results as preliminary and states that M_K and finite-volume effects are neglected.","tokens_in":9771,"tokens_out":6245,"duration_ms":63756,"significance":"If the continuum and chiral extrapolations are controlled, this is a first nonperturbative calculation of the Weinberg and mixed susceptibilities in the gradient-flow scheme, and it would provide an important ingredient for future EDM calculations. The paper has several strengths: the susceptibilities are directly measured lattice quantities rather than fitted parameters; the topological susceptibility is checked against prior determinations; the gradient-flow renormalization framework is well motivated; and the authors are explicit that the theta-induced estimate is not yet useful. The significance is limited by the absence of final numbers for chi_W and chi_M and by the unquantified systematic uncertainties in the extrapolations.","major_comments":[{"comment":"The fits displayed in Figs. 1 and 2 use the three-parameter a^2 ansatz of Eq. (7). For chi_Q, the AIC strongly favors the four-parameter a ansatz of Eq. (8): AIC 106.234 with chi^2/dof 1.382, versus AIC 142.890 with chi^2/dof 2.189 for the displayed form. The choice is justified by the theoretical expectation that chi_Q vanishes in the chiral-continuum limit, not by the data. Since chi_W and chi_M are new results and the fits are the main product of the paper, it is load-bearing to show how the extrapolated values, and the ratio in Eq. (9), change under all four fit forms or under an AIC-based model average. Without that, the central claim that the continuum and chiral extrapolations are under control is not supported.","section":"Section 7.3, Table 2"},{"comment":"The text states that the fits exclude C13, D220, D5L, and D6 to avoid coarse lattice spacings, heavy pion masses, or small volumes. This criterion does not obviously apply to D6, which has a=0.0914 fm, M_pi=175 MeV, and L=48a ~ 4.4 fm: it is neither the coarsest nor the heaviest ensemble, and its volume is only modestly smaller than that of D7, which is included. If D6 is excluded because of finite-volume effects, the volume criterion should be quantified (e.g., M_pi L values) and the stability of the fits with and without D6 should be demonstrated. The exclusion of four of eleven ensembles, particularly a near-physical-pion ensemble, can shift the extrapolated central values by more than the quoted statistical errors.","section":"Section 7.2, Table 1"},{"comment":"The fit model (Eqs. 7 and 8) contains no M_K dependence, and Section 7.2 states that contributions from M_K being heavier than physical are 'assumed negligible'. However, Table 1 shows M_K varies from 476 MeV (C13) to 575 MeV (E5) across the ensembles, with several values more than 10% above the physical kaon mass. Since the included ensembles have unphysical strange-quark masses, the omission of an M_K^2 term (or an equivalent constraint to the physical kaon mass) can bias the intercept of the M_pi^2 extrapolation. The paper should provide an estimate of this systematic effect, for example by adding a term proportional to (M_K^2 - M_K^2_phys) to the fit or by showing that the seven included ensembles have sufficiently physical M_K values.","section":"Section 7.2, Table 1"},{"comment":"The paper reports that chi_W has a strong flow-time dependence, yet the continuum extrapolation is performed at fixed flow time with the spline coefficients s_i(t_gf). For a scheme defined by the gradient flow, the physical reference scale is set by t_gf, and the continuum limit should be taken at fixed physical flow time with t_gf/a^2 -> infinity. The paper does not specify the flow-time window over which the continuum extrapolation is stable, nor does it show the extrapolated values at a chosen reference flow time with their uncertainties. The right panel of Fig. 2 shows the ratio chi_M/chi_Q without any error band, which makes it difficult to verify the statement that the statistical uncertainties are too large for a useful estimate. This point is load-bearing for the central claim that the method produces controlled extrapolations of these susceptibilities.","section":"Section 8, Figs. 1 and 2"}],"minor_comments":[{"comment":"The definitions of chi_W and chi_M are written as products of integrals without expectation values; they should be defined as variances and covariances, for example chi_W = <Q_W^2> and chi_M = <Q Q_W>, where Q_W is the volume integral of the Weinberg operator.","section":"Section 6"},{"comment":"There is a typo in 'for both both valence and sea quarks'; the doubled word should be removed.","section":"Section 7.1"},{"comment":"The word 'vaccuum' should be 'vacuum'.","section":"Section 6"},{"comment":"The phrase 'an useful estimate' should be 'a useful estimate'.","section":"Section 8"},{"comment":"The entries labeled 1/(tau_gf a)^2 appear to have inconsistent dimensions for a mixing coefficient between a dimension-four and a dimension-six operator; the notation relating tau_gf and t_gf should be checked and clarified.","section":"Eq. (5)"},{"comment":"For chi_W, the AIC values of the four fit forms differ by less than 0.3; the text should state whether the differences are considered significant or treated as statistically equivalent.","section":"Table 2"}],"recommendation":"major_revision","confidential_remarks":"This is a LATTICE2024 proceedings contribution and the authors are appropriately cautious, repeatedly labeling the results as preliminary and the theta-induced estimate as not yet useful. The main issue is that the displayed continuum and chiral extrapolations rest on fit-form choices, ensemble exclusions, and neglected M_K and finite-volume terms that could shift the new chi_W and chi_M results by more than the quoted statistical errors. The paper would be acceptable after a revision that adds robustness checks, quantitative estimates of the neglected systematics, and a clearer specification of the flow-time window for the continuum limit. I see no circularity issue: the susceptibilities are directly measured quantities, and the topological-susceptibility check against prior work is a genuine validation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First: this is a genuinely new lattice calculation — the Weinberg and mixed CP-violating susceptibilities in the gradient-flow scheme, with continuum/chiral extrapolations attempted. It is also an honest proceedings paper: it says plainly that chi_W has strong flow-time dependence and that chi_M/chi_Q is too noisy for a useful PQ-theta estimate. The topological sector agrees with earlier determinations, which is a good sanity check. The citation pattern is normal; close citations to the group's own theta-term work are appropriate because this is a direct extension.\n\nThe soft spots are in the fit/systematics. The extrapolated curves rest on low-order polynomial ansaetze in a and M_pi^2, with coefficients as cubic splines in flow time, fit to 7 of 11 ensembles. Four ensembles are dropped \"to avoid coarse lattice spacing, heavy pion masses or small volumes\"; D6 (a ~ 0.09 fm, M_pi ~ 175 MeV, L ~ 4.4 fm) does not obviously fit those categories, and the paper does not quantify what changes if it is included. M_K and finite-volume terms are 'assumed negligible' without an estimate. For chi_Q, the displayed fit is not the AIC-best one (Table 2: Eq. 8 with a has AIC 106 vs. 143 for the shown Eq. 7 a^2 form); the choice is defended by the theoretical expectation that chi_Q vanishes in the chiral-continuum limit, not by data. That is a defensible physics choice, but it is exactly where a skeptic wants sensitivity checks.\n\nNone of this is fatal to the central claim, because the central claim is 'this route is feasible,' not 'here is the final number.' The paper itself says the interesting ratio is not yet determinable. The problem is that the extrapolated chi_W and chi_M values are the main new products, and each of the above choices could shift them by more than the quoted statistical error. So the right reading is: first step, not final result.\n\nI would send this to a referee. The calculation is new, the literature is engaged, and the limitations are acknowledged in-text. A competent referee should ask for an ensemble-inclusion study and a bound on M_K/FV terms, not reject the paper. For a LATTICE proceedings it clears the bar; for a full journal paper it would need those systematics worked out.","headline":"First lattice results for the Weinberg and mixed susceptibilities: a genuine but explicitly preliminary calculation whose extrapolated numbers should not be used until the fit-selection and systematics are nailed down.","tokens_in":10395,"tokens_out":4044,"would_cite":true,"duration_ms":39691,"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":"A gradient-flow lattice calculation of the Weinberg operator's CP-violating susceptibilities reports chi_Q ~ (71 MeV)^4 at physical quark masses, but the ratio chi_M/chi_Q that sets the Peccei-Quinn-induced theta remains too noisy to quote.","keywords":["Weinberg operator","gradient flow","lattice QCD","CP violation","topological susceptibility","Peccei-Quinn mechanism","electric dipole moment","chiral extrapolation"],"falsifier":"Produce a new ensemble with lattice spacing smaller than the smallest used here ($a<0.05$ fm), pion mass within a few MeV of the physical value, a larger spatial volume, and a near-physical kaon mass, then compare $\\chi_Q$ and $\\chi_W$ measured at flow times $\\tau_{\\mathrm{gf}} \\simeq 0.4$--$0.6$ fm with the error band of the paper's extrapolation; a mismatch beyond the combined uncertainties would falsify the extrapolation.","tokens_in":9200,"feed_emoji":"⚛️","tokens_out":14220,"duration_ms":111947,"temperature":0.7,"pith_summary":"This paper tries to establish that the CP-violating Weinberg three-gluon operator can be studied nonperturbatively in lattice QCD using the gradient-flow scheme, where flowed composite operators need no extra renormalization. From variances and covariances of the flowed topological and Weinberg charges, the authors compute the topological, Weinberg, and mixed susceptibilities on eleven ensembles, model their flow-time dependence with cubic splines, and extrapolate to the continuum at physical quark masses. They find a topological susceptibility of about $(71\\,\\mathrm{MeV})^4$ at the physical point, consistent with earlier determinations, and a strong flow-time dependence for the Weinberg susceptibility. The ratio of the mixed to the topological susceptibility, which would fix the $\\theta$ angle induced by the Weinberg operator in Peccei-Quinn theories, remains too imprecise for a useful estimate. These are the nonperturbative inputs needed, alongside nucleon matrix elements, to convert EDM measurements into constraints on beyond-standard-model CP violation.","feed_headline":"Weinberg-operator response computed; key ratio still noisy","feed_subtitle":"The Peccei-Quinn-induced theta remains out of reach, but the lattice method for the Weinberg operator is now in place.","key_machinery":"The central machinery is the gradient-flow scheme, which smears gauge fields to a scale $\\tau_{\\mathrm{gf}} = \\sqrt{8t_{\\mathrm{gf}}}\\,a$ and lets the lattice cutoff be removed while composite operators remain finite, so the flowed charges need no extra operator renormalization. The susceptibilities are then formed as variances and covariances of the volume integrals $Q = \\int d^4x\\, G\\cdot\\tilde G$ and $W = \\int d^4x\\, G\\cdot\\tilde G\\cdot G$, which defines renormalized quantities directly. The extrapolation is carried by a fitting ansatz whose coefficients are low-order polynomials in the lattice spacing and $M_\\pi^2$, with each coefficient modeled as a cubic spline in flow time and the spline knots chosen by an information criterion.","core_discovery":"On the paper's own terms, the central discovery is that the gradient-flow scheme gives a workable nonperturbative definition of the renormalized Weinberg operator: the three susceptibilities $\\chi_Q$, $\\chi_W$, and $\\chi_M$ can be obtained as second derivatives of the effective action, equivalently as variances and covariances of the integrated topological charge and the Weinberg charge, and their combined dependence on flow time, lattice spacing, and pion mass can be extrapolated to the physical continuum. The extrapolated $\\chi_Q$ is consistent with zero in the chiral-continuum limit and is approximately $(71\\,\\mathrm{MeV})^4$ at physical quark masses. $\\chi_W$ depends strongly on flow time, and the paper reports that the statistical uncertainties in the ratio $\\chi_M/\\chi_Q$, which would set $\\theta_{\\mathrm{induced}} = -w\\,\\chi_M/\\chi_Q$ in Peccei-Quinn theories, are still too large to yield a useful number.","pith_inferences":["Because the gradient-flow scheme avoids hard-cutoff power-divergence subtractions, the same flowed operators could be reused for nucleon matrix elements of the Weinberg operator, giving a single consistent scheme for the full EDM calculation; the paper motivates but does not carry out that step.","The strong flow-time dependence of $\\chi_W$ may be dominated by mixing with the topological charge and by contact terms; checking whether the ratio $\\chi_M/\\chi_Q$ is more flow-time stable than either susceptibility alone could reveal how much of that dependence cancels.","A direct test of the fitting strategy would be to compute the ratio on a single ensemble with a physical pion mass and a large volume before any extrapolation; if it does not stabilize with flow time there, the continuum limit of the induced theta may be slow to approach.","If the $(71\\,\\mathrm{MeV})^4$ value and the large uncertainty in the ratio persist with higher statistics, near-term EDM searches would remain more sensitive to the topological term than to the Weinberg-induced theta, steering lattice effort toward the nucleon-side matrix elements."],"forward_implications":["If the extrapolated susceptibilities are correct, the topological susceptibility at the physical point is fixed near $(71\\,\\mathrm{MeV})^4$, giving a nonperturbative anchor for the QCD vacuum's response to a CP-violating theta term.","The chiral-continuum vanishing of $\\chi_Q$ is reproduced, consistent with the expectation that the theta term's effects disappear as quark masses go to zero.","Once the precision of $\\chi_M/\\chi_Q$ improves, the Peccei-Quinn-induced theta from the Weinberg operator can be computed directly from lattice QCD rather than estimated by models.","Together with nucleon matrix elements of quark-bilinear operators, these susceptibilities enable first-principles estimates of the neutron and proton electric dipole moments and CP-violating pion-nucleon couplings induced by the Weinberg operator, sharpening beyond-standard-model constraints."],"supporting_citations":[{"why":"Defines the gradient-flow renormalization scheme and establishes that flowed composite operators have finite matrix elements, the basis for computing the susceptibilities.","marker":"[24]"},{"why":"Supplies the one-loop matching of the CP-odd three-gluon operator to the gradient-flow scheme, fixing the renormalization constants Z_W and Z_M used in the conversion matrix.","marker":"[26]"},{"why":"Earlier theta-term EDM study whose chiral-continuum expectation for chi_Q motivates the choice among fit forms.","marker":"[31]"},{"why":"Prior lattice determination of the topological susceptibility that the paper's (71 MeV)^4 result is compared with.","marker":"[35]"},{"why":"Derives the relation theta_induced = -w chi_M/chi_Q for the Peccei-Quinn minimum, giving the physical target of the measured ratio.","marker":"[17]"},{"why":"Introduces the Peccei-Quinn mechanism that dynamically relaxes theta, the scenario in which the Weinberg-induced theta matters.","marker":"[18]"},{"why":"Provides the ensemble parameters, lattice spacings, pion and kaon masses, and volumes, that enter the fits.","marker":"[30]"}],"fun_headline_variants":["Gradient flow tames Weinberg operator, theta still elusive","Lattice maps Weinberg susceptibilities, ratio too noisy for theta","Weinberg response computed via flow, theta estimate out of reach","New lattice method for Weinberg operator, but ratio stays noisy","Flow scheme yields Weinberg susceptibilities, theta ratio still noisy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole extrapolation rests on the assumption that the low-order polynomial ansatz in lattice spacing and pion mass, with flow-time dependence represented by cubic splines and fitted to only seven of the eleven ensembles while taking kaon-mass and finite-volume effects to be negligible, correctly captures the approach to the physical continuum limit.","fun_headline_variants_meta":{"raw":{"variants":["Gradient flow tames Weinberg operator, theta still elusive","Lattice maps Weinberg susceptibilities, ratio too noisy for theta","Weinberg response computed via flow, theta estimate out of reach","New lattice method for Weinberg operator, but ratio stays noisy","Flow scheme yields Weinberg susceptibilities, theta ratio still noisy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000217,"raw_usage":{"total_tokens":1394,"prompt_tokens":863,"completion_tokens":531,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":479,"completion_tokens_details":{"reasoning_tokens":443}},"tokens_in":479,"tokens_out":531,"duration_ms":5442,"temperature":1.0,"reasoning_tokens":443,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T18:54:37.070258+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Produce a new ensemble with lattice spacing smaller than the smallest used here ($a<0.05$ fm), pion mass within a few MeV of the physical value, a larger spatial volume, and a near-physical kaon mass, then compare $\\chi_Q$ and $\\chi_W$ measured at flow times $\\tau_{\\mathrm{gf}} \\simeq 0.4$--$0.6$ fm with the error band of the paper's extrapolation; a mismatch beyond the combined uncertainties would falsify the extrapolation.","supporting_citations":[{"cited_title":"One-loop matching of the $CP$-odd three-gluon operator to the gradient flow","cited_arxiv_id":"2308.16221","evidence_quote":"Supplies the one-loop matching of the CP-odd three-gluon operator to the gradient-flow scheme, fixing the renormalization constants Z_W and Z_M used in the conversion matrix."},{"cited_title":"nEDM from the theta-term and chromoEDM operators","cited_arxiv_id":"2301.08161","evidence_quote":"Prior lattice determination of the topological susceptibility that the paper's (71 MeV)^4 result is compared with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the Peccei-Quinn mechanism that dynamically relaxes theta, the scenario in which the Weinberg-induced theta matters."},{"cited_title":"Hadron spectrum and lattice scale on 2+1-flavor lattice QCD ensembles","cited_arxiv_id":null,"evidence_quote":"Provides the ensemble parameters, lattice spacings, pion and kaon masses, and volumes, that enter the fits."}],"review_version":1}