Silicon 3D Microdetectors for Microdosimetry in Hadron Therapy

Assessing radiation dosimetry for microorganisms in naturally radioactive mineral springs using GATE and Geant4-DNA Monte Carlo simulations

Mineral springs in Massif Central, France can be characterized by higher levels of natural radioactivity in comparison to the background. The biota in these waters is constantly under radiation exposure mainly from the α-emitters of the natural decay chains, with 226Ra in sediments ranging from 21 Bq/g to 43 Bq/g and 222Rn activity concentrations in water up to 4600 Bq/L. This study couples for the first time micro- and nanodosimetric approaches to radioecology by combining GATE and Geant4-DNA to assess the dose rates and DNA damages to microorganisms living in these naturally radioactive ecosystems. It focuses on unicellular eukaryotic microalgae (diatoms) which display an exceptional abundance of teratological forms in the most radioactive mineral springs in Auvergne. Using spherical geometries for the microorganisms and based on γ-spectrometric analyses, we evaluate the impact of the external exposure to 1000 Bq/L 222Rn dissolved in the water and 30 Bq/g 226Ra in the sediments. Our results show that the external dose rates for diatoms are significant (9.7 μGy/h) and comparable to the threshold (10 μGy/h) for the protection of the ecosystems suggested by the literature. In a first attempt of simulating the radiation induced DNA damage on this species, the rate of DNA Double Strand Breaks per day is estimated to 1.11E-04. Our study confirms the significant mutational pressure from natural radioactivity to which microbial biodiversity has been exposed since Earth origin in hydrothermal springs.

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.

First experimental measurements of 2D microdosimetry maps in proton therapy

BACKGROUND

Empirical data in proton therapy indicate that relative biological effectiveness (RBE) is not constant, and it is directly related to the linear energy transfer (LET). The experimental assessment of LET with high resolution would be a powerful tool for minimizing the LET hot spots in intensity-modulated proton therapy, RBE- or LET-guided evaluation and optimization to achieve biologically optimized proton plans, verifying the theoretical predictions of variable proton RBE models, and so on. This could impact clinical outcomes by reducing toxicities in organs at risk.

PURPOSE

The present work shows the first 2D LET maps obtained at a proton therapy facility using the double scattering delivery mode in clinical conditions by means of new silicon 3D-cylindrical microdetectors.

METHODS

The device consists of a matrix of 121 independent silicon-based detectors that have 3D-cylindrical electrodes of 25-µm diameter and 20-µm depth, resulting each one of them in a well-defined micrometric radiation sensitive volume etched inside the silicon. They have been specifically designed for a hadron therapy, improving the performance of current silicon-based microdosimeters. Microdosimetry spectra were obtained at different positions of the Bragg curve by using a water-equivalent phantom along an 89-MeV pristine proton beam generated in the Y1 proton passive scattering beamline of the Orsay Proton Therapy Centre (Institut Curie, France).

RESULTS

Microdosimetry 2D-maps showing the variation of the lineal energy with depth in the three dimensions were obtained in situ during irradiation at clinical fluence rates (∼10⁸  s^-1  cm^-2 ) for the first time with a spatial resolution of 200 µm, the highest achieved in the transverse plane so far. The experimental results were cross-checked with Monte Carlo simulations and a good agreement between the spectra shapes was found. The experimental frequency-mean lineal energy values in silicon were 1.858 ± 0.019 keV µm^-1 at the entrance, 2.61 ± 0.03 keV µm^-1 at the proximal distance, 4.97 ± 0.05 keV µm^-1 close to the Bragg peak, and 8.6 ± 0.1 keV µm^-1 at the distal edge. They are in good agreement with the expected trends in the literature in clinical proton beams.

CONCLUSIONS

We present the first 2D microdosimetry maps obtained in situ during irradiation at clinical fluence rates in proton therapy. Our results show that the arrays of 3D-cylindrical microdetectors are a reliable microdosimeter to evaluate LET maps not only in the longitudinal axis of the beam, but also in the transverse plane allowing for LET characterization in three dimensions. This work is a proof of principle showing the capacity of our system to deliver LET 2D maps. This kind of experimental data is needed to validate variable proton RBE models and to optimize LET-guided plans.

