skip to main content


Title: A list-mode multi-energy window low-count SPECT reconstruction method for isotopes with multiple emission peaks
Abstract Background

Single-photon emission computed tomography (SPECT) provides a mechanism to perform absorbed-dose quantification tasks for$$\alpha$$α-particle radiopharmaceutical therapies ($$\alpha$$α-RPTs). However, quantitative SPECT for$$\alpha$$α-RPT is challenging due to the low number of detected counts, the complex emission spectrum, and other image-degrading artifacts. Towards addressing these challenges, we propose a low-count quantitative SPECT reconstruction method for isotopes with multiple emission peaks.

Methods

Given the low-count setting, it is important that the reconstruction method extracts the maximal possible information from each detected photon. Processing data over multiple energy windows and in list-mode (LM) format provide mechanisms to achieve that objective. Towards this goal, we propose a list-mode multi energy window (LM-MEW) ordered-subsets expectation–maximization-based SPECT reconstruction method that uses data from multiple energy windows in LM format and include the energy attribute of each detected photon. For computational efficiency, we developed a multi-GPU-based implementation of this method. The method was evaluated using 2-D SPECT simulation studies in a single-scatter setting conducted in the context of imaging [$$^{223}$$223Ra]RaCl$${_2}$$2, an FDA-approved RPT for metastatic prostate cancer.

Results

The proposed method yielded improved performance on the task of estimating activity uptake within known regions of interest in comparison to approaches that use a single energy window or use binned data. The improved performance was observed in terms of both accuracy and precision and for different sizes of the region of interest.

Conclusions

