Pith. sign in

REVIEW 4 major objections 5 minor 13 references

Using krypton-83m to determine the neutrino-mass bias caused by a non-constant electric potential of the KATRIN source

T0 review · 4 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read KATRIN's neutrino-mass bias from a non-constant source potential is set by shifts between scattered and unscattered spectra, not by Gaussian broadening.

desk verdict The paper's real contribution is reframing the KATRIN source-potential bias as a first-moment problem; the quantitative size of the effect is still open, but the qualitative argument is solid. read the letter →

arxiv 2506.13829 v1 pith:GOK4W376 submitted 2025-06-15 physics.ins-det nucl-exphysics.data-an

Moritz Machatschek This is my paper · ORCID
classification physics.ins-detnucl-exphysics.data-an
keywords KATRINneutrinomasssourcepotentialkrypton-83menergy-lossshiftstarting-potentialdistributionbetaspectrumendpointMAC-Efilter
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper tries to establish that the leading-order neutrino-mass bias from an inhomogeneous electric potential in the KATRIN tritium source comes from first moments of the starting-potential distribution, specifically the mutual shifts between spectra of different electron-scattering multiplicities, rather than from the Gaussian broadening of the potential. It derives a perturbative expansion in the statistical moments of the starting-potential distribution and shows that the allowed parameter space of these shifts is an ellipsoid constrained by a single measured standard deviation. If correct, this means the krypton-83m calibration spectrum, where the once-scattered line sits about 13 eV below the unscattered line, is the right instrument for the dominant source-potential systematic: the distance between the two lines directly measures the shift. The paper estimates that the squared neutrino-mass bias scales roughly as a 1 eV coefficient times the shift, turning a tens-of-millivolt potential asymmetry into a substantial systematic, and projects that measuring the shift will reduce this uncertainty by a factor of two to three. This matters because the experiment's sensitivity goal depends on controlling every spectrum-shaping systematic near the endpoint.

What carries the argument

The central object is the starting-potential distribution $\mathrm{SPD}_i(\nu) = \langle \delta(\nu - V(z)) P_i(z) \rangle$, which, because convolution is commutative and associative, enters the KATRIN response exactly like an additional energy-loss distribution $f'_i = \mathrm{SPD}_i \otimes f_i$. The argument is carried by the first two moments of these distributions, especially the energy-loss shifts $\Delta_{i0}[V] = \mathrm{Cov}_0[(P_i - P_0)/P_0, V]$, the bound $|\Delta_{i0}[V]| \le \kappa_i \sigma_0[V]$, and the ellipsoidal shape-operator space $\vec{\rho}^\mathsf{T} P^{-1} \vec{\rho} = 1$. The paper then linearizes the response-function moments to obtain the susceptibility coefficients $a_i$ and $\varepsilon_i$ that convert the shifts into endpoint and squared-mass biases.

What would settle it

Run the full KATRIN simulation with a known antisymmetric potential proportional to the boundary shape $Q_1(z)$ and fit the generated spectrum with the zero-potential response over the standard analysis range; if the recovered $\Delta m^2$ deviates from $-\varepsilon_1 \Delta_{10}$ by more than the statistical precision, the leading-order coefficient is wrong. A more direct check is to compare the predicted location of the maximum of $a_1$ with the full simulation, since the paper reports a visible mismatch there.

Watch

Extended reading notes

Core claim

The paper's central claim is that, once electron scattering in the KATRIN source is taken into account, the leading-order neutrino-mass bias from an inhomogeneous source potential is not the familiar Gaussian broadening $\sigma_0[V]$ but the energy-loss shifts $\Delta_{i0}[V] = \langle V \rangle_i - \langle V \rangle_0$, which are the first moments of the starting-potential distributions for each scattering multiplicity. Because scattered electrons come preferentially from the rear of the source and unscattered electrons from the front, even a small antisymmetric potential produces such shifts, and the squared neutrino mass is biased approximately by $\Delta m^2 \propto 1\,\mathrm{eV} \cdot \Delta_{10}$. The paper proves these shifts are constrained by the measured standard deviation through $|\Delta_{i0}[V]| \le \kappa_i \sigma_0[V]$, with the full vector of shifts lying on an ellipsoid, and it identifies the $^{83\mathrm{m}}$Kr conversion-electron lines as the observable that directly measures $\Delta_{10}$.

