-
An analytic, moment-based method to estimate orthopositronium lifetimes in positron annihilation lifetime spectroscopy measurements
Authors:
Lucas Berens,
Isaac Hsu,
Chin-Tu Chen,
Howard Halpern,
Chien-Min Kao
Abstract:
The presence of tumor hypoxia is known to correlate with poor patient prognosis. Measurement of tissue oxygen concentration can be challenging, but recent advancements using positron annihilation lifetime spectroscopy (PALS) in three-dimensional positron emission tomography (PET) scans have shown promise for hypoxia detection. In this work, a novel method for estimating the orthopositronium lifeti…
▽ More
The presence of tumor hypoxia is known to correlate with poor patient prognosis. Measurement of tissue oxygen concentration can be challenging, but recent advancements using positron annihilation lifetime spectroscopy (PALS) in three-dimensional positron emission tomography (PET) scans have shown promise for hypoxia detection. In this work, a novel method for estimating the orthopositronium lifetime in PALS is presented. This method is analytical and uses moments of the time-difference histogram from photon arrival times. For sufficient statistical power, the method produces monotonic, stable estimates. For cases with a lower number of photon counts, the method was characterized and solutions are presented to correct for bias and estimation variability.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
Enhancing Positronium Lifetime Imaging through Two-Component Reconstruction in Time-of-Flight Positron Emission Tomography
Authors:
Zhuo Chen,
Chien-Min Kao,
Hsin-Hsiung Huang,
Lingling An
Abstract:
Positron Emission Tomography (PET) is a crucial tool in medical imaging, particularly for diagnosing diseases like cancer and Alzheimer's. The advent of Positronium Lifetime Imaging (PLI) has opened new avenues for assessing the tissue micro-environment, which is vital for early-stage disease detection. In this study, we introduce a two-component reconstruction model for PLI in Time-of-Flight (TOF…
▽ More
Positron Emission Tomography (PET) is a crucial tool in medical imaging, particularly for diagnosing diseases like cancer and Alzheimer's. The advent of Positronium Lifetime Imaging (PLI) has opened new avenues for assessing the tissue micro-environment, which is vital for early-stage disease detection. In this study, we introduce a two-component reconstruction model for PLI in Time-of-Flight (TOF) PET, incorporating both ortho-positronium and para-positronium decays. Our model enhances the accuracy of positronium imaging by providing a more detailed representation of the tissue environment. We conducted simulation studies to evaluate the performance of our model and compared it with existing single-component models. The results demonstrate the superiority of the two-component model in capturing the intricacies of the tissue micro-environment, thus paving the way for more precise and informative PET diagnostics.
△ Less
Submitted 22 March, 2024;
originally announced March 2024.
-
Simultaneous activity and attenuation estimation in TOF-PET with TV-constrained nonconvex optimization
Authors:
Zhimei Ren,
Emil Y. Sidky,
Rina Foygel Barber,
Chien-Min Kao,
Xiaochuan Pan
Abstract:
An alternating direction method of multipliers (ADMM) framework is developed for nonsmooth biconvex optimization for inverse problems in imaging. In particular, the simultaneous estimation of activity and attenuation (SAA) problem in time-of-flight positron emission tomography (TOF-PET) has such a structure when maximum likelihood estimation (MLE) is employed. The ADMM framework is applied to MLE…
▽ More
An alternating direction method of multipliers (ADMM) framework is developed for nonsmooth biconvex optimization for inverse problems in imaging. In particular, the simultaneous estimation of activity and attenuation (SAA) problem in time-of-flight positron emission tomography (TOF-PET) has such a structure when maximum likelihood estimation (MLE) is employed. The ADMM framework is applied to MLE for SAA in TOF-PET, resulting in the ADMM-SAA algorithm. This algorithm is extended by imposing total variation (TV) constraints on both the activity and attenuation map, resulting in the ADMM-TVSAA algorithm. The performance of this algorithm is illustrated using the penalized maximum likelihood activity and attenuation estimation (P-MLAA) algorithm as a reference. Additional results on step-size tuning and on the use of unconstrained ADMM-SAA are presented in the previous arXiv submission: arXiv:2303.17042v1.
△ Less
Submitted 9 February, 2024; v1 submitted 29 March, 2023;
originally announced March 2023.
-
A statistical reconstruction algorithm for positronium lifetime imaging using time-of-flight positron emission tomography
Authors:
Hsin-Hsiung Huang,
Zheyuan Zhu,
Slun Booppasiri,
Zhuo Chen,
Shuo Pang,
Chien-Min Kao
Abstract:
Positron emission tomography (PET) is an important modality for diagnosing diseases such as cancer and Alzheimer's disease, capable of revealing the uptake of radiolabeled molecules that target specific pathological markers of the diseases. Recently, positronium lifetime imaging (PLI) that adds to traditional PET the ability to explore properties of the tissue microenvironment beyond tracer uptake…
▽ More
Positron emission tomography (PET) is an important modality for diagnosing diseases such as cancer and Alzheimer's disease, capable of revealing the uptake of radiolabeled molecules that target specific pathological markers of the diseases. Recently, positronium lifetime imaging (PLI) that adds to traditional PET the ability to explore properties of the tissue microenvironment beyond tracer uptake has been demonstrated with time-of-flight (TOF) PET and the use of non-pure positron emitters. However, achieving accurate reconstruction of lifetime images from data acquired by systems having a finite TOF resolution still presents a challenge. This paper focuses on the two-dimensional PLI, introducing a maximum likelihood estimation (MLE) method that employs an exponentially modified Gaussian (EMG) probability distribution that describes the positronium lifetime data produced by TOF PET. We evaluate the performance of our EMG-based MLE method against \st{traditional} approaches using exponential likelihood functions and penalized surrogate methods. Results from computer-simulated data reveal that the proposed EMG-MLE method can yield quantitatively accurate lifetime images. We also demonstrate that the proposed MLE formulation can be extended to handle PLI data containing multiple positron populations.
△ Less
Submitted 9 May, 2024; v1 submitted 13 June, 2022;
originally announced June 2022.
-
Performance portable ice-sheet modeling with MALI
Authors:
Jerry Watkins,
Max Carlson,
Kyle Shan,
Irina Tezaur,
Mauro Perego,
Luca Bertagna,
Carolyn Kao,
Matthew J. Hoffman,
Stephen F. Price
Abstract:
High resolution simulations of polar ice-sheets play a crucial role in the ongoing effort to develop more accurate and reliable Earth-system models for probabilistic sea-level projections. These simulations often require a massive amount of memory and computation from large supercomputing clusters to provide sufficient accuracy and resolution. The latest exascale machines poised to come online con…
▽ More
High resolution simulations of polar ice-sheets play a crucial role in the ongoing effort to develop more accurate and reliable Earth-system models for probabilistic sea-level projections. These simulations often require a massive amount of memory and computation from large supercomputing clusters to provide sufficient accuracy and resolution. The latest exascale machines poised to come online contain a diverse set of computing architectures. In an effort to avoid architecture specific programming and maintain productivity across platforms, the ice-sheet modeling code known as MALI uses high level abstractions to integrate Trilinos libraries and the Kokkos programming model for performance portable code across a variety of different architectures. In this paper, we analyze the performance portable features of MALI via a performance analysis on current CPU-based and GPU-based supercomputers. The analysis highlights performance portable improvements made in finite element assembly and multigrid preconditioning within MALI with speedups between 1.26-1.82x across CPU and GPU architectures but also identifies the need to further improve performance in software coupling and preconditioning on GPUs. We also perform a weak scalability study and show that simulations on GPU-based machines perform 1.24-1.92x faster when utilizing the GPUs. The best performance is found in finite element assembly which achieved a speedup of up to 8.65x and a weak scaling efficiency of 82.9% with GPUs. We additionally describe an automated performance testing framework developed for this code base using a changepoint detection method. The framework is used to make actionable decisions about performance within MALI. We provide several concrete examples of scenarios in which the framework has identified performance regressions, improvements, and algorithm differences over the course of two years of development.
△ Less
Submitted 8 April, 2022;
originally announced April 2022.
-
Development of a PET/EPRI combined imaging system for assessing tumor hypoxia
Authors:
Heejong Kim,
Boris Epel,
Subramanian Sundramoorthy,
Hsiu-Ming Tsai,
Eugene Barth,
Inna Gertsenshteyn,
Howard Halpern,
Yuexuan Hua,
Qingguo Xie,
Chin-Tu Chen,
Chien-Min Kao
Abstract:
Precise quantitative delineation of tumor hypoxia is essential in radiation therapy treatment planning to improve the treatment efficacy by targeting hypoxic sub-volumes. We developed a combined imaging system of positron emission tomography (PET) and electron para-magnetic resonance imaging (EPRI) of molecular oxygen to investigate the accuracy of PET imaging in assessing tumor hypoxia. The PET/E…
▽ More
Precise quantitative delineation of tumor hypoxia is essential in radiation therapy treatment planning to improve the treatment efficacy by targeting hypoxic sub-volumes. We developed a combined imaging system of positron emission tomography (PET) and electron para-magnetic resonance imaging (EPRI) of molecular oxygen to investigate the accuracy of PET imaging in assessing tumor hypoxia. The PET/EPRI combined imaging system aims to use EPRI to precisely measure the oxygen partial pressure in tissues. This will evaluate the validity of PET hypoxic tumor imaging by (near) simultaneously acquired EPRI as ground truth. The combined imaging system was constructed by integrating a small animal PET scanner (inner ring diameter 62 mm and axial field of view 25.6 mm) and an EPRI subsystem (field strength 25 mT and resonant frequency 700 MHz). The compatibility between the PET and EPRI subsystems were tested with both phantom and animal imaging. Hypoxic imaging on a tumor mouse model using $^{18}$F-fluoromisonidazole radio-tracer was conducted with the developed PET/EPRI system. We report the development and initial imaging results obtained from the PET/EPRI combined imaging system.
△ Less
Submitted 27 January, 2021;
originally announced January 2021.
-
Disparate exciton-phonon couplings for zone center and boundary phonons in solid-state graphite
Authors:
Xuefei Feng,
Shawn Sallis,
Yu-Cheng Shao,
Ruimin Qiao,
Yi-Sheng Liu,
Li Cheng Kao,
Anton Tremsin,
Zahid Hussain,
Wanli Yang,
Jinghua Guo,
Yi-De Chuang
Abstract:
The exciton-phonon coupling in highly oriented pyrolytic graphite is studied using resonant inelastic X-ray scattering (RIXS) spectroscopy. With ~ 70 meV energy resolution, multiple low energy excitations associated with coupling to phonons can be clearly resolved in RIXS spectra. Using resonance dependence and the closed form for RIXS cross-section without considering the intermediate state mixin…
▽ More
The exciton-phonon coupling in highly oriented pyrolytic graphite is studied using resonant inelastic X-ray scattering (RIXS) spectroscopy. With ~ 70 meV energy resolution, multiple low energy excitations associated with coupling to phonons can be clearly resolved in RIXS spectra. Using resonance dependence and the closed form for RIXS cross-section without considering the intermediate state mixing of phonon modes, the dimensionless coupling constant g is determined to be 5 and 0.4, corresponding to the coupling strength of 0.42 eV +/- 40 meV and 0.21 eV +/- 30 meV, for zone center and boundary phonons respectively. The reduced g value for zone-boundary phonon may be related to its double resonance nature.
△ Less
Submitted 5 January, 2020;
originally announced January 2020.
-
Lifetime and population of the 2S state in muonic protium and deuterium
Authors:
Marc Diepold,
Fernando D. Amaro,
Aldo Antognini,
François Biraben,
João M. R. Cardoso,
Daniel S. Covita,
Andreas Dax,
Satish Dhawan,
Luis M. P. Fernandes,
Adolf Giesen,
Andrea L. Gouvea,
Thomas Graf,
Theodor W. Hänsch,
Paul Indelicato,
Lucile Julien,
Cheng-Yang Kao,
Paul Knowles,
Franz Kottmann,
Eric-Olivier Le Bigot,
Yi-Wei Liu,
José A. M. Lopes,
Livia Ludhova,
Cristina M. B. Monteiro,
Françoise Mulhauser,
Tobias Nebel
, et al. (10 additional authors not shown)
Abstract:
Radiative deexcitation (RD) of the metastable 2S state of muonic protium and deuterium atoms has been observed. In muonic protium, we improve the precision on lifetime and population (formation probability) values for the short-lived μp(2S) component, and give an upper limit for RD of long-lived μp(2S) atoms. In muonic deuterium at 1 hPa, 3.1 +-0.3 % of all stopped muons form μd(2S) atoms. The sho…
▽ More
Radiative deexcitation (RD) of the metastable 2S state of muonic protium and deuterium atoms has been observed. In muonic protium, we improve the precision on lifetime and population (formation probability) values for the short-lived μp(2S) component, and give an upper limit for RD of long-lived μp(2S) atoms. In muonic deuterium at 1 hPa, 3.1 +-0.3 % of all stopped muons form μd(2S) atoms. The short-lived 2S component has a population of 1.35 +0.57 -0.33 % and a lifetime of τ_short(μd) = 138 +32 -34 ns. We see evidence for RD of long-lived μd(2S) with a lifetime of τ_long(μd) = 1.15 +0.75 -0.53 μs. This is interpreted as formation and decay of excited muonic molecules.
△ Less
Submitted 10 October, 2013; v1 submitted 25 July, 2013;
originally announced July 2013.
-
Detecting high density ultracold molecules using atom-molecule collision
Authors:
Jun-Ren Chen,
Cheng-Yang Kao,
Hung-Bin Chen,
Yi-Wei Liu
Abstract:
Utilizing single-photon photoassociation, we have achieved ultracold rubidium molecules with a high number density that provides a new efficient approach toward molecular quantum degeneracy. A new detection mechanism for ultracold molecule utilizing the inelastic atom-molecule collision is demonstrated. The resonant coupling effect on the formation of the ${\rm X^1Σ^+_g}$ ground state…
▽ More
Utilizing single-photon photoassociation, we have achieved ultracold rubidium molecules with a high number density that provides a new efficient approach toward molecular quantum degeneracy. A new detection mechanism for ultracold molecule utilizing the inelastic atom-molecule collision is demonstrated. The resonant coupling effect on the formation of the ${\rm X^1Σ^+_g}$ ground state ${\rm ^{85}Rb_2}$ allows for a sufficient number of more deeply bound ultracold molecules, which induced an additional trap loss and heating of the co-existing atoms owing to the inelastic atom-molecule collision. Therefore, after photoassociation process, the ultracold molecules can be investigated using the absorption image of the ultracold rubidium atoms mixed with the molecules in a crossed optical dipole trap. The existence of the ultracold molecules was then verified, and the amount of the accumulated molecules was measured. This method is to detect the final produced ultracold molecules, and hence distinct from the conventional trap loss experiments, which is used to study the association resonance. It is composed of measurements of the time evolution of atomic cloud and a decay model, by which the number density of the ultracold ${\rm ^{85}Rb_2}$ molecules in the optical trap was estimated to be ${\rm > 5.2\times10^{11} cm^{-3}}$.
△ Less
Submitted 7 March, 2013; v1 submitted 28 November, 2012;
originally announced November 2012.
-
Patterned coating by suspensions
Authors:
Justin C. T. Kao,
A. E. Hosoi
Abstract:
This is a fluid dynamics video submitted to the Gallery of Fluid Motion of the 63rd Annual Meeting of the American Physical Society Division of Fluid Dynamics.
We examine the Landau-Levich flow of a suspension of neutrally buoyant particles within a rotating glass cylinder. This differs from previously investigated systems in that we deal with macroscopic particles which immediately exceed the t…
▽ More
This is a fluid dynamics video submitted to the Gallery of Fluid Motion of the 63rd Annual Meeting of the American Physical Society Division of Fluid Dynamics.
We examine the Landau-Levich flow of a suspension of neutrally buoyant particles within a rotating glass cylinder. This differs from previously investigated systems in that we deal with macroscopic particles which immediately exceed the thickness of the deposited liquid film. Thus, the presence of particles on the liquid-coated wall creates excess surface area, leading to a capillary attraction between particles. In turn, this capillary attraction induces phase separation of the suspension coating into regions of either dense particle clustering or clear liquid.
△ Less
Submitted 15 October, 2010;
originally announced October 2010.
-
Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT
Authors:
Emil Y. Sidky,
Chien-Min Kao,
Xiaochuan Pan
Abstract:
In practical applications of tomographic imaging, there are often challenges for image reconstruction due to under-sampling and insufficient data. In computed tomography (CT), for example, image reconstruction from few views would enable rapid scanning with a reduced x-ray dose delivered to the patient. Limited-angle problems are also of practical significance in CT. In this work, we develop and…
▽ More
In practical applications of tomographic imaging, there are often challenges for image reconstruction due to under-sampling and insufficient data. In computed tomography (CT), for example, image reconstruction from few views would enable rapid scanning with a reduced x-ray dose delivered to the patient. Limited-angle problems are also of practical significance in CT. In this work, we develop and investigate an iterative image reconstruction algorithm based on the minimization of the image total variation (TV) that applies to divergent-beam CT. Numerical demonstrations of our TV algorithm are performed with various insufficient data problems in fan-beam CT. The TV algorithm can be generalized to cone-beam CT as well as other tomographic imaging modalities.
△ Less
Submitted 28 April, 2009;
originally announced April 2009.