Results of our studies show that the use of multiple energy windows and processing data in LM format with the proposed LM-MEW method led to improved quantification performance in low-count SPECT of isotopes with multiple emission peaks. These results motivate further development and validation of the LM-MEW method for such imaging applications, including for$$\alpha$$α-RPT SPECT.

 
more » « less
Award ID(s):
2239707
NSF-PAR ID:
10424901
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
EJNMMI Physics
Volume:
10
Issue:
1
ISSN:
2197-7364
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Ultra-pure NaI(Tl) crystals are the key element for a model-independent verification of the long standing DAMA result and a powerful means to search for the annual modulation signature of dark matter interactions. The SABRE collaboration has been developing cutting-edge techniques for the reduction of intrinsic backgrounds over several years. In this paper we report the first characterization of a 3.4 kg crystal, named NaI-33, performed in an underground passive shielding setup at LNGS. NaI-33 has a record low$$^{39}$$39K contamination of 4.3 ± 0.2 ppb as determined by mass spectrometry. We measured a light yield of 11.1 ± 0.2 photoelectrons/keV and an energy resolution of 13.2% (FWHM/E) at 59.5 keV. We evaluated the activities of$$^{226}$$226Ra and$$^{228}$$228Th inside the crystal to be$$5.9\pm 0.6~\upmu $$5.9±0.6μBq/kg and$$1.6\pm 0.3~\upmu $$1.6±0.3μBq/kg, respectively, which would indicate a contamination from$$^{238}$$238U and$$^{232}$$232Th at part-per-trillion level. We measured an activity of 0.51 ± 0.02 mBq/kg due to$$^{210}$$210Pb out of equilibrium and a$$\alpha $$αquenching factor of 0.63 ± 0.01 at 5304 keV. We illustrate the analyses techniques developed to reject electronic noise in the lower part of the energy spectrum. A cut-based strategy and a multivariate approach indicated a rate, attributed to the intrinsic radioactivity of the crystal, of$$\sim $$1 count/day/kg/keV in the [5–20] keV region.

     
    more » « less
  2. Abstract

    We develop a two-timing perturbation analysis to provide quantitative insights on the existence of temporal ratchets in an exemplary system of a particle moving in a tank of fluid in response to an external vibration of the tank. We consider two-mode vibrations with angular frequencies$$\omega $$ωand$$\alpha \omega $$αω, where$$\alpha $$αis a rational number. If$$\alpha $$αis a ratio of odd and even integers (e.g.,$$\tfrac{2}{1},\,\tfrac{3}{2},\,\tfrac{4}{3}$$21,32,43), the system yields a net response: here, a nonzero time-average particle velocity. Our first-order perturbation solution predicts the existence of temporal ratchets for$$\alpha =2$$α=2. Furthermore, we demonstrate, for a reduced model, that the temporal ratcheting effect for$$\alpha =\tfrac{3}{2}$$α=32and$$\tfrac{4}{3}$$43appears at the third-order perturbation solution. More importantly, we find closed-form formulas for the magnitude and direction of the induced net velocities for these$$\alpha $$αvalues. On a broader scale, our methodology offers a new mathematical approach to study the complicated nature of temporal ratchets in physical systems.

    Graphic abstract 
    more » « less
  3. Abstract

    We introduce a family of Finsler metrics, called the$$L^p$$Lp-Fisher–Rao metrics$$F_p$$Fp, for$$p\in (1,\infty )$$p(1,), which generalizes the classical Fisher–Rao metric$$F_2$$F2, both on the space of densities$${\text {Dens}}_+(M)$$Dens+(M)and probability densities$${\text {Prob}}(M)$$Prob(M). We then study their relations to the Amari–C̆encov$$\alpha $$α-connections$$\nabla ^{(\alpha )}$$(α)from information geometry: on$${\text {Dens}}_+(M)$$Dens+(M), the geodesic equations of$$F_p$$Fpand$$\nabla ^{(\alpha )}$$(α)coincide, for$$p = 2/(1-\alpha )$$p=2/(1-α). Both are pullbacks of canonical constructions on$$L^p(M)$$Lp(M), in which geodesics are simply straight lines. In particular, this gives a new variational interpretation of$$\alpha $$α-geodesics as being energy minimizing curves. On$${\text {Prob}}(M)$$Prob(M), the$$F_p$$Fpand$$\nabla ^{(\alpha )}$$(α)geodesics can still be thought as pullbacks of natural operations on the unit sphere in$$L^p(M)$$Lp(M), but in this case they no longer coincide unless$$p=2$$p=2. Using this transformation, we solve the geodesic equation of the$$\alpha $$α-connection by showing that the geodesic are pullbacks of projections of straight lines onto the unit sphere, and they always cease to exists after finite time when they leave the positive part of the sphere. This unveils the geometric structure of solutions to the generalized Proudman–Johnson equations, and generalizes them to higher dimensions. In addition, we calculate the associate tensors of$$F_p$$Fp, and study their relation to$$\nabla ^{(\alpha )}$$(α).

     
    more » « less
  4. Abstract

    Complete theoretical understanding of the most complex superconductors requires a detailed knowledge of the symmetry of the superconducting energy-gap$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkα, for all momentakon the Fermi surface of every bandα. While there are a variety of techniques for determining$$|{\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha |$$Δkα, no general method existed to measure the signed values of$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkα. Recently, however, a technique based on phase-resolved visualization of superconducting quasiparticle interference (QPI) patterns, centered on a single non-magnetic impurity atom, was introduced. In principle, energy-resolved and phase-resolved Fourier analysis of these images identifies wavevectors connecting allk-space regions where$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkαhas the same or opposite sign. But use of a single isolated impurity atom, from whose precise location the spatial phase of the scattering interference pattern must be measured, is technically difficult. Here we introduce a generalization of this approach for use with multiple impurity atoms, and demonstrate its validity by comparing the$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkαit generates to the$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkαdetermined from single-atom scattering in FeSe where s±energy-gap symmetry is established. Finally, to exemplify utility, we use the multi-atom technique on LiFeAs and find scattering interference between the hole-like and electron-like pockets as predicted for$${\mathrm{{\Delta}}}_{\mathbf{k}}^\alpha$$Δkαof opposite sign.

     
    more » « less
  5. Abstract

    A type-II InAs/AlAs$$_{0.16}$$0.16Sb$$_{0.84}$$0.84multiple-quantum well sample is investigated for the photoexcited carrier dynamics as a function of excitation photon energy and lattice temperature. Time-resolved measurements are performed using a near-infrared pump pulse, with photon energies near to and above the band gap, probed with a terahertz probe pulse. The transient terahertz absorption is characterized by a multi-rise, multi-decay function that captures long-lived decay times and a metastable state for an excess-photon energy of$$>100$$>100meV. For sufficient excess-photon energy, excitation of the metastable state is followed by a transition to the long-lived states. Excitation dependence of the long-lived states map onto a nearly-direct band gap ($$E{_g}$$Eg) density of states with an Urbach tail below$$E{_g}$$Eg. As temperature increases, the long-lived decay times increase$$<Eg, due to the increased phonon interaction of the unintentional defect states, and by phonon stabilization of the hot carriers$$>E{_g}$$>Eg. Additionally, Auger (and/or trap-assisted Auger) scattering above the onset of the plateau may also contribute to longer hot-carrier lifetimes. Meanwhile, the initial decay component shows strong dependence on excitation energy and temperature, reflecting the complicated initial transfer of energy between valence-band and defect states, indicating methods to further prolong hot carriers for technological applications.

     
    more » « less