Comparative studies of the sensitivities of sparse and full geometries of Total-Body PET scanners built from crystals and plastic scintillators.

BACKGROUND: Alongside the benefits of Total-Body imaging modalities, such as higher sensitivity, single-bed position, low dose imaging, etc., their final construction cost prevents worldwide utilization. The main aim of this study is to present a simulation-based comparison of the sensitivities of existing and currently developed tomographs to introduce a cost-efficient solution for constructing a Total-Body PET scanner based on plastic scintillators. METHODS: For the case of this study, eight tomographs based on the uEXPLORER configuration with different scintillator materials (BGO, LYSO), axial field-of-view (97.4 cm and 194.8 cm), and detector configurations (full and sparse) were simulated. In addition, 8 J-PET scanners with different configurations, such as various axial field-of-view (200 cm and 250 cm), different cross sections of plastic scintillator, and multiple numbers of plastic scintillator layers (2, 3, and 4), based on J-PET technology have been simulated by GATE software. Furthermore, Siemens' Biograph Vision has been simulated to compare the results with standard PET scans. Two types of simulations have been performed. The first one with a centrally located source with a diameter of 1 mm and a length of 250 cm, and the second one with the same source inside a water-filled cylindrical phantom with a diameter of 20 cm and a length of 183 cm. RESULTS: With regards to sensitivity, among all the proposed scanners, the ones constructed with BGO crystals give the best performance ([Formula: see text] 350 cps/kBq at the center). The utilization of sparse geometry or LYSO crystals significantly lowers the achievable sensitivity of such systems. The J-PET design gives a similar sensitivity to the sparse LYSO crystal-based detectors while having full detector coverage over the body. Moreover, it provides uniform sensitivity over the body with additional gain on its sides and provides the possibility for high-quality brain imaging. CONCLUSION: Taking into account not only the sensitivity but also the price of Total-Body PET tomographs, which till now was one of the main obstacles in their widespread clinical availability, the J-PET tomography system based on plastic scintillators could be a cost-efficient alternative for Total-Body PET scanners.

Efficiency determination of J-PET: first plastic scintillators-based PET scanner.

BACKGROUND: The Jagiellonian Positron Emission Tomograph is the 3-layer prototype of the first scanner based on plastic scintillators, consisting of 192 half-metre-long strips with readouts at both ends. Compared to crystal-based detectors, plastic scintillators are several times cheaper and could be considered as a more economical alternative to crystal scintillators in future PETs. JPET is also a first multi-photon PET prototype. For the development of multi-photon detection, with photon characterized by the continuous energy spectrum, it is important to estimate the efficiency of J-PET as a function of energy deposition. The aim of this work is to determine the registration efficiency of the J-PET tomograph as a function of energy deposition by incident photons and the intrinsic efficiency of the J-PET scanner in detecting photons of different incident energies. In this study, 3-hit events are investigated, where 2-hits are caused by 511 keV photons emitted in [Formula: see text] annihilations, while the third hit is caused by one of the scattered photons. The scattered photon is used to accurately measure the scattering angle and thus the energy deposition. Two hits by a primary and a scattered photon are sufficient to calculate the scattering angle of a photon, while the third hit ensures the precise labeling of the 511 keV photons. RESULTS: By comparing experimental and simulated energy distribution spectra, the registration efficiency of the J-PET scanner was determined in the energy deposition range of 70-270 keV, where it varies between 20 and 100[Formula: see text]. In addition, the intrinsic efficiency of the J-PET was also determined as a function of the energy of the incident photons. CONCLUSION: A method for determining registration efficiency as a function of energy deposition and intrinsic efficiency as a function of incident photon energy of the J-PET scanner was demonstrated. This study is crucial for evaluating the performance of the scanner based on plastic scintillators and its applications as a standard and multi-photon PET systems. The method may be also used in the calibration of Compton-cameras developed for the ion-beam therapy monitoring and simultaneous multi-radionuclide imaging in nuclear medicine.

