{"id":"5c996bd0-5981-4b38-8b6d-186d3cf97047","arxiv_id":"1909.01523","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"New high-resolution spectra of tellurium-125 conversion electrons give penetration parameter λ=-1.2(6), mixing ratio |δ|=0.015(2), and a shake probability of about 50%, above the 20% prediction.","lead":"This paper measures the electrons emitted when the medical isotope iodine-125 decays, at higher energy resolution than before. It finds new values for two nuclear parameters and reports that a side effect called electron shake happens about twice as often as theory predicts.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The common-tail line-shape model is the load-bearing assumption: the 50% shake fraction is read directly from tail intensities with no error budget, and any shell dependence of the tails biases the peak-area ratios used for λ and |δ|.","rationale":"The reader's weakest assumption identifies exactly this point, and I agree. The 50% shake fraction is a direct output of the phenomenological tail model in Table I, and the paper explicitly lists tail/main overlap as a possible overestimate (Sec. V.B). Because the same tail template is imposed on all lines, any shell dependence of shake or extrinsic loss would bias the peak areas used in the λ–δ fit (Eq. 8) in a way that the current natural-width checks do not test. The proposed re-fit with a physical asymmetric tail profile and free tail intensities would show whether the ≈50% fraction is an artifact. If the tail fraction drops toward the 20% prediction, the shake discrepancy disappears; if it remains ~50%, the claim is supported. The literature-data exclusions are a second concern, but the tail/shape issue is more fundamental because it affects the experimental input itself. I therefore see no ground to move from the reader's CONDITIONAL verdict; the paper is careful and the central nuclear parameters are plausible, but the shake claim and, to a lesser extent, λ and δ need this robustness check.","tokens_in":16555,"tokens_out":8345,"duration_ms":86612,"concrete_test":"Refit the high-resolution L1 spectrum (Fig. 4) with the tail intensities as free parameters and an asymmetric tail profile (e.g., an exponential tail or the shake-energy distribution from Krause-Carlson [24]) instead of the fixed four-Gaussian set, allowing the Shirley background to float. Then propagate the best-fit L1 tail parameters into a refit of all low-resolution L/M/N spectra without enforcing identical relative tail intensities. If the total tail fraction (currently 1.1x the main peak) changes by more than ~0.2, or if the relative tail intensities change by more than 2σ, the 50% shake probability and the affected area ratios lack support.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantities depend on the assumption that every conversion-electron line has the same tail distribution, fixed from the high-resolution L1 line (Sec. III, Table I: tails at -5, -18, -42 eV with relative intensities 0.4, 0.5, 0.2). The paper's shake probability of 'about 50%' is literally the ratio of these tail intensities to the total peak area, quoted with no uncertainty. Section V.B itself concedes that (i) the large overlap between the tail and main peak may have overestimated the tail intensity, and (ii) the symmetric Gaussian tail shape is not the expected asymmetric shake distribution (Ref. [57]). If the tail fractions are not identical across subshells, or if the L1 fit is not a good template, the extracted areas for the weaker lines (L2, L3, M2, M3, N1) shift, and those ratios enter the least-squares fit of λ and |δ| through Eq. (8). The paper's checks vary the natural-width database but never vary the tail shape or the common-tail constraint, so this systematic is not bracketed. A wrong shake attribution (extrinsic energy loss rather than intrinsic shake) would further change the absolute tail strength and the 50% value.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports high-resolution conversion-electron spectroscopy of the electron-capture decay of 125I, using an electrostatic spectrometer, to determine nuclear parameters of the 35.5-keV M1+E2 transition in 125Te and to study electron shake. The authors fit the present L1:L2:L3, M1:M2:M3, and M1:N1 intensity ratios together with literature ICC data and angular-correlation data using Eqs. (7)-(9), obtaining λ = -1.2(6) and |δ| = 0.015(2) with reduced χ² = 1.2. They also estimate the shake probability from the tail intensity of the conversion lines as about 50%, roughly 2.5 times the calculated value of about 20%, and compare the measured λ with particle-vibrational and other nuclear models.","tokens_in":16842,"tokens_out":8857,"duration_ms":77694,"significance":"If correct, the results provide a precise, small negative penetration parameter for an l-forbidden M1 transition, which is in tension with core-polarization predictions, and they suggest a large enhancement of shake probability in an open-shell atom. The work also underpins Auger-yield determinations for the medical isotope 125I. The paper is commendable for a transparent fitting procedure, a reduced χ² of 1.2, a transmission check using the L1:M1 ratio that is insensitive to λ and δ, and for exploring the effect of different natural-width databases in the line-shape fits. However, the central claims rest on two assumptions that need further scrutiny: the common-tail line-shape model and the post-hoc rejection of discrepant literature data.","major_comments":[{"comment":"The least-squares fit excludes all literature values that are \"more than two standard deviations away from the corresponding fitted values,\" including the four most precise angular-correlation determinations of δ (rows marked \"c\": +0.09(1), +0.095(25), +0.078(12), +0.08(3)). Because the fit is itself used to define the reference values, this exclusion is post-hoc and circular; it removes precisely the data that disagree with the fitted |δ| = 0.015(2). The paper should report the fit with these points included, or apply a pre-defined robust rejection criterion, and should discuss the physical tension between the ICC ratios and the angular-correlation δ values. This is load-bearing because |δ| is a central claimed result.","section":"Section V.A, Table III, Eq. (9)"},{"comment":"All conversion-electron lines are fitted with tail parameters fixed to the L1 line, based on the assumption that all conversion lines share the same tail distribution. The extracted peak areas for L2, L3, M2, M3, and N1 therefore depend on the L1 tail shape and on the assumption that shake probabilities are equal across subshells. These areas enter the determination of λ and |δ| via Eq. (8). The paper varies the natural-width database but never varies the tail shape or relaxes the common-tail constraint. A test that fits the tail parameters for at least the stronger lines independently, or a systematic variation of the tail shape (e.g., an asymmetric tail), is needed to bracket this systematic. Without it, the quoted uncertainties on λ and |δ| are not complete.","section":"Section III, Table I, Section V.B, Eq. (8)"},{"comment":"The 50% shake probability has no uncertainty and no systematic budget. It is computed directly as (0.4+0.5+0.2)/(1+0.4+0.5+0.2) from the fixed tail intensities. The paper itself acknowledges that (i) the large overlap between the broad tails and the main peak may overestimate the tail intensity and (ii) the symmetric Gaussian tail shape is not the expected asymmetric shake distribution (Ref. [57]). These caveats need to be quantified. In addition, the comparison value of about 20% is taken from Ref. [23], a co-author's PhD thesis; the prediction should be documented in a peer-reviewed source or reproduced in the paper, and the measured shake probability should be quoted with an error bar.","section":"Section V.B"}],"minor_comments":[{"comment":"The reduced χ² = 1.2 is quoted without the number of degrees of freedom; please provide it to allow the reader to judge the fit quality.","section":"Section V.A"},{"comment":"The caption states \"The reduced χ2 of the this fit is 1.9\"; \"the this\" is a typo.","section":"Fig. 3 caption"},{"comment":"The text compares the tail shifts in Table I with the outer-shell binding energies, but tail #2 at -18 eV does not correspond directly to any of the listed N4/N5/O1/O2/O3 binding energies; a brief discussion of how the three discrete tail components map onto the expected shake distribution would be helpful.","section":"Section V.B"},{"comment":"The column headers \"Calculated λ=+2.4a |δ|=0.029|δ|=0.015\" are visually ambiguous; please clarify which column corresponds to which parameter set.","section":"Table III"},{"comment":"The phrase \"the deviation was fairy constant\" should read \"fairly constant\".","section":"Section II"}],"recommendation":"major_revision","confidential_remarks":"The most serious issue is the post-hoc exclusion of the precise angular-correlation δ measurements, which are exactly the data that disagree with the fitted |δ|. If the authors cannot provide a fit that includes these data and a frank discussion of the discrepancy, the paper's central claim is not supported. The reliance on a co-author's PhD thesis for the shake prediction should also be disclosed and ideally replaced by an independent reference."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThe short version: this is a careful experimental paper that gives new values for λ and |δ| in the 35.5-keV transition in 125Te, and it reports a shake probability about 2.5 times larger than single-configuration theory. The nuclear-parameter part is credible. The shake number is a rough estimate and should not be quoted without an uncertainty.\n\nWhat's actually new: high-resolution L1 line shape, new intensity ratios for L, M, N lines, and a re-evaluation that moves λ from +2.4 (1982) to −1.2(6) and |δ| from 0.029 to 0.015(2). The fit is transparent, reduced χ²=1.2, and the M1/L1 transmission check is a good idea. The PV-model calculation with the semi-empirical sign analysis is a reasonable attempt to interpret the negative sign, though it depends on calculated signs of matrix elements.\n\nSoft spots, in proportion:\n\n1. The common-tail line-shape model is load-bearing. All lines are fit with the same three tails taken from the high-resolution L1 line. The paper justifies this with similar predicted shake probabilities, but never varies the tail shape or tests shell-dependent tails. If the tails differ across subshells, the peak areas shift and λ and |δ| move. This is a real systematic and it isn't bracketed.\n\n2. The 50% shake probability is just the tail-to-total area ratio with no error budget. The paper itself concedes the overlap with the main peak may overestimate the tail intensity and that the symmetric Gaussian shape is not the expected asymmetric shake. So treat it as an order-of-magnitude estimate, not a measurement.\n\n3. The fit excludes data more than 2σ from the fitted curve. That includes four of the ten δ measurements and a few ICC ratios (not the entire angular-correlation set, as one summary claimed). The procedure can bias the result; a quick include/exclude sensitivity check would have settled it.\n\nThe citation pattern is fine; the co-author thesis is a comparison, not an input to the fit.\n\nBottom line: this deserves serious referee time. It is a solid update for nuclear data evaluations and for anyone using 125I conversion yields. The shake part needs more work before it can be used quantitatively.\n\nRecommendation: I would accept it with minor revisions, asking for a sensitivity analysis on the excluded data and an uncertainty on the shake fraction.\n\nBest,\n[Your name]","headline":"A careful new measurement that credibly updates λ and |δ| for 125Te; the 50% shake fraction is a rough estimate that needs an error budget.","tokens_in":17406,"tokens_out":5609,"would_cite":true,"duration_ms":48005,"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":"High-resolution spectroscopy of the 125I decay finds a 50% shake probability for its conversion electrons, 2.5 times the single-configuration prediction, and yields nuclear parameters $\\lambda=-1.2(6)$ and $|\\delta|=0.015(2)$ for the…","keywords":["125I electron-capture decay","conversion electron spectroscopy","penetration parameter","internal conversion coefficients","E2/M1 mixing ratio","shake probability","125Te","particle-vibrational model"],"falsifier":"A coincidence experiment that detects a shake electron at the same time as the 35.5-keV conversion electron would separate true shake from surface-plasmon energy loss: if the low-energy tail does not require a second electron, the 50% shake assignment collapses. Likewise, if high-resolution spectra of the $L_1$, $M_1$, and $N_1$ lines, recorded with better statistics, demand different tail parameters for different shells, the common-tail assumption, and with it the fitted $\\lambda$ and $|\\delta|$, would need revision.","tokens_in":16363,"feed_emoji":"⚡️","tokens_out":7376,"duration_ms":65474,"temperature":0.7,"pith_summary":"The paper reports high-resolution measurements of the conversion electrons emitted when the 35.5-keV excited state of $^{125}$Te is populated by electron capture in $^{125}$I. By fitting the measured $L$, $M$, and $N$ sub-shell intensities together with literature values, it extracts a penetration parameter $\\lambda = -1.2(6)$ and an E2/M1 mixing ratio $|\\delta| = 0.015(2)$ for the 35.5-keV transition. The small negative $\\lambda$ means the nuclear-penetration correction to the internal conversion coefficients is under 4%. The line-shape analysis also finds that about 50% of the conversion-electron intensity sits in low-energy shake tails, more than twice the 20% single-configuration prediction. A sympathetic reader would care because the result sharpens a long-standing test of how nuclear currents enter internal conversion and because $^{125}$I is a widely used medical isotope whose Auger and conversion-electron yields are needed for dosimetry.","feed_headline":"Electron shake in 125I decay is 2.5 times predicted","feed_subtitle":"High-resolution spectra also pin the 35.5-keV transition's penetration parameter at -1.2(6), a small negative anomaly.","key_machinery":"The argument runs through two identities. First, the conversion coefficient of each sub-shell is written as $\\alpha_i = \\alpha_i(M1)(1+b_1(i)\\lambda+b_2(i)\\lambda^2)$ for the $M1$ part, combined with the $E2$ part through $|\\delta|$; the coefficients $b_1,b_2$ come from Dirac-Hartree-Fock-Slater wavefunctions. Second, every conversion line is fitted with a Lorentzian core plus four Gaussians, with the three low-energy tail components' shifts, widths, and intensities fixed from the high-resolution $L_1$ line and then applied to all other lines. The tail areas are identified with shake electrons, whose energy distribution is tied to outer-shell binding energies by the Krause-Carlson relation. The fit uses the measured sub-shell intensity ratios plus literature conversion coefficients to constrain $\\lambda$ and $\\delta$ in a least-squares sense.","core_discovery":"The paper's central claim is that the high-resolution conversion-electron spectrum of the 35.5-keV $M1+E2$ transition in $^{125}$Te can be described by a common line shape whose low-energy tail is mostly atomic shake, and that with this line shape the measured sub-shell ratios, combined with earlier data, determine $\\lambda = -1.2(6)$ and $|\\delta| = 0.015(2)$. The fitted tail intensity corresponds to a shake probability of about 50%, which is 2.5 times the predicted value of 20% from single-configuration calculations. The magnitude of $\\lambda$ is close to the value obtained by combining the particle-vibrational model's allowed penetration matrix elements with the experimental forbidden $M1$ gamma matrix element, while its negative sign contradicts the positive sign predicted by core-polarization theory.","pith_inferences":["If the common-tail assumption is checked and survives, the 50% shake fraction would be a direct experimental benchmark for correlation-inclusive shake calculations on open-shell atoms, not just for Te.","The negative sign of $\\lambda$ in $^{125}$Te, taken with the signs reported for $^{121}$Te and $^{123}$Te, suggests that core-polarization theories giving a positive sign for these $\\nu d_{3/2}\\to\\nu s_{1/2}$ transitions should be revisited; a high-resolution remeasurement of the other two isotopes at similar precision would test that pattern.","Because the published Auger yields for $^{125}$I are normalized through the conversion-electron intensities used in this paper, the shift from the older $\\lambda=+2.4$, $|\\delta|=0.029$ evaluation to the new values changes those absolute yields at the few-percent level, which matters for dosimetry models of this medical isotope."],"forward_implications":["The 35.5-keV transition's conversion coefficients deviate from the no-penetration values by less than 4%, so the earlier evaluation with a larger positive penetration anomaly is not supported.","The measured $|\\delta|=0.015(2)$ agrees with the Kisslinger-Sörensen-model prediction, while the sign of $\\delta$ remains undetermined by conversion-coefficient data alone.","The measured $L_1:M_1$ and $M_1:N_1$ ratios agree with theoretical internal conversion coefficients, confirming that the electrostatic spectrometer's transmission is effectively energy-independent over the measured range.","A shake probability near 50% means the low-energy tails carry as much intensity as the main conversion lines, so any absolute conversion-electron or Auger-yield determination for $^{125}$I must include this tail intensity.","For l-forbidden $M1$ transitions generally, the paper's analysis suggests that the magnitude of the penetration parameter can be estimated by combining a theoretical calculation of the allowed penetration matrix elements with the experimental forbidden $M1$ gamma matrix element."],"supporting_citations":[{"why":"Introduces the penetration effect and the l-forbidden mechanism that motivates the measurement.","marker":"[1]"},{"why":"Supplies the natural (Lorentzian) widths used to fit each conversion line.","marker":"[16]"},{"why":"Provides the previous evaluated $\\lambda$ and $|\\delta|$ and the magnetic-spectrometer ratios this work compares and partly excludes.","marker":"[18]"},{"why":"Gives the single-configuration shake-probability prediction of about 20% that the measured 50% is compared against.","marker":"[23]"},{"why":"Provides the Krause-Carlson methodology used to compute shake probabilities and the energy relation for shake electrons.","marker":"[24]"},{"why":"Supplies BrIcc conversion coefficients and binding energies used in the fitting and in the theoretical ICC values.","marker":"[25]"},{"why":"The CATAR code, modified, is used to calculate the penetration coefficients $b_1$ and $b_2$.","marker":"[32]"},{"why":"The least-squares minimizer used for the extraction of $\\lambda$ and $|\\delta|$.","marker":"[33]"}],"fun_headline_variants":["Electron shake in 125I decay is 2.5× predicted","125I spectra show electron shake 2.5× theory","Conversion electron shake 2.5× expected in 125I decay","125I decay electron shake exceeds prediction by 2.5×","High-resolution 125I data: shake 2.5× predicted"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes that all conversion-electron lines share exactly the same low-energy tail shape, fixed from the $L_1$ line, and that this tail is mostly atomic shake rather than electrons that lost energy while leaving the source.","fun_headline_variants_meta":{"raw":{"variants":["Electron shake in 125I decay is 2.5× predicted","125I spectra show electron shake 2.5× theory","Conversion electron shake 2.5× expected in 125I decay","125I decay electron shake exceeds prediction by 2.5×","High-resolution 125I data: shake 2.5× predicted"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000343,"raw_usage":{"total_tokens":1820,"prompt_tokens":812,"completion_tokens":1008,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":428,"completion_tokens_details":{"reasoning_tokens":914}},"tokens_in":428,"tokens_out":1008,"duration_ms":9588,"temperature":1.0,"reasoning_tokens":914,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:15:46.726961+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A coincidence experiment that detects a shake electron at the same time as the 35.5-keV conversion electron would separate true shake from surface-plasmon energy loss: if the low-energy tail does not require a second electron, the 50% shake assignment collapses. Likewise, if high-resolution spectra of the $L_1$, $M_1$, and $N_1$ lines, recorded with better statistics, demand different tail parameters for different shells, the common-tail assumption, and with it the fitted $\\lambda$ and $|\\delta|$, would need revision.","supporting_citations":[{"cited_title":"Church and J","cited_arxiv_id":null,"evidence_quote":"Introduces the penetration effect and the l-forbidden mechanism that motivates the measurement."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the natural (Lorentzian) widths used to fit each conversion line."},{"cited_title":"[19], where both studies used mag- netic spectrometers","cited_arxiv_id":null,"evidence_quote":"Provides the previous evaluated $\\lambda$ and $|\\delta|$ and the magnetic-spectrometer ratios this work compares and partly excludes."},{"cited_title":"Brabec, M","cited_arxiv_id":null,"evidence_quote":"Gives the single-configuration shake-probability prediction of about 20% that the measured 50% is compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Krause-Carlson methodology used to compute shake probabilities and the energy relation for shake electrons."},{"cited_title":"Kib´ edi, T","cited_arxiv_id":null,"evidence_quote":"Supplies BrIcc conversion coefficients and binding energies used in the fitting and in the theoretical ICC values."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The CATAR code, modified, is used to calculate the penetration coefficients $b_1$ and $b_2$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The least-squares minimizer used for the extraction of $\\lambda$ and $|\\delta|$."}],"review_version":1}