Characterization of176Lu background in LSO-based PET scanners

Quality control performance evaluation of uMI550 digital PET-CT by utilising LYSO crystal natural radioactivity of 176Lu

PURPOSE

Quality Control (QC) of the Positron Emission Tomography-Computed Tomography (PET-CT) system must be performed prior to the PET-CT acquisition to ensure the reproducibility as per the manufacturer recommendation. In this study we have evaluated the performance of daily PET QC test by utilising lutetium yttrium oxyorthosilicate (LYSO) scintillation crystal natural radioactivity of 176Lu as a source of radiation to perform the PET uCare.iQC with uMI550 digital PET-CT system. This was also compared with existing radioactive external source-based QC test with other manufacturer PET-CT systems.

METHOD

This radioactive source free daily QC study was performed on uMI550 digital PET-CT system. The daily QC data report was captured and interpreted. This PET-CT system has unique feature that utilises the inherent property of LYSO crystal that is 176Lu with natural radioactivity abundance of 2.6%. The Lutetium-176 (176Lu) radioactivity is used to perform the daily QC in PET in place of external radioactive source of Germanium-68 (68Ge). This feature work automatically in preschedule manner to complete the daily QC at preset time in the morning and system get ready after the QC test.

RESULTS

Over 120 automatic PET daily uCare.iQC test were performed. The daily PET QC test was prescheduling setup for 6:00 am in every morning. No failure on daily QC test were observed. The QC parameters and system parameters consistency was observed.

CONCLUSION

It was concluded that the daily PET QC can be performed by utilising LYSO crystal inherent natural radioactivity of 176Lu as a source of radiation to perform the test as replacement of external 68Ge radioactive source.

IMPLICATION FOR PRACTICE

PET-CT daily QC by utilizing the 176Lu radioactivity of LYSO crystal results in reducing the radiation exposure to operation staff and reducing operational cost by elimination 68Ge shield source Phantom.

Ultrasensitive and multiplexed tracking of single cells using whole-body PET/CT

In vivo molecular imaging tools are crucially important for elucidating how cells move through complex biological systems; however, achieving single-cell sensitivity over the entire body remains challenging. Here, we report a highly sensitive and multiplexed approach for tracking upward of 20 single cells simultaneously in the same subject using positron emission tomography (PET). The method relies on a statistical tracking algorithm (PEPT-EM) to achieve a sensitivity of 4 becquerel per cell and a streamlined workflow to reliably label single cells with over 50 becquerel per cell of 18F-fluorodeoxyglucose (FDG). To demonstrate the potential of the method, we tracked the fate of more than 70 melanoma cells after intracardiac injection and found they primarily arrested in the small capillaries of the pulmonary, musculoskeletal, and digestive organ systems. This study bolsters the evolving potential of PET in offering unmatched insights into the earliest phases of cell trafficking in physiological and pathological processes and in cell-based therapies.

Efficient and multiplexed tracking of single cells using whole-body PET/CT

In vivo molecular imaging tools are crucially important for elucidating how cells move through complex biological systems, however, achieving single-cell sensitivity over the entire body remains challenging. Here, we report a highly sensitive and multiplexed approach for tracking upwards of 20 single cells simultaneously in the same subject using positron emission tomography (PET). The method relies on a new tracking algorithm (PEPT-EM) to push the cellular detection threshold to below 4 Bq/cell, and a streamlined workflow to reliably label single cells with over 50 Bq/cell of 18F-fluorodeoxyglucose (FDG). To demonstrate the potential of method, we tracked the fate of over 70 melanoma cells after intracardiac injection and found they primarily arrested in the small capillaries of the pulmonary, musculoskeletal, and digestive organ systems. This study bolsters the evolving potential of PET in offering unmatched insights into the earliest phases of cell trafficking in physiological and pathological processes and in cell-based therapies.

Technical opportunities and challenges in developing total-body PET scanners for mice and rats

