Title: A list-mode multi-energy window low-count SPECT reconstruction method for isotopes with multiple emission peaks
Abstract BackgroundSingle-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. MethodsGiven 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}$$ 223 Ra]RaCl$${_2}$$ 2 , an FDA-approved RPT for metastatic prostate cancer. ResultsThe 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. ConclusionsResults 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
PAR ID:
10482419
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Pubmed
Date Published:
Journal Name:
EJNMMI Physics
Volume:
10
Issue:
1
ISSN:
2197-7364
Subject(s) / Keyword(s):
List mode SPECT alpha-particle therapy Quantitative SPECT Low-count SPECT
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. AbstractWe 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}$$ α= 32 and$$\tfrac{4}{3}$$ 43 appears 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
  2. Abstract By the Aharonov–Casher theorem, the Pauli operatorPhas no zero eigenvalue when the normalized magnetic flux$$\alpha $$ α satisfies$$|\alpha |<1$$ |α|<1 , but it does have a zero energy resonance. We prove that in this case a Lieb–Thirring inequality for the$$\gamma $$ γ -th moment of the eigenvalues of$$P+V$$ P+V is valid under the optimal restrictions$$\gamma \ge |\alpha |$$ γ≥|α| and$$\gamma >0$$ γ>0 . Besides the usual semiclassical integral, the right side of our inequality involves an integral where the zero energy resonance state appears explicitly. Our inequality improves earlier works that were restricted to moments of order$$\gamma \ge 1$$ γ≥1 . 
    more » « less
  3. Abstract Low-lying states in$$^{54}$$ 54 Cr have been investigated via the$$\alpha $$ α -transfer reaction$$^{50}$$ 50 Ti($$^{7}$$ 7 Li,t) at a bombarding energy of 20 MeV. The exclusive$$\alpha $$ α -transfer channel is separated from other reaction channels through the appropriate energy gate on the complementary particle, triton. Levels of$$^{54}$$ 54 Cr populated exclusively by the$$\alpha $$ α -transfer process could be identified up to$$\approx $$ ≈ 5 MeV excitation energy and angular momentum up to$$(8)^{+}$$ (8)+ , by identifying the corresponding known$$\gamma $$ γ -rays. These include multiple low-lying non-yrast 2$$^+$$ + and 4$$^+$$ + states, which would otherwise be unfavorable via fusion evaporation reactions. The feeding-subtracted$$\gamma $$ γ -ray yields have been extracted to estimate the population of various excited states through the transfer process. The measured integrated transfer cross sections for all the observed yrast and non-yrast states are compared with Coupled Channels calculations usingfrescoto extract the$$\alpha $$ α +$$^{50}$$ 50 Ti core spectroscopic factors. For the yrast states, a higher$$\alpha $$ α +core overlap is seen for the$$2^+$$ 2+ and$$4^+$$ 4+ states, while it is seen to be less favorable for the$$6^+$$ 6+ and$$(8)^+$$ (8)+ states when$$\alpha $$ α -transfer is considered to occur predominantly as a direct one-step process to the$$^{50}$$ 50 Ti core ground state. The yrast$$2^+$$ 2+ , and$$4^+$$ 4+ states are predominantly populated by single-step transfer, while for the states with spin$$\ge $$ ≥ 5, the possibility of core excitation followed by$$\alpha $$ α -transfer shows a larger$$\alpha $$ α -core overlap. For the non-yrast$$0^+$$ 0+ ,$$2^+$$ 2+ , and$$4^+$$ 4+ states, single-step transfer shows moderate to small$$\alpha $$ α -core overlap. No higher spin non-yrast states are observed. 
    more » « less
  4. Abstract Large-scale organic liquid scintillator detectors are highly efficient in the detection of MeV-scale electron antineutrinos. These signal events can be detected through inverse beta decay on protons, which produce a positron accompanied by a neutron. A noteworthy background for antineutrinos coming from nuclear power reactors and from the depths of the Earth (geoneutrinos) is generated by ($$\alpha ,\,n$$ α,n ) reactions. In organic liquid scintillator detectors,$$\alpha $$ α particles emitted from intrinsic contaminants such as$$^{238}$$ 238 U,$$^{232}$$ 232 Th, and$$^{210}$$ 210 Pb/$$^{210}$$ 210 Po, can be captured on$$^{13}$$ 13 C nuclei, followed by the emission of a MeV-scale neutron. Three distinct interaction mechanisms can produce prompt energy depositions preceding the delayed neutron capture, leading to a pair of events correlated in space and time within the detector. Thus, ($$\alpha ,\,n$$ α,n ) reactions represent an indistinguishable background in liquid scintillator-based antineutrino detectors, where their expected rate and energy spectrum are typically evaluated via Monte Carlo simulations. This work presents results from the open-source SaG4n software, used to calculate the expected energy depositions from the neutron and any associated de-excitation products. Also simulated is a detailed detector response to these interactions, using a dedicated Geant4-based simulation software from the JUNO experiment. An expected measurable$$^{13}$$ 13 C$$(\alpha ,\,n)^{16}$$ (α,n) 16 O event rate and reconstructed prompt energy spectrum with associated uncertainties, are presented in the context of JUNO, however, the methods and results are applicable and relevant to other organic liquid scintillator neutrino detectors. 
    more » « less
  5. Abstract The evolutionary path of massive stars begins at helium burning. Energy production for this phase of stellar evolution is dominated by the reaction path 3$$\alpha \rightarrow ^{12}$$ α→12 C$$(\alpha ,\gamma )^{16}$$ (α,γ)16 O and also determines the ratio of$$^{12}$$ 12 C/$$^{16}$$ 16 O in the stellar core. This ratio then sets the evolutionary trajectory as the star evolves towards a white dwarf, neutron star or black hole. Although the reaction rate of the 3$$\alpha $$ α process is relatively well known, since it proceeds mainly through a single narrow resonance in$$^{12}$$ 12 C, that of the$$^{12}$$ 12 C$$(\alpha ,\gamma )^{16}$$ (α,γ)16 O reaction remains uncertain since it is the result of a more difficult to pin down, slowly-varying, portion of the cross section over a strong interference region between the high-energy tails of subthreshold resonances, the low-energy tails of higher-energy broad resonances and direct capture. Experimental measurements of this cross section require herculean efforts, since even at higher energies the cross section remains small and large background sources are often present that require the use of very sensitive experimental methods. Since the$$^{12}$$ 12 C$$(\alpha ,\gamma )^{16}$$ (α,γ)16 O reaction has such a strong influence on many different stellar objects, it is also interesting to try to back calculate the required rate needed to match astrophysical observations. This has become increasingly tempting, as the accuracy and precision of observational data has been steadily improving. Yet, the pitfall to this approach lies in the intermediary steps of modeling, where other uncertainties needed to model a star’s internal behavior remain highly uncertain. 
    more » « less