Testing CPT symmetry in ortho-positronium decays with positronium annihilation tomography.

Charged lepton system symmetry under combined charge, parity, and time-reversal transformation (CPT) remains scarcely tested. Despite stringent quantum-electrodynamic limits, discrepancies in predictions for the electron-positron bound state (positronium atom) motivate further investigation, including fundamental symmetry tests. While CPT noninvariance effects could be manifested in non-vanishing angular correlations between final-state photons and spin of annihilating positronium, measurements were previously limited by knowledge of the latter. Here, we demonstrate tomographic reconstruction techniques applied to three-photon annihilations of ortho-positronium atoms to estimate their spin polarisation without magnetic field or polarised positronium source. We use a plastic-scintillator-based positron-emission-tomography scanner to record ortho-positronium (o-Ps) annihilations with single-event estimation of o-Ps spin and determine the complete spectrum of an angular correlation operator sensitive to CPT-violating effects. We find no violation at the precision level of 10-4, with an over threefold improvement on the previous measurement.

Optimisation of the event-based TOF filtered back-projection for online imaging in total-body J-PET.

We perform a parametric study of the newly developed time-of-flight (TOF) image reconstruction algorithm, proposed for the real-time imaging in total-body Jagiellonian PET (J-PET) scanners. The asymmetric 3D filtering kernel is applied at each most likely position of electron-positron annihilation, estimated from the emissions of back-to-back γ-photons. The optimisation of its parameters is studied using Monte Carlo simulations of a 1-mm spherical source, NEMA IEC and XCAT phantoms inside the ideal J-PET scanner. The combination of high-pass filters which included the TOF filtered back-projection (FBP), resulted in spatial resolution, 1.5 times higher in the axial direction than for the conventional 3D FBP. For realistic 10-minute scans of NEMA IEC and XCAT, which require a trade-off between the noise and spatial resolution, the need for Gaussian TOF kernel components, coupled with median post-filtering, is demonstrated. The best sets of 3D filter parameters were obtained by the Nelder-Mead minimisation of the mean squared error between the resulting and reference images. The approach allows training the reconstruction algorithm for custom scans, using the IEC phantom, when the temporal resolution is below 50 ps. The image quality parameters, estimated for the best outcomes, were systematically better than for the non-TOF FBP.

Simulating NEMA characteristics of the modular total-body J-PET scanner-an economic total-body PET from plastic scintillators.

The purpose of the presented research is estimation of the performance characteristics of the economic total-body Jagiellonian-PET system (TB-J-PET) constructed from plastic scintillators. The characteristics are estimated according to the NEMA NU-2-2018 standards utilizing the GATE package. The simulated detector consists of 24 modules, each built out of 32 plastic scintillator strips (each with cross section of 6 mm times 30 mm and length of 140 or 200 cm) arranged in two layers in regular 24-sided polygon circumscribing a circle with the diameter of 78.6 cm. For the TB-J-PET with an axial field-of-view (AFOV) of 200 cm, a spatial resolutions (SRs) of 3.7 mm (transversal) and 4.9 mm (axial) are achieved. The noise equivalent count rate (NECR) peak of 630 kcps is expected at 30 kBq cc-1. Activity concentration and the sensitivity at the center amounts to 38 cps kBq-1. The scatter fraction (SF) is estimated to 36.2 %. The values of SF and SR are comparable to those obtained for the state-of-the-art clinical PET scanners and the first total-body tomographs: uExplorer and PennPET. With respect to the standard PET systems with AFOV in the range from 16 to 26 cm, the TB-J-PET is characterized by an increase in NECR approximately by factor of 4 and by the increase of the whole-body sensitivity by factor of 12.6 to 38. The time-of-flight resolution for the TB-J-PET is expected to be at the level of CRT = 240 ps full width at half maximum. For the TB-J-PET with an AFOV of 140 cm, an image quality of the reconstructed images of a NEMA IEC phantom was presented with a contrast recovery coefficient and a background variability parameters. The increase of the whole-body sensitivity and NECR estimated for the TB-J-PET with respect to current commercial PET systems makes the TB-J-PET a promising cost-effective solution for the broad clinical applications of total-body PET scanners. TB-J-PET may constitute an economic alternative for the crystal TB-PET scanners, since plastic scintillators are much cheaper than BGO or LYSO crystals and axial arrangement of the strips significantly reduces the costs of readout electronics and SiPMs.