Positron emission tomography (PET) is the most sensitive in vivo molecular imaging technique available. Small animal PET has been widely used in studying pharmaceutical biodistribution and disease progression over time by imaging a wide range of biological processes. However, it remains true that almost all small animal PET studies using mouse or rat as preclinical models are either limited by the spatial resolution or the sensitivity (especially for dynamic studies), or both, reducing the quantitative accuracy and quantitative precision of the results. Total-body small animal PET scanners, which have axial lengths longer than the nose-to-anus length of the mouse/rat and can provide high sensitivity across the entire body of mouse/rat, can realize new opportunities for small animal PET. This article aims to discuss the technical opportunities and challenges in developing total-body small animal PET scanners for mice and rats.

Validation of a computational chain from PET Monte Carlo simulations to reconstructed images

The study aimed to create a pipeline from Monte Carlo simulated projections of a Gate PET system to reconstructed images. The PET system was modelled after the GE Discovery MI (DMI) PET/CT, and the simulated projections were reconstructed with the stand-alone reconstruction software CASToR. Attenuation correction, normalisation calibration, random estimation, and scatter estimation for the simulations were computed with in-house programs. The pipeline was compared in both projection and image space with data acquired on a clinical DMI and reconstructed with GE's off-line PET reconstruction software (PET Toolbox) and CASToR. The simulated and measured data were compared for the number of prompt coincidences, scatter fraction, contrast recovery coefficient (CRC), signal-to-noise ratio (SNR), background variability, residual lung error, and image profiles. A slight discrepancy was noted in the projection space, but good agreements were generally achieved in image space between simulated and measured data. The CRC was found to be 81 % for Gate – CASToR, 84 % for GE – CASToR, and 84 % for GE - PET Toolbox for the largest sphere of the NEMA image quality (IQ) phantom, and the SNR was found to be 98 for Gate – CASToR, 91 for GE – CASToR, and 93 for GE – PET Toolbox. Profiles drawn over the spheres for the NEMA IQ phantom and the Data Spectrum (DS) phantom show a good match between measurement and simulation. The results indicate feasibility to utilise the pipeline as a tool for off-line simulation-based studies. A complete pipeline introduces possibilities to study the impact of single parameters in the whole chain from simulation to reconstructed images.

Precision dosimetry in yttrium-90 radioembolization through CT imaging of radiopaque microspheres in a rabbit liver model

PURPOSE

To perform precision dosimetry in yttrium-90 radioembolization through CT imaging of radiopaque microspheres in a rabbit liver model and to compare extracted dose metrics to those produced from conventional PET-based dosimetry.

MATERIALS AND METHODS

A CT calibration phantom was designed containing posts with nominal microsphere concentrations of 0.5 mg/mL, 5.0 mg/mL, and 25.0 mg/mL. The mean Hounsfield unit was extracted from the post volumes to generate a calibration curve to relate Hounsfield units to microsphere concentration. A nominal bolus of 40 mg of microspheres was administered to the livers of eight rabbits, followed by PET/CT imaging. A CT-based activity distribution was calculated through the application of the calibration curve to the CT liver volume. Post-treatment dosimetry was performed through the convolution of yttrium-90 dose-voxel kernels and the PET- and CT-based cumulated activity distributions. The mean dose to the liver in PET- and CT-based dose distributions was compared through linear regression, ANOVA, and Bland-Altman analysis.

RESULTS

A linear least-squares fit to the average Hounsfield unit and microsphere concentration data from the calibration phantom confirmed a strong correlation (r2 > 0.999) with a slope of 14.13 HU/mg/mL. A poor correlation was found between the mean dose derived from CT and PET (r2 = 0.374), while the ANOVA analysis revealed statistically significant differences (p < 10-12) between the MIRD-derived mean dose and the PET- and CT-derived mean dose. Bland-Altman analysis predicted an offset of 15.0 Gy between the mean dose in CT and PET. The dose within the liver was shown to be more heterogeneous in CT than in PET with an average coefficient of variation equal to 1.99 and 1.02, respectively.

CONCLUSION

The benefits of a CT-based approach to post-treatment dosimetry in yttrium-90 radioembolization include improved visualization of the dose distribution, reduced partial volume effects, a better representation of dose heterogeneity, and the mitigation of respiratory motion effects. Post-treatment CT imaging of radiopaque microspheres in yttrium-90 radioembolization provides the means to perform precision dosimetry and extract accurate dose metrics used to refine the understanding of the dose-response relationship, which could ultimately improve future patient outcomes.