Load-bearing premise

The result depends on the response function depending only on surplus energy $\epsilon = E - qU$ and on evaluating the susceptibility coefficients with a simplified response model and unweighted averaging over the retarding-energy scan; if the real energy-dependent scattering or the actual measurement-time weighting breaks that form, the numerical coefficients, including the claimed $\Delta m^2 \propto 1\,\mathrm{eV} \cdot \Delta_{10}$ scaling, would change, and the paper itself notes that the maximum of $a_1$ does not match the full simulation.

Editorial extensions

If this is right

  • The dominant source-potential systematic can no longer be treated as a Gaussian broadening; future KATRIN analyses must include the first-moment shifts as separate parameters.
  • Measuring $\Delta_{10}$ from the $^{83\mathrm{m}}$Kr spectrum, with the projected few-millivolt uncertainty, can shrink the neutrino-mass systematic by roughly a factor of two to three.
  • The ellipsoidal parameter space provides a practical way to constrain all $\Delta_{i0}$ simultaneously from a single measured $\sigma_0$.
  • The effective single-parameter treatment of the shifts is discouraged because it corrects the mass bias but leaves residual structure in other observables; the full ellipsoid parametrization is recommended instead.
  • A measured $^{83\mathrm{m}}$Kr $\Delta_{10}$ near zero would only weakly constrain the potential, because antisymmetric potentials dominate $\Delta_{10}$ while symmetric ones contribute mostly to $\sigma_0$.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the first-moment logic is correct, other gaseous-source tritium endpoint experiments using MAC-E filters should re-examine their source-potential systematics in terms of scattering-multiplicity shifts, not just variance broadening.
  • The same formalism suggests that calibration lines with higher scattering multiplicities, beyond the once-scattered line, offer additional handles on the potential shape if their widths can be resolved.
  • A direct experimental test would be to compare the fitted neutrino-mass bias in simulated KATRIN data with a known asymmetric potential against the prediction $-\varepsilon_1 \Delta_{10}$; the paper already notes a mismatch in the location of the $a_1$ maximum, making such a test decisive.
  • The distinction between $\sigma_0$ and $\Delta_{10}$ implies that the source potential should be monitored by line-position shifts between scattered and unscattered calibration lines, not only by line broadening.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper analyzes how a non-constant electric potential in the KATRIN tritium source biases the measured neutrino mass and beta-spectrum endpoint. It introduces starting-potential distributions for each scattering multiplicity, defines their moments (mean, standard deviation, and mutual shifts Δ_i0), and derives a covariance bound |Δ_i0| ≤ κ_i σ_0 together with an ellipsoidal parameter space for the shifts (Section 2). In Section 3 it argues that, because of electron scattering, the leading-order neutrino-mass bias is caused by the first-moment shifts Δ_i0 rather than by the Gaussian broadening σ_0, and it derives susceptibility coefficients a_i and ε_i (Eq. 52) that relate ΔE_0 and Δm² to Δ_i0. The paper concludes that Δm² scales roughly as 1 eV × Δ_10, and that the krypton-83m measurement of Δ_10 can constrain this systematic.

Significance. If the central qualitative claim holds, the paper corrects a potentially important systematic treatment for KATRIN: the dominant source-potential bias would come from shifts between spectra of different scattering multiplicities, not from a simple Gaussian broadening. The covariance and ellipsoid mathematics in Section 2.5 are clean and appear correct: the bound |Δ_i0| ≤ κ_i σ_0 follows directly from Cauchy-Schwarz, and the ellipsoid parameterization is a sound geometric result. The distinction between shift-type and broadening-type observables is conceptually valuable and likely to influence how source-potential systematics are parametrized and constrained by krypton-83m calibration data. However, the quantitative susceptibility coefficients in Eq. (52) are not yet validated: they rest on a simplified response model and unweighted ΔqU averaging, and Fig. 9 reports a visible mismatch between a_1 and the full Asimov simulation of reference [6]. The paper itself repeatedly describes the calculation as qualitative, so the strength of the contribution lies in the qualitative hierarchy and the parameter-space geometry, not in the specific 1 eV coefficient.

