Measurement of carbon ion microdosimetric distributions with ultrathin 3D silicon diodes

State-of-the-art and potential of experimental microdosimetry in ion-beam therapy

In radiotherapy, radiation-quality should be an expression of the biological and physical characteristics of ionizing radiation such as spatial distribution of ionization or energy deposition. Linear energy transfer (LET) and lineal energy (y) are two descriptors used to quantify the radiation quality. These two quantities are connected and exhibit similar features. In ion-beam therapy (IBT), lineal energy can be measured with microdosimeters, which are specifically designed to cope with the high fluence of particles in clinical beams, while the quantification of LET is generally based on calculations. In pre-clinical studies, microdosimetric spectra are used for the indirect determination of relative biological effectiveness (RBE), e.g., using the microdosimetric kinetic model (MKM) or biophysical response functions. In this context it is important to consider saturation effects, which occur when the highest values of y become less biologically relevant compared to the relative contribution they make to the physical dose. Recent clinical data suggests that local tumor control and normal tissue effects can be linked to macroscopic and microscopic dosimetry parameters. In particular, positive clinical outcomes have been correlated to the highest LET values in the density distribution, and there is no evident link to the saturation discussed above. A systematic collection of microdosimetric information in combination with clinical data in retrospective studies may clarify the role of radiation quality at the highest LET. In the clinical setting, microdosimetry is not widely used yet, despite its potential to be linked with LET by experimentally-determined y values. Through this connection, both play an important role in complex therapy techniques such as intensity modulated particle therapy (IMPT), LET-painting and multi-ion optimization. This review summarizes the current state of microdosimetry for IBT and its potential, as well as research and development needed to make experimental microdosimetry a mature procedure in a clinical context.

Evaluation of a therapeutic carbon beam using a G2000 glass scintillator

A G2000 glass scintillator (G2000-SC) was used to determine the carbon profile and range of a 290-MeV/n carbon beam used in heavy-ion therapy because it was sensitive enough to detect single-ion hits at hundreds of mega electron Volts. An electron-multiplying charge-coupled device camera was used to detect the ion luminescence generated during the irradiation of G2000-SC with the beam. The resulting image showed that the position of the Bragg peak can be determined. The beam passes through the 112-mm-thick water phantom and stops 5.73 ± 0.03 mm from the incident side to the G2000-SC. Additionally, the location of the Bragg peak was simulated when irradiating G2000-SC with the beam using the Monte Carlo code particle and heavy ion transport system (PHITS). Simulation results show that the incident beam stops at 5.60 mm after entering G2000-SC. The beam stop location obtained from images and the PHITS code is defined at 80% distal fall-off from the Bragg peak position. Consequently, G2000-SC provided effective profile measurements of therapeutic carbon beams.

Predicting the Biological Effects of Human Salivary Gland Tumour Cells for Scanned 4He-, 12C-, 16O-, and 20Ne-Ion Beams Using an SOI Microdosimeter

Experimental microdosimetry along with the microdosimetric kinetic (MK) model can be utilized to predict the biological effects of ions. To predict the relative biological effectiveness (RBE) of ions and the survival fraction (SF) of human salivary gland tumour (HSGc-C5) cells, microdosimetric quantities measured by a silicon-on-insulator (SOI) MicroPlus-mushroom microdosimeter along the spread-out Bragg peak (SOBP) delivered by pencil beam scanning of 4He, 12C, 16O, and 20Ne ions were used. The MK model parameters of HSGc-C5 cells were obtained from the best fit of the calculated SF for the different linear energy transfer (LET) of these ions and the formerly reported in vitro SF for the same LET and ions used for calculations. For a cube-shaped target of 10 × 10 × 6 cm3, treatment plans for 4He, 12C, 16O, and 20Ne ions were produced with proprietary treatment planning software (TPS) aiming for 10% SF of HSGc-C5 cells over the target volume and were delivered to a polymethyl methacrylate (PMMA) phantom. Afterwards, the saturation-corrected dose-mean lineal energy derived based on the measured microdosimetry spectra, along with the physical dose at various depths in PMMA phantoms, was used for the estimation of the SF, RBE, and RBE-weighted dose using the MK model. The predicted SF, RBE, and the RBE-weighted dose agreed with what was planned by the TPS within 3% at most depths for these ions.

Microdosimetry in low energy proton beam at therapeutic-equivalent fluence rate with silicon 3D-cylindrical microdetectors