3D TOF-PET image reconstruction using total variation regularization.

In this paper we introduce a semi-analytic algorithm for 3-dimensional image reconstruction for positron emission tomography (PET). The method consists of the back-projection of the acquired data into the most likely image voxel according to time-of-flight (TOF) information, followed by the filtering step in the image space using an iterative optimization algorithm with a total variation (TV) regularization. TV regularization in image space is more computationally efficient than usual iterative optimization methods for PET reconstruction with full system matrix that use TV regularization. The efficiency comes from the one-time TOF back-projection step that might also be described as a reformatting of the acquired data. An important aspect of our work concerns the evaluation of the filter operator of the linear transform mapping an original radioactive tracer distribution into the TOF back-projected image. We obtain concise, closed-form analytical formula for the filter operator. The proposed method is validated with the Monte Carlo simulations of the NEMA IEC phantom using a one-layer, 50 cm-long cylindrical device called Jagiellonian PET scanner. The results show a better image quality compared with the reference TOF maximum likelihood expectation maximization algorithm.

Performance assessment of the 2 γpositronium imaging with the total-body PET scanners.

PURPOSE: In living organisms, the positron-electron annihilation (occurring during the PET imaging) proceeds in about 30% via creation of a metastable ortho-positronium atom. In the tissue, due to the pick-off and conversion processes, over 98% of ortho-positronia annihilate into two 511 keV photons. In this article, we assess the feasibility for reconstruction of the mean ortho-positronium lifetime image based on annihilations into two photons. The main objectives of this work include the (i) estimation of the sensitivity of the total-body PET scanners for the ortho-positronium mean lifetime imaging using 2γ annihilations and (ii) estimation of the spatial and time resolution of the ortho-positronium image as a function of the coincidence resolving time (CRT) of the scanner. METHODS: Simulations are conducted assuming that radiopharmaceutical is labeled with 44Sc isotope emitting one positron and one prompt gamma. The image is reconstructed on the basis of triple coincidence events. The ortho-positronium lifetime spectrum is determined for each voxel of the image. Calculations were performed for cases of total-body detectors build of (i) LYSO scintillators as used in the EXPLORER PET and (ii) plastic scintillators as anticipated for the cost-effective total-body J-PET scanner. To assess the spatial and time resolution, the four cases were considered assuming that CRT is equal to 500 ps, 140 ps, 50 ps, and 10 ps. RESULTS: The estimated total-body PET sensitivity for the registration and selection of image forming triple coincidences (2γ+γprompt) is larger by a factor of 13.5 (for LYSO PET) and by factor of 5.2 (for plastic PET) with respect to the sensitivity for the standard 2γ imaging by LYSO PET scanners with AFOV = 20 cm. The spatial resolution of the ortho-positronium image is comparable with the resolution achievable when using TOF-FBP algorithms already for CRT = 50 ps. For the 20-min scan, the resolution better than 20 ps is expected for the mean ortho-positronium lifetime image determination. CONCLUSIONS: Ortho-positronium mean lifetime imaging based on the annihilations into two photons and prompt gamma is shown to be feasible with the advent of the high sensitivity total-body PET systems and time resolution of the order of tens of picoseconds.

Estimating relationship between the time over threshold and energy loss by photons in plastic scintillators used in the J-PET scanner.

