{"id":"38d03ebe-141f-47d5-9fb6-6de8a30c65a7","arxiv_id":"2607.10606","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Using background electric fields, local topological charge, and non-Hermitian GEVP, lattice QCD yields dn = -0.0050(4)stat(8)sys θ-bar e fm at the physical point.","lead":"Lattice QCD calculation of the neutron's electric dipole moment induced by the QCD theta term, using a background electric field and local topological charge. The result constrains strong CP violation and improves statistical precision over prior form-factor methods.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"Chiral extrapolation from mπ ≥ 340 MeV with a two-parameter ansatz is the load-bearing step for the physical-point claim, and the quoted systematics do not cover its model dependence.","rationale":"The Reader correctly isolates the chiral extrapolation (together with the incomplete systematic budget) as the weakest assumption supporting the strongest claim. The lattice methodology—local topological charge + non-Hermitian GEVP isolation of the deformed ground state—is carefully cross-checked across operators, topological definitions and field strengths, so the three lattice points themselves are reliable. The load-bearing risk therefore sits entirely in the two-parameter fit that maps those points to the physical pion mass. Because the paper already flags conventional systematics as future work, the appropriate verdict remains CONDITIONAL; the concrete refit test simply quantifies how much of the present (8)sys is an under-estimate of model dependence. No stronger internal inconsistency is present.","tokens_in":49068,"tokens_out":591,"duration_ms":7697,"concrete_test":"Refit the three |nz|=1 points of Table III with the pure linear ansatz dn = c0 mπ^{2} (and, separately, with an added c2 mπ^{4} term). If the physical-point central value shifts by more than the present 0.0008 systematic, the quoted (8)sys under-covers the chiral-model uncertainty and the headline number must be re-qualified.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The headline number dn/θ-bar = -0.0050(4)stat(8)sys e fm (abstract, Eq. (80)) is obtained by fitting the three lattice points (mπ = 340, 420, 576 MeV) to the two-parameter form dn = c0 mπ^{2} + c1 mπ^{2} log(mπ^{2}/mN,phy^{2}) (Eq. (79), Sec. IV E, Fig. 25). The (8)sys is assembled only from Euclidean-time fit-window variation and the |nz|=1 vs |nz|=2 difference (Table III); the paper explicitly defers discretization, finite-volume and full chiral systematics. With only three points above 340 MeV and no higher-order chiral terms or continuum/volume data, the physical-point value is controlled by the assumed functional form rather than by the lattice matrix elements themselves. The multi-operator GEVP and local-topology improvements are solid and make the lattice points trustworthy, but they do not protect the extrapolation step that produces the quoted physical result.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper computes the QCD θ-term contribution to the neutron EDM on three 2+1-flavor domain-wall ensembles at fixed a≈0.11 fm (mπ=340, 420, 576 MeV) by measuring the energy shift of the neutron in a uniform Euclidean background electric field. Motivated by the Feynman–Hellmann theorem, the authors replace the noisy global topological charge with the forward matrix element of the local topological charge density between ground-state nucleons that have been deformed by the electric field; those states are isolated by a non-Hermitian GEVP that mixes positive- and negative-parity interpolators. After multi-operator consistency checks, ABJ-anomaly validation of two topological-charge definitions, and a two-parameter chiral extrapolation, they quote dn=−0.0050(4)stat(8)sys θ-bar e fm at the physical point, with the quoted systematic covering only Euclidean-time fit windows and |nz|=1 versus |nz|=2 differences. Conventional discretization, finite-volume and full chiral systematics are deferred.","tokens_in":49378,"tokens_out":1148,"duration_ms":30092,"significance":"If the physical-point number survives a complete error budget, it supplies one of the more precise lattice determinations of the θ-induced nEDM and tightens the bound on the strong-CP angle. Independently of the final number, the methodological advances—local rather than global topology, non-Hermitian GEVP isolation of the parity-mixed ground state, and explicit demonstration that conventional positive-parity operators suffer large, sign-changing contamination—are substantial and immediately reusable for other CP-odd operators (Weinberg three-gluon, four-quark, etc.). The multi-operator and ABJ cross-checks give high that the lattice matrix elements themselves are under control.","major_comments":[{"comment":"Sec. IV E, Eqs. (79)–(80), Fig. 25 and Table III: the headline physical-point value and its (8)sys are obtained from a two-parameter fit of three points (mπ≥340 MeV) to the chiral form dn=c0 mπ^{2} + c1 mπ^{2} log(mπ^{2}/mN,phy^{2}). The quoted systematic contains only fit-window and |nz| variation; ansatz dependence, higher-order chiral terms and the large lever arm to the physical point are explicitly deferred. Because the abstract and Eq. (80) present this number as the principal result, the model dependence of the extrapolation must be quantified (e.g., pure mπ^{2}, inclusion of mπ^{4} or Nπ continuum terms) and folded into the error, or the physical value must be demoted to a secondary, illustrative extrapolation while the three lattice points are emphasized.","section":null},{"comment":"Sec. III and the discussion surrounding Eq. (42): the calculation is electro-quenched (background field couples only to valence quarks). While the authors cite a ~1 % sea contribution to the magnetic moment, the same argument does not automatically apply to a CP-odd, topology-sensitive matrix element. An estimate or a controlled test of the quenching error is needed before the physical number can be regarded as complete.","section":null},{"comment":"Throughout Sec. IV and the abstract: only a single lattice spacing is used. Domain-wall fermions are automatically O(a)-improved, yet residual O(a^{2}) and finite-volume effects remain uncontrolled for the quoted physical result. The paper correctly flags these as future work, but the present error budget therefore understates the total uncertainty attached to dn/θ-bar at the physical point.","section":null}],"minor_comments":[{"comment":"Fig. 7 and the accompanying discussion of κest: the large operator dependence and lack of plateau are used to motivate the GEVP, but a short quantitative estimate of the polarizability bias (Appendix D) already appears in the text; cross-referencing it more prominently would help the reader.","section":null},{"comment":"Notation for Euclidean versus Minkowski quantities (Appendix A) is careful, yet the main text occasionally drops the E subscript after declaring it the default; a single clarifying sentence at the start of Sec. II would remove any residual ambiguity.","section":null},{"comment":"Table III and Fig. 24: the choice of fit windows for the two |nz|=2 cases that yield poorer χ^{2} is explained in the text, but the table caption itself should note which windows were adopted so that the numbers can be read without hunting through the prose.","section":null},{"comment":"Several multi-panel figures (e.g., Figs. 13–15, 22–23) pack four operators × three ensembles; increasing the font size of the operator labels or adding a common legend would improve readability.","section":null}],"recommendation":"major_revision","confidential_remarks":"The methodological core is strong enough that, once the chiral and continuum systematics are under better control, the paper will be a high-impact contribution. The present tension with some form-factor results (Fig. 25) is interesting and should be discussed more explicitly once the error budgets are comparable; I would not hold the paper for that discussion, only for a more honest treatment of the extrapolation uncertainty that currently dominates the physical claim."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The real news here is the method, not the final number. They extract the θ-induced nEDM from the forward matrix element of single-time-slice topological charge between nucleon ground states deformed by a Euclidean background electric field, isolating those states with a non-Hermitian GEVP that mixes positive- and negative-parity operators. That combination is new, and it works: statistical noise drops dramatically relative to global Qtop, and the multi-operator GEVP removes the large excited-state contamination that flips the sign of the usual positive-parity estimator. Cross-checks with gluonic vs ABJ topology (after axial correction), four interpolators, two field strengths, and several fit windows all line up. The lattice points themselves look trustworthy.\n\nWhat they do well is the technical control of the ground-state isolation and the honest documentation of every intermediate ratio. The appendices on anti-Hermitian perturbation theory and polarizability bias are useful. The comparison to Dragos, Alexandrou, Bhattacharya and Liang is fair and shows why their points sit higher than some earlier form-factor results.\n\nThe soft spot is exactly where the stress-test says: the headline dn/θ-bar = -0.0050(4)stat(8)sys is produced by fitting three points (340, 420, 576 MeV) to the two-parameter chiral form of Crewther et al. The (8)sys only covers Euclidean-time windows and |nz|=1 vs 2; discretization, volume and full chiral model dependence are deferred. With only one spacing and no higher-order terms, the physical-point value is still ansatz-driven. That does not invalidate the lattice matrix elements, but it does mean the quoted number is not yet a final hadronic input.\n\nThis is for lattice practitioners who care about CP-odd matrix elements and for phenomenologists who need to know the current status of the θ contribution. It deserves a serious referee. I would cite the method and the intermediate points; I would not yet treat the physical-point central value as definitive. Engage.","headline":"Solid methodological advance on θ-induced nEDM with local topology + non-Hermitian GEVP; the physical-point number is still controlled by a two-parameter chiral fit from mπ ≥ 340 MeV that the quoted systematics do not cover.","tokens_in":49964,"tokens_out":549,"would_cite":true,"duration_ms":8340,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"A lattice calculation extracts the theta-induced neutron EDM as -0.0050(4)(8) e fm by sampling local topology on a deformed nucleon ground state.","keywords":["neutron electric dipole moment","lattice QCD","theta term","background electric field","topological charge","generalized eigenvalue problem","CP violation","domain-wall fermions"],"falsifier":"A calculation on the same ensembles that isolates the ground-state matrix element with an independent multi-state fit (or a larger GEVP basis) and finds a central value outside the quoted -0.0050(4)(8) band, or a continuum-limit study that moves the physical-point result by more than the present systematic.","tokens_in":50003,"feed_emoji":"⚛️","tokens_out":963,"duration_ms":10940,"temperature":0.7,"pith_summary":"Permanent electric dipole moments of the neutron are a precision probe of CP violation. The only CP-odd operator of dimension four allowed by QCD is the theta term, but its contribution to the neutron EDM has been hard to compute cleanly on the lattice: global topological charge produces large noise, and ordinary positive-parity nucleon operators suffer severe excited-state contamination once a background electric field is turned on. This paper shows that the EDM can be read off as the forward matrix element of the local (single-time-slice) topological charge density between ground-state nucleons that have been deformed by a weak Euclidean electric field. The deformed ground state is isolated by a non-Hermitian generalized eigenvalue problem that mixes positive- and negative-parity operators. After a chiral extrapolation the authors report dn = -0.0050(4)stat(8)sys theta-bar e fm. The result is free of the form-factor extrapolation that earlier methods required and is consistent across several choices of interpolating field once the ground state is properly projected.","feed_headline":"Lattice QCD pins neutron EDM at -0.0050 theta e fm","feed_subtitle":"Local topology plus a deformed-nucleon GEVP yields a clean physical-point result","key_machinery":"Non-Hermitian generalized eigenvalue problem (GEVP) for the 2\times2 (or multi-operator) correlator matrix of positive- and negative-parity nucleon operators in a constant Euclidean electric field, combined with Feynman–Hellmann sampling of the single-time-slice topological charge density.","core_discovery":"The theta-induced neutron EDM is equal to the forward matrix element of the local topological charge density evaluated between the left and right ground-state nucleon eigenvectors that solve a non-Hermitian GEVP in a background Euclidean electric field. After multi-operator GEVP isolation and chiral extrapolation the authors obtain dn/theta-bar = -0.0050(4)stat(8)sys e fm at the physical point.","pith_inferences":["If residual Nπ contamination still hides inside the present GEVP, the physical-point central value could shift once true multi-hadron operators are added to the basis.","The method’s statistical power on local topology suggests it could become the default route for lattice EDM calculations once continuum and volume systematics are under control.","A parallel calculation of the proton EDM would immediately test whether the same deformed-ground-state machinery survives acceleration effects."],"forward_implications":["The same local-topology + non-Hermitian GEVP method can be applied without modification to other CP-odd operators (Weinberg three-gluon, four-quark, chromo-EDM).","Once lighter pion masses and finer lattices become available, the remaining conventional systematics can be quantified and the physical-point error budget reduced.","Form-factor extrapolations and global-topology noise can be bypassed for any CP-odd nucleon observable that couples to a local density.","Consistency between gluonic and (ABJ-corrected) fermionic definitions of topology supplies an internal cross-check for future lattice EDM programs."],"fun_headline_variants":["Lattice QCD sets nEDM to -0.0050 θ e fm via background E-field GEVP","Local topology plus deformed-nucleon GEVP yields nEDM -0.0050 θ e fm","Electric-field matrix element of top charge pins nEDM at -0.0050 θ e fm","Non-Hermitian GEVP extracts physical nEDM of -0.0050(4)(8) e fm","θ-induced neutron EDM equals -0.0050 e fm after lattice extrapolation"],"cache_read_input_tokens":32768,"weakest_assumption_plain":"The chiral extrapolation formula and the claim that residual excited-state and higher-order electric-field effects are fully captured by the quoted systematic error, even though only one lattice spacing is used and conventional discretization and finite-volume uncertainties are left for later work.","fun_headline_variants_meta":{"raw":{"variants":["Lattice QCD sets nEDM to -0.0050 θ e fm via background E-field GEVP","Local topology plus deformed-nucleon GEVP yields nEDM -0.0050 θ e fm","Electric-field matrix element of top charge pins nEDM at -0.0050 θ e fm","Non-Hermitian GEVP extracts physical nEDM of -0.0050(4)(8) e fm","θ-induced neutron EDM equals -0.0050 e fm after lattice extrapolation"]},"model":"grok-4.5","effort":"low","cost_usd":0.00756,"raw_usage":{"total_tokens":1944,"prompt_tokens":930,"num_sources_used":0,"completion_tokens":134,"cost_in_usd_ticks":75600000,"prompt_tokens_details":{"text_tokens":930,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":880,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":930,"tokens_out":134,"duration_ms":8434,"temperature":1.0,"reasoning_tokens":880,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T10:30:53.694009+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A calculation on the same ensembles that isolates the ground-state matrix element with an independent multi-state fit (or a larger GEVP basis) and finds a central value outside the quoted -0.0050(4)(8) band, or a continuum-limit study that moves the physical-point result by more than the present systematic.","supporting_citations":[],"review_version":1}