{"id":"19b82f02-2e3b-494a-a35e-c46071083788","arxiv_id":"2504.19515","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"For Mg+, MCDHF/RCI and relativistic coupled-cluster calculations agree on energies, lifetimes, and polarizability, while coupled-cluster matches measured hyperfine and isotope-shift constants more closely.","lead":"This paper compares two advanced computational methods for calculating the properties of a magnesium ion, and shows they mostly agree, with coupled-cluster matching measurements slightly better for hyperfine and isotope-shift quantities. It also proposes more precise lifetimes for excited states than current experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The RCC-over-MCDHF preference for IS constants rests on SMS values whose convergence is not established: triples shift the 3s-4s SMS by roughly 53 GHz u and the D-line SMS by about 31 GHz u, yet quoted uncertainties are 5-15 GHz u, and the paper itself says quadruples are non-negligible.","rationale":"The paper has real strengths: the layer-by-layer MCDHF/RCI convergence tables, the independent RCC sequence, consistent lifetime estimates from both methods, and an RCC ground-state hyperfine constant that matches experiment. Those parts support the paper's broader comparative claim. The concentrated weakness is the IS/SMS recommendation. The authors themselves flag quadruple substitutions as non-negligible for SMS constants, yet the final SMS uncertainties in Table XV are smaller than the observed triple-excitation shifts in the same quantity. The only experimental SMS benchmark, the D1/D2 lines, is also the basis for the statement that RCC IS constants agree better with measurements. Because the omitted corrections are acknowledged to be potentially significant and are not included in the quoted uncertainties, the Table XV recommended values and the RCC-over-MCDHF preference should be treated as provisional. This is an internal convergence issue, not a disagreement with consensus, and it is directly testable with a higher-order calculation. The reader's CONDITIONAL verdict already captures this appropriately; the concern sharpens the condition but does not change the verdict.","tokens_in":28019,"tokens_out":7322,"duration_ms":78253,"concrete_test":"Recompute the differential SMS constants of Table XV at the RCCSDT level with h and i orbitals included in the triple-excitation space (or with a perturbative quadruples correction), using the same finite-field protocol. If the 3s-4s SMS value changes by more than about 10 GHz u, or if either D-line SMS value shifts by more than 15 GHz u relative to the quoted 356(15)/362(15) GHz u, then the neglected higher-order terms are not small and the RCC-over-MCDHF IS conclusion is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that RCC IS constants agree with measurements better than the MCDHF/RCI values is carried by the SMS differential constants in Table XV, but these are the least converged quantities in the paper. In the RCC sequence (Table XII), adding triples changes the 3s-4s SMS constant from +55.7 to +3.1 GHz u and changes the D-line SMS constants by about -31 GHz u, while the uncertainties quoted in Table XV are only 5 GHz u and 15 GHz u, respectively. Section IV.D explicitly states that quadruple substitutions are 'expected to be non-negligible' for SMS constants, and the Summary repeats that quadruple excitations may further improve the RCC results. The RCCSDT calculations are also truncated to g symmetry, with the h/i corrections to Table XV computed only at the RCCSD level. Thus the load-bearing assumption that neglected higher-order correlation is smaller than the method differences is not established for exactly the quantity used to benchmark RCC against MCDHF. If quadruple or high-l triple effects shift the SMS by an amount comparable to the observed triple-excitation shift, the RCC D-line values (356 and 362 GHz u) would move away from the experimental values (369.3 and 367.7 GHz u), and the sign of the recommended 3s-4s SMS could flip. The MCDHF side is similarly exposed: triple substitutions are restricted to the {12s,10p,5d} active space, and the paper notes that further core-core correlation layers are needed for IS convergence.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a systematic comparison of two many-body methods, MCDHF/RCI and relativistic coupled-cluster (RCCSD/RCCSDT), applied to low-lying states of Mg+. For each method it computes excitation energies, E1 matrix elements, lifetimes, ground-state electric dipole polarizability, hyperfine A constants, and isotope-shift constants (field shift, normal mass shift, and specific mass shift), with layer-by-layer or method-level convergence analysis. The MCDHF/RCI results are taken from layer 7 with triple-excitation corrections, while the RCC results include basis, Breit, and QED corrections. The central claim is that the two methods agree well for most properties, and that the RCC values for hyperfine constants and isotope-shift constants agree with experiment better than the MCDHF/RCI values, leading to recommended lifetimes and differential IS constants.","tokens_in":28352,"tokens_out":6053,"duration_ms":65838,"significance":"If the results hold, the paper provides useful benchmarks for Mg+ — in particular excited-state lifetimes that are more precise than the currently available experimental values — and general guidance on when MCDHF/RCI and RCC results can be trusted. The paper has clear strengths: systematic layer-by-layer expansions, use of three independent approaches to isotope-shift constants (finite-field, expectation-value, and analytic response), and no property fitted to the benchmark data. The cross-method consistency for E1 matrix elements, lifetimes, polarizability, and field-shift constants is convincing. The central recommendation, however, is only as strong as the convergence of the specific mass shift (SMS) constants, and here the evidence presented is incomplete.","major_comments":[{"comment":"The load-bearing claim that RCC IS constants agree better with experiment rests on differential SMS constants that are the least converged quantities in the paper. In Table XII, the RCCSD-to-RCCSDT shifts in the finite-field SMS constants are large: the ground-state 3s value changes from 111.34 to 46.03 GHz u, the 3p 2P°_1/2 value changes from -276.20 to -310.19 GHz u, and the 3p 2P°_3/2 value changes from -280.93 to -314.74 GHz u. The final differential values in Table XV (356(15) and 362(15) GHz u for the D lines) inherit these triple-excitation shifts, yet the quoted uncertainties are smaller than the triples shift, and Sec. IV.D explicitly states that quadruple substitutions are expected to be non-negligible for SMS constants. In addition, the RCCSDT calculation is truncated at g symmetry, and the h/i corrections in Table XIV are computed only at the RCCSD level. If neglected quadruple or high-l triple effects shift the SMS constants by an amount comparable to the observed triple-excitation shifts, the RCC D-line values would move away from the experimental values, and the sign of the 3s-4s SMS constant could flip. The claimed preference for RCC over MCDHF for the D-line SMS values is therefore not established without a quantitative estimate of these neglected effects.","section":"Sec. IV.D, Tables XII and XV"},{"comment":"The MCDHF/RCI side of the comparison is similarly affected by truncation. Triple substitutions are restricted to the {12s,10p,5d} active space, and the SD-to-SDT changes in Table XIII are substantial: K_SMS for 3s-4s changes from -75 to -54 GHz u, and for the D lines from 304 to 321 GHz u. The text notes that further core-core correlation layers are needed for isotope-shift convergence. Since the central value of the RCC 3s-4s SMS constant has the opposite sign to the MCDHF/RCI value and is reported as 3(5) GHz u, the current calculations do not establish which method is more reliable for the SMS constants. At minimum, the comparison should be accompanied by a sensitivity estimate based on the observed SD-to-SDT shifts and the known truncation of the triple-excitation space.","section":"Sec. IV.D and Table XIII"},{"comment":"The quoted uncertainties in Table XV are not derived or explained in the text. For K_SMS, uncertainties of 5 GHz u and 10-15 GHz u are assigned to transitions for which the triples contribution alone is tens of GHz u and quadruples are stated to be non-negligible. The paper should state explicitly how the uncertainties were obtained, and demonstrate that they account for truncation of the RCCSDT excitation space to g symmetry, treatment of h/i corrections only at RCCSD, and neglect of quadruple substitutions. Without this, the uncertainties give a false impression of the reliability of the recommended SMS constants.","section":"Table XV"}],"minor_comments":[{"comment":"The table heading states that E1 transition rates are given \"in s\", but rates are in s^-1; please correct the units.","section":"Table IV"},{"comment":"The text contains an unresolved citation \"[ ? ]\" after presenting the relativistic NMS and SMS operators; this should be replaced with the proper reference, presumably Ref. [38].","section":"Sec. II, after Eq. (15)"},{"comment":"The definition of the mean level deviation and the energy shift ES is unclear: it should be stated explicitly whether ES is added or subtracted, and the quoted values (ES = -183 cm^-1 and ES = -149 cm^-1) should be interpreted in terms of the ground-state binding-energy imbalance.","section":"Eq. (42) and surrounding text"},{"comment":"The column labeled \"Scaling [36]\" for the NMS constants is not defined in the text; please state how the scaled values are obtained from the NIST ASD ionization potentials.","section":"Table XII"},{"comment":"There is a typo in \"Babaushkin\" (should be Babushkin); similar minor typographical errors appear elsewhere, including \"bench-marked\" in the abstract, and would benefit from a careful proofreading pass.","section":"Sec. IV.A"},{"comment":"References [3] and [20] are identical (Nataraj et al.); please consolidate or cite different relevant works if duplication was unintended.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal's scope and the comparative data set will be useful to the atomic-physics community, but the abstract's claim that RCC IS constants agree with measurements better than MCDHF/RCI is stronger than the convergence evidence for the SMS constants. I would advise the editor that the recommended SMS constants and the RCC-vs-MCDHF ranking should either be supported by additional higher-order calculations (for example, RCCSDT with h/i orbitals and a quadruple-excitation estimate) or be presented with substantially more conservative uncertainties. No concerns about attribution or novelty."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Things to know: this is a competent, honest benchmark paper, not a breakthrough. It runs MCDHF/RCI and RCC on Mg+ for energies, E1 amplitudes, hyperfine constants, isotope-shift constants, and derived lifetimes and polarizability, with no property fitted to experiment. The layer-by-layer convergence study is the real value: it shows how different operators converge differently in the two methods.\n\nWhat is good: the recommended lifetimes from the two methods agree and are more precise than the available measurements; the RCC ground-state hyperfine constant lands within 0.4 MHz of experiment, with QED corrections helping; the FS constants agree between methods and with the one experimental value; the RCC NMS values track the experimental-energy scaling. The paper is also frank about its own gaps, which is rare.\n\nThe soft spot is the headline IS claim. The abstract's statement that RCC IS constants agree with measurements better is carried by the SMS differentials in Table XV. But the triples contribution to those SMS differentials is roughly 52 GHz u for the 3s-4s line and 31 GHz u for the D-lines, while the quoted uncertainties are 5 and 15 GHz u. The paper itself states that quadruple substitutions are expected to be non-negligible for SMS. So the RCC-over-MCDHF preference rests on exactly the quantity that is least converged. The authors flag this, but it means the recommended SMS constants and their error bars are provisional, not final. A shift comparable to the triples effect would move the D-line values away from experiment and could flip the sign of the 3s-4s SMS.\n\nOther, smaller issues: the MCDHF values come without explicit uncertainties; no code or input data are provided (minor, since the codes are established); there is a dangling reference placeholder in the NMS/SMS operator definition and a missing citation for the earlier MCDHF calculation that produced -595.262 MHz for the ground-state A; and the claim that RCC hyperfine agrees better with experiment relies on a single precisely measured state.\n\nNone of this is disqualifying. The comparative benchmark itself is useful and the recommended lifetimes and FS constants are likely fine. But the IS preference for RCC should be framed as preliminary, with wider error bars, until SMS convergence is demonstrated.\n\nI would send this to peer review: it deserves a serious referee. The revision should address the SMS convergence and uncertainty question and clean up the references, but the core work is legitimate and worth publishing. Reading group: maybe; cite: yes.","headline":"Solid, honest MCDHF-vs-RCC benchmark on Mg+ that is useful for lifetimes and FS constants, but its headline claim about RCC isotope-shift agreement is underpinned by SMS values whose convergence is not yet established.","tokens_in":28851,"tokens_out":8617,"would_cite":true,"duration_ms":79683,"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 paper benchmarks Mg+ properties with two many-body methods and finds that relativistic coupled-cluster matches measured hyperfine and isotope-shift constants better than MCDHF.","keywords":["multiconfiguration Dirac-Hartree-Fock","relativistic coupled-cluster","magnesium ion Mg+","hyperfine structure constants","isotope shift constants","E1 transition matrix elements","atomic lifetimes","dipole polarizability"],"falsifier":"Measure the 3s→4s isotope shift in 24Mg+–26Mg+ (or 25Mg+–24Mg+) with enough precision to fix the sign and magnitude of the specific mass-shift constant; the two methods predict opposite signs (-54 GHz amu from MCDHF/RCI versus +3(5) GHz amu from RCC), so the measurement would immediately show which treatment of two-body correlation is right. A second check would be a remeasurement of the ground-state polarizability to see whether it falls near the calculated 34.9–35.2 a.u. or near the old experimental 33.8 a.u.","tokens_in":27827,"feed_emoji":"⚛️","tokens_out":6592,"duration_ms":62701,"temperature":0.7,"pith_summary":"The paper sets out to see how two standard atomic many-body methods, multiconfiguration Dirac-Hartree-Fock with relativistic configuration interaction (MCDHF/RCI) and relativistic coupled-cluster (RCC), behave as correlation is added layer by layer for singly ionized magnesium. It benchmarks excitation energies, E1 transition matrix elements, magnetic-dipole hyperfine constants, isotope-shift constants, and derived lifetimes and polarizability against measured values. The central finding is that the two methods agree with each other for most properties, but for hyperfine constants and isotope shifts the RCC results track the measurements more closely. The paper therefore recommends the RCC-based lifetimes and differential isotope-shift constants as the more reliable inputs for future experiments and for benchmarking other atomic systems.","feed_headline":"Coupled-cluster matches measured Mg+ hyperfine and isotope shifts","feed_subtitle":"Head-to-head benchmark of two atomic many-body methods says where they agree—and which values to trust for experiments.","key_machinery":"The argument is carried by two complementary wave-function machines and a set of operator-specific evaluation protocols. The MCDHF/RCI side builds correlation in seven orbital layers with singles and doubles substitutions plus a restricted set of triple substitutions, and estimates transition uncertainties by scanning the full gauge parameter rather than just comparing two gauges. The RCC side uses a Fock-space coupled-cluster ansatz, $e^T\\{1+S_v\\}$, with singles-doubles-triples amplitudes, finite-field, expectation-value, and analytical-response routes to isotope-shift constants, plus separate basis, Breit, and QED corrections. The isotope-shift operators—field shift, normal mass shift, and specific mass shift—have different radial sensitivities, and the paper tracks which approximation level stabilizes each, which is the mechanism that lets it attribute the RCC advantage to a better treatment of two-body correlation.","core_discovery":"On the paper's own terms, the discovery is a calibration: in Mg+, RCCSDT with basis, Breit, and QED corrections reproduces the measured ground-state hyperfine constant, -596.6(8) MHz versus -596.2542487(42) MHz, and gives differential isotope-shift constants for the D1 and D2 lines whose specific mass shifts (356(15) and 362(15) GHz amu) sit close to experiment, while the MCDHF/RCI values (321 and 321 GHz amu) deviate more. The E1 amplitudes from the two methods agree to the second decimal place, leading to essentially identical lifetimes and a ground-state polarizability of 35.16 (MCDHF/RCI) versus 34.92(7) (RCC), both above the single available measurement. Where the methods differ, the difference is systematic: correlation trends for one-body operators (energies, E1) converge smoothly, while two-body isotope-shift operators, especially the specific mass shift, remain sensitive to the treatment of triples and beyond.","pith_inferences":["Editorial inference: the same head-to-head protocol applied to heavier alkaline-earth ions such as Ca+, Sr+, or Ba+ would probably show larger MCDHF/RCI disagreements, because relativistic and core-correlation effects grow with Z; the Mg+ agreement should not be read as a general license.","Editorial inference: the opposite signs predicted for the 3s→4s specific mass shift make that transition a sharper experimental test than the D lines, where both methods already sit closer to experiment.","Editorial inference: since the field-shift constants are robust but the specific mass shift is not, future method work on isotope shifts should concentrate on two-body correlation treatments (full connected triples, or explicit quadruples) rather than on enlarging one-particle basis sets.","Editorial inference: the gauge-scan approach implies that for other ions, relying on length-velocity gauge agreement to certify E1 accuracy can be misleading; the full gauge-parameter dependence should be checked before quoting sub-percent uncertainties."],"forward_implications":["The recommended lifetimes for the 3p, 4s, 3d, and 4p states of Mg+ are more precise than any single available measurement, so they can serve as reference values for beam or trap lifetime experiments.","Because both methods put the ground-state polarizability near 35 a.u., the existing experimental value near 33.8 a.u. is called into question and should be re-examined.","For isotope-shift work, the field-shift constants agree across methods and with experiment, so nuclear-charge-radius differences extracted from Mg+ data are on solid ground; the specific mass shift is the part that needs the RCC treatment.","The layer-by-layer comparison provides a template for judging where MCDHF/RCI and RCC can be trusted in heavier ions, where no measurements exist for quantities like parity-violation or time-reversal-violation enhancement factors."],"supporting_citations":[{"why":"Supplies the MCDHF theory and computational framework used for the layer-by-layer wave functions.","marker":"[12]"},{"why":"Documents the MCDHF/RCI implementation that produces the energies and E1 matrix elements.","marker":"[13]"},{"why":"Establishes the coupled-cluster formalism that the RCC energy and property calculations are built on.","marker":"[18]"},{"why":"Introduces the analytical-response approach used as one of the three routes to isotope-shift constants.","marker":"[32]"},{"why":"Provides the measured excitation energies and the scaling values used to judge the NMS constants.","marker":"[36]"},{"why":"Supplies the program used to evaluate isotope-shift parameters in the MCDHF/RCI calculations.","marker":"[37]"},{"why":"High-precision measured ground-state hyperfine constant against which RCC and MCDHF/RCI values are judged.","marker":"[60]"},{"why":"Independent high-precision measurement of the ground-state hyperfine constant used for the same comparison.","marker":"[61]"},{"why":"Measured field-shift constant for the D1 line used to validate the recommended FS values.","marker":"[71]"},{"why":"Measured specific-mass-shift constants of the D1/D2 lines used to show RCC agrees better with experiment.","marker":"[72]"}],"fun_headline_variants":["RCC edges out MCDHF on Mg+ hyperfine and isotope shifts","For Mg+ hyperfine and shifts, RCC outperforms MCDHF","Mg+ atomic benchmark: RCC closer to measured hyperfine","Coupled-cluster better for Mg+ isotope shifts than MCDHF","Mg+ hyperfine and isotope constants: RCC ahead of MCDHF"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The comparison assumes that correlation effects left out of both methods—especially triple excitations beyond a small orbital set in MCDHF/RCI, and quadruple excitations in RCC—are smaller than the method-to-method differences under discussion; the paper itself notes that quadruple substitutions are expected to be non-negligible for specific mass shifts.","fun_headline_variants_meta":{"raw":{"variants":["RCC edges out MCDHF on Mg+ hyperfine and isotope shifts","For Mg+ hyperfine and shifts, RCC outperforms MCDHF","Mg+ atomic benchmark: RCC closer to measured hyperfine","Coupled-cluster better for Mg+ isotope shifts than MCDHF","Mg+ hyperfine and isotope constants: RCC ahead of MCDHF"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000788,"raw_usage":{"total_tokens":3488,"prompt_tokens":973,"completion_tokens":2515,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":589,"completion_tokens_details":{"reasoning_tokens":2418}},"tokens_in":589,"tokens_out":2515,"duration_ms":20065,"temperature":1.0,"reasoning_tokens":2418,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:50:37.971763+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the 3s→4s isotope shift in 24Mg+–26Mg+ (or 25Mg+–24Mg+) with enough precision to fix the sign and magnitude of the specific mass-shift constant; the two methods predict opposite signs (-54 GHz amu from MCDHF/RCI versus +3(5) GHz amu from RCC), so the measurement would immediately show which treatment of two-body correlation is right. A second check would be a remeasurement of the ground-state polarizability to see whether it falls near the calculated 34.9–35.2 a.u. or near the old experimental 33.8 a.u.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the MCDHF theory and computational framework used for the layer-by-layer wave functions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the coupled-cluster formalism that the RCC energy and property calculations are built on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the analytical-response approach used as one of the three routes to isotope-shift constants."},{"cited_title":"J¨ onsson, G","cited_arxiv_id":null,"evidence_quote":"Provides the measured excitation energies and the scaling values used to judge the NMS constants."},{"cited_title":"Lindgren and J","cited_arxiv_id":null,"evidence_quote":"Supplies the program used to evaluate isotope-shift parameters in the MCDHF/RCI calculations."},{"cited_title":"Lindgren and D","cited_arxiv_id":null,"evidence_quote":"High-precision measured ground-state hyperfine constant against which RCC and MCDHF/RCI values are judged."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Independent high-precision measurement of the ground-state hyperfine constant used for the same comparison."},{"cited_title":"Schiffmann, J","cited_arxiv_id":null,"evidence_quote":"Measured field-shift constant for the D1 line used to validate the recommended FS values."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Measured specific-mass-shift constants of the D1/D2 lines used to show RCC agrees better with experiment."}],"review_version":1}