pith. machine review for the scientific record. sign in

arxiv: 2605.00554 · v1 · submitted 2026-05-01 · ⚛️ physics.comp-ph · nucl-th· physics.atom-ph

Recognition: unknown

MuDirac 1.3.0: A Sustainable Software Tool for Calculating Ground State Nuclear Properties Using Muonic X-Ray Measurements

Adrian Hillier, Albert Bartok, Leandro Liborio, Martin Plummer, Milan Kumar, Philip Jones, Subindev Devadasan

Pith reviewed 2026-05-09 15:09 UTC · model grok-4.3

classification ⚛️ physics.comp-ph nucl-thphysics.atom-ph
keywords muonic X-raynuclear charge radiusFermi distributionsoftware toolnuclear propertiesmuon physicsatomic physicscomputational physics
0
0 comments X

The pith

MuDirac 1.3.0 provides an open software tool to calculate nuclear charge radii from muonic X-ray measurements using a two-parameter Fermi distribution.

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

The paper introduces MuDirac version 1.3.0 as a publicly available software package for the negative muon community. It processes experimental muonic X-ray transition energies to derive ground-state nuclear properties such as the nuclear charge radius. The calculations rely on modeling the nuclear charge density with a two-parameter Fermi distribution to achieve both accuracy and computational speed. This setup removes the need for users to build custom simulation code from scratch.

Core claim

With MuDirac (1.3.0), the community will be able to accurately and efficiently estimate nuclear properties, such as the nuclear charge radius, by assuming a 2-parameter Fermi distribution of the nuclear charge. The software is open, publicly available, sustainable and computationally efficient for calculating ground state nuclear properties using muonic X-ray measurements.

What carries the argument

The two-parameter Fermi distribution model of nuclear charge density, which converts measured muonic X-ray transition energies into extracted values for nuclear radii and related ground-state properties.

If this is right

  • Researchers gain a ready method to extract nuclear charge radii directly from muonic X-ray data without writing new simulation code.
  • The tool supports routine, low-cost computations that fit within standard laboratory workflows.
  • Public release of the code enables community inspection, reuse, and incremental improvement of the fitting procedures.
  • Ground-state nuclear properties for nuclei accessible to negative muons become easier to obtain and compare across experiments.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Adoption of the tool could encourage more uniform reporting standards for nuclear radii derived from muonic data.
  • Later versions might add options for alternative charge-density models to handle nuclei where the Fermi form is known to be inadequate.
  • Wider use may accelerate cross-checks between muonic results and data from other probes such as laser spectroscopy of muonic atoms.

Load-bearing premise

Modeling the nuclear charge density with a two-parameter Fermi distribution is sufficient to produce accurate nuclear radii from the muonic X-ray transition energies.

What would settle it

A direct comparison in which nuclear charge radii obtained from MuDirac on a set of nuclei differ markedly from independent determinations of the same radii by electron scattering or other established techniques.

Figures

Figures reproduced from arXiv: 2605.00554 by Adrian Hillier, Albert Bartok, Leandro Liborio, Martin Plummer, Milan Kumar, Philip Jones, Subindev Devadasan.

