Energy imparted and ionisation yield at the nanometre scale: results for extended beams.

The metrological problem of interpreting ionisation-based micro- and nanodosimetric measurements in terms of quantities proportional to energy imparted becomes particularly relevant when the sensitive volume (SV) size is in the nanometre range. At these scales, a constant W-value cannot be assumed, and the stochastics of the energy transfer per single collision could play a more important role. This problem was recently analysed by our group by means of track-structure Monte Carlo simulations with the Geant4-DNA code, finding a strong correlation between the energy imparted and ionisation yield also for SV diameters of 1 nm. As the previous study was limited to primary beams of radius zero crossing the sensitive sphere along its diameter, it is the aim of the present work to extend the analysis to beams with a radius larger than the dimensions of the SV, to better assess the role played by secondary electrons.

A Geant4 shielding design for the first US carbon multi-ion hybrid synchrotron facility.

The Mayo Clinic Florida Integrated Oncology Building will be the home of the first spot-scanning only carbon/proton hybrid therapy system by Hitachi, Ltd. It will provide proton beams up to kinetic energies of 230 MeV and carbon beams up to 430 MeV n-1for clinical deployment. To provide adequate radiation protection, the Geant4 (v10.6) Monte Carlo toolkit was utilized to quantify the ambient dose equivalent at a 10 mm depth (H*(10)) for photons and neutrons. To perform accurate calculations of the ambient dose equivalent, three-dimensional computer-aided design files of the entire planned facility were imported into Geant4, as well as certain particle system components such as the bending magnets, fast Faraday cup, and gantry. Particle fluence was scored using 60 cm diameter spheres, which were strategically placed throughout areas of interests. Analytical calculations were performed as first-pass design checks. Major shielding slabs were optimized using Geant4 simulations iteratively, with more than 20 alternative designs evaluated within Geant4. The 430 MeV n-1carbon beams played the most significant role in concrete thickness Requirements. The primary wall thickness for the carbon fixed beam room is 4 meters. The presence of the proton gantry structure in the simulation caused the ambient dose equivalent to increase by around 67% at the maze entrance, but a decrease in the high energy beam transport corridor. All shielding primary and secondary goals for clinical operations were met per state regulation and national guidelines.

Development of a more accurate Geant4 quantum molecular dynamics model for hadron therapy.

Objective. Although in heavy-ion therapy, the quantum molecular dynamics (QMD) model is one of the most fundamental physics models providing an accurate daughter-ion production yield in the final state, there are still non-negligible differences with the experimental results. The aim of this study is to improve fragment production in water phantoms by developing a more accurate QMD model in Geant4.Approach. A QMD model was developed by implementing modern Skyrme interaction parameter sets, as well as by incorporating with an ad hocα-cluster model in the initial nuclear state. Two adjusting parameters were selected that can significantly affect the fragment productions in the QMD model: the radius to discriminate a cluster to which nucleons belong after the nucleus-nucleus reaction, denoted byR, and the squared standard deviation of the Gaussian packet, denoted byL. Squared Mahalanobis's distance of fragment yields and angular distributions with 1, 2, and the higher atomic number for the produced fragments were employed as objective functions, and multi-objective optimization (MOO), which make it possible to compare quantitatively the simulated production yields with the reference experimental data, was performed.Main results. The MOO analysis showed that the QMD model with modern Skyrme parameters coupled with the proposedα-cluster model, denoted as SkM*α, can drastically improve light fragments yields in water. In addition, the proposed model reproduced the kinetic energy distribution of the fragments accurately. The optimizedLin SkM*αwas confirmed to be realistic by the charge radii analysis in the ground state formation.Significance. The proposed framework using MOO was demonstrated to be very useful in judging the superiority of the proposed nuclear model. The optimized QMD model is expected to improve the accuracy of heavy-ion therapy dosimetry.

Estimating the biological effects of helium, carbon, oxygen, and neon ion beams using 3D silicon microdosimeters.

