Validation of the GATE Monte Carlo simulation platform for modelling a CsI(Tl) scintillation camera dedicated to small-animal imaging

Monte Carlo simulations are increasingly used in scintigraphic imaging to model imaging systems and to develop and assess tomographic reconstruction algorithms and correction methods for improved image quantitation. GATE (GEANT4 application for tomographic emission) is a new Monte Carlo simulation platform based on GEANT4 dedicated to nuclear imaging applications. This paper describes the GATE simulation of a prototype of scintillation camera dedicated to small-animal imaging and consisting of a CsI(Tl) crystal array coupled to a position-sensitive photomultiplier tube. The relevance of GATE to model the camera prototype was assessed by comparing simulated 99mTc point spread functions, energy spectra, sensitivities, scatter fractions and image of a capillary phantom with the corresponding experimental measurements. Results showed an excellent agreement between simulated and experimental data: experimental spatial resolutions were predicted with an error less than 100 microns. The difference between experimental and simulated system sensitivities for different source-to-collimator distances was within 2%. Simulated and experimental scatter fractions in a [98-182 keV] energy window differed by less than 2% for sources located in water. Simulated and experimental energy spectra agreed very well between 40 and 180 keV. These results demonstrate the ability and flexibility of GATE for simulating original detector designs. The main weakness of GATE concerns the long computation time it requires: this issue is currently under investigation by the GEANT4 and the GATE collaborations.     

Dosimetric evaluation of a dedicated stereotactic linear accelerator using measurement and Monte Carlo simulation

Dosimetric parameters of a dedicated stereotactic linear accelerator have been investigated using measurements and Monte Carlo simulations. This linac has a unique built in multileaf collimation (MLC) system with the maximum opening of 16 x 21 cm2 and 4 mm leaf width at the isocenter and has successfully been modeled for the first time using the Monte Carlo simulation. The high resolution MLC, combined with its relatively large maximum field size, opens up a new opportunity for expanding stereotactic radiation treatment techniques from traditionally treating smaller targets to larger ones for both cranial and extracranial lesions. Dosimetric parameters of this linac such as accuracy of leaf positioning and field shaping, leakage and transmission, percentage depth doses, off-axes dose profiles, and dose penumbras were measured and calculated for different field sizes, depths, and source to surface distances. In addition, the ability of the linac in accurate dose delivery of several treatment plans, including intensity modulated radiation therapy (IMRT), performed on phantom and patients was determined. Ionization chamber, photon diode detector, films, several solid water phantoms, and a water tank were used for the measurements. The MLC leaf positioning to any particular point in the maximum aperture was accurate with a standard deviation of 0.29 mm. Maximum and average leakages were 1.7% and 1.1% for the reference field of 10.4 x 9.6 cm2. Measured penumbra widths (80%-20%) for this field at source axis distance (SAD) of 100 cm at a depth of 1.5 cm (dmax) were 3.2 and 4 mm for the leaf-sides and leaf-ends, respectively. The corresponding results at 10 cm depth and SAD =100 cm were 5.4 and 6.3 mm. Monte Carlo results generally agreed with the measurements to within 1% and or 1 mm, with respective uncertainties of 0.5% and 0.2 mm. The linac accuracy in delivering non-IMRT treatment plans was better than 1%. Ionization chamber dosimetry results for a phantom IMRT plan in the high dose and low dose regions were -0.5% and +3.6%, respectively. Dosimetry results at isocenter for three patients' IMRT plans were measured to be within 3% of their corresponding treatment plans. Film dosimetry was also used to compare dose distributions of IMRT treatment plans and delivered cumulative doses at different cross sectional planes.     

Characterization of dual layer phoswich detector performance for small animal PET using Monte Carlo simulation