In this work we show the first microdosimetry measurements on a low energy proton beam with therapeutic-equivalent fluence rates by using the second generation of 3D-cylindrical microdetectors. The sensors belong to an improved version of a novel silicon-based 3D-microdetector design with electrodes etched inside silicon, which were manufactured at the National Microelectronics Centre (IMB-CNM, CSIC) in Spain. A new microtechnology has been employed using quasi-toroid electrodes of 25μm diameter and a depth of 20μm within the silicon bulk, resulting in a well-defined cylindrical radiation sensitive volume. These detectors were tested at the 18 MeV proton beamline of the cyclotron at the National Accelerator Centre (CNA, Spain). They were assembled into an in-house low-noise readout electronics to assess their performance at a therapeutic-equivalent fluence rate. Microdosimetry spectra of lineal energy were recorded at several proton energies starting from 18 MeV by adding 50μm thick tungsten foils gradually at the exit-window of the cyclotron external beamline, which corresponds to different depths along the Bragg curve. The experimentalyF¯values in silicon cover from (5.7 ± 0.9) to (8.5 ± 0.4) keV μm^-1in the entrance to (27.4 ± 2.3) keV μm^-1in the distal edge. Pulse height energy spectra were crosschecked with Monte Carlo simulations and an excellent agreement was obtained. This work demonstrates the capability of the second generation 3D-microdetectors to assess accurate microdosimetric distributions at fluence rates as high as those used in clinical centers in proton therapy.

Silicon 3D Microdetectors for Microdosimetry in Hadron Therapy

The present overview describes the evolution of new microdosimeters developed in the National Microelectronics Center in Spain (IMB-CNM, CSIC), ranging from the first ultra-thin 3D diodes (U3DTHINs) to the advanced 3D-cylindrical microdetectors, which have been developed over the last 10 years. In this work, we summarize the design, main manufacture processes, and electrical characterization of these devices. These sensors were specifically customized for use in particle therapy and overcame some of the technological challenges in this domain, namely the low noise capability, well-defined sensitive volume, high spatial resolution, and pile-up robustness. Likewise, both architectures reduce the loss of charge carriers due to trapping effects, the charge collection time, and the voltage required for full depletion compared to planar silicon detectors. In particular, a 3D‒cylindrical architecture with electrodes inserted into the silicon bulk and with a very well‒delimited sensitive volume (SV) mimicked a cell array with shapes and sizes similar to those of mammalian cells for the first time. Experimental tests of the carbon beamlines at the Grand Accélérateur National d’Lourds (GANIL, France) and Centro Nazionale Adroterapia Oncologica (CNAO, Italy) showed the feasibility of the U3DTHINs in hadron therapy beams and the good performance of the 3D‒cylindrical microdetectors for assessing linear energy distributions of clinical beams, with clinical fluence rates of 5 × 10⁷ s⁻¹cm⁻² without saturation. The dose-averaged lineal energies showed a generally good agreement with Monte Carlo simulations. The results indicated that these devices can be used to characterize the microdosimetric properties in hadron therapy, even though the charge collection efficiency (CCE) and electronic noise may pose limitations on their performance, which is studied and discussed herein. In the last 3D‒cylindrical microdetector generation, we considerably improved the CCE due to the microfabrication enhancements, which have led to shallower and steeper dopant profiles. We also summarize the successive microdosimetric characterizations performed with both devices in proton and carbon beamlines.

Experimental validation of an analytical microdosimetric model based on Geant4-DNA simulations by using a silicon-based microdosimeter

PURPOSE

To study the agreement between proton microdosimetric distributions measured with a silicon-based cylindrical microdosimeter and a previously published analytical microdosimetric model based on Geant4-DNA in-water Monte Carlo simulations for low energy proton beams.

METHODS AND MATERIAL

Distributions for lineal energy (y) are measured for four proton monoenergetic beams with nominal energies from 2.0 MeV to 4.5 MeV, with a tissue equivalent proportional counter (TEPC) and a silicon-based microdosimeter. The actual energy for protons traversing the silicon-based microdosimeter is simulated with SRIM. Monoenergetic beams with these energies are simulated with Geant4-DNA code by simulating a water cylinder site of dimensions equal to those of the microdosimeter. The microdosimeter response is calibrated by using the distribution peaks obtained from the TEPC. Analytical calculations fory ¯ F andy ¯ D using our methodology based on spherical sites are also performed choosing the equivalent sphere to be checked against experimental results.

RESULTS

Distributions for y at silicon are converted into tissue equivalent and compared to the Geant4-DNA simulated, yielding maximum deviations of 1.03% fory ¯ F and 1.17% fory ¯ D . Our analytical method generates maximum deviations of 1.29% and 3.33%, respectively, with respect to experimental results.

CONCLUSION

Simulations in Geant4-DNA with ideal cylindrical sites in liquid water produce similar results to the measurements in an actual silicon-based cylindrical microdosimeter properly calibrated. The found agreement suggests the possibility to experimentally verify the calculated clinicaly ¯ D with our analytical method.

Impact of charge collection efficiency and electronic noise on the performance of solid-state 3D microdetectors