major comments (4)
  1. [§3.3, Eq. (52)] The central quantitative claim, Δm² ≈ −ε_1 Δ_10 with ε_1 on the 1 eV scale, rests entirely on the coefficients a_i and ε_i derived in Eq. (52). The derivation is only sketched ("only the basic ideas are described"), uses a simplified response model, and assumes an unweighted ΔqU average. Figure 9 shows that the maximum of a_1 deviates visibly from the full Asimov simulation in [6]; the authors attribute this to the approximations. As presented, ε_i is an unvalidated estimate, not an established result. Please either derive Eq. (52) from stated assumptions in an appendix and validate both a_1 and ε_1 against the full simulation, or explicitly downgrade the 1 eV scaling to a motivated heuristic and remove it from the abstract and conclusions.
  2. [§3.2, Eq. (45) and footnote 2] Equations (45) define the bias as an average over all retarding potentials U with a weight that must contain the measurement-time distribution, and the footnote adds that a χ²-minimization would modify the result. The evaluation in Fig. 9 instead uses an unweighted ΔqU average. This is a different weighting from the actual KATRIN analysis, so the numerical values of a_i and ε_i are not tied to a real measurement procedure. Please state the relationship between the unweighted average and the experimental analysis weighting, and estimate the resulting uncertainty on the coefficients, or restrict the conclusions to the unweighted model.
  3. [§3.3, Normalization] The normalization step is calibrated with the constant-potential case ΔE_0 = q⟨V⟩_0 and Δm² = 0, adding a correction term proportional to the endpoint. This is an ad-hoc convention that enters the derivation of Eq. (52). It is not shown that the final coefficients are independent of this calibration or that the 1 eV scale is robust under a different normalization choice. Please make the calibration dependence explicit and demonstrate that the quoted scale is not an artifact of this convention.
  4. [§2.3, suppression of higher moments] The perturbative truncation at second order relies on the statement that moments of order n ≥ 3 are suppressed by (σ_i/w)^n, justified only by a reference to Hölder's inequality without showing the argument. Since this truncation underlies the entire moment-based expansion, please provide the derivation or a precise citation; otherwise the stated suppression factor is unsubstantiated.
minor comments (5)
  1. [Throughout] There are several typographical errors, including "througout" (abstract or Section 2) and "previous believe" in the conclusions; these should be corrected.
  2. [Fig. 5 caption] The caption refers to an "N123 spectrum" and lines "N1", "N2", "N3" without defining these labels in the text; please add a brief definition.
  3. [§3.3, Fig. 9] The figure caption says "good qualitative agreement" for ε_1 but "the position of its maximum deviates visibly" for a_1; please quantify the deviations (e.g., shift in eV and relative amplitude) so the reader can judge the level of agreement.
  4. [§2.3, Eq. (31)] The definition of the shape operators ρ̂_i uses σ_0[V] both in the numerator denominator; consider clarifying that σ_0 is the standard deviation of the unscattered electron distribution, as stated around Eq. (25).
  5. [§4.3, Eq. (67)] The propagation of the krypton Δ_10 constraint to tritium via Eq. (67) would benefit from a brief derivation or a reference to the covariance calculation, as the expression with the ± sign is not immediately obvious.

Circularity Check

0 steps flagged · score 2.0 of 10

No circular reduction; the moment-based derivation is self-contained, with [6] cited only as a comparison benchmark.

full rationale

I find no circular step in the derivation chain. The central objects are the moments of the starting-potential distributions (Eqs. 20, 24, 25), and the claim that the leading-order neutrino-mass bias is linear in the shifts Delta_i0 rather than quadratic in sigma0 follows from the Taylor expansion (Eq. 48) and the explicit first-order result (Eq. 51). No fitted constant is renamed as a prediction. The bound |Delta_i0| <= kappa_i sigma0 (Eq. 32) and the ellipsoid parametrization (Eq. 41) are derived in-text from covariance identities using Cauchy-Schwarz; they are not imported from prior work. The susceptibility coefficients a_i and epsilon_i (Eq. 52) are analytic expansions obtained under stated approximations, and are not fit to the bias they are used to predict. Reference [11] supplies only the independently measured sigma0, while reference [6], the author's own thesis, is used as a simulation benchmark for qualitative comparison in Fig. 9 and is not an input to the central equations. The paper explicitly labels the analytic results as 'only qualitative' and reports that the maximum of a1 does not match the full simulation, which is a validation gap rather than a circularity. The only mild flag is the presence of self-citation [6], but it is not load-bearing, so the score is 2 rather than 0.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

