Physics study of microbeam radiation therapy with PSI-version of Monte Carlo code GEANT as a new computational tool

A high resolution dose calculation engine for x-ray microbeams radiation therapy

BACKGROUND

Microbeam radiation therapy (MRT) is a treatment modality based on spatial fractionation of synchrotron generated x-rays into parallel, high dose, microbeams of a few microns width. MRT is still an under-development radiosurgery technique for which, promising preclinical results on brain tumors and epilepsy encourages its clinical transfer.

PURPOSE

A safe clinical transfer of MRT needs a specific treatment planning system (TPS) that provides accurate dose calculations in human patients, taking into account the MRT beams properties (high dose gradients, spatial fractionation, polarization effects). So far, the most advanced MRT treatment planning system, based on a hybrid dose calculation algorithm, is limited to a macroscopic rendering of the dose and does not account for the complex dose distribution inherent to MRT if delivered as conformal irradiations with multiple incidences. For overcoming these limitations, a multi-scale full Monte-Carlo calculation engine called penMRT has been developed and benchmarked against two general purpose Monte Carlo codes: penmain based on PENELOPE and Gate based on Geant4.

METHODS

PenMRT, is based on the PENELOPE (2018) Monte Carlo (MC) code, modified to take into account the voxelized geometry of the patients (CT-scans) and offering an adaptive micrometric dose calculation grid independent to the CT size, location and orientation. The implementation of the dynamic memory allocation in penMRT, makes the simulations feasible within a huge number of dose scoring bins. The possibility of using a source replication approach to simulate arrays of microbeams, and the parallelization using OpenMPI have been added to penMRT in order to increase the calculation speed for clinical usages. This engine can be implemented in a TPS as a dose calculation core.

RESULTS

The performance tests highlight the reliability of penMRT to be used for complex irradiation conditions in MRT. The benchmarking against a standard PENELOPE code did not show any significant difference for calculations in centimetric beams, for a single microbeam and for a microbeam array. The comparisons between penMRT and Gate as an independent MC code did not show any difference in the beam paths, whereas in valley regions, relative differences between the two codes rank from 1 to 7.5% which are probably due to the differences in physics lists that are used in these two codes. The reliability of the source replication approach has also been tested and validated with an underestimation of no more than 0.6% in low dose areas.

CONCLUSIONS

Good agreements (a relative difference between 0 to 8%) were found when comparing calculated peak to valley dose ratio (PVDR) values using penMRT, for irradiations with a full microbeam array, with calculated values in the literature. The high-resolution calculated dose maps obtained with penMRT are used to extract differential and cumulative dose-volume histograms (DVHs) and analyze treatment plans with much finer metrics regarding the irradiation complexity. To our knowledge, these are the first high-resolution dose maps and associated DVHs ever obtained for cross-fired microbeams irradiation, which is bringing a significant added value to the field of treatment planning in spatially fractionated radiation therapy. This article is protected by copyright. All rights reserved.

Investigating the Dosimetric Characteristics of Microbeam Radiation Treatment

BACKGROUND

High-radiation therapeutic gain could be achieved by the modern technique of microbeam radiation treatment (MRT). The aim of this study was to investigate the dosimetric properties of MRT.

METHODS

The EGSnrc Monte Carlo (MC) code system was used to transport photons and electrons in MRT. The mono-energetic beams (1 cm × 1 cm array) of 50, 100, and 150 keV and the spectrum photon beam (European Synchrotron Radiation Facility [ESRF]) were modeled to transport through multislit collimators with the aperture's widths of 25 and 50 μm and the center-to-center (c-t-c) distance between two adjacent microbeams (MBs) of 200 μm. The calculated phase spaces at the upper surface of water phantom (1 cm × 1 cm) were implemented in DOSXYZnrc code to calculate the percentage depth dose (PDD), the dose profile curves (in depths of 0-1, 1-2, and 3-4 cm), and the peak-to-valley dose ratios (PVDRs).

RESULTS

The PDD, dose profile curves, and PVDRs were calculated for different effective parameters. The more flatness of lateral dose profile was obtained for the ESRF spectrum MB. With constant c-t-c distance, an increase in the MB size increased the peak and valley dose; simultaneously, the PVDR was larger for the 25 μm MB (33.5) compared to 50 μm MB (21.9) beam, due to the decreased scattering photons followed to the lower overlapping of the adjacent MBs. An increase in the depth decreased the PVDRs (i.e., 54.9 in depth of 0-1 cm).