In this study, the survival fraction (SF) and relative biological effectiveness (RBE) of pancreatic cancer cells exposed to spread-out Bragg peak helium, carbon, oxygen, and neon ion beams are estimated from the measured microdosimetric spectra using a microdosimeter and the application of the microdosimetric kinetic (MK) model. To measure the microdosimetric spectra, a 3D mushroom silicon-on-insulator microdosimeter connected to low noise readout electronics (MicroPlus probe) was used. The parameters of the MK model were determined for pancreatic cancer cells such that the calculated SFs reproduced previously reported in vitro SF data. For a cuboid target of 10 × 10 × 6 cm3, treatment plans of helium, carbon, oxygen, and neon ion beams were designed using in-house treatment planning software (TPS) to achieve a 10% SF of pancreatic cancer cells throughout the target. The physical doses and microdosimetric spectra of the planned fields were measured at different depths in polymethyl methacrylate phantoms. The biological effects, such as SF, RBE, and RBE-weighted dose at different depths along the fields were predicted from the measurements. The predicted SFs at the target region were generally in good agreement with the planned SF from the TPS in most cases.

Response of SOI microdosimeter in fast neutron beams: experiment and Monte Carlo simulations.

In this study, Monte Carlo codes, Geant4 and MCNP6, were used to characterize the fast neutron therapeutic beam produced at iThemba LABS in South Africa. Experimental and simulation results were compared using the latest generation of Silicon on Insulator (SOI) microdosimeters from the Centre for Medical Radiation Physics (CMRP). Geant4 and MCNP6 were able to successfully model the neutron gantry and simulate the expected neutron energy spectrum produced from the reaction by protons bombarding a 9Be target. The neutron beam was simulated in a water phantom and its characteristics recorded by the silicon microdosimeters; bare and covered by a 10B enriched boron carbide converter, at different positions. The microdosimetric quantities calculated using Geant4 and MCNP6 are in agreement with experimental measurements. The thermal neutron sensitivity and production of 10B capture products in the p+ boron-implanted dopant regions of the Bridge microdosimeter is investigated. The obtained results are useful for the future development of dedicated SOI microdosimeters for Boron Neutron Capture Therapy (BNCT). This paper provides a benchmark comparison of Geant4 and MCNP6 capabilities in the context of further applications of these codes for neutron microdosimetry.

A benchmarking study of Geant4 for Auger electrons emitted by medical radioisotopes.

Auger emitting radioisotopes are of great interest in targeted radiotherapy because, once internalised in the tumour cells, they can deliver dose locally to the radiation sensitive targets, while not affecting surrounding cells. Geant4 is a Monte Carlo code widely used to characterise the physics mechanism at the basis of targeted radiotherapy. In this work, we benchmarked the modelling of the emission of Auger electrons in Geant4 deriving from the decay of 123I, 124I, 125I radionuclides against existing theoretical approaches. We also compared Geant4 against reference data in the case of 131Cs, which is of interest for brachytherapy. In the case of 125I and 131Cs, the simulation results are compared to experimental measurements as well. Good agreement was found between Geant4 and the reference data. As far as we know, this is the first study aimed to benchmark against experimental measurements the emission of Auger electrons in Geant4 for radiotherapy applications.

In-field and out-of-field microdosimetric characterisation of a 62 MeV proton beam at CATANA.

