A GATE/Geant4 beam model for the MedAustron non-isocentric proton treatment plans quality assurance

PurposeTo present a reference Monte Carlo (MC) beam model developed in GATE/Geant4 for the MedAustron fixed beam line. The proposed model includes an absolute dose calibration in Dose-Area-Product (DAP) and it has been validated within clinical tolerances for non-isocentric treatments as routinely performed at MedAustron.Material and MethodsThe proton beam model was parametrized at the nozzle entrance considering optic and energy properties of the pencil beam. The calibration in terms of absorbed dose to water was performed exploiting the relationship between number of particles and DAP by mean of a recent formalism. Typical longitudinal dose distribution parameters (range, distal penumbra and modulation) and transverse dose distribution parameters (spot sizes, field sizes and lateral penumbra) were evaluated. The model was validated in water, considering regular-shaped dose distribution as well as clinical plans delivered in non-isocentric conditions.ResultsSimulated parameters agree with measurements within the clinical requirements at different air gaps. The agreement of distal and longitudinal dose distribution parameters is mostly better than 1 mm. The dose difference in reference conditions and for 3D dose delivery in water is within 0.5% and 1.2%, respectively. Clinical plans were reproduced within 3%.ConclusionA full nozzle beam model for active scanning proton pencil beam is described using GATE/Geant4. Absolute dose calibration based on DAP formalism was implemented. The beam model is fully validated in water over a wide range of clinical scenarios and will be inserted as a reference tool for research and for independent dose calculation in the clinical routine.     

Organ doses from a proton gantry-mounted cone-beam computed tomography system characterized with MCNP6 and GATE

PurposeTo determine organ doses from a proton gantry-mounted cone-beam computed tomography (CBCT) system using two Monte Carlo codes and to study the influence on organ doses from different acquisition modes and repeated imaging.MethodsThe CBCT system was characterized with MCNP6 and GATE using measurements of depth doses in water and spatial profiles in air. The beam models were validated against absolute dose measurements and used to simulate organ doses from CBCT imaging with head, thorax and pelvis protocols. Anterior and posterior 190° scans were simulated and the resulting organ doses per mAs were compared to those from 360° scans. The influence on organ doses from repeated imaging with different imaging schedules was also investigated.ResultsThe agreement between MCNP6, GATE and measurements with regard to depth doses and beam profiles was within 4% for all protocols and the corresponding average agreement in absolute dose validation was 4%. Absorbed doses for in-field organs from 360° scans ranged between 6 and 8 mGy, 15–17 mGy and 24–54 mGy for the head, thorax and pelvis protocols, respectively. Cumulative organ doses from repeated CBCT imaging ranged between 0.04 and 0.32 Gy for weekly imaging and 0.2–1.6 Gy for daily imaging. The anterior scans resulted in an average increase in dose per mAs of 24% to the organs of interest relative to the 360° scan, while the posterior scan showed a 37% decrease.ConclusionsA proton gantry-mounted CBCT system was accurately characterized with MCNP6 and GATE. Organ doses varied greatly depending on acquisition mode, favoring posterior scans.     

Prompt-gamma emission in GEANT4 revisited and confronted with experiment

PurposeThe field of online monitoring of the beam range is one of the most researched topics in proton therapy over the last decade. The development of detectors that can be used for beam range verification under clinical conditions is a challenging task. One promising possible solution are modalities that record prompt-gamma radiation produced by the interactions of the proton beam with the target tissue. A good understanding of the energy spectra of the prompt gammas and the yields in certain energy regions is crucial for a successful design of a prompt-gamma detector. Monte-Carlo simulations are an important tool in development and testing of detector concepts, thus the proper modelling of the prompt-gamma emission in those simulations are of vital importance. In this paper, we confront a number of GEANT4 simulations of prompt-gamma emission, performed with different versions of the package and different physics lists, with experimental data obtained from a phantom irradiation with proton beams of four different energies in the range 70–230 MeV.MethodsThe comparison is made on different levels: features of the prompt-gamma energy spectrum, gamma emission depth profiles for discrete transitions and the width of the distal fall-off in those profiles.ResultsThe best agreement between the measurements and the simulations is found for the GEANT4 version 10.4.2 and the reference physics list QGSP_BIC_HP.ConclusionsModifications to prompt-gamma emission modelling in higher versions of the software increase the discrepancy between the simulation results and the experimental data.     

