The GEANT4 toolkit capability in the hadron therapy field: simulation of a transport beam line

At Laboratori Nazionali del Sud of the Instituto Nazionale di Fisica Nucleare of Catania (Sicily, Italy), the first Italian hadron therapy facility named CATANA (Centro di AdroTerapia ed Applicazioni Nucleari Avanzate) has been realized. Inside CATANA 62 MeV proton beams, accelerated by a superconducting cyclotron, are used for the radiotherapeutic treatments of some types of ocular tumours. Therapy with hadron beams still represents a pioneer technique, and only a few centers worldwide can provide this advanced specialized cancer treatment. On the basis of the experience so far gained, and considering the future hadron-therapy facilities to be developed (Rinecker, Munich Germany, Heidelberg/GSI, Darmstadt, Germany, PSI Villigen, Switzerland, CNAO, Pavia, Italy, Centro di Adroterapia, Catania, Italy) we decided to develop a Monte Carlo application based on the GEANT4 toolkit, for the design, the realization and the optimization of a proton-therapy beam line. Another feature of our project is to provide a general tool able to study the interactions of hadrons with the human tissue and to test the analytical-based treatment planning systems actually used in the routine practice. All the typical elements of a hadron-therapy line, such as diffusers, range shifters, collimators and detectors were modelled. In particular, we simulated the Markus type ionization chamber and a Gaf Chromic film as dosimeters to reconstruct the depth (Bragg peak and Spread Out Bragg Peak) and lateral dose distributions, respectively. We validated our simulated detectors comparing the results with the experimental data available in our facility

Correction factors to convert microdosimetry measurements in silicon to tissue in \u3csup\u3e12\u3c/sup\u3eC ion therapy

Silicon microdosimetry is a promising technology for heavy ion therapy (HIT) quality assurance, because of its sub-mm spatial resolution and capability to determine radiation effects at a cellular level in a mixed radiation field. A drawback of silicon is not being tissue-equivalent, thus the need to convert the detector response obtained in silicon to tissue. This paper presents a method for converting silicon microdosimetric spectra to tissue for a therapeutic 12C beam, based on Monte Carlo simulations. The energy deposition spectra in a 10 μm sized silicon cylindrical sensitive volume (SV) were found to be equivalent to those measured in a tissue SV, with the same shape, but with dimensions scaled by a factor κ equal to 0.57 and 0.54 for muscle and water, respectively. A low energy correction factor was determined to account for the enhanced response in silicon at low energy depositions, produced by electrons. The concept of the mean path length (lPath) to calculate the lineal energy was introduced as an alternative to the mean chord length (l) because it was found that adopting Cauchy\u27s formula for the (l) was not appropriate for the radiation field typical of HIT as it is very directional (lPath) can be determined based on the peak of the lineal energy distribution produced by the incident carbon beam. Furthermore it was demonstrated that the thickness of the SV along the direction of the incident 12C ion beam can be adopted as (lPath). The tissue equivalence conversion method and (lPath) were adopted to determine the RBE10, calculated using a modified microdosimetric kinetic model, applied to the microdosimetric spectra resulting from the simulation study. Comparison of the RBE10 along the Bragg peak to experimental TEPC measurements at HIMAC, NIRS, showed good agreement. Such agreement demonstrates the validity of the developed tissue equivalence correction factors and of the determination of (lPath)

Comparison of phantom materials for use in quality assurance of microbeam radiation therapy

Microbeam radiation therapy (MRT) is a promising radiotherapy modality that uses arrays of spatially fractionated micrometre-sized beams of synchrotron radiation to irradiate tumours. Routine dosimetry quality assurance (QA) prior to treatment is necessary to identify any changes in beam condition from the treatment plan, and is undertaken using solid homogeneous phantoms. Solid phantoms are designed for, and routinely used in, megavoltage X-ray beam radiation therapy. These solid phantoms are not necessarily designed to be water-equivalent at low X-ray energies, and therefore may not be suitable for MRT QA. This work quantitatively determines the most appropriate solid phantom to use in dosimetric MRT QA. Simulated dose profiles of various phantom materials were compared with those calculated in water under the same conditions. The phantoms under consideration were RMI457 Solid Water (Gammex-RMI, Middleton, WI, USA), Plastic Water (CIRS, Norfolk, VA, USA), Plastic Water DT (CIRS, Norfolk, VA, USA), PAGAT (CIRS, Norfolk, VA, USA), RW3 Solid Phantom (PTW Freiburg, Freiburg, Germany), PMMA, Virtual Water (Med-Cal, Verona, WI, USA) and Perspex. RMI457 Solid Water and Virtual Water were found to be the best approximations for water in MRT dosimetry (within ±3% deviation in peak and 6% in valley). RW3 and Plastic Water DT approximate the relative dose distribution in water (within ±3% deviation in the peak and 5% in the valley). PAGAT, PMMA, Perspex and Plastic Water are not recommended to be used as phantoms for MRT QA, due to dosimetric discrepancies greater than 5%