Optimization of GFAG crystal surface treatment for SiPM based TOF PET detector

Coincidence timing resolution (CTR) is an important parameter in clinical positron emission tomography (PET) scanners to increase the signal-to-noise ratio of PET images by using time-of-flight (TOF) information. Lutetium (Lu) based scintillators are often used for TOF-PET systems. However, the self-radiation of Lu-based scintillators may influence the image quality for ultra-low activity PET imaging. Recently, a gadolinium fine aluminum gallate (Ce:GFAG) scintillation crystal that features a fast decay time (~55 ns) and no self-radiation was developed. The present study aimed at optimizing the GFAG crystal surface treatment to enhance both CTR and energy resolution (ER). The TOF-PET detector consisted of a GFAG crystal (3.0 × 3.0 × 20 mm3) and a SiPM with an effective area of 3.0 × 3.0 mm2. The timing and energy signals were extracted using a high-frequency SiPM readout circuit and then were digitized using a CAMAC DAQ system. The CTR and ER were evaluated with nine different crystal surface treatments such as partial saw-cut and chemical polishing and the 1-side saw-cut was the best choice among the treatments. The respective CTR and ER of 202±2 ps and 9.5±0.1% were obtained with the 1-side saw-cut; the other 5-side mechanically polished GFAG crystals had respective values which were 18 ps (9.0%) and 1.3% better than those of the all-side mechanically polished GFAG crystal. The chemically polished GFAG crystals also offered enhanced CTR and ER of about 17 ps (8.2%) and 2.1%, respectively, over the mechanically polished GFAG crystals.

Lutetium background radiation in total-body PET-A simulation study on opportunities and challenges in PET attenuation correction

The current generation of total-body positron emission tomography (PET) scanners offer significant sensitivity increase with an extended axial imaging extent. With the large volume of lutetium-based scintillation crystals that are used as detector elements in these scanners, there is an increased flux of background radiation originating from ¹⁷⁶Lu decay in the crystals and higher sensitivity for detecting it. Combined with the ability of scanning the entire body in a single bed position, this allows more effective utilization of the lutetium background as a transmission source for estimating 511 keV attenuation coefficients. In this study, utilization of the lutetium background radiation for attenuation correction in total-body PET was studied using Monte Carlo simulations of a 3D whole-body XCAT phantom in the uEXPLORER PET scanner, with particular focus on ultralow-dose PET scans that are now made possible with these scanners. Effects of an increased acceptance angle, reduced scan durations, and Compton scattering on PET quantification were studied. Furthermore, quantification accuracy of lutetium-based attenuation correction was compared for a 20-min scan of the whole body on the uEXPLORER, a one-meter-long, and a conventional 24-cm-long scanner. Quantification and lesion contrast were minimally affected in both long axial field-of-view scanners and in a whole-body 20-min scan, the mean bias in all analyzed organs of interest were within a ±10% range compared to ground-truth activity maps. Quantification was affected in certain organs, when scan duration was reduced to 5 min or a reduced acceptance angle of 17° was used. Analysis of the Compton scattered events suggests that implementing a scatter correction method for the transmission data will be required, and increasing the energy threshold from 250 keV to 290 keV can reduce the computational costs and data rates, with negligible effects on PET quantification. Finally, the current results can serve as groundwork for transferring lutetium-based attenuation correction into research and clinical practice.

A CT-less approach to quantitative PET imaging using the LSO intrinsic radiation for long-axial FOV PET scanners

Purpose

Long‐axial field‐of‐view (FOV) positron emission tomography (PET) scanners have gained a lot of interest in the recent years. Such scanners provide increased sensitivity and enable unique imaging opportunities that were not previously feasible. Benefiting from the high sensitivity of a long‐axial FOV PET scanner, we studied a computed tomography (CT)–less reconstruction algorithm for the Siemens Biograph Vision Quadra with an axial FOV of 106 cm.

Methods

In this work, the background radiation from radioisotope lutetium‐176 in the scintillators was used to create an initial estimate of the attenuation maps. Then, joint activity and attenuation reconstruction algorithms were used to create an improved attenuation map of the object. The final attenuation maps were then used to reconstruct quantitative PET images, which were compared against CT‐based PET images. The proposed method was evaluated on data from three patients who underwent a flurodeoxyglucouse PET scan.