Investigations of line scanning proton therapy with dynamic multi-leaf collimator

PurposeScanning proton therapy has dosimetric advantage over passive treatment, but has a large penumbra in low-energy region. This study investigates the penumbra reduction when multi-leaf collimators (MLCs) are used for line scanning proton beams and secondary neutron production from MLCs.MethodsScanning beam plans with and without MLC shaping were devised. Line scanning proton plan of 36 energy layers between 71.2 and 155.2 MeV was generated. The MLCs were shaped according to the cross-sectional target shape for each energy layer. The two-dimensional doses were measured through an ion-chamber array, depending on the presence of MLC field, and Monte Carlo (MC) simulations were performed. The plan, measurement, and MC data, with and without MLC, were compared at each depth. The secondary neutron dose was simulated with MC. Ambient neutron dose equivalents were computed for the line scanning with 10 × 10 × 5 cm³ volume and maximum proton energy of 150 MeV, with and without MLCs, at lateral distances of 25–200 cm from the isocenter. The neutron dose for a wobbling plan with 10 × 10 × 5 cm³ volume was also evaluated.ResultsThe lateral penumbra width using MLC was reduced by 23.2% on average, up to a maximum of 32.2%, over the four depths evaluated. The ambient neutron dose equivalent was 18.52% of that of the wobbling beam but was 353.1% larger than the scanning open field.ConclusionsMLC field shaping with line scanning reduced the lateral penumbra and should be effective in sparing normal tissue. However, it is important to investigate the increase in neutron dose.     

Photon beam energy dependent single-arc volumetric modulated arc optimization

PurposeThe purpose of this work was to present a new single-arc mixed photon (6&18MV) VMAT (SAMP) optimization framework that concurrently optimizes for two photon energies with corresponding partial arc lengths.Methods and materialsOwing to simultaneous optimization of energy dependent intensity maps and corresponding arc locations, the proposed model poses nonlinearity. Unique relaxation constraints based on McCormick approximations were introduced for linearization. Energy dependent intensity maps were then decomposed to generate apertures. Feasibility of the proposed framework was tested on a sample of ten prostate cancer cases with lateral separation ranging from 34 cm (case no.1) to 52 cm (case no.6). The SAMP plans were compared against single energy (6MV) VMAT (SE) plans through dose volume histograms (DVHs) and radiobiological parameters including normal tissue complication probability (NTCP) and equivalent uniform dose (EUD).ResultsThe contribution of higher energy photon beam optimized by the algorithm demonstrated an increase for cases with a lateral separation >40 cm. SAMP–VMAT notably improved bladder and rectum sparing in large size cases. Compared to single energy, SAMP–VMAT plans reduced bladder and rectum NTCP in cases with large lateral separation. With the exception of one case, SAMP–VMAT either improved or maintained femoral heads compared to SE–VMAT. SAMP–VMAT reduced the nontarget tissue integral dose in all ten cases.ConclusionsA single-arc VMAT optimization framework comprising mixed photon energy partial arcs was presented. Overall results underline the feasibility and potential of the proposed approach for improving OAR sparing in large size patients without compromising the target homogeneity and coverage.     

Out-of-field doses from secondary radiation produced in proton therapy and the associated risk of radiation-induced cancer from a brain tumor treatment

