{"id":"f9579174-c550-4ccd-a7d8-b983628b4ee8","arxiv_id":"2502.06225","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"General-order relativistic coupled-cluster calculations with triple excitations yield electric-field response properties for Ca+ and Yb+ clock states, including new quadrupole polarizabilities.","lead":"This paper computes how the clock states of calcium and ytterbium ions respond to electric fields, giving values for polarizabilities and quadrupole moments used to correct optical clock errors. The new calculations include triple-electron correlations, which the authors find are key to matching experiments for the ytterbium clock.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Final Yb+ 2D3/2 values rest on an untested 2ζ-triples additivity; the paper's own ground-state data show the triple correction is strongly basis-dependent.","rationale":"The reader's weakest assumption is the load-bearing issue, and the arithmetic sharpens it: the procedure actually used for Table III is not the additivity described in Section III, and the two differ by 9.25 a.u. for αS_d. Since the summary's benchmark claims all rest on these numbers, this is central. The proposed 3ζ SDT calculation is feasible—the authors already performed it for the Yb+ ground state—and it directly tests whether the 2ζ triple correction is transferable. Given the paper's own ground-state data show a factor-of-four basis dependence for the triple correction, this test is mandatory before the recommended values are used as clock-systematic benchmarks. Encouraging agreement with experiment and the useful Ca+ analysis do not substitute for direct validation in Yb+. This does not move the reader's CONDITIONAL verdict; it specifies the condition that must be met.","tokens_in":16752,"tokens_out":13604,"duration_ms":111268,"concrete_test":"Run an e23-SDT<10 RCCSDT calculation for the 4f14 5d 2D3/2 state of Yb+ with the dyall.cv3z basis, the same calculation already reported for the Yb+ ground state in Table II. Form the final estimate as αS_d(4ζ, SD) + [αS_d(3ζ, SDT) − αS_d(3ζ, SD)], and the analogous estimates for ΔαS_d and Θ. If these values differ from the quoted 101(4), 38(4), and 1.973(31) by more than the reported uncertainties, the 2ζ-based additivity fails and the error budget must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section IV's final D3/2 values do not implement the additivity announced in Section III. Section III says to take large-basis RCCSD and add (RCCSDT−RCCSD) from a small basis; for Table III that would give αS_d = 108.36 + (103.31 − 119.71) = 91.96 a.u. Instead the text takes the 2ζ RCCSDT value 103.31 and adds the SD-only 3ζ→4ζ basis shift (108.36−110.46 = −2.10), yielding 101.21 a.u., algebraically equivalent to adding (SDT_2ζ − SD_3ζ) = −7.15 to SD_4ζ. The two procedures differ by 9.25 a.u., far outside the quoted 4 a.u. uncertainty. The paper's own ground-state data (Table II) show triple corrections are strongly basis-dependent: ΔPT(2ζ)=−2.27 vs ΔPT(3ζ)=−0.53 for αS_d(Yb+), a fourfold change. No RCCSDT result exists for the 2D3/2 state beyond 2ζ, so the central benchmark values (αS_d=101(4), ΔαS_d=38(4), and by the same scheme Θ=1.973(31)) rest on an untested mixed-basis additivity. The quoted uncertainties are formed from incremental corrections, not from this additivity assumption, so their coverage is unknown.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports general-order relativistic coupled-cluster (RCC) calculations, at the RCCSD and RCCSDT levels, of static electric dipole polarizabilities, quadrupole moments, and quadrupole polarizabilities for clock-relevant states of Ca+ and Yb+. A finite-field approach is used to extract these properties from energy shifts. The authors present recommended values for Yb+ (e.g., αS_d = 101(4) a.u. for the 5d 2D_3/2 state, ΔαS_d = 38(4) a.u. for the 2S_1/2 → 2D_3/2 transition, and Θ = 1.973(31) a.u. for the 2D_3/2 state) and argue that triple excitations are decisive in bringing the calculated values into agreement with recent experiments. The Ca+ calculations are used as a lighter proxy to validate the role of triples and to guide basis-set choices for Yb+.","tokens_in":16943,"tokens_out":7411,"duration_ms":58158,"significance":"If correct, the recommended values would be a valuable ab initio benchmark for systematic-effect corrections in Ca+ and Yb+ optical clocks. The work uses a genuine general-order RCC implementation with triples and a finite-field method that avoids property-evaluation issues of some earlier approaches; it is also free of experimental input in the Hamiltonian. The Ca+ analysis is a useful methodological demonstration. However, the final Yb+ numbers rely on a composite-scheme assumption that is not validated and that is implemented inconsistently with the scheme stated in Sec. III. Because the main claim is to provide benchmark-quality values with reliable uncertainties, this issue must be resolved before the numbers can be adopted with confidence.","major_comments":[{"comment":"The composite scheme actually used for the 2D_3/2 state differs from the scheme stated in Sec. III. Sec. III defines the final value as RCCSD(large basis) + [RCCSDT(small basis) − RCCSD(small basis)]. For Table III this would give αS_d = 108.36 + (103.31 − 119.71) = 91.96 a.u. Instead, the text combines the 2ζ RCCSDT value with the 3ζ→4ζ basis shift obtained from RCCSD: 103.31 + (108.36 − 110.46) = 101.21 a.u. The two prescriptions differ by 9.25 a.u., which is more than twice the quoted uncertainty of 4 a.u. This directly affects the recommended αS_d = 101(4) a.u. and ΔαS_d = 38(4) a.u., and therefore the claimed agreement with the experimental result of Ref. [13]. The authors must either reconcile the two schemes or justify why the second one is the correct composite estimate.","section":"Sec. IV, Table III"},{"comment":"The additivity assumption that triple-excitation corrections are nearly independent of basis size is contradicted by the paper's own ground-state data. For αS_d of Yb+, Table II gives ΔPT = −2.27 a.u. at 2ζ but −0.53 a.u. at 3ζ, a factor-of-four change. No RCCSDT result is reported for the 2D_3/2 state beyond 2ζ, so the triple correction in Table III is taken at the smallest basis, where basis error is largest. The uncertainty budget for the final 2D_3/2 values is built from Δbasis, Δ4d, Δvirt, and Δsaug, but it does not include the uncertainty in the additivity assumption itself. The error bars therefore have unknown coverage, and the recommended values could be biased beyond the quoted uncertainties.","section":"Sec. IV, Tables II and III"},{"comment":"For the 5d 2D_5/2 state, no RCCSDT calculation is reported. The final values are obtained from RCCSD only, with an uncertainty estimated by multiplying Δbasis by a factor of 2 'to be on safer side.' This is an ad hoc uncertainty that is not based on an actual treatment of triple excitations, and it is not justified by the Ca+ analysis, where triples contribute at the level of several percent to the 2D-state polarizabilities. As the paper recommends αS_d = 89(6) a.u., αT_d = −74(2) a.u., and Θ = 3.06(7) a.u. for this state, the reliability of these benchmark values is not established.","section":"Sec. IV, Table V"}],"minor_comments":[{"comment":"The discussion refers to 'Exp-2015' for Ref. [12], but Ref. [12] is Schneider et al. (2005); this appears to be a typographical error.","section":"Sec. IV, text near Table III"},{"comment":"The final Θ is quoted as 1.973(31) a.u., but adding the listed corrections (1.985 − 0.025 + 0.017 − 0.004 − 0.005) gives 1.968 a.u.; the manuscript should clarify how 1.973 is obtained.","section":"Sec. IV, Table IV"},{"comment":"The notation 'e32-CCSD' and '5D5/2' appears in the table and accompanying text; these should be 'e23-CCSD' and '2D_5/2' or '5d 2D_5/2'.","section":"Sec. IV, Table V and text"},{"comment":"The text states that the ground-state polarizability used in the final ΔαS_d is αS_d = 65.53(35) a.u., but the recommended ground-state value in Table II and the Summary is 63.53(35) a.u.; this appears to be a typographical error.","section":"Sec. IV, Table III"}],"recommendation":"major_revision","confidential_remarks":"The central issue is the unvalidated composite-scheme assumption. In my reading this is not a presentation problem: the two composite procedures described in the paper differ by about 9 a.u. for αS_d of the 2D_3/2 state, which is larger than the quoted uncertainty. If the authors can provide an RCCSDT result for the 2D_3/2 state in a larger basis (or otherwise demonstrate that the additivity is valid), the paper could become publishable after revision. If not, the recommended values and their uncertainties would need to be substantially revised or withdrawn."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Good to see this one. It's a serious ab initio effort: first finite-field RCCSDT calculation of dipole and quadrupole polarizabilities and quadrupole moments for Ca+ and Yb+ clock states, and the first alpha_q values for the Yb+ 2D3/2 state. The final Yb+ D3/2 numbers land close to recent measurements, and the Ca+ proxy study is a reasonable strategy.\n\nThe soft spot is in how the Yb+ D3/2 values are assembled. Section III describes a standard composite: RCCSD at the largest basis plus the (RCCSDT−RCCSD) correction from a smaller basis. For Table III, that gives 108.36 + (103.31 − 119.71) = 91.96 a.u. What the paper actually does is take the 2ζ RCCSDT value 103.31 and add an SD-level basis shift from 3ζ to 4ζ of −2.10, getting 101.21 a.u. The two procedures differ by 9.25 a.u., far outside the quoted 4 a.u. uncertainty, and the text doesn't acknowledge the switch. The same mismatch appears in Table II, where the final ground-state polarizability is shifted by Δbasis even though the text says those corrections go into the uncertainty budget.\n\nThe bigger worry is that the additivity is unvalidated for the target states. The paper's own ground-state data show the triple correction is strongly basis dependent: ΔPT is −2.27 a.u. at 2ζ and −0.53 a.u. at 3ζ for αS_d(Yb+). For the 2D3/2 state, no RCCSDT result exists beyond 2ζ, so the large triple correction (16.4 a.u.) could change significantly with basis. The quoted ±4 a.u. uncertainty is built from incremental SD-level corrections, not from the additivity assumption, so its coverage is unknown. And the statement that a prior RCCSD calculation by the same group cannot be reproduced due to 'implementation differences' is a reproducibility flag that needs an explanation.\n\nWhat I credit: the calculations are genuinely ab initio, no experimental fitting, and the experiment comparisons are retrospective validation. The work gives the clock community new numbers for alpha_q and a plausible resolution of the theory-experiment discrepancy for ΔαS_d of the E2 transition. The Ca+ results are on firmer ground because triples are included in larger bases there.\n\nThis deserves a serious referee, but not acceptance as written. The authors should either compute the 3ζ RCCSDT for the D3/2 state or justify why the 2ζ triple correction is adequate, reconcile the method description with the executed arithmetic, and explain the Ref. [35] discrepancy. As is, I'd treat the final Yb+ D3/2 values as preliminary rather than benchmarks.","headline":"Solid Ca+ results and plausible Yb+ numbers, but the Yb+ D3/2 benchmark rests on an untested additivity that contradicts the stated protocol.","tokens_in":17616,"tokens_out":9390,"would_cite":false,"duration_ms":66911,"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":"Relativistic coupled-cluster calculations with triple excitations fix the electric-field response properties of the Yb+ optical-clock states, yielding a differential polarizability of 38(4) a.u. in accord with measurement.","keywords":["relativistic coupled-cluster","dipole polarizability","quadrupole moment","optical clock","ytterbium ion","calcium ion","finite-field method","triple excitations"],"falsifier":"A full RCCSDT calculation in the 4ζ basis for the Yb+ 5d 2D3/2 state (or an equivalent large-basis treatment of triples) would settle whether the additivity assumption holds; if the computed Δαd shifts by more than the quoted 4 a.u., the recommended value is biased. Alternatively, a more precise experimental measurement of Δαd for the Yb+ clock transition with uncertainty below 2 a.u. would directly test the central value.","tokens_in":16368,"feed_emoji":"⚛️","tokens_out":5023,"duration_ms":40860,"temperature":0.7,"pith_summary":"This paper tries to establish that the electric-field response properties of the Yb+ ion's optical clock states—the differential scalar dipole polarizability, the quadrupole moment, and the quadrupole polarizability—can be computed accurately enough to serve as benchmarks for systematic-effect corrections in clocks. Using general-order relativistic coupled-cluster theory through singles, doubles, and triples (RCCSDT) in a finite-field scheme, the authors obtain a differential scalar dipole polarizability of 38(4) a.u. for the 6s 2S1/2 → 5d 2D3/2 clock transition, matching a 2018 measurement, and a quadrupole moment of 1.973(31) a.u. for the 5d 2D3/2 state, matching a 2020 measurement. The central lesson is that triple excitations, previously ignored in most RCC property calculations, are decisive for these quantities. They also validate the method on the lighter Ca+ ion, whose similar clock states allow a more complete correlation treatment.","feed_headline":"Triple excitations fix Yb+ clock polarizability at 38(4) a.u.","feed_subtitle":"Relativistic coupled-cluster results resolve a theory-experiment clash for the ytterbium ion's E2 clock transition.","key_machinery":"The central object is the general-order relativistic coupled-cluster expansion with single, double, and triple excitations (RCCSDT), applied through a finite-field approach: the energy of each clock state is computed under a series of static electric fields (and field gradients), and the dipole polarizability, quadrupole moment, and quadrupole polarizability are read off from the coefficients of the polynomial energy shift versus field strength. This bypasses direct property-evaluation expressions, and the authors argue it satisfies the Hellmann-Feynman theorem in the RCC framework. The strategy of first benchmarking on Ca+—whose similar 4s–3d clock transitions allow nearly complete correlation—is used to decide basis-set and triple-excitation handling before transferring the method to Yb+.","core_discovery":"The paper's central claim is that the previously discrepant theoretical and experimental values for the Yb+ E2 clock transition's electric-field response arise primarily from missing triple-excitation correlations. Computing energies under external electric fields and field gradients with RCCSDT wavefunctions, and extracting polarizabilities and quadrupole moments from polynomial fits of those energies, the authors obtain ground-state αd = 63.53(35) a.u., excited-state αd = 101(4) a.u. and αT = −72(5) a.u. for the 5d 2D3/2 state, differential Δαd = 38(4) a.u., and quadrupole moment Θ = 1.973(31) a.u. These numbers agree with the most recent experiments and support using them as benchmarks, while also explaining why older calculations and one of the two experimental groups sat systematically higher.","pith_inferences":["One could test whether the same RCCSDT-plus-basis-correction strategy applies to other heavy-ion clocks such as Sr+ or Ba+, where triple excitations may be underestimated as well.","A direct extension would be to compute the E3 clock state (4f13 6s2 2F7/2) polarizabilities with this method; the paper does not attempt that transition.","The finite-field extraction of αq from third-order energy terms may be sensitive to the fitting range, so future work should report stability across field strengths.","The explanation for older RCCSD results suggests that earlier sum-over-states analyses may need to revisit triple-excitation effects in their property evaluations."],"forward_implications":["The recommended values give clock-systematics corrections with uncertainties small enough to improve the Yb+ E2 clock systematic budget.","The result resolves the conflict between the two Yb+ differential-polarizability measurements, siding with the 2018 value.","The quadrupole moment agrees with recent experiment, implying the quadrupole shift can be reduced in Yb+ E2 clock operations.","The method provides first reported values of αq for the Yb+ clock states, useful for estimating gradient-induced shifts.","The improved Ca+ values also strengthen benchmarks for the Ca+ optical clock."],"supporting_citations":[{"why":"Provides the experimental differential scalar dipole polarizability (35.7(1.8) a.u.) that the central Yb+ result is compared against.","marker":"[13]"},{"why":"Cited as the recent experimental source of the quadrupole moment that the recommended Θ aligns with.","marker":"[8]"},{"why":"Reports the 2020 measured quadrupole moment of 1.95(1) a.u. used as the primary comparison for Θ in the full text.","marker":"[14]"},{"why":"Previous RCCSD calculation of ground and excited state polarizabilities that the paper's larger-basis RCCSDT results supersede.","marker":"[19]"},{"why":"Earlier RCC calculation of the quadrupole moment that overestimates Θ relative to the recent experiment.","marker":"[22]"},{"why":"Fock-space coupled-cluster calculation of Θ used as a benchmark for the 5d 2D3/2 and 5d 2D5/2 states.","marker":"[23]"},{"why":"Configuration-interaction plus core-polarization calculation providing a recent alternative polarizability value for comparison.","marker":"[24]"},{"why":"Earlier experimental values for both the differential polarizability and quadrupole moment that show the previous discrepancy with theory.","marker":"[12]"},{"why":"Previous RCCSD calculation in Ca+ whose values the paper reproduces only after adding triple-excitation corrections.","marker":"[35]"}],"fun_headline_variants":["Triple excitations resolve Yb+ clock polarizability clash","RCCSDT yields Yb+ polarizability 38(4) a.u. matching expt","Coupled-cluster triples settle Yb+ quadrupole moment","General-order RCC fixes Yb+ electric-field response"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The final Yb+ values rely on the assumption that triple-excitation contributions computed with a smaller (2ζ or 3ζ) basis can be added as a correction to RCCSD results from a larger (4ζ) basis; if triple corrections depend strongly on basis size, the recommended values could be biased by more than the quoted uncertainties.","fun_headline_variants_meta":{"raw":{"variants":["Triple excitations resolve Yb+ clock polarizability clash","RCCSDT yields Yb+ polarizability 38(4) a.u. matching expt","Coupled-cluster triples settle Yb+ quadrupole moment","General-order RCC fixes Yb+ electric-field response"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000444,"raw_usage":{"total_tokens":2261,"prompt_tokens":974,"completion_tokens":1287,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":590,"completion_tokens_details":{"reasoning_tokens":1206}},"tokens_in":590,"tokens_out":1287,"duration_ms":9039,"temperature":1.0,"reasoning_tokens":1206,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T16:21:30.326539+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A full RCCSDT calculation in the 4ζ basis for the Yb+ 5d 2D3/2 state (or an equivalent large-basis treatment of triples) would settle whether the additivity assumption holds; if the computed Δαd shifts by more than the quoted 4 a.u., the recommended value is biased. Alternatively, a more precise experimental measurement of Δαd for the Yb+ clock transition with uncertainty below 2 a.u. would directly test the central value.","supporting_citations":[{"cited_title":"Lange, N","cited_arxiv_id":null,"evidence_quote":"Reports the 2020 measured quadrupole moment of 1.95(1) a.u. used as the primary comparison for Θ in the full text."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Previous RCCSD calculation of ground and excited state polarizabilities that the paper's larger-basis RCCSDT results supersede."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier RCC calculation of the quadrupole moment that overestimates Θ relative to the recent experiment."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Fock-space coupled-cluster calculation of Θ used as a benchmark for the 5d 2D3/2 and 5d 2D5/2 states."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Configuration-interaction plus core-polarization calculation providing a recent alternative polarizability value for comparison."},{"cited_title":"Schneider, E","cited_arxiv_id":null,"evidence_quote":"Earlier experimental values for both the differential polarizability and quadrupole moment that show the previous discrepancy with theory."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Previous RCCSD calculation in Ca+ whose values the paper reproduces only after adding triple-excitation corrections."}],"review_version":1}