skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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:
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. 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}$$ 2 1 , 3 2 , 4 3 ), 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}$$ α = 3 2 and$$\tfrac{4}{3}$$ 4 3 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 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
  3. 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
  4. 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
  5. Abstract Let$$\lambda $$ λ denote the Liouville function. We show that the logarithmic mean of$$\lambda (\lfloor \alpha _1n\rfloor )\lambda (\lfloor \alpha _2n\rfloor )$$ λ ( α 1 n ) λ ( α 2 n ) is 0 whenever$$\alpha _1,\alpha _2$$ α 1 , α 2 are positive reals with$$\alpha _1/\alpha _2$$ α 1 / α 2 irrational. We also show that for$$k\geqslant 3$$ k 3 the logarithmic mean of$$\lambda (\lfloor \alpha _1n\rfloor )\cdots \lambda (\lfloor \alpha _kn\rfloor )$$ λ ( α 1 n ) λ ( α k n ) has some nontrivial amount of cancellation, under certain rational independence assumptions on the real numbers$$\alpha _i.$$ α i . Our results for the Liouville function generalise to produce independence statements for general bounded real-valued multiplicative functions evaluated at Beatty sequences. These results answer the two-point case of a conjecture of Frantzikinakis (and provide some progress on the higher order cases), generalising a recent result of Crnčević–Hernández–Rizk–Sereesuchart–Tao. As an ingredient in our proofs, we establish bounds for the logarithmic correlations of the Liouville function along Bohr sets. 
    more » « less