PurposeTo determine out-of-field doses produced in proton pencil beam scanning (PBS) therapy using Monte Carlo simulations and to estimate the associated risk of radiation-induced second cancer from a brain tumor treatment.MethodsSimulations of out-of-field absorbed doses were performed with MCNP6 and benchmarked against measurements with tissue-equivalent proportional counters (TEPC) for three irradiation setups: two irradiations of a water phantom using proton energies of 78–147 MeV and 177–223 MeV, and one brain tumor irradiation of a whole-body phantom. Out-of-field absorbed and equivalent doses to organs in a whole-body phantom following a brain tumor treatment were subsequently simulated and used to estimate the risk of radiation-induced cancer. Additionally, the contribution of absorbed dose originating from radiation produced in the nozzle was calculated from simulations.ResultsOut-of-field absorbed doses to the TEPC ranged from 0.4 to 135 µGy/Gy. The average deviation between simulations and measurements of the water phantom irradiations was about 17%. The absorbed dose contribution from radiation produced in the nozzle ranged between 0 and 70% of the total dose; the contribution was however small in absolute terms. The absorbed and equivalent doses to the organs ranged between 0.2 and 60 µGy/Gy and 0.5–151 µSv/Gy. The estimated lifetime risk of radiation-induced second cancer was approximately 0.01%.ConclusionsThe agreement of out-of-field absorbed doses between measurements and simulations was good given the sources of uncertainties. Calculations of out-of-field organ doses following a brain tumor treatment indicated that proton PBS therapy of brain tumors is associated with a low risk of radiation-induced cancer.     

Design of commissioning process for Halcyon™ linac with a new rigid board: A clinical experience

PurposeThis study performed the accurate measurements of beam profiles with a new rigid board, which was consistent with the supplied reference beam profiles (RBPs) for clinical Halcyon model.MethodsPercentage depth doses (PDDs), lateral and diagonal dose profiles were measured and compared with RBPs. A water tank was set on the rigid board bridged Halcyon bore without sagging and source-to-surface distance was 90.0 cm. Field sizes were from 2.0 to 28.0 cm squares and depths of lateral and diagonal dose profiles were 1.3, 5.0, 10.0, and 20.0 cm. For the PDD, the depth of maximum dose (d_max), PDD value at depth of 10.0 cm (PDD₁₀), and absolute dose difference (DD) between RBP and measured beam profiles (MBP) were evaluated. For lateral and diagonal dose profiles, DDs for the whole and divided areas (central, shoulder, and extended areas) defined by third derivative, and distance-to-agreement (DTA) in the penumbra area were evaluated.ResultsFor PDDs, the differences of d_max and PDD₁₀ and DD beyond the d_max were within 1.0 mm, 0.3%, and 1.0%, respectively. For lateral and diagonal dose profiles, the DDs reached approximately 5.0% in the whole area because of penumbra area, while the DDs in the central, shoulder, and extended areas were within 1.0%, 2.0%, and 1.0%, respectively. The DTAs in the penumbra area were within 0.8 mm.ConclusionsThe supplied RBPs can be used clinically owing to the good agreement with the accurate MBPs with rigid board.     

FRoG dose computation meets Monte Carlo accuracy for proton therapy dose calculation in lung

PurposeTo benchmark and evaluate the clinical viability of novel analytical GPU-accelerated and CPU-based Monte Carlo (MC) dose-engines for spot-scanning intensity-modulated-proton-therapy (IMPT) towards the improvement of lung cancer treatment.MethodsNine patient cases were collected from the CNAO clinical experience and The Cancer Imaging Archive-4D-Lung-Database for in-silico study. All plans were optimized with 2 orthogonal beams in RayStation (RS) v.8. Forward calculations were performed with FRoG, an independent dose calculation system using a fast robust approach to the pencil beam algorithm (PBA), RS-MC (CPU for v.8) and general-purpose MC (gp-MC). Dosimetric benchmarks were acquired via irradiation of a lung-like phantom and ionization chambers for both a single-field-uniform-dose (SFUD) and IMPT plans. Dose-volume-histograms, dose-difference and γ-analyses were conducted.ResultsWith respect to reference gp-MC, the average dose to the GTV was 1.8% and 2.3% larger for FRoG and the RS-MC treatment planning system (TPS). FRoG and RS-MC showed a local γ-passing rate of ~96% and ~93%. Phantom measurements confirmed FRoG’s high accuracy with a deviation < 0.1%.ConclusionsDose calculation performance using the GPU-accelerated analytical PBA, MC-TPS and gp-MC code were well within clinical tolerances. FRoG predictions were in good agreement with both the full gp-MC and experimental data for proton beams optimized for thoracic dose calculations. GPU-accelerated dose-engines like FRoG may alleviate current issues related to deficiencies in current commercial analytical proton beam models. The novel approach to the PBA implemented in FRoG is suitable for either clinical TPS or as an auxiliary dose-engine to support clinical activity for lung patients.     

