Benchmarking and validation of a Geant4–SHADOW Monte Carlo simulation for dose calculations in microbeam radiation therapy

Microbeam radiation therapy (MRT) is a synchrotron‐based radiotherapy modality that uses high‐intensity beams of spatially fractionated radiation to treat tumours. The rapid evolution of MRT towards clinical trials demands accurate treatment planning systems (TPS), as well as independent tools for the verification of TPS calculated dose distributions in order to ensure patient safety and treatment efficacy. Monte Carlo computer simulation represents the most accurate method of dose calculation in patient geometries and is best suited for the purpose of TPS verification. A Monte Carlo model of the ID17 biomedical beamline at the European Synchrotron Radiation Facility has been developed, including recent modifications, using the Geant4 Monte Carlo toolkit interfaced with the SHADOW X‐ray optics and ray‐tracing libraries. The code was benchmarked by simulating dose profiles in water‐equivalent phantoms subject to irradiation by broad‐beam (without spatial fractionation) and microbeam (with spatial fractionation) fields, and comparing against those calculated with a previous model of the beamline developed using the PENELOPE code. Validation against additional experimental dose profiles in water‐equivalent phantoms subject to broad‐beam irradiation was also performed. Good agreement between codes was observed, with the exception of out‐of‐field doses and toward the field edge for larger field sizes. Microbeam results showed good agreement between both codes and experimental results within uncertainties. Results of the experimental validation showed agreement for different beamline configurations. The asymmetry in the out‐of‐field dose profiles due to polarization effects was also investigated, yielding important information for the treatment planning process in MRT. This work represents an important step in the development of a Monte Carlo‐based independent verification tool for treatment planning in MRT.     

First trial of spatial and temporal fractionations of the delivered dose using synchrotron microbeam radiation therapy

The technical feasibility of temporal and spatial fractionations of the radiation dose has been evaluated using synchrotron microbeam radiation therapy for brain tumors in rats. A significant increase in lifespan (216%, p < 0.0001) resulted when three fractions of microbeam irradiation were applied to the tumor through three different ports, orthogonal to each other, at 24 h intervals. However, there were no long‐term survivors, and immunohistological studies revealed that 9 L tumors were not entirely ablated.     

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%.     

Synchrotron X‐ray microbeam dosimetry with a 20 micrometre resolution scintillator fibre‐optic dosimeter

Cancer is one of the leading causes of death worldwide. External beam radiation therapy is one of the most important modalities for the treatment of cancers. Synchrotron microbeam radiation therapy (MRT) is a novel pre‐clinical therapy that uses highly spatially fractionated X‐ray beams to target tumours, allowing doses much higher than conventional radiotherapies to be delivered. A dosimeter with a high spatial resolution is required to provide the appropriate quality assurance for MRT. This work presents a plastic scintillator fibre optic dosimeter with a one‐dimensional spatial resolution of 20 µm, an improvement on the dosimeter with a resolution of 50 µm that was demonstrated in previous work. The ability of this probe to resolve microbeams of width 50 µm has been demonstrated. The major limitations of this method were identified, most notably the low‐light signal resulting from the small sensitive volume, which made valley dose measurements very challenging. A titanium‐based reflective paint was used as a coating on the probe to improve the light collection, but a possible effect of the high‐Z material on the probes water‐equivalence has been identified. The effect of the reflective paint was a 28.5 ± 4.6% increase in the total light collected; it did not affect the shape of the depth‐dose profile, nor did it explain an over‐response observed when used to probe at low depths, when compared with an ionization chamber. With improvements to the data acquisition, this probe design has the potential to provide a water‐equivalent, inexpensive dosimetry tool for MRT.     

Photon dosimetry methods outside the target volume in radiation therapy: Optically stimulated luminescence (OSL), thermoluminescence (TL) and radiophotoluminescence (RPL) dosimetry

Dosimetry methods outside the target volume are still not well established in radiotherapy. Luminescence detectors due to their small dimensions, very good sensitivity, well known dose and energy response are considered as an interesting approach in verification of doses outside the treated region. The physical processes of thermoluminescence (TL), radiophotoluminescence (RPL) and optically stimulated luminescence (OSL) are very similar and can be described in terms of the energy band model of electron-hole production following irradiation.This work is a review of the main dosimetric characteristics of luminescence detectors which were used in experiments performed by EURADOS Working Group 9 for in-phantom measurements of secondary radiation (scattered and leakage photons). TL LiF:Mg,Ti detectors type MTS-7 (IFJ PAN, Poland), types TLD-100 and TLD-700 (Harshaw), OSL Al₂O₃:C detectors type nanoDot™ (Landauer Inc.) and RPL rod glass elements type GD-352M (Asahi Techno Glass Coorporation) are described. The main characteristics are discussed, together with the readout and calibration procedures which lead to a determination of absorbed dose to water.All dosimeter types used show very good uniformity, batch reproducibility and homogeneity. For improved accuracy, individual sensitivity correction factors should be applied for TL and OSL dosimeters while for RPL dosimeters there is no need for individual sensitivity corrections.The dose response of all dosimeters is linear for a wide range of doses.The energy response of GD-352M type dosimeters (with Sn filter) used for out-of-field measurements is flat for medium and low energy X-rays.The energy dependence for TLDs is low across the range of photon energies used and the energy correction was neglected. A significant over response of Al₂O₃:C OSLDs irradiated in kilovoltage photon beams was taken into account. The energy correction factor f_en was calculated by using the 2006 PENELOPE Monte Carlo code.With suitable calibration, all dosimeter types are appropriate for out-of-field dose measurements as well as for the in-phantom measurements of radiotherapy MV X-rays beams.     

