{"id":"bf5d4d2e-50e3-423b-9902-24ca363e6019","arxiv_id":"2412.05932","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Improved potassium isotope-shift factors and updated nuclear polarizability corrections give a more precise radius for 38mK and an isospin-breaking test in the A=38 isotriplet that is consistent with zero.","lead":"This paper computes extremely precise atomic factors for potassium and uses them to extract the nuclear size of a rare isotope, 38mK, then tests whether the charge radii of three related nuclei obey isospin symmetry. The test finds isospin-breaking effects consistent with zero, offering a benchmark for nuclear calculations that feed the precision extraction of the Cabibbo-Kobayashi-Maskawa matrix element Vud.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'most stringent' isospin-breaking test depends on an assumed mirror charge-radius relation for 38Ar; if that relation fails, the test loses its claimed directness.","rationale":"The reader's weakest_assumption identifies the same load-bearing concern: the stringent ISB test assumes the empirical mirror charge-radius relation supplies the poorly measured 38Ar radius. This is indeed the most load-bearing point for the central claim. The paper has independent support for the atomic part: energies agree with experiment, hyperfine benchmarks are at the percent level, and the FF results are cross-checked against semi-empirical mass shifts from muonic data. The updated nuclear polarization correction is promised in a separate publication and does affect absolute radii, but footnote [75] notes that it shifts all members of the isotriplet nearly identically, so it is not decisive for ΔM_B^(1). The mirror relation, by contrast, directly enters the construction of the 'stringent' test. If the direct 38Ar radius were used, the uncertainty would grow and the phrase 'most stringent test' would no longer be justified. Because the reader already conditioned the verdict on this and other missing pieces, my read does not change the verdict; it sharpens the specific condition that must be checked. The proposed test replaces the assumed radius with the direct value and asks whether the qualitative conclusion survives; that settles whether the concern is decisive or merely a framing issue.","tokens_in":35085,"tokens_out":3600,"duration_ms":36194,"concrete_test":"Recompute ΔM_B^(1)(38) using the previously recommended direct value r(38Ar) = 3.402(6) fm from Ref. [29] in place of the mirror-relation value 3.3973(36) fm, propagating all correlations as in Ref. [75]. Compare the resulting central value and uncertainty with the model band -0.42 to -0.04 fm^2. If the direct-radius result remains inside the band with an uncertainty around 0.7 fm^2, the qualitative ISB conclusion survives but the 'most stringent test' claim must be revised; if it moves outside the band, the mirror assumption is decisive and the claimed stringent test is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline ISB result ΔM_B^(1)(38) = -0.48(63) fm^2 is presented as the most stringent charge-radius test of isospin symmetry breaking. But the stringent version is not built from a directly measured 38Ar radius. Section VII states, 'To perform a more stringent test, we assume that the mirror fit holds in this region of the nuclear chart,' giving r_N(38Ca) = r_N(38Ar) + 0.0727(18) fm. The direct 38Ar radius quoted from Ref. [29], r(38Ar) = 3.402(6) fm, is replaced by the mirror-relation value 3.3973(36) fm. The quoted -0.48(63) fm^2 therefore contains a nuclear-model assumption. If the mirror relation is inaccurate at A=38, the central value shifts and the uncertainty underestimates the model dependence; the result would no longer be a direct experimental measurement of ISB, and the 'stringent benchmark' framing would need tempering. The paper is transparent about the assumption, but the abstract and summary present the resulting ΔM_B as the main finding without carrying the caveat. This is not an internal inconsistency; it is a limitation on the load-bearing interpretation of the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports relativistic coupled-cluster (RCC) calculations of isotope-shift factors for seven low-lying states of potassium, comparing finite-field, expectation-value, and analytical-response approaches. The finite-field results are selected as the most reliable and are benchmarked against experimental energies, hyperfine constants, muonic radius differences, and semi-empirical mass shifts. Using the recommended field-shift and mass-shift factors, the authors extract the mean-square radius difference between 38mK and 39K, update the absolute radius of 38mK with a revised nuclear-polarization correction, and combine this with an updated 38Ca radius to obtain ΔM_B^(1)(38) = -0.48(63) fm², which they interpret as the most stringent charge-radius-based test of isospin symmetry breaking to date.","tokens_in":35212,"tokens_out":5760,"duration_ms":61360,"significance":"If the central result stands, it provides a valuable experimental benchmark for the isospin-symmetry-breaking corrections needed in superallowed beta-decay studies and for the extraction of V_ud. The paper's strengths include a systematic comparison of three many-body approaches, careful benchmarking against hyperfine and muonic data, explicit treatment of triples, basis, Breit, and QED corrections, and a transparent decomposition of uncertainties. The FF-based mass-shift factors are cross-checked against semi-empirical values, which lends credibility to the atomic theory. However, the headline isospin-breaking test rests on an assumed mirror charge-radius relation for 38Ar and an unpublished nuclear-polarization calculation; these caveats need to be carried explicitly in the abstract and summary, and the sensitivity of the result to those assumptions should be quantified.","major_comments":[{"comment":"The stringent test uses r_N(38Ar) not from the direct measurement r(38Ar)=3.402(6) fm quoted from Ref. [29], but from the assumed mirror relation r_N(38Ca)=r_N(38Ar)+0.0727(18) fm, as stated in the sentence 'To perform a more stringent test, we assume that the mirror fit holds in this region of the nuclear chart.' The quoted ΔM_B^(1)(38) = -0.48(63) fm² therefore contains an untested nuclear-model assumption, and the abstract and summary present this as the 'most stringent test' without carrying that caveat. Since this is the paper's central claim, please quantify the sensitivity of ΔM_B to the mirror-relation slope and to plausible local deviations at A=38, and rephrase the abstract and summary to distinguish the direct from the assumption-dependent determination.","section":"§VII, Eq. (53)"},{"comment":"The updated absolute radii depend on nuclear-polarization shifts that are asserted but not derived in the manuscript: ΔE_NP(39K)=156(47) eV and ΔE_NP(40Ca)=177 eV, attributed to an additional 'nucleon polarization' effect, with the text stating that details 'will be published elsewhere.' Because these shifts change r_N(39K), r_N(38mK), and r_N(38Ca), the reported radius updates are not reproducible from the information provided. At minimum, include the calculation in an appendix or clearly state which reported quantities are preliminary pending the separate publication. Footnote [75] indicates that ΔM_B itself is only slightly affected because the shift is common to the isotriplet, but the absolute radii are also primary outputs of the paper.","section":"§VI and §VII"},{"comment":"The semi-empirical total mass-shift factor K_tot_SE is constructed using the measured isotope shifts, the muonic δ⟨r²⟩^{39,41}, and this paper's own F factors from Table IV. The weighted-average K_tot_WA then combines K_tot_SE with K_tot_FF, which share the same atomic wavefunctions and, in particular, the same field-shift factor; the two quantities are therefore not fully independent. Please state the assumed correlation structure and, if the values are simply weighted by their quoted uncertainties, justify that treatment, since it affects the uncertainty on ⟨r²⟩^{38m,39} and hence on r_N(38mK).","section":"§V, Table VIII"}],"minor_comments":[{"comment":"The header for the field-shift block reads 'F values (in MHz/fm^{-1})'; the units should be MHz/fm^{-2}.","section":"Table VI"},{"comment":"The subscripts Z_{+1}, Z_0, and Z_{-1} should be defined before Eq. (53) is used, or the isospin assignments T3 = ±1, 0 should be stated explicitly.","section":"Eq. (53)"},{"comment":"The phrase 'most stringent test' should be qualified whenever it is used, since the stringent form of the test relies on the assumed mirror relation for 38Ar.","section":"Abstract and §IX"},{"comment":"The paper states that second-order contributions in the finite-field method are 'typically considered negligible' and that this assumption might differ for heavier elements; this limitation should be restated in the conclusions or in the uncertainty discussion, as it is worth carrying forward to applications beyond potassium.","section":"§IV B and Table IV"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Re: arXiv:2412.05932. The thing to know: this is the most careful potassium isotope-shift calculation to date, and it delivers a more accurate radius for 38mK plus an isospin-breaking test consistent with zero. The atomic physics is the real content; the beta-decay framing is the payoff but also where the caveats live.\n\nWhat's new: first IS factors for seven low-lying K states at RCCSDT with triples, and a genuine three-method comparison (finite-field, expectation-value, analytical response). They benchmark energies against NIST to ~0.1% and hyperfine constants to ~1%, show the FF method matches semi-empirical mass shifts best, and explain why via orbital relaxation. The weighted radius difference <r^2>^{38m,39} = -0.019(14) fm^2 is twice as accurate as the previous value. That part is solid.\n\nSoft spots, in proportion. First, the absolute radii rely on an updated nuclear-polarization correction that is not derived here; it is a separate publication in preparation. You cannot reproduce r(38mK)=3.4396(37) fm from this manuscript alone. That is a real gap for a paper whose abstract leads with the radius. Second, the \"most stringent test\" claim deserves tempering. The stringent version replaces the poorly measured 38Ar radius with a mirror-fit relation from Ref. [29]; the paper is transparent about this in Section VII, but the abstract and summary do not carry the caveat. If that relation is inaccurate at A=38, both the central value and uncertainty shift, and the result becomes a benchmark for the relation as much as a direct measurement. The quoted uncertainty (0.63 fm^2) also exceeds the model spread it is meant to discriminate, so \"stringent\" overstates the discriminating power. Third and minor: the semi-empirical K_tot_SE is built with the paper's own F, so the weighted average of SE and FF factors is not fully independent. Given the cross-checks, this is likely a small effect, but it should be stated.\n\nNothing here is load-bearing wrong. The atomic calculations are extensive, the uncertainty estimates are transparent, and the central ISB result is consistent with zero, which is what the models expect. A referee should ask for the nuclear-polarization derivation to be included or linked to a public preprint, and for the abstract to note that the stringent test assumes the mirror radius relation.\n\nWho for: atomic physicists working on isotope shifts, and the superallowed-beta-decay community. It deserves a serious referee; I would send it out. I would cite it for the K factors, and it is a good reading-group paper for the three-method comparison alone.","headline":"Solid RCCSDT isotope-shift factors for K with an honest three-method comparison, but the headline ISB benchmark rests on an assumed mirror-radius relation and a nuclear-polarization calculation deferred to another paper.","tokens_in":35957,"tokens_out":2544,"would_cite":true,"duration_ms":27135,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Potassium isotope-shift calculation puts isospin-breaking test at –0.48(63) fm², consistent with zero.","keywords":["isotope shift","potassium","relativistic coupled-cluster theory","finite-field method","nuclear charge radius","isospin symmetry breaking","superallowed beta decay","Vud"],"falsifier":"Measure the charge radius of 38Ar directly, via electron scattering on a radioactive 38Ar target or muonic-atom x-ray spectroscopy, and compare the extracted r_N(38Ar) with the mirror-relation prediction; a deviation beyond the quoted 0.0018 fm uncertainty would change ΔM_B^(1)(38) by more than its stated error.","tokens_in":34723,"feed_emoji":"⚛️","tokens_out":4732,"duration_ms":42918,"temperature":0.7,"pith_summary":"This paper aims to sharpen the extraction of nuclear charge radii for the A = 38 isotriplet (38Ca, 38mK, 38Ar) by calculating isotope-shift factors of potassium with relativistic coupled-cluster theory, and to use these radii as a test of isospin symmetry breaking. The authors argue that the finite-field approach gives the most reliable isotope-shift factors, because it naturally includes orbital relaxation effects. Combining new isotope-shift factors with muonic-atom radii and an updated nuclear-polarization correction, they obtain the isospin-breaking quantity ΔM_B^(1)(38) = −0.48(63) fm², which is compatible with zero. If correct, this is the most stringent charge-radius-based test of isospin symmetry breaking to date and a benchmark for the nuclear-model corrections needed to extract the CKM matrix element Vud from superallowed $\\beta$ decays.","feed_headline":"New K radii: isospin breaking in A=38 is zero within error","feed_subtitle":"Precise isotope-shift calculations benchmark the corrections needed to extract Vud from superallowed beta decays.","key_machinery":"The load-bearing object is the isospin-breaking matrix element ΔM_B^(1) ≈ ½(Z+1 r²_{N,+1} + Z-1 r²_{N,-1}) − Z0 r²_{N,0}, built from charge radii of the isotriplet. The argument also rests on isotope-shift factors F, K_NMS, K_SMS computed with the relativistic coupled-cluster method in singles-doubles-triples approximation, evaluated via finite-field, expectation-value, and analytical-response approaches; the finite-field method is selected because it captures orbital relaxation. The updated nuclear-polarization correction to the muonic reference radius of 39K and the mirror relation r_N(38Ca) = r_N(38Ar) + 0.0727(18) fm carry the final step.","core_discovery":"The central claim is that the isospin-symmetry-breaking combination of the A = 38 triplet charge radii, ΔM_B^(1)(38) = −0.48(63) fm², is consistent with zero within uncertainty, and that this constitutes the most stringent test of isospin symmetry breaking using charge radii. The paper also claims that the finite-field method yields isotope-shift factors for potassium that agree with semi-empirical data, and that the previous nuclear-polarization correction to the 39K reference radius was underestimated by about 37 eV due to a previously overlooked nucleon-polarization contribution.","pith_inferences":["If the mirror relation holds generally, the same method could be applied to other isobaric triplets to map isospin breaking across the nuclear chart.","The overlooked nucleon-polarization contribution may shift absolute radii in other medium-mass muonic atoms, not just potassium, suggesting that previously published reference radii in that region could need revision.","A direct measurement of the 38Ar charge radius (e.g., electron scattering or a muonic-atom measurement) would remove the dependence on the mirror-fit assumption and convert the stringent test into a model-independent one.","The technique of combining finite-field isotope-shift factors with semi-empirical mass shifts could be extended to other alkali atoms with limited stable isotopes."],"forward_implications":["The result provides a rigorous experimental benchmark for nuclear-model calculations of the isospin-breaking correction δC in superallowed beta decays, a crucial input to Vud.","The compatibility of ΔM_B^(1)(38) with zero indicates that the charge radii of the A = 38 isotriplet do not show evidence of significant isospin breaking at current precision.","The updated r_N(38mK) = 3.4396(37) fm differs by about one standard error from the previously recommended value, affecting the statistical rate function f for 38mK → 38Ar decay.","The discrepancy for the A = 26 isotriplet (4.5σ) motivates a reanalysis of the 26mAl charge radius.","The improved isotope-shift factors are twice as accurate as previous values for the 38mK–39K radius difference, opening the way for more precise optical measurements."],"supporting_citations":[{"why":"Provides the mirror charge-radius fit and the A = 38 radii used in the final isospin-breaking test.","marker":"[29]"},{"why":"Muonic-atom x-ray measurements supply the reference charge radii of 39K and 41K.","marker":"[50]"},{"why":"Magneto-optical trap measurement of the 38mK–39K isotope shift for the 4S–4P3/2 transition.","marker":"[57]"},{"why":"Collinear laser spectroscopy of the 38mK–39K isotope shift, the previous radius extraction.","marker":"[58]"},{"why":"Collinear laser spectroscopy measurement of the 38Ca–40Ca isotope shift used to update the 38Ca radius.","marker":"[74]"},{"why":"Defines the isospin-breaking quantity ΔM_B^(1) and the benchmark idea connecting charge radii to δC.","marker":"[16]"},{"why":"Provides the baseline nuclear-polarization correction and the reference radius of 39K that the paper updates.","marker":"[62]"},{"why":"Earlier isotope-shift factors and radius extraction for exotic potassium isotopes, providing the comparison point for accuracy.","marker":"[28]"}],"fun_headline_variants":["Isospin breaking in A=38 compatible with zero","Potassium isotope shifts: isospin symmetry holds","Most precise isospin test from K radii","Finite-field method yields accurate K IS factors","Zero isospin breaking in A=38 from new calculations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The stringent isospin-breaking test assumes that the mirror charge-radius relation r_N(38Ca) = r_N(38Ar) + 0.0727(18) fm, which comes from a fit to other mirror pairs, holds at A = 38; if the 38Ar radius deviates from this relation, the quoted ΔM_B value is not a direct measurement.","fun_headline_variants_meta":{"raw":{"variants":["Isospin breaking in A=38 compatible with zero","Potassium isotope shifts: isospin symmetry holds","Most precise isospin test from K radii","Finite-field method yields accurate K IS factors","Zero isospin breaking in A=38 from new calculations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000209,"raw_usage":{"total_tokens":1400,"prompt_tokens":931,"completion_tokens":469,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":547,"completion_tokens_details":{"reasoning_tokens":392}},"tokens_in":547,"tokens_out":469,"duration_ms":4852,"temperature":1.0,"reasoning_tokens":392,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T20:11:51.743897+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the charge radius of 38Ar directly, via electron scattering on a radioactive 38Ar target or muonic-atom x-ray spectroscopy, and compare the extracted r_N(38Ar) with the mirror-relation prediction; a deviation beyond the quoted 0.0018 fm uncertainty would change ΔM_B^(1)(38) by more than its stated error.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the mirror charge-radius fit and the A = 38 radii used in the final isospin-breaking test."},{"cited_title":"Ohayon, R","cited_arxiv_id":null,"evidence_quote":"Muonic-atom x-ray measurements supply the reference charge radii of 39K and 41K."},{"cited_title":"Falke, E","cited_arxiv_id":null,"evidence_quote":"Magneto-optical trap measurement of the 38mK–39K isotope shift for the 4S–4P3/2 transition."},{"cited_title":"Halloran, S","cited_arxiv_id":null,"evidence_quote":"Collinear laser spectroscopy of the 38mK–39K isotope shift, the previous radius extraction."},{"cited_title":"Satula, J","cited_arxiv_id":null,"evidence_quote":"Collinear laser spectroscopy measurement of the 38Ca–40Ca isotope shift used to update the 38Ca radius."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the isospin-breaking quantity ΔM_B^(1) and the benchmark idea connecting charge radii to δC."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the baseline nuclear-polarization correction and the reference radius of 39K that the paper updates."},{"cited_title":"Koszorus, X","cited_arxiv_id":null,"evidence_quote":"Earlier isotope-shift factors and radius extraction for exotic potassium isotopes, providing the comparison point for accuracy."}],"review_version":1}