3D Simulation, Electrical Characteristics and Customized Manufacturing Method for a Hemispherical Electrode Detector

The theoretical basis of a hypothetical spherical electrode detector was investigated in our previous work. It was found that the proposed detector has very good electrical characteristics, such as greatly reduced full depletion voltage, small capacitance and ultra-fast collection time. However, due to the limitations of current technology, spherical electrode detectors cannot be made. Therefore, in order to use existing CMOS technology to realize the fabrication of the detector, a hemispherical electrode detector is proposed. In this work, 3D modeling and simulation including potential and electric field distribution and hole concentration distribution are carried out using the TCAD simulation tools. In addition, the electrical characteristics, such as I-V, C-V, induced current and charge collection efficiency (CCE) with different radiation fluences, are studied to predict the radiation hardness property of the device. Furthermore, a customized manufacturing method is proposed and simulated with the TCAD-SPROCESS simulation tool. The key is to reasonably set the aspect ratio of the deep trench in the multi-step repetitive process and optimize parameters such as the angle, energy, and dose of ion implantation to realize the connection of the heavily doped region of the near-hemispherical electrode. Finally, the electrical characteristics of the process simulation are compared with the device simulation results to verify its feasibility.

Microdosimetry performance of the first multi-arrays of 3D-cylindrical microdetectors

The present work reports on the microdosimetry measurements performed with the two first multi-arrays of microdosimeters with the highest radiation sensitive surface covered so far. The sensors are based on new silicon-based radiation detectors with a novel 3D cylindrical architecture. Each system consists of arrays of independent microdetectors covering 2 mm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document}×2 mm and 0.4 mm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document}×12 cm radiation sensitive areas, the sensor distributions are arranged in layouts of 11\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document}×11 microdetectors and 3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document}×3 multi-arrays, respectively. We have performed proton irradiations at several energies to compare the microdosimetry performance of the two systems, which have different spatial resolution and detection surface. The unitcell of both arrays is a 3D cylindrical diode with a 25 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μm diameter and a 20 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μm depth that results in a welldefined and isolated radiation sensitive micro-volume etched inside a silicon wafer. Measurements were carried out at the Accélérateur Linéaire et Tandem à Orsay (ALTO) facility by irradiating the two detection systems with monoenergetic proton beams from 6 to 20 MeV at clinical-equivalent fluence rates. The microdosimetry quantities were obtained with a spatial resolution of 200 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μm and 600 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μm for the 11\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document}×11 system and for the 3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document}×3 multi-array system, respectively. Experimental results were compared with Monte Carlo simulations and an overall good agreement was found. The good performance of both microdetector arrays demonstrates that this architecture and both configurations can be used clinically as microdosimeters for measuring the lineal energy distributions and, thus, for RBE optimization of hadron therapy treatments. Likewise, the results have shown that the devices can be also employed as a multipurpose device for beam monitoring in particle accelerators.

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.

Microdosimetric characterization of a clinical proton therapy beam: comparison between simulated lineal energy distributions in spherical water targets and experimental measurements with a silicon detector

Objective Treatment planning based on computer simulations were proposed to account for the increase in the relative biological effectiveness (RBE) of proton radiotherapy beams near to the edges of the irradiated volume. Since silicon detectors could be used to validate the results of these simulations, it is important to explore the limitations of this comparison. Approach Microdosimetric measurements with a MicroPlus Bridge V2 silicon detector (thickness = 10 µm) were performed along the Bragg peak of a clinical proton beam. The lineal energy distributions, the dose mean values, and the RBE calculated with a biological weighting function were compared with simulations with PHITS (microdosimetric target = 1 µm water sphere), and published clonogenic survival in vitro RBE data for the V79 cell line. The effect of the silicon-to-water conversion was also investigated by comparing three different methodologies (conversion based on a single value, novel bin-to-bin conversions based on SRIM and PSTAR). Main results Mainly due to differences in the microdosimetric targets, the experimental dose-mean lineal energy and RBE values at the distal edge were respectively up to 53% and 28% lower than the simulated ones. Furthermore, the methodology chosen for the silicon-to-water conversion was proven to affect the dose mean lineal energy and the RBE10 up to 32% and 11% respectively. The best methodology to compensate for this underestimation was the bin-to-bin silicon-to-water conversion based on PSTAR. Significance This work represents the first comparison between PHITS-simulated lineal energy distributions in water targets and corresponding experimental spectra measured with silicon detectors. Furthermore, the effect of the silicon-to-water conversion on the RBE was explored for the first time. The proposed methodology based on the PSTAR bin-to-bin conversion appears to provide superior results with respect to commonly used single scaling factors and is recommended for future studies.

