{"id":"c6cc5f33-c993-4b69-b31d-95bd14bcbef7","arxiv_id":"2502.10312","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Relativistic coupled-cluster linear-response calculations produce recommended static dipole scalar and tensor polarizabilities for six low-lying states of Li, Na, and K.","lead":"This paper computes electric dipole polarizabilities for six low-lying states of lithium, sodium, and potassium using several relativistic many-body methods, and recommends coupled-cluster values. These numbers help reduce uncertainties in atomic clocks and cold-atom experiments, where stray electric fields shift the energy levels.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Uncertainty budget in Table V is internally inconsistent: Li 3D values disagree with Ref. [30] by roughly 5 sigma, and Na 3D 'agreement' with Ref. [48] fails by 4–8 sigma, so the accuracy claim is not supported.","rationale":"The paper's central claim is that Table V provides accurate ab initio polarizabilities with realistic uncertainties. The method (relativistic linear-response RCCSD) is well established, and the ground-state and low-lying P-state results match experiments. The weakest link is the uncertainty budget and the agreement claims in Table V. The Li 3D3/2,5/2 states show deviations of 78–82 a.u. from the high-precision calculations of Ref. [30], roughly 4.5 sigma with the quoted combined errors. The paper attributes this to deficiencies of sum-over-states, but no evidence is given that those deficiencies are this large. On the contrary, for Na 3D states the text claims agreement with Ref. [48] while the scalar polarizabilities differ by about 3.6 sigma and 5.6 sigma respectively, an internal inconsistency. Since the quantitative error bars are central to the 'accurate values' claim, this is a load-bearing concern. It does not invalidate the entire paper but requires the authors to provide triple-excitation contributions, basis-convergence checks, and a quantitative reconciliation of the D-state discrepancies. The reader's CONDITIONAL verdict remains appropriate, so no change in verdict is recommended.","tokens_in":16821,"tokens_out":7119,"duration_ms":69760,"concrete_test":"Recompute the Li 3D3/2 scalar polarizability with the same RCCSD approach but using a systematically larger GTO basis (e.g., roughly 30% more functions per symmetry) and including explicit perturbative triples (RCCSD(T)); report the shift from the Table V value of -15003. If the combined shift exceeds the quoted 16 a.u. uncertainty, the error budget is incomplete. If it remains small, the discrepancy with Ref. [30] must be demonstrated to arise from sum-over-states errors, which the paper should then justify quantitatively.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that Table V reports accurate ab initio polarizabilities with uncertainties derived from leading triple excitations. This is not supported for several D states. For Li 3D3/2, the recommended value is -15003(16) a.u., while the high-precision sum-over-states result of Ref. [30] is -14925(8) a.u.; the 78 a.u. difference is about 4.4 sigma under the quoted error bars. For Li 3D5/2 the difference is 82 a.u., about 4.7 sigma. The paper attributes the discrepancy to limitations of sum-over-states, but provides no calculation demonstrating that those limitations produce a shift of this size; it could equally be basis-set incompleteness or omitted triples in RCCSD. Additionally, the text states that Na 3D3/2,5/2 'are in agreement with Ref. [48]', yet the differences are 14.6 a.u. (about 3.6 sigma) and 16.9 a.u. (about 5.6 sigma) respectively. The uncertainty estimate is therefore internally inconsistent with the paper's own comparisons. Because the error bars in Table V are the basis for the 'accurate values' claim, this is load-bearing: it affects the central conclusion, not just presentation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports ab initio calculations of the static electric dipole scalar and tensor polarizabilities of six low-lying states of Li, Na, and K, using Dirac-Hartree-Fock (DHF), third-order many-body perturbation theory (MBPT(3)), random-phase approximation (RPA), and relativistic coupled-cluster singles and doubles (RCCSD). The RCCSD values are recommended and assigned uncertainties estimated from selected perturbative triple-excitation diagrams, and the final numbers are compared with earlier calculations and experiments. The paper also gives a term-by-term decomposition of core, core-valence, and valence contributions for each method.","tokens_in":17111,"tokens_out":4653,"duration_ms":44102,"significance":"If the quoted uncertainties were reliable, the recommended RCCSD values would be a useful systematic dataset for these three alkali-metal atoms: the calculations are genuinely ab initio, no parameters are fitted to polarizability data, and the comparison across DHF, MBPT(3), RPA, and RCCSD gives a clear picture of how correlation effects enter. The extensive tabulations of core, core-valence, and valence contributions, and of individual coupled-cluster terms, are a useful asset. However, the central accuracy claim is not supported for several D states because the quoted uncertainties are inconsistent with the paper's own comparisons to high-precision references. The uncertainty budget is incomplete (no basis-set convergence study and no estimate of several neglected contributions), so the reported 'accurate ab initio values' go beyond what the evidence establishes.","major_comments":[{"comment":"The recommended Li 3D values are not consistent with the uncertainties quoted. For 3D3/2 the recommended value -15003(16) differs from Ref. [30] (-14925(8)) by 78 a.u., which is about 4.4σ when the two quoted errors are combined; for 3D5/2 the difference is 82 a.u., about 4.7σ. The text dismisses this as a limitation of the sum-over-states approach, but no calculation is shown that demonstrates that sum-over-states truncations actually produce a shift of this size. The discrepancy could equally be due to basis-set incompleteness or omitted triple excitations in the present RCCSD calculation. Because the error bars in Table V are the basis of the 'accurate values' claim, this issue is load-bearing and needs to be addressed quantitatively.","section":"§IV, Table V (Li 3D states)"},{"comment":"The statement that the Na 3D3/2,5/2 scalar polarizabilities 'are in agreement with Ref. [48]' is contradicted by the numbers in the same table. The scalar differences are 14.6 a.u. (about 3.6σ) for 3D3/2 and 16.9 a.u. (about 5.6σ) for 3D5/2, using the quoted uncertainties. For the tensor polarizability of 3D5/2 the difference is 10.3 a.u. against a quoted error of 1.3 a.u., about 7.9σ. If the quoted uncertainties are standard deviations, these are not agreements. The comparison claims and the uncertainty estimates cannot both stand as written; either the uncertainties must be enlarged to cover the spread of benchmark results, or the text must be revised to report the discrepancies honestly.","section":"§IV, Table V (Na 3D states)"},{"comment":"The uncertainty budget is incomplete in a way that directly affects the central claim. The text says the uncertainties are 'derived from the leading order triple excitations' and shows only selected diagrams in Fig. 2, but there is no basis-set convergence study in the manuscript. Since polarizabilities of diffuse excited states are sensitive to the Gaussian basis, tests with systematically varied numbers of GTOs and/or GTO parameters are needed. In addition, contributions from omitted triple-excitation diagrams, higher excitations, and Breit/QED effects are not estimated or added to the quoted errors. Without such a convergence study, the reported error bars cannot be taken as reflecting the actual accuracy of the RCCSD method.","section":"§IV (uncertainty estimation) and Table V"}],"minor_comments":[{"comment":"'Block equation' should be 'Bloch equation'; the same misspelling appears in the sentences introducing Eqs. (12).","section":"§III, Eq. (12) and surrounding text"},{"comment":"The publisher location is misspelled: 'Signapore' should be 'Singapore'.","section":"Reference [1]"},{"comment":"Several entries in the K rows contain malformed spacing, such as '0 .02' and '0 .03'; please fix the formatting.","section":"Table III"},{"comment":"The abstract states that core-valence contributions are negligibly small in all methods and atoms, but §IV says that for the K ground state the core-valence contribution 'cannot be disregarded for accurate calculations.' These statements should be reconciled.","section":"Abstract and §IV (K core-valence contribution)"},{"comment":"In the Li 'Others' row, the tensor polarizability entry for 2P3/2 appears as '-3.54 3.54', which looks like a formatting or alignment error; please verify the entry.","section":"Table IV"},{"comment":"The caption reads 'Ratios of scalar and tensor polarizability values from different many-body methods and their DHF values'; this should be '... to their DHF values'.","section":"Fig. 1 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within the journal's scope and the underlying coupled-cluster machinery is standard. The main weakness is uncertainty quantification: the manuscript currently makes 'accurate ab initio values' claims that its own Table V comparisons contradict. I believe this is fixable in a revision, provided the authors add a basis-set convergence study, expand the uncertainty budget, and revise the 'agreement' statements for the D states. I would not recommend rejection because the formalism and the extensive term-by-term analysis are useful; the issue is that the conclusions currently outrun the evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is a systematic RCCSD linear-response calculation of static dipole polarizabilities for the ground and five excited states of Li, Na, and K, with DHF, MBPT(3), and RPA values alongside for comparison. That comparison, plus the breakdown into core, core-valence, and valence contributions, is genuinely useful and well presented. The ground-state and P-state results are solid and agree with experiment.\n\nThe soft spot is the uncertainty budget, and it is load-bearing. The final error bars come only from perturbative triple-excitation shifts; there is no basis-set convergence study and no estimate of higher-order or Breit/QED effects. More tellingly, Table V contradicts its own 'agreement' claims. For Li 3D3/2 and 3D5/2, the recommended values sit roughly 4.5 sigma away from the high-precision sum-over-states results of Ref. [30]. For Na 3D, the text says the values agree with Ref. [48], but the differences are 3.6 and 5.6 times the quoted error bars. The paper attributes the Li discrepancy to limitations of the sum-over-states approach, but nothing in the paper shows that those limitations actually produce a shift of that size; it could just as easily be basis-set incompleteness or omitted triples in the RCCSD calculation. Either way, the quoted uncertainties understate the real error for the D states, and the 'accurate ab initio values' claim rests on those error bars.\n\nWhat is new: the multi-method comparison and the D-state scalar and tensor values are not in the earlier literature, so the table is a useful resource even if some entries need revision. The method itself follows the group's earlier work, so the novelty is modest.\n\nOverall, the paper deserves a serious referee. I would not desk-reject it. But I would send it back for major revision: the authors need to add a basis-set convergence study, expand the uncertainty estimate, and either reconcile the D-state discrepancies or soften the accuracy claims. The agreement statements for Na 3D also need to be corrected.\n\nFor a reading group, it's a decent illustration of how correlation contributions behave across alkali atoms, but the unresolved discrepancy makes it a 'maybe' rather than a 'yes'.\n\nRecommendation: conditional acceptance after major revision.","headline":"A useful systematic RCCSD polarizability table for three alkalis, but the quoted error bars don't cover the D-state discrepancies the paper itself shows, so the 'accurate' claim is overstated.","tokens_in":17624,"tokens_out":2673,"would_cite":false,"duration_ms":25798,"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 theory supplies ab initio scalar and tensor polarizabilities for six low-lying states each of Li, Na, and K.","keywords":["electric dipole polarizability","relativistic coupled-cluster theory","linear response","alkali atoms","lithium","sodium","potassium","scalar and tensor polarizabilities"],"falsifier":"Take the same RCCSD calculation, add full iterative triple excitations in a systematically enlarged Gaussian basis, and check whether the scalar polarizabilities move by more than the quoted uncertainties; a shift larger than the error bars would falsify the uncertainty prescription. A direct Stark-shift measurement of the Li 3D states, where this work gives about $-15003(16)$ a.u. for $3D_{3/2}$ versus $-14925(8)$ a.u. in an earlier high-precision calculation, would also decide whether the one-percent discrepancy is real.","tokens_in":16625,"feed_emoji":"⚛️","tokens_out":8073,"duration_ms":71979,"temperature":0.7,"pith_summary":"This paper sets out to show that a purely ab initio relativistic coupled-cluster calculation, using no experimental energy input, can determine the static electric dipole polarizabilities of lithium, sodium, and potassium for the ground state and five low-lying excited states. The claim matters because these quantities control Stark shifts and blackbody-radiation shifts in atomic clocks and cold-atom experiments, and most earlier excited-state values relied on semi-empirical sum-over-states input with core and continuum contributions added separately. Across all four methods used, the authors find that pair-correlation effects, not core polarization, dominate the correlation corrections, and they quote recommended singles-and-doubles coupled-cluster values with uncertainties estimated from the leading triple-excitation contributions.","feed_headline":"Coupled-cluster theory pins down alkali-atom polarizabilities","feed_subtitle":"Pair correlations, not core polarization, set the scalar and tensor response in Li, Na, and K.","key_machinery":"The central object is the linear-response relativistic coupled-cluster singles-and-doubles (RCCSD) formula, which expresses the polarizability as $$\\$\\alpha$ = 2\\,\\frac{\\langle\\Phi_v|\\{1+$S_v^{{\\dagger}}$\\}\\,\\bar{\\tilde D}\\,\\{$T^{{(1)}}$(1+S_v)+$S_v^{{(1)}}$\\}|\\Phi_v\\rangle}{\\langle\\Phi_v|\\{$S_v^{{\\dagger}}$+1\\}\\,\\bar N\\,\\{1+S_v\\}|\\Phi_v\\rangle},$$ with $\\bar{\\tilde D}=e^{T^{\\dagger}}\\tilde D e^{T}$ the dressed dipole operator. Here $\\tilde D$ encodes the scalar and tensor angular factors from the reduced-matrix-element expression, so one calculation yields both $\\alpha^S_d$ and $\\alpha^T_d$. The machinery solves the first-order perturbation equations for the core ($T^{(1)}$) and valence ($S_v^{(1)}$) cluster amplitudes against the Dirac-Coulomb Hamiltonian, so that core, valence, and continuum intermediate states enter on an equal footing without a sum-over-states truncation. Triple excitations are not included in the wave operators; a selected set of leading triple diagrams is added perturbatively, and the difference from the RCCSD result is used as the uncertainty estimate.","core_discovery":"On its own terms, the paper's central claim is that the recommended values from the relativistic coupled-cluster singles-and-doubles (RCCSD) method in Table V are accurate ab initio scalar and tensor static dipole polarizabilities for the ground and five low-lying excited states of Li, Na, and K. The comparison across methods shows a consistent pattern: Dirac-Hartree-Fock results sit close to random-phase-approximation results, while third-order many-body perturbation theory results sit close to RCCSD results, which the paper reads as evidence that pair correlations, not core polarization, carry the important many-body physics for these atoms. The reported uncertainties come from comparing RCCSD with a perturbative treatment of the leading triple excitations, and the final values agree with available experiments for the ground and low-lying P states. For the Li 3D states, the paper's values differ from earlier sum-over-states results by about one percent, which it attributes to the more complete treatment of core and continuum intermediate states in the ab initio approach.","pith_inferences":["Editorial inference: the same linear-response RCC machinery should become more necessary for the heavier alkalis Rb, Cs, and Fr, where the paper's own core and core-valence contributions grow with atomic number and will be harder to treat by simplified methods.","Editorial inference: the paper's finding that core-valence contributions are negligible for light alkalis suggests that model potentials tuned to reproduce these polarizabilities can safely omit core-valence coupling for Li and Na, but not for K; this is a testable prediction for effective-potential builders.","Editorial inference: the uncertainty protocol, which uses only selected triple-excitation diagrams, could be stress-tested by computing the same quantities with a different basis set and with the Breit interaction; agreement would confirm the quoted error bars, while disagreement would show the error budget is too narrow."],"forward_implications":["The reported values give clock and cold-atom experiments direct ab initio input for Stark-shift and blackbody-radiation systematic corrections in Li, Na, and K.","For the Li $3D_{3/2,5/2}$ states, the paper predicts polarizabilities about one percent larger in magnitude than earlier sum-over-states values, so a future measurement can discriminate between the two treatments.","For the K $4P$ and $3D$ states, the calculation narrows the uncertainty far below existing experiments, giving specific targets for new measurements.","The method-level pattern (DHF close to RPA, MBPT(3) close to RCCSD) implies that any accurate alkali polarizability calculation must include pair correlations beyond the random-phase approximation; core polarization alone is not enough."],"supporting_citations":[{"why":"Supplies the linear-response RCC formulation for electric dipole polarizabilities that this work implements and extends.","marker":"[17]"},{"why":"Provides the Dalgarno-Lewis approach for computing polarizabilities without explicit sum over intermediate states, the conceptual basis of the linear-response treatment.","marker":"[34]"},{"why":"Gives prior relativistic coupled-cluster polarizability values for Li that serve as a central comparison for the new ground and excited-state results.","marker":"[14]"},{"why":"Provides high-precision sum-over-states scalar and tensor polarizabilities for Li and Na against which the RCCSD values are compared.","marker":"[30]"},{"why":"Supplies the relativistic sum-over-states K polarizability values used as the main comparison for the K excited states.","marker":"[32]"},{"why":"Gives semi-empirical DHF-plus-core-polarization polarizabilities for Li, Na, and K that the paper compares with its ab initio numbers.","marker":"[48]"},{"why":"Provides recent many-body K polarizability values and error bars that bracket the paper's recommended K results.","marker":"[58]"},{"why":"Provides the classic beam-deflection measurements of Li, Na, and K ground-state polarizabilities used to validate the recommended values.","marker":"[26]"},{"why":"Supplies the high-precision Li ground-state measurement used to benchmark the Li $2S_{1/2}$ ab initio value.","marker":"[29]"}],"fun_headline_variants":["Pair effects, not core, fix alkali polarizabilities","Coupled-cluster theory yields Li, Na, K polarizabilities","Alkali polarizabilities: pair correlations matter most","Ab initio alkali polarizabilities via RCCSD method"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quoted error bars assume that the neglected electron-correlation terms, especially terms that move three electrons at once, plus basis-set incompleteness and small relativistic or QED effects, are all smaller than the estimated triple-excitation shift; if any of those neglected errors is comparable in size, the reported RCCSD values overstate their accuracy.","fun_headline_variants_meta":{"raw":{"variants":["Pair effects, not core, fix alkali polarizabilities","Coupled-cluster theory yields Li, Na, K polarizabilities","Alkali polarizabilities: pair correlations matter most","Ab initio alkali polarizabilities via RCCSD method"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000244,"raw_usage":{"total_tokens":1552,"prompt_tokens":985,"completion_tokens":567,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":601,"completion_tokens_details":{"reasoning_tokens":496}},"tokens_in":601,"tokens_out":567,"duration_ms":5764,"temperature":1.0,"reasoning_tokens":496,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T18:34:48.811295+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the same RCCSD calculation, add full iterative triple excitations in a systematically enlarged Gaussian basis, and check whether the scalar polarizabilities move by more than the quoted uncertainties; a shift larger than the error bars would falsify the uncertainty prescription. A direct Stark-shift measurement of the Li 3D states, where this work gives about $-15003(16)$ a.u. for $3D_{3/2}$ versus $-14925(8)$ a.u. in an earlier high-precision calculation, would also decide whether the one-percent discrepancy is real.","supporting_citations":[{"cited_title":"Monin and G","cited_arxiv_id":null,"evidence_quote":"Supplies the linear-response RCC formulation for electric dipole polarizabilities that this work implements and extends."},{"cited_title":"Marrus and J","cited_arxiv_id":null,"evidence_quote":"Provides the Dalgarno-Lewis approach for computing polarizabilities without explicit sum over intermediate states, the conceptual basis of the linear-response treatment."},{"cited_title":"[30] applied a sum- over-states method using a mix of experimental and the- oretical energy values","cited_arxiv_id":null,"evidence_quote":"Gives prior relativistic coupled-cluster polarizability values for Li that serve as a central comparison for the new ground and excited-state results."},{"cited_title":"Miffre, M","cited_arxiv_id":null,"evidence_quote":"Provides high-precision sum-over-states scalar and tensor polarizabilities for Li and Na against which the RCCSD values are compared."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the relativistic sum-over-states K polarizability values used as the main comparison for the K excited states."},{"cited_title":"Mitroy, Phys","cited_arxiv_id":null,"evidence_quote":"Gives semi-empirical DHF-plus-core-polarization polarizabilities for Li, Na, and K that the paper compares with its ab initio numbers."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides recent many-body K polarizability values and error bars that bracket the paper's recommended K results."},{"cited_title":"Joachim, J","cited_arxiv_id":null,"evidence_quote":"Provides the classic beam-deflection measurements of Li, Na, and K ground-state polarizabilities used to validate the recommended values."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the high-precision Li ground-state measurement used to benchmark the Li $2S_{1/2}$ ab initio value."}],"review_version":1}