Figure 1
Figure 1. Figure 1: Charge density in the spherically symmetric Uniform model and 2pF model view at source ↗
Figure 2
Figure 2. Figure 2: 2pF nuclear model parameter bands which reproduce valid muonic X-ray view at source ↗
Figure 3
Figure 3. Figure 3: 2pF nuclear model parameter bands, which reproduce valid muonic X view at source ↗
Figure 4
Figure 4. Figure 4: χ 2 contours for 206Pb in polar coordinates. χ 2 is calculated using the experimental energies and the MuDirac simulated energies, for the same 4 transitions, over a polar coordinate grid with a resolution of 4 decimal places. The contours mark the χ 2 minimiser and lines enclosing regions of χ 2 min + n for each nσ level of uncertainty. The χ 2 min = 0.48 is found at Rrms = 5.487, θ = 0.324. efficiency fo… view at source ↗
Figure 5
Figure 5. Figure 5: Optimisation accuracy of different optimisation algorithms in view at source ↗
Figure 6
Figure 6. Figure 6: Optimisation runtime of different optimisation algorithms in view at source ↗
Figure 7
Figure 7. Figure 7: Optimisation accuracy and runtime of MuDirac using different optimisation algorithms and coordinate systems using 24 muonic atom datasets. The axes of the plots follow figs. 5 and 6 using equations (23) and (24). The accuracy plot (left) shows the true χ 2 min is found most frequently by the global optimiser as expected, shortly followed by the lm optimisers. The runtime plot (right) shows that the lm opti… view at source ↗
Figure 8
Figure 8. Figure 8: χ 2 (c, t) for a subset of 24 muonic atoms [15]. The χ 2 (c, t) is calculated across 4 transitions for each atom using the default 2pF values (blue) and fitted 2pF values (orange) in MuDirac. The horizontal red dashed line at χ 2 (c, t) = 4 marks the χ 2 (c, t)max where the simulated energies can still be within the uncertainty of the measured energies. For clarity of the plot, data from separate sources o… view at source ↗
Figure 9
Figure 9. Figure 9: The 70JP1 label at the top refers to the reference for the experimental view at source ↗
Figure 10
Figure 10. Figure 10: The 69LB1 label at the top refers to the reference for the experimental view at source ↗
Figure 11
Figure 11. Figure 11: The 66AB1 label at the top refers to the reference for the experimental view at source ↗
Figure 12
Figure 12. Figure 12: The 66AB1, 69AH2 and 68PO1 labels at the top refer to the reference view at source ↗
Figure 13
Figure 13. Figure 13: 2pF nuclear model parameter bands, which reproduce valid muonic X-ray view at source ↗
Figure 14
Figure 14. Figure 14: 2pF nuclear model parameter bands, which reproduce valid muonic X-ray view at source ↗
Figure 15
Figure 15. Figure 15: 2pF nuclear model parameter bands, which reproduce valid muonic X-ray view at source ↗
Figure 16
Figure 16. Figure 16: 2pF nuclear model parameter bands, which reproduce valid muonic X-ray view at source ↗
read the original abstract

The nuclear charge radius is one of the most fundamental quantities of the atomic nucleus. It can be deduced from a combination of experimental measurements of muonicX-raytransitionenergieswithmodellingofthoseX-raytransitionenergies. In thisworkwepresentMuDirac (1.3.0), whichisanopen, publiclyavailable, sustainable and computationally efficient software tool that will be at put the disposal of the negative muon community. With MuDirac (1.3.0), the community will be able to accurately and efficiently estimate nuclear properties, such as the nuclear charge radius, by assuming a 2-parameter Fermi distribution of the nuclear charge.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript presents MuDirac 1.3.0, an open-source software tool for computing ground-state nuclear properties (primarily the nuclear charge radius) from measured muonic X-ray transition energies. It solves the Dirac equation for a two-parameter Fermi nuclear charge distribution and is positioned as a sustainable, publicly available, and computationally efficient resource for the muonic-atom community.

Significance. If the accuracy claim is substantiated, the work would deliver a reusable, open-source implementation of a standard forward-modeling approach, lowering barriers for extracting nuclear radii from muonic data while addressing software sustainability. The emphasis on public availability and efficiency constitutes a concrete contribution to the computational-physics toolkit in this niche.

major comments (2)
  1. [Abstract] Abstract: the repeated claim that the tool enables 'accurate' estimation of nuclear charge radii (and other properties) is unsupported by any quantitative validation. No comparison of extracted radii to independent electron-scattering values, no benchmark against alternative density models (3-parameter Fermi or otherwise), and no error-propagation analysis for the two free Fermi parameters appear in the manuscript. This directly underpins the central utility assertion.
  2. [Method/Implementation] Method/Implementation sections: the forward-modeling procedure is standard, yet the manuscript supplies no test cases, tables, or figures demonstrating numerical convergence, sensitivity to the Fermi radius and diffuseness parameters, or reproduction of published muonic radii for even a single nucleus.