Electrically active defect annealing in neutron and in ion-damaged Si

Hall effect and electrical conductivity measurements of defect annealing in 1 ohm-cm n-type and 2 ohm-cm p-type silicon were made following neutron irradiation at ～50°C. Measurements were also made following 400-keV B¹¹ ion implantation into a 100 ohm-cm n-type Si substrate. As the neutron fluence is increased the electrical effects of the damage eventually outweigh those of the chemical dopants, and further changes in the electrical properties become small. Conversely, significant electrical recovery upon annealing begins only when the electrical effects of the remaining damage become comparable to those of the chemical dopants. This condition will occur at higher anneal temperatures for higher fluence irradiations. The neutron fluence dependence of the damage and the annealing is interpreted in terms of the neutron energy per cm³. E, spent in atomic processes divided by the number/cm³, N, of electrically active dopants. When E/N ≤ 0.5 keV the electrical measurements show that the predominant defect annealing occurs below 400°C. However, when E/N > 0.5 keV electrical measurements emphasize the annealing at temperatures > 400°C. After 500°C annealing, energy levels in neutron damaged Si are observed at E_v +0.1 and E_v +0.15 eV in p-type and at E_c -0.33 eV in n-type Si. Application of the E/N criteria to room temperature implant-doped Si predicts that the electrical effects will be dominated by lattice damage even if all the implanted ions are substitutional.     

Egyptian limestone for gamma dosimetry: an electron paramagnetic resonance study

The electron paramagnetic resonance (EPR) properties of limestone from a certain Egyptian site were investigated in order to propose an efficient and low-cost gamma dosimeter. Radiation-induced free radicals were of one type which was produced in the limestone samples at g=2.0066 after exposure to gamma radiation (⁶⁰Co). EPR spectrum was recorded and analyzed. The microwave power saturation curve and the effect of changing modulation amplitude on peak-to- peak signal height were investigated. The response of limestone to different radiation doses (0.5–20 kGy) was studied. Except for the decrease in signal intensities during the first five hours following irradiation, over the period of two months fair stabilities of signal intensities were noticed. From the current results, it is possible to conclude that natural limestone may be a suitable material for radiation dosimetry in the range of irradiation processing.     

Characterization of Gafchromic EBT3 films for ultraviolet radiation dosimetry

A newly released processless Gafchromic EBT3 film was investigated for its response to ultraviolet radiation (UVR) at three different wavelengths: 365 (UVA), 302 (UVB) and 254 nm (UVC). The change in the film's optical density was dependent on the radiance exposure and UV wavelength and was found to correlate with the UV intensity. The dynamic range of the EBT3 films is ～5 to ～60 J/cm² for UVA and UVB and ～5 to ～300 J/cm² for UVC. The dose sensitivity of EBT3 films for UVA and UVB radiation was approximately 10 times higher than that for UVC radiation. The results indicate the suitability of EBT3 films for quantitative and qualitative measurements of UVR exposure.     

The electrical activity of neutron damage centres in silicon diodes

In this paper, a study is made of radiation damage in silicon p n diodes produced by low dose neutron irradiation, and a comparison is made with the results of earlier work on carbon irradiated diodes. Particular emphasis is placed on the effects of damage on the I/V characteristics of the diodes. The physical properties of the damage are studied by means of thermally stimulated current measurements, and further information is obtained from the temperature dependence of the reverse leakage.     

Aus der Praxis der Ferrosulfat-(Fricke-)Dosimetrie

Unter den chmischen Dosimetersystem genügt das Ferrosulfat-(Fricke-) Dosimeter in fast “idealer” Weise einer Reihe von Grundforderungen der praktischen Dosimetrie. Mit Hitfe einer 10^?4?10^?2 m Fe(II)- sulfatlösung in 0,4?2,0 m H₂SO₄, gesättigt mit Luft, Kann es als Vergleichsstandard zwischen verschiedenen Laboratorien bis zu einer Dosis von 5·10⁴ rad verwendet werden. Das Ferrosulfat-Dosimetersystem dient als sekundärer Stand zur Eichung anderer chemischer Dosimetersysteme.     