Validation of a Monte Carlo simulation for Microbeam Radiation Therapy on the Imaging and Medical Beamline at the Australian Synchrotron

Microbeam Radiation Therapy (MRT) is an emerging cancer treatment modality characterised by the use of high-intensity synchrotron-generated x-rays, spatially fractionated by a multi-slit collimator (MSC), to ablate target tumours. The implementation of an accurate treatment planning system, coupled with simulation tools that allow for independent verification of calculated dose distributions are required to ensure optimal treatment outcomes via reliable dose delivery. In this article we present data from the first Geant4 Monte Carlo radiation transport model of the Imaging and Medical Beamline at the Australian Synchrotron. We have developed the model for use as an independent verification tool for experiments in one of three MRT delivery rooms and therefore compare simulation results with equivalent experimental data. The normalised x-ray spectra produced by the Geant4 model and a previously validated analytical model, SPEC, showed very good agreement using wiggler magnetic field strengths of 2 and 3 T. However, the validity of absolute photon flux at the plane of the Phase Space File (PSF) for a fixed number of simulated electrons was unable to be established. This work shows a possible limitation of the G4SynchrotronRadiation process to model synchrotron radiation when using a variable magnetic field. To account for this limitation, experimentally derived normalisation factors for each wiggler field strength determined under reference conditions were implemented. Experimentally measured broadbeam and microbeam dose distributions within a Gammex RMI457 Solid Water® phantom were compared to simulated distributions generated by the Geant4 model. Simulated and measured broadbeam dose distributions agreed within 3% for all investigated configurations and measured depths. Agreement between the simulated and measured microbeam dose distributions agreed within 5% for all investigated configurations and measured depths

BrachyView, A novel inbody imaging system for HDR prostate brachytherapy: Design and Monte Carlo feasibility study

Purpose: High dose rate (HDR) brachytherapy is a form of radiation therapy for treating prostate cancer whereby a high activity radiation source is moved between predefined positions inside applicators inserted within the treatment volume. Accurate positioning of the source is essential in delivering the desired dose to the target area while avoiding radiation injury to the surrounding tissue. In this paper, HDR BrachyView, a novel inbody dosimetric imaging system for real time monitoring and verification of the radioactive seed position in HDR prostate brachytherapy treatment is introduced. The current prototype consists of a 15 × 60 mm2 silicon pixel detector with a multipinhole tungsten collimator placed 6.5 mm above the detector. Seven identical pinholes allow full imaging coverage of the entire treatment volume. The combined pinhole and pixel sensor arrangement is geometrically designed to be able to resolve the three-dimensional location of the source. The probe may be rotated to keep the whole prostate within the transverse plane. The purpose of this paper is to demonstrate the efficacy of the design through computer simulation, and to estimate the accuracy in resolving the source position (in detector plane and in 3D space) as part of the feasibility study for the BrachyView project. Methods: Monte Carlo simulations were performed using the GEANT4 radiation transport model, with a 192Ir source placed in different locations within a prostate phantom. A geometrically accurate model of the detector and collimator were constructed. Simulations were conducted with a single pinhole to evaluate the pinhole design and the signal to background ratio obtained. Second, a pair of adjacent pinholes were simulated to evaluate the error in calculated source location. Results: Simulation results show that accurate determination of the true source position is easily obtainable within the typical one second source dwell time. The maximum error in the estimated projection position was found to be 0.95 mm in the imaging (detector) plane, resulting in a maximum source positioning estimation error of 1.48 mm. Conclusions: HDR BrachyView is a feasible design for real-time source tracking in HDR prostate brachytherapy. It is capable of resolving the source position within a subsecond dwell time. In combination with anatomical information obtained from transrectal ultrasound imaging, HDR BrachyView adds a significant quality assurance capability to HDR brachytherapy treatment systems. © 2013 American Association of Physicists in Medicine

Thin Silicon Microdosimeter utilizing 3D MEMS Fabrication Technology: Charge Collection Study and its application in mixed radiation fields

