No Arabic abstract
In this paper we investigate, with a detailed Monte-Carlo simulation based on Geant4, the novel approach [Nucl. Instrum. Methods A588 (2008) 457] to 3D imaging with photon scattering. A monochromatic and well collimated gamma beam is used to illuminate the object to be imaged and the photons Compton scattered are detected by means of a surrounding germanium strip detector. The impact position and the energy of the photons are measured with high precision and the scattering position along the beam axis is calculated. We study as an application of this technique the case of brain imaging but the results can be applied as well to situations where a lighter object, with localized variations of density, is embedded in a denser container. We report here the attainable sensitivity in the detection of density variations as a function of the beam energy, the depth inside the object and size and density of the inclusions. Using a 600 keV gamma beam, for an inclusion with a density increase of 30% with respect to the sorrounding tissue and thickness along the beam of 5 mm, we obtain at midbrain position a resolution of about 2 mm and a contrast of 12%. In addition the simulation indicates that for the same gamma beam energy a complete brain scan would result in an effective dose of about 1 mSv.
The Monte Carlo simulation of the electron transport through thin slabs is studied with five general purpose codes: PENELOPE, GEANT3, GEANT4, EGSnrc and MCNPX. The different material foils analyzed in the old experiments of Kulchitsky and Latyshev [Phys. Rev. 61 (1942) 254-266] and Hanson et al. [Phys. Rev. 84 (1951) 634-637] are used to perform the comparison between the Monte Carlo codes. Non-negligible differences are observed in the angular distributions of the transmitted electrons obtained with the some of the codes. The experimental data are reasonably well described by EGSnrc, PENELOPE (v. 2005) and GEANT4. A general good agreement is found for EGSnrc and GEANT4 in all the cases analyzed.
Charged Particle Therapy is a technique for cancer treatment that exploits hadron beams, mostly protons and carbons. A critical issue is the monitoring of the dose released by the beam to the tumor and to the surrounding tissues. We present the design of a new tracking device for monitoring on-line the dose in ion therapy through the detection of secondary charged particles produced by the beam interactions in the patient tissues. In fact, the charged particle emission shape can be correlated with the spatial dose release and the Bragg peak position. The detector uses the information provided by 12 layers of scintillating fibers followed by a plastic scintillator and a small calorimeter made of a pixelated Lutetium Fine Silicate crystal. Simulations have been performed to evaluate the achievable spatial resolution and a possible application of the device for the monitoring of the dose profile in a real treatment is presented.
Liquid water has been proved to be an excellent medium for specimen structure imaging by a scanning electron microscope. Knowledge of electron-water interaction physics and particularly the secondary electron yield is essential to the interpretation of the imaging contrast. However, very little is known up to now experimentally on the low energy electron interaction with liquid water because of certain practical limitations. It is then important to gain some useful information about electron emission from water by a Monte Carlo (MC) simulation technique that can numerically model electron transport trajectories in water. In this study, we have performed MC simulations of electron emission from liquid water in the primary energy range of 50 eV-30 keV by using two different codes, i.e. a classical MC (CMC) code developed in our laboratory and the Geant4-DNA (G4DNA) code. The calculated secondary electron yield and electron backscattering coefficient are compared with experimental results wherever applicable to verify the validity of physical models for the electron-water interaction. The secondary electron yield vs. primary energy curves calculated by the two codes present the same generic curve shape as that of metals but in rather different absolute values. G4DNA yields the underestimated absolute values due to the application of one step thermalization model by setting a cutoff energy at 7.4 eV so that the low energy losses due to phonon excitations are omitted. Our CMC calculation of secondary electron yield is closer to the experimental data and the energy distribution is reasonable. It is concluded that a full dielectric function data at low energy loss values below 7.4 eV shall be employed in G4DNA model for the modeling of low energy electrons.
Image reconstruction in Single Photon Emission Computed Tomography (SPECT) is affected by physical effects such as photon attenuation, Compton scatter and detector response. These effects can be compensated for by modeling the corresponding spread of photons in 3D within the system matrix used for tomographic reconstruction. The fully 3D Monte Carlo (F3DMC) reconstruction technique consists in calculating this system matrix using Monte Carlo simulations. The inverse problem of tomographic reconstruction is then solved using conventional iterative algorithms such as maximum likelihood expectation maximization (MLEM). Although F3DMC has already shown promising results, its use is currently limited by two major issues: huge size of the fully 3D system matrix and long computation time required for calculating a robust and accurate system matrix. To address these two issues, we propose to calculate the F3DMC system matrix using a spatial sampling matching the functional regions to be reconstructed. In this approach, different regions of interest can be reconstructed with different spatial sampling. For instance, a single value is reconstructed for a functional region assumed to contain uniform activity. To assess the value of this approach, Monte Carlo simulations have been performed using GATE. Results suggest that F3DMC reconstruction using functional regions improves quantitative accuracy compared to the F3DMC reconstruction method proposed so far. In addition, it considerably reduces disk space requirement and duration of the simulations needed to estimate the system matrix. The concept of functional regions might therefore make F3DMC reconstruction practically feasible.
Precise quantitative delineation of tumor hypoxia is essential in radiation therapy treatment planning to improve the treatment efficacy by targeting hypoxic sub-volumes. We developed a combined imaging system of positron emission tomography (PET) and electron para-magnetic resonance imaging (EPRI) of molecular oxygen to investigate the accuracy of PET imaging in assessing tumor hypoxia. The PET/EPRI combined imaging system aims to use EPRI to precisely measure the oxygen partial pressure in tissues. This will evaluate the validity of PET hypoxic tumor imaging by (near) simultaneously acquired EPRI as ground truth. The combined imaging system was constructed by integrating a small animal PET scanner (inner ring diameter 62 mm and axial field of view 25.6 mm) and an EPRI subsystem (field strength 25 mT and resonant frequency 700 MHz). The compatibility between the PET and EPRI subsystems were tested with both phantom and animal imaging. Hypoxic imaging on a tumor mouse model using $^{18}$F-fluoromisonidazole radio-tracer was conducted with the developed PET/EPRI system. We report the development and initial imaging results obtained from the PET/EPRI combined imaging system.