CONCLUSION

Our MC model of MRT successfully calculated the effect of dosimetric parameters including photon's energy, beam width, and depth to estimate the dose distribution.

Technical advances in x-ray microbeam radiation therapy

In the last 25 years microbeam radiation therapy (MRT) has emerged as a promising alternative to conventional radiation therapy at large, third generation synchrotrons. In MRT, a multi-slit collimator modulates a kilovoltage x-ray beam on a micrometer scale, creating peak dose areas with unconventionally high doses of several hundred Grays separated by low dose valley regions, where the dose remains well below the tissue tolerance level. Pre-clinical evidence demonstrates that such beam geometries lead to substantially reduced damage to normal tissue at equal tumour control rates and hence drastically increase the therapeutic window. Although the mechanisms behind MRT are still to be elucidated, previous studies indicate that immune response, tumour microenvironment, and the microvasculature may play a crucial role. Beyond tumour therapy, MRT has also been suggested as a microsurgical tool in neurological disorders and as a primer for drug delivery. The physical properties of MRT demand innovative medical physics and engineering solutions for safe treatment delivery. This article reviews technical developments in MRT and discusses existing solutions for dosimetric validation, reliable treatment planning and safety. Instrumentation at synchrotron facilities, including beam production, collimators and patient positioning systems, is also discussed. Specific solutions reviewed in this article include: dosimetry techniques that can cope with high spatial resolution, low photon energies and extremely high dose rates of up to 15 000 Gy s-1, dose calculation algorithms-apart from pure Monte Carlo Simulations-to overcome the challenge of small voxel sizes and a wide dynamic dose-range, and the use of dose-enhancing nanoparticles to combat the limited penetrability of a kilovoltage energy spectrum. Finally, concepts for alternative compact microbeam sources are presented, such as inverse Compton scattering set-ups and carbon nanotube x-ray tubes, that may facilitate the transfer of MRT into a hospital-based clinical environment. Intensive research in recent years has resulted in practical solutions to most of the technical challenges in MRT. Treatment planning, dosimetry and patient safety systems at synchrotrons have matured to a point that first veterinary and clinical studies in MRT are within reach. Should these studies confirm the promising results of pre-clinical studies, the authors are confident that MRT will become an effective new radiotherapy option for certain patients.

Validation of a Monte Carlo simulation for Microbeam Radiation Therapy on the Imaging and Medical Beamline at the Australian Synchrotron

Microbeam Radiation Therapy (MRT) is an emerging cancer treatment modality characterised by the use of high-intensity synchrotron-generated x-rays, spatially fractionated by a multi-slit collimator (MSC), to ablate target tumours. The implementation of an accurate treatment planning system, coupled with simulation tools that allow for independent verification of calculated dose distributions are required to ensure optimal treatment outcomes via reliable dose delivery. In this article we present data from the first Geant4 Monte Carlo radiation transport model of the Imaging and Medical Beamline at the Australian Synchrotron. We have developed the model for use as an independent verification tool for experiments in one of three MRT delivery rooms and therefore compare simulation results with equivalent experimental data. The normalised x-ray spectra produced by the Geant4 model and a previously validated analytical model, SPEC, showed very good agreement using wiggler magnetic field strengths of 2 and 3 T. However, the validity of absolute photon flux at the plane of the Phase Space File (PSF) for a fixed number of simulated electrons was unable to be established. This work shows a possible limitation of the G4SynchrotronRadiation process to model synchrotron radiation when using a variable magnetic field. To account for this limitation, experimentally derived normalisation factors for each wiggler field strength determined under reference conditions were implemented. Experimentally measured broadbeam and microbeam dose distributions within a Gammex RMI457 Solid Water® phantom were compared to simulated distributions generated by the Geant4 model. Simulated and measured broadbeam dose distributions agreed within 3% for all investigated configurations and measured depths. Agreement between the simulated and measured microbeam dose distributions agreed within 5% for all investigated configurations and measured depths.

Hybrid dose calculation: a dose calculation algorithm for microbeam radiation therapy