A positron emission tomograph dedicated to small animal imaging should have high spatial resolution and sensitivity, and dual layer scintillators have been developed for this purpose. In this study, simulations were performed to optimize the order and the length of each crystal of a dual layer phoswich detector, and to evaluate the possibility of measuring signals from each layer of the phoswich detector. A simulation tool GATE was used to estimate the sensitivity and resolution of a small PET scanner. The proposed scanner is based on dual layer phoswich detector modules arranged in a ring of 10 cm diameter. Each module is composed of 8 x 8 arrays of phoswich detectors consisting of LSO and LuYAP with a 2 mm x 2 mm sensitive area coupled to a Hamamatsu R7600-00-M64 PSPMT. The length of the front layer of the phoswich detector varied from 0 to 10 mm at 1 mm intervals, and the total length (LSO + LuYAP) was fixed at 20 mm. The order of the crystal layers of the phoswich detector was also changed. Radial resolutions were kept below 3.4 mm and 3.7 mm over 8 cm FOV, and sensitivities were 7.4% and 8.0% for LSO 5 mm-LuYAP 15 mm, and LuYAP 6 mm-LSO 14 mm phoswich detectors, respectively. Whereas, high and uniform resolutions were achieved by using the LSO front layer, higher sensitivities were obtained by changing the crystal order. The feasibilities for applying crystal identification methods to phoswich detectors consisting of LSO and LuYAP were investigated using simulation and experimentally derived measurements of the light outputs from each layer of the phoswich detector. In this study, the optimal order and lengths of the dual layer phoswich detector were derived in order to achieve high sensitivity and high and uniform radial resolution.     

Feasibility of a novel design of high resolution parallax-free Compton enhanced PET scanner dedicated to brain research

A novel concept for a positron emission tomography (PET) camera module is proposed, which provides full 3D reconstruction with high resolution over the total detector volume, free of parallax errors. The key components are a matrix of long scintillator crystals and hybrid photon detectors (HPDs) with matched segmentation and integrated readout electronics. The HPDs read out the two ends of the scintillator package. Both excellent spatial (x, y, z) and energy resolution are obtained. The concept allows enhancing the detection efficiency by reconstructing a significant fraction of events which underwent Compton scattering in the crystals. The proof of concept will first be demonstrated with yttrium orthoaluminate perovskite (YAP):Ce crystals, but the final design will rely on other scintillators more adequate for PET applications (e.g. LSO:Ce or LaBr3:Ce). A promising application of the proposed camera module, which is currently under development, is a high resolution 3D brain PET camera with an axial field-of-view of approximately 15 cm dedicated to brain research. The design philosophy and performance predictions based on analytical calculations and Monte Carlo simulations are presented. Image correction and reconstruction tools required to operate this transmissionless device in a research environment are also discussed. Better or similar performance parameters were obtained compared to other known designs at lower fabrication cost. The axial geometrical concept also seems to be promising for applications such as positron emission mammography.     

Monte Carlo simulation of the spatial resolution and depth sensitivity of two-dimensional optical imaging of the brain

Absorption or fluorescence-based two-dimensional (2-D) optical imaging is widely employed in functional brain imaging. The image is a weighted sum of the real signal from the tissue at different depths. This weighting function is defined as "depth sensitivity." Characterizing depth sensitivity and spatial resolution is important to better interpret the functional imaging data. However, due to light scattering and absorption in biological tissues, our knowledge of these is incomplete. We use Monte Carlo simulations to carry out a systematic study of spatial resolution and depth sensitivity for 2-D optical imaging methods with configurations typically encountered in functional brain imaging. We found the following: (i) the spatial resolution is <200 μm for NA≤0.2 or focal plane depth≤300 μm. (ii) More than 97% of the signal comes from the top 500 μm of the tissue. (iii) For activated columns with lateral size larger than spatial resolution, changing numerical aperature (NA) and focal plane depth does not affect depth sensitivity. (iv) For either smaller columns or large columns covered by surface vessels, increasing NA and/or focal plane depth may improve depth sensitivity at deeper layers. Our results provide valuable guidance for the optimization of optical imaging systems and data interpretation.     

