{"id":"e83c9835-81a4-462d-9b4b-2348b923a8b9","arxiv_id":"2412.04044","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"New MCDHF-CI energy levels and M1 transition rates are reported for I8+ through I12+, including first published data for I10+, I11+, and I12+.","lead":"This paper calculates the energy levels and magnetic transition rates for five highly charged iodine ions using a standard atomic-structure method. It provides the first published data for the three most stripped ions, which had no such values in the literature.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 3% accuracy class for unmeasured I11+ and I12+ levels rests on an extrapolation that the paper's own I9+ comparison shows is unreliable, since convergence-based errors (0.3%) are nearly ten times smaller than the observed theory–experiment deviations (up to 2.6%).","rationale":"The reader's weakest_assumption correctly identifies the extrapolation of the 3% accuracy class from I8+/I9+ to the unmeasured I10+/I12+ ions. My stress-test confirms this is the single most load-bearing concern: the paper's own convergence-based uncertainty (Section 3.2) is inconsistent with its experiment comparison (Section 3.3), showing that the dominant error is systematic rather than captured by AS2→AS3 convergence. Because all I9+ deviations are positive, a systematic missing contribution (e.g., core-valence correlation or higher-order QED) may scale with configuration complexity, making the 3% bound for 4d⁶ and 4d⁵ configurations fragile. The proposed independent calculation or isoelectronic benchmark would settle whether the extrapolation holds. The paper is otherwise a standard, transparent MCDHF-CI study with useful tabulated data; the issue is the over-strong accuracy claim for the new ions, not the underlying calculation. Therefore the conditional verdict remains appropriate, with no adjustment needed.","tokens_in":22493,"tokens_out":4420,"duration_ms":47916,"concrete_test":"Perform independent relativistic many-body calculations for the first 10 low-lying levels of I11+ and I12+ using a different methodology (e.g., CI+MBPT or all-order coupled-cluster, or a second MCDHF implementation with different core-valence correlation treatment) and compare the resulting level energies with Tables 4 and 5. If any level differs by more than 3% from the paper's value, the extrapolated A accuracy class is not supported and the central claim of high accuracy for the unmeasured ions would need to be softened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that these MCDHF-CI results are high-accuracy reference data for I10+ through I12+ depends on the uncertainty estimate in Section 3.3: 'we roughly estimate the uncertainties for energy levels for I9+ through I12+ to be 3% (A accuracy class).' This estimate is not supported by the paper's internal evidence. In Section 3.2, the convergence-based uncertainty δE = |E_AS3 − E_AS2| is reported as up to 0.3% for I9+ level energies. Yet Section 3.3 and Table D show that the same I9+ levels deviate from experiment by 0.1% to 2.6%, with every theoretical value systematically above the measured one (e.g., level 5 at 30638.93 cm⁻¹ vs 29860 cm⁻¹, level 9 at 83632.17 cm⁻¹ vs 81670 cm⁻¹). Thus the convergence metric underestimates the true error by roughly an order of magnitude and does not capture a systematic offset. The 3% class is then obtained by simply taking the maximum observed I8+/I9+ deviation and assuming it bounds the errors for I10+ through I12+, for which no measured spectra exist. This is the load-bearing step: if the systematic error grows with open-shell complexity or with charge state, the claimed A accuracy class fails for exactly the ions the paper advertises as new data. The paper is transparent about the lack of reference values, but the assumption that the I8+/I9+ error envelope transfers to the unmeasured ions is not independently tested.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports GRASP2018-based Multi-Configuration Dirac-Hartree-Fock with Configuration Interaction (MCDHF-CI) calculations of energy levels and M1 transition rates for the [Kr]4d^n (n=5-9) valence configurations of I8+ through I12+, with the explicit goal of providing reference data for I10+ through I12+, for which no measured spectra are available. The authors document active-space convergence (AS0 through AS3), estimate theoretical uncertainties from the AS2-to-AS3 difference, benchmark the I8+ and I9+ level energies against published experimental values, and compare selected M1 rates with FAC calculations. The main advertised result is a set of level-energy tables and transition-rate tables for the five iodine ions, with an A accuracy class (about 3%) claimed for the unmeasured ions I9+-I12+.","tokens_in":22856,"tokens_out":2785,"duration_ms":31290,"significance":"If the accuracy claims hold, the tables would be a useful reference for plasma spectroscopy, astrophysical modeling, and studies of forbidden transitions in open-d-shell highly charged ions, filling a genuine gap in the literature for I10+-I12+. The paper's strengths are that the calculations follow a systematic active-space protocol, the convergence pattern is reported for every ion, the I8+ and I9+ comparisons are genuine external benchmarks, and no parameter is fitted to the target data. The data tables are extensive and clearly organized. The main weakness is that the uncertainty classification for the unmeasured ions rests on an extrapolation that the paper's own benchmark data do not fully support; this weakens the central 'high-accuracy' claim unless it is corrected or recast.","major_comments":[{"comment":"The assignment of a 3% (A accuracy class) uncertainty to I9+-I12+ level energies is not supported by the paper's internal evidence. The convergence-based uncertainty in §3.2 is up to 0.3% for I9+, yet Table D shows theory-experiment deviations for the same ion ranging from 0.1% to 2.6%, with all theoretical values systematically above the experimental ones (e.g., level 5 at 30638.93 cm^-1 vs 29860 cm^-1 and level 9 at 83632.17 cm^-1 vs 81670 cm^-1). Thus the AS2/AS3 convergence metric underestimates the true error by roughly an order of magnitude and does not capture the systematic offset. Taking the maximum observed I8+/I9+ deviation as the uncertainty for I10+, I11+, and I12+ is an extrapolation that needs either an explicit physical justification (e.g., why the error envelope should not grow with open-shell complexity or charge state) or a substantially more conservative accuracy statement. This is load-bearing because the abstract and conclusions advertise high-accuracy data for exactly the ions with no measured spectra.","section":"§3.3 and Table D"},{"comment":"The error budget is incomplete in a way that matters for the stated accuracy class. The paper accounts for active-space convergence (AS2 vs AS3) and, for transition rates, a GRASP-FAC comparison, but it does not quantify the combined contribution of the Breit interaction treatment, the QED model (screened hydrogenic self-energy and Fullerton-Rinker vacuum polarization), or the effect of restricting substitutions to 4s, 4p, and 4d subshells while keeping the n=1-3 core frozen. Because the I9+ benchmark shows a systematic offset much larger than the convergence estimate, one or more of these omitted terms likely dominate the error. I request a transparent discussion of these contributions, a level- or term-dependent error estimate, or a downgrade of the accuracy class for the unmeasured ions to a more defensible statement such as 'estimated to be within a few percent based on the I8+/I9+ envelope, but with an unquantified systematic component.'","section":"§3.2-§3.3"}],"minor_comments":[{"comment":"The column header 'Upper level No. Upper level No.' appears in each transition table; the second column should be labeled 'Lower level No.' or 'Lower level' for clarity.","section":"Tables 6-10"},{"comment":"The phrase 'Basing on this finding' should read 'Based on this finding'.","section":"§3.3"},{"comment":"The sentence 'Recent study suggested that HCIs with d(4−6) configurations might host suitable forbidden transitions' has a number and tense error ('Recent study suggested'); it should be 'A recent study suggested' or 'Recent studies have suggested'.","section":"§1"},{"comment":"The definition of the transition-rate uncertainty (selecting transitions with intensity above 10% of the strongest and extending that uncertainty to all transitions) should be justified for weak transitions, because the fractional uncertainty of a weak M1 rate can be much larger than that of the strong transitions used for calibration.","section":"§3.2"},{"comment":"The notation |∆| = |AGRASP − AFAC|/AGRASP should be defined in the table caption, and the values are reported without explicit signs; specifying the sign would help identify systematic biases in the FAC comparison.","section":"Table B"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid data-production effort with a clear protocol, and the I8+/I9+ benchmarks are honest. The central issue is that the 3% A accuracy class for the unmeasured I10+-I12+ levels is an extrapolation that is not supported by the paper's own benchmark data; this needs either a quantitative justification or a more cautious accuracy statement. As written, the claim of 'high-accuracy' reference data for those ions is not yet established. If the authors recalibrate the uncertainty statement and add an explicit error budget, the manuscript would be suitable for publication as a data reference."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a straightforward, well-executed MCDHF-CI calculation that fills a real gap — there were no published level energies for I10+, I11+, I12+. The tables will be useful to anyone modeling iodine plasmas or planning spectroscopy on these ions. The benchmark against measured I8+/I9+ energies, with agreement within 2.6%, gives me reasonable confidence in the numbers.\n\nWhat the paper does well: it uses the standard GRASP2018 protocol with a sensible active-space expansion (AS0–AS3) and reports both convergence behavior and explicit comparison with experiment. The transition-rate tables come with a GRASP–FAC cross-check, which is exactly the right kind of independent sanity check for M1 rates. The paper is transparent about what is new (I10+–I12+) and what is benchmarked (I8+, I9+).\n\nThe soft spots are real but not disqualifying. The biggest one is the accuracy-class assignment in Section 3.3. The paper states \"we roughly estimate the uncertainties for energy levels for I9+ through I12+ to be 3% (A accuracy class).\" That 3% is taken as the maximum observed theory–experiment deviation for I8+/I9+ and carried over to the unmeasured ions. The stress-test note is right that the convergence-based δE (0.3–0.5%) is not the error budget — but the paper doesn't use δE for the final claim; it uses the experiment comparison. The real question is whether a 2.6% error in I9+ bounds the error in I11+/I12+, which have more open-shell electrons and no measured spectra. That's an extrapolation, and the paper should say so more explicitly rather than calling it an estimate. I'd ask the authors to soften the language to \"expected accuracy\" and add a sentence acknowledging the possibility of larger errors for the more complex ions.\n\nA second, smaller issue: the transition-rate uncertainty is derived from the high-intensity subset and then extended to all listed rates. For weak M1 rates, cancellations can make the error much larger. The GRASP–FAC comparison partially mitigates this, but it only covers three transitions per ion. Worth a comment in the paper.\n\nAlso, in Table D every theoretical I9+ level is above experiment — a systematic offset. That's not a fatal problem, but it suggests a missing correlation contribution, and the paper doesn't discuss it. A referee should ask about it.\n\nOverall: this is a legitimate reference-data paper. It doesn't introduce new methodology, but it doesn't need to. I'd send it to peer review, likely accept after minor revision. The tables will get cited.","headline":"Solid, useful reference data for iodine ions; the 3% accuracy-class claim for unmeasured ions is a fair extrapolation but needs to be flagged as such.","tokens_in":23329,"tokens_out":3499,"would_cite":false,"duration_ms":33938,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["32.30.-r","31.15.Ar","32.70.Cs"],"model":"deepseek-v4-flash","headline":"High-accuracy MCDHF-CI calculations give reliable level energies and M1 transition rates for the low-lying $[\\mathrm{Kr}]4d^n$ states in I8+ through I12+, and the I10+ through I12+ results fill a gap in the literature.","keywords":["iodine ions","highly charged ions","energy levels","M1 transition rates","MCDHF-CI","atomic data","4d open-shell configurations","forbidden transitions"],"falsifier":"A clean test would be an electron-beam ion trap or beam-foil measurement of a predicted M1 line in I10+, I11+, or I12+: for example, an I10+ M1 transition predicted at 295.20 nm (level 12 to level 2) or the I12+ line predicted at 265.87 nm (level 2 to level 1). A measured wavelength that differs from the predicted value by more than the assigned 3% would show that the uncertainty estimate underestimates the true error, while agreement within the stated class would validate the tables as reference data.","tokens_in":22313,"feed_emoji":"⚛️","tokens_out":9648,"duration_ms":80434,"temperature":0.7,"pith_summary":"The paper reports high-accuracy Multi-Configuration Dirac–Hartree–Fock with Configuration Interaction (MCDHF-CI) calculations of the low-lying energy levels and magnetic-dipole transition rates of the five iodine ions I8+ through I12+, whose valence configurations are $[\\mathrm{Kr}]4d^n$ for $n=5$–$9$. For I10+ through I12+, the authors argue, these tables fill a genuine gap in the published data, since no detailed analysis of the I XI, I XII, and I XIII spectra exists in the literature. The results are argued to be converged at the AS3 active-space stage, with level energies for I9+ through I12+ assigned roughly 3% (A accuracy class) uncertainties and the strongest M1 rates assigned B or C+ classes. The evidence is convergence between successive active spaces, agreement with measured I8+ and I9+ spectra to within 2.6%, and cross-checks of the strongest transitions against an independent relativistic code.","feed_headline":"New energies fill data gap for iodine ions I8+ to I12+","feed_subtitle":"Calculated levels and M1 rates cover I10+ to I12+, which had no published reference data.","key_machinery":"The machinery is the MCDHF-CI expansion of atomic state functions as linear combinations of configuration state functions, with virtual orbitals included through principal quantum number $n=7$ and orbital angular momentum $l=5$. The calculation steps through four active spaces, AS0 through AS3, with AS3 containing up to about $2\\times10^6$ CSFs for the heavier ions; the Breit interaction is added perturbatively and QED self-energy and vacuum-polarization corrections are included. The difference between the AS3 and AS2 energies serves as the estimated level-energy uncertainty, transition-rate uncertainties follow the Kramida procedure, and the independent FAC code provides a cross-check on the strongest M1 rates. This machinery is what supports the claim that the AS3 results are converged and that the I10+ through I12+ tables can stand as reference data.","core_discovery":"The central claim is that a consistent set of MCDHF-CI calculations, carried out with the GRASP 2018 code, gives reliable level energies and M1 transition rates for the $[\\mathrm{Kr}]4d^n$ ($n=5$–$9$) configurations of I8+ through I12+. For the two lowest-charge ions the computed levels reproduce experimentally deduced spectra within 2.6%, while for I10+, I11+, and I12+ the present results are offered as the missing reference data. The authors assign an A accuracy class (roughly 3% uncertainties) to the I9+ through I12+ level energies and B or C+ classes to the tabulated M1 transition rates, based on active-space convergence and on differences between the GRASP 2018 and FAC calculations.","pith_inferences":["Applying the same AS2/AS3 convergence procedure to neighbouring tellurium and xenon ions would likely show similar accuracy, since the open-$d$-shell correlation structure varies smoothly along isoelectronic sequences; this is an inference, not a claim in the paper.","A systematic GRASP-versus-FAC comparison across all tabulated transitions, rather than only the strongest ones, would turn the overall B/C+ rate classification into a per-transition reliability map, which the paper does not provide.","Extending the active space beyond $n=7$ would indicate how much of the residual level-energy uncertainty comes from basis-set truncation rather than from physical effects such as higher-order QED or core-valence correlation."],"forward_implications":["The tables give the first computed energy-level and M1 transition-rate data for I10+, I11+, and I12+ in the $[\\mathrm{Kr}]4d^n$ configurations, filling the gap noted in the NIST Atomic Spectra Database.","The I8+ and I9+ comparisons against measured spectra provide a benchmark for using the same method on neighbouring ions where experimental data are missing.","The listed M1 transitions, especially the strongest lines, give plasma diagnosticians and astrophysical modelers direct input for line identifications in iodine spectra.","The I11+ ($4d^6$, Mo-like) and I12+ ($4d^5$, Nb-like) data connect to the search for forbidden transitions suitable as optical-clock frequency references, a direction the introduction explicitly raises."],"supporting_citations":[{"why":"Provides the GRASP 2018 MCDHF code used to compute all level energies and transition rates.","marker":"[8]"},{"why":"Supplies the Kramida procedure used to estimate transition-rate uncertainties from the spread of high-intensity transitions.","marker":"[16]"},{"why":"Provides the FAC code whose independent results are compared with GRASP 2018 to assign B and C+ accuracy classes.","marker":"[17]"},{"why":"Supplies the experimental I8+ splitting used to benchmark the theory.","marker":"[4]"},{"why":"Supplies the experimentally deduced I9+ energy levels used to benchmark the $4d^8$ calculations.","marker":"[5]"},{"why":"Supplies an additional experimental value for the I8+ fine-structure splitting used in the comparison.","marker":"[18]"},{"why":"Documents the lack of published iodine-ion energy-level data in the NIST Atomic Spectra Database that motivates the gap-filling claim.","marker":"[3]"}],"fun_headline_variants":["Iodine ion energies fill missing data for I10+ to I12+","Energies for iodine ions I8+ to I12+ now complete","MCDHF-CI data fills iodine ion gap from I10+ to I12+","New iodine ion energies cover I10+ to I12+ gap"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"To assign 3% uncertainties to the unmeasured I10+ through I12+ levels, the paper assumes that the difference between the AS3 and AS2 active-space results bounds the true computational error, and that the accuracy measured for I8+ and I9+ against experiment carries over to the ions for which no measured spectra exist.","fun_headline_variants_meta":{"raw":{"variants":["Iodine ion energies fill missing data for I10+ to I12+","Energies for iodine ions I8+ to I12+ now complete","MCDHF-CI data fills iodine ion gap from I10+ to I12+","New iodine ion energies cover I10+ to I12+ gap"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000431,"raw_usage":{"total_tokens":2118,"prompt_tokens":783,"completion_tokens":1335,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":399,"completion_tokens_details":{"reasoning_tokens":1251}},"tokens_in":399,"tokens_out":1335,"duration_ms":10476,"temperature":1.0,"reasoning_tokens":1251,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T21:48:47.886981+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A clean test would be an electron-beam ion trap or beam-foil measurement of a predicted M1 line in I10+, I11+, or I12+: for example, an I10+ M1 transition predicted at 295.20 nm (level 12 to level 2) or the I12+ line predicted at 265.87 nm (level 2 to level 1). A measured wavelength that differs from the predicted value by more than the assigned 3% would show that the uncertainty estimate underestimates the true error, while agreement within the stated class would validate the tables as reference data.","supporting_citations":[{"cited_title":"Froese Fischer, G","cited_arxiv_id":null,"evidence_quote":"Provides the GRASP 2018 MCDHF code used to compute all level energies and transition rates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Kramida procedure used to estimate transition-rate uncertainties from the spread of high-intensity transitions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the experimental I8+ splitting used to benchmark the theory."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the experimentally deduced I9+ energy levels used to benchmark the $4d^8$ calculations."},{"cited_title":"Kimura and N","cited_arxiv_id":null,"evidence_quote":"Supplies an additional experimental value for the I8+ fine-structure splitting used in the comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the lack of published iodine-ion energy-level data in the NIST Atomic Spectra Database that motivates the gap-filling claim."}],"review_version":1}