Microbeam radiation therapy (MRT) is still a preclinical approach in radiation oncology that uses planar micrometre wide beamlets with extremely high peak doses, separated by a few hundred micrometre wide low dose regions. Abundant preclinical evidence demonstrates that MRT spares normal tissue more effectively than conventional radiation therapy, at equivalent tumour control. In order to launch first clinical trials, accurate and efficient dose calculation methods are an inevitable prerequisite. In this work a hybrid dose calculation approach is presented that is based on a combination of Monte Carlo and kernel based dose calculation. In various examples the performance of the algorithm is compared to purely Monte Carlo and purely kernel based dose calculations. The accuracy of the developed algorithm is comparable to conventional pure Monte Carlo calculations. In particular for inhomogeneous materials the hybrid dose calculation algorithm out-performs purely convolution based dose calculation approaches. It is demonstrated that the hybrid algorithm can efficiently calculate even complicated pencil beam and cross firing beam geometries. The required calculation times are substantially lower than for pure Monte Carlo calculations.

Experimental optimisation of the X-ray energy in microbeam radiation therapy

Microbeam radiation therapy has demonstrated superior normal tissue sparing properties compared to broadbeam radiation fields. The ratio of the microbeam peak dose to the valley dose (PVDR), which is dependent on the X-ray energy/spectrum and geometry, should be maximised for an optimal therapeutic ratio. Simulation studies in the literature report the optimal energy for MRT based on the PVDR. However, most of these studies have considered different microbeam geometries to that at the Imaging and Medical Beamline (50 μm beam width with a spacing of 400 μm). We present the first fully experimental investigation of the energy dependence of PVDR and microbeam penumbra. Using monochromatic X-ray energies in the range 40-120 keV the PVDR was shown to increase with increasing energy up to 100 keV before plateauing. PVDRs measured for pink beams were consistently higher than those for monochromatic energies similar or equivalent to the average energy of the spectrum. The highest PVDR was found for a pink beam average energy of 124 keV. Conversely, the microbeam penumbra decreased with increasing energy before plateauing for energies above 90 keV. The effect of bone on the PVDR was investigated at energies 60, 95 and 120 keV. At depths greater than 20 mm beyond the bone/water interface there was almost no effect on the PVDR. In conclusion, the optimal energy range for MRT at IMBL is 90-120 keV, however when considering the IMBL flux at different energies, a spectrum with 95 keV weighted average energy was found to be the best compromise.

Sparing of tissue by using micro-slit-beam radiation therapy reduces neurotoxicity compared with broad-beam radiation therapy

Micro-slit-beam radiation therapy (MRT) using synchrotron-generated X-ray beams allows for extremely high-dose irradiation. However, the toxicity of MRT in central nervous system (CNS) use is still unknown. To gather baseline toxicological data, we evaluated mortality in normal mice following CNS-targeted MRT. Male C57BL/6 J mice were head-fixed in a stereotaxic frame. Synchrotron X-ray-beam radiation was provided by the SPring-8 BL28B2 beam-line. For MRT, radiation was delivered to groups of mice in a 10 × 12 mm unidirectional array consisting of 25-μm-wide beams spaced 100, 200 or 300 μm apart; another group of mice received the equivalent broad-beam radiation therapy (BRT) for comparison. Peak and valley dose rates of the MRT were 120 and 0.7 Gy/s, respectively. Delivered doses were 96-960 Gy for MRT, and 24-120 Gy for BRT. Mortality was monitored for 90 days post-irradiation. Brain tissue was stained using hematoxylin and eosin to evaluate neural structure. Demyelination was evaluated by Klüver-Barrera staining. The LD₅₀ and LD₁₀₀ when using MRT were 600 Gy and 720 Gy, respectively, and when using BRT they were 80 Gy and 96 Gy, respectively. In MRT, mortality decreased as the center-to-center beam spacing increased from 100 μm to 300 μm. Cortical architecture was well preserved in MRT, whereas BRT induced various degrees of cerebral hemorrhage and demyelination. MRT was able to deliver extremely high doses of radiation, while still minimizing neuronal death. The valley doses, influenced by beam spacing and irradiated dose, could represent important survival factors for MRT.

Assessment and Implications of Scattered Microbeam and Broadbeam Synchrotron Radiation for Bystander Effect Studies

