Search

Microdosimetry is a particularly powerful method to estimate the relative biological effectiveness (RBE) of any mixed radiation field. This is particularly convenient for therapeutic heavy ion therapy (HIT) beams, referring to ions larger than protons, where the RBE of the beam can vary significantly along the Bragg curve. Additionally, due to the sharp dose gradients at the end of the Bragg peak (BP), or spread out BP, to make accurate measurements and estimations of the biological properties of a beam a high spatial resolution is required, less than a millimetre. This requirement makes silicon microdosimetry particularly attractive due to the thicknesses of the sensitive volumes commonly being  ∼10 [Formula: see text]m or less. Monte Carlo (MC) codes are widely used to study the complex mixed HIT radiation field as well as to model the response of novel microdosimeter detectors when irradiated with HIT beams. Therefore it is essential to validate MC codes against experimental measurements. This work compares measurements performed with a silicon microdosimeter in mono-energetic [Formula: see text], [Formula: see text] and [Formula: see text] ion beams of therapeutic energies, against simulation results calculated with the Geant4 toolkit. Experimental and simulation results were compared in terms of microdosimetric spectra (dose lineal energy, [Formula: see text]), the dose mean lineal energy, y  D and the RBE 10 , as estimated by the microdosimetric kinetic model (MKM). Overall Geant4 showed reasonable agreement with experimental measurements. Before the distal edge of the BP, simulation and experiment agreed within  ∼10% for y  D and  ∼2% for RBE 10 . Downstream of the BP less agreement was observed between simulation and experiment, particularly for the [Formula: see text] and [Formula: see text] beams. Simulation results downstream of the BP had lower values of y  D and RBE 10 compared to the experiment due to a higher contribution from lighter fragments compared to heavier fragments.

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., (© The Author(s) 2019. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Using the CMRP 'bridge' μ+ probe, microdosimetric measurements were undertaken out-of-field using a therapeutic scanning proton pencil beam and in-field using a 12C ion therapy field. These measurements were undertaken at Mayo Clinic, Rochester, USA and at HIMAC, Chiba, Japan, respectively. For a typical proton field used in the treatment of deep-seated tumors, we observed dose-equivalent values ranging from 0.62 to 0.99 mSv/Gy at locations downstream of the distal edge. Lateral measurements at depths close to the entrance and along the SOBP plateau were found to reach maximum values of 3.1 mSv/Gy and 5.3 mSv/Gy at 10 mm from the field edge, respectively, and decreased to ~0.04 mSv/Gy 120 mm from the field edge. The ability to measure the dose-equivalent with high spatial resolution is particularly relevant to healthy tissue dose calculations in hadron therapy treatments. We have also shown qualitatively and quantitively the effects critical organ motion would have in treatment using microdosimetric spectra. Large differences in spectra and RBE10 were observed for treatments where miscalculations of 12C ion range would result in critical structures being irradiated, showing the importance of motion management.

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 RBE 10 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., (© 2018 American Association of Physicists in Medicine.)

Purpose: The measurement of depth dose profiles for range and energy verification of heavy ion beams is an important aspect of quality assurance procedures for heavy ion therapy facilities. The steep dose gradients in the Bragg peak region of these profiles require the use of detectors with high spatial resolution. The aim of this work is to characterize a one dimensional monolithic silicon detector array called the "serial Dose Magnifying Glass" (sDMG) as an independent ion beam energy and range verification system used for quality assurance conducted for ion beams used in heavy ion therapy., Methods: The sDMG detector consists of two linear arrays of 128 silicon sensitive volumes each with an effective size of 2mm × 50μm × 100μm fabricated on a p-type substrate at a pitch of 200 μm along a single axis of detection. The detector was characterized for beam energy and range verification by measuring the response of the detector when irradiated with a 290 MeV/u 12 C ion broad beam incident along the single axis of the detector embedded in a PMMA phantom. The energy of the 12 C ion beam incident on the detector and the residual energy of an ion beam incident on the phantom was determined from the measured Bragg peak position in the sDMG. Ad hoc Monte Carlo simulations of the experimental setup were also performed to give further insight into the detector response., Results: The relative response profiles along the single axis measured with the sDMG detector were found to have good agreement between experiment and simulation with the position of the Bragg peak determined to fall within 0.2 mm or 1.1% of the range in the detector for the two cases. The energy of the beam incident on the detector was found to vary less than 1% between experiment and simulation. The beam energy incident on the phantom was determined to be (280.9 ± 0.8) MeV/u from the experimental and (280.9 ± 0.2) MeV/u from the simulated profiles. These values coincide with the expected energy of 281 MeV/u., Conclusions: The sDMG detector response was studied experimentally and characterized using a Monte Carlo simulation. The sDMG detector was found to accurately determine the 12 C beam energy and is suited for fast energy and range verification quality assurance. It is proposed that the sDMG is also applicable for verification of treatment planning systems that rely on particle range., (© 2017 American Association of Physicists in Medicine.)

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 12 C 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 12 C ion beam can be adopted as [Formula: see text]. The tissue equivalence conversion method and [Formula: see text] were adopted to determine the RBE 10 , calculated using a modified microdosimetric kinetic model, applied to the microdosimetric spectra resulting from the simulation study. Comparison of the RBE 10 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].