Results

Segmentation of the PET images of the three studied patients showed an average quantitative error of 6.5%–8.3% across all studied organs when using attenuation maps from maximum likelihood estimation of attenuation and activity and 5.3%–6.6% when using attenuation maps from maximum likelihood estimation of activity and attenuation correction coefficients.

Conclusions

Benefiting from the background radiation of lutetium‐based scintillators, a quantitative CT‐less PET imaging technique was evaluated in this work.

Discrete Convolution-Based Energy Spectrum Configuring Method for the Analysis of the Intrinsic Radiation of 176Lu

There has been considerable interest in inorganic scintillators based on lutetium due to their favorable physical properties. Despite their advantages, lutetium-based scintillators could face issues because of the natural occurring radioisotope of ¹⁷⁶Lu that is contained in natural lutetium. In order to mitigate its potential shortcomings, previous works have studied to understand the energy spectrum of the intrinsic radiation of ¹⁷⁶Lu (IRL). However, few studies have focused on the various principal types of photon interactions with matter; in other words, only the full-energy peak according to the photoelectric effect or internal conversion have been considered for understanding the energy spectrum of IRL. Thus, the approach we have used in this study considers other principal types of photon interactions by convoluting each energy spectrum with combinations for generating the spectrum of the intrinsic radiation of ¹⁷⁶Lu. From the results, we confirm that the method provides good agreement with the experiment. A significant contribution of this study is the provision of a new approach to process energy spectra induced by mutually independent radiation interactions as a single spectrum.

Yttrium-90 TOF-PET-Based EUD Predicts Response Post Liver Radioembolizations Using Recommended Manufacturer FDG Reconstruction Parameters

Purpose

Explaining why ⁹⁰Y TOF-PET based equivalent uniform dose (EUD) using recommended manufacturer FDG reconstruction parameters has been shown to predict response.

Methods

The hot rods insert of a Jaszczak deluxe phantom was partially filled with a 2.65 GBq ⁹⁰Y - 300ml DTPA water solution resulting in a 100 Gy mean absorbed dose in the 6 sectors. A two bed 20min/position acquisition was performed on a 550ps- and on a 320ps- TOF-PET/CT and reconstructed with recommended manufacturer FDG reconstruction parameters, without and with additional filtering. The whole procedure was repeated on both PET after adding 300ml of water (50Gy setup). The phantom was acquired again after decay by a factor of 10 (5Gy setup), but with 200min per bed position. For comparison, the phantom was also acquired with ¹⁸F activity corresponding to a clinical FDG whole body acquisition.

Results

The 100Gy-setup provided a hot rod sectors image almost as good as the ¹⁸F phantom. However, despite acquisition time compensation, the 5Gy-setup provides much lower quality imaging. TOF-PET based sectors EUDs for the three large rod sectors agreed with the actual EUDs computed with a radiosensitivity of 0.021Gy^-1 well in the range observed in external beam radiotherapy (EBRT), i.e. 0.01-0.04Gy^-1. This agreement explains the reunification of the dose-response relationships of the glass and resin spheres in HCC using the TOF-PET based EUD. Additional filtering reduced the EUDs agreement quality.

Conclusions

Recommended manufacturer FDG reconstruction parameters are suitable in TOF-PET post ⁹⁰Y liver radioembolization for accurate tumour EUD computation. The present results rule out the use of low specific activity phantom studies to optimize reconstruction parameters.

Energy spectra due to the intrinsic radiation of LYSO/LSO scintillators for a wide range of crystal sizes

PURPOSE

Most detectors in current positron emission tomography (PET) scanners and prototypes use lutetium oxyorthosilicate (LSO) or lutetium yttrium oxyorthosilicate (LYSO) scintillators. The aim of this work is to provide a complete set of background energy spectra, due to the scintillator intrinsic radioactivity, for a wide range of crystal sizes.

METHODS