Synchrotron radiation is an excellent tool for investigating bystander effects in cell and animal models because of the well-defined and controllable configuration of the beam. Although synchrotron radiation has many advantages for such studies compared to conventional radiation, the contribution of dose exposure from scattered radiation nevertheless remains a source of concern. Therefore, the influence of scattered radiation on the detection of bystander effects induced by synchrotron radiation in biological in vitro models was evaluated. Radiochromic XRQA2 film-based dosimetry was employed to measure the absorbed dose of scattered radiation in cultured cells at various distances from a field exposed to microbeam radiotherapy and broadbeam X-ray radiation. The level of scattered radiation was dependent on the distance, dose in the target zone and beam mode. The number of γ-H2AX foci in cells positioned at the same target distances was measured and used as a biodosimeter to evaluate the absorbed dose. A correlation of absorbed dose values measured by the physical and biological methods was identified. The γ-H2AX assay successfully quantitated the scattered radiation in the range starting from 10 mGy and its contribution to the observed radiation-induced bystander effect.

Medical physics aspects of the synchrotron radiation therapies: Microbeam radiation therapy (MRT) and synchrotron stereotactic radiotherapy (SSRT)

Stereotactic Synchrotron Radiotherapy (SSRT) and Microbeam Radiation Therapy (MRT) are both novel approaches to treat brain tumor and potentially other tumors using synchrotron radiation. Although the techniques differ by their principles, SSRT and MRT share certain common aspects with the possibility of combining their advantages in the future. For MRT, the technique uses highly collimated, quasi-parallel arrays of X-ray microbeams between 50 and 600 keV. Important features of highly brilliant Synchrotron sources are a very small beam divergence and an extremely high dose rate. The minimal beam divergence allows the insertion of so called Multi Slit Collimators (MSC) to produce spatially fractionated beams of typically ∼25-75 micron-wide microplanar beams separated by wider (100-400 microns center-to-center(ctc)) spaces with a very sharp penumbra. Peak entrance doses of several hundreds of Gy are extremely well tolerated by normal tissues and at the same time provide a higher therapeutic index for various tumor models in rodents. The hypothesis of a selective radio-vulnerability of the tumor vasculature versus normal blood vessels by MRT was recently more solidified. SSRT (Synchrotron Stereotactic Radiotherapy) is based on a local drug uptake of high-Z elements in tumors followed by stereotactic irradiation with 80 keV photons to enhance the dose deposition only within the tumor. With SSRT already in its clinical trial stage at the ESRF, most medical physics problems are already solved and the implemented solutions are briefly described, while the medical physics aspects in MRT will be discussed in more detail in this paper.

Pilot study for compact microbeam radiation therapy using a carbon nanotube field emission micro-CT scanner

PURPOSE

Microbeam radiation therapy (MRT) is defined as the use of parallel, microplanar x-ray beams with an energy spectrum between 50 and 300 keV for cancer treatment and brain radiosurgery. Up until now, the possibilities of MRT have mainly been studied using synchrotron sources due to their high flux (100s Gy/s) and approximately parallel x-ray paths. The authors have proposed a compact x-ray based MRT system capable of delivering MRT dose distributions at a high dose rate. This system would employ carbon nanotube (CNT) field emission technology to create an x-ray source array that surrounds the target of irradiation. Using such a geometry, multiple collimators would shape the irradiation from this array into multiple microbeams that would then overlap or interlace in the target region. This pilot study demonstrates the feasibility of attaining a high dose rate and parallel microbeam beams using such a system.

METHODS

The microbeam dose distribution was generated by our CNT micro-CT scanner (100 μm focal spot) and a custom-made microbeam collimator. An alignment assembly was fabricated and attached to the scanner in order to collimate and superimpose beams coming from different gantry positions. The MRT dose distribution was measured using two orthogonal radiochromic films embedded inside a cylindrical phantom. This target was irradiated with microbeams incident from 44 different gantry angles to simulate an array of x-ray sources as in the proposed compact CNT-based MRT system. Finally, phantom translation in a direction perpendicular to the microplanar beams was used to simulate the use of multiple parallel microbeams.

RESULTS