The authors attempt to establish the relative biological effectiveness (RBE) calculation for designing therapeutic proton beams on the basis of microdosimetry. The tissue-equivalent proportional counter (TEPC) was used to measure microdosimetric lineal energy spectra for proton beams at various depths in a water phantom. An RBE-weighted absorbed dose is defined as an absorbed dose multiplied by an RBE for cell death of human salivary gland (HSG) tumor cells in this study. The RBE values were calculated by a modified microdosimetric kinetic model using the biological parameters for HSG tumor cells. The calculated RBE distributions showed a gradual increase to about 1cm short of a beam range and a steep increase around the beam range for both the mono-energetic and spread-out Bragg peak (SOBP) proton beams. The calculated RBE values were partially compared with a biological experiment in which the HSG tumor cells were irradiated by the SOBP beam except around the distal end. The RBE-weighted absorbed dose distribution for the SOBP beam was derived from the measured spectra for the mono-energetic beam by a mixing calculation, and it was confirmed that it agreed well with that directly derived from the microdosimetric spectra measured in the SOBP beam. The absorbed dose distributions to planarize the RBE-weighted absorbed dose were calculated in consideration of the RBE dependence on the prescribed absorbed dose and cellular radio-sensitivity. The results show that the microdosimetric measurement for the mono-energetic proton beam is also useful for designing RBE-weighted absorbed dose distributions for range-modulated proton beams.

Background and Purpose: Characterization of the out-of-field dose profile following irradiation of the target with a 3D treatment plan delivered with modern techniques., Methods: An anthropomorphic RANDO phantom was irradiated with a treatment plan designed for a simulated 5 × 2 × 5 cm(3) tumor volume located in the center of the head. The experiment was repeated with all most common radiation treatment types (photons, protons and carbon ions) and delivery techniques (Intensity Modulated Radiation Therapy, passive modulation and spot scanning). The measurements were performed with active diamond detector and passive thermoluminescence (TLD) detectors to investigate the out-of-field dose both inside and outside the phantom., Results: The highest out-of-field dose values both on the surface and inside the phantom were measured during the treatment with 25 MV photons. In the proximity of the Planned Target Volume (PTV), the lowest lateral dose profile was observed for passively modulated protons mainly because of the presence of the collimator in combination with the chosen volume shape. In the far out-of-field region (above 100mm from the PTV), passively modulated ions were characterized by a less pronounced dose fall-off in comparison with scanned beams. Overall, the treatment with scanned carbon ions delivered the lowest dose outside the target volume., Conclusions: For the selected PTV, the use of the collimator in proton therapy drastically reduced the dose deposited by ions or photons nearby the tumor. Scanning modulation represents the optimal technique for achieving the highest dose reduction far-out-of-field., (Copyright © 2012 Elsevier Ireland Ltd. All rights reserved.)

The RBE-weighted absorbed dose, called "biological dose," has been routinely used for carbon-ion treatment planning in Japan to formulate dose prescriptions for treatment protocols. This paper presents a microdosimetric approach to measuring the biological dose, which was redefined to be derived from microdosimetric quantities measured by a tissue-equivalent proportional counter (TEPC). The TEPC was calibrated in (60)Co gamma rays to assure a traceability of the TEPC measurement to Japanese standards and to eliminate the discrepancies among matching counters. The absorbed doses measured by the TEPC were reasonably coincident with those measured by a reference ionization chamber. The RBE value was calculated from the microdosimetric spectrum on the basis of the microdosimetric kinetic model. The biological doses obtained by the TEPC were compared with those prescribed in the carbon-ion treatment planning system. We found that it was reasonable for the measured biological doses to decrease with depth around the rear SOBP region because of beam divergence, scattering effect, and fragmentation reaction. These results demonstrate that the TEPC can be an effective tool to assure the radiation quality in carbon-ion radiotherapy.