PURPOSE: The time-over-threshold (TOT) technique is being used widely due to itsimplications in developing the multi-channel readouts, mainly when fast signal processing is required. Using the TOT technique, as a measure of energy loss instead of charge integration methods, significantly reduces the signal readout costs by combining the time and energy information. Therefore, this approach can potentially be utilized in J-PET tomograph which is built from plastic scintillators characterized by fast light signals. The drawback in adopting this technique lies in the non-linear correlation between input energy loss and TOT of the signal. The main motivation behind this work is to develop the relationship between TOT and energy loss and validate it by the J-PET tomograph setup. METHODS: The experiment was performed using a 22Na beta emitter source placed in the center of the J-PET tomograph. This isotope produces photons of two different energies: 511 keV photons from the positron annihilation (direct annihilation or through the formation of a para-positronium atom or pick-off process of ortho-positronium atoms) and a 1275 keV prompt photon. This allows the study of the correlation between TOT values and energy loss for energy ranges up to 1000 keV. Since the photon interacts predominantly via Compton scattering inside the plastic scintillator, there is no direct information of the energy deposition. However, using the J-PET geometry, one can measure the scattering angle of the interacting photon. Since the 22Na source emits photons of two different energies, it is necessary to know unambiguously the energy of incident photons and their corresponding scattering angles in order to estimate energy deposition. In summary, this work presents a dedicated algorithm developed to tag photons of different energies and studying their scattering angles to calculate the energy deposition by the interacting photons. RESULTS: A new method was elaborated to measure the energy loss by photons interacting with plastic scintillators used in the J-PET tomograph. We find the relationship between the energy loss and TOT is non-linear and can be described by the functions TOT = A0 + A1 * ln(E dep + A2) + A3 * (ln(E dep + A2))2 and TOT = A0 - A1 * A2[Formula: see text]. In addition, we also introduced a theoretical model to calculate the TOT as a function of energy loss in plastic scintillators. CONCLUSIONS: A relationship between TOT and energy loss by photons interacting inside the plastic scintillators used in J-PET scanner is established for a deposited energy range of 100-1000 keV.

Directed, Elliptic, and Higher Order Flow Harmonics of Protons, Deuterons, and Tritons in Au+Au Collisions at sqrt[s_{NN}]=2.4 GeV.

Flow coefficients v_{n} of the orders n=1-6 are measured with the High-Acceptance DiElectron Spectrometer (HADES) at GSI for protons, deuterons, and tritons as a function of centrality, transverse momentum, and rapidity in Au+Au collisions at sqrt[s_{NN}]=2.4 GeV. Combining the information from the flow coefficients of all orders allows us to construct for the first time, at collision energies of a few GeV, a multidifferential picture of the angular emission pattern of these particles. It reflects the complicated interplay between the effect of the central fireball pressure on the emission of particles and their subsequent interaction with spectator matter. The high precision information on higher order flow coefficients is a major step forward in constraining the equation of state of dense baryonic matter.

Strong Absorption of Hadrons with Hidden and Open Strangeness in Nuclear Matter.

We present the first observation of K^{-} and ϕ absorption within nuclear matter by means of π^{-}-induced reactions on C and W targets at an incident beam momentum of 1.7 GeV/c studied with HADES at SIS18/GSI. The double ratio (K^{-}/K^{+})_{W}/(K^{-}/K^{+})_{C} is found to be 0.319±0.009(stat)_{-0.012}^{+0.014}(syst) indicating a larger absorption of K^{-} in heavier targets as compared to lighter ones. The measured ϕ/K^{-} ratios in π^{-}+C and π^{-}+W reactions within the HADES acceptance are found to be equal to 0.55±0.04(stat)_{-0.07}^{+0.06}(syst) and to 0.63±0.06(stat)_{-0.11}^{+0.11}(syst), respectively. The similar ratios measured in the two different reactions demonstrate for the first time experimentally that the dynamics of the ϕ meson in nuclear medium is strongly coupled to the K^{-} dynamics. The large difference in the ϕ production off C and W nuclei is discussed in terms of a strong ϕN in-medium coupling. These results are relevant for the description of heavy-ion collisions and the structure of neutron stars.

Feasibility study of the positronium imaging with the J-PET tomograph.