PURPOSE: A 5 and 10 µm thin silicon on insulator (SOI) 3D mushroom microdosimeter was used to characterize both the in-field and out-of-field of a 62 MeV proton beam. METHODS: The SOI mushroom microdosimeter consisted of an array of cylindrical sensitive volumes (SVs), developed by the Centre for Medical Radiation Physics, University of Wollongong, was irradiated with 62 MeV protons at the CATANA (Centro di AdroTerapia Applicazioni Nucleari Avanzate) facility in Catania, Italy, a facility dedicated to the radiation treatment of ocular melanomas. Dose mean lineal energy, ( y D ¯ ), values were obtained at various depths in PMMA along a pristine and spread out Bragg peak (SOBP). The measured microdosimetric spectra at each position were then used as inputs into the modified Microdosimetric Kinetic Model (MKM) to derive the RBE for absorbed dose in a middle of the SOBP 2Gy (RBED ). Microdosimetric spectra were obtained with both the 5 and 10 µm 3D SOI microdosimeters, with a focus on the distal part of the BP. The in-field and out-of-field measurement configurations along the Bragg curve were modeled in Geant4 for comparison with experimental results. Lateral out-of-field measurements were performed to study secondary particles' contribution to normal tissue's dose, up to 12 mm from the edge of the beam field, and quality factor and dose equivalent results were obtained. RESULTS: Comparison between experimental and simulation results showed good agreement between one another for both the pristine and SOBP beams in terms of y D ¯ and RBED. Though a small discrepancy between experiment and simulation was seen at the entrance of the Bragg curve, where experimental results were slightly lower than Geant4. The dose equivalent value measured 12 mm from the edge of the target volume was 1.27 ± 0.15 mSv/Gy with a Q ¯ value of 2.52 ± 0.30, both of which agree within uncertainty with Geant4 simulation. CONCLUSIONS: These results demonstrate that SOI microdosimeters are an effective tool to predict RBED in-field as well as dose equivalent monitoring out-of-field to provide insight to probability of second cancer generation.

Advances in modelling gold nanoparticle radiosensitization using new Geant4-DNA physics models.

Gold nanoparticles have demonstrated significant radiosensitization of cancer treatment with x-ray radiotherapy. To understand the mechanisms at the basis of nanoparticle radiosensitization, Monte Carlo simulations are used to investigate the dose enhancement, given a certain nanoparticle concentration and distribution in the biological medium. Earlier studies have ordinarily used condensed history physics models to predict nanoscale dose enhancement with nanoparticles. This study uses Geant4-DNA complemented with novel track structure physics models to accurately describe electron interactions in gold and to calculate the dose surrounding gold nanoparticle structures at nanoscale level. The computed dose in silico due to a clinical kilovoltage beam and the presence of gold nanoparticles was related to in vitro brain cancer cell survival using the local effect model. The comparison of the simulation results with radiobiological experimental measurements shows that Geant4-DNA and local effect model can be used to predict cell survival in silico in the case of x-ray kilovoltage beams.

Validation of linear energy transfer computed in a Monte Carlo dose engine of a commercial treatment planning system.

The relative biological effectiveness (RBE) of protons is highly variable and difficult to quantify. However, RBE is related to the local ionization density, which can be related to the physical measurable dose weighted linear energy transfer (LETD). The aim of this study was to validate the LETD calculations for proton therapy beams implemented in a commercially available treatment planning system (TPS) using microdosimetry measurements and independent LETD calculations (Open-MCsquare (MCS)). The TPS (RayStation v6R) was used to generate treatment plans on the CIRS-731-HN anthropomorphic phantom for three anatomical sites (brain, nasopharynx, neck) for a spherical target (Ø = 5 cm) with uniform target dose to calculate the LETD distribution. Measurements were performed at the University Medical Center Groningen proton therapy center (Proteus Plus, IBA) using a µ +-probe utilizing silicon on insulator microdosimeters capable of detecting lineal energies as low as 0.15 keV µm-1 in tissue. Dose averaged mean lineal energy [Formula: see text] depth-profiles were measured for 70 and 130 MeV spots in water and for the three treatment plans in water and an anthropomorphic phantom. The [Formula: see text] measurements were compared to the LETD calculated in the TPS and MCS independent dose calculation engine. D · [Formula: see text] was compared to D · LETD in terms of a gamma-index with a distance-to-agreement criteria of 2 mm and increasing dose difference criteria to determine the criteria for which a 90% pass rate was accomplished. Measurements of D · [Formula: see text] were in good agreement with the D · LETD calculated in the TPS and MCS. The 90% passing rate threshold was reached at different D · LETD difference criteria for single spots (TPS: 1% MCS: 1%), treatment plans in water (TPS: 3% MCS: 6%) and treatment plans in an anthropomorphic phantom (TPS: 6% MCS: 1%). We conclude that D · LETD calculations accuracy in the RayStation TPS and open MCSquare are within 6%, and sufficient for clinical D · LETD evaluation and optimization. These findings remove an important obstacle in the road towards clinical implementation of D · LETD evaluation and optimization of proton therapy treatment plans. Novelty and significance The dose weighed linear energy transfer (LETD) distribution can be calculated for proton therapy treatment plans by Monte Carlo dose engines. The relative biological effectiveness (RBE) of protons is known to vary with the LETD distribution. Therefore, there exists a need for accurate calculation of clinical LETD distributions. Previous LETD validations have focused on general purpose Monte Carlo dose engines which are typically not used clinically. We present the first validation of mean lineal energy [Formula: see text] measurements of the LETD against calculations by the Monte Carlo dose engines of the Raystation treatment planning system and open MCSquare.