Microbeams delivered from 44 gantry angles were superimposed to form a single microbeam dose distribution in the phantom with a FWHM of 300 μm (calculated value was 290 μm). Also, during the multiple beam simulation, a peak to valley dose ratio of ~10 was found when the phantom translation distance was roughly 4x the beam width. The first prototype CNT-based x-ray tube dedicated to the development of compact MRT technology development was proposed and planned based on the preliminary experimental results presented here and the previous corresponding Monte Carlo simulations.

CONCLUSIONS

The authors have demonstrated the feasibility of creating microbeam dose distributions at a high dose rate using a proposed compact MRT system. The flexibility of CNT field emission x-ray sources could possibly bring compact and low cost MRT devices to the larger research community and assist in the translational research of this promising new approach to radiation therapy.

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.

Monte Carlo-based dose calculation engine for minibeam radiation therapy

Minibeam radiation therapy (MBRT) is an innovative radiotherapy approach based on the well-established tissue sparing effect of arrays of quasi-parallel micrometre-sized beams. In order to guide the preclinical trials in progress at the European Synchrotron Radiation Facility (ESRF), a Monte Carlo-based dose calculation engine has been developed and successfully benchmarked with experimental data in anthropomorphic phantoms. Additionally, a realistic example of treatment plan is presented. Despite the micron scale of the voxels used to tally dose distributions in MBRT, the combination of several efficiency optimisation methods allowed to achieve acceptable computation times for clinical settings (approximately 2 h). The calculation engine can be easily adapted with little or no programming effort to other synchrotron sources or for dose calculations in presence of contrast agents.

Monte Carlo simulation of a compact microbeam radiotherapy system based on carbon nanotube field emission technology

PURPOSE

Microbeam radiation therapy (MRT) is an experimental radiotherapy technique that has shown potent antitumor effects with minimal damage to normal tissue in animal studies. This unique form of radiation is currently only produced in a few large synchrotron accelerator research facilities in the world. To promote widespread translational research on this promising treatment technology we have proposed and are in the initial development stages of a compact MRT system that is based on carbon nanotube field emission x-ray technology. We report on a Monte Carlo based feasibility study of the compact MRT system design.

METHODS

Monte Carlo calculations were performed using EGSnrc-based codes. The proposed small animal research MRT device design includes carbon nanotube cathodes shaped to match the corresponding MRT collimator apertures, a common reflection anode with filter, and a MRT collimator. Each collimator aperture is sized to deliver a beam width ranging from 30 to 200 μm at 18.6 cm source-to-axis distance. Design parameters studied with Monte Carlo include electron energy, cathode design, anode angle, filtration, and collimator design. Calculations were performed for single and multibeam configurations.

RESULTS

Increasing the energy from 100 kVp to 160 kVp increased the photon fluence through the collimator by a factor of 1.7. Both energies produced a largely uniform fluence along the long dimension of the microbeam, with 5% decreases in intensity near the edges. The isocentric dose rate for 160 kVp was calculated to be 700 Gy∕min∕A in the center of a 3 cm diameter target. Scatter contributions resulting from collimator size were found to produce only small (<7%) changes in the dose rate for field widths greater than 50 μm. Dose vs depth was weakly dependent on filtration material. The peak-to-valley ratio varied from 10 to 100 as the separation between adjacent microbeams varies from 150 to 1000 μm.

CONCLUSIONS

Monte Carlo simulations demonstrate that the proposed compact MRT system design is capable of delivering a sufficient dose rate and peak-to-valley ratio for small animal MRT studies.

Monte carlo simulation of an X-ray pixel beam microirradiation system

Monte Carlo simulations are used in the development of a nanotechnology-based multi-pixel beam array small animal microirradiation system. The microirradiation system uses carbon nanotube field emission technology to generate arrays of individually controllable X-ray pixel beams that electronically form irregular irradiation fields having intensity and temporal modulation without any mechanical motion. The microirradiation system, once developed, will be incorporated with the micro-CT system already developed that is based on the same nanotechnology to form an integrated image-guided and intensity-modulated microirradiation system for high-temporal-resolution small animal research. Prospective microirradiation designs were evaluated based on dosimetry calculated using EGSnrc-based Monte Carlo simulations. Design aspects studied included X-ray anode design, collimator design, and dosimetric considerations such as beam energy, dose rate, inhomogeneity correction, and the microirradiation treatment planning strategies. The dosimetric properties of beam energies between 80-400 kVp with varying filtration were studied, producing a pixel beam dose rate per current of 0.35-13 Gy per min per mA at the microirradiation isocenter. Using opposing multi-pixel-beam array pairs reduces the dose inhomogeneity between adjacent pixel beams to negligible levels near the isocenter and 20% near the mouse surface.