Purpose: Successful results in carbon-ion and proton radiotherapies can extend patients' lives and thus present a treatment option for younger patients; however, the undesired exposure to normal tissues outside the treatment volume is a concern. Organ-specific information on the absorbed dose and the biological effectiveness in the patient is essential for assessing the risk, but experimental dose assessment has seldom been done. In this study, absorbed doses, quality factors, and dose equivalents in water phantom outside of the irradiation field were determined based on lineal energy distributions measured with a commercial tissue equivalent proportional counter (TEPC) at passive carbon-ion and proton radiotherapy facilities., Methods: Measurements at eight positions in the water phantom were carried out at the Heavy-Ion Medical Accelerator in Chiba of the National Institute of Radiological Sciences for 400 and 290 MeV/u carbon beams and at the National Cancer Center Hospital East for a 235 MeV proton beam., Results: The dose equivalent per treatment absorbed dose at the center of the range-modulated region H/Dt, decreased as the position became farther from the beam axis and farther from the phantom surface. The values of H/Dt ranged from 6.7 to 0.16 mSv/Gy for the 400 MeV/u carbon beam, from 1.3 to 0.055 mSv/Gy for the 290 MeV/u carbon beam, and from 4.7 to 0.24 mSv/GV for the 235 MeV proton beam. The values of the dose-averaged quality factor QD ranged from 2.4 to 4.6 for the 400 MeV/u beam, from 2.8 to 5.3 for the 290 MeV/u beam, and from 5.1 to 8.2 for the proton beam. The authors also observed differences in the distributions of H/Dt and QD between the carbon and proton beams., Conclusions: The authors experimentally obtained absorbed doses, dose-averaged quality factors, and dose equivalents in water phantom outside of the irradiation field in passive carbon-ion and proton radiotherapies with TEPC. These data are very useful for estimating the risk of secondary cancer after receiving passive radiotherapies and for verifying Monte Carlo calculations.

The features of relativistic carbon-ion beams are attractive from the viewpoint of radiotherapy. They exhibit not only a superior physical dose distribution but also an increase in biological efficiency with depth, because energy loss of the beams increases as they penetrate the body. This paper reviews clinical aspects of carbon-beam radiotherapy using the experience at the National Institute of Radiological Sciences. The paper also outlines the dosimetry related to carbon-beam radiotherapy, including absolute dosimetry of the carbon beam, neutron measurements and radiation protection measurements.

Secondary neutron ambient dose equivalents per the treatment absorbed dose in passive carbon-ion and proton radiotherapies were measured using a rem meter, WENDI-II at two carbon-ion radiotherapy facilities and four proton radiotherapy facilities in Japan. Our measured results showed that (1) neutron ambient dose equivalent in carbon-ion radiotherapy is lower than that in proton radiotherapy, and (2) the difference to the measured neutron ambient dose equivalents among the facilities is within a factor of 3 depending on the operational beam setting used at the facility and the arrangement of the beam line, regardless of the method for making a laterally uniform irradiation field: the double scattering method or the single-ring wobbling method. The reoptimization of the beam line in passive particle radiotherapy is an effective way to reduce the risk of secondary cancer because installing an adjustable precollimator and designing the beam line devices with consideration of their material, thickness and location, etc., can significantly reduce the neutron exposure. It was also found that the neutron ambient dose equivalent in passive particle radiotherapy is equal to or less than that in the photon radiotherapy. This result means that not only scanning particle radiotherapy but also passive particle radiotherapy can provide reduced exposure to normal tissues around the target volume without an accompanied increase in total body dose.

In carbon therapy, doses at center of spread-out Bragg peaks depend on field size. For a small field of 5 x 5 cm2, the central dose reduces to 96% of the central dose for the open field in case of 400 MeV/n carbon beam. Assuming the broad beam injected to the water phantom is made up of many pencil beams, the transverse dose distribution can be reconstructed by summing the dose distribution of the pencil beams. We estimated dose profiles of this pencil beam through measurements of dose distributions of broad uniform beams blocked half of the irradiation fields. The dose at a distance of a few cm from the edge of the irradiation field reaches up to a few percent of the central dose. From radiation quality measurements of this penumbra, the large-angle scattered particles were found to be secondary fragments which have lower LET than primary carbon beams. Carbon ions break up in beam modifying devices or in water phantom through nuclear interaction with target nuclei. The angular distributions of these fragmented nuclei are much broader than those of primary carbon particles. The transverse dose distribution of the pencil beam can be approximated by a function of the three-Gaussian form. For a simplest case of mono-energetic beam, contributions of the Gaussian components which have large mean deviations become larger as the depth in the water phantom increases.

Treatment planning of heavy-ion radiotherapy involves predictive calculation of not only the physical dose but also the biological dose in a patient body. The biological dose is defined as the product of the physical dose and the relative biological effectiveness (RBE). In carbon-ion radiotherapy at National Institute of Radiological Sciences, the RBE value has been defined as the ratio of the 10% survival dose of 200 kVp x-rays to that of the radiation of interest for in vitro human salivary gland tumour cells. In this note, the physical and biological dose distributions of a typical therapeutic carbon-ion beam are calculated using the GEANT4 Monte Carlo simulation toolkit in comparison with those with the biological dose estimate system based on the one-dimensional beam model currently used in treatment planning. The results differed between the GEANT4 simulation and the one-dimensional beam model, indicating the physical limitations in the beam model. This study demonstrates that the Monte Carlo physics simulation technique can be applied to improve the accuracy of the biological dose distribution in treatment planning of heavy-ion radiotherapy.

