Positron range effects of 66Ga in small-animal PET imaging.

PURPOSE: Gallium-66 is a non-conventional positron emitter that stands out not only for its high potential to label peptides, proteins and antibodies, but also because it can provide spatio-temporal information of relatively slow physiological processes in the body due to its conveniently long half-life of 9.5 h. However, 66Ga emits the most energetic positrons for PET imaging. The lack of information of the positron range effect on spatial resolution for this positron emitter is an issue, particularly in preclinical imaging. METHODS: The line spread function (LSF) in tissue-equivalent materials with densities between 0.2 and 1.93 g/cm3 was obtained with 66Ga and 18F. A complementary study with the NEMA NU 4-2008 image quality phantom is also included. RESULTS: High-energy positrons moving in lower density materials produce far-reaching activity distributions. The LSFs were characterized with Lorentzian-Gaussian fits, with spatial resolution (FWHM) in the 2.14-3.2 mm range, and long tails extending a few tens of mm depending on the material type and density. A narrowing of the LSF was observed for lung-equivalent materials, indicating the lack of enough material for the positron annihilation to take place. The NEMA NU 4-2008 image quality phantom produced blurred images, notoriously observed in the hot and cold cylinders used for evaluation of recovery coefficients (RC) and spill-over ratios (SOR), producing very low RC and very large SOR. CONCLUSIONS: Quantitative PET imaging with the non-conventional 66Ga is hampered due to the large range of its high-energy positrons affecting both spatial resolution and activity concentration quantification.

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 176Lu which decays beta to excited states of 176Hf 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 176Lu from zero to the corresponding energy value for all combinations of the isomeric transitions of 176Hf (γ-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.

Optimization of a large-area detector-block based on SiPM and pixelated LYSO crystal arrays.

We present the performance evaluation of a large-area detector module based on the ArrayC-60035-64P, an 8×8 array of tileable, 7.2mm pitch, silicon photomultipliers (SiPM) by SensL, covering a total area of 57.4mm×57.4mm. We characterized the ArrayC-60035-64P, operating at room temperature, using LYSO pixelated crystal arrays of different pitch sizes (1.075, 1.430, 1.683, 2.080 and 2.280mm) to determine the resolvable crystal size. After an optimization process, a 7mm thick coupling light guide was used for all crystal pitches. To identify the interaction position a 16-channel (8 columns, 8 rows) symmetric charge division (SCD) readout board together with a center-of-gravity algorithm was used. Based on this, we assembled the detector modules using a 40×40 LYSO, 1.43mm pitch array, covering the total detector area. Calibration was performed using a 137Cs source resulting in excellent crystal maps with minor geometric distortion, a mean 4.1 peak-to-valley ratio and 9.6% mean energy resolution for 662keV photons in the central region. The resolvability index was calculated in the x and y directions with values under 0.42 in all cases. We show that these large area SiPM arrays, combined with a 16-channel SCD readout board, can offer high spatial resolution, without processing a big number of signals, attaining excellent energy resolution and detector uniformity.

Positron range in tissue-equivalent materials: experimental microPET studies.

In this work an experimental investigation was carried out to study the effect that positron range has over positron emission tomography (PET) scans through measurements of the line spread function (LSF) in tissue-equivalent materials. Line-sources consisted of thin capillary tubes filled with (18)F, (13)N or (68)Ga water-solution inserted along the axis of symmetry of cylindrical phantoms constructed with the tissue-equivalent materials: lung (inhale and exhale), adipose tissue, solid water, trabecular and cortical bone. PET scans were performed with a commercial small-animal PET scanner and image reconstruction was carried out with filtered-backprojection. Line-source distributions were analyzed using radial profiles taken on axial slices from which the spatial resolution was determined through the full-width at half-maximum, tenth-maximum, twentieth-maximum and fiftieth-maximum. A double-Gaussian model of the LSFs was used to fit experimental data which can be incorporated into iterative reconstruction methods. In addition, the maximum activity concentration in the line-sources was determined from reconstructed images and compared to the known values for each case. The experimental data indicates that positron range in different materials has a strong effect on both spatial resolution and activity concentration quantification in PET scans. Consequently, extra care should be taken when computing standard-uptake values in PET scans, in particular when the radiopharmaceutical is taken up by different tissues in the body, and more even so with high-energy positron emitters.

Effect of correction methods of radiochromic EBT2 films on the accuracy of IMRT QA.

Radiochromic films are dosimeters with more favorable characteristics than other two-dimensional (2D) radiation detectors. Transmission film scanners using a linear charge-coupled device (CCD) array have a drawback of variation in response along the detector array, which may result in a nonuniform transmission of signal over the scanned image. This study uses 2D gamma index analysis to compare two methods for correcting the nonuniform response of EBT2: the first method was based on the new red:blue method proposed by the manufacturer (to compensate for small nonuniformities in the film coating) and the second method, proposed by Menegotti et al. (2008), was based on dose-dependent matrix correction factors. The gamma index analysis shows that both the methods are comparably accurate for all the criteria values used for evaluation (1 mm/1%, 2 mm/2%, 3 mm/3%). Centers around the world use both the methods to correct EBT2 local heterogeneities, but it is important to note that the former method has several advantages such as less time consumption and easy implementation.

Effect of dosimeter type for commissioning small photon beams on calculated dose distribution in stereotactic radiosurgery.

PURPOSE: To assess the impact of the detector used to commission small photon beams on the calculated dose distribution in stereotactic radiosurgery (SRS). METHODS: In this study, six types of detectors were used to characterize small photon beams: three diodes [a silicon stereotactic field diode SFD, a silicon diode SRS, and a silicon diode E], an ionization chamber CC01, and two types of radiochromic film models EBT and EBT2. These detectors were used to characterize circular collimated beams that were generated by a Novalis linear accelerator. This study was conducted in two parts. First, the following dosimetric data, which are of particular interest in SRS, were compared for the different detectors: the total scatter factor (TSF), the tissue phantom ratios (TPRs), and the off-axis ratios (OARs). Second, the commissioned data sets were incorporated into the treatment planning system (TPS) to compare the calculated dose distributions and the dose volume histograms (DVHs) that were obtained using the different detectors. RESULTS: The TSFs data measured by all of the detectors were in good agreement with each other within the respective statistical uncertainties: two exceptions, where the data were systematically below those obtained for the other detectors, were the CC01 results for all of the circular collimators and the EBT2 film results for circular collimators with diameters below 10.0 mm. The OAR results obtained for all of the detectors were in excellent agreement for all of the circular collimators. This observation was supported by the gamma-index test. The largest difference in the TPR data was found for the 4.0 mm circular collimator, followed by the 10.0 and 20.0 mm circular collimators. The results for the calculated dose distributions showed that all of the detectors passed the gamma-index test at 100% for the 3 mm/3% criteria. The aforementioned observation was true regardless of the size of the calculation grid for all of the circular collimators. Finally, the dose volume histogram results were independent of the size of the calculation grid used. CONCLUSIONS: The results of this study showed that all of the studied detectors produced similar commissioned data sets for the TPS dose calculations. However, this result only validated the dose distribution calculation in the TPS and could not be used to assess the dose delivery to the target in which the TFS data were used to calculate the monitor units (the TFS data were not used in the TPS dose distribution calculation). Therefore, this study could not be used to determine the most accurate detector commissioning data set; however, all of the detectors exhibited superior performance for the relative dosimetry of small photon beams.

Response of LiF:Mg,Ti to low energy carbon and oxygen ions.

The extended track interaction model (ETIM) has been formulated to explain the TL-fluence response for peak 5 to heavy ions using radial dose distributions produced by the ions in LiF and their luminescent centre occupation probability distributions. In this work, an experimental study of the TLD-100 fluence response to carbon and oxygen ions and its interpretation in terms of a Monte Carlo simulation of ETIM applied to peak 5 are presented. Irradiations were performed with 7.34 and 10.3 MeV 12C and 8.34 MeV 16O ions in the fluence interval between 2 x 10(7) and 2 x 10(11) cm(-2). Individual glow curve responses show the expected increase of supralinearity as the peak temperature increases. Data for peak 5 show a weak dependence with energy. These latter results are difficult to understand when one considers the differences in the expected radial occupancies for different ion energies.

Search for ionisation density effects in the radiation absorption stage in LiF:Mg,Ti.

Optical absorption (OA) dose-response of LiF:Mg,Ti (TLD-100) is studied as a function of electron energy (ionisation density) and irradiation dose. Contrary to the situation in thermoluminescence dose-response where the supralinearity is strongly energy-dependent, no dependence of the OA dose filling constants on energy is observed. This result is interpreted as indicating a lack of competitive process in the radiation absorption stage. The lack of an energy dependence of the dose filling constant also suggests that the charge carrier migration distances are sufficiently large to smear out the differences in the non-uniform distribution of ionisation events created by the impinging gamma/electron radiation of various energies.

Radiochromic dye film studies for brachytherapy applications.

Commercial radiochromic dye films have been used in recent years to quantify absorbed dose in several medical applications. In this study we present the characterisation of the GafChromic MD-55-2 dye film, a double sensitive layer film suitable for photon irradiation in brachytherapy applications. Dose measurements were carried out with a low dose rate 137Cs brachytherapy source, which produces very steep dose gradients in its vicinity, and therefore requires the capability of producing high spatial resolution isodose curves. Quantification of the dose rate in water per unit air kerma strength was obtained using a high-resolution transmission commercial scanner (Agfa DuoScan T1200 with the capability of digitising up to 600 x 1200 pixels per inch using 36 bits per pixel, together with optical density measurements. The Monte Carlo calculations and experimental measurements compared well in the 0-50 Gy dose interval used in this study.

Thermoluminescence induced by heavy charged particles.

Thermoluminescent materials display complex features when exposed to heavy charged particles: the TL response shows a linear-supralinear saturation dependence on the fluence and the efficiency is a strong function of the particle's velocity. Results are presented from measurements performed with low energy ions incident on TLD-100. The measured supralinearity can be reproduced by a Monte Carlo track interaction model which incorporates saturation effects; however, the model parameters which describe the ion's radii do not display the expected dependence on incident energy. The efficiency data support the choice of reference radiation having a secondary electron spectrum 'matched' to that of the heavy charged particles. The systematics of the measured glow curve shapes and maximum supralinearity unifies results from light ions and photons based on characteristics of the secondary electron spectra.

On the correct measurement of relative heavy charged particles to gamma thermoluminescent efficiencies.

In order to better understand the most important experimental aspects for performing correct measurements of relative thermoluminescent (TL) efficiencies, an investigation has been carried out to quantify the effect of using different experimental procedures in the evaluation of 3 MeV proton-to-gamma relative efficiency (etap,gamma) of LiF:Mg,Ti. Variations in batch, presentation, annealing and reader have been studied. When the same protocol is used to measure proton and gamma TL response, efficiency values obtained range from 0.36 to 0.59 for peak 5 and from 0.44 to 0.79 for the total signal. The use of different annealings and different batches leads to 20% and 10% differences in etap,gamma respectively. Large differences (40%) are found between efficiency values measured with TLD-100 chips and those obtained using TLD-100 microcubes. When 'mixed' procedures are used to measure the proton and the gamma response, differences in etap,gamma may increase even more. The main conclusion of this work is to stress the importance of measuring an entire series of experiments in the same laboratory with a carefully defined protocol and using dosemeters from the same batch to obtain heavy charged particle TL response and gamma TL response with identical annealing and readout procedures.

Supralinearity of peaks 5a, 5 and 5b in TLD-100 following 6.8 MeV and 2.6 MeV He ion irradiation: the extended track interaction model.

The extended track interaction model (ETIM) is applied to the TL fluence response of peaks 5, 5b and 5a following 6.8 MeV and 2.6 MeV He ion irradiation of LiF:Mg,Ti (TLD-100).

The TL fluence response to heavy charged particles using the track interaction model and track structure information.

The extended track interaction model, ETIM, has recently being proposed to explain the TLD-100 fluence response of peak 5 to heavy ions. This model includes the track structure information through the use of the luminescent centre occupation probability obtained from radial dose distributions produced by the ions as they travel through the dosemeter. In this work an implementation of ETIM using Monte Carlo techniques is presented. The simulation was applied to calculate the response of peak 5 of both sensitised and normal TLD-100 crystals to 2.6 and 6.8 MeV 4He ions. The simulation shows that the TL-fluence response has a strong dependence on ion energy, in disagreement with experimental observations. In spite of this, good agreement between the simulated TL-fluence response calculated for the 6.8 MeV 4He radial distributions and the experimental data for the two energies was achieved.

The use of a reflective scanner to study radiochromic film response.

GafChromic MD-55-2 film response was studied using a flatbed, reflective Microtek ScanMaker E3 scanner, under different conditions to optimize its use. Irradiations were performed using 60Co gamma rays in a 0-300 Gy dose interval to produce a calibrated step wedge. A 24-bit colour-scale mode along with image splitting into its red, green and blue components is suggested as an improved dosimetry method over a 256 grey level (8-bit) mode, extending the dynamic range for this film. Diverse film orientation and positioning, a black or a white background, and individual step film scanning were evaluated. Unwanted normalization is overcome by adding reference black and white steps adjacent to the radiochromic film, ensuring reproducibility. The use of a red filter was found to be equivalent to the use of the red component of the image after image colour splitting. The useful range for MD-55-2 film is extended up to 300 Gy if colour components, other than red, are used to evaluate the response. Comparisons with optical density measurements show that inexpensive commercial scanners might be a good alternative to densitometers.

Thermoluminescent response and relative efficiency of TLD-100 exposed to low-energy x-rays.

The dose-response of LiF:Mg,Ti (TLD-100) exposed to 15 and 35 kVp (8.0 +/- 0.1 and 8.1 +/- 0.1 keV effective energy respectively) x-rays and 60Co gamma-rays has been measured in the dose interval from (1.2-5.4) x 10(3) Gy for x-rays, and from 0.14 to 850 Gy for gamma-rays. In both cases the total TL signal and glow curve peaks 3 to 9 show supralinearity. The supralinearity function f(D) is similar for both x-ray beams, except for peak 8, where a 30% difference is observed. The maxima of f(D) for the total TL signal and peaks 5 to 8 are 2.1, 1.7, 6.4, 3.3 and 7.5 respectively for 8.1 keV x-rays and 3.7, 3.1, 13.6, 9.9 and 11.0 for gamma-rays. The measured relative efficiencies for x-rays with respect to 60Co, for the total TL signal and peaks 5 and 7, were 1.04, 0.97 and 3.2 respectively.

Monte Carlo simulation of depth-dose distributions in TLD-100 under 90Sr-90Y irradiation.

In this work the depth-dose distribution in TLD-100 dosimeters under beta irradiation from a 90Sr-90Y source was investigated using the Monte Carlo method. Comparisons between the simulated data and experimental results showed that the depth-dose distribution is strongly affected by the different components of both the source and dosimeter holders due to the large number of electron scattering events.

Comparative assessment of nine scatter correction methods based on spectral analysis using Monte Carlo simulations.

UNLABELLED: We compared nine scatter correction methods based on spectral analysis which process SPECT projections. METHODS: Monte Carlo simulation was used to generate histories of photons emitted from a realistic 99mTc phantom. A particular projection was considered. Information regarding the history, location and energy of the photons detected in this projection was analyzed to test the assumptions underlying each scatter correction method. Relative and absolute quantification and signal-to-noise ratio were assessed for each scatter corrected image. RESULTS: For the simulated data, two methods do not enable activity quantification. Among the methods requiring some parameters to be calibrated, the dual-energy window method shows the best compromise between accuracy and ease of implementation but introduces a bias in relative quantification. In this respect, a triple-energy window technique is more accurate than the dual-window method. A factor analysis approach results in more stable quantitative accuracy (error approximately 10%) for a wide range of activity but requires a more sophisticated acquisition mode (30 energy windows). CONCLUSION: These results show that a scatter correction method using spectral analysis can be used to substantially improve accurate quantification.