Monte Carlo simulation of electron beams generated by a 12 MeV dedicated mobile IORT accelerator

The aim of this study was to investigate the dosimetric characteristics of the electron beams generated by the light intraoperative accelerator, Liac? (SORDINA, Italy), using Monte Carlo (MC) calculations. Moreover we investigated the possibility of characterizing the Liac? dosimetry with a minimal set of dosimetric data. In fact accelerator commissioning requires measurements of both percentage depth doses (PDDs) and off-axis profiles for all the possible combinations of energy, applicator diameter and bevelled angle. The Liac? geometry and water phantom were simulated in a typical measurement setup, using the MC code EGSnrc/BEAMnrc. A simulated annealing optimization algorithm was used in order to find the optimal non-monoenergetic spectrum of the initial electron beam that minimizes the differences between calculated and measured PDDs. We have concluded that, for each investigated nominal energy beam, only the PDDs of applicators with diameters of 30, 70 and 100 mm and the PDD without an applicator were needed to find the optimal spectra. Finally, the output factors of the entire set of applicator diameters/bevelled angles were calculated. The differences between calculated and experimental output factors were better than 2%, with the exception of the smallest applicator which gave differences between 3% and 4% for all energies. The code turned out to be useful for checking the experimental data from various Liac? beams and will be the basis for developing a tool based on MC simulation to support the medical physicist in the commissioning phase.     

Study of the point spread function (PSF) for 123I SPECT imaging using Monte Carlo simulation

The iterative reconstruction algorithms employed in brain single-photon emission computed tomography (SPECT) allow some quantitative parameters of the image to be improved. These algorithms require accurate modelling of the so-called point spread function (PSF). Nowadays, most in vivo neurotransmitter SPECT studies employ pharmaceuticals radiolabelled with 123I. In addition to an intense line at 159 keV, the decay scheme of this radioisotope includes some higher energy gammas which may have a non-negligible contribution to the PSF. The aim of this work is to study this contribution for two low-energy high-resolution collimator configurations, namely, the parallel and the fan beam. The transport of radiation through the material system is simulated with the Monte Carlo code PENELOPE. We have developed a main program that deals with the intricacies associated with tracking photon trajectories through the geometry of the collimator and detection systems. The simulated PSFs are partly validated with a set of experimental measurements that use the 511 keV annihilation photons emitted by a 18F source. Sensitivity and spatial resolution have been studied, showing that a significant fraction of the detection events in the energy window centred at 159 keV (up to approximately 49% for the parallel collimator) are originated by higher energy gamma rays, which contribute to the spatial profile of the PSF mostly outside the 'geometrical' region dominated by the low-energy photons. Therefore, these high-energy counts are to be considered as noise, a fact that should be taken into account when modelling PSFs for reconstruction algorithms. We also show that the fan beam collimator gives higher signal-to-noise ratios than the parallel collimator for all the source positions analysed.     

PETSIM: Monte Carlo simulation of all sensitivity and resolution parameters of cylindrical positron imaging systems.

Monte Carlo simulation techniques are applied to track the annihilation photons from positron decay, and store the photon histories. Reasonably realistic models of the isotope distribution in the brain and heart during typical PET studies, as well as the traditional phantoms used for measuring PET scanner performance can be built out of up to 10 hollow or solid cylinders. Separate programs model the source distribution and its attenuation characteristics, the collimators and the detectors. These modules are connected by compact gamma history files which are stored on disc or tape. Over 50 million gamma ray histories can be saved on a 1 Gbyte disc, representing the decay of several billion atoms. This allows for good precision even for single thin slices in scanners with wide axial acceptance. The simulation results include spectrum analysis, sensitivity to true coincident events, scattered coincident and single rays, and the effects on these parameters of detector dead time. The storage of intermediate results on tape reduces simulation time, since most common source geometries need be generated only once. The sensitivities in multi-slice systems are presented as matrices of coincident crystal planes. The matrix shows the true count sensitivity and the scatter fraction together for each valid combination of planes. This presentation is very useful for assessing the effects of various degrees of inter-plane collimation. The spatial resolution analysis includes the effects of positron range, non-collinearity of the gamma rays, multiple interaction within the detectors, and the effects of quantization into single crystals in multiple-crystal block detectors. Each of these effects can be turned on or off without repeating the simulation. Both in-plane and axial resolutions are calculated as a function of location of the positron-emitting nucleus and the angle of incidence of gamma rays on the crystals. Single crystals, blocks and crystals with depth of interaction encoding can be specified, as can the method of backprojection (planar, or 3D), so that the detector geometry can be optimized.     