Characterization of the Mixed Radiation Field Produced by Carbon and Oxygen Ion Beams of Therapeutic Energy: A Monte Carlo Simulation Study.

PURPOSE: The main advantages of charged particle radiotherapy compared to conventional X-ray external beam radiotherapy are a better tumor conformality coupled with the capability of treating deep-seated radio-resistant tumors. This work investigates the possibility to use oxygen beams for hadron therapy, as an alternative to carbon ions. MATERIALS AND METHODS: Oxygen ions have the advantage of a higher relative biological effectiveness (RBE) and better conformality to the tumor target. This work describes the mixed radiation field produced by an oxygen beam in water and compares it to the one produced by a therapeutic carbon ion beam. The study has been performed using Geant4 simulations. The dose is calculated for incident carbon ions with energies of 162 MeV/u and 290 MeV/u, and oxygen ions with energies of 192 MeV/u and 245 MeV/u, and hence that the range of the primary oxygen ions projectiles in water was located at the same depth as the carbon ions. RESULTS: The results show that the benefits of oxygen ions are more pronounced when using lower energies because of a slightly higher peak-to-entrance ratio, which allows either providing higher dose in tumor target or reducing it in the surrounding healthy tissues. It is observed that, per incident particle, oxygen ions deliver higher doses than carbon ions. CONCLUSIONS: This result coupled with the higher RBE shows that it may be possible to use a lower fluence of oxygen ions to achieve the same therapeutic dose in the patient as that obtained with carbon ion therapy.

Characterization of prompt gamma ray emission for in vivo range verification in particle therapy: A simulation study.

In this paper we investigate the emission and detection characteristics of prompt gamma (PG) rays for in vivo range verification during hadron therapy, using Geant4 simulations. Proton, 4He and 12C beams of varying energy are incident on water phantoms. The PG production yield, energy spectral characteristics and spatial correlation with the Bragg Peak (BP) have been quantified. Further, the angular distributions for PG detection with respect to a point-of-reference on the phantom surface have been explored. The temporal properties of PG emission and time-of-flight (TOF) of PG detection have also been investigated in correlation with the changing particle beam range. Our results show that the primary PG rays from nuclear interactions of the primary beam exhibit the closest correlation to the beam range but its signal is significantly masked by the concurrent secondary PG rays, particularly for heavier ions such as carbon ion beams. The PG TOF spectroscopy encodes the essential information of the beam range but requires high time resolution measurements to retrieve it. A hybrid PG detection system to utilize the energy, timing and spatial characteristics of PG rays is desirable for BP tracking in real-time.

INVESTIGATING VARIABLE RBE IN A 12C MINIBEAM FIELD WITH MICRODOSIMETRY AND GEANT4.

