Report on G4-Med, a Geant4 benchmarking system for medical physics applications developed by the Geant4 Medical Simulation Benchmarking Group.

BACKGROUND: Geant4 is a Monte Carlo code extensively used in medical physics for a wide range of applications, such as dosimetry, micro- and nanodosimetry, imaging, radiation protection, and nuclear medicine. Geant4 is continuously evolving, so it is crucial to have a system that benchmarks this Monte Carlo code for medical physics against reference data and to perform regression testing. AIMS: To respond to these needs, we developed G4-Med, a benchmarking and regression testing system of Geant4 for medical physics. MATERIALS AND METHODS: G4-Med currently includes 18 tests. They range from the benchmarking of fundamental physics quantities to the testing of Monte Carlo simulation setups typical of medical physics applications. Both electromagnetic and hadronic physics processes and models within the prebuilt Geant4 physics lists are tested. The tests included in G4-Med are executed on the CERN computing infrastructure via the use of the geant-val web application, developed at CERN for Geant4 testing. The physical observables can be compared to reference data for benchmarking and to results of previous Geant4 versions for regression testing purposes. RESULTS: This paper describes the tests included in G4-Med and shows the results derived from the benchmarking of Geant4 10.5 against reference data. DISCUSSION: Our results indicate that the Geant4 electromagnetic physics constructor G4EmStandardPhysics_option4 gives a good agreement with the reference data for all the tests. The QGSP_BIC_HP physics list provided an overall adequate description of the physics involved in hadron therapy, including proton and carbon ion therapy. New tests should be included in the next stage of the project to extend the benchmarking to other physical quantities and application scenarios of interest for medical physics. CONCLUSION: The results presented and discussed in this paper will aid users in tailoring physics lists to their particular application.

Glandular dose in breast computed tomography with synchrotron radiation.

The purpose of this work is to provide an evaluation of the mean glandular dose (MGD) for breast computed tomography (CT) with synchrotron radiation in an axial scanning configuration with a partial or total organ volume irradiation, for the in vivo program of breast CT ongoing at the ELETTRA facility (Trieste, Italy). A Geant4 Monte Carlo code was implemented, simulating the photon irradiation from a synchrotron radiation source in the energetic range from 8 to 50 keV with 1 keV intervals, to evaluate the MGD. The code was validated with literature data, in terms of mammographic normalized glandular dose coefficients (DgN) and with ad hoc experimental data, in terms of computed tomography dose index (CTDI). Simulated cylindrical phantoms of different sizes (diameter at phantom base 8, 10, 12, 14 or 16 cm, axial length 1.5 times the radius) and glandular fraction by weight (0%, 14.3%, 25%, 50%, 75% and 100%) were implemented into the code. The validation of the code shows an excellent agreement both with previously published work and in terms of DgN and CDTI measurements. The implemented simulations show a dependence of the glandular dose estimate on the vertical dimension of the irradiated zone when a partial organ irradiation was implemented. Specific normalized coefficients for calculating the MGD to the whole breast or to the single irradiated slice were reported.

Energy response of GR-200A thermoluminescence dosemeters to 60Co and to monoenergetic synchrotron radiation in the energy range 28-40 keV.

The response of LiF:Mg,Cu,P thermoluminescence dosemeters (type GR-200A) to monoenergetic radiation of energy 28, 35, 38 and 40 keV was evaluated with respect to irradiation with a calibrated (60)Co gamma-ray source. High-precision measurements of the relative air kerma response performed at the SYRMEP beamline of the ELETTRA synchrotron radiation facility (Trieste, Italy) showed a significant deviation of the average response to low-energy X-rays from that to (60)Co, with an over-response from 6 % (at 28 keV) to 22 % (at 40 keV). These data are not consistent with literature data for these dosemeters, where model predictions gave deviation from unity of the relative air kerma response of about 10 %. The authors conclude for the need of additional determinations of the low-energy relative response of GR-200A dosemeters, covering a wider range of monoenergetic energies sampled at a fine energy step, as planned in future experiments by their group at the ELETTRA facility.

Imaging performance of phase-contrast breast computed tomography with synchrotron radiation and a CdTe photon-counting detector.