The ledger is dominated by domain assumptions about the KATRIN response model rather than fitted constants. The central formulas are derived from those assumptions, with measured sigma0 entering from [11] and simulation comparisons from [6]. No new physical entities are introduced.

free parameters (3)
  • Simplified response model = not specified
    Used in Section 3.3 and Figure 9 to evaluate a_i and epsilon_i; the paper does not specify its form, and one coefficient's maximum disagrees with full simulation [6].
  • Unweighted Delta-qU averaging = uniform
    The coefficients in Figure 9 assume unweighted averaging over retarding potentials, whereas the actual bias depends on the measurement-time distribution (Section 3.2).
  • Analysis range = 40 eV
    The susceptibility values and conclusions depend on the chosen 40 eV analysis window; the paper plots coefficients versus this range, so it functions as a parameter of the calculation.
assumptions (6)
  • domain assumption The starting-potential effect is fully described by starting-potential distributions SPD_i; different z-dependent potentials with equal SPDs are indistinguishable in measured electron spectra.
    Section 2.2, Eqs. 11-12. This is the modeling foundation that reduces the potential to its moments.
  • domain assumption Energy and angle dependence of the starting-potential distributions can be neglected, SPD_i(nu) approximately equals <delta(nu-V) P_i(z)>.
    Eq. 13 in Section 2.2; the paper estimates the angle effect via Eqs. 58-60 and argues it is subleading, but it remains an approximation.
  • domain assumption The response function depends only on surplus energy epsilon = E - qU.
    Section 3.3, 'Energy dependency'. Needed to reduce the response integrals to functions of Delta-qU and derive Eq. 52; not proven against energy-dependent scattering cross sections.
  • domain assumption A second-order Taylor expansion of the measured spectrum around the potential moments is sufficient; orders n >= 3 are suppressed by (sigma_i/w)^n.
    Section 2.3, Hoelder argument; uses sigma0 around 30 mV and w around 1 eV. The paper recommends further studies to confirm the suppression.
  • domain assumption Scattering probabilities follow a Poisson model with effective column density N_eff per Eq. 3-4.
    Section 2.1, based on [5]; this is the standard KATRIN response model adopted without re-derivation.
  • ad hoc to paper The normalization correction for the response-moment calculation is calibrated with the constant-potential case Delta-E0 = q <V>0 and Delta-m2 = 0.
    Section 3.3, 'Normalization'; this calibration chooses the boundary condition for the endpoint correction and is not derived from first principles.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Using krypton-83m to determine the neutrino-mass bias caused by a non-constant electric potential of the KATRIN source." pith.science (2026). https://pith.science/paper/GOK4W376

@misc{pith2026250613829,
  author       = {Pith},
  title        = {Pith review of: Using krypton-83m to determine the neutrino-mass bias caused by a non-constant electric potential of the KATRIN source},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/GOK4W376}},
  note         = {Machine review of arXiv:2506.13829}
}
abstract

Precision spectroscopy of the electron spectrum of the tritium $\beta$ decay near the kinematic endpoint is a direct method to determine the effective electron antineutrino mass. The KArlsruhe TRItium Neutrino (KATRIN) experiment aims to determine this quantity with a sensitivity of better than $0.3\,\mathrm{eV}$ ($90\,\%$~C.L.). An inhomogeneous electric potential in the tritium source of KATRIN leads to a distortion of the $\beta$ spectrum, which directly impacts the neutrino-mass observable. This effect can be quantified through precision spectroscopy of the conversion electrons of co-circulated metastable $^{83\mathrm{m}}Kr$. This work perturbatively describes the effect of the source potential on the recorded spectra, and thus establishes the leading-order observables of this effect.

Figures

Figures reproduced from arXiv: 2506.13829 by the authors.