Development, Characterization and Valuable Use of Novel Dosimeter Film Based on PVA Polymer Doped Nitro Blue Tetrazolium Dye and AgNO3 for the Accurate Detection of Low X-ray Doses

Currently, the uncontrolled exposure of individuals to X-rays during medical examinations represents a substantial danger that threatens both medical professionals and patients. Therefore, radiation dosimetry for low X-ray doses is a very important control of radiation practice in medical diagnostic radiology. In line with this, the current study proposes a valuable dosimeter-based PVA thin film doubly doped with silver nitrate salt and nitro blue tetrazolium dye. The nanocomposite film was prepared via a simple casting method and the different processing parameters were optimized. The performance of radiation detection was evaluated according to optical, chromic, chemical and structural changes after exposure to variable low X-ray doses (0, 2, 4, 10 and 20 mGy). The different film labels exhibited an excellent stability behavior in dark and light upon 30 days of storage. The UV-Vis spectrophotometric study showed a gradual increase in the maximum absorbance as a function of the dose and the corresponding response curve confirmed this linear variation (R = 0.998). A clear structural modification was recorded via X-ray diffraction (XRD) analysis revealing the increase in crystallinity with the level of the dose received by the nanocomposite films. Microscopic surface analysis via SEM assessments revealed a significant morphological change in PVA/Ag⁺/NBT films exposed to increased radiation doses and typical dendrites growing in needle- or tree-like microstructures appeared with a high X-ray dose. Finally, the nanocomposite films before and after irradiation were evaluated via a spectrocolorimetric study and the different CIELab coordinates, the color difference, as well as the color strength, showed a linear correlation with the intensity of the applied dose. This new dosimeter design could, therefore, provide a promising and efficient alternative for prompt and accurate detection of low X-rays doses in diagnostic radiology.

What can space radiation protection learn from radiation oncology?

Protection from cosmic radiation of crews of long-term space missions is now becoming an urgent requirement to allow a safe colonization of the moon and Mars. Epidemiology provides little help to quantify the risk, because the astronaut group is small and as yet mostly involved in low-Earth orbit mission, whilst the usual cohorts used for radiation protection on Earth (e.g. atomic bomb survivors) were exposed to a radiation quality substantially different from the energetic charged particle field found in space. However, there are over 260,000 patients treated with accelerated protons or heavier ions for different types of cancer, and this cohort may be useful for quantifying the effects of space-like radiation in humans. Space radiation protection and particle therapy research also share the same tools and devices, such as accelerators and detectors, as well as several research topics, from nuclear fragmentation cross sections to the radiobiology of densely ionizing radiation. The transfer of the information from the cancer radiotherapy field to space is manifestly complicated, yet the two field should strengthen their relationship and exchange methods and data.

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.

Dosimetry with gafchromic films based on a new micro-opto-electro-mechanical system

This work presents the first tests performed with radiochromic films and a new Micro‒Opto‒Electro-Mechanical system (MOEMS) for in situ dosimetry evaluation in radiotherapy in real time. We present a new device and methodology that overcomes the traditional limitation of time-delay in radiochromic film analysis by turning a passive detector into an active sensor. The proposed system consists mainly of an optical sensor based on light emitting diodes and photodetectors controlled by both customized electronic circuit and graphical user interface, which enables optical measurements directly. We show the first trials performed in a low‒energy proton cyclotron with this MOEMS by using gafchromic EBT3 films. Results show the feasibility of using this system for in situ dose evaluations. Further adaptation is ongoing to develop a full real‒time active detector by integrating MOEM multi‒arrays and films in flexible printed circuits. Hence, we point to improve the clinical application of radiochromic films with the aim to optimize radiotherapy treatment verifications.