The treatment planning approach used within the heavy ion tumour therapy project at GSI Darmstadt includes a biological optimisation, which is based on a biophysical model, the Local Effect Model (LEM). Here we show that the predictions of the LEM are in good agreement with clinical data obtained at the HIMAC in Chiba for the treatment of non-small-cell lung cancer, and the steep dose response for carbon ions is reproduced correctly. This steeper increase corresponds to an increasing RBE with increasing dose, which apparently is in contradiction to the systematics observed in general for in vitro measurements. A possible explanation of this discrepancy is based on the interindividual variation of photon sensitivity.

The production of projectile fragments is one of the most important, but not yet perfectly understood, problems to be considered when planning for the utilization of high-energy heavy charged particles for radiotherapy. This paper reports our investigation of the fragments' fluence and linear energy transfer (LET) spectra produced from various incident ions using an experimental approach to reveal these physical qualities of the beams. Polymethyl methacrylate, as a substitute for the human body, was used as a target. A deltaE-E counter telescope with a plastic scintillator and a BGO scintillator made it possible to identify the species of fragments based on differences of various elements. By combining a gas-flow proportional counter with a counter telescope system, LET spectra as well as fluence spectra of the fragments were derived for each element down from the primary particles to hydrogen. Among them, the information on hydrogen and helium fragments was derived for the first time. The result revealed that the number of light fragments, such as hydrogen and helium, became larger than the number of primaries in the vicinity of the range end. However, the greater part of the dose delivered to a cell was still governed by the primaries. The calculated result of a simulation used for heavy-ion radiotherapy indicated room for improving the reaction model.

A method to establish the relationship between CT number and effective density for therapeutic radiations is proposed. We approximated body tissues to mixtures of muscle, air, fat and bone. Consequently, the relationship can be calibrated only with a CT scan of their substitutes, for which we chose water, air, ethanol and potassium phosphate solution, respectively. With simple and specific corrections for non-equivalencies of the substitutes, a calibration accuracy of 1% will be achieved. We tested the calibration method with some biological materials to verify that the proposed method would offer the accuracy, simplicity and specificity required for a standard in radiotherapy treatment planning, in particular with heavy charged particles.

In heavy-ion radiotherapy, an accelerated beam is modified to realize a desired dose distribution in patients. The setup of the beam-modifying devices in the irradiation system is changed according to the patient, and it is important to check the depth dose distributions in the patient. In order to measure dose distributions realized by an irradiation system for heavy-ion radiotherapy, a multi-layer ionization chamber(MLIC) was developed. The MLIC consists of 64 ionization chambers, which are stacked mutually. The interval between each ionization chamber is about 4.1 mm water. There are signal and high voltage plates in the MILC, which are used as electrodes of the ionization chambers and phantom. Depth dose distribution from 5.09 mm to 261.92 mm water can be measured in about 30 seconds using this MLIC. Thus, it is possible to check beam quality in a short amount of time.

For exact determination of absorbed dose in heavy-ion irradiation fields which are used in radiation therapy and biological experiments, ionization chambers have been characterized with defined heavy-ion beams and correction factors. The LET (linear energy transfer) dependence of columnar recombination in a parallel-plate ionization chamber has been examined. Using 135 MeV/u carbon and neon beams, the ion collection efficiency was measured for several gases (air, carbon dioxide, argon and tissue-equivalent gas). 95 MeV/u argon beams and 90 MeV/u iron beams were also used for measurements of columnar recombination in air. As expected by Jaffe theory, the inverse of the ratio of the ionization charge to the saturated ionization charge had a linear relationship with the inverse of the electric field strength in the region below 0.002 V(-1) cm. The gradient of the line increases as the LET of the heavy ions increases. A strong LET dependence of the gradient was observed in air and carbon dioxide. The LET dependence was not observed in tissue-equivalent gas, nitrogen or argon. The exact depth-dose distribution of the heavy-ion beam was obtained by this correction of the initial recombination effect for the collected ionization charge. The columnar recombination in air was analysed using Jaffe theory; the obtained parameter b (a track radius) should be in the range between 0.001 cm and 0.005 cm, whereas the value obtained by Jaffe is 0.00179 cm. The value of the parameter b should increase as the LET of the heavy-ion beam increases in order to reproduce the experimental values of the initial recombination.