Figure 1
Figure 1. Schematic of the KATRIN beamline with FPD. [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Longitudinal density profile in the source. [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 4
Figure 4. Transmission function. This plot illustrates equa￾tion 11. The transmission condition is a step function in the surplus energy, corresponding to the energy required for electrons at a given angle to be transmitted. It is modified by the potential distribution, here modeled as a Gaussian, which may also depend on the angle due to the weighting with the signal electron distributions. For unscattered elec￾trons, the tr… view at source ↗
Figures from the paper (4 more)
Figure 5
Figure 5. Figure 5: Krypton-83m N23 spectrum with energy loss. The orange line shows a fit of the measured N123 spectrum. One￾time scattered electrons are visible as a rate increase approx￾imately 13 eV below the main peak (purple). Changes of this distance ∆10[V] are an observable of the…
Figure 7
Figure 7. Figure 7: Parameter spaces. The left plot shows the high correlation of the ⃗ρ parameter space: the major axis of the ellipsoid is oriented almost perfectly along the diagonal, such that high (low) values of one ρˆi correspond to high (low) values of the other ρˆj . In the right…
Figure 8
Figure 8. Figure 8: Integrations over the tritium spectrum. In the con￾volution with an energy fluctuation, the spectrum is aver￾aged around the selected energy E, weighted with the fluc￾tuation. Since the width of the fluctuation is small, the cut by the β spectrum is only relevant for e…
Figure 10
Figure 10. Figure 10: Dependency of the unscattered electron distribu [PITH_FULL_IMAGE:figures/full_fig_p011_10.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

13 extracted references · 7 canonical work pages

  1. [6]

    M. B. Machatschek. A Phenomenological Theory of KATRIN Source Potential Systematics and its Applica- tion in Krypton-83m Calibration Measurements . PhD thesis, KIT, 2021. URL https://doi.org/10.5445/ IR/1000132391

    arXiv 2021

  2. [1]

    Fukuda et al

    Y . Fukuda et al. Phys. Rev. Lett., 82:2644–2648, 1999. doi:10.1103/PhysRevLett.82.2644

    doi:10.1103/physrevlett.82.2644 1999

  3. [2]

    Q. R. Ahmad et al. Phys. Rev. Lett., 89(1):011301, 2002. doi:10.1103/PhysRevLett.89.011301

    doi:10.1103/physrevlett.89.011301 2002

  4. [3]

    Science, 388(6743):180–185,

    KATRIN Collaboration. Science, 388(6743):180–185,

    work page

  5. [4]

    Aker et al

    M. Aker et al. J. Instrum. , 16(08):T08015, 2021. doi:10.1088/1748-0221/16/08/T08015

    doi:10.1088/1748-0221/16/08/t08015 2021

  6. [5]

    Kleesiek et al

    M. Kleesiek et al. Eur. Phys. J. C , 79(3):204, 2019. doi:10.1140/epjc/s10052-019-6686-7

    work page doi:10.1140/epjc/s10052-019-6686-7 2019

  7. [7]

    Kalempa and F

    D. Kalempa and F. Sharipov. Separation phenomenon in the windowless gaseous tritium source of katrin exper- iment. ternary mixtures, 2010. URL http://fisica. ufpr.br/sharipov/KATRIN/report5.pdf. Ac- cessed 17 April 2025

    work page 2010

  8. [8]

    M. Hötzel. Simulation and analysis of source-related effects for KATRIN . PhD thesis, KIT, 2012. URL https://doi.org/10.5445/IR/1000031259

    work page arXiv 2012

Show all 13 references
  1. [9]

    S. Groh. Modeling of the response function and mea- surement of transmission properties of the KATRIN ex- periment. PhD thesis, KIT, 2015. URL https://doi. org/10.5445/IR/1000046546

    2015

  2. [10]

    M. Slezák. Monitoring of the energy scale in the KA- TRIN neutrino experiment . PhD thesis, Prague, Inst. Phys., 2016. URL https://cds.cern.ch/record/ 2243754

    2016

  3. [11]

    Acharya et al

    H. Acharya et al. arXiv:2503.13221, 2025. doi:https://doi.org/10.48550/arXiv.2503.13221

    2025 doi

  4. [12]

    R. G. H. Robertson and D. A. Knapp. Annu. Rev. Nucl. Part. Sci. , 38:185–215, 1988. doi:10.1146/annurev.ns.38.120188.001153

    1988

  5. [2025]

    doi:10.1126/science.adq9592

    doi

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.