Z. Francis

THE Geant4-DNA project

The Geant4-DNA project proposes to develop an open-source simulation software based and fully included in the general-purpose Geant4 Monte-Carlo simulation toolkit. The main objective of this software is to simulate biological damages induced by ionizing radiations at the cellular and sub-cellular scale. This project was originally initiated by the European Space Agency for the prediction of the deleterious effects of radiations that may affect astronauts during future long duration space exploration missions. In this paper, the Geant4-DNA collaboration presents an overview of the whole on-going project, including its most recent developments that are available in the Geant4 toolkit since December 2009 (release 9.3), as well as an illustration example simulating the direct irradiation of a biological chromatin fiber. Expected extensions involving several research domains, such as particle physics, chemistry and cellular and molecular biology, within a fully interdisciplinary activity of the Geant4 collaboration are also discussed.

Geant4 Physics Processes for Microdosimetry Simulation: Design Foundation and Implementation of the First Set of Models

New physical processes specific for microdosimetry simulation are under development in the Geant4 Low Energy Electromagnetic package. The first set of models implemented for this purpose cover the interactions of electrons, protons and light ions in liquid water; they address a physics domain relevant to the simulation of radiation effects in biological systems, where water represents an important component. The design developed for effectively handling particle interactions down to a low energy scale and the physics models implemented in the first public release of the software are described.

Track structure modeling in liquid water: A review of the Geant4-DNA very low energy extension of the Geant4 Monte Carlo simulation toolkit

Understanding the fundamental mechanisms involved in the induction of biological damage by ionizing radiation remains a major challenge of todays radiobiology research. The Monte Carlo simulation of physical, physicochemical and chemical processes involved may provide a powerful tool for the simulation of early damage induction. The Geant4-DNA extension of the general purpose Monte Carlo Geant4 simulation toolkit aims to provide the scientific community with an open source access platform for the mechanistic simulation of such early damage. This paper presents the most recent review of the Geant4-DNA extension, as available to Geant4 users since June 2015 (release 10.2 Beta). In particular, the review includes the description of new physical models for the description of electron elastic and inelastic interactions in liquid water, as well as new examples dedicated to the simulation of physicochemical and chemical stages of water radiolysis. Several implementations of geometrical models of biological targets are presented as well, and the list of Geant4-DNA examples is described.

Molecular scale track structure simulations in liquid water using the Geant4-DNA Monte-Carlo processes

This paper presents a study of energy deposits induced by ionising particles in liquid water at the molecular scale. Particles track structures were generated using the Geant4-DNA processes of the Geant4 Monte-Carlo toolkit. These processes cover electrons (0.025 eV-1 MeV), protons (1 keV-100 MeV), hydrogen atoms (1 keV-100 MeV) and alpha particles (10 keV-40 MeV) including their different charge states. Electron ranges and lineal energies for protons were calculated in nanometric and micrometric volumes.

Physical models implemented in the GEANT4-DNA extension of the GEANT-4 toolkit for calculating initial radiation damage at the molecular level.

The ROSIRIS project aims to study the radiobiology of integrated systems for medical treatment optimisation using ionising radiations and evaluate the associated risk. In the framework of this project, one research focus is the interpretation of the initial radio-induced damage in DNA created by ionising radiation (and detected by γ-H2AX foci analysis) from the track structure of the incident particles. In order to calculate the track structure of ionising particles at a nanometric level, the Geant4 Monte Carlo toolkit was used. Geant4 (Object Oriented Programming Architecture in C++) offers a common platform, available free to all users and relatively easy to use. Nevertheless, the current low-energy threshold for electromagnetic processes in GEANT4 is set to 1 keV (250 eV using the Livermore processes), which is an unsuitable value for nanometric applications. To lower this energy threshold, the necessary interaction processes and models were identified, and the corresponding available cross sections collected from the literature. They are mostly based on the plane-wave Born approximation (first Born approximation, or FBA) for inelastic interactions and on semi-empirical models for energies where the FBA fails (at low energies). In this paper, the extensions that have been introduced into the 9.3 release of the Geant4 toolkit are described, the so-called Geant4-DNA extension, including a set of processes and models adapted in this study and permitting the simulation of electron (8 eV-1 MeV), proton (100 eV-100 MeV) and alpha particle (1 keV-10 MeV) interactions in liquid water.

Monte Carlo simulation of energy-deposit clustering for ions of the same LET in liquid water.