Microdosimetry has been traditionally performed through gaseous proportional counters, although in recent years different solid-state microdosimeters have been proposed and constructed for this task. In this paper, we analyze the response of solid-state devices of micrometric size with no intrinsic gain developed by CNM-CSIC (Spain). There are two major aspects of the operation of these devices that affect the reconstruction of the probability distributions and momenta of stochastic quantities related to microdosimetry. For micrometric volumes, the drift and diffusion of the charge carriers gives rise to a partial charge collection efficiency in the peripheral region of the depleted volume. This effect produces a perturbation of the reconstructed pulse height (i.e. imparted energy) distributions with respect to the actual microdosimetric distributions. The relevance of this deviation depends on the size, geometry and operating conditions of the device. On the other hand, the electronic noise from the single-event readout set-up poses a limit on the minimum detectable lineal energy when the microdosimeter size is reduced. This article addresses these issues to provide a framework on the physical constraints for the design and operation of solid-state microdosimeters.

Study on the RBE estimation for carbon beam scanning irradiation using a solid-state microdosimeter

PURPOSE

The purpose of this study was to study the field size effect on the estimated Relative Biological Effectiveness (RBE) for carbon scanning beam irradiation.

METHODS

A silicon-on-insulator (SOI) microdosimeter system developed by the Centre for Medical Radiation Physics, University of Wollongong, Australia, was used for lineal-energy measurements (microdosimetric quantity). The RBE values were derived based on the modified microdosimetric kinetic model (MKM) at different depths in a water phantom in the scanning carbon beam for various scanned areas.

RESULTS

Our study shows that the difference in RBE values derived from the SOI microdosimeter measurements with the MKM model and from the Treatment Planning System (TPS). The difference of the RBE values is within 6.5 % at the peak point of the spread-out Bragg Peak (SOBP) region. Compared to the spot-beam, RBE values obtained in the scanned-beam with a larger scanned area of 1.0 × 1.0 cm² have better agreement with which estimated by the TPS.

CONCLUSIONS

This study shows the possibility of using the SOI microdosimeter system as a quality assurance (QA) tool for RBE evaluation in carbon-pencil beam scanning radiotherapy.

Report of a National Cancer Institute special panel: Characterization of the physical parameters of particle beams for biological research

PURPOSE

To define the physical parameters needed to characterize a particle beam in order to allow intercomparison of different experiments performed using different ions at the same facility and using the same ion at different facilities.

METHODS

At the request of the National Cancer Institute (NCI), a special panel was convened to review the current status of the field and to provide suggested metrics for reporting the physical parameters of particle beams to be used for biological research. A set of physical parameters and measurements that should be performed by facilities and understood and reported by researchers supported by NCI to perform pre-clinical radiobiology and medical physics of heavy ions were generated.

RESULTS

Standard measures such as radiation delivery technique, beam modifiers used, nominal energy, field size, physical dose and dose rate should all be reported. However, more advanced physical measurements, including detailed characterization of beam quality by microdosimetric spectrum and fragmentation spectra, should also be established and reported. Details regarding how such data should be incorporated into Monte Carlo simulations and the proper reporting of simulation details are also discussed.

CONCLUSIONS

In order to allow for a clear relation of physical parameters to biological effects, facilities and researchers should establish and report detailed physical characteristics of the irradiation beams utilized including both standard and advanced measures. Biological researchers are encouraged to actively engage facility staff and physicists in the design and conduct of experiments. Modeling individual experimental setups will allow for the reporting of the uncertainties in the measurement or calculation of physical parameters which should be routinely reported.

Linear energy transfer (LET) spectra and survival fraction distribution based on the CR-39 plastic charged-particle detector in a spread-out Bragg peak irradiation by a 12C beam

Facilities for heavy ion therapies are steadily increasing in number worldwide. One of the advantages of heavy ions is their high relative biological effect (RBE). In a model used at NIRS (National Institute of Radiological Sciences), linear energy transfer (LET) spectra are required to estimate biological dose (physical dose  ×  RBE). The CR-39 plastic charged-particle detector (CR-39) is suitable for measurement of LET. For the present study, done at the Gunma University Heavy Ion Medical Center (GHMC), we measured LET spectra at 11 depths in spread-out Bragg peak (SOBP) irradiation by a ¹²C beam of 380 MeV/u. The lower threshold of the CR-39 to measure LET was about 5 keV µm^-1 due to poor sensitivity for low LET. Then we calculated biological dose and survival fraction distributions and compared them with treatment planning results at GHMC. We used Monte Carlo simulation (Geant4) to calculate LET spectra. The simulation results were in good agreement with the experimental spectra. Moreover, the biological dose and survival fraction distributions estimated from the CR-39 reproduced the treatment planning. The CR-39 is suitable for estimating biological dose in carbon ion therapy.