Determination of inflection points of CyberKnife dose profiles within acceptability criteria of deviations in measurements

AimThe aim of this study was to determine the Inflection Points (IPs) of flattening filter free (FFF) CyberKnife dose profiles for cone-based streotactic radiotherapy. In addition, dosimetric field sizes were determined.BackgroundThe increased need for treatment in the early stages of cancer necessitated the treatment of smaller tumors. However, efforts in that direction required the modeling accuracy of the beam. Removal of the flattening filter (FF) from the path of x-ray beam has provided the solution to those efforts, but required a different normalization approach for the beam to ensure the delivery of the dose accurately. As a solution, researchers proposed a normalization factor based on IPs.Materials and methodsMeasurements using microDiamond (PTW 60019), Diode SRS (PTW 60018) and Monte Carlo (MC) calculations of dose profiles were completed at SAD 80 cm and 5 cm depth for 15–60 mm cones. Performance analysis of detectors with respect to MC calculation was carried out. Gamma evaluation method was used to determine achievable acceptability criteria for FFF CyberKnife beams.ResultsAcceptability within (3%–0.5 mm) was found to be anachievable criterion for all dose profile measurements of the cone beams used in this study. To determine the IP, the first and second derivatives of the dose profile were determined via the cubic spline interpolation technique.ConclusionDerivatives of the interpolated profiles showed that locations of IPs and 50% isodose points coincide.     

Integration of the M6 Cyberknife in the Moderato Monte Carlo platform and prediction of beam parameters using machine learning

PurposeThis work describes the integration of the M6 Cyberknife in the Moderato Monte Carlo platform, and introduces a machine learning method to accelerate the modelling of a linac.MethodsThe MLC-equipped M6 Cyberknife was modelled and integrated in Moderato, our in-house platform offering independent verification of radiotherapy dose distributions. The model was validated by comparing TPS dose distributions with Moderato and by film measurements. Using this model, a machine learning algorithm was trained to find electron beam parameters for other M6 devices, by simulating dose curves with varying spot size and energy. The algorithm was optimized using cross-validation and tested with measurements from other institutions equipped with a M6 Cyberknife.ResultsOptimal agreement in the Monte Carlo model was reached for a monoenergetic electron beam of 6.75 MeV with Gaussian spatial distribution of 2.4 mm FWHM. Clinical plan dose distributions from Moderato agreed within 2% with the TPS, and film measurements confirmed the accuracy of the model. Cross-validation of the prediction algorithm produced mean absolute errors of 0.1 MeV and 0.3 mm for beam energy and spot size respectively. Prediction-based simulated dose curves for other centres agreed within 3% with measurements, except for one device where differences up to 6% were detected.ConclusionsThe M6 Cyberknife was integrated in Moderato and validated through dose re-calculations and film measurements. The prediction algorithm was successfully applied to obtain electron beam parameters for other M6 devices. This method would prove useful to speed up modelling of new machines in Monte Carlo systems.     

Dose calculation validation of a convolution algorithm in a solid water phantom

PurposeThe dose calculated using a convolution algorithm should be validated in a simple homogeneous water-equivalent phantom before clinical use. The dose calculation accuracy within a solid water phantom was investigated.MethodsThe specific Gamma knife design requires a dose rate calibration within a spherical solid water phantom. The TMR10 algorithm, which approximates the phantom material as liquid water, correctly computes the absolute dose in water. The convolution algorithm, which considers electron density miscalculates the dose in water as the phantom Hounsfield units were converted into higher electron density when the original CT calibration curve was used. To address this issue, the electron density of liquid water was affected by modifying the CT calibration curve. The absolute dose calculated using the convolution algorithm was compared with that computed by the TMR10. The measured depth dose profiles were also compared to those computed by the convolution and TMR10 algorithms. A patient treatment was recalculated in the solid-water phantom and the delivery quality assurance was checked.ResultsThe convolution algorithm and the TMR10 calculate an absolute dose within 1% when using the modified CT calibration curve. The dose depth profile calculated using the convolution algorithms was superimposed on the TMR10 and measured dose profiles when the modified CT calibration curve was applied. The Gamma index was better than 93%.ConclusionsDose calculation algorithms, which consider electron density, require a CT calibration curve adapted to the phantom material to correctly compute the dose in water.     