minor comments (2)
  1. [Abstract] Abstract contains typographical and grammatical errors (e.g., 'muonicX-raytransitionenergieswithmodellingofthoseX-raytransitionenergies' and 'will be at put the disposal').
  2. [Method] Notation for the two Fermi parameters is introduced without an explicit equation or table defining their symbols and units.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments, which highlight opportunities to strengthen the manuscript's claims about the tool's utility. We address each major point below and will incorporate revisions to provide the requested validation and demonstrations.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the repeated claim that the tool enables 'accurate' estimation of nuclear charge radii (and other properties) is unsupported by any quantitative validation. No comparison of extracted radii to independent electron-scattering values, no benchmark against alternative density models (3-parameter Fermi or otherwise), and no error-propagation analysis for the two free Fermi parameters appear in the manuscript. This directly underpins the central utility assertion.

    Authors: We agree that the abstract's assertion of 'accurate' estimation requires substantiation to fully support the tool's claimed utility. The two-parameter Fermi model is a standard choice in muonic-atom literature, and MuDirac implements the established Dirac-equation solution for it, but the current manuscript does not include explicit benchmarks. In the revised version we will add a dedicated validation subsection that (i) compares extracted charge radii for several nuclei against independent electron-scattering values, (ii) briefly contrasts results with a 3-parameter Fermi model where literature data exist, and (iii) includes a short discussion of parameter sensitivity and basic error propagation for the two Fermi parameters. We will also moderate the abstract wording to reflect the added evidence. revision: yes

  2. Referee: [Method/Implementation] Method/Implementation sections: the forward-modeling procedure is standard, yet the manuscript supplies no test cases, tables, or figures demonstrating numerical convergence, sensitivity to the Fermi radius and diffuseness parameters, or reproduction of published muonic radii for even a single nucleus.

    Authors: We accept that the manuscript would be improved by explicit demonstrations of the implementation's numerical behavior. In the revised manuscript we will insert a new 'Validation and Test Cases' subsection containing: (i) tables or figures showing convergence of transition energies with respect to radial grid size and integration tolerances, (ii) sensitivity plots or tables illustrating how extracted radii vary with the Fermi radius and diffuseness parameters, and (iii) reproduction of published muonic X-ray derived radii for at least two nuclei, with direct comparison to literature values. These additions will be placed in the Methods/Implementation section or as a supplementary results subsection. revision: yes

Circularity Check

0 steps flagged

No circularity: standard forward-modeling tool with external data dependence

full rationale

The manuscript describes MuDirac as an open software implementation of the Dirac equation solved for a 2-parameter Fermi nuclear charge density to match measured muonic X-ray transition energies and thereby extract the nuclear charge radius. No derivation chain, uniqueness theorem, or fitted parameter is presented that reduces by construction to the paper's own inputs or prior self-citations. The central procedure is a conventional numerical forward model whose outputs are determined by external experimental spectra rather than internal redefinition or ansatz smuggling. The accuracy qualifier is an empirical claim (not demonstrated here) but does not create circularity in the reported method.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central functionality rests on the standard Dirac-equation treatment of a muon in a Fermi-distributed nuclear potential; the two Fermi parameters are the only adjustable quantities introduced by the modeling choice.

free parameters (2)
  • Fermi radius parameter
    One of the two free parameters in the nuclear charge distribution model that is adjusted to match observed transition energies.
  • Fermi diffuseness parameter
    Second free parameter controlling the surface thickness of the nuclear charge distribution.
axioms (1)
  • domain assumption The muonic atom can be accurately described by solving the Dirac equation with a static, spherically symmetric nuclear charge distribution given by the two-parameter Fermi form.
    Standard modeling assumption in muonic-atom physics invoked to connect measured X-ray energies to nuclear radius.

pith-pipeline@v0.9.0 · 5437 in / 1387 out tokens · 38105 ms · 2026-05-09T15:09:51.991522+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