PURPOSE: Within the SYRMA-CT collaboration based at the ELETTRA synchrotron radiation (SR) facility the authors investigated the imaging performance of the phase-contrast computed tomography (CT) system dedicated to monochromatic in vivo 3D imaging of the female breast, for breast cancer diagnosis. METHODS: Test objects were imaged at 38keV using monochromatic SR and a high-resolution CdTe photon-counting detector. Signal and noise performance were evaluated using modulation transfer function (MTF) and noise power spectrum. The analysis was performed on the images obtained with the application of a phase retrieval algorithm as well as on those obtained without phase retrieval. The contrast to noise ratio (CNR) and the capability of detecting test microcalcification clusters and soft masses were investigated. RESULTS: For a voxel size of (60µm)(3), images without phase retrieval showed higher spatial resolution (6.7mm(-1) at 10% MTF) than corresponding images with phase retrieval (2.5mm(-1)). Phase retrieval produced a reduction of the noise level and an increase of the CNR by more than one order of magnitude, compared to raw phase-contrast images. Microcalcifications with a diameter down to 130µm could be detected in both types of images. CONCLUSIONS: The investigation on test objects indicates that breast CT with a monochromatic SR source is technically feasible in terms of spatial resolution, image noise and contrast, for in vivo 3D imaging with a dose comparable to that of two-view mammography. Images obtained with the phase retrieval algorithm showed the best performance in the trade-off between spatial resolution and image noise.

Towards breast tomography with synchrotron radiation at Elettra: first images.

The aim of the SYRMA-CT collaboration is to set-up the first clinical trial of phase-contrast breast CT with synchrotron radiation (SR). In order to combine high image quality and low delivered dose a number of innovative elements are merged: a CdTe single photon counting detector, state-of-the-art CT reconstruction and phase retrieval algorithms. To facilitate an accurate exam optimization, a Monte Carlo model was developed for dose calculation using GEANT4. In this study, high isotropic spatial resolution (120 µm)(3) CT scans of objects with dimensions and attenuation similar to a human breast were acquired, delivering mean glandular doses in the range of those delivered in clinical breast CT (5-25 mGy). Due to the spatial coherence of the SR beam and the long distance between sample and detector, the images contain, not only absorption, but also phase information from the samples. The application of a phase-retrieval procedure increases the contrast-to-noise ratio of the tomographic images, while the contrast remains almost constant. After applying the simultaneous algebraic reconstruction technique to low-dose phase-retrieved data sets (about 5 mGy) with a reduced number of projections, the spatial resolution was found to be equal to filtered back projection utilizing a four fold higher dose, while the contrast-to-noise ratio was reduced by 30%. These first results indicate the feasibility of clinical breast CT with SR.

GEANT4 for breast dosimetry: parameters optimization study.

Mean glandular dose (MGD) is the main dosimetric quantity in mammography. MGD evaluation is obtained by multiplying the entrance skin air kerma (ESAK) by normalized glandular dose (DgN) coefficients. While ESAK is an empirical quantity, DgN coefficients can only be estimated with Monte Carlo (MC) methods. Thus, a MC parameters benchmark is needed for effectively evaluating DgN coefficients. GEANT4 is a MC toolkit suitable for medical purposes that offers to the users several computational choices. In this work we investigate the GEANT4 performances testing the main PhysicsLists for medical applications. Four electromagnetic PhysicsLists were implemented: the linear attenuation coefficients were calculated for breast glandularity 0%, 50%, 100% in the energetic range 8-50 keV and DgN coefficients were evaluated. The results were compared with published data. Fit equations for the estimation of the G-factor parameter, introduced by the literature for converting the dose delivered in the heterogeneous medium to that in the glandular tissue, are proposed and the application of this parameter interaction-by-interaction or retrospectively is discussed. G4EmLivermorePhysicsList shows the best agreement for the linear attenuation coefficients both with theoretical values and published data. Moreover, excellent correlation factor (r2>0.99) is found for the DgN coefficients with the literature. The final goal of this study is to identify, for the first time, a benchmark of parameters that could be useful for future breast dosimetry studies with GEANT4.