REVIEW 2 major objections 5 minor 15 references
Comprehensive characterization of a YAG:Ce scintillator: light yield, alpha quenching and pulse-shape discrimination
T0 review · 2 major / 5 minor · reviewed 2026-07-15 · grok-4.5
Pith's one-line read YAG:Ce alpha quenching falls from 0.17 to 0.10 between 5.5 and 1 MeV, while pulse shape still separates alphas from gammas.
desk verdict Solid, usable characterization of commercial YAG:Ce: continuous α QF curve 1–6 MeV, τ_long temperature dependence, and quantitative PSD numbers that groups can actually cite. 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
Reconstruction of measured waveforms as random sums of single-SiPM-cell pulses drawn from a double-exponential emission law, which simultaneously yields the short and long decay constants, their relative weights, and the absolute number of detected photons used for light-yield and quenching calculations.
What would settle it
Repeat the pressure-scanned alpha measurements with an independent electron-beam or Compton-edge calibration that covers the same 1–6 MeV window and check whether the extracted quenching factors remain inside the reported 0.10–0.17 band.
Extended reading notes
Core claim
For a 1 cm YAG:Ce cube read out by a SiPM, the alpha quenching factor decreases smoothly from approximately 0.17 to 0.10 as the energy deposited by alphas falls from 5.49 MeV to about 1 MeV, while a single partial-charge pulse-shape parameter separates those alphas from equal-amplitude gammas at the 1–2 σ level and full-energy alphas from 511 keV gammas at about 2.3 σ.
Load-bearing premise
The light yield for electrons and gammas is taken as constant above roughly 100 keV and is therefore fixed solely by the 511 keV peak; any real non-linearity in that range would stretch or compress the entire quenching-factor curve.
Editorial extensions
If this is right
- Energy reconstruction for alphas in YAG:Ce must apply an energy-dependent correction that reaches a factor of ten near 1 MeV.
- Particle identification remains usable down to amplitudes matching the 511 keV gamma peak, albeit with reduced separation power.
- Cooling to -50 °C lengthens the slow component by a factor of two without spoiling light yield, so timing cuts must be retuned with temperature.
- The same vacuum-chamber method can be extended below 1 MeV by thinning the crystal to suppress the competing 60 keV gamma line.
Reading between the lines
- Because quenching continues to fall toward 1 MeV, applications that rely on low-energy alphas (e.g., surface contamination monitors) will need an explicit energy-dependent quenching model rather than a single constant.
- The observed temperature-independent light yield combined with a slowing decay suggests that YAG:Ce may remain competitive for outdoor or space instruments where power budgets forbid active cooling of the crystal itself.
- If the partial-charge separation scales with photon statistics, higher-light-yield YAG:Ce variants already reported in the literature should push the discrimination figure of merit well above 3 σ at the same energies.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports a laboratory characterization of a 10 imes10 imes10 mm^{3} YAG:Ce crystal read out by a Hamamatsu S13360 SiPM. Using a double-exponential reconstruction of physics pulses from single-cell templates, the authors extract short and long scintillation time constants (~63–67 ns and ~245–273 ns) and a light yield of ~19 000 photons MeV⁻¹. By varying gas pressure, gas species (air, Ar, He) and source–crystal distance they map the α quenching factor from ~0.17 at 5.49 MeV down to ~0.10 near 1 MeV. Temperature scans from room temperature to –50 °C show that τ_long roughly doubles while light output remains stable. Finally, a partial-charge pulse-shape parameter separates full-energy α events from 511 keV γ events at ~2.3 σ (~1.3 σ when amplitudes are matched).
Significance. The work supplies a self-consistent, multi-configuration data set for a commercially available YAG:Ce crystal that is already used in beam diagnostics, space instrumentation and low-background applications. The energy-dependent QF curve, the temperature evolution of the slow component, and the quantitative PSD figures of merit are directly usable for detector design and Monte-Carlo modeling. Strengths include the explicit cross-check of QF with three gases and two geometries, the transparent single-cell reconstruction method (Fig. 1), and the clear statement of the linearity assumption used for the β/γ calibration. The results are incremental rather than foundational, yet they fill a practical gap in the literature for this material.
major comments (2)
- Sec. II B and Fig. 2: the entire QF curve is scaled by the light yield fixed solely from the 511 keV ²²Na peak under the assumption of energy-independent β/γ response above ~100 keV. While this is conventional and acknowledged by the authors, a single-point calibration leaves an overall systematic that is not quantified. A brief cross-check with an additional γ line (e.g., 1.27 MeV ²²Na or an external ¹³⁷Cs source) or a short discussion of the expected non-linearity from literature would strengthen the absolute scale of the reported QF values.
- Sec. III: the PSD performance is quoted as ~2.3 σ (full-energy α) and ~1.3 σ (amplitude-matched). The text does not specify how σ is computed (Gaussian fit, RMS of the partial-charge distributions, or clustering purity). Clarifying the definition and, if possible, providing the figure-of-merit (e.g., FoM = |μ_α – μ_γ| / (σ_α + σ_γ)) would make the discrimination claim fully reproducible.
minor comments (5)
- Table I: the fractional short-component weight F_short is listed as 0.10 (α) and 0.15 (γ); a short sentence on how these fractions were obtained (free fit or constrained) would help readers reproduce the reconstruction.
- Fig. 3: the vertical axis is labeled “Norm. τ_long(T)” but the caption and text state that τ_long increases by a factor of two; confirming whether the plotted quantity is τ_long(T)/τ_long(RT) or the inverse would remove ambiguity.
- Introduction and Sec. II A: the light-yield value of ~19 000 photons MeV⁻¹ is lower than some recent literature reports (up to 3.5×10⁴). A brief remark on crystal quality, Ce concentration or light-collection geometry would place the number in context.
- Throughout: the abstract and body use both “YAG:Ce” and “\YAG”; consistent notation would improve readability.
- References: Ref. [10] is cited for the reconstruction method and experimental setup; ensuring that the arXiv or journal version is publicly available at the time of publication would aid reproducibility.
Circularity Check
No circularity: all reported quantities (decay constants, light yield, QF, PSD) are extracted from new experimental data under standard phenomenological models.
full rationale
This is a pure experimental characterization paper. Scintillation parameters (τ_short, τ_long, F) are obtained by fitting a conventional double-exponential model to newly recorded waveforms; light yield is counted from the number of single-cell pulses needed to reconstruct the 511 keV peak; the quenching factor is the ratio of measured α light output to the β/γ light yield calibrated at 511 keV, with deposited energy supplied by independent Geant4 simulations of the gas path; temperature dependence and partial-charge PSD separations are likewise direct measurements. None of these steps is definitional, none re-uses a fitted parameter as a subsequent prediction, and the only self-citation (Ref. [10]) supplies experimental methodology rather than a uniqueness theorem or ansatz that forces the present results. The double-exponential form is a standard phenomenological ansatz, not a circular definition of the measured times. The paper is therefore self-contained against external benchmarks and scores 0.
Assumptions & free parameters
free parameters (2)
- τ_short, τ_long, F_short =
α: 63 ns / 273 ns / 0.10; γ: 67 ns / 245 ns / 0.15
- light yield (photons MeV⁻¹) =
~19 000
assumptions (3)
- domain assumption Scintillation light for β/γ events is linear with deposited energy above ~100 keV, so a single 511 keV calibration point suffices for the MeV range.
- domain assumption A two-component exponential adequately describes the scintillation kinetics; any ~1 µs component can be neglected.
- domain assumption Geant4 correctly computes the mean energy lost by α particles in air, argon and helium at the stated pressures and distances.
Cite this review
Pith. "Pith review of Comprehensive characterization of a YAG:Ce scintillator: light yield, alpha quenching and pulse-shape discrimination." pith.science (2026). https://pith.science/paper/J4ABXV3W
@misc{pith2026260306429,
author = {Pith},
title = {Pith review of: Comprehensive characterization of a YAG:Ce scintillator: light yield, alpha quenching and pulse-shape discrimination},
year = {2026},
howpublished = {\url{https://pith.science/paper/J4ABXV3W}},
note = {Machine review of arXiv:2603.06429}
}
abstract
Solid-state scintillators are widely used in particle and applied physics due to their versatility and resistance to diverse environments and operating conditions. This broad range of applications calls for thorough characterization of scintillating crystals. Among these materials, cerium-doped yttrium aluminum garnet (YAG:Ce) is a promising scintillator owing to its favorable timing characteristics, high light yield, good mechanical properties, and chemical stability. In this work, we report a comprehensive experimental characterization of a YAG:Ce crystal exposed to both $\gamma$ and $alpha$ radiation. We extract the scintillation decay time and light yield, and study their evolution from room temperature down to approximately $-50 ^\circ$ C. We perform a detailed investigation of the quenching factor for \al particles in the energy range from about $6$ MeV down to $1$ MeV, finding a value that decreases from approximately $0.17$ to $0.10$. We also explore the possibility of pulse-shape discrimination based on the different signal evolution depending on the interaction type, demonstrating strong classification capabilities. These results provide a detailed assessment of the performance of \YAG for radiation-detection applications and offer insight into its potential use in environments requiring reliable particle identification and stable response across a wide range of operating conditions.
Reference graph
Works this paper leans on
-
[1]
J. B. Birks, Proc. Phys. Soc. A64, 874 (1951)
1951
-
[2]
van Eijk, J
C. van Eijk, J. Andriessen, P. Dorenbos, and R. Visser, Nucl. Instrum. Methods A348, 546 (1994)
1994
-
[3]
Moszy´ nskiet al., Nucl
M. Moszy´ nskiet al., Nucl. Instrum. Meth. A345, 461 (1994)
1994
-
[4]
Zapadl´ ıket al., Opt
O. Zapadl´ ıket al., Opt. Mater. X15, 100165 (2022)
2022
-
[5]
Lucchiniet al., IEEE Trans
M. Lucchiniet al., IEEE Trans. Nucl. Sci.63, 586 (2016)
2016
-
[6]
Ludziejewskiet al., Nucl
T. Ludziejewskiet al., Nucl. Instrum. Meth. A398, 287 (1997)
1997
-
[7]
M. Kim, H. J. Kim, J. Y. Cho, S. Kaewjaeng, and J. Kaewkhao, Appl. Rad. Iso.145, 126 (2019)
2019
-
[8]
Epic crystal,
Epic Crystal, “Epic crystal,” (2026)
2026
Show all 15 references
-
[9]
Hamamatsu Photonics,S13360 Series MPPC (Multi- Pixel Photon Counter) Datasheet, Hamamatsu Photon- ics (2021),https://www.hamamatsu.com/content/dam/ hamamatsu-photonics/sites/documents/99_SALES_ LIBRARY/ssd/s13360_series_kapd1052e.pdf
2021
-
[10]
Gironiet al., The European Physical Journal Plus 140, 415 (2025)
L. Gironiet al., The European Physical Journal Plus 140, 415 (2025)
2025
-
[11]
Agostinelliet al.(GEANT4 Group), Nucl
S. Agostinelliet al.(GEANT4 Group), Nucl. Instrum. Meth. A506, 250 (2003)
2003
-
[12]
Scintillation Detectors for Charged Particles and Photons,
P. Lecoq, “Scintillation Detectors for Charged Particles and Photons,” inParticle Physics Reference Library. Vol- ume 2: Detectors for Particles and Radiation, edited by C. W. Fabjan and H. Schopper (Springer, 2020) pp. 45– 89
2020
-
[13]
V. I. Tretyak, Astropart. Phys.33, 40 (2010)
2010
-
[14]
Beretta, S
M. Beretta, S. Capelli, L. Gironi, E. Previtali, and M. Sisti, J. Instrum.12, P04007 (2017)
2017
-
[15]
H. M. Qadr and A. M. Hamad, Res. Rev. J. Pure Appl. Phys.5, 23 (2017)
2017
Reviewed July 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.