{"id":"d4e0a84b-10e8-4cf5-81c6-915d4e3cf6c1","arxiv_id":"2508.18562","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Including subthreshold atomic hole states in the spectral function of 159Dy electron capture raises the predicted rate in the neutrino endpoint region by more than an order of magnitude.","lead":"Electron capture in 159Dy is usually modeled with only energetically allowed atomic hole states. This paper adds the tails of energetically forbidden subthreshold hole states and finds the predicted decay rate near the endpoint rises by more than an order of magnitude, which could improve the case for a neutrino mass experiment with this isotope.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Strongest quantitative claim depends on the unvalidated Lorentzian tail of Eq. (8) at tens-to-hundreds of widths from resonance; the paper itself flags this approximation as requiring improvement.","rationale":"The reader's weakest-assumption analysis identifies exactly the same load-bearing point: the quantitative enhancement rests on the Lorentzian tail of Eq. (8) far from resonance. My independent reading agrees. The paper's arithmetic is internally consistent and the qualitative idea of subthreshold-state contributions is plausible; the self-reported limitation on Eq. (8) is an in-scope warning that the quantitative claim is not yet secured. This justifies the CONDITIONAL verdict already reached, so I do not recommend changing it. The concern is not that the model is internally contradictory or that the authors acted improperly; rather, the empirical or many-body validation of the spectral function far from resonance is missing, and the paper itself acknowledges this. The proposed test using the measured 163Ho spectrum is a practical way to settle whether the tail shape assumed in Eq. (8) survives at the relevant offset.","tokens_in":7026,"tokens_out":5369,"duration_ms":60767,"concrete_test":"Use the high-resolution 163Ho electron-capture spectrum of Velte et al. (2019) to extract the empirical spectral function on the side of the M1 resonance corresponding to the same reduced offset as in 159Dy, i.e., E_ex - eps_M1 about 0.8 keV (roughly 60 natural widths for the M1 line). Fit that region without imposing the Lorentzian form of Eq. (8), then insert the extracted tail shape into the 159Dy calculation and recompute the two r values in Table II. If the empirical tail at this offset differs from the Lorentzian prediction by more than a factor of about 3, the declared order-of-magnitude enhancement is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is the factor ~10 enhancement of the endpoint-region EC rate for 159Dy in Table II (r = 3.3e-11 vs 2.9e-12 for m_nu = 1 eV; 8.4e-12 vs 7.3e-13 for m_nu = 0.1 eV). These numbers follow from evaluating P_x(E_ex) in Eq. (8) as a Lorentzian with an energy-independent width at large offsets from the resonance centers: for Q - Delta_nucl = 1.18 keV, the M1 line at 1.968 keV is 0.788 keV (about 60 times its 13 eV width) away, and the M2 line at 1.768 keV is 0.588 keV (about 100 times its 5.8 eV width) away. For the older Q value 1.7 keV, the M2 offset is only 68 eV. At such offsets, the Lorentzian tail is not independently established: Auger decay channels, Fano-type interference, or an underlying non-resonant continuum could change P(E_ex) substantially. The authors explicitly state that the estimation of the spectral function using Eq. (8) requires improvements, particularly the resonance energies and the energy-independent parameterization of the width. Because the claimed enhancement is an order of magnitude, a line-shape correction of even a factor of a few could weaken or erase the headline result.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript studies the electron-capture (EC) spectrum of 159Dy and argues that atomic hole states whose resonance energies lie above the decay endpoint (subthreshold states) contribute through the tails of their line shapes to the partial rate in the neutrino-mass-sensitive endpoint region. Using a Lorentzian parameterization of the spectral function with parameters taken from atomic compilations, the authors compute the fraction r of the total EC rate in a 1 eV window above the minimum neutrino energy for the proposed Q value of 1.18 keV and for an older value of 1.7 keV. They report an enhancement of r by more than an order of magnitude when subthreshold M states are included, and they argue that the spectral function can be determined experimentally. They also mention a similar enhancement for 111In.","tokens_in":7323,"tokens_out":8142,"duration_ms":73725,"significance":"The idea that subthreshold atomic states should not be discarded in EC endpoint analyses is physically interesting and, if quantitatively correct, would improve the prospects of 159Dy as a neutrino-mass candidate by increasing the endpoint-event rate. The authors are to be credited for using external atomic parameters rather than fitting the 159Dy data, which avoids circularity, and for transparently flagging the approximations in Eq. (8). The central numerical claim, however, rests entirely on the unvalidated Lorentzian tail of the M1/M2 lines at energy offsets of tens to hundreds of line widths, and the paper does not yet provide the sensitivity analysis needed to establish the claimed enhancement. The manuscript is clearly written and the formalism is standard, but the quantitative results require substantial strengthening before publication.","major_comments":[{"comment":"The headline enhancement of r by factors of about 11 and 300 (Table II) is obtained by evaluating Eq. (8) at E_ex = 1.18 keV (or 1.7 keV), which lies 0.79 keV below the M1 line at 1.968 keV and 0.59 keV below the M2 line at 1.768 keV (for the 1.7 keV Q, the M2 offset is only 68 eV). These offsets correspond to about 60 and 120 times the stated widths, where the Lorentzian shape has no empirical or ab-initio support. The authors explicitly write that 'the estimation of the spectral function using Eq. (8) requires improvements, particularly regarding the resonance energies ... and the energy-independent parameterization of the width.' Since the enhancement factor is directly proportional to the tail value, a modest change in the line shape (e.g., exponential falloff or an additional non-resonant background) could reduce or erase the claimed effect. I request a quantitative sensitivity study that varies the tail parameterization and shows the resulting range of r in Table II.","section":"Eq. (8), Table II"},{"comment":"The only empirical validation offered is the statement that the measured 163Ho spectrum shows that P(E_ex) is 'flat' near the endpoint. No comparison plot or fit is shown, and the Ho endpoint lies essentially at the M-line complex, whereas the 159Dy endpoint is 0.6–0.8 keV below the M lines. The Ho data, as presented, cannot test the far tail that drives the 159Dy enhancement. Please quantify the agreement of the Lorentzian model with the Ho data over an energy range comparable to the offsets relevant for 159Dy, or otherwise qualify the validation claim.","section":"Sec. IV (validation with 163Ho)"},{"comment":"The numerical results in Table II and Fig. 2 require the full spectral function including continuum N, O, and P states (the dashed curves), but the parameters for the O and P shells (and for the K and L shells shown in Fig. 1) are not provided. Only M1, M2, N1, N2 are listed in Table I, and the text refers to L states in the 111In discussion without giving their parameters. Without a complete parameter table, the computation is not reproducible and the central values of r cannot be verified. Please include all parameters used, or make the calculation script available.","section":"Table I and Fig. 2"}],"minor_comments":[{"comment":"The caption states 'The parameters of ϵx and Γx correspond to Tb, while βx pertains to Dy.' Since the spectral function is for the daughter atom (Tb) but the captured electron is initially in Dy, one sentence explaining why the parent and daughter parameters are mixed would help readers.","section":"Table I caption"},{"comment":"The factor 2 multiplying n_x is not defined. If n_x is the occupation number of the subshell, the factor 2 is redundant; if it is the number of electrons, it should be stated. This also affects the normalization of P(E_ex).","section":"Eq. (8)"},{"comment":"Refs. [19, 20] are the AME 2020 mass evaluation, not a direct measurement; the phrase 'a previous measurement of the Q value' is misleading and should be rephrased.","section":"References [19, 20]"},{"comment":"The claim that 'r is substantially enhanced by approximately 10^4' for 111In is not quantified in Table II or any other table, and the parameters for the L2 state are not given. Please provide the calculation details or move this to future work.","section":"111In discussion"},{"comment":"The statement 'no adjustments are required for previous findings concerning 163Ho' is too strong; the model makes a prediction for 163Ho that could be tested against the high-resolution data of ref. [10], and a comparison would strengthen the paper.","section":"Conclusions"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a short Letter-style paper. The referees should ask the authors to supply the missing parameter table and a sensitivity study of the line-shape model; these are necessary to support the central claim. The archival value of the paper would be increased if the calculation code were released."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this if you care about future neutrino-mass experiments with low-Q EC isotopes. The new thing here is simple: the authors include energetically forbidden (subthreshold) atomic hole states in the spectral function for 159Dy EC, and find the rate in the last eV near the endpoint is boosted by more than an order of magnitude. That is a genuinely new quantitative prediction, and it runs contrary to the usual assumption that only continuum states matter.\n\nWhat the paper does well: the calculation is transparent and reproducible from the tabulated M- and N-shell parameters; the arithmetic is consistent; the comparison with 163Ho data is a sensible, if indirect, consistency check; and the suggestion that P(E) can be mapped out using different Q values is a constructive experimental idea. The authors also state openly that Eq. (8) is an approximation that needs improvement.\n\nThe soft spot is the one they flag themselves. The enhancement factors 11 and 302 come from evaluating a Lorentzian tail at offsets of ~68 eV to ~790 eV from the M-line centers, which is 10 to 100 times the line widths. That tail is not independently established; Auger channels, Fano interference, or a smooth background could change the shape and hence the enhancement by a factor of several. Without any uncertainty on the central numbers, the reader cannot tell how seriously to take the factor-of-ten claim. The 111In estimate (10^4 enhancement) is even more sensitive to the same assumption. Still, this is a feasibility calculation, not a precision measurement, and the qualitative conclusion—that subthreshold states can dominate the endpoint window—is physically reasonable and worth testing.\n\nNo circularity: the parameters come from external compilations, and the 163Ho comparison is a check, not a fit. So the core logic is sound even if the quantitative output is soft.\n\nWho is this for? Experimental groups planning calorimetric EC experiments, and theorists working on atomic effects in neutrino mass searches. It deserves a serious referee: a good referee should ask for a discussion of line-shape uncertainties, a more realistic spectral function, or at least a clear caveat that the quoted factors are ansatz-dependent estimates. I would send it to peer review.","headline":"A transparent, well-caveated estimate that subthreshold bound states can boost the 159Dy endpoint EC rate by an order of magnitude, but the quantitative factors rest on an unvalidated Lorentzian tail.","tokens_in":7895,"tokens_out":2736,"would_cite":false,"duration_ms":26256,"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":"By adding energetically forbidden atomic hole states to the spectral function, this paper predicts that the electron-capture rate of 159Dy near the neutrino endpoint is more than an order of magnitude larger than previously estimated.","keywords":["electron capture","subthreshold atomic states","spectral function","159Dy","neutrino mass","ultra-low Q value electron capture","Lorentzian line shape","endpoint spectrum"],"falsifier":"Measure the 159Dy EC spectrum with a TES microcalorimeter and count events in the neutrino-energy window $m_\\nu<E_\\nu<m_\\nu+1$ eV at the endpoint; if the rate per decay matches the continuum-only prediction ($r\\sim2.9\\times10^{-12}$) rather than the subthreshold-enhanced prediction ($r\\sim3.3\\times10^{-11}$) for $m_\\nu=1$ eV, the Lorentzian-tail contribution is ruled out. Alternatively, a high-statistics fit of the 163Ho spectrum that shows the spectral function falling faster than a Lorentzian at 60-800 eV offsets would falsify the assumed tail shape.","tokens_in":6785,"feed_emoji":"⚛️","tokens_out":10827,"duration_ms":97096,"temperature":0.7,"pith_summary":"This paper argues that the electron-capture (EC) spectrum of 159Dy has been underestimated near the neutrino endpoint because conventional analyses include only atomic hole states whose excitation energy is kinematically allowed. The authors add 'subthreshold' M-shell hole states, whose binding energy exceeds the available endpoint energy, and show that the tails of their Lorentzian line shapes reach into the allowed region. Including them enhances the partial EC rate near zero-momentum neutrino emission by more than an order of magnitude (r rises from $2.9\\times10^{-12}$ to $3.3\\times10^{-11}$ for a 1 eV window at $m_\\nu=1$ eV). Because the endpoint region is where a neutrino mass would show up, this makes 159Dy a more promising ultra-low-Q EC candidate and suggests that experimentally determining $P(E_{ex})$ can resolve endpoint ambiguities.","feed_headline":"Forbidden atomic states multiply 159Dy's neutrino-mass signal tenfold","feed_subtitle":"For 159Dy, this makes an ultra-low Q-value electron-capture candidate viable for neutrino-mass searches.","key_machinery":"The load-bearing object is the atomic spectral function $P_x(E_{ex})$ of Eq. (8), written for each hole state $x=(n,l,j)$ as a Breit-Wigner/Lorentzian line $$P_x(E_{ex}) \\propto \\frac{2 n_x B_x \\$beta_x^{2}$}{4\\pi}\\,\\frac{\\Gamma_x/(2\\pi)}{(E_{ex}-\\epsilon_x)^2+\\$Gamma_x^{2}$/4},$$ with $\\epsilon_x$ the binding energy of the hole, $\\Gamma_x$ its atomic width, $\\beta_x^2/(4\\pi)$ the electron density at the nucleus, and $n_x$, $B_x$ occupation and exchange-overlap factors. The argument works because this shape does not vanish when $E_{ex}<\\epsilon_x$: the M-shell peaks at 1.97 and 1.77 keV extend tails into the allowed region around 1.18 keV, with the M2 center only 68 eV (about 23 widths) from the 1.7 keV endpoint. The rate formula, $d\\lambda_{EC}/dE_\\nu = (G_\\beta^2/2\\pi)\\,p_\\nu E_\\nu\\, C\\, P(E_{ex})$, then multiplies the phase space by this tail-enhanced spectral function, producing the order-of-magnitude endpoint enhancement.","core_discovery":"The central claim is that the residual-energy spectral function $P(E_{ex})$ of electron capture in $^{159}$Dy is dominated near the neutrino endpoint by contributions from energetically forbidden subthreshold atomic hole states, specifically the M-shell holes at 1.768 keV (M2) and 1.968 keV (M1). Because each hole state contributes a Lorentzian tail of width $\\Gamma_x$ centered at its binding energy, these states still populate excitation energies below the endpoint $Q-\\Delta_{\\rm nucl}\\simeq1.18$ keV, even though their centers lie above it. The authors demonstrate that including these states raises the fraction $r$ of decays in the endpoint window $m_\\nu<E_\\nu<m_\\nu+\\delta$ ($\\delta=1$ eV) to $3.3\\times10^{-11}$ at $m_\\nu=1$ eV and $8.4\\times10^{-12}$ at $m_\\nu=0.1$ eV, roughly an order of magnitude above the continuum-only values. They further find that a larger $Q$ value, such as the older 1.7(12) keV estimate, moves the endpoint closer to the M2 resonance and increases the enhancement to factors of $1.13\\times10^{1}$ and $3.02\\times10^{2}$, opposite to the usual expectation that larger $Q$ dilutes the endpoint fraction. The paper concludes that $P(E_{ex})$ can be extracted experimentally, making the endpoint region of $^{159}$Dy accessible for neutrino-mass studies.","pith_inferences":["A natural extension of the paper's logic is that any EC or beta-decay candidate whose endpoint lies within a few hundred eV of a deep atomic hole state should be screened for a similar subthreshold boost; proximity to atomic resonances may matter as much as the bare $Q$ value.","If later atomic-structure calculations or $^{163}$Ho data show the line shape falls faster than a Lorentzian, the enhancement factor for $^{159}$Dy could change by orders of magnitude; a data-driven $P(E_{ex})$ extracted from the measured spectrum would make the neutrino-mass analysis more robust than relying on parameterized tails.","The same 'virtual capture' reasoning might apply to other decay modes and to daughter charge states that shift binding energies, so the size of the boost could be tuned by choosing the initial atomic configuration, a possibility the paper does not explore.","A dedicated $^{159}$Dy source measured with a transition-edge sensor array in the gap between hole states could map $P(E_{ex})$ directly and test whether the tail is the smooth Lorentzian assumed here."],"forward_implications":["The ultra-low-Q electron-capture candidate $^{159}$Dy becomes competitive for neutrino-mass searches: its endpoint-window fraction $r$ reaches $3\\times10^{-11}$, larger than tritium's $1.5\\times10^{-12}$ for the same 1 eV window at $m_\\nu=1$ eV.","Experiments do not need to pin down the $Q$ value precisely: the endpoint fraction depends only weakly on $Q$, and a larger $Q$ value (up to 1.7 keV) actually increases the count rate near the endpoint.","The spectral function $P(E_{ex})$ can be treated as an experimentally determinable parameter, so the same measurement that searches for the neutrino mass can also calibrate the atomic response.","The same subthreshold-tail enhancement applies to other EC candidates: for $^{111}$In, whose effective $Q$ value lies within 0.04 keV of the L2 hole energy, including the subthreshold state enhances the endpoint rate by roughly $10^4$.","The measured $^{163}$Ho EC spectrum can serve as a test of the Lorentzian-tail parametrization, supporting the use of this parametrization for $^{159}$Dy after correcting for phase space."],"supporting_citations":[{"why":"Supplies the 159Dy electron-capture candidate, its effective Q value Q−Δnucl = 1.18(19) keV, and the 5/2− excited state of 159Tb that define the endpoint region.","marker":"[13]"},{"why":"Gives the standard particle-hole form of the atomic spectral function P(Eex), including the Lorentzian line shape of Eq. (8) that the paper extends to subthreshold states.","marker":"[15]"},{"why":"Supplies atomic binding energies used for the hole-state centers in Eq. (8).","marker":"[16]"},{"why":"Supplies the atomic level widths Γx for the M and N shells used in the Lorentzian tails of Table I.","marker":"[17]"},{"why":"Provides the measured 163Ho electron-capture spectrum used to argue that P(Eex) is experimentally accessible and flat near the endpoint.","marker":"[10]"},{"why":"Supplies the older mass-evaluation result Q−Δnucl = 1.7(12) keV, used to show the enhancement grows with larger Q.","marker":"[19, 20]"},{"why":"Offers a more sophisticated ab initio treatment of the EC spectrum used to motivate the caveat that the simple Lorentzian parametrization needs improvement.","marker":"[21]"},{"why":"Provides the 111In Q value whose proximity to the L2 subthreshold hole state supports the extension of the enhancement mechanism to another EC candidate.","marker":"[25]"}],"fun_headline_variants":["Forbidden atomic states boost 159Dy neutrino endpoint by 10x","Subthreshold holes multiply 159Dy's neutrino-mass signal","Off-shell states make 159Dy viable for neutrino-mass hunt","159Dy: forbidden states lift neutrino endpoint fraction tenfold","Energetically forbidden states enhance 159Dy EC rate"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole enhancement rests on the assumption that the atomic line shape remains a Lorentzian tail many widths away from resonance, for example with the M2 hole center about 23 widths from the 1.7 keV endpoint, and that no steeper cutoff or unresolved background changes the tail at those offsets.","fun_headline_variants_meta":{"raw":{"variants":["Forbidden atomic states boost 159Dy neutrino endpoint by 10x","Subthreshold holes multiply 159Dy's neutrino-mass signal","Off-shell states make 159Dy viable for neutrino-mass hunt","159Dy: forbidden states lift neutrino endpoint fraction tenfold","Energetically forbidden states enhance 159Dy EC rate"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000582,"raw_usage":{"total_tokens":2756,"prompt_tokens":975,"completion_tokens":1781,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":591,"completion_tokens_details":{"reasoning_tokens":1695}},"tokens_in":591,"tokens_out":1781,"duration_ms":14671,"temperature":1.0,"reasoning_tokens":1695,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:56:48.842524+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the 159Dy EC spectrum with a TES microcalorimeter and count events in the neutrino-energy window $m_\\nu<E_\\nu<m_\\nu+1$ eV at the endpoint; if the rate per decay matches the continuum-only prediction ($r\\sim2.9\\times10^{-12}$) rather than the subthreshold-enhanced prediction ($r\\sim3.3\\times10^{-11}$) for $m_\\nu=1$ eV, the Lorentzian-tail contribution is ruled out. Alternatively, a high-statistics fit of the 163Ho spectrum that shows the spectral function falling faster than a Lorentzian at 60-800 eV offsets would falsify the assumed tail shape.","supporting_citations":[{"cited_title":"$^{159}$Dy electron-capture: a strong new candidate for neutrino mass determination","cited_arxiv_id":"2106.06626","evidence_quote":"Supplies the 159Dy electron-capture candidate, its effective Q value Q−Δnucl = 1.18(19) keV, and the 5/2− excited state of 159Tb that define the endpoint region."},{"cited_title":"Orbital electron capture by the nucleus,","cited_arxiv_id":null,"evidence_quote":"Gives the standard particle-hole form of the atomic spectral function P(Eex), including the Lorentzian line shape of Eq. (8) that the paper extends to subthreshold states."},{"cited_title":"X-ray data booklet (2009),","cited_arxiv_id":null,"evidence_quote":"Supplies atomic binding energies used for the hole-state centers in Eq. (8)."},{"cited_title":"Widths of the Atomic K – N 7 LEVELS,","cited_arxiv_id":null,"evidence_quote":"Supplies the atomic level widths Γx for the M and N shells used in the Lorentzian tails of Table I."},{"cited_title":"High-resolution and low-background 163Ho spectrum: interpretation of the resonance tails,","cited_arxiv_id":null,"evidence_quote":"Provides the measured 163Ho electron-capture spectrum used to argue that P(Eex) is experimentally accessible and flat near the endpoint."},{"cited_title":"Ab initio calculation of the electron capture spectrum of 163ho: Auger–meitner decay into continuum states,","cited_arxiv_id":null,"evidence_quote":"Offers a more sophisticated ab initio treatment of the EC spectrum used to motivate the caveat that the simple Lorentzian parametrization needs improvement."},{"cited_title":"High-precision electron-capture Q value measurement of 111In for electron-neutrino mass determination,","cited_arxiv_id":null,"evidence_quote":"Provides the 111In Q value whose proximity to the L2 subthreshold hole state supports the extension of the enhancement mechanism to another EC candidate."}],"review_version":2}