A detection system of the conventional PET tomograph is set-up to record data from [Formula: see text] annihilation into two photons with energy of 511 keV, and it gives information on the density distribution of a radiopharmaceutical in the body of the object. In this paper we explore the possibility of performing the three gamma photons imaging based on ortho-positronium annihilation, as well as the possibility of positronium mean lifetime imaging with the J-PET tomograph constructed from plastic scintillators. For this purposes simulations of the ortho-positronium formation and its annihilation into three photons were performed taking into account distributions of photons' momenta as predicted by the theory of quantum electrodynamics and the response of the J-PET tomograph. In order to test the proposed ortho-positronium lifetime image reconstruction method, we concentrate on the decay of the ortho-positronium into three photons and applications of radiopharmaceuticals labeled with isotopes emitting a prompt gamma. The proposed method of imaging is based on the determination of hit-times and hit-positions of registered photons which enables the reconstruction of the time and position of the annihilation point as well as the lifetime of the ortho-positronium on an event-by-event basis. We have simulated the production of the positronium in point-like sources and in a cylindrical phantom composed of a set of different materials in which the ortho-positronium lifetime varied from 2.0 ns to 3.0 ns, as expected for ortho-positronium created in the human body. The presented reconstruction method for total-body J-PET like detector allows to achieve a mean lifetime resolution of ∼40 ps. Recent positron annihilation lifetime spectroscopy measurements of cancerous and healthy uterine tissues show that this sensitivity may allow to study the morphological changes in cell structures.

Estimating the NEMA characteristics of the J-PET tomograph using the GATE package.

A novel whole-body positron emission tomography (PET) system based on plastic scintillators is developed by the J-PET Collaboration. It consists of plastic scintillator strips arranged axially in the form of a cylinder, allowing the cost-effective construction of the total-body PET system. In order to determine the properties of the scanner prototype and optimize its geometry, advanced computer simulations were performed using the GATE (Geant4 application for tomographic emission) software. The spatial resolution, sensitivity, scatter fraction and noise equivalent count rate were estimated according to the National Electrical Manufacturers Association norm, as a function of the length of the tomograph, the number of detection layers, the diameter of the tomographic chamber and for various types of applied readout. For the single-layer geometry with a diameter of 85 cm, a strip length of 100 cm, a cross-section of 4 mm × 20 mm and silicon photomultipliers with an additional layer of wavelength shifter as the readout, the spatial resolution (full width at half maximum) in the centre of the scanner is equal to 3 mm (radial, tangential) and 6 mm (axial). For the analogous double-layer geometry with the same readout, diameter and scintillator length, with a strip cross-section of 7 mm × 20 mm, a noise equivalent count rate peak of 300 kcps was reached at 40 kBq cc-1 activity concentration, the scatter fraction is estimated to be about 35% and the sensitivity at the centre amounts to 14.9 cps kBq-1. Sensitivity profiles were also determined.

Evaluation of Single-Chip, Real-Time Tomographic Data Processing on FPGA SoC Devices.

A novel approach to tomographic data processing has been developed and evaluated using the Jagiellonian positron emission tomography scanner as an example. We propose a system in which there is no need for powerful, local to the scanner processing facility, capable to reconstruct images on the fly. Instead, we introduce a field programmable gate array system-on-chip platform connected directly to data streams coming from the scanner, which can perform event building, filtering, coincidence search, and region-of-response reconstruction by the programmable logic and visualization by the integrated processors. The platform significantly reduces data volume converting raw data to a list-mode representation, while generating visualization on the fly.

Feasibility studies of the polarization of photons beyond the optical wavelength regime with the J-PET detector.