31 extracted references · 22 canonical work pages

  1. [1]

    High-Energy Electron Scattering from Lithium-6 , author =. Phys. Rev. , volume =. 1958 , month =. doi:10.1103/PhysRev.112.1282 , url =

    work page doi:10.1103/physrev.112.1282 1958

  2. [2]

    Precise Measurement of the Muonic X Rays in the Lead Isotopes , author =. Phys. Rev. , volume =. 1969 , month =. doi:10.1103/PhysRev.187.1565 , url =

    work page doi:10.1103/physrev.187.1565 1969

  3. [3]

    Hyperfine Effects in Muonic Bi and Pb , author =. Phys. Rev. , volume =. 1968 , month =. doi:10.1103/PhysRev.169.1 , url =

    work page doi:10.1103/physrev.169.1 1968

  4. [4]

    Lee, W. Y. and Bernow, S. and Chen, M. Y. and Cheng, S. C. and Hitlin, D. and Kast, J. W. and Macagno, E. R. and Rushton, A. M. and Wu, C. S. and Budick, B. , journal =. Finite Distributions of Nuclear M1 and E2 Moments in Muonic ^. 1969 , month =. doi:10.1103/PhysRevLett.23.648 , url =

    work page doi:10.1103/physrevlett.23.648 1969

  5. [5]

    and Powers, R

    Jenkins, David A. and Powers, R. J. and Kunselman, A. R. Muonic Interactions in Zn-66 and Zn-68. Phys. Rev. C. 1970. doi:10.1103/PhysRevC.2.458

    work page doi:10.1103/physrevc.2.458 1970

  6. [6]

    Nucleic Acids Research , author =

    The. Nucleic Acids Research , author =. 2022 , pages =. doi:10.1093/nar/gkac247 , abstract =

    work page doi:10.1093/nar/gkac247 2022

  7. [8]

    Acker and G

    H.L. Acker and G. Backenstoss and C. Daum and J.C. Sens and S.A. Measurement and analysis of muonic X-ray spectra in spherical nuclei , journal =. 1966 , issn =. doi:https://doi.org/10.1016/0029-5582(66)90355-5 , url =

    work page doi:10.1016/0029-5582(66)90355-5 1966

  8. [9]

    , title =

    Marquardt, Donald W. , title =. Journal of the Society for Industrial and Applied Mathematics , volume =. 1963 , doi =

    1963

  9. [10]

    2023 , month =

    Agarwal, Sameer and Mierle, Keir and The Ceres Solver Team , title =. 2023 , month =

    2023

  10. [11]

    and Haberland, Matt and Reddy, Tyler and Cournapeau, David and Burovski, Evgeni and Peterson, Pearu and Weckesser, Warren and Bright, Jonathan and

    Virtanen, Pauli and Gommers, Ralf and Oliphant, Travis E. and Haberland, Matt and Reddy, Tyler and Cournapeau, David and Burovski, Evgeni and Peterson, Pearu and Weckesser, Warren and Bright, Jonathan and. Nature Methods , year =

  11. [12]

    Benchmarking optimization software with performance profiles

    Dolan, Elizabeth D and Mor \'e , Jorge J. Benchmarking optimization software with performance profiles. Math. Program

  12. [13]

    Verbeek, Benjamin , title =

  13. [14]

    Muon Nuclear Data , doi=

    Megumi Niikura and Shinichiro Abe and Shoichiro Kawase and Teiichiro Matsuzaki and Futoshi Minato and Rurie Mizuno and Yukinobu Watanabe and Yuji Yamaguchi , booktitle =. Muon Nuclear Data , doi=

  14. [15]

    2024 , eprint=

    Muon Nuclear Data , author=. 2024 , eprint=

    2024

  15. [16]

    K. A. Beyer and T. E. Cocolios and C. Costache and M. Deseyn and P. Demol and A. Doinaki and O. Eizenberg and M. Gorshteyn and M. Heines and A. Herzáň and P. Indelicato and K. Kirch and A. Knecht and R. Lica and V. Matousek and E. A. Maugeri and B. Ohayon and N. S. Oreshkina and W. W. M. M. Phyo and R. Pohl and S. Rathi and W. Ryssens and A. Turturica and...

    work page arXiv

  16. [17]

    Mizuno and M

    R. Mizuno and M. Niikura and T.Y. Saito and T. Matsuzaki and H. Sakurai and A. Amato and S. Asari and S. Biswas and I. Chiu and L. Gerchow and Z. Guguchia and G. Janka and K. Ninomiya and N. Ritjoho and A. Sato and K. Development of wide range photon detection system for muonic X-ray spectroscopy , journal =. 2024 , issn =. doi:https://doi.org/10.1016/j.n...

    work page doi:10.1016/j.nima.2023.169029 2024

  17. [18]

    Muonic x-ray measurement for the nuclear charge distribution: The case of stable palladium isotopes , author =. Phys. Rev. C , volume =. 2025 , month =. doi:10.1103/PhysRevC.111.034313 , url =

    work page doi:10.1103/physrevc.111.034313 2025

  18. [19]

    The MONUMENT experiment: ordinary muon capture studies for decay , author =. Eur. Phys. J. C , volume =. 2024 , publisher =. doi:https://doi.org/10.1140/epjc/s10052-024-13470-6 , url =

    work page doi:10.1140/epjc/s10052-024-13470-6 2024

  19. [20]

    Theoretical prediction of the fine and hyperfine structure of heavy muonic atoms , author =. Phys. Rev. A , volume =. 2017 , month =. doi:10.1103/PhysRevA.96.032510 , url =

    work page doi:10.1103/physreva.96.032510 2017

  20. [21]

    , TITLE =

    Cataldo, Matteo and Clemenza, Massimiliano and Ishida, Katsuiko and Hillier, Adrian D. , TITLE =. Applied Sciences , VOLUME =. 2022 , NUMBER =

    2022

  21. [22]

    Dutt and V

    Nuclear Ground State Charge Radii from Electromagnetic Interactions , journal =. 1995 , issn =. doi:https://doi.org/10.1006/adnd.1995.1007 , url =

    work page doi:10.1006/adnd.1995.1007 1995

  22. [23]

    arXiv: Nuclear Experiment , year=

    Electronic catalogue of muonic X-rays , author=. arXiv: Nuclear Experiment , year=

  23. [24]

    X-Ray Spectrometry , volume =

    Sturniolo, Simone and Hillier, Adrian , title =. X-Ray Spectrometry , volume =. doi:https://doi.org/10.1002/xrs.3212 , url =. https://analyticalsciencejournals.onlinelibrary.wiley.com/doi/pdf/10.1002/xrs.3212 , year =

    work page doi:10.1002/xrs.3212

  24. [25]

    Angeli and K.P

    Table of experimental nuclear ground state charge radii: An update , journal =. 2013 , issn =. doi:https://doi.org/10.1016/j.adt.2011.12.006 , url =

    work page doi:10.1016/j.adt.2011.12.006 2013

  25. [26]

    2017 , issn =

    Prospects for electron scattering on unstable, exotic nuclei , journal =. 2017 , issn =. doi:https://doi.org/10.1016/j.ppnp.2017.04.002 , url =

    work page doi:10.1016/j.ppnp.2017.04.002 2017

  26. [27]

    1974 , note =

    Charge-distribution parameters, isotope shifts, isomer shifts, and magnetic hyperfine constants from muonic atoms , journal =. 1974 , note =. doi:https://doi.org/10.1016/S0092-640X(74)80003-3 , url =

    work page doi:10.1016/s0092-640x(74)80003-3 1974

  27. [28]

    Mu-Mesonic X-Rays and the Shape of the Nuclear Charge Distribution , author =. Phys. Rev. , volume =. 1954 , month =. doi:10.1103/PhysRev.94.1617 , url =

    work page doi:10.1103/physrev.94.1617 1954

  28. [29]

    1974 , month =

    R C Barrett , title =. 1974 , month =. doi:10.1088/0034-4885/37/1/001 , url =

    work page doi:10.1088/0034-4885/37/1/001 1974

  29. [30]

    1970 , issn =

    Model-independent parameters of the nuclear charge distribution from muonic X-rays , journal =. 1970 , issn =. doi:https://doi.org/10.1016/0370-2693(70)90611-8 , url =

    work page doi:10.1016/0370-2693(70)90611-8 1970

  30. [31]

    Comparison of muonic-atom and electronic x-ray isotope shifts in lead , author =. Phys. Rev. C , volume =. 1974 , month =. doi:10.1103/PhysRevC.9.2444 , url =

    work page doi:10.1103/physrevc.9.2444 1974

  31. [32]

    2025 , issn =

    Critical evaluation of reference charge radii and applications in mirror nuclei , journal =. 2025 , issn =. doi:https://doi.org/10.1016/j.adt.2025.101732 , url =

    work page doi:10.1016/j.adt.2025.101732 2025