Solvent Effects are Important in Elucidating Radical-Scavenging Mechanisms of Antioxidants. A Case Study on Genistein

Abstract

In two recent researches on antioxidant mechanisms of flavonoids and isoflavonoids, quantum chemical method was employed to calculate the proton dissociation energies, the bond dissociation energies and the ionization potentials for the phenols and derived radicals to help determine the radical-scavenging mechanisms. As the solvent effect was left out of consideration, the conclusion drawn from the calculation incurred some controversies. In the current study, we re-calculated the parameters for genistein and its anion by employing a B3LYP/6–311+G(d,p)//6–31G(d,p) method with solvent effect. Accordingly, a more reasonable explanation on the experimentally observed behavior of genistein as a radical scavenger was obtained. Therefore, solvent effect should be considered in the investigation of radical-scavenging mechanisms of antioxidants in polar solvents.     

Scan path optimization with/without clustering for active beam delivery in charged particle therapy

PurposeTo compare different algorithms to optimize the scanning path in charged particle therapy with quasidiscrete scanning. We implemented a Hybrid Genetic Algorithm with Heuristics (HyGA) and combined it with clustering techniques. The performance was compared to Simulated Annealing (SA) and to commercially available treatment planning system (TPS).MethodsPerformance and clinical implications were assessed using data from 10 patients treated at CNAO (Centro Nazionale di Adroterapia Oncologica). Clinical treatments are performed relying on beam deflection, avoiding irradiation for transitions between adjacent spots larger than 2 cm. A clustering method was implemented with HyGA (HyGA_Cl), which assumes beam deflection during transition between clusters. Clinical performance was determined as the total number of particles delivered during spot transitions and the number of particles wasted due to beam deflection. Results were compared to scan paths obtained with CNAO TPS.ResultsSA and HyGA produced on average shorter paths compared to the currently available TPS. This did not result in a reduction of transit particles, due to the concomitant effect of beam deflection out of the extraction line. HyGA_Cl achieved 2% average reduction in transit particles when compared to CNAO TPS. As a drawback, wasted particles increased, due to more frequent use of beam deflection. Both the SA and HyGA algorithms reduced the number of wasted particles.ConclusionSA and HyGA proved to be the most cost-effective methods in reducing wasted particles, with benefits in terms of shorter scan paths. A decrease in transit particles delivered with beam deflection can be achieved using HyGA_Cl.     

RBE variation in prostate carcinoma cells in active scanning proton beams: In-vitro measurements in comparison with phenomenological models

PurposeIn-vitro radiobiological studies are essential for modelling the relative biological effectiveness (RBE) in proton therapy. The purpose of this study was to experimentally determine the RBE values in proton beams along the beam path for human prostate carcinoma cells (Du-145). RBE-dose and RBE-LET_d (dose-averaged linear energy transfer) dependencies were investigated and three phenomenological RBE models, i.e. McNamara, Rørvik and Wilkens were benchmarked for this cell line.MethodsCells were placed at multiple positions along the beam path, employing an in-house developed solid phantom. The experimental setup reflected the clinical prostate treatment scenario in terms of field size, depth, and required proton energies (127.2–180.1 MeV) and the physical doses from 0.5 to 6 Gy were delivered. The reference irradiation was performed with 200 kV X-ray beams. Respective (α/β) values were determined using the linear quadratic model and LET_d was derived from the treatment planning system at the exact location of cells.Results and ConclusionIndependent of the cell survival level, all experimental RBE values were consistently higher in the target than the generic clinical RBE value of 1.1; with the lowest RBE value of 1.28 obtained at the beginning of the SOBP. A systematic RBE decrease with increasing dose was observed for the investigated dose range. The RBE values from all three applied models were considerably smaller than the experimental values. A clear increase of experimental RBE values with LET_d parameter suggests that proton LET must be taken into consideration for this low (α/β) tissue.     

Determination of initial electron parameters by means of Monte Carlo simulations for the Siemens Artiste Linac 6 MV photon beam