An analytical model, developed and validated in a previous work, was used to obtain the background energy spectra of square base cuboids of different dimensions. The model uses the photon absorption probabilities of the three gamma rays (88, 202, and 307 keV) emitted following the beta decay of ¹⁷⁶ Lu to ¹⁷⁶ Hf excited states. These probabilities were obtained for each crystal size considered in this work from Monte Carlo simulations using the PENELOPE code. The probabilities are then used to normalize and shift the beta spectrum to the corresponding energy value of the simultaneous detection of one, two, or three gamma rays in the scintillator. The simulated cuboids had side lengths of 5, 10, 20, 30, 40, 50, and 60 mm and crystal thickness T = 2.5, 5, 10, 15, and 20 mm. From these results a complete set of energy spectra, including intermediate dimensions, were obtained. In addition, LYSO and LSO were compared in terms of their analytical background energy spectra for two crystal sizes. The analytical spectra were convolved using a variable Gaussian kernel to account for the energy resolution of a typical detector. A parameterization of the photon absorption probabilities for each gamma ray energy as a function of the cuboid volume to surface area ratio was obtained.

RESULTS

A data set of L(Y)SO background energy spectra was obtained and is available for the reader as 2D histograms. The model accurately predicts the structure of the energy spectra including the relative peak and valley intensities. The data allow visualizing how the structure evolves with increasing crystal length and thickness. Lutetium yttrium oxyorthosilicate and LSO present very similar background energy spectra for the range of sizes studied in this work and therefore the data generated can be confidently used for both scintillator materials. The filtered spectra showed a variable shift in the main peaks, depending on crystal size, alerting that to achieve a correct detector calibration using the background spectrum is not straight forward and requires precise data analysis and measurements. In addition, we found that square base L(Y)SO cuboids with same volume to surface area ratio have background spectra with the same structure.

CONCLUSIONS

We present the energy spectra of L(Y)SO crystal of different sizes which will be very useful for industry and research groups developing and simulating detectors for positron imaging applications in terms of calibration, quality assurance, crystal maps, detector fine gain tuning, background reduction and other applications using the long-lived ¹⁷⁶ Lu source. We analyzed the data produced in this work and found that crystal cuboids with equal volume to surface area ratio produce the same background energy spectra, a conclusion that simplifies its calculation and clarifies why the same energy spectrum is observed under different experimental setups.

Experimental validation of the ANTS2 code for modelling optical photon transport in monolithic LYSO crystals

In this work the scintillation energy spectra originating from the background radioactivity from polished monolithic lutetium yttrium oxyorthosilicate coupled to position-sensitive silicon photomultipliers (SiPM) was studied using the open source Monte Carlo simulation package ANTS2. Two crystal sizes, fully and partially covering the photosensor area, three surface crystal wrappings (black, specular or diffuse) and the full signal formation process in the photosensor were considered. The simulation results were validated with experimental data acquired under the same geometric and detector operating conditions. In all cases ANTS2 simulated spectra have very good agreement with experimental results, reproducing the expected shape, with correct onset and end at 88 and 1190 keV, respectively, as well as sharp edges at the reference energies of 88, 88 + 202, 88 + 307 and 88 + 202 + 307 keV. The normalized root-mean square error between simulated and measured spectra varied between 4.3% and 10.4%.

Coincidence energy spectra due to the intrinsic radioactivity of LYSO scintillation crystals

Background

Lutetium oxyorthosilicate or lutetium yttrium oxyorthosilicate (LYSO) scintillation crystals used in most current PET scanner detectors contain ¹⁷⁶Lu, which decays by beta emission to excited states of ¹⁷⁶Hf accompanied by the emission of prompt gamma rays or internal conversion electrons. This intrinsic radioactivity can be self-detected in singles mode as a constant background signal that has an energy spectrum whose structure has been explained previously. In this work, we studied the energy spectrum due to the intrinsic radioactivity of LYSO scintillation crystals of two opposing detectors working in coincidence mode. The investigation included experimental data, Monte Carlo simulations and an analytical model.

Results

The structure of the energy spectrum was completely understood and is the result of the self-detection of beta particles from ¹⁷⁶Lu in one crystal and the detection of one or more prompt gamma rays detected in coincidence by the opposing crystal. The most probable coincidence detection involves the gamma rays of 202 and 307 keV, which result in two narrow photopeaks, superimposed on a continuous energy distribution due to the beta particle energy deposition. The relative intensities of the gamma ray peaks depend on crystal size and detector separation distance, as is explained by the analytical model and verified through the Monte Carlo simulations and experiments.