This work presents a Monte Carlo study of energy depositions due to protons, alpha particles and carbon ions of the same linear-energy-transfer (LET) in liquid water. The corresponding track structures were generated using the Geant4-DNA toolkit, and the energy deposition spatial distributions were analyzed using an adapted version of the DBSCAN clustering algorithm. Combining the Geant4 simulations and the clustering algorithm it was possible to compare the quality of the different radiation types. The ratios of clustered and single energy depositions are shown versus particle LET and frequency-mean lineal energies. The estimated effect of these types of radiation on biological tissues is then discussed by comparing the results obtained for different particles with the same LET.

Technical Note: Improvements in GEANT4 energy-loss model and the effect on low-energy electron transport in liquid water

PURPOSE The geant4-DNA physics models are upgraded by a more accurate set of electron cross sections for ionization and excitation in liquid water. The impact of the new developments on low-energy electron transport simulations by the geant4 Monte Carlo toolkit is examined for improving its performance in dosimetry applications at the subcellular and nanometer level. METHODS The authors provide an algorithm for an improved implementation of the Emfietzoglou model dielectric response function of liquid water used in the geant4-DNA existing model. The algorithm redistributes the imaginary part of the dielectric function to ensure a physically motivated behavior at the binding energies, while retaining all the advantages of the original formulation, e.g., the analytic properties and the fulfillment of the f-sum-rule. In addition, refinements in the exchange and perturbation corrections to the Born approximation used in the geant4-DNA existing model are also made. RESULTS The new ionization and excitation cross sections are significantly different from those of the geant4-DNA existing model. In particular, excitations are strongly enhanced relative to ionizations, resulting in higher W-values and less diffusive dose-point-kernels at sub-keV electron energies. CONCLUSIONS An improved energy-loss model for the excitation and ionization of liquid water by low-energy electrons has been implemented in geant4-DNA. The suspiciously low W-values and the unphysical long tail in the dose-point-kernel have been corrected owing to a different partitioning of the dielectric function.

A study on the physicochemical parameters for Penicillium expansum growth and patulin production: effect of temperature, pH, and water activity

Abstract Penicillium expansum is among the most ubiquitous fungi disseminated worldwide, that could threaten the fruit sector by secreting patulin, a toxic secondary metabolite. Nevertheless, we lack sufficient data regarding the growth and the toxigenesis conditions of this species. This work enables a clear differentiation between the favorable conditions to the P. expansum growth and those promising for patulin production. A mathematical model allowing the estimation of the P. expansum growth rate according to temperature, a W, and pH, was also developed. An optimal growth rate of 0.92 cm/day was predicted at 24°C with pH level of 5.1 and high a W level of 0.99. The models predictive capability was tested successfully on artificial contaminated apples. This model could be exploited by apple growers and the industrialists of fruit juices in order to predict the development of P. expansum during storage and apple processing.

Carbon ion fragmentation effects on the nanometric level behind the Bragg peak depth

In this study, fragmentation yields of carbon therapy beams are estimated using the Geant4 simulation toolkit version 9.5. Simulations are carried out in a step-by-step mode using the Geant4-DNA processes for each of the major contributing fragments. The energy of the initial beam is taken 400 MeV amu(-1) as this is the highest energy, which is used for medical accelerators and this would show the integral role of secondary contributions in radiotherapy irradiations. The obtained results showed that 64% of the global dose deposition is initiated by carbon ions, while up to 36% is initiated by the produced fragments including all their isotopes. The energy deposition clustering yields of each of the simulated fragments are then estimated using the DBSCAN clustering algorithm and they are compared to the yields of the incident primary beam.

The Influence of DNA Configuration on the Direct Strand Break Yield

Purpose. To study the influence of DNA configuration on the direct damage yield. No indirect effect has been accounted for. Methods. The GEANT4-DNA code was used to simulate the interactions of protons and alpha particles with geometrical models of the A-, B-, and Z-DNA configurations. The direct total, single, and double strand break yields and site-hit probabilities were determined. Certain features of the energy deposition process were also studied. Results. A slight increase of the site-hit probability as a function of the incident particle linear energy transfer was found for each DNA configuration. Each DNA form presents a well-defined site-hit probability, independently of the particle linear energy transfer. Approximately 70% of the inelastic collisions and ~60% of the absorbed dose are due to secondary electrons. These fractions are slightly higher for protons than for alpha particles at the same incident energy. Conclusions. The total direct strand break yield for a given DNA form depends weakly on DNA conformation topology. This yield is practically determined by the target volume of the DNA configuration. However, the double strand break yield increases with the packing ratio of the DNA double helix; thus, it depends on the DNA conformation.