AimIn this study, we investigated initial electron parameters of Siemens Artiste Linac with 6 MV photon beam using the Monte Carlo method.BackgroundIt is essential to define all the characteristics of initial electrons hitting the target, i.e. mean energy and full width of half maximum (FWHM) of the spatial distribution intensity, which is needed to run Monte Carlo simulations. The Monte Carlo is the most accurate method for simulation of radiotherapy treatments.Materials and methodsLinac head geometry was modeled using the BEAMnrc code. The phase space files were used as input file to DOSXYZnrc simulation to determine the dose distribution in a water phantom. We obtained percent depth dose curves and the lateral dose profile. All the results were obtained at 100 cm of SSD and for a 10 × 10 cm² field.ResultsWe concluded that there existed a good conformity between Monte Carlo simulation and measurement data when we used electron mean energy of 6.3 MeV and 0.30 cm FWHM value as initial parameters. We observed that FWHM values had very little effect on PDD and we found that the electron mean energy and FWHM values affected the lateral dose profile. However, these effects are between tolerance values.ConclusionsThe initial parameters especially depend on components of a linac head. The phase space file which was obtained from Monte Carlo Simulation for a linac can be used as calculation of scattering, MLC leakage, to compare dose distribution on patients and in various studies.     

Prompt gamma ray imaging for verification of proton boron fusion therapy: A Monte Carlo study

PurposeThe purpose of this study was to verify acquisition feasibility of a single photon emission computed tomography image using prompt gamma rays for proton boron fusion therapy (PBFT) and to confirm an enhanced therapeutic effect of PBFT by comparison with conventional proton therapy without use of boron.MethodsMonte Carlo simulation was performed to acquire reconstructed image during PBFT. We acquired percentage depth dose (PDD) of the proton beams in a water phantom, energy spectrum of the prompt gamma rays, and tomographic images, including the boron uptake region (BUR; target). The prompt gamma ray image was reconstructed using maximum likelihood expectation maximisation (MLEM) with 64 projection raw data. To verify the reconstructed image, both an image profile and contrast analysis according to the iteration number were conducted. In addition, the physical distance between two BURs in the region of interest of each BUR was measured.ResultsThe PDD of the proton beam from the water phantom including the BURs shows more efficient than that of conventional proton therapy on tumour region. A 719 keV prompt gamma ray peak was clearly observed in the prompt gamma ray energy spectrum. The prompt gamma ray image was reconstructed successfully using 64 projections. Different image profiles including two BURs were acquired from the reconstructed image according to the iteration number.ConclusionWe confirmed successful acquisition of a prompt gamma ray image during PBFT. In addition, the quantitative image analysis results showed relatively good performance for further study.     

Effects of shielding on pelvic and abdominal IORT dose distributions

PurposeTo study the impact of shielding elements in the proximity of Intra-Operative Radiation Therapy (IORT) irradiation fields, and to generate graphical and quantitative information to assist radiation oncologists in the design of optimal shielding during pelvic and abdominal IORT.MethodAn IORT system was modeled with BEAMnrc and EGS++ Monte Carlo codes. The model was validated in reference conditions by gamma index analysis against an experimental data set of different beam energies, applicator diameters, and bevel angles. The reliability of the IORT model was further tested considering shielding layers inserted in the radiation beam. Further simulations were performed introducing a bone-like layer embedded in the water phantom. The dose distributions were calculated as 3D dose maps.ResultsThe analysis of the resulting 2D dose maps parallel to the clinical axis shows that the bevel angle of the applicator and its position relative to the shielding have a major influence on the dose distribution. When insufficient shielding is used, a hotspot nearby the shield appears near the surface. At greater depths, lateral scatter limits the dose reduction attainable with shielding, although the presence of bone-like structures in the phantom reduces the impact of this effect.ConclusionsDose distributions in shielded IORT procedures are affected by distinct contributions when considering the regions near the shielding and deeper in tissue: insufficient shielding may lead to residual dose and hotspots, and the scattering effects may enlarge the beam in depth. These effects must be carefully considered when planning an IORT treatment with shielding.     

The determination of a dose deposited in reference medium due to (p,n) reaction occurring during proton therapy