Monte Carlo simulation of the image formation process in portal imaging

We have written Monte Carlo programs to simulate the formation of radiological images. Our code is used to propagate a simulated x-ray fluence through each component of an existing video-based portal imaging system. This simulated fluence consists of a 512 x 512 pixel image containing both contrast-detail patterns as well as checker patterns to assess spatial resolution of the simulated portal imager. All of the components of the portal imaging system were modeled as a cascade of eight linear stages. Using this code, one can assess the visual impact of changing components in the imaging chain by changing the appropriate probability density function. Virtual experiments were performed to assess the visual impact of replacing the lens and TV camera by an amorphous silicon array, and the effect of scattered radiation on portal images.     

Monte Carlo simulations of a scintillation camera using GATE: validation and application modelling

Geant4 application for tomographic emission (GATE) is a recently developed simulation platform based on Geant4, specifically designed for PET and SPECT studies. In this paper we present validation results of GATE based on the comparison of simulations against experimental data, acquired with a standard SPECT camera. The most important components of the scintillation camera were modelled. The photoelectric effect. Compton and Rayleigh scatter are included in the gamma transport process. Special attention was paid to the processes involved in the collimator: scatter, penetration and lead fluorescence. A LEHR and a MEGP collimator were modelled as closely as possible to their shape and dimensions. In the validation study, we compared the simulated and measured energy spectra of different isotopes: 99mTc, 22Na, 57Co and 67Ga. The sensitivity was evaluated by using sources at varying distances from the detector surface. Scatter component analysis was performed in different energy windows at different distances from the detector and for different attenuation geometries. Spatial resolution was evaluated using a 99mTc source at various distances. Overall results showed very good agreement between the acquisitions and the simulations. The clinical usefulness of GATE depends on its ability to use voxelized datasets. Therefore, a clinical extension was written so that digital patient data can be read in by the simulator as a source distribution or as an attenuating geometry. Following this validation we modelled two additional camera designs: the Beacon transmission device for attenuation correction and the Solstice scanner prototype with a rotating collimator. For the first setup a scatter analysis was performed and for the latter design. the simulated sensitivity results were compared against theoretical predictions. Both case studies demonstrated the flexibility and accuracy of GATE and exemplified its potential benefits in protocol optimization and in system design.     

Numerical study of reflectance imaging using a parallel Monte Carlo method

Reflectance imaging of biological tissues with visible and near-infrared light has the significant potential to provide a noninvasive and safe imaging modality for diagnosis of dysplastic and malignant lesions in the superficial tissue layers. The difficulty in the extraction of optical and structural parameters lies in the lack of efficient methods for accurate modeling of light scattering in biological tissues of turbid nature. We present a parallel Monte Carlo method for accurate and efficient modeling of reflectance images from turbid tissue phantoms. A parallel Monte Carlo code has been developed with the message passing interface and evaluated on a computing cluster with 16 processing elements. The code was validated against the solutions of the radiative transfer equation on the bidirectional reflection and transmission functions. With this code we investigated numerically the dependence of reflectance image on the imaging system and phantom parameters. The contrasts of reflectance images were found to be nearly independent of the numerical aperture (NA) of the imaging camera despite the fact that reflectance depends on the NA. This enables efficient simulations of the reflectance images using an NA at 1.00. Using heterogeneous tissue phantoms with an embedded region simulating a lesion, we investigated the correlation between the reflectance image profile or contrast and the phantom parameters. It has been shown that the image contrast approaches 0 when the single-scattering albedos of the two regions in the heterogeneous phantoms become matched. Furthermore, a zone of detection has been demonstrated for determination of the thickness of the embedded region and optical parameters from the reflectance image profile and contrast. Therefore, the utility of the reflectance imaging method with visible and near-infrared light has been firmly established. We conclude from these results that the optical parameters of the embedded region can be determined inversely from reflectance images acquired with full-field illumination at multiple incident angles or multiple wavelengths.     