An experimental and simulation-based study was performed on a 12C ion minibeam radiation therapy (MBRT) field produced with a clinical broad beam and a brass multi-slit collimator (MSC). Silicon-on-insulator (SOI) microdosimeters developed at the Centre for Medical Radiation Physics (CMRP) with micron sized sensitive volumes were used to measure the microdosimetric spectra at varying positions throughout the MBRT field and the corresponding dose-mean lineal energies and RBE for 10% cell survival (RBE10) were calculated using the modified Microdosimetric Kinetic Model (MKM). An increase in the average RBE10 of ∼30% and 10% was observed in the plateau region compared to broad beam for experimental and simulation values, respectively. The experimental collimator misalignment was determined to be 0.7° by comparison between measured and simulated microdosimetric spectra at varying collimator angles. The simulated dose-mean lineal energies in the valley region between minibeams were found to be higher on average than in the minibeams due to higher LET particles being produced in these regions from the MSC. This work presents the first experimental microdosimetry measurements and characterisation of the local biological effectiveness in a MBRT field.

Monte Carlo investigation of the characteristics of radioactive beams for heavy ion therapy.

This work presents a simulation study evaluating relative biological effectiveness at 10% survival fraction (RBE10) of several different positron-emitting radionuclides in heavy ion treatment systems, and comparing these to the RBE10s of their non-radioactive counterparts. RBE10 is evaluated as a function of depth for three positron-emitting radioactive ion beams (10C, 11C and 15O) and two stable ion beams (12C and 16O) using the modified microdosimetric kinetic model (MKM) in a heterogeneous skull phantom subject to a rectangular 50 mm × 50 mm × 60 mm spread out Bragg peak. We demonstrate that the RBE10 of the positron-emitting radioactive beams is almost identical to the corresponding stable isotopes. The potential improvement in PET quality assurance image quality which is obtained when using radioactive beams is evaluated by comparing the signal to background ratios of positron annihilations at different intra- and post-irradiation time points. Finally, the incidental dose to the patient resulting from the use of radioactive beams is also quantified and shown to be negligible.

A novel methodology to assess linear energy transfer and relative biological effectiveness in proton therapy using pairs of differently doped thermoluminescent detectors.

A new methodology for assessing linear energy transfer (LET) and relative biological effectiveness (RBE) in proton therapy beams using thermoluminescent detectors is presented. The method is based on the different LET response of two different lithium fluoride thermoluminescent detectors (LiF:Mg,Ti and LiF:Mg,Cu,P) for measuring charged particles. The relative efficiency of the two detector types was predicted using the recently developed Microdosimetric d(z) Model in combination with the Monte Carlo code PHITS. Afterwards, the calculated ratio of the expected response of the two detector types was correlated with the fluence- and dose- mean values of the unrestricted proton LET. Using the obtained proton dose mean LET as input, the RBE was assessed using a phenomenological biophysical model of cell survival. The aforementioned methodology was benchmarked by exposing the detectors at different depths within the spread out Bragg peak (SOBP) of a clinical proton beam at iThemba LABS. The assessed LET values were found to be in good agreement with the results of radiation transport computer simulations performed using the Monte Carlo code GEANT4. Furthermore, the estimated RBE values were compared with the RBE values experimentally determined by performing colony survival measurements with Chinese Hamster Ovary (CHO) cells during the same experimental run. A very good agreement was found between the results of the proposed methodology and the results of the in vitro study.

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

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.

Evaluation of early radiation DNA damage in a fractal cell nucleus model using Geant4-DNA.

The advancement of multidisciplinary research fields dealing with ionising radiation induced biological damage - radiobiology, radiation physics, radiation protection and, in particular, medical physics - requires a clear mechanistic understanding of how cellular damage is induced by ionising radiation. Monte Carlo (MC) simulations provide a promising approach for the mechanistic simulation of radiation transport and radiation chemistry, towards the in silico simulation of early biological damage. We have recently developed a fully integrated MC simulation that calculates early single strand breaks (SSBs) and double strand breaks (DSBs) in a fractal chromatin based human cell nucleus model. The results of this simulation are almost equivalent to past MC simulations when considering direct/indirect strand break fraction, DSB yields and fragment distribution. The simulation results agree with experimental data on DSB yields within 13.6% on average and fragment distributions agree within an average of 34.8%.

