REVIEW 3 major objections 4 minor 34 references
Depth profiling the elemental composition with negative muons: Monte Carlo based tools for improved data analysis
T0 review · 3 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read This paper claims that coupling muonic X-ray depth profiles with Monte Carlo stopping simulations determines the thickness of a surface gold layer on copper alloys to about ±1 µm, in agreement with SEM measurements.
desk verdict A useful proof-of-concept for muonic depth profiling of gold layers, but the EM2 discrepancy and tuned parameters mean the ±1 µm claim needs a stronger uncertainty treatment. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the momentum-scan depth profile: a series of muonic X-ray spectra taken at increasing beam momenta, in which each element's fitted peak area is plotted against momentum. The comparison engine is the reduced chi-square between measured and simulated profiles, computed over all momentum runs including zero-intensity runs. The simulation chain uses GEANT4/ARBY, a detailed-geometry Monte Carlo transport code, and SRIM/TRIM, a layered stopping-and-range code run with a hydrogen ion of one-ninth the muon mass, to output the number of stopped muons per layer; the paper assumes this number is proportional to the measured X-ray intensity. The muon itself is the physical enabler: at 207 times the electron mass, its cascade X-rays reach hundreds of keV to MeV, escaping from deep inside metals and giving the technique its non-destructive penetration.
What would settle it
Measure a multilayer standard of known composition (for example gold on nickel on copper) with a precisely machined gold thickness, run the full momentum scan, and fit it. If the best-fit simulated thickness disagrees with the known thickness by more than the claimed ±1 µm, or if switching to a substrate with very different muon capture properties shifts the fitted thickness while SEM confirms a constant gold thickness, the proportionality assumption is wrong. A more targeted check is to measure the actual beam momentum spread at the sample position; if it is 4% while simulations require 5%, the method is compensating for an unmodelled effect rather than measuring thickness cleanly.
Extended reading notes
Core claim
The paper's central claim is that a layer's thickness can be recovered by matching a simulated muon-stopping depth profile to a measured one. Because a negative muon stopped in a layer produces that layer's characteristic muonic X-rays, the authors assume the full-energy peak intensity from a layer is proportional to the number of muons stopped there. They scan the muon beam momentum so that muons stop at increasing depths, build momentum-dependent profiles for gold, nickel, copper, aluminium, and nitrogen, then vary the simulated gold thickness until the reduced chi-square between simulated and measured profiles is minimized. For electroplated samples, the best-fit gold thicknesses are 3.5±1, 4.5±1, and 7.5±1 µm against SEM values of 3.3±0.2, 4.6±0.6, and 7.3±0.8 µm; for the amalgam-gilded brass sample the best fit is 11±1 µm versus 11±1 µm by SEM, and for the bronze sample the best fit is 5±1 µm with a 20% reduced gold density to mimic surface unevenness and air pockets. Both simulation codes agree, with GEANT4/ARBY handling thin layers better than SRIM/TRIM.
Load-bearing premise
The whole method assumes that the intensity of a layer's muonic X-ray peak is directly proportional to the number of muons stopped in that layer; if capture probabilities, cascade yields, or muon transfer between layers vary, the fitted thickness will be biased.
Editorial extensions
If this is right
- Artworks with metal gildings, patinas, or corrosion crusts of ten to a few tens of microns can be measured non-destructively, recovering the thickness of each identifiable layer with about a micrometre of uncertainty.
- The same profile-matching protocol applies to any layered material in which each layer has a distinct muonic X-ray line, including buried layers beneath thick outer shells that XRF and PIXE cannot see.
- A quick open-source simulation route (SRIM/TRIM) gives usable thickness estimates, while a detailed-geometry code (GEANT4/ARBY) improves accuracy for thin or irregular layers.
- Because the method reads the full layer sequence in one momentum scan, it yields not just the coating thickness but also information about the substrate and intermediate layers, such as nickel under gold.
- For irregular or porous layers, the fit can absorb thickness variations as an effective density reduction, allowing the method to flag uneven gilding even when the nominal thickness is unknown.
Reading between the lines
- If the stopped-muon-to-X-ray proportionality is verified layer by layer, the method could be extended from thickness fitting to quantitative per-layer composition, because the same machinery would then connect absolute peak intensities to elemental concentrations.
- The reported ±1 µm uncertainty is set by the momentum spread and by the chi-square step size; a direct measurement of the beam's momentum distribution at the sample would likely tighten the error and remove the need to tune spread and source distance by hand.
- The effective-density trick used for the uneven amalgam-gilded bronze suggests a testable extension: porosity of a gilding could be estimated from the density reduction needed to reach a best fit, giving conservators a quantitative measure of layer quality.
- The same reduced-chi-square matching could be transferred to other depth-sensitive probes, but the muon's centimetre reach makes it uniquely suited to layered metal artefacts, where competing surface techniques lose signal.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a data-analysis protocol for muonic atom X-ray emission spectroscopy (μ-XES) in which Monte Carlo simulations (GEANT4/ARBY and SRIM/TRIM) are matched to measured muon stopping depth profiles to determine the thickness of a superficial gold layer. The method is demonstrated on two sets of gilded copper-alloy standards: three electroplated foils (EPA, EPB, EPC) and two amalgam-gilded foils (SM3, EM2). The simulated depth profiles are compared with measured normalised peak areas via reduced chi-square, with the gold thickness as the fitted parameter. The authors report best-fit thicknesses of 3.5±1, 4.5±1, and 7.5±1 µm for EPA, EPB, and EPC (SEM: 3.3±0.2, 4.6±0.6, 7.3±0.8 µm), 11±1 µm for SM3 (SEM: 11±1 µm), and 5±1 µm for EM2, the last with no SEM value listed in Table 2. The central claim is that the two simulation tools, despite their different approaches, yield consistent thickness estimates that agree with the SEM characterisation to within about ±1 µm.
Significance. If the ±1 µm accuracy claim holds, the protocol would be a valuable non-destructive depth-profiling tool for heritage science, where gilding and alteration layers of tens of microns are otherwise difficult to assess. The paper's strengths include the use of two independent simulation codes, a straightforward reduced-chi-square comparison criterion, open data links for the ISIS beamtime, and a clear validation design against SEM on multiple standard samples. The approach is falsifiable in principle: each sample's fitted thickness can be checked against an independent reference. However, the significance is substantially reduced by the incomplete validation record for EM2 and by the post-hoc adjustment of several physical inputs, which means the reported ±1 µm uncertainty is not yet established as a general systematic error budget.
major comments (3)
- [§2.2 and Table 2]
- [§2.1 and Discussion]
- [§1.1 and §1.4]
minor comments (4)
- [Figure 4 and Figure 5]
- [§2.2]
- [Figure 4d]
- [§2.1]
Circularity Check
Validation is partially circular: the gold-thickness search range is centered on the SEM values later used as the reference, so the claimed agreement is partly built into the scan.
-
fitted input called prediction
[Section 2.1, first paragraph (Electroplated samples)]
"For each sample, simulations with different gold layer sizes were performed, in a range of ±0.5 µm from the SEM average value."
The fitted thickness is not searched freely; it is confined to a 1 µm window centered on the SEM thickness that is then invoked as the independent validation. Because the paper states that the reduced chi-square is almost equal for 0.5 µm differences, any best fit inside this window is automatically consistent with SEM at the reported level. Thus the agreement between the simulated thickness and the SEM thickness is guaranteed by the choice of search range, not demonstrated by the simulation. This is a fitted input presented as an independent confirmation.
-
fitted input called prediction
[Section 2.2, amalgam samples, SM3 paragraph]
"Here, only size was changed, in a range from 8 to 12 μm. The best fit was reached with a thickness of 11 ± 1 µm µm and standard gold density, in agreement with the SEM measurement"
For SM3, the simulation search range is again centered on the SEM value of 11 ± 1 µm. The fitted result of 11 ± 1 µm is therefore not an independent determination; the allowed thickness interval was chosen around the reference value before fitting. The subsequent statement that the simulation 'is in agreement with SEM' is a consequence of the prior choice of the scan interval, making this validation step circular by construction.
full rationale
The core workflow is an inverse fit: measured muonic X-ray depth profiles are compared with simulated stopping profiles for assumed gold thicknesses, and the thickness is chosen by minimizing reduced chi-square. Inverse fitting is not itself circular. However, the paper's validation claim is significantly weakened because, for EPA, EPB, EPC, and SM3, the thickness search range is anchored to the SEM thickness that is later cited as the reference. Since chi-square is reported to be nearly flat over 0.5 µm steps, any fit value in the SEM-centered scan is within the stated ±1 µm uncertainty, so the 'agreement with SEM' is largely built into the search grid. Additional tuned inputs, including raising the momentum spread from 4% to 5% and adjusting source-to-sample distance for EP_C, are fit to the same measured profiles; the paper admits the spread 'was not extensively investigated' and that positioning represented a source of error. The ±1 µm thickness uncertainty is conditional on these tuned parameters and excludes their systematic effect. EM2 provides a further unaddressed inconsistency: the fitted thickness is 5 ± 1 µm while the SEM value is 11 ± 1 µm, and Table 2 omits the SEM entry for EM2, leaving the general claim of agreement unsupported. No load-bearing self-citation or imported uniqueness theorem is present; ARBY and TRIM are externally developed tools, and the cited prior uses are not used to justify the thickness result. The method has independent content, but the headline validation is partially circular because the reference values are used to define the fitting search space.
Assumptions & free parameters
free parameters (4)
- Gold layer thickness (per sample) =
EPA 3.5, EPB 4.5, EPC 7.5, SM3 11, EM2 5 um
- Beam momentum spread (dp/p) =
5% (increased from nominal 4%)
- Source-to-sample distance (EPC) =
9.7 cm (nominal 10 cm)
- Gold layer density (EM2) =
80% of nominal (15.46 g/cm3)
assumptions (4)
- domain assumption Muon stopping in matter is accurately modeled by GEANT4 and by TRIM using a hydrogen ion with one-ninth the proton mass (effectively muon mass).
- domain assumption The number of stopped muons in a layer is proportional to the emitted X-ray intensity.
- domain assumption Detection efficiency is approximately constant over the 102-134 keV range, so no efficiency correction is applied.
- domain assumption The samples consist of uniform, flat layers with sharp interfaces and the nominal compositions.
Cite this review
Pith. "Pith review of Depth profiling the elemental composition with negative muons: Monte Carlo based tools for improved data analysis." pith.science (2026). https://pith.science/paper/3DPIGNW7
@misc{pith2026250504679,
author = {Pith},
title = {Pith review of: Depth profiling the elemental composition with negative muons: Monte Carlo based tools for improved data analysis},
year = {2026},
howpublished = {\url{https://pith.science/paper/3DPIGNW7}},
note = {Machine review of arXiv:2505.04679}
}
read the original abstract
Gildings, patinas and alteration crusts are common features of many heritage artefacts, especially for metals. Their size depends on many factors, like the manufacturing method for gildings or the conservation state for alteration crusts: in some cases, it can be in the scale of the tens of microns. Such thickness would be difficult to investigate with classical non-destructive methods and would prevent getting information from the bulk of the sample. This work proposes an innovative approach for the study of multi-layered materials with the Muonic atom X-ray Emission Spectroscopy technique ({\mu}-XES). Based on the detection of the high-energy X-rays emitted after the muon capture by the atom, this method is characterised by a remarkable penetration depth (from microns to cm). From the surface to the bulk, this technique can evaluate the variation of the elemental composition as a function of depth. The paper focuses on providing an improved interpretation of {\mu}-XES data by coupling the analysis with the use of two Monte Carlo simulation software, GEANT4/ARBY and SRIM/TRIM. With these two software, it is possible to replicate the negative muon experiments and compare the experimental and simulated outputs to address the size of a given layer. To validate this approach, a set of standard gilded bronze and brass foils were measured at the ISIS Neutron and Muon source. From simulations, it was possible to evaluate the thickness of the superficial gold layer, with results in agreement with the preliminary SEM characterisation of the samples.
Reference graph
Works this paper leans on
-
[1]
External Beam IBA Measurements for Cultural Heritage,
M. Chiari, “External Beam IBA Measurements for Cultural Heritage,” Applied Sciences (Switzerland), vol. 13, no. 5, Mar. 2023, doi: 10.3390/app13053366
-
[2]
C. Bouvier, S. Van Nuffel, P. Walter, and A. Brunelle, “Time-of-flight secondary ion mass spectrometry imaging in cultural heritage: A focus on old paintings,” Journal of Mass Spectrometry, vol. 57, no. 1, p. e4803, 2022, doi: https://doi.org/10.1002/jms.4803
-
[3]
R. Escobar Galindo, R. Gago, D. Duday, and C. Palacio, “Towards nanometric resolution in multilayer depth profiling: A comparative study of RBS, SIMS, XPS and GDOES,” 2010, Springer Verlag. doi: 10.1007/s00216-009-3339-y
-
[4]
Comparison of depth profiling techniques using ion sputtering from the practical point of view,
S. Oswald and S. Baunack, “Comparison of depth profiling techniques using ion sputtering from the practical point of view,” 2003
work page 2003
-
[5]
B. Kanngießer, W. Malzer, I. Mantouvalou, D. Sokaras, and A. G. Karydas, “A deep view in cultural heritage-confocal micro X-ray spectroscopy for depth resolved elemental analysis,” Appl Phys A Mater Sci Process, vol. 106, no. 2, pp. 325–338, Feb. 2012, doi: 10.1007/s00339- 011-6698-0
doi:10.1007/s00339- 2012
-
[6]
The Role of PIXE and XRF in Heritage Science: The INFN-CHNet LABEC Experience,
L. Sottili et al., “The Role of PIXE and XRF in Heritage Science: The INFN-CHNet LABEC Experience,” Applied Sciences, vol. 12, no. 13, 2022, doi: 10.3390/app12136585
- [7]
-
[8]
A. D. Hillier et al., “Muon spin spectroscopy,” Nature Reviews Methods Primers, vol. 2, no. 1, p. 4, 2022, doi: 10.1038/s43586-021-00089-0
Show all 34 references
-
[9]
S. J. Blundell, R. De Rienzi, T. Lancaster, and F. L. Pratt, Muon Spectroscopy: an Introduction. Oxford University press, 2021
2021
-
[10]
A Novel Non-Destructive Technique for Cultural Heritage: Depth Profiling and Elemental Analysis Underneath the Surface with Negative Muons,
M. Cataldo, M. Clemenza, K. Ishida, and A. D. Hillier, “A Novel Non-Destructive Technique for Cultural Heritage: Depth Profiling and Elemental Analysis Underneath the Surface with Negative Muons,” May 01, 2022, MDPI. doi: 10.3390/app12094237
2022 doi
-
[11]
Development of elemental analysis by muonic X-ray measurement in J- PARC,
K. Ninomiya et al., “Development of elemental analysis by muonic X-ray measurement in J- PARC,” J Phys Conf Ser, vol. 225, pp. 1–3, 2010, doi: 10.1088/1742-6596/225/1/012040
2010 doi
-
[12]
Elemental Analysis of Bronze Artifacts by Muonic X-ray Spectroscopy,
K. Ninomiya et al., “Elemental Analysis of Bronze Artifacts by Muonic X-ray Spectroscopy,” Physical Society of Japan, Sep. 2015. doi: 10.7566/jpscp.8.033005
2015 doi
-
[13]
Understanding Roman Gold Coinage Inside Out,
G. A. Green, K. Ishida, B. V. Hampshire, K. Butcher, A. M. Pollard, and A. D. Hillier, “Understanding Roman Gold Coinage Inside Out,” J Archaeol Sci, vol. 134, Oct. 2021, doi: 10.1016/j.jas.2021.105470
2021
-
[14]
The non-destructive investigation of a late antique knob bow fibula (Bügelknopffibel) from Kaiseraugst/CH using Muon Induced X-ray Emission (MIXE),
S. Biswas et al., “The non-destructive investigation of a late antique knob bow fibula (Bügelknopffibel) from Kaiseraugst/CH using Muon Induced X-ray Emission (MIXE),” Herit Sci, vol. 11, no. 1, Dec. 2023, doi: 10.1186/s40494-023-00880-0. © 2025. This manuscript version is mad...
2023 doi
-
[15]
Use of Muonic X Rays for Nondestructive Analysis of Bulk Samples for Low Z Constituents,
J. J. Reidy, R. L. Hutson, H. Daniel1, and K. Springer, “Use of Muonic X Rays for Nondestructive Analysis of Bulk Samples for Low Z Constituents,” Anal Chem, vol. 50, no. 1, pp. 40–44, 1978, doi: 10.1021/ac50023a015
1978 doi
-
[16]
Application of muonic X-ray techniques to the elemental analysis of archeological objects,
E. Köhler, R. Bergmann, H. Daniel, P. Ehrhart, and F. J. Hartmann, “Application of muonic X-ray techniques to the elemental analysis of archeological objects,” Nuclear Instruments and Methods, vol. 187, no. 2–3, pp. 563–568, 1981, doi: 10.1016/0029-554X(81)90389-X
1981 doi
-
[17]
FLUKA: A multi-particle transport code (Program version 2005),
A. Ferrari, P. R. Sala, A. Fasso, and J. Ranft, “FLUKA: A multi-particle transport code (Program version 2005),” Jul. 2005, doi: 10.2172/877507
2005 doi
-
[18]
GEANT4 - A simulation toolkit,
S. Agostinelli et al., “GEANT4 - A simulation toolkit,” Nucl Instrum Methods Phys Res A, vol. 506, no. 3, pp. 250–303, 2003, doi: 10.1016/S0168-9002(03)01368-8
2003 doi
-
[19]
MCNP User Manual . Code Versione 6.2,
C. J. Werner, “MCNP User Manual . Code Versione 6.2,” Los Alamos National Laboratory, 2017
2017
-
[20]
The projected background for the CUORE experiment,
C. Alduino et al., “The projected background for the CUORE experiment,” The European Physical Journal C, vol. 77, no. 8, p. 543, 2017, doi: 10.1140/epjc/s10052-017-5080-6
2017 doi
-
[21]
Measurement of the two-neutrino double-beta decay half-life of $$^{130}$$Te with the CUORE-0 experiment,
C. Alduino et al., “Measurement of the two-neutrino double-beta decay half-life of $$^{130}$$Te with the CUORE-0 experiment,” The European Physical Journal C, vol. 77, no. 1, p. 13, 2017, doi: 10.1140/epjc/s10052-016-4498-6
2017 doi
-
[22]
SRIM – The stopping and range of ions in matter (2010),
J. F. Ziegler, M. D. Ziegler, and J. P. Biersack, “SRIM – The stopping and range of ions in matter (2010),” Nucl Instrum Methods Phys Res B, vol. 268, no. 11–12, pp. 1818–1823, Jun. 2010, doi: 10.1016/J.NIMB.2010.02.091
2010 doi
-
[23]
Thickness determination of metal multilayers by ED-XRF multivariate analysis using Monte Carlo simulated standards,
W. Giurlani, E. Berretti, A. Lavacchi, and M. Innocenti, “Thickness determination of metal multilayers by ED-XRF multivariate analysis using Monte Carlo simulated standards,” Anal Chim Acta, vol. 1130, pp. 72–79, 2020, doi: https://doi.org/10.1016/j.aca.2020.07.047
2020 doi
-
[24]
Monte Carlo simulations of XRF intensities in non- homogeneous matrices,
J. A. Helsen and B. A. R. Vrebos, “Monte Carlo simulations of XRF intensities in non- homogeneous matrices,” Spectrochim Acta Part B At Spectrosc, vol. 39, no. 6, pp. 751–759, 1984, doi: https://doi.org/10.1016/0584-8547(84)80083-X
1984 doi
-
[25]
The Implementation of MuDirac in Geant4: A Preliminary Approach to the Improvement of the Simulation of the Muonic Atom Cascade Process,
M. Cataldo et al., “The Implementation of MuDirac in Geant4: A Preliminary Approach to the Improvement of the Simulation of the Muonic Atom Cascade Process,” Condens Matter, vol. 8, no. 4, Dec. 2023, doi: 10.3390/condmat8040101
2023 doi
-
[26]
The practice and characterization of historic fire gilding techniques,
K. Anheuser, “The practice and characterization of historic fire gilding techniques,” JOM, vol. 49, no. 11, pp. 58–62, 1997, doi: 10.1007/s11837-997-0015-6
1997 doi
-
[27]
Gilding through the ages,
A. Oddy, “Gilding through the ages,” Gold Bull, vol. 14, no. 2, pp. 75–79, 1981, doi: 10.1007/BF03214601
1981 doi
-
[28]
Plating and Surface Treatments on Ancient Metalwork,
A. Giumlia-Mair, “Plating and Surface Treatments on Ancient Metalwork,” Advances in Archaeomaterials, vol. 1, no. 1, pp. 1–26, Dec. 2020, doi: 10.1016/j.aia.2020.10.001
2020 doi
-
[29]
Electroplating for decorative applications: Recent trends in research and development,
W. Giurlani et al., “Electroplating for decorative applications: Recent trends in research and development,” Aug. 01, 2018, MDPI AG. doi: 10.3390/coatings8080260
2018 doi
-
[30]
Muons at ISIS,
A. D. Hillier, J. S. Lord, K. Ishida, and C. Rogers, “Muons at ISIS,” Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 377, no. 2137, pp. 1–7, 2019, doi: 10.1098/rsta.2018.0064. © 2025. This manuscript version is made ava...
2019
-
[31]
The RIKEN-RAL pulsed Muon Facility,
T. Matsuzaki, K. Ishida, K. Nagamine, I. Watanabe, G. H. Eaton, and W. G. Williams, “The RIKEN-RAL pulsed Muon Facility,” 2001
2001
-
[32]
Electronic catalogue of muonic X-rays,
Zinatulina Daniya et al., “Electronic catalogue of muonic X-rays,” EPJ Web Conf., vol. 177, p. 3006, 2018, doi: 10.1051/epjconf/201817703006
2018
-
[33]
Charge-distribution parameters, isotope shifts, isomer shifts, and magnetic hyperfine constanst from muonic atoms,
R. Engfer, H. Schneuwly, J. L. Vuilleumier, H. K. Walter, and A. Zehnder, “Charge-distribution parameters, isotope shifts, isomer shifts, and magnetic hyperfine constanst from muonic atoms,” At Data Nucl Data Tables, vol. 14, pp. 509–597, 1974
1974
-
[34]
Mudirac: A Dirac equation solver for elemental analysis with muonic X-rays,
S. Sturniolo and A. Hillier, “Mudirac: A Dirac equation solver for elemental analysis with muonic X-rays,” X-Ray Spectrometry, no. June, pp. 1–17, 2020, doi: 10.1002/xrs.3212
2020 doi
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.