蒙特卡罗技术在放射诊断剂量快速计算中的应用

采用蒙特卡罗（MC）技术研究放射性核素131I的放疗特性。给出了131I衰变过程的MC抽样表达式,实现了含电子和光子的混合抽样模拟。对水客体中131I衰变的MC模拟结果表明:131I能够有效地治疗甲状腺肿瘤而几乎不影响正常组织。变密度客体的MC模拟结果表明:沉积能量与面密度之间的关系与密度无关,满足卷积核自适应的要求;并且MC模拟能够提供剂量卷积核函数,为快速剂量计算提供数据。     

蒙特卡罗技术在放射诊断剂量快速计算中的应用

采用蒙特卡罗（MC）技术研究放射性核素131I的放疗特性。给出了131I衰变过程的MC抽样表达式,实现了含电子和光子的混合抽样模拟。对水客体中131I衰变的MC模拟结果表明：131I能够有效地治疗甲状腺肿瘤而几乎不影响正常组织。变密度客体的MC模拟结果表明：沉积能量与面密度之间的关系与密度无关,满足卷积核自适应的要求;并且MC模拟能够提供剂量卷积核函数,为快速剂量计算提供数据。     

A new method of creating minibeam patterns for synchrotron radiation therapy: a feasibility study

Several synchrotrons around the world are currently developing innovative radiotherapy techniques with the aim of palliating and possibly curing human brain tumors. Amongst them, microbeam radiation therapy (MRT) and, more recently, minibeam radiation therapy (MBRT) have shown promising results. In MBRT the beam thickness ranges from 500 to 700 µm with a separation between two adjacent minibeams of the same value, whilst in MRT the thickness is of the order of 25–50 µm with a distance between adjacent microbeams of the order of 200 µm. An original method has been developed and tested at the ESRF ID17 biomedical beamline to produce the minibeam patterns. It utilizes a specially developed high‐energy white‐beam chopper whose action is synchronized with the vertical motion of the target moving at constant speed. Each opening of the chopper generates a horizontal beam print. The method described here has the advantage of being simple and reliable, and it allows for an easy control of the patient safety in future clinical trials. To study the feasibility of the method, dosimetric measurements have been performed using Gafchromic HD‐810 films and compared with Monte Carlo simulations. The results of this comparison are discussed.     

Semi empirical formula for exposure buildup factors

The nuclear data of photon buildup factor is an important concept that must be considered in nuclear safety aspects such as radiation shielding and dosimetry. The buildup factor is a coefficient that represents the contribution of collided photons with the target medium. Present work formulated a semi empirical formulae for exposure buildup factors (EBF) in the energy region 0.015–15?MeV, atomic number range 1?≤?Z?≤?92 and for mean free path up to 40?mfp. The EBFs produced by the present formula are compared with that of data available in the literature. It is found that present work agree with literature. This formula is first of its kind to calculate EBFs without using geometric progression fitting parameters. This formula may also use to calculate EBFs for compounds/mixtures/Biological samples. The present formula is useful in producing EBFs for elements and mixtures quickly. This semi empirical formula finds importance in the calculations of EBFs which intern helps in the radiation protection and dosimetry.     

Evaluation of a metaphase chromosome finder: Potential application to chromosome-based radiation dosimetry

An automated metaphase chromosome finder is described which combines a microscope, state-of-the-art computer technology and a simple decision-making algorithm. A microscope slide is systematically scanned under computer control and the location of each positive ‘signal’ placed into memory for later recall and review by a human operator. The software identifies two events, positives (the presence of a ‘signal’) and negatives (the absence of a ‘signal’). The performance of the metaphase finder was evaluated using receiver operating characteristic curve analysis. At the optimum decision threshold, the detection rates for true positives (metaphase spreads) was about 74%, false positives (type I error) about 6%, and false negatives (type II error) about 26%. The overall accuracy, which accounts for differences in the sensitivity of the detector to positive and negative events, was 89.4% (±0.01%; standard error of the mean, n = 8). Potential applications to radiation dosimetry are discussed.     

A thermoluminescence study of mineral silicates extracted from herbs in the dose range 0.5–5 Gy

The presence of silicates in many personal objects suggests their potential use at low dose as fortuitous dosimeter in an accidental radiological exposure, when conventional dosimetry is not available. The goal of the present work is the dosimetric characterization of mineral silicates extracted from the plant Hibiscus Sabdariffa L, known as Jamaica flower, in the dose range 0.5–5 Gy. By studying the radiation-induced signal in time, the temperature integration region between 210 °C and 250 °C was found to be the most stable and also reduced the effects of thermal fading in the dose reconstruction process; the dose response curve was linear between 0.5 Gy and 5 Gy. By checking the change in sensitivity after repeated exposures to ionizing radiations and to high temperature heating, no variation in the glow curve shape or peak intensities were detected. To eliminate a pre-existing background signal, all the characterization measurements were performed with aliquots “annealed” by a preliminary readout of the TL.     