New 10-μm-thick silicon microdosimeters utilizing 3-D technology have been developed and investigated in this paper. The TCAD simulations were carried out to understand the electrical properties of the microdosimeters\u27 design. A charge collection study of the devices was performed using 5.5-MeV He2+ ions which were raster scanned over the surface of the detectors and the charge collection median energy maps were obtained and the detection yield was also evaluated. The devices were tested in a 290 MeV/u carbon ion beam at the Heavy Ion Medical Accelerator in Chiba (HIMAC) in Japan. Based on the microdosimetric measurements, the quality factor and dose equivalent out of field were obtained in a mixed radiation field mimicking the radiation environment for spacecraft in deep space

Comparative study of alternative Geant4 hadronic ion inelastic physics models for prediction of positron-emitting radionuclide production in carbon and oxygen ion therapy

© 2019 Commonwealth of Australia, Australian Nuclear Science and Technology Organisation, ANSTO.. The distribution of fragmentation products predicted by Monte Carlo simulations of heavy ion therapy depend on the hadronic physics model chosen in the simulation. This work aims to evaluate three alternative hadronic inelastic fragmentation physics options available in the Geant4 Monte Carlo radiation physics simulation framework to determine which model most accurately predicts the production of positron-emitting fragmentation products observable using in-beam PET imaging. Fragment distributions obtained with the BIC, QMD, and INCL + + physics models in Geant4 version 10.2.p03 are compared to experimental data obtained at the HIMAC heavy-ion treatment facility at NIRS in Chiba, Japan. For both simulations and experiments, monoenergetic beams are applied to three different block phantoms composed of gelatin, poly(methyl methacrylate) and polyethylene. The yields of the positron-emitting nuclei 11C, 10C and 15O obtained from simulations conducted with each model are compared to the experimental yields estimated by fitting a multi-exponential radioactive decay model to dynamic PET images using the normalised mean square error metric in the entrance, build up/Bragg peak and tail regions. Significant differences in positron-emitting fragment yield are observed among the three physics models with the best overall fit to experimental 12C and 16O beam measurements obtained with the BIC physics model

GEANT4 : a simulation toolkit

Abstract Geant4 is a toolkit for simulating the passage of particles through matter. It includes a complete range of functionality including tracking, geometry, physics models and hits. The physics processes offered cover a comprehensive range, including electromagnetic, hadronic and optical processes, a large set of long-lived particles, materials and elements, over a wide energy range starting, in some cases, from 250 eV and extending in others to the TeV energy range. It has been designed and constructed to expose the physics models utilised, to handle complex geometries, and to enable its easy adaptation for optimal use in different sets of applications. The toolkit is the result of a worldwide collaboration of physicists and software engineers. It has been created exploiting software engineering and object-oriented technology and implemented in the C++ programming language. It has been used in applications in particle physics, nuclear physics, accelerator design, space engineering and medical physics. PACS: 07.05.Tp; 13; 2

Effects of Alpha Interferon Treatment on Intrinsic Anti-HIV-1 Immunity In Vivo

Alpha interferon (IFN-α) suppresses human immunodeficiency virus type 1 (HIV-1) replication in vitro by inducing cell-intrinsic retroviral restriction mechanisms. We investigated the effects of IFN-α/ribavirin (IFN-α/riba) treatment on 34 anti-HIV-1 restriction factors in vivo. Expression of several anti-HIV-1 restriction factors was significantly induced by IFN-α/riba in HIV/hepatitis C virus (HCV)-coinfected individuals. Fold induction of cumulative restriction factor expression in CD4+ T cells was significantly correlated with viral load reduction during IFN-α/riba treatment (r2 = 0.649; P < 0.016). Exogenous IFN-α induces supraphysiologic restriction factor expression associated with a pronounced decrease in HIV-1 viremia

Virtual histological assessment of the prenatal life history and age at death of the Upper Paleolithic fetus from Ostuni (Italy)

The fetal remains from the Ostuni 1 burial (Italy, ca 27 ka) represent a unique opportunity to explore the prenatal biological parameters, and to reconstruct the possible patho-biography, of a fetus (and its mother) in an Upper Paleolithic context. Phase-contrast synchrotron X-ray microtomography imaging of two deciduous tooth crowns and microfocus CT measurements of the right hemimandible of the Ostuni 1b fetus were performed at the SYRMEP beamline and at the TomoLab station of the Elettra - Sincrotrone laboratory (Trieste, Italy) in order to refne age at death and to report the enamel developmental history and dental tissue volumes for this fetal individual. The virtual histology allowed to estimate the age at death of the fetus at 31–33 gestational weeks. Three severe physiological stress episodes were also identifed in the prenatal enamel. These stress episodes occurred during the last two months and half of pregnancy and may relate to the death of both individuals. Compared with modern prenatal standards, Os1b’s skeletal development was advanced. This cautions against the use of modern skeletal and dental references for archaeological fnds and emphasizes the need for more studies on prenatal archaeological skeletal samples