The relative biological effectiveness for carbon, nitrogen, and oxygen ion beams using passive and scanning techniques evaluated with fully 3D silicon microdosimeters.

BACKGROUND: The aim of this study was to measure the microdosimetric distributions of a carbon pencil beam scanning (PBS) and passive scattering system as well as to evaluate the relative biological effectiveness (RBE) of different ions, namely 12 C, 14 N, and 16 O, using a silicon-on-insulator (SOI) microdosimeter with well-defined 3D-sensitive volumes (SV). Geant4 simulations were performed with the same experimental setup and results were compared to the experimental results for benchmarking. METHOD: Two different silicon microdosimeters with rectangular parallelepiped and cylindrical shaped SVs, both 10 µm in thickness were used in this study. The microdosimeters were connected to low noise electronics which allowed for the detection of lineal energies as low as 0.15 keV/µm in tissue. The silicon microdosimeters provide extremely high spatial resolution and can be used for in-field and out-of-field measurements in both passive scattering and PBS deliveries. The response of the microdosimeters was studied in 290 MeV/u 12 C, 180 MeV/u 14 N, 400 MeV/u 16 O passive ion beams, and 290 MeV/u 12 C scanning carbon therapy beam at heavy ion medical accelerator in Chiba (HIMAC) and Gunma University Heavy Ion Medical Center (GHMC), Japan, respectively. The microdosimeters were placed at various depths in a water phantom along the central axis of the ion beam, and at the distal part of the Spread Out Bragg Peak (SOBP) in 0.5 mm increments. The RBE values of the pristine Bragg peak (BP) and SOBP were derived using the microdosimetric lineal energy spectra and the modified microdosimetric kinetic model (MKM), using MKM input parameters corresponding to human salivary gland (HSG) tumor cells. Geant4 simulations were performed in order to verify the calculated depth-dose distribution from the treatment planning system (TPS) and to compare the simulated dose-mean lineal energy to the experimental results. RESULTS: For a 180 MeV/u 14 N pristine BP, the dose-mean lineal energy yD¯ obtained with two types of silicon microdosimeters started from approximately 29 keV/µm at the entrance to 92 keV/µm at the BP, with a maximum value in the range of 412 to 438 keV/µm at the distal edge. For 400 MeV/u 16 O ions, the dose-mean lineal energy yD¯ started from about 24 keV/µm at the entrance to 106 keV/µm at the BP, with a maximum value of approximately 381 keV/µm at the distal edge. The maximum derived RBE10 values for 14 N and 16 O ions were found to be 3.10 ± 0.47 and 2.93 ± 0.45, respectively. Silicon microdosimetry measurements using pencilbeam scanning 12 C ions were also compared to the passive scattering beam. CONCLUSIONS: These SOI microdosimeters with well-defined three-dimensional (3D) SVs have applicability in characterizing heavy ion radiation fields and measuring lineal energy deposition with sub-millimeter spatial resolution. It has been shown that the dose-mean lineal energy increased significantly at the distal part of the BP and SOBP due to very high LET particles. Good agreement was observed for the experimental and simulation results obtained with silicon microdosimeters in 14 N and 16 O ion beams, confirming the potential application of SOI microdosimeter with 3D SV for quality assurance in charged particle therapy.

Characterization of prompt gamma-ray emission with respect to the Bragg peak for proton beam range verification: A Monte Carlo study.