Dosimetric studies of micropencil X-ray beam interacting with labelled tissues by Au and Gd agents using Geant4

The aim of microbeam radiation therapy is to deliver a high dose to tumours while sparing adjacent healthy tissues. Recovery of normal tissues injured by the beam irradiation and ablation of tumour are dependent on the dose distribution generated by the incident microbeams. Using microbeams has the advantage that the areas outside the beams' trajectories (valley region) are poorly irradiated by the radiation scattered inside the tissues. Thus, the normal tissues not directly irradiated are adequately preserved, resulting in a rapid regeneration of blood vessels in the directly irradiated areas (peak region). The goal of this work was to study the effects of using gold (Au) and gadolinium (Gd) as dose enhancement factors on the radial dose distribution when target tissue is irradiated by a micropencil X-ray beam. The Monte Carlo Geant4 simulation program was used to evaluate dose distribution in the phantom in two phases. In phase 1, validity of this model based on Geant4 was evaluated by comparing the obtained results with those of the published reports. In phase 2 of this simulation, Au and Gd were introduced to the assumed cancerous cylindrical shell-shaped region both on the surface (i.e. in the 0-1 cm depth of phantom) and in the depth (i.e. in the 4-5 cm depth of phantom). Then the phantom was exposed to a micropencil beam mimicking the typical conditions used at the European synchrotron radiation facility in the simulated model. The simulated dose profiles indicate that introducing high Z elements considerably enhances the absorbed dose both in the beam path and in the surrounding region. However, this enhancement is more effective for Au in the beam path and for Gd in the surrounding region. This approach of introducing high Z elements leading to their accumulation in cancerous tissue could hopefully prepare new treatment planning of preclinical trials.

Calculation of dose profiles in stereotactic Synchrotron microplanar beam radiotherapy in a tissue-lung phantom

Synchrotron x-ray beams with high fluence rate and highly collimated may be used in stereotactic radiotherapy of lung tumours. A bundle of converging monochromatic x-ray beams having uniform microscopic thickness i.e. (microplanar beams) are directed to the center of the tumour, delivering lethal dose to the target volume while sparing normal cells. The proposed technique takes advantage of the hypothesised repair mechanism of capillary cells between alternate microbeam zones, which regenerate the lethally irradiated endothelial cells. The sharply dropping lateral dose of a microbeam provides low scattered dose to the off-target interbeam volume. In the target volume the converging bundle of beams are closely spaced, and relatively high primary and secondary electron doses overlap and produce a high dose region between the beams. This higher and lower dose margins in the target volume allows precise targeting. The advantages of stereotactic microbeam radiotherapy will be lost as the dose between microbeams exceeds the tolerance dose of the dose limiting tissues. Therefore, it is essential to optimize the interbeam doses in off-target volume. The lateral and depth doses of 100 keV microplanar beams are investigated for a single beam and an array of converging microplanar beams in a tissue, lung and tissue-lung phantoms. The EGS5 Monte Carlo code is used to calculate dose profiles at different depths and bundles of beams. The maximum dose on the beam axis (peak) and the minimum interbeam dose (valley) are compared at different energies, depths, bundle sizes, heights, widths and beam spacings. The interbeam dose is calculated at different depths and an isodose map of the phantom is obtained. An acceptable energy region is found for tissue and lung microbeam radiotherapy and a stereotactic microbeam radiotherapy model is proposed for a 4 cm diameter and 1 cm thick tumour on the lung phantom.

Determination of dosimetrical quantities used in microbeam radiation therapy (MRT) with Monte Carlo simulations

