{"id":"847f73de-0872-4625-886b-1b1837a7ac49","arxiv_id":"1908.04151","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"The authors use two relativistic coupled-cluster methods to compute xenon-129's EDM sensitivity to nuclear Schiff moments and electron-nucleus tensor-pseudotensor interactions, with claimed uncertainties below 1%.","lead":"The paper reports improved theoretical calculations of how a xenon-129 atom responds to two time-reversal-violating effects that could generate an electric dipole moment. These response coefficients will be needed to interpret future xenon EDM experiments searching for physics beyond the Standard Model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed 0.2% (T-PT) and 0.7% (NSM) uncertainties are inconsistent with the 2% RCCSD(SC)-RNCCSD difference in eta and the 6% difference from the earlier RCCSD zeta, so the central accuracy claim is not yet supported.","rationale":"The paper is a serious precision calculation with genuine independent support: the alpha_d values from both RCCSD(SC) and RNCCSD are close to experiment, and the DF and CPDF results match previous calculations. That support is real but not decisive for the EDM coefficients, because alpha_d is an electric-dipole response of the same rank and parity, not a direct measure of P,T-odd operator convergence. The most direct evidence about the accuracy claim is internal: Table I contains related many-body results whose spread is 2% for eta (0.50 from RCCSD, 0.49 from RNCCSD, 0.48 from RCCSD(SC)) and 6% for zeta (0.34 from previous RCCSD versus 0.32 here). The stated errors (0.2%, 0.7%) are smaller than the smallest inter-method difference. The authors assert rapid convergence, but the only quantitative error indicators they provide, partial triples and Breit, contribute far less than the observed spread. This makes the central claim of accuracy to 0.2%/0.7% unsubstantiated. However, the numerical coefficients themselves are close to previous work, so the paper remains useful for extracting limits if the error bars are revised upward. Thus the appropriate verdict is conditional acceptance, which is unchanged from the reader's verdict.","tokens_in":9468,"tokens_out":7001,"duration_ms":71027,"concrete_test":"Recompute eta and zeta with the present RNCCSD code using the GTO basis of Ref. [28] and with the Breit interaction included in the T1/T2 amplitude equations. If either eta deviates from 0.49 by more than 0.001 or zeta deviates from 0.32 by more than 0.002, the claimed 0.2%/0.7% uncertainties are contradicted by basis/method sensitivity.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing concern is the uncertainty budget. The key result is that RNCCSD yields eta = 0.49 and zeta = 0.32 with estimated errors of 0.2% and 0.7%, obtained by adding partial-triples corrections (3.9e-5 in eta, 1.3e-4 in zeta) and Breit corrections (1.1e-3 in eta, 2.1e-3 in zeta). But Table I shows RCCSD(SC) gives eta = 0.48, a 2% difference from RNCCSD, and the earlier RCCSD calculation of Ref. [28] gives zeta = 0.34, about 6% above the adopted 0.32. These method/basis spreads are an order of magnitude larger than the quoted errors. The paragraph beginning 'We have evaluated the numerical error...' and the final conclusion treat triples and Breit as if they exhausted the omitted effects, but no calculation is shown demonstrating that the RCCSD(SC)-RNCCSD or prior-RCCSD differences are explained by these corrections. Since the extraction of limits on S and C_T uses these coefficients directly, the derived limits |S| < 4.7e-10 and |C_T| < 6.1e-7 inherit the underestimated uncertainty.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports relativistic coupled-cluster calculations of the 129Xe EDM sensitivity coefficients eta = d_a/(<sigma_N> C_T 10^20 |e| cm) and zeta = d_a/(S 10^17 |e| cm/(|e| fm^3)), using both the self-consistent RCCSD method and the normal RCCSD (RNCCSD) method with Gaussian-type orbitals. It reports RNCCSD values eta = 0.49 and zeta = 0.32, adds corrections for the Breit interaction and partial triple excitations, and estimates final uncertainties of 0.2% for the T-PT channel and 0.7% for the NSM channel. Combining these coefficients with the recent 129Xe EDM limit from Ref. [19] yields |C_T| < 6.1e-7 and |S| < 4.7e-10 |e| fm^3. The paper also highlights the good agreement of its computed electric dipole polarizability with experiment and argues that 129Xe is a particularly promising EDM probe.","tokens_in":9802,"tokens_out":3551,"duration_ms":37931,"significance":"If the quoted accuracy is justified, the calculated coefficients are valuable: they are ab initio in the sense that no EDM data are used to adjust parameters, and the polarizability benchmark from Ref. [29] provides an independent check of the many-body treatment. The comparison between RCCSD(SC) and RNCCSD, and the explicit tabulation of leading correlation terms, are also useful. However, the central claim is the quoted 0.2% and 0.7% uncertainties, because the limits in Eqs. (25)-(26) scale directly with eta and zeta. The manuscript does not currently reconcile those error bars with the 2% method spread for eta and the 6% difference from the earlier RCCSD zeta in Table I.","major_comments":[{"comment":"The central accuracy claim is not supported by the spread of methods reported in Table I. RCCSD(SC) gives eta = 0.48 versus RNCCSD eta = 0.49, a 2% difference, while the earlier RCCSD result zeta = 0.34 from Ref. [28] differs from the adopted RNCCSD value zeta = 0.32 by about 6%. The error paragraph attributes the omitted effects entirely to partial triples and the Breit interaction, but no calculation is shown demonstrating that these two corrections account for either of those differences. The quoted uncertainties of 0.2% and 0.7% are roughly an order of magnitude smaller than the observed method/basis spread, so the derived limits in Eqs. (25) and (26) inherit an underestimation of uncertainty. The authors should either demonstrate convergence across the two methods with a detailed accounting of the RCCSD(SC)-RNCCSD and Ref. [28] differences, or quote more conservative uncertainties that reflect the actual spread.","section":"Table I and the paragraph beginning 'We have evaluated the numerical error in our RCC calculations...'"},{"comment":"The text states that the Breit contributions were 0.6% and 0.9% of the total Dirac-Coulomb contributions in the CPDF and RCCSD approximations, respectively, but Table I lists only Delta_CPDF^Breit = -0.001 for eta and -0.002 for zeta; the RCCSD-level Breit values are neither tabulated nor sufficiently described. Since the Breit estimate enters the final error budget directly, the missing values and the definition of the percentages should be provided, and the relation between these percentages and the final absolute uncertainties in eta and zeta should be made explicit.","section":"Paragraph beginning 'In this work, the Breit interaction contributions were found to be...'"},{"comment":"The paper uses the quoted uncertainties to assert that the results are 'more accurate and reliable' than those for 199Hg and to argue that 129Xe is the most promising stable diamagnetic EDM probe. This conclusion is load-bearing for the paper's stated implications. Because the uncertainty budget is not yet reconciled with the method spread, the statement that the estimated errors are 0.2% and 0.7% is premature, and the derived bounds on S and C_T should be presented with error bars that reflect the unresolved spread until the convergence question is settled.","section":"Concluding discussion and Eqs. (25)-(26)"}],"minor_comments":[{"comment":"There is a typo in the title as displayed: 'Standard Mo del' should be 'Standard Model'; also 'Schiff moment' appears with a nonstandard spacing in several places.","section":"Title and abstract"},{"comment":"The triples corrections are listed as '~ 0' in Table I, while the text reports absolute contributions of 3.9e-5 for eta and 1.3e-4 for zeta; these numbers should be included in the table and their relative sizes compared explicitly with the final errors.","section":"Table I"},{"comment":"The notation in Eq. (19), '<Phi_0(1 + ~T)O|Phi_0>_c', is missing the bra bar and should be written as '<Phi_0|(1 + ~T)O|Phi_0>_c' for consistency with the other equations.","section":"Equation (19)"},{"comment":"The phrase 'improving by factors of one-and-half and five times' should read 'by factors of 1.5 and 5'.","section":"Introduction"},{"comment":"The caption refers to 'RCCSD(SC)' but the table lists terms with 'h.c.'; it would be helpful to state explicitly which terms are the hermitian conjugates and whether the listed numerical values include them or not.","section":"Table II caption"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this is a genuine new calculation: the authors apply RCCSD(SC) and RNCCSD to the 129Xe EDM and give a transparent term-by-term breakdown. The central coefficients, eta=0.49, zeta=0.32, sit close to earlier values, which is reassuring. Second, don't take the error bars at face value. The 0.2% and 0.7% estimates come from adding partial triples and Breit corrections only, and they are an order of magnitude smaller than the internal spread: RCCSD(SC) gives eta=0.48, about 2% below the RNCCSD value, and the earlier RCCSD of Ref. [28] gives zeta=0.34, about 6% above the adopted 0.32. That doesn't make the calculation wrong, but it does make the precision claim as written unsupported.\n\nWhat is genuinely good: the work is methodologically careful. The equations follow standard coupled-cluster practice, the DF and CPDF results reproduce prior numbers, and the polarizability benchmark (27.51 vs experimental 27.815) is an independent check on the basis and correlation treatment. No EDM data or fitted constants enter the calculation; the coefficients are ab initio. The term decompositions in Tables II and III are useful and allow readers to see where the physics sits.\n\nThe soft spot is exactly the error budget. The paragraph beginning 'We have evaluated the numerical error...' treats triples and Breit as if they exhausted the omitted effects, but no calculation shows that these corrections account for the 2% eta gap or the 6% zeta gap. The earlier RCCSD value of 0.34 is described as being in good agreement, which is fair at the 6% level, but that sits in tension with a claimed 0.7% uncertainty. Part of the spread may be basis-set incompleteness, since the earlier calculations used different GTO bases, but the paper does not estimate that. The derived limits (|S|<4.7e-10, |C_T|<6.1e-7) inherit the understated uncertainty.\n\nFor a reader in the EDM subfield, this is a useful independent determination of Xe sensitivity coefficients, and the broader case that Xe could rival Hg is reasonable, though not new. But the headline accuracy should be revised from sub-percent to a few percent unless the authors can reconcile the method spread. I would send this to peer review; it's a credible calculation with a fixable but real discrepancy in the error analysis.","headline":"A credible, careful RCC calculation of 129Xe EDM coefficients whose central values are fine, but the claimed 0.2%/0.7% uncertainties are not supported by the 2-6% spread across methods and earlier work.","tokens_in":10286,"tokens_out":3657,"would_cite":true,"duration_ms":32691,"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":"The 129Xe EDM sensitivity coefficients are computed to sub-percent accuracy.","keywords":["xenon-129","electric dipole moment","Schiff moment","tensor-pseudotensor interaction","relativistic coupled cluster","CP violation","time-reversal violation","diamagnetic atom"],"falsifier":"Perform an independent all-order calculation with a different basis set that includes all connected triple excitations and the Breit interaction; if $\\eta$ moves by more than about 0.001 (0.2% of 0.49) or $\\zeta$ by more than about 0.002 (0.7% of 0.32), the claimed error budget is incomplete.","tokens_in":9304,"feed_emoji":"⚛️","tokens_out":9766,"duration_ms":90983,"temperature":0.7,"pith_summary":"The paper aims to establish that the electric dipole moment (EDM) of 129Xe can be calculated accurately enough from relativistic many-body theory that future measurements become sharp tests of time-reversal and parity violation. It reports sensitivity coefficients for the two dominant sources—the nuclear Schiff moment and the electron–nucleus tensor-pseudotensor interaction—using self-consistent and normal coupled-cluster variants, and estimates the remaining theoretical error at 0.7% and 0.2%. If those uncertainties hold, combining the coefficients with improved Xe EDM measurements would tighten limits on hadronic CP violation and help identify new physics beyond the Standard Model. The authors also argue that 129Xe, unlike 199Hg, has both reliable atomic and nuclear theory, making it a leading candidate among stable diamagnetic atoms.","feed_headline":"Xenon-129 EDM theory sharpened to sub-percent precision","feed_subtitle":"At 0.2–0.7% uncertainty, the new Xe-129 sensitivity coefficients would turn future EDM runs into sharper tests of CP violation.","key_machinery":"The load-bearing object is the relativistic normal coupled-cluster singles-and-doubles (RNCCSD) expectation value, in which the bra state is replaced by $\\langle\\widetilde{\\Psi}| = \\langle\\Phi_0|(1+\\widetilde{T})e^{-T}$, so the EDM expression terminates and is stationary with respect to the bra amplitudes. This gives the finite sum $d_a/\\lambda = \\langle\\Phi_0|\\widetilde{T}^{(1)}H_{\\mathrm{PTV}} + (1+\\widetilde{T}^{(0)})H_{\\mathrm{PTV}}T^{(1)}|\\Phi_0\\rangle_c$, avoiding the non-terminating series of ordinary coupled cluster. The companion RCCSD(SC) calculation sums powers of $T^{(0)}$ and $T^{(0)\\dagger}$ self-consistently; agreement between the two variants is the paper's evidence that higher-order many-body effects have converged. The Gaussian-type orbital basis and the perturbed-triples and Breit-interaction error estimates come from the authors' earlier polarizability study of the same atom.","core_discovery":"On the paper's own terms, the central discovery is that the P,T-odd sensitivity of 129Xe is converged at the relativistic coupled-cluster singles-and-doubles level: the self-consistent and normal (RNCCSD) methods give $\\eta = d_a/(\\langle\\sigma_N\\rangle C_T) \\times 10^{20}\\,|e|\\mathrm{cm}$ values of 0.48 and 0.49, respectively, and identical $\\zeta = d_a/S \\times 10^{17}\\,|e|\\mathrm{cm}/(|e|\\,\\mathrm{fm}^3)$ values of 0.32. Including estimated corrections from partial triple excitations and the Breit interaction, the paper quotes $d_a = 0.49\\times10^{-20}\\langle\\sigma_N\\rangle C_T\\,e\\,\\mathrm{cm}$ for the T-PT channel and $d_a = 0.32\\times10^{-17} S/(|e|\\,\\mathrm{fm}^3)\\,|e|\\,\\mathrm{cm}$ for the Schiff-moment channel, with estimated errors of 0.2% and 0.7%. Combining these with the current experimental limit $|d_a|<1.5\\times10^{-27}\\,|e|\\mathrm{cm}$, and assuming a single source, yields $|S|<4.7\\times10^{-10}\\,|e|\\mathrm{fm}^3$ and $|C_T|<6.1\\times10^{-7}$.","pith_inferences":["The quoted 0.2% and 0.7% are atomic-structure errors; the nuclear Schiff-moment calculation carries its own uncertainty, so the total physics reach may be set by nuclear theory rather than these atomic coefficients.","The same terminating expectation-value machinery, already validated on the electric dipole polarizability, could be applied to other noble-gas EDMs to cross-check systematics.","The single-source bounds from the paper would need a two-parameter analysis if both T-PT and Schiff contributions contribute comparably to a future measured Xe EDM."],"forward_implications":["If the computed coefficients are correct, an improved 129Xe EDM measurement directly bounds the tensor-pseudotensor coupling $C_T$ and the nuclear Schiff moment, assuming one dominant source.","Planned 129Xe runs aiming for sensitivity near $10^{-30}\\,e\\,\\mathrm{cm}$ would bring xenon's reach to the level of, or beyond, the current mercury limit.","The small 2% spread between RCCSD(SC) and RNCCSD for the T-PT channel, and exact agreement for the Schiff channel, indicates rapid convergence of higher-order many-body effects.","With QCD input, sharper Schiff-moment and T-PT limits translate into tighter constraints on the QCD theta term and quark chromo-EDMs."],"supporting_citations":[{"why":"Introduces the relativistic normal coupled-cluster method whose terminating expectation value expression is used for the EDM.","marker":"[26]"},{"why":"Earlier RCCSD calculation of the 129Xe EDM that supplies the baseline values extended here and the comparison with RCCSD(SC).","marker":"[28]"},{"why":"Prior polarizability study of the same atom that provides the Gaussian-type orbital basis and the procedures for estimating triples and Breit corrections.","marker":"[29]"},{"why":"Supplies the CPDF comparison values and the forms of the Schiff-moment and T-PT Hamiltonians used in the calculation.","marker":"[32]"},{"why":"Latest 129Xe EDM experimental limit, combined with the computed coefficients to derive the bounds on the Schiff moment and $C_T$.","marker":"[19]"},{"why":"Best current mercury EDM limit, used as the benchmark for sensitivity that the authors argue xenon can match or surpass.","marker":"[10]"},{"why":"199Hg EDM calculation whose larger many-body spread is contrasted with the rapidly converging 129Xe results.","marker":"[45]"},{"why":"Nuclear Schiff-moment calculation that underpins the claim that nuclear theory for 129Xe is more reliable than for 199Hg.","marker":"[6]"}],"fun_headline_variants":["Sub-percent Xe-129 EDM theory tightens CP violation constraints","Xe-129 EDM theory narrows uncertainty to 0.2-0.7% for new physics","Sharper Xe-129 EDM: 0.7% Schiff-moment accuracy for CP tests","Relativistic many-body theory refines Xe-129 EDM to sub-percent"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that partial triple excitations plus the Breit interaction account for all omitted correlation and relativity, so the RNCCSD results are converged to 0.2% and 0.7%; this is asserted even though the two coupled-cluster variants differ by 2% for the T-PT coefficient.","fun_headline_variants_meta":{"raw":{"variants":["Sub-percent Xe-129 EDM theory tightens CP violation constraints","Xe-129 EDM theory narrows uncertainty to 0.2-0.7% for new physics","Sharper Xe-129 EDM: 0.7% Schiff-moment accuracy for CP tests","Relativistic many-body theory refines Xe-129 EDM to sub-percent"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00083,"raw_usage":{"total_tokens":3675,"prompt_tokens":1047,"completion_tokens":2628,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":663,"completion_tokens_details":{"reasoning_tokens":2528}},"tokens_in":663,"tokens_out":2628,"duration_ms":18806,"temperature":1.0,"reasoning_tokens":2528,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:49:59.603964+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform an independent all-order calculation with a different basis set that includes all connected triple excitations and the Breit interaction; if $\\eta$ moves by more than about 0.001 (0.2% of 0.49) or $\\zeta$ by more than about 0.002 (0.7% of 0.32), the claimed error budget is incomplete.","supporting_citations":[{"cited_title":"Flambaum, I.B","cited_arxiv_id":null,"evidence_quote":"Introduces the relativistic normal coupled-cluster method whose terminating expectation value expression is used for the EDM."},{"cited_title":"9 × 10−5 and 1","cited_arxiv_id":null,"evidence_quote":"Prior polarizability study of the same atom that provides the Gaussian-type orbital basis and the procedures for estimating triples and Breit corrections."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the CPDF comparison values and the forms of the Schiff-moment and T-PT Hamiltonians used in the calculation."},{"cited_title":"Kuchler et al., Hyperfine Interact","cited_arxiv_id":null,"evidence_quote":"Latest 129Xe EDM experimental limit, combined with the computed coefficients to derive the bounds on the Schiff moment and $C_T$."},{"cited_title":"Chupp, M","cited_arxiv_id":null,"evidence_quote":"199Hg EDM calculation whose larger many-body spread is contrasted with the rapidly converging 129Xe results."},{"cited_title":"Results of this calculation are of the same sign and of the same order of magnitude as the previous calculation, unlike the case of 199Hg [8, 11]","cited_arxiv_id":null,"evidence_quote":"Nuclear Schiff-moment calculation that underpins the claim that nuclear theory for 129Xe is more reliable than for 199Hg."}],"review_version":1}