Multilayer imaging and compositional analysis of human male breast by laser reflectometry and Monte Carlo simulation

The laser backscattering from biological tissues depends on their composition and blood flow. The onset of abnormalities in tissues is associated with the change in composition at a specific location which may affect laser backscattering. The objective of this work is to study the point-to-point compositional variation of male breast tissues as this site has been prone to cancer development. The normalized backscattered intensity (NBI) profiles at various locations of human chest region of five subjects by multi-probe laser reflectometer are obtained. Based on these data the images of tissue composition, showing the point-to-point changes at various depths from the tissue surface, are reconstructed. The analysis of data shows that the maximum NBI variation is at the pectoralis major muscle and minimum variation is observed at the sternum. The optical parameters, based on the NBI data obtained for five human subjects, show the maximum increase in absorption (p < 0.0001) and minimum change in scattering (p < 0.0001) coefficients compared to that as observed at the sternum. Also the minimum absorption and maximum scattering coefficients are observed at the pectoralis major muscles. The regional variations of NBI and optical parameters further support these findings. The variations in the NBI and optical parameters may indicate the compositional change in tissues, which could be used for diagnostic and therapeutic applications of laser.     

Characterization of a novel micro-irradiator using Monte Carlo radiation transport simulations

Small animals are highly valuable resources for radiobiology research. While rodents have been widely used for decades, zebrafish embryos have recently become a very popular research model. However, unlike rodents, zebrafish embryos lack appropriate irradiation tools and methodologies. Therefore, the main purpose of this work is to use Monte Carlo radiation transport simulations to characterize dosimetric parameters, determine dosimetric sensitivity and help with the design of a new micro-irradiator capable of delivering irradiation fields as small as 1.0 mm in diameter. The system is based on a miniature x-ray source enclosed in a brass collimator with 3 cm diameter and 3 cm length. A pinhole of 1.0 mm diameter along the central axis of the collimator is used to produce a narrow photon beam. The MCNP5, Monte Carlo code, is used to study the beam energy spectrum, percentage depth dose curves, penumbra and effective field size, dose rate and radiation levels at 50 cm from the source. The results obtained from Monte Carlo simulations show that a beam produced by the miniature x-ray and the collimator system is adequate to totally or partially irradiate zebrafish embryos, cell cultures and other small specimens used in radiobiology research.     

Optimization of accelerator target and detector for portal imaging using Monte Carlo simulation and experiment

Megavoltage portal images suffer from poor quality compared to those produced with kilovoltage x-rays. Several authors have shown that the image quality can be improved by modifying the linear accelerator to generate more low-energy photons. This work addresses the problem of using Monte Carlo simulation and experiment to optimize the beam and detector combination to maximize image quality for a given patient thickness. A simple model of the whole imaging chain was developed for investigation of the effect of the target parameters on the quality of the image. The optimum targets (6 mm thick aluminium and 1.6 mm copper) were installed in an Elekta SL25 accelerator. The first beam will be referred to as A16 and the second as Cu1.6. A tissue-equivalent contrast phantom was imaged with the 6 MV standard photon beam and the experimental beams with standard radiotherapy and mammography film/screen systems. The arrangement with a thin Al target/mammography system improved the contrast from 1.4 cm bone in 5 cm water to 19% compared with 2% for the standard arrangement of a thick, high-Z target/radiotherapy verification system. The linac/phantom/detector system was simulated with the BEAM/EGS4 Monte Carlo code. Contrast calculated from the predicted images was in good agreement with the experiment (to within 2.5%). The use of MC techniques to predict images accurately, taking into account the whole imaging system, is a powerful new method for portal imaging system design optimization.     

