Pith. sign in

REVIEW 4 major objections 5 minor 14 references

Positron Annihilation Lifetime Spectroscopy Using Fast Scintillators and Digital Electronics

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

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

Pith's one-line read This paper reports a fully digital positron-annihilation lifetime spectroscopy setup based on BC-418 plastic scintillators and 500 MS/s digitization with sinc interpolation and digital constant-fraction discrimination, reaching 198.3 ±…

desk verdict The 198 ps timing result is plausible and well-measured, but the PALS demonstration is under-validated: the source background is uncorrected and the component intensities disagree with the cited quartz reference. read the letter →

arxiv 1908.04434 v1 pith:E2RD2JZI submitted 2019-08-12 physics.ins-det nucl-exphysics.data-an

Ming Fang , Nathan Bartholomew , Angela Di Fulvio This is my paper · ORCID
classification physics.ins-detnucl-exphysics.data-an PACS 29.40.Mc78.70.Bj
keywords positronannihilationlifetimespectroscopyorganicscintillatorsBC-418digitalpulseprocessingconstantfractiondiscriminationsincinterpolationtimeresolutionsingle-crystalquartz
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

Positron annihilation lifetime spectroscopy (PALS) times the gamma ray that announces a positron's birth against the gamma ray from its annihilation, and that time encodes how the positronium interacts with vacancies and pores in a material. The paper aims to show that this timing can be done digitally, without analog timing modules, using fast organic scintillators and a 500 MS/s digitizer. With the fastest detector pair, BC-418 plastic scintillators, the setup reaches a time resolution of 198.3 ± 0.8 ps FWHM ($\sigma = 84.2 \pm 0.3$ ps). Applied to single-crystal quartz, it fits positronium lifetimes of $\tau_1 = 159 \pm 9$ ps and $\tau_2 = 366 \pm 22$ ps, in agreement with reference values of 156 ps and 358 ps for that material. The result implies that material-defect studies can get picosecond timing from small plastic detectors and digital pulse processing rather than bulky analog timing electronics.

What carries the argument

The load-bearing element is the arrival-time estimator: a Gaussian-windowed truncated sinc interpolation (Eqs. 4–6) reconstructs each pulse between 2-ns samples, and a digital constant-fraction discriminator (Eq. 7) forms $F \cdot S(i) - S(i-\Delta)$ and takes its zero crossing as the timestamp. The two free parameters, fraction $F$ and delay $\Delta$, are optimized for each detector pair by Gaussian-fitting the 60Co coincidence time-difference spectrum and minimizing its FWHM. Interpolation does the specific job of recovering the rising edge and true peak that sparse sampling misses; the paper shows this reduces the time-difference spread and removes the skewness introduced by too-small $F$ values.

What would settle it

Record 60Co coincidence pulses with a 20 GS/s digitizer, down-sample the same waveforms to 500 MS/s, and run this paper's interpolation plus CFD on the down-sampled version. If the recovered arrival-time distribution has a sigma above about 84 ps or shows a systematic shift relative to the 20 GS/s timestamps, the unbiased-reconstruction premise is wrong.

Watch

Extended reading notes

Core claim

The paper establishes that a positron-annihilation lifetime spectrometer built from two BC-418 fast plastic scintillators, a 14-bit 500 MS/s digitizer, and a two-step digital timing algorithm can reach 198.3 ± 0.8 ps FWHM time resolution, sufficient to separate the short para-positronium component from the longer ortho-positronium and free-positron component in a defect-free material. Using this setup on a 22Na source sandwiched between two single-crystal quartz samples, the fitted lifetimes are $\tau_1 = 159\pm9$ ps and $\tau_2 = 366\pm22$ ps, consistent with the 156 ps and 358 ps values reported for the same material. The implication is that the digital pipeline itself drives timing performance: interpolation narrows the measured time-difference distribution, and the CFD parameters change the spread substantially across the tested detector pairs.

Load-bearing premise

The whole timing accuracy rests on the premise that voltages sampled every 2 nanoseconds can be faithfully reconstructed into a smooth pulse, so that the constant-fraction zero crossing is an unbiased arrival time; the reconstruction's window size is not specified, so the reader cannot check its error from the paper alone.

Editorial extensions