Conclusions

The analytical model used in this work accurately explains the general features of the coincidence energy spectrum due to the presence of ¹⁷⁶Lu in the scintillation crystals, as observed experimentally and with Monte Carlo simulations. This work will be useful to those research studies aimed at using the intrinsic radioactivity of LYSO crystals for transmission scans and detector calibration in coincidence mode.

Convenient method for improving the light output of scintillators by using buffer layers coated with photonic crystals

The low light-extraction efficiency of scintillators is due to total internal reflection and has led to the extensive use of photonic crystals to improve the light output. However, in some applications, photonic crystals cannot be fabricated directly on scintillators. Here, we demonstrate a promising method to improve the light output of scintillators by using a buffer layer coated with photonic crystals and then fixed to the scintillator. Through both numerical simulations and experiments, we investigate how the refractive indexes of the buffer layer and photonic crystal affect the light output from scintillators. The experimental results indicate that the light output of (Lu,Y)₂SiO₅:Ce scintillators is enhanced 1.9 times by using a sapphire buffer layer coated with an array of polystyrene nanospheres. This method can be used to improve the detection efficiency of radiation-detection systems when photonic crystals cannot be fabricated directly on the scintillator.

Comparison of the Biograph Vision and Biograph mCT for quantitative 90Y PET/CT imaging for radioembolisation

BACKGROUND

New digital PET scanners with improved time of flight timing and extended axial field of view such as the Siemens Biograph Vision have come on the market and are expected to replace current generation photomultiplier tube (PMT)-based systems such as the Siemens Biograph mCT. These replacements warrant a direct comparison between the systems, so that a smooth transition in clinical practice and research is guaranteed, especially when quantitative values are used for dosimetry-based treatment guidance. The new generation digital PET scanners offer increased sensitivity. This could particularly benefit 90Y imaging, which tends to be very noisy owing to the small positron branching ratio and high random fraction of 90Y. This study aims to determine the ideal reconstruction settings for the digital Vision for quantitative 90Y imaging and to evaluate the image quality and quantification of the digital Vision in comparison with its predecessor, the PMT-based mCT, for 90Y imaging in radioembolisation procedures.

METHODS

The NEMA image quality phantom was scanned to determine the ideal reconstruction settings for the Vision. In addition, an anthropomorphic phantom was scanned with both the Vision and the mCT, mimicking a radioembolisation patient with lung, liver, tumour, and extrahepatic deposition inserts. Image quantification of the anthropomorphic phantom was assessed by the lung shunt fraction, the tumour to non-tumour ratio, the parenchymal dose, and the contrast to noise ratio of extrahepatic depositions.

RESULTS

For the Vision, a reconstruction with 3 iterations, 5 subsets, and no post-reconstruction filter is recommended for quantitative 90Y imaging, based on the convergence of the recovery coefficient. Comparing both systems showed that the noise level of the Vision is significantly lower than that of the mCT (background variability of 14% for the Vision and 25% for the mCT at 2.5·103 MBq for the 37 mm sphere size). For quantitative 90Y measures, such as needed in radioembolisation, both systems perform similarly.

CONCLUSIONS

We recommend to reconstruct 90Y images acquired on the Vision with 3 iterations, 5 subsets, and no post-reconstruction filter for quantitative imaging. The Vision provides a reduced noise level, but similar quantitative accuracy as compared with its predecessor the mCT.

Simultaneous emission and attenuation reconstruction in time-of-flight PET using a reference object

Background

Simultaneous reconstruction of emission and attenuation images in time-of-flight (TOF) positron emission tomography (PET) does not provide a unique solution. In this study, we propose to solve this limitation by including additional information given by a reference object with known attenuation placed outside the patient. Different configurations of the reference object were studied including geometry, material composition, and activity, and an optimal configuration was defined. In addition, this configuration was tested for different timing resolutions and noise levels.

Results

