Monte Carlo simulation is an essential tool in emission tomography that can assist in the design of new medical imaging devices, the optimization of acquisition protocols and the development or assessment of image reconstruction algorithms and correction techniques. GATE, the Geant4 Application for Tomographic Emission, encapsulates the Geant4 libraries to achieve a modular, versatile, scripted simulation toolkit adapted to the field of nuclear medicine. In particular, GATE allows the description of time-dependent phenomena such as source or detector movement, and source decay kinetics. This feature makes it possible to simulate time curves under realistic acquisition conditions and to test dynamic reconstruction algorithms. This paper gives a detailed description of the design and development of GATE by the OpenGATE collaboration, whose continuing objective is to improve, document and validate GATE by simulating commercially available imaging systems for PET and SPECT. Large effort is also invested in the ability and the flexibility to model novel detection systems or systems still under design. A public release of GATE licensed under the GNU Lesser General Public License can be downloaded at http:/www-lphe.epfl.ch/GATE/. Two benchmarks developed for PET and SPECT to test the installation of GATE and to serve as a tutorial for the users are presented. Extensive validation of the GATE simulation platform has been started, comparing simulations and measurements on commercially available acquisition systems. References to those results are listed. The future prospects towards the gridification of GATE and its extension to other domains such as dosimetry are also discussed.

The feasibility and limits in performing tomographic bioluminescence imaging with a combined optical-PET (OPET) system was explored by simulating its image formation process. A micro-MRI based virtual mouse phantom was assigned appropriate tissue optical properties to each of its segmented internal organs at wavelengths spanning the emission spectrum of the firefly luciferase at 37 °C. The TOAST finite-element code was employed to simulate the diffuse transport of photons emitted from bioluminescence sources in the mouse. OPET measurements were simulated for singlepoint, two-point and distributed bioluminescence sources located in different organs such as the liver, the kidneys and the gut. An expectation maximization code was employed to recover the intensity and location of these simulated sources. It was found that spectrally resolved measurements were necessary in order to perform tomographic bioluminescence imaging. The true location of emission sources could be recovered if the mouse background optical properties were known a priori. Assumption of a homogeneous optical property background proved inadequate for describing photon transport in optically heterogeneous tissues and lead to inaccurate source localization in the reconstructed images. The simulation results pointed out specific methodological challenges that need to be addressed before a practical implementation of OPET-based bioluminescence tomography is achieved.

Attenuation correction is one of the important corrections required for quantitative positron emission tomography (PET). This work will compare the quantitative accuracy of attenuation correction using a simple global scale factor with traditional transmission-based methods acquired either with a small animal PET or a small animal x-ray computed tomography (CT) scanner. Two phantoms (one mouse-sized and one rat-sized) and two animal subjects (one mouse and one rat) were scanned in CTI Concorde Microsystem's microPET R FocusTM for emission and transmission data and in ImTek's MicroCAT TM II for transmission data. PET emission image values were calibrated against a scintillation well counter. Results indicate that the scale factor method of attenuation correction places the average measured activity concentration about the expected value, without correcting for the cupping artefact from attenuation. Noise analysis in the phantom studies with the PET-based method shows that noise in the transmission data increases the noise in the corrected emission data. The CT-based method was accurate and delivered low-noise images suitable for both PET data correction and PET tracer localization.

Optical/IR images of transition disks (TDs) have revealed deep intensity decrements in the rings of HAeBes HD 142527 and HD 100453, that can be interpreted as shadowing from sharply tilted inner disks, such that the outer disks are directly exposed to stellar light. Here we report similar dips in SPHERE+IRDIS differential polarized imaging (DPI) of TTauri DoAr 44. With a fairly axially symmetric ring in the sub mm radio continuum, DoAr 44 is likely also a warped system. We constrain the warp geometry by comparing radiative transfer predictions with the DPI data in H band (Q φ (H)) and with a re-processing of archival 336 GHz ALMA observations. The observed DPI shadows have coincident radio counterparts, but the intensity drops are much deeper in Q φ (H) (∼88%), compared to the shallow drops at 336 GHz (∼24%). Radiative transfer predictions with an inner disk tilt of ∼ 30 ± 5 deg approximately account for the observations. ALMA long-baseline observations should allow the observation of the warped gas kinematics inside the cavity of DoAr 44.