EPR dosimeter material properties of potassium tartrate hemihydrate

The ability of potassium tartrate hemihydrate as a radiation sensitive material for electron paramagnetic resonance (EPR) dosimetry was investigated. The samples were subjected to different doses, in the range of 1–9 Gy of ⁶⁰Co gamma rays at room temperature. The EPR spectra were investigated through variation of signal intensity with respect to absorbed dose, magnetic field modulation amplitude, microwave power and time stability. The results indicate that the sensitivity of potassium tartrate hemihydrate is about 30% higher than that of alanine. However, the EPR signal is timely less stable within the first two weeks after irradiation.     

Possible Use of Polymethylmethacrylate Doped with Hg Dithizonate as a Dosimeter in the Mrad Range

It is the found that the 3 μ thin film polymethylmethacrylate (PMMA) doped with Hg dithizonate can be used for dosimetry in a range of adsorbed dose from about 1 … 4 Mrad from which the decrease of optical density of the main absorption band at 485 … 490 nm decreases exponentially by increasing absorbed dose. The effects of post storage irradiation and storage temperature on the radiation induced visible spectrum were investigated.     

Increase of lifespan for glioma‐bearing rats by using minibeam radiation therapy

This feasibility work assesses the therapeutic effectiveness of minibeam radiation therapy, a new synchrotron radiotherapy technique. In this new approach the irradiation is performed on 9L gliosarcoma‐bearing rats with arrays of parallel beams of width 500–700 µm. Two irradiation configurations were compared: a lateral unidirectional irradiation and two orthogonal arrays interlacing at the target. A dose escalation study was performed. A factor of three gain in the mean survival time obtained for some animals paves the way for further exploration of the different possibilities of this technique and its further optimization.     

Influence of cosmic radiation spectrum and its variation on the relative efficiency of LiF thermoluminescent detectors – Calculations and measurements

Lithium fluoride thermoluminescent detectors (TLD) were used for cosmic radiation dosimetry already in early 1960s. Since that time they have been constantly applied in numerous space missions for personal dosimetry, area monitoring, phantom measurements and dosimetry for biological experiments. The relative efficiency of TLDs, defined as the ratio of their response to a given radiation and to a reference radiation, is not constant, but depends on ionization density. This raises a question about the relative efficiency of TLDs exposed to the complex cosmic radiation spectrum encountered in Earth's orbit, which consists of a variety of particles, including heavy ions, the spectrum of which covers an extremely broad energy range. The present work is an attempt to find an answer to this question.The particle energy spectra were calculated for realistic flight conditions of the International Space Station (ISS). The calculation of the Galactic Cosmic Ray (GCR) component was based on the input spectra generated with the DLR model for solar minimum (2009) and solar maximum (2000) conditions. Contributions of trapped protons were estimated based on the AP8 model for solar minimum and maximum taking into account the altitude variations of the ISS. The interactions of the primary particles with the ISS were simulated with GEANT4 using a shielding geometry derived from the mass distribution of the Columbus Laboratory of the ISS and several constant aluminum shieldings. The calculated spectra were convoluted with the experimental data on the relative TL efficiency measured for ions ranging from H to Xe at various particle accelerators for two commonly applied TL-materials, namely LiF:Mg,Ti and LiF:Mg,Cu,P.The results showed the differences in the average TL-efficiency for these two TL-materials. For LiF:Mg,Ti the relative efficiency is within a few percent from unity for any of the analyzed values of shielding, altitude and solar cycle conditions. This means that one can assume cosmic radiation doses measured in Low Earth Orbit (LEO) with LiF:Mg,Ti detectors to be correct within such uncertainty. LiF:Mg,Cu,P underestimates the cosmic radiation doses by more than 15% in all cases. Altitude and solar cycle were found to have a very weak influence on the TL efficiency. In contrast, the influence of shielding thickness is quite significant. The reason for this is a change of contributions of radiation field components: trapped protons dominate at low shielding (97% of dose at 1 g/cm²), but are negligible above 60 g/cm², as well as changes within GCR spectrum (increase of dose due to lower LET secondaries for higher shielding). Shielding thickness affects both TLD types in different ways: the efficiency of LiF:Mg,Cu,P increases with increasing shielding thickness, while the efficiency of LiF:Mg,Ti shows some fluctuations, with a weak minimum for 60 g/cm². The response ratio of these TLDs decreases monotonically with the shielding thickness and could be used as an indicator for the average shielding conditions in which the TLDs were exposed.