Microbeam radiation therapy (MRT) is being performed by using an array of narrow rectangular x-ray beams (typical beam sizes 25 microm X 1 cm), positioned close to each other (typically 200 microm separation), to irradiate a target tissue. The ratio of peak-to-valley doses (PVDR's) in the composite dose distribution has been found to be strongly correlated with the normal tissue tolerance and the therapeutic effect of MRT. In this work a Monte Carlo (MC) study of the depth- and lateral-dose profiles in water for single x-ray microbeams of different shapes and energies has been performed with the MC code PENELOPE. The contributions to the dose deposition from different interaction types have been determined at different distances from the center of the microbeam. The dependence of the peak dose, in a water phantom, on the microbeam field size used in the preclinical trials, has been demonstrated. Composite dose distributions for an array of microbeams were obtained using superposition algorithms and PVDR's were determined and compared with literature results obtained with other Monte Carlo codes. The dependence of the PVDR's on microbeam width, x-ray energy used, and on the separation between adjacent microbeams has been studied in detail.

Interlaced x-ray microplanar beams: a radiosurgery approach with clinical potential

Studies have shown that x-rays delivered as arrays of parallel microplanar beams (microbeams), 25- to 90-microm thick and spaced 100-300 microm on-center, respectively, spare normal tissues including the central nervous system (CNS) and preferentially damage tumors. However, such thin microbeams can only be produced by synchrotron sources and have other practical limitations to clinical implementation. To approach this problem, we first studied CNS tolerance to much thicker beams. Three of four rats whose spinal cords were exposed transaxially to four 400-Gy, 0.68-mm microbeams, spaced 4 mm, and all four rats irradiated to their brains with large, 170-Gy arrays of such beams spaced 1.36 mm, all observed for 7 months, showed no paralysis or behavioral changes. We then used an interlacing geometry in which two such arrays at a 90-degree angle produced the equivalent of a contiguous beam in the target volume only. By using this approach, we produced 90-, 120-, and 150-Gy 3.4 x 3.4 x 3.4 mm(3) exposures in the rat brain. MRIs performed 6 months later revealed focal damage within the target volume at the 120- and 150-Gy doses but no apparent damage elsewhere at 120 Gy. Monte Carlo calculations indicated a 30-microm dose falloff (80-20%) at the edge of the target, which is much less than the 2- to 5-mm value for conventional radiotherapy and radiosurgery. These findings strongly suggest potential application of interlaced microbeams to treat tumors or to ablate nontumorous abnormalities with minimal damage to surrounding normal tissue.

In vivo two-photon microscopy study of short-term effects of microbeam irradiation on normal mouse brain microvasculature

PURPOSE

The purpose of this study was to assess the early effects of microbeam irradiation on the vascular permeability and volume in the parietal cortex of normal nude mice using two-photon microscopy and immunohistochemistry.

METHODS AND MATERIALS

The upper part of the left hemisphere of 55 mice was irradiated anteroposteriorly using 18 vertically oriented beams (width 25 microm, interdistance 211 microm; peak entrance doses: 312 or 1000 Gy). At different times after microbeam exposure, the microvasculature in the cortex was analyzed using intravital two-photon microscopy after intravascular injection of fluorescein isothiocyanate (FITC)-dextrans and sulforhodamine B (SRB). Changes of the vascular volume were observed at the FITC wavelength over a maximum depth of 650 mum from the dura. The vascular permeability was detected as extravasations of SRB.

RESULTS

For all times (12 h to 1 month) after microbeam irradiation and for both doses, the FITC-dextran remained in the vessels. No significant change in vascular volume was observed between 12 h and 3 months after irradiation. Diffusion of SRB was observed in microbeam irradiated regions from 12 h until 12 days only after a 1000 Gy exposure.

CONCLUSION

No radiation damage to the microvasculature was detected in normal brain tissue after a 312 Gy microbeam irradiation. This dose would be more appropriate than 1000 Gy for the treatment of brain tumors using crossfired microbeams.

Dose distribution from x‐ray microbeam arrays applied to radiation therapy: An EGS4 Monte Carlo study

We present EGS4 Monte Carlo calculations of the spatial distribution of the dose deposited by a single x-ray pencil beam, a planar microbeam, and an array of parallel planar microbeams as used in radiation therapy research. The profiles of the absorbed dose distribution in a phantom, including the peak-to-valley ratio of the dose distribution from microbeam arrays, were calculated at micrometer resolution. We determined the dependence of the findings on the main parameters of photon and electron transport. The results illustrate the dependence of the electron range and the deposited in-beam dose on the cut-off energy, of the electron transport, as well as the effects on the dose profiles of the beam energy, the array size, and the beam spacing. The effect of beam polarization also was studied for a single pencil beam and for an array of parallel planar microbeams. The results show that although the polarization effect on the dose distribution from a 3 cm x 3 cm microbeam array inside a water phantom is large enough to be measured at the outer side of the array (16% difference of the deposited dose for x-ray beams of 200 keV), it is not detectable at the array's center, thus being irrelevant for the radiation therapy purposes. Finally we show that to properly compare the dose profiles determined with a metal oxide semiconductor field emission transistor detector with the computational method predictions, it is important to simulate adequately the size and the material of the device's Si active element.

New irradiation geometry for microbeam radiation therapy

Microbeam radiation therapy (MRT) has the potential to treat infantile brain tumours when other kinds of radiotherapy would be excessively toxic to the developing normal brain. MRT uses extraordinarily high doses of x-rays but provides unusual resistance to radioneurotoxicity, presumably from the migration of endothelial cells from 'valleys' into 'peaks', i.e., into directly irradiated microslices of tissues. We present a novel irradiation geometry which results in a tolerable valley dose for the normal tissue and a decreased peak-to-valley dose ratio (PVDR) in the tumour area by applying an innovative cross-firing technique. We propose an MRT technique to orthogonally crossfire two arrays of parallel, nonintersecting, mutually interspersed microbeams that produces tumouricidal doses with small PVDRs where the arrays meet and tolerable radiation doses to normal tissues between the microbeams proximal and distal to the tumour in the paths of the arrays.

Photon penetration and scatter in micro-pinhole imaging: a Monte Carlo investigation

Pinhole SPECT is rapidly gaining popularity for imaging laboratory animals using gamma-emitting molecules. Penetration and scattering of gamma radiation in the pinhole edge material can account for a significant fraction of the total number of photons detected, particularly if the pinholes have small diameters. This study characterizes the effects of penetration and scatter with micro-pinholes made of lead, tungsten, gold and platinum. Monte Carlo simulations are performed for 1-125 (27-35 keV) and Tc-99m (140 keV) point sources with pinhole diameters ranging from 50 to 500 microm. The simulations account for the effects of photo-electric interaction, Rayleigh scattering, Compton scattering, ionization, bremsstrahlung and electron multiple scattering. As a typical example, in the case of a Tc-99m point source and pinholes with a diameter of 300 microm in gold or platinum, approximately 55% of the photons detected resulted from penetration and approximately 3% from scatter. For pinhole diameters ranging from 100 to 500 microm, the penetration fraction for tungsten and lead was approx a factor of 1.0 to 1.6 higher and the scatter fraction was 1.0 to 1.8 times higher than in case of gold or platinum. Using I-125 instead of Tc-99m decreases the penetration fraction by a factor ranging from 3 to 11 and the scatter fraction by a factor ranging from 12 to 40. For all materials studied, the total amounts of penetrated and scattered photons changed approximately linearly with respect to the pinhole diameter.

MOSFET dosimetry for microbeam radiation therapy at the European Synchrotron Radiation Facility

Preclinical experiments are carried out with approximately 20-30 microm wide, approximately 10 mm high parallel microbeams of hard, broad-"white"-spectrum x rays (approximately 50-600 keV) to investigate microbeam radiation therapy (MRT) of brain tumors in infants for whom other kinds of radiotherapy are inadequate and/or unsafe. Novel physical microdosimetry (implemented with MOSFET chips in the "edge-on" mode) and Monte Carlo computer-simulated dosimetry are described here for selected points in the peak and valley regions of a microbeam-irradiated tissue-equivalent phantom. Such microbeam irradiation causes minimal damage to normal tissues, possible because of rapid repair of their microscopic lesions. Radiation damage from an array of parallel microbeams tends to correlate with the range of peak-valley dose ratios (PVDR). This paper summarizes comparisons of our dosimetric MOSFET measurements with Monte Carlo calculations. Peak doses at depths <22 mm are 18% less than Monte Carlo values, whereas those depths >22 mm and valley doses at all depths investigated (2 mm-62 mm) are within 2-13% of the Monte Carlo values. These results lend credence to the use of MOSFET detector systems in edge-on mode for microplanar irradiation dosimetry.