In this paper we report a Geant4 simulation study to investigate the characteristic prompt gamma (PG) emission in a water phantom for real-time monitoring of the Bragg peak (BP) during proton beam irradiation. The PG production, emission spatial correlation with the BP, and position preference for detection with respect to the BP have been quantified in different PG energy windows as a function of proton pencil-beam energy from 100 to 200MeV. The PG response to small BP shifts was evaluated using a 2cm-thick slab with different human body materials embedded in a water phantom. Our results show that the prominent characteristic PG emissions of 4.44, 5.21 and 6.13MeV exhibit distinctive correlation with the dose deposition curve. The accuracy in BP position identification using these characteristic PG rays is highly consistent as the beam energy increases from 100 to 200MeV. There exists a position preference for PG detection with respect to the BP position, which has a strong dependence on the proton beam energy and PG energies. It was also observed that a submillimeter shift of the BP position can be realized by using PG signals. These results indicate that the characteristic PG signal is sensitive and reliable for BP tracking. Although the maximization of the PG measurement associated with the BP is difficult, it can be optimized with energy and detection position preferences.

Correction factors to convert microdosimetry measurements in silicon to tissue in 12C 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 [Formula: see text] to calculate the lineal energy was introduced as an alternative to the mean chord length [Formula: see text] because it was found that adopting Cauchy's formula for the [Formula: see text] was not appropriate for the radiation field typical of HIT as it is very directional. [Formula: see text] 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 [Formula: see text]. The tissue equivalence conversion method and [Formula: see text] 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 [Formula: see text].

Characterization of proton pencil beam scanning and passive beam using a high spatial resolution solid-state microdosimeter.

PURPOSE: This work aims to characterize a proton pencil beam scanning (PBS) and passive double scattering (DS) systems as well as to measure parameters relevant to the relative biological effectiveness (RBE) of the beam using a silicon on insulator (SOI) microdosimeter with well-defined 3D sensitive volumes (SV). The dose equivalent downstream and laterally outside of a clinical PBS treatment field was assessed and compared to that of a DS beam. METHODS: A novel silicon microdosimeter with well-defined 3D SVs was used in this study. It was connected to low noise electronics, allowing for detection of lineal energies as low as 0.15 keV/µm. The microdosimeter was placed at various depths in a water phantom along the central axis of the proton beam, and at the distal part of the spread-out Bragg peak (SOBP) in 0.5 mm increments. The RBE values of the pristine Bragg peak (BP) and SOBP were derived using the measured microdosimetric lineal energy spectra as inputs to the modified microdosimetric kinetic model (MKM). Geant4 simulations were performed in order to verify the calculated depth-dose distribution from the treatment planning system (TPS) and to compare the simulated dose-mean lineal energy to the experimental results. RESULTS: For a 131 MeV PBS spot (124.6 mm R90 range in water), the measured dose-mean lineal energy yD¯ increased from 2 keV/µm at the entrance to 8 keV/µm in the BP, with a maximum value of 10 keV/µm at the distal edge. The derived RBE distribution for the PBS beam slowly increased from 0.97 ± 0.14 at the entrance to 1.04 ± 0.09 proximal to the BP, then to 1.1 ± 0.08 in the BP, and steeply rose to 1.57 ± 0.19 at the distal part of the BP. The RBE distribution for the DS SOBP beam was approximately 0.96 ± 0.16 to 1.01 ± 0.16 at shallow depths, and 1.01 ± 0.16 to 1.28 ± 0.17 within the SOBP. The RBE significantly increased from 1.29 ± 0.17 to 1.43 ± 0.18 at the distal edge of the SOBP. CONCLUSIONS: The SOI microdosimeter with its well-defined 3D SV has applicability in characterizing proton radiation fields and can measure relevant physical parameters to model the RBE with submillimeter spatial resolution. It has been shown that for a physical dose of 1.82 Gy at the BP, the derived RBE based on the MKM model increased from 1.14 to 1.6 in the BP and its distal part. Good agreement was observed between the experimental and simulation results, confirming the potential application of SOI microdosimeter with 3D SV for quality assurance in proton therapy.