The proposed strategy was tested in 2D simulations obtained by forward projection of available PET/CT data and noise was included using Monte Carlo techniques. Obtained results suggest that the optimal configuration corresponds to a water cylinder inserted in the patient table and filled with activity. In that case, mean differences between reconstructed and true images were below 10%. However, better results can be obtained by increasing the activity of the reference object.

Conclusion

This study shows promising results that might allow to obtain an accurate attenuation map from pure TOF-PET data without prior knowledge obtained from CT, MRI, or transmission scans.

Understanding the intrinsic radioactivity energy spectrum from 176Lu in LYSO/LSO scintillation crystals

Lutetium oxyorthosilicate (LSO) or lutetium yttrium oxyorthosilicate (LYSO) are the scintillator materials most widely used today in PET detectors due to their convenient physical properties for the detection of 511 keV annihilation photons. Natural lutetium contains 2.6% of ¹⁷⁶Lu which decays beta to excited states of ¹⁷⁶Hf producing a constant background signal. Although previous works have studied the background activity from LSO/LYSO, the shape of the spectrum, resulting from β-particle and γ radiation self-detection, has not been fully explained. The present work examines the contribution of the different β-particle and γ-ray interactions to provide a fuller comprehension of this background spectrum and to explain the differences observed when using crystals of different sizes. To this purpose we have shifted the continuous β-particle energy spectrum of ¹⁷⁶Lu from zero to the corresponding energy value for all combinations of the isomeric transitions of ¹⁷⁶Hf (γ-rays/internal conversion). The area of each shifted β-spectrum was normalized to reflect the probability of occurrence. To account for the probability of the γ-rays escaping from the crystal, Monte Carlo simulations using PENELOPE were performed in which point-like sources of monoenergetic photons were generated, inside LYSO square base prisms (all 1 cm thick) of different sizes: 1.0 cm to 5.74 cm. The analytic distributions were convolved using a varying Gaussian function to account for the measured energy resolution. The calculated spectra were compared to those obtained experimentally using monolithic crystals of the same dimensions coupled to SiPM arrays. Our results are in very good agreement with the experiment, and even explain the differences observed due to crystal size. This work may prove useful to calibrate and assess detector performance, and to measure energy resolution at different energy values.

The origin and reduction of spurious extrahepatic counts observed in 90Y non-TOF PET imaging post radioembolization

Our literature survey revealed a physical effect unknown to the nuclear medicine community, i.e. internal bremsstrahlung emission, and also the existence of long energy resolution tails in crystal scintillation. None of these effects has ever been modelled in PET Monte Carlo (MC) simulations. This study investigates whether these two effects could be at the origin of two unexplained observations in ⁹⁰Y imaging by PET: the increasing tails in the radial profile of true coincidences, and the presence of spurious extrahepatic counts post radioembolization in non-TOF PET and their absence in TOF PET. These spurious extrahepatic counts hamper the microsphere delivery check in liver radioembolization. An acquisition of a ³²P vial was performed on a GSO PET system. This is the ideal setup to study the impact of bremsstrahlung x-rays on the true coincidence rate when no positron emission and no crystal radioactivity are present. A MC simulation of the acquisition was performed using Gate-Geant4. MC simulations of non-TOF PET and TOF-PET imaging of a synthetic ⁹⁰Y human liver radioembolization phantom were also performed. Internal bremsstrahlung and long energy resolution tails inclusion in MC simulations quantitatively predict the increasing tails in the radial profile. In addition, internal bremsstrahlung explains the discrepancy previously observed in bremsstrahlung SPECT between the measure of the ⁹⁰Y bremsstrahlung spectrum and its simulation with Gate-Geant4. However the spurious extrahepatic counts in non-TOF PET mainly result from the failure of conventional random correction methods in such low count rate studies and poor robustness versus emission-transmission inconsistency. A novel proposed random correction method succeeds in cleaning the spurious extrahepatic counts in non-TOF PET. Two physical effects not considered up to now in nuclear medicine were identified to be at the origin of the unusual ⁹⁰Y true coincidences radial profile. TOF reconstruction removing of the spurious extrahepatic counts was theoretically explained by a better robustness against emission-transmission inconsistency. A novel random correction method was proposed to overcome the issue in non-TOF PET. Further studies are needed to assess the novel random correction method robustness.