AimThe aim of the investigation was to determine the undesirable dose coming from neutrons produced in reactions (p,n) in irradiated tissues represented by water.BackgroundProduction of neutrons in the system of beam collimators and in irradiated tissues is the undesirable phenomenon related to the application of protons in radiotherapy. It makes that proton beams are contaminated by neutrons and patients receive the undesirable neutron dose.Materials and methodsThe investigation was based on the Monte Carlo simulations (GEANT4 code). The calculations were performed for five energies of protons: 50 MeV, 55 MeV, 60 MeV, 65 MeV and 75 MeV. The neutron doses were calculated on the basis of the neutron fluence and neutron energy spectra derived from simulations and by means of the neutron fluence–dose conversion coefficients taken from the ICRP dosimetry protocol no. 74 for the antero-posterior irradiation geometry.ResultsThe obtained neutron doses are much less than the proton ones. They do not exceed 0.1%, 0.4%, 0.5%, 0.6% and 0.7% of the total dose at a given depth for the primary protons with energy of 50 MeV, 55 MeV, 60 MeV, 65 MeV and 70 MeV, respectively.ConclusionsThe neutron production takes place mainly along the central axis of the beam. The maximum neutron dose appears at about a half of the depth of the maximum proton dose (Bragg peak), i.e. in the volume of a healthy tissue. The doses of neutrons produced in the irradiated medium (water) are about two orders of magnitude less than the proton doses for the considered range of energy of protons.     

Quality assurance of carbon ion and proton beams: A feasibility study for using the 2D MatriXX detector

PurposeThe quality assurance (QA) procedures in particle therapy centers with active beam scanning make extensive use of films, which do not provide immediate results. The purpose of this work is to verify whether the 2D MatriXX detector by IBA Dosimetry has enough sensitivity to replace films in some of the measurements.MethodsMatriXX is a commercial detector composed of 32 × 32 parallel plate ionization chambers designed for pre-treatment dose verification in conventional radiation therapy. The detector and GAFCHROMIC® films were exposed simultaneously to a 131.44 MeV proton and a 221.45 MeV/u carbon-ion therapeutic beam at the CNAO therapy center of Pavia – Italy, and the results were analyzed and compared.ResultsThe sensitivity MatriXX on the beam position, beam width and field flatness was investigated. For the first two quantities, a method for correcting systematic uncertainties, dependent on the beam size, was developed allowing to achieve a position resolution equal to 230 μm for carbon ions and less than 100 μm for protons. The beam size and the field flatness measured using MatriXX were compared with the same quantities measured with the irradiated film, showing a good agreement.ConclusionsThe results indicate that a 2D detector such as MatriXX can be used to measure several parameters of a scanned ion beam quickly and precisely and suggest that the QA would benefit from a new protocol where the MatriXX detector is added to the existing systems.     

Optimization of treatment planning for hypoxic tumours and re-modulation of radiation intensity in heavy-ion radiotherapy

AimThe purpose of this study is to optimize treatment planning in carbon ion radiotherapy, taking into account the effect of tumour hypoxia.BackgroundIn conventional hadron therapy, the goal is to create a homogenous dose in the tumour area and, thus, achieve a uniform cell survival level. Since the induction of a specific damage to cells is directly influenced by the level of hypoxia in the tissue, the varying oxygen pressure in the different regions of hypoxic tumours would disrupt the uniformity of the cell survival level.Materials and methodsUsing the Geant4 Monte Carlo Code, the physical dose profile and dose-averaged linear energy transfer were calculated in the tumour. Then, the oxygen enhancement ratio in different areas of the tumour were compared with different pressures.ResultsModulations of radiation intensities as well as energies of ion beams were calculated, both considering and disregarding the effect of hypoxia, and the required dose profiles were compared with each other. Cell survival levels were also compared between the two methods. An equation was obtained for re-modulating the beams in the presence of hypoxia, and radiation weighting factors were extracted for the beam intensities.ConclusionThe results show that taking the effect of hypoxia into account would cause the reduction of average doses delivered to the tumour tissues up to 1.54 times. In this regard, the required dose is reduced by 1.63 times in the healthy tissues before the tumour. This will result in an effective protection of healthy tissues around the tumour.