J-PET is a detector optimized for registration of photons from the electron-positron annihilation via plastic scintillators where photons interact predominantly via Compton scattering. Registration of both primary and scattered photons enables to determinate the linear polarization of the primary photon on the event by event basis with a certain probability. Here we present quantitative results on the feasibility of such polarization measurements of photons from the decay of positronium with the J-PET and explore the physical limitations for the resolution of the polarization determination of 511 keV photons via Compton scattering. For scattering angles of about 82 ∘ (where the best contrast for polarization measurement is theoretically predicted) we find that the single event resolution for the determination of the polarization is about 40 ∘ (predominantly due to properties of the Compton effect). However, for samples larger than ten thousand events the J-PET is capable of determining relative average polarization of these photons with the precision of about few degrees. The obtained results open new perspectives for studies of various physics phenomena such as quantum entanglement and tests of discrete symmetries in decays of positronium and extend the energy range of polarization measurements by five orders of magnitude beyond the optical wavelength regime.

Calculation of the time resolution of the J-PET tomograph using kernel density estimation.

In this paper we estimate the time resolution of the J-PET scanner built from plastic scintillators. We incorporate the method of signal processing using the Tikhonov regularization framework and the kernel density estimation method. We obtain simple, closed-form analytical formulae for time resolution. The proposed method is validated using signals registered by means of the single detection unit of the J-PET tomograph built from a 30 cm long plastic scintillator strip. It is shown that the experimental and theoretical results obtained for the J-PET scanner equipped with vacuum tube photomultipliers are consistent.

A feasibility study of ortho-positronium decays measurement with the J-PET scanner based on plastic scintillators.

We present a study of the application of the Jagiellonian positron emission tomograph (J-PET) for the registration of gamma quanta from decays of ortho-positronium (o-Ps). The J-PET is the first positron emission tomography scanner based on organic scintillators in contrast to all current PET scanners based on inorganic crystals. Monte Carlo simulations show that the J-PET as an axially symmetric and high acceptance scanner can be used as a multi-purpose detector well suited to pursue research including e.g. tests of discrete symmetries in decays of ortho-positronium in addition to the medical imaging. The gamma quanta originating from o-Ps decay interact in the plastic scintillators predominantly via the Compton effect, making the direct measurement of their energy impossible. Nevertheless, it is shown in this paper that the J-PET scanner will enable studies of the [Formula: see text] decays with angular and energy resolution equal to [Formula: see text] and [Formula: see text], respectively. An order of magnitude shorter decay time of signals from plastic scintillators with respect to the inorganic crystals results not only in better timing properties crucial for the reduction of physical and instrumental background, but also suppresses significantly the pile-ups, thus enabling compensation of the lower efficiency of the plastic scintillators by performing measurements with higher positron source activities.

Time resolution of the plastic scintillator strips with matrix photomultiplier readout for J-PET tomograph.

Recent tests of a single module of the Jagiellonian Positron Emission Tomography system (J-PET) consisting of 30 cm long plastic scintillator strips have proven its applicability for the detection of annihilation quanta (0.511 MeV) with a coincidence resolving time (CRT) of 0.266 ns. The achieved resolution is almost by a factor of two better with respect to the current TOF-PET detectors and it can still be improved since, as it is shown in this article, the intrinsic limit of time resolution for the determination of time of the interaction of 0.511 MeV gamma quanta in plastic scintillators is much lower. As the major point of the article, a method allowing to record timestamps of several photons, at two ends of the scintillator strip, by means of matrix of silicon photomultipliers (SiPM) is introduced. As a result of simulations, conducted with the number of SiPM varying from 4 to 42, it is shown that the improvement of timing resolution saturates with the growing number of photomultipliers, and that the [Formula: see text] configuration at two ends allowing to read twenty timestamps, constitutes an optimal solution. The conducted simulations accounted for the emission time distribution, photon transport and absorption inside the scintillator, as well as quantum efficiency and transit time spread of photosensors, and were checked based on the experimental results. Application of the [Formula: see text] matrix of SiPM allows for achieving the coincidence resolving time in positron emission tomography of [Formula: see text]0.170 ns for 15 cm axial field-of-view (AFOV) and [Formula: see text]0.365 ns for 100 cm AFOV. The results open perspectives for construction of a cost-effective TOF-PET scanner with significantly better TOF resolution and larger AFOV with respect to the current TOF-PET modalities.