Fast modelling of the collimator-detector response in Monte Carlo simulation of SPECT imaging using the angular response function

Interactions of incident photons with the collimator and detector, including septal penetration, scatter and x-ray fluorescence, are significant sources of image degradation in applications of SPECT including dual isotope imaging and imaging using radioisotopes that emit high- or medium-energy photons. Modelling these interactions using full Monte Carlo (MC) simulations is computationally very demanding. We present a new method based on the use of angular response functions (ARFs). The ARF is a function of the incident photon's direction and energy and represents the probability that a photon will either interact with or pass through the collimator, and be detected at the intersection of the photon's direction vector and the detection plane in an energy window of interest. The ARFs were pre-computed using full MC simulations of point sources that include propagation through the collimator-detector system. We have implemented the ARF method for use in conjunction with the SimSET/PHG MC code to provide fast modelling of both interactions in the patient and in the collimator-detector system. Validation results in the three cases studied show that there was good agreement between the projections generated using the ARF method and those from previously validated full MC simulations, but with hundred to thousand fold reductions in simulation time.     

Computer phantom study of brain PET glucose metabolism imaging using a rotating SPECT/PET camera

Positron emission tomography (PET) with [18F] fluoro-deoxy-glucose (FDG) provides information about glucose metabolism and is used to measure tissue glucose kinetics in the brain. The recent interest in hybrid SPECT/PET systems emerged as a practical approach to reduce the high cost of purchasing a dedicated ring-detector PET system. We have implemented interpolation methods for processing the projection data that could potentially reduce artifacts when reconstructing a dynamic imaging sequence in a PET study from a dual-head rotating SPECT/PET system. The computer simulations predict that parameter estimates from the dedicated PET system will be superior to results using the rotating camera system. However, the rotating camera system using projection interpolation may approach the accuracy of the dedicated PET system if the data noise is below 20%.     

Monte Carlo simulation of a protontherapy platform devoted to ocular melanoma

Patients with ocular melanoma have been treated since June 1991 at the medical cyclotron of the Centre Antoine Lacassagne (CAL). Positions and sizes of the ocular nozzle elements were initially defined based on experimental work, taking as a pattern functional existing facilities. Nowadays Monte Carlo (MC) calculation offers a tool to refine this geometry by adjusting size and place of beam modeling devices. Moreover, the MC tool is a useful way to calculate the dose and to evaluate the impact of secondary particles in the field of radiotherapy or radiation protection. Both LINAC and cyclotron producing x rays, electrons, protons, and neutrons are available in CAL, which suggests choosing MCNPX for its particle versatility. As a first step, the existing installation was input in MCNPX to check its aptitude to reproduce experimentally measured depth-dose profile, lateral profile, output-factor (OF), and absolute dose. The geometry was defined precisely and described from the last achromatic bending magnet of our proton beam line to the position of treated eyes. Relative comparisons of percentage depth-dose and lateral profiles, performed between measured data and simulations, show an agreement of the order of 2% in dose and 0.1 mm in range accuracy. These comparisons, carried out with and without beam-modifying device, yield results compatible to the required precision in ocular melanoma treatments, as long as adequate choices are made on MCNPX input decks for physics card. Absolute dose and OF issued from calculations and measurements were also compared. Results obtained for these two kinds of data, carried out in the simplified situation of an unmodulated beam, indicate that MC calculation could effectively complement measurements. These encouraging results are a large source of motivation to promote further studies, first in a new design of the ocular nozzle, and second in the analysis of the influence of beam-modifying devices attached to the final patient collimator, such as wedge or compensators, on dose values.