If this is right

  • The BC-418 digital setup resolves two positronium lifetime components in quartz, so vacancy and defect studies can be done without analog timing electronics.
  • The 198 ps FWHM makes it possible to discriminate the spin-singlet para-positronium component, something earlier systems with roughly 330 ps resolution could not do for this material.
  • The same interpolation-plus-CFD algorithm transfers to any fast-timing application, including nuclear medicine and radiation imaging, where timestamp accuracy is the limiting step.
  • Interpolation alone improves the timing resolution by about 33 ps at the BC-418 operating point and also makes skewed time-difference histograms symmetric.
  • The authors state the optimized setup is intended next for analyzing vacancies and damage in radiation detectors exposed to high fluence.

Reading between the lines

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

  • If the reconstruction is unbiased, the same 500 MS/s digitizer plus interpolation pipeline could be ported to positron emission tomography, where roughly 200 ps coincidence timing is a practical target and commodity digitizers could replace analog constant-fraction discriminators.
  • Because BC-418 does not provide pulse-shape discrimination, the current readout separates only two lifetime components; pairing the same timing algorithm with EJ-309 or EJ-276, which do support pulse-shape discrimination, would let users reject scattered events and resolve additional components.
  • A direct stress test would be to record waveforms at 20 GS/s, down-sample them to 500 MS/s, and compare the paper's timestamps with the original ones; the reported $\sigma = 84.2$ ps would be falsified by any systematic offset larger than that on synthetic pulses of known arrival time.
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 presents a digital positron annihilation lifetime spectroscopy (PALS) setup built from organic scintillators, a 500 MS/s digitizer, and a pulse-processing chain that interpolates waveforms and applies a constant-fraction discrimination (CFD) algorithm. The authors test three detector pairs (BC-418 plastic, EJ-309 liquid, EJ-276 plastic), optimize the CFD parameters (F and Δ) and an energy threshold for each pair, and report a best time resolution of 198.3 ± 0.8 ps FWHM (σ = 84.2 ± 0.3 ps) with BC-418 detectors after rejecting pulses below 600 keVee. Using this optimized pair, they measure positron lifetimes in single-crystal quartz and report τ1 = 159 ± 9 ps and τ2 = 366 ± 22 ps, which they claim are in good agreement with literature values. The central claims are the sub-200 ps timing resolution and the demonstration that the setup can resolve the p-Ps and o-Ps components in quartz.

Significance. If the timing result is reliable, the paper provides a useful, cost-effective digital timing method that achieves competitive time resolution with fast plastic scintillators and commercially available digitizers. The open-source ROOT-based pulse-processing code is a practical contribution that could benefit PALS and other fast-timing applications. However, the PALS demonstration is not yet convincing: the fitted component intensities disagree strongly with the cited reference, and the fitting model omits the known source/Kapton contribution. The timing measurement itself is well described, with explicit parameter optimization and Gaussian fits, but the interpolation parameters and Nyquist assumptions are not fully specified, leaving some uncertainty in the timing method. Overall, the instrumentation part is promising, while the material-characterization claim needs substantial additional support.

major comments (4)
  1. [Sec. 2.3 and Eq. (8)] The 22Na source is sealed between two Kapton foils, yet the PALS fitting model in Eq. (8) contains no source/background term and no source-only (blank) measurement is reported. Kapton has a characteristic positron lifetime near 380 ps, so even a modest fraction of annihilations in the source assembly would inflate the fitted τ2 and alter the intensity balance. To support the assignment of τ2 to quartz o-Ps pick-off, the authors must either measure and subtract the source contribution, include a source term in the fit, or demonstrate with a blank spectrum that the source contribution is negligible.
  2. [Table 3 and Sec. 3.2] The fitted intensities I1/I2 = (60 ± 6)/(40 ± 6) disagree with the cited reference values (84.2 ± 0.3)/(15.8 ± 0.3) by more than four standard deviations, even though the lifetimes agree. Because the intensities directly indicate which fraction of annihilations belongs to each component, this discrepancy means the two fitted components are not cleanly identified with quartz p-Ps and o-Ps/free positron. The statement in Sec. 3.2 that the lifetimes are 'in good agreement' is therefore incomplete; the authors need to explain the intensity mismatch or substantially temper the claim that the setup can separate the positronium states.
  3. [Sec. 3.2 and Fig. 8] The two-component fit is adopted because the authors assert that the third (trapped o-Ps) component is undetectable in their high-purity quartz, but no statistical goodness-of-fit test or sensitivity analysis is reported. A chi-square value, residuals plot, or equivalent would allow the reader to judge whether the two-exponential model with no background term is adequate. Without such a test, the reported lifetimes cannot be distinguished from an arbitrary decomposition of a spectrum that may include source-related contributions.
  4. [Sec. 2.2 and Eqs. (4)-(6)] The interpolation parameters L (window width) and T (Gaussian decay constant) in Eqs. (4)-(6) are never specified, and the Nyquist condition for the 500 MS/s sampling of sub-nanosecond BC-418 pulses is not verified. Since the CFD time stamp depends on the interpolated rising edge, the missing parameters make the timing algorithm unreproducible and leave open the possibility of systematic timing bias. Although the 60Co-based timing resolution is measured directly, this issue affects the credibility of the headline timing result and should be addressed.
minor comments (5)
  1. [Abstract vs. Sec. 3.2/Table 3] The abstract reports τ1 = 156 ± 9 ps, while the text and Table 3 report 159 ± 9 ps; these numbers should be made consistent.
  2. [Sec. 3.2] The text says 'We fitted Eq. (6) to the positron lifetime spectrum' but Eq. (6) is the interpolation formula; the intended reference is the PALS fit model Eq. (8).
  3. [Fig. 4 caption] The caption contains the typo 'Normailzed Counts'; it should read 'Normalized Counts'.
  4. [Sec. 3.1] The phrase 'The minimum 195.7 ps σ (293.4 ps FWHM)' refers to the result before the 600 keVee energy cut; this should be stated explicitly to avoid confusion with the final 84.2 ps σ in Table 2.
  5. [Sec. 3.1] The term 'DIACFD' appears in the text; this is likely a typo for 'CFD'.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the timing resolution is measured with 60Co, and the quartz lifetimes are fitted to data and compared with external literature, not derived from the inputs.

full rationale

The paper's derivation chain is self-contained and empirically grounded. The claimed time resolution of 198.3 ± 0.8 ps FWHM is obtained by Gaussian fitting of a 60Co coincidence time-difference distribution, which is an external reference measurement independent of the PALS analysis. The processing parameters F and Δ are optimized against this 60Co resolution, not against the quartz lifetimes, so no fitted parameter is renamed as a prediction. The quartz lifetimes τ1 = 159 ± 9 ps and τ2 = 366 ± 22 ps are fitted to the measured PALS spectrum using the LT10 program and compared with literature values from refs. [3] and [6]; they are not derived from any quantity defined in terms of the claimed result. The interpolation algorithm (Eqs. 4–6) is a standard sinc-based reconstruction cited to Shannon and to Warburton and Hennig, and it is applied to the measured pulses rather than used to force any outcome. The only potentially self-authored citation is ref. [13] (which includes A. Di Fulvio) for the digital CFD algorithm, but this is a methodological citation and not load-bearing for the central claims. Issues such as the unsubtracted Kapton/source background, the intensity mismatch with the reference, or the Eq. (6) vs Eq. (8) typo are correctness or presentation risks, not circular reasoning.

Assumptions & free parameters 5 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new entities. Its central claims rest on a handful of fitted tuning parameters (F, Delta, energy threshold, interpolation constants) and on standard signal-processing and spectral-modeling assumptions that are not all verified.

free parameters (5)
  • CFD fraction F (per detector pair) = BC-418: 0.4; EJ-309: 0.2; EJ-276: 0.05
    Chosen in Sec. 2.2 and Sec. 3.1 by scanning to minimize the Gaussian FWHM of the 60Co time-difference distribution; the headline 198.3 ps result depends on these choices.
  • CFD delay Delta (per detector pair) = BC-418: 4 ns; EJ-309: 10 ns; EJ-276: 6 ns
    A tuned parameter in Eq. (7), optimized along with F; reported in Table 2.
  • Energy threshold for timing = 600 keVee
    Applied post hoc in Sec. 3.1 to reject low-energy pulses and reduce timing tails; the best resolution in Table 2 is conditional on this cut.
  • Interpolation window L and Gaussian decay T = not specified
    Constants in Eqs. (5)-(6) that control the truncated sinc kernel; no values are given, so the timing resolution cannot be exactly reproduced.
  • PALS spectral fit parameters = τ1=159 ps, τ2=366 ps, I1=60%, I2=40%
    Obtained with LT10 in Sec. 3.2; these are the measured lifetimes and intensities, but they are free parameters of the fit, not independent predictions.
assumptions (5)
  • domain assumption Digitized pulses satisfy the Nyquist condition and can be exactly reconstructed by sinc interpolation.
    Invoked in Sec. 2.2 Eqs. (1)-(6); 500 MS/s sampling of sub-nanosecond scintillator pulses may not meet Nyquist, and truncation with a Gaussian window is an approximation from Ref. [12].
  • domain assumption A single Gaussian resolution function describes the system response at both 1.17/1.33 MeV (60Co) and 511 keV/1.27 MeV (PALS).
    Sec. 2.1 and Eq. (8); the resolution is measured at 60Co energies and applied to PALS energies without an energy-dependent correction.
  • ad hoc to paper The quartz PALS spectrum contains exactly two exponential components; the third (trapped o-Ps) component is undetectable.
    Sec. 3.2 states 'we believe the third component was actually undetectable' with no statistical test; this choice changes the fitted τ1 and I1.
  • domain assumption Positrons annihilating in the Kapton foils and source support are a negligible component of the spectrum.
    Not stated or measured; the 22Na source is sealed in Kapton (Sec. 2.3), which itself has a characteristic positronium lifetime.
  • domain assumption Constant-fraction discrimination zero crossing is an unbiased estimator of gamma arrival time.
    Sec. 2.2 Eq. (7); relies on the standard CFD model with no verification of walk or amplitude dependence for the three scintillators.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Positron Annihilation Lifetime Spectroscopy Using Fast Scintillators and Digital Electronics." pith.science (2026). https://pith.science/paper/E2RD2JZI

@misc{pith2026190804434,
  author       = {Pith},
  title        = {Pith review of: Positron Annihilation Lifetime Spectroscopy Using Fast Scintillators and Digital Electronics},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/E2RD2JZI}},
  note         = {Machine review of arXiv:1908.04434}
}
read the original abstract

Positron Annihilation Lifetime Spectroscopy (PALS) is a non-destructive radiological technique widely used in material science studies. PALS typically relies on an analog coincidence measurement setup and allows the estimate of the positron lifetime in a material sample under investigation. The positronium trapping at vacancies in the material results in an increased lifetime. In this work, we have developed and optimized a PALS experimental setup using organic scintillators, fast digitizers, and advanced pulse processing algorithms. We tested three pairs of different organic scintillation detectors: EJ-309 liquid, EJ-276 newly developed plastic, and BC-418 plastic, and optimized the data processing parameters for each pair separately. Our high-throughput data analysis method is based on single-pulse interpolation and a constant fraction discrimination (CFD) algorithm. The setup based on the BC-418 detector achieved the best time resolution of 198.3 +- 0.8 ps. We used such optimized setup to analyze two single-crystal quartz samples and found lifetimes of 156 +- 9 ps and 366 +- 22 ps, in good agreement with the characteristic time constants of this material. The proposed experimental set up achieve an excellent time resolution, which makes it possible to accurately characterize material vacancies by discriminating between the lifetimes of either the spin singlet or triplet states of positronium. The optimized data processing algorithms are relevant to all the applications where fast timing is important, such as nuclear medicine and radiation imaging.

Figures

Figures reproduced from arXiv: 1908.04434 by the authors.

Figure 1
Figure 1. Detector time resolution measurement setup. [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Schematic diagram of the PALS measurement setup using two BC-418 detectors [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. An example pulse before and after interpolation. [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: The time difference distribution before interpolation and after interpolation. Mea [PITH_FULL_IMAGE:figures/full_fig_p008_4.png]
Figure 6
Figure 6. Figure 6: Optimized time resolution We can further reduce the FWHM by rejecting the low energy pulses. These pulses have small amplitudes and the sampling values could be easily affected 9 [PITH_FULL_IMAGE:figures/full_fig_p009_6.png]
Figure 7
Figure 7. Figure 7: shows the comparison of the positron lifetime spectrum in quartz and the distribution of arrival times obtained using the 60Co. The positron lifetime spectrum shows a longer tail due to longer lifetime, as expected. ● ● ●● ● ● ●● ● ●● ● ● ●● ● ●●● ● ● ●● ● ● ●● ● ● ● ●…
Figure 8
Figure 8. Figure 8: Fit the PAL spectrum of the single-crystal quartz sample. [PITH_FULL_IMAGE:figures/full_fig_p011_8.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

14 extracted references · 11 canonical work pages

  1. [1]

    Brandt, S

    W. Brandt, S. Berko, W. W. Walker, Positronium decay in molecular sub- stances, physical review 120 (4) (1960) 1289. doi:10.1103/PhysRev.120. 1289

    doi:10.1103/physrev.120 1960

  2. [2]

    Gidley, A

    D. Gidley, A. Rich, E. Sweetman, D. West, New precision measurements of the decay rates of singlet and triplet positronium, Physical Review Letters 49 (8) (1982) 525. doi:10.1103/PhysRevLett.49.525

    work page doi:10.1103/physrevlett.49.525 1982

  3. [3]

    Saito, T

    H. Saito, T. Hyodo, Direct measurement of the parapositronium lifetime in α- s i o 2, Physical review letters 90 (19) (2003) 193401. doi:10.1103/ PhysRevLett.90.193401

    work page 2003

  4. [4]

    Kansy, K

    J. Kansy, K. Mroczka, J. Dutkiewicz, Pals determination of defect density within friction stir welded joints of aluminium alloys, in: Journal of Physics: Conference Series, Vol. 265, IOP Publishing, 2011, p. 012010. doi:10. 1088/1742-6596/265/1/012010

    work page 2011

  5. [5]

    Gidley, W

    D. Gidley, W. Frieze, T. Dull, A. Yee, E. Ryan, H.-M. Ho, Positronium annihilation in mesoporous thin films, Physical Review B 60 (8) (1999) R5157. doi:10.1103/PhysRevB.60.R5157. 12

    work page doi:10.1103/physrevb.60.r5157 1999

  6. [6]

    J. D. Van Horn, F. Wu, G. Corsiglia, Y. C. Jean, Asymmetric positron in- teractions with chiral quartz crystals?, in: Defect and Diffusion Forum, Vol. 373, Trans Tech Publ, 2016, pp. 221–226. doi:10.4028/www.scientific. net/DDF.373.221

    doi:10.4028/www.scientific 2016

  7. [7]

    Hodges, B

    C. Hodges, B. McKee, W. Triftsh¨ auser, A. Stewart, Umklapp annihilation of positronium in crystals, Canadian Journal of Physics 50 (2) (1972) 103–

    work page 1972

  8. [8]

    Ryts¨ ol¨ a, J

    K. Ryts¨ ol¨ a, J. Nissil¨ a, J. Kokkonen, A. Laakso, R. Aavikko, K. Saarinen, Digital measurement of positron lifetime, Applied Surface Science 194 (1-4) (2002) 260–263. doi:10.1016/S0169-4332(02)00128-9

    work page doi:10.1016/s0169-4332(02)00128-9 2002

Show all 14 references
  1. [9]

    Beˇ cv´ aˇ r, J.ˇC´ ıˇ zek, I

    F. Beˇ cv´ aˇ r, J.ˇC´ ıˇ zek, I. Prochazka, High-resolution positron lifetime mea- surement using ultra fast digitizers acqiris dc211, Applied Surface Science 255 (1) (2008) 111–114. doi:10.1016/j.apsusc.2008.05.184

    2008 doi

  2. [10]

    S.pA., CoMPASS Multiparametric DAQ Software for Physics Applica- tions, CAEN S.pA

    C. S.pA., CoMPASS Multiparametric DAQ Software for Physics Applica- tions, CAEN S.pA

  3. [11]

    C. E. Shannon, Communication in the presence of noise, Proceedings of the IEEE 86 (2) (1998) 447–457. doi:10.1109/JPROC.1998.659497

    1998

  4. [12]

    W. K. Warburton, W. Hennig, New algorithms for improved digital pulse arrival timing with sub-gsps adcs, IEEE Transactions on Nuclear Science 64 (12) (2017) 2938–2950. doi:10.1109/TNS.2017.2766074

    2017

  5. [13]

    Steinberger, M

    W. Steinberger, M. Ruch, A. Di-Fulvio, S. Clarke, S. Pozzi, Timing perfor- mance of organic scintillators coupled to silicon photomultipliers, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 922 (20...

    2019 doi

  6. [14]

    Kansy, D

    J. Kansy, D. Giebel, Study of defect structure with new software for nu- merical analysis of pal spectra, in: Journal of Physics: Conference Series, 13 Vol. 265, IOP Publishing, 2011, p. 012030. doi:10.1088/1742-6596/265/ 1/012030. 14

    2011 doi

Pith tools

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