Simulation and experimental benchmarking of a proton pencil beam scanning nozzle model for development of MR-integrated proton therapy

BACKGROUND

MR-integrated proton therapy is under development. It consists of the unique challenge of integrating a proton pencil beam scanning (PBS) beam line nozzle with an magnetic resonance imaging (MRI) scanner. The magnetic interaction between these two components is deemed high risk as the MR images can be degraded if there is cross-talk during beam delivery and image acquisition.

PURPOSE

To create and benchmark a self-consistent proton PBS nozzle model for empowering the next stages of MR-integrated proton therapy development, namely exploring and de-risking complete integrated prototype system designs including magnetic shielding of the PBS nozzle.

MATERIALS AND METHODS

Magnetic field (COMSOL Multiphysics ${\text{Multiphysics}}$ ) and radiation transport (Geant4) models of a proton PBS nozzle located at OncoRay (Dresden, Germany) were developed according to the manufacturers specifications. Geant4 simulations of the PBS process were performed by using magnetic field data generated by the COMSOL Multiphysics ${\text{Multiphysics}}$ simulations. In total 315 spots were simulated which consisted of a40 × 30 cm 2 $40\times 30\,{\text{cm}}^{2}$ scan pattern with 5 cm spot spacings and for proton energies of 70, 100, 150, 200, and 220 MeV. Analysis of the simulated deflection at the beam isocenter plane was performed to determine the self-consistency of the model. The magnetic fringe field from a sub selection of 24 of the 315 spot simulations were directly compared with high precision magnetometer measurements. These focused on the maximum scanning setting of ± $\pm$  20 cm beam deflection as generated from the second scanning magnet in the PBS for a proton beam energy of 220 MeV. Locations along the beam line central axis (CAX) were measured at beam isocenter and downstream of 22, 47, 72, 97, and 122 cm. Horizontal off-axis positions were measured at 22 cm downstream of isocenter ( ± $\pm$  50, ± $\pm$  100, and ± $\pm$  150 cm from CAX).

RESULTS

The proton PBS simulations had good spatial agreement to the theoretical values in all 315 spots examined at the beam line isocenter plane (0-2.9 mm differences or within 1.5 % of the local spot deflection amount). Careful analysis of the experimental measurements were able to isolate the changes in magnetic fields due solely to the scanning magnet contribution, and showed 1.9  ± $\pm$  1.2μ T $\bf{\mu} {\text{T}}$ -9.4 ± $\pm$  1.2μ T $\bf{\mu} {\text{T}}$ changes over the range of measurement locations. Direct comparison with the equivalent simulations matched within the measurement apparatus and setup uncertainty in all but one measurement point.

CONCLUSIONS

For the first time a robust, accurate and self-consistent model of a proton PBS nozzle assembly has been created and successfully benchmarked for the purposes of advancing MR-integrated proton therapy research. The model will enable confidence in further simulation based work on fully integrated designs including MRI scanners and PBS nozzle magnetic shielding in order to de-risk and realize the full potential of MR-integrated proton therapy.

Development of an algorithm for proton dose calculation in magnetic fields

BACKGROUND

The advantages of proton therapy can be further enhanced with online magnetic resonance imaging (MRI) guidance. One of the challenges in the realization of MRI-guided proton therapy (MRPT) is accurately calculating the radiation dose in the presence of magnetic fields.

PURPOSE

This study aims to develop an efficient and accurate proton dose calculation algorithm adapted to the presence of magnetic fields.

METHODS

An analytical-numerical radiation dose calculation algorithm, Proton and Ion Dose Engine (PRIDE), was developed. The algorithm combines the pencil beam algorithm (PBA) with a novel iterative voxel-based ray-tracing algorithm. The new ray-tracing method uses fewer assumptions and ensures broader applicability for proton beam trajectory prediction in magnetic fields, and has been compared to Wolf's method and Schellhammer's method. The accuracy of PRIDE algorithm was validated on three phantoms and two practical plans (one single-field water plan and one prostate tumor plan) in different magnetic field strengths up to 3.0 T. The validation was performed by comparing the results against the Monte Carlo (MC) simulations, using the global gamma index criteria of 2%/2 mm and 3%/3 mm with a 10% threshold.

RESULTS

PRIDE showed good agreement with MC in homogeneous and slab heterogeneous phantom, achieving gamma passing rates (%GPs) above 99% for 2%/2 mm criteria when magnetic field strength is not greater than 1.5 T. Although the agreement decreased for scenarios involving high proton energy (240 MeV) and strong magnetic field (3.0 T), the 2%/2 mm %GPs still remained above 98%. In lateral heterogeneous phantom, the accuracy of PRIDE decreased due to the PBA's limitation. For the two practical plans in different magnetic fields, %GPs exceeded 98% and 99% for 2%/2 mm and 3%/3 mm criteria, respectively.

CONCLUSIONS

PRIDE can perform efficient and accurate proton dose calculation in magnetic fields up to 3.0 T, and is expected to work as a useful tool for proton dose calculation in MRPT.

The impact of motion on onboard MRI-guided pencil beam scanned proton therapy treatments

Objective.Online magnetic resonance imaging (MRI) guidance could be especially beneficial for pencil beam scanned (PBS) proton therapy of tumours affected by respiratory motion. For the first time to our knowledge, we investigate the dosimetric impact of respiratory motion on MRI-guided proton therapy compared to the scenario without magnetic field.Approach.A previously developed analytical proton dose calculation algorithm accounting for perpendicular magnetic fields was extended to enable 4D dose calculations. For two geometrical phantoms and three liver and two lung patient cases, static treatment plans were optimised with and without magnetic field (0, 0.5 and 1.5 T). Furthermore, plans were optimised using gantry angle corrections (0.5 T +5° and 1.5 T +15°) to reproduce similar beam trajectories compared to the 0 T reference plans. The effect of motion was then considered using 4D dose calculations without any motion mitigation and simulating 8-times volumetric rescanning, with motion for the patient cases provided by 4DCT(MRI) data sets. Each 4D dose calculation was performed for different starting phases and the CTV dose coverageV95%and homogeneityD5%-D95%were analysed.Main results.For the geometrical phantoms with rigid motion perpendicular to the beam and parallel to the magnetic field, a comparable dosimetric effect was observed independent of the magnetic field. Also for the five 4DCT(MRI) cases, the influence of motion was comparable for all magnetic field strengths with and without gantry angle correction. On average, the motion-induced decrease in CTVV95%from the static plan was 17.0% and 18.9% for 1.5 T and 0.5 T, respectively, and 19.9% without magnetic field.Significance.For the first time, this study investigates the combined impact of magnetic fields and respiratory motion on MR-guided proton therapy. The comparable dosimetric effects irrespective of magnetic field strength indicate that the effects of motion for future MR-guided proton therapy may not be worse than for conventional PBS proton therapy.

Proton dosimetry in a magnetic field: Measurement and calculation of magnetic field correction factors for a plane-parallel ionization chamber

BACKGROUND

The combination of magnetic resonance imaging and proton therapy offers the potential to improve cancer treatment. The magnetic field (MF)-dependent change in the dosage of ionization chambers in magnetic resonance imaging-integrated proton therapy (MRiPT) is considered by the correction factork B ⃗ , M , Q $k_{\vec{B},M,Q}$ , which needs to be determined experimentally or computed via Monte Carlo (MC) simulations.

PURPOSE

In this study,k B ⃗ , M , Q $k_{\vec{B},M,Q}$ was both measured and simulated with high accuracy for a plane-parallel ionization chamber at different clinical relevant proton energies and MF strengths.

MATERIAL AND METHODS

The dose-response of the Advanced Markus chamber (TM34045, PTW, Freiburg, Germany) irradiated with homogeneous 10 × $\times$ 10 cm2 $^2$ quasi mono-energetic fields, using 103.3, 128.4, 153.1, 223.1, and 252.7 MeV proton beams was measured in a water phantom placed in the MF of an electromagnet with MF strengths of 0.32, 0.5, and 1 T. The detector was positioned at a depth of 2 g/cm2 $^2$ in water, with chamber electrodes parallel to the MF lines and perpendicular to the proton beam incidence direction. The measurements were compared with TOPAS MC simulations utilizing COMSOL-calculated 0.32, 0.5, and 1 T MF maps of the electromagnet.k B ⃗ , M , Q $k_{\vec{B},M,Q}$ was calculated for the measurements for all energies and MF strengths based on the equation:k B ⃗ , M , Q = M Q M Q B ⃗ $k_{\vec{B},M,Q}=\frac{M_\mathrm{Q}}{M_\mathrm{Q}^{\vec{B}}}$ , whereM Q B ⃗ $M_\mathrm{Q}^{\vec{B}}$ andM Q $M_\mathrm{Q}$ were the temperature and air-pressure corrected detector readings with and without the MF, respectively. MC-based correction factors were determined ask B ⃗ , M , Q = D det D det B ⃗ $k_{\vec{B},M,Q}=\frac{D_\mathrm{det}}{D_\mathrm{det}^{\vec{B}}}$ , whereD det B ⃗ $D_\mathrm{det}^{\vec{B}}$ andD det $D_\mathrm{det}$ were the doses deposited in the air cavity of the ionization chamber model with and without the MF, respectively. Furthermore, MF effects on the chamber dosimetry are studied using MC simulations, examining the impact on the absorbed dose-to-water (D W $D_{W}$ ) and the shift in depth of the Bragg peak.

RESULTS

The detector showed a reduced dose-response for all measured energies and MF strengths, resulting in experimentally determinedk B ⃗ , M , Q $k_{\vec{B},M,Q}$ values larger than unity. For all energies and MF strengths examined,k B ⃗ , M , Q $k_{\vec{B},M,Q}$ ranged between 1.0065 and 1.0205. The dependence on the energy and the MF strength was found to be non-linear with a maximum at 1 T and 252.7 MeV. The MC simulatedk B ⃗ , M , Q $k_{\vec{B},M,Q}$ values agreed with the experimentally determined correction factors within their standard deviations with a maximum difference of 0.6%. The MC calculated impact onD W $D_{W}$ was smaller 0.2 %.

CONCLUSION

For the first time, measurements and simulations were compared for proton dosimetry within MFs using an Advanced Markus chamber. Good agreement ofk B ⃗ , M , Q $k_{\vec{B},M,Q}$ was found between experimentally determined and MC calculated values. The performed benchmarking of the MC code allows for calculatingk B ⃗ , M , Q $k_{\vec{B},M,Q}$ for various ionization chamber models, MF strengths and proton energies to generate the data needed for a proton dosimetry protocol within MFs and is, therefore, a step towards MRiPT.

MR-guided ion therapy: Detector response in magnetic fields during carbon ion irradiation

BACKGROUND

Combining carbon ion therapy with on-bed MR imaging has the potential to bring particle therapy to a new level of precision. However, the introduction of magnetic fields brings challenges for dosimetry and quality assurance. For protons, a small, but significant change in detector response was shown in the presence of magnetic fields previously. For carbon ion beams, so far no such experiments have been performed.

PURPOSE

To investigate the influence of external magnetic fields on the response of air-filled ionization chambers.

METHODS

Four commercially available ionization chambers, three thimble type (Farmer, Semiflex, and PinPoint), and a plane parallel (Bragg peak) detector were investigated. Detectors were aligned in water such that their effective point of measurement was located at 2 cm depth. Irradiations were performed using10 × 10 cm 2 $10\times 10\nobreakspace \mathrm{cm}^2$ square fields for carbon ions of 186.1, 272.5, and 402.8 MeV/u employing magnetic field strengths of 0, 0.25, 0.5, and 1 T. In addition, the detector response for protons and carbon ions was compared taking into account the secondary electron spectra and employing protons of 252.7 MeV for comparison.

RESULTS

For all four detectors, a statistically significant change in detector response, dependent on the magnetic field strength, was found. The effect was more pronounced for higher energies. The highest effects were found at 0.5 T for the PinPoint detector with a change in detector response of 1.1%. The response of different detector types appeared to be related to the cavity diameter. For proton and carbon ion irradiation with similar secondary electron spectra, the change in detector response was larger for carbon ions compared to protons.

CONCLUSION

A small, but significant dependence of the detector response was found for carbon ion irradiation in a magnetic field. The effect was found to be larger for smaller cavity diameters and at medium magnetic field strengths. Changes in detector response were more pronounced for carbon ions compared to protons.

A fast analytical dose calculation approach for MRI-guided proton therapy

Objective.Magnetic resonance (MR) is an innovative technology for online image guidance in conventional radiotherapy and is also starting to be considered for proton therapy as well. For MR-guided therapy, particularly for online plan adaptations, fast dose calculation is essential. Monte Carlo (MC) simulations, however, which are considered the gold standard for proton dose calculations, are very time-consuming. To address the need for an efficient dose calculation approach for MRI-guided proton therapy, we have developed a fast GPU-based modification of an analytical dose calculation algorithm incorporating beam deflections caused by magnetic fields.Approach.Proton beams (70-229 MeV) in orthogonal magnetic fields (0.5/1.5 T) were simulated using TOPAS-MC and central beam trajectories were extracted to generate look-up tables (LUTs) of incremental rotation angles as a function of water-equivalent depth. Beam trajectories are then reconstructed using these LUTs for the modified ray casting dose calculation. The algorithm was validated against MC in water, different materials and for four example patient cases, whereby it has also been fully incorporated into a treatment plan optimisation regime.Main results.Excellent agreement between analytical and MC dose distributions could be observed with sub-millimetre range deviations and differences in lateral shifts <2 mm even for high densities (1000 HU). 2%/2 mm gamma pass rates were comparable to the 0 T scenario and above 94.5% apart for the lung case. Further, comparable treatment plan quality could be achieved regardless of magnetic field strength.Significance.A new method for accurate and fast proton dose calculation in magnetic fields has been developed and successfully implemented for treatment plan optimisation.

Commissioning a beam line for MR-guided particle therapy assisted by in silico methods

BACKGROUND

Radiation therapy is continuously moving towards more precise dose delivery. The combination of online MR imaging and particle therapy, for example, radiation therapy using protons or carbon ions, could enable the next level of precision in radiotherapy. In particle therapy, research towards a combination of MR and particle therapy is well underway, but still far from clinical systems. The combination of high magnetic fields with particle therapy delivery poses several challenges for treatment planning, treatment workflow, dose delivery, and dosimetry.

PURPOSE

To present a workflow for commissioning of a light ion beam line with an integrated dipole magnet to perform MR-guided particle therapy (MRgPT) research, producing not only basic beam data but also magnetic field maps for accurate dose calculation. Accurate dose calculation in magnetic field environments requires high-quality magnetic field maps to compensate for magnetic-field-dependent trajectory changes and dose perturbations.

METHODS

The research beam line at MedAustron was coupled with a resistive dipole magnet positioned at the isocenter. Beam data were measured for proton and carbon ions with and without an applied magnetic field of 1 T. Laterally integrated depth-dose curves (IDC) as well as beam profiles were measured in water while beam trajectories were measured in air. Based on manufacturer data, an in silico model of the magnet was created, allowing to extract high-quality 3D magnetic field data. An existing GATE/Geant4 Monte Carlo (MC) model of the beam line was extended with the generated magnetic field data and benchmarked against experimental data.

RESULTS

A 3D magnetic field volume covering fringe fields until 50 mT was found to be sufficient for an accurate beam trajectory modeling. The effect on particle range retraction was found to be 2.3 and 0.3 mm for protons and carbon ions, respectively. Measured lateral beam offsets in water agreed within 0.4 and -0.5 mm with MC simulations for protons and carbon ions, respectively. Experimentally determined in-air beam trajectories agreed within 0.4 mm in the homogeneous magnetic field area.

CONCLUSION

The presented approach based on in silico modeling and measurements allows to commission a beam line for MRgPT while providing benchmarking data for the magnetic field modeling, required for state-of-the art dose calculation methods.

An external perpendicular magnetic field does not influence survival and DNA damage after proton and carbon ion irradiation in human cancer cells

BACKGROUND AND PURPOSE

Magnetic field effects on the radiobiological effectiveness during treatment of magnetic resonance (MRI) guided particle therapy are being debated. This study aims at assessing the influence of a perpendicular magnetic field on the biological effects in two human cancer cell lines irradiated with proton or carbon ions.

METHODS AND MATERIALS

In vitro cell irradiations were performed in water inside a perpendicular magnetic field of 0 and 1T for both protons and carbon ions. Samples were located in the center of a spread-out Bragg peak at 8cm water equivalent depth with a dose averaged linear energy transfer (LETd) of 4.2 or 83.4keV/μm for protons and carbon ions, respectively. Physical dose levels of 0, 0.5, 1, 2, 4 and 6Gy were employed. The irradiation field was shifted and laterally enlarged, to compensate for the beam deflection due to the magnetic field and ensure consistent and homogenous irradiations of the flasks. The human cancer cell lines SKMel (Melanoma) and SW1353 (chondrosarcoma) were selected which represent a high and a low (α/β)x ratio cell type. Cell survival curves were generated applying a linear-quadratic curve fit. DNA damage and DNA damage clearance were assessed via γH2AX foci quantification at 1 and 24h post radiation treatment.

RESULTS

Without a magnetic field, RBE10 values of 1.04±0.03 (SW1353) and 1.51±0.06 (SKMel) as well as RBE80 values of 0.93±0.15 (SW1353) and 2.28±0.40 (SKMel) were calculated for protons. Carbon treatments yielded RBE10 values of 1.68±0.04 (SW1353) and 2.30±0.07 (SKMel) and RBE80 values of 2.19±0.24 (SW1353) and 4.06±0.33 (SKMel). For a field strength of B=1T, RBE10 values of 1.06±0.03 (SW1353) and 1.47±0.06 (SKMel) resulted from protons, while RBE10 values of 1.70±0.05 (SW1353) and 2.37±0.08 (SKMel) were obtained for carbon ions. RBE80 values were calculated to be 1.06±0.12 (SW1353) and 2.33±0.40 (SKMel) following protons and 2.13±0.25 (SW1353) and 4.29±0.35 (SKMel) following carbon treatments. Substantially increased γH2AX foci per nucleus were found in both cell lines 1h after radiation with both ion species. At the 24h time point only carbon treated samples of both cell lines showed increased γH2AX levels. The presence of the magnetic field did neither influence the survival parameters of either cell line, nor initial DNA damage and DNA damage clearance.

CONCLUSIONS

Applying a perpendicular magnetic field did not influence the cell survival, DNA repair, nor the biological effectiveness of protons or carbon ions in two human cancer cell lines.

Integrated MRI-guided proton therapy planning: accounting for the full MRI field in a perpendicular system

PURPOSE

To present a first study on the treatment planning feasibility in perpendicular field MRI-integrated proton therapy which considers the full transport of protons from the pencil beam scanning assembly to the patient inside the MRI scanner.

METHODS

A generic proton pencil beam scanning (PBS) gantry was modelled as being integrated with a realistic split-bore MRI system in the perpendicular orientation. MRI field strengths were modeled as 0.5 T, 1 T, and 1.5 T. The PBS beam delivery and dose calculation was modeled using the TOPAS Monte Carlo toolkit coupled with matRad as the optimizer engine. A water phantom, liver and prostate plans were evaluated and optimized in the presence of the full MRI field distribution. A simple combination of gantry angle offset and small PBS nozzle skew was used to direct the proton beams along a path that closely follows the reference planning scenario, i.e. without magnetic field.

RESULTS

All planning metrics could be successfully achieved with the inclusion of gantry angle offsets in the range of 8°-29° when coupled with a PBS nozzle skew of 1.6°-4.4°. These two hardware based corrections were selected to minimize the average Euclidean distance (AED) in the beam path enabling the proton beams to travel inside the patient in a path that is close to the original path (AED smaller than 3 mm at 1.5 T). Final dose optimization, performed through further changes in the pencil beam scanning delivery, was then shown to be feasible for our selection of plans studied yielding comparable plan quality metrics to reference conditions.

CONCLUSIONS

For the first time, we have shown a robust method to account for the full proton beam deflection in a perpendicular orientation MRI-integrated proton therapy. These results support the ongoing development of the current prototype systems. This article is protected by copyright. All rights reserved.

Toward MR-integrated proton therapy: modeling the potential benefits for liver tumors

Magnetic resonance imaging (MRI)-integrated proton therapy (MRiPT) is envisioned to improve treatment quality for many cancer patients. However, given the availability of alternative image-guided strategies, its clinical need is yet to be justified. This study aims to compare the expected clinical outcomes of MRiPT with standard of practice cone-beam CT (CBCT)-guided PT, and other MR-guided methods, i.e. offline MR-guided PT and MR-linac, for treatment of liver tumors. Clinical outcomes were assessed by quantifying the dosimetric and biological impact of target margin reduction enabled by each image-guided approach. Planning target volume (PTV) margins were calculated using random and systematic setup, delineation and motion uncertainties, which were quantified by analyzing longitudinal MRI data for 10 patients with liver tumors. Proton treatment plans were created using appropriate PTV margins for each image-guided PT method. Photon plans with margins equivalent to MRiPT were generated to represent MR-linac. Normal tissue complication probabilities (NTCP) of the uninvolved liver were compared. We found that PTV margin can be reduced by 20% and 40% for offline MR-guided PT and MRiPT, respectively, compared with CBCT-guided PT. Furthermore, clinical target volume expansion could be largely alleviated when delineating on MRI rather than CT. Dosimetric implications included decreased equivalent mean dose of the uninvolved liver, i.e. up to 24.4 Gy and 27.3 Gy for offline MR-guided PT and MRiPT compared to CBCT-guided PT, respectively. Considering Child-Pugh score increase as endpoint, NTCP of the uninvolved liver was significantly decreased for MRiPT compared to CBCT-guided PT (up to 48.4%,p < 0.01), offline MR-guided PT (up to 12.9%,p < 0.01) and MR-linac (up to 30.8%,p < 0.05). Target underdose was possible in the absence of MRI-guidance (D90 reduction up to 4.2 Gy in 20% of cases). In conclusion, MRiPT has the potential to significantly reduce healthy liver toxicities in patients with liver tumors. It is superior to other image-guided techniques currently available.

Dosimetric Deviations of Bragg-Peak Position Shifts in Uniform Magnetic Fields for Magnetic Resonance Imaging-Guiding Proton Radiotherapy: A Monte Carlo Study

Objective: To investigate dosimetric deviations in scanning protons for Bragg-peak position shifts, which were caused by proton spiral tracks in an ideal uniform field of magnetic resonance (MRI) imaging-guided proton radiotherapy (MRI-IGPRT).

Methods: The FLUKA Monte-Carlo (MC) code was used to simulate the spiral tracks of protons penetrating water with initial energies of 70–270 MeV under the influence of field strength of 0.0–3.0 Tesla in commercial MRI systems. Two indexes, lateral shift (marked as WD) perpendicular to the field and a penetration-depth shift (marked as ΔDD) along the beam path, were employed for the Bragg-peak position of spiral proton track analysis. A comparison was performed between MC and classical analytical model to check the simulation results. The shape of the 2D/3D dose distribution of proton spots at the depth of Bragg-Peak was also investigated. The ratio of Gaussian-fit value between longitudinal and transverse major axes was used to indicate the asymmetric index. The skewness of asymmetry was evaluated at various dose levels by the radius ratio of circumscribed and inscribed circles by fitting a semi-ellipse circle of 2D distribution.

Results: The maximum of WD deflection is 2.82 cm while the maximum of shortening ΔDD is 0.44 cm for proton at 270 MeV/u under a magnetic field of 3.0 Tesla. The trend of WD and ΔDD from MC simulation was consistent with the analytical model, which means the reverse equation of the analytical model can be applied to determine the proper field strength of the magnet and the initial energy of the proton for the planned dose. The asymmetry of 2D/3D dose distribution under the influence of a magnetic field was increased with higher energy, and the skewness of asymmetry for one proton energy at various dose levels was also increased with a larger radius, i.e., a lower dose level.

Conclusions: The trend of the spiral proton track under a uniform magnetic field was obtained in this study using either MC simulation or the analytical model, which can provide an optimized and planned dose of the proton beam in the clinical application of MRI-IGPRT.

MR-guided proton therapy: Impact of magnetic fields on the detector response

PURPOSE

To investigate the response of detectors for proton dosimetry in the presence of magnetic fields.

MATERIAL AND METHODS

Four ionization chambers (ICs), two thimble-type and two plane-parallel-type, and a diamond detector were investigated. All detectors were irradiated with homogeneous single-energy-layer fields, using 252.7 MeV proton beams. A Farmer IC was additionally irradiated in the same geometrical configuration, but with a lower nominal energy of 97.4 MeV. The beams were subjected to magnetic field strengths of 0, 0.25, 0.5, 0.75, and 1 T produced by a research dipole magnet placed at the room's isocenter. Detectors were positioned at 2 cm water equivalent depth, with their stem perpendicular to both the magnetic field lines and the proton beam's central axis, in the direction of the Lorentz force. Normality and two sample statistical Student's t tests were performed to assess the influence of the magnetic field on the detectors' responses.

RESULTS

For all detectors, a small but significant magnetic field-dependent change of their response was found. Observed differences compared to the no magnetic field case ranged from +0.5% to -0.7%. The magnetic field dependence was found to be nonlinear and highest between 0.25 and 0.5 T for 252.7 MeV proton beams. A different variation of the Farmer chamber response with magnetic field strength was observed for irradiations using lower energy (97.4 MeV) protons. The largest magnetic field effects were observed for plane-parallel ionization chambers.

CONCLUSION

Small magnetic field-dependent changes in the detector response were identified, which should be corrected for dosimetric applications.

Technical Note: Design and commissioning of a water phantom for proton dosimetry in magnetic fields

Purpose

To design and commission a water phantom suitable for constrained environments and magnetic fields for magnetic resonance (MR)‐guided proton therapy.

Methods

A phantom was designed, to enable precise, remote controlled detector positioning in water within the constrained environment of a magnet for MR‐guided proton therapy. The phantom consists of a PMMA enclosure whose outer dimensions of 81×40×12.5cm3 were chosen to optimize space usage inside the 13.5‐cm bore gap of the magnet. The moving mechanism is based on a low‐height H‐shaped non‐ferromagnetic belt drive, driven by stepper motors located outside of the magnetic field. The control system and the associated electronics were designed in house, with similar features as available in commercial water phantoms. Reproducibility as well as accuracy of the phantom positioning were tested using a high‐precision Leica AT 402 laser tracker. Laterally integrated depth dose curves and lateral beam profiles at three depths were acquired repeatedly for a 148.2 MeV proton beam in water.

Results

The phantom was successfully operated with and without applied magnetic fields. For complex movements, a positioning uncertainty within 0.16 mm was found with an absolute accuracy typically below 0.3 mm. Laterally integrated depth dose curves agreed within 0.1 mm with data taken using a commercial water phantom. The lateral beam offset determined from beam profile measurements agreed well with data from Monte Carlo simulations.

Conclusion

The phantom is optimally suited for detector positioning and dosimetric experiments within constrained environments in high magnetic fields.

Characterization of magnetic interference and image artefacts during simultaneous in-beam MR imaging and proton pencil beam scanning

For the first time, a low-field open magnetic resonance (MR) scanner was combined with a proton pencil beam scanning (PBS) research beamline. The aim of this study was to characterize the magnetic fringe fields produced by the PBS system and measure their effects on MR image quality during simultaneous PBS irradiation and image acquisition. A magnetic field camera measured the change in central resonance frequency (Δf _res) and magnetic field homogeneity (ΔMFH) of the B₀ field of the MR scanner during operation of the beam transport and scanning magnets. The beam energy was varied between 70 - 220 MeV and beam scanning was performed along the central horizontal and vertical axis of a 48 × 24 cm² radiation field. The time structure of the scanning magnets' fringe fields was simultaneously recorded by a tri-axial Hall probe. MR imaging experiments were conducted using the ACR (American College of Radiology) Small MRI Phantom and a spoiled gradient echo pulse sequence during simultaneous volumetric irradiation. Computer simulations were performed to predict the effects of B ₀ field perturbations due to PBS irradiation on MR image formation in k-space. Setting the beam transport magnets, horizontal and vertical scanning magnets resulted in a maximum Δf _res of 50, 235 and 4 Hz, respectively. The ΔMFH was less than 3 parts per million for all measurements. MR images acquired during beam energy variation and vertical beam scanning showed no visual loss in image quality. However, MR images acquired during horizontal beam scanning showed severe coherent ghosting artefacts in phase encoding direction. Both simulated and measured k-space phase maps prove that these artefacts are caused by phase-offsets. This study shows first experimental evidence that simultaneous in-beam MR imaging during proton PBS irradiation is subject to severe loss of image quality in the absence of magnetic decoupling between the PBS and MR system.

Plasmid-loadable magnetic/ultrasound-responsive nanodroplets with a SPIO-NP dispersed perfluoropentane core and lipid shell for tumor-targeted intracellular plasmid delivery

Using ultrasound activating contrast agents to induce sonoporation is a potential strategy for effective lesion-targeted gene delivery. Previous reports have proven that submicron nanodroplets have a better advantage than microbubbles in that they can pass through tumor vasculature endothelial gaps by passive targeting; however, they cannot achieve an adequate dose in tumors to facilitate ultrasound-enhanced gene delivery. Additionally, a few studies focused on delivering macromolecular genetic materials (i.e. overexpression plasmid and CRISPR plasmid) have presented more unique advantages than small-molecular genetic materials (i.e. miRNA mimics, siRNA and shRNA etc.), such as enhancing the expression of target genes with long-term effectiveness. Thereby, we constructed novel plasmid-loadable magnetic/ultrasound-responsive nanodroplets, where superparamagnetic iron oxide nanoparticle dispersed perfluoropentane was encapsulated with lipids to which plasmids could be adhered, and branched polyethylenimine was used to protect the plasmids from enzymolysis. Furthermore, in vitro and in vivo studies were performed to verify the magnetic tumor-targeting ability of the plasmid-loadable magnetic/ultrasound-responsive nanodroplets and focused ultrasound enhanced intracellular plasmid delivery. The plasmid-loadable magnetic/ultrasound-responsive nanodroplets, carrying 16-19 plasmids per droplet, had desirable diameters less than 300 nm, and integrated the merits of excellent magnetic targeting capabilities and phase transition sensitivity to focused ultrasound. Under programmable focused ultrasound exposure, the plasmid-loadable magnetic/ultrasound-responsive nanodroplets underwent a phase-transition into echogenic microbubbles and the subsequent inertial cavitation of the microbubbles achieved an ∼40% in vitro plasmid delivery efficiency. Following intravenous administration, T2-weighted magnet resonance imaging, scanning electron microscopy and inductively coupled plasma optical emission spectrometry of the tumors showed significantly enhanced intratumoral accumulation of the plasmid-loadable magnetic/ultrasound-responsive nanodroplets under an external magnetic field. And a GFP ELISA assay and immunofluorescence staining indicated that focused ultrasound-induced inertial cavitation of the plasmid-loadable magnetic/ultrasound-responsive nanodroplets significantly enhanced the intracellular delivery of plasmids within the tumor after magnet-assisted accumulation, while only lower GFP levels were observed in the tumors on applying focused ultrasound or an external magnet alone. Taken together, utilizing the excellent plasmid-loadable magnetic/ultrasound-responsive nanodroplets combined with magnetism and ultrasound could efficiently deliver plasmids to cancer cells, which could be a potential strategy for macromolecular genetic material delivery in the clinic to treat cancer.

Multipotent miRNA Sponge-Loaded Magnetic Nanodroplets with Ultrasound/Magnet-Assisted Delivery for Hepatocellular Carcinoma Therapy

Gene therapy is likely to be the most promising way to tackle cancer, while defects in molecular strategies and delivery systems have led to an impasse in clinical application. Here, it is found that onco-miRNAs of the miR-515 and -449 families were upregulated in hepatocellular carcinoma (HCC), and the sponge targeting miR-515 family had a significant probability to suppress cancer cell proliferation. Then, we constructed non-toxic sponge-loaded magnetic nanodroplets containing 20% C6F14 (SLMNDs-20%) that are incorporated with fluorinated superparamagnetic iron oxide nanoparticles enhancing external magnetism-assisted targeting and enabling a direct visualization of SLMNDs-20% distribution in vivo via magnetic resonance imaging monitoring. SLMNDs-20% could be vaporized by programmable focused ultrasound (FUS) activation, achieving ∼45% in vitro sponge delivery efficiency and significantly enhancing in vivo sponge delivery without a clear apoptosis. Moreover, the sponge-1-carrying SLMNDs-20% could effectively suppress proliferation of xenograft HCC after FUS exposure because sponge-1-suppressing onco-miR-515 enhanced the expression of anti-oncogenes (P21, CD22, TIMP1, NFKB, and E-cadherin) in cancer cells. The current results indicated that ultrasonic cavitation-inducing sonoporation enhanced the intracellular delivery of sponge-1 using SLMNDs-20% after magnetic-assisted accumulation, which was a therapeutic approach to inhibit HCC progression.

MR-guided proton therapy: a review and a preview

Background

The targeting accuracy of proton therapy (PT) for moving soft-tissue tumours is expected to greatly improve by real-time magnetic resonance imaging (MRI) guidance. The integration of MRI and PT at the treatment isocenter would offer the opportunity of combining the unparalleled soft-tissue contrast and real-time imaging capabilities of MRI with the most conformal dose distribution and best dose steering capability provided by modern PT. However, hybrid systems for MR-integrated PT (MRiPT) have not been realized so far due to a number of hitherto open technological challenges. In recent years, various research groups have started addressing these challenges and exploring the technical feasibility and clinical potential of MRiPT. The aim of this contribution is to review the different aspects of MRiPT, to report on the status quo and to identify important future research topics.

Methods

Four aspects currently under study and their future directions are discussed: modelling and experimental investigations of electromagnetic interactions between the MRI and PT systems, integration of MRiPT workflows in clinical facilities, proton dose calculation algorithms in magnetic fields, and MRI-only based proton treatment planning approaches.

Conclusions

Although MRiPT is still in its infancy, significant progress on all four aspects has been made, showing promising results that justify further efforts for research and development to be undertaken. First non-clinical research solutions have recently been realized and are being thoroughly characterized. The prospect that first prototype MRiPT systems for clinical use will likely exist within the next 5 to 10 years seems realistic, but requires significant work to be performed by collaborative efforts of research groups and industrial partners.

The effects of a transverse magnetic field on the dose enhancement of nanoparticles in a proton beam: a Monte Carlo simulation

High-Z nanoparticles (NP) as radio-sensitization agents provide the feasibility of dose localization within the tumor in radiotherapy. Dose enhancement of NPs in the presence of a magnetic field (MF) could be challenged when magnetic resonance imaging (MRI) systems are used as an image-guided system. The MF can influence dose enhancement of NPs at their interfaces and surrounding medium and affect their dose deposition behavior. In the TOPAS Monte Carlo code, gold nanoparticle (GNP) and superparamagnetic iron oxide nanoparticle (SPION) were irradiated using 70 and 150 MeV proton beams, in presence of transverse MF strengths with 0, 1, 3, and 7 T. The changes in the liberated secondary electrons from NPs and their dose enhancement ratio (DER), magnetic dose enhancement ratio (MDER), and angular dose distribution in 10 nm shell thicknesses up to 500 nanometers from their centers were measured. The central plane of NPs was considered as a scorer. Its thickness was 2 nm and divided into 6-degree sectors with 10 nm radial length. The dose deposition in this voxelated scorer was calculated. The values of the deposited doses around NPs decrease rapidly while the DERs resulted from the secondary electrons are increased. MDERs are changed within [Formula: see text] and [Formula: see text] for 20 and 50 nm radius NPs, respectively. The variation in the angular dose distribution around a singular NP was not considerable when different MF strengths were applied. The dose values in the voxelated central plane show very similar results for the same NPs types in the different MF strengths. The typically used MF in the MRI systems would not considerably affect the energy deposition behavior of the secondary electrons produced in the interaction of proton beam with NPs, at least in the near vicinity of NPs. The DERs of NPs in a water medium resulted from emerged secondary electrons, experience a low degree of perturbation in the presence of an MF. The results of this study show that the NPs as dose enhancement agents can also be used in an MF without pronounced modification in their efficacy.

Benchmarking a GATE/Geant4 Monte Carlo model for proton beams in magnetic fields

Purpose

Magnetic resonance guidance in proton therapy (MRPT) is expected to improve its current performance. The combination of magnetic fields with clinical proton beam lines poses several challenges for dosimetry, treatment planning and dose delivery. Proton beams are deflected by magnetic fields causing considerable changes in beam trajectories and also a retraction of the Bragg peak positions. A proper prediction and compensation of these effects is essential to ensure accurate dose calculations. This work aims to develop and benchmark a Monte Carlo (MC) beam model for dose calculation of MRPT for static magnetic fields up to 1 T.

Methods

Proton beam interactions with magnetic fields were simulated using the GATE/Geant4 toolkit. The transport of charged particle in custom 3D magnetic field maps was implemented for the first time in GATE. Validation experiments were done using a horizontal proton pencil beam scanning system with energies between 62.4 and 252.7 MeV and a large gap dipole magnet (B = 0–1 T), positioned at the isocenter and creating magnetic fields transverse to the beam direction. Dose was measured with Gafchromic EBT3 films within a homogeneous PMMA phantom without and with bone and tissue equivalent material slab inserts. Linear energy transfer (LET) quenching of EBT3 films was corrected using a linear model on dose‐averaged LET method to ensure a realistic dosimetric comparison between simulations and experiments. Planar dose distributions were measured with the films in two different configurations: parallel and transverse to the beam direction using single energy fields and spread‐out Bragg peaks. The MC model was benchmarked against lateral deflections and spot sizes in air of single beams measured with a Lynx PT detector, as well as dose distributions using EBT3 films. Experimental and calculated dose distributions were compared to test the accuracy of the model.

Results

Measured proton beam deflections in air at distances of 465, 665, and 1155 mm behind the isocenter after passing the magnetic field region agreed with MC‐predicted values within 4 mm. Differences between calculated and measured beam full width at half maximum (FWHM) were lower than 2 mm. For the homogeneous phantom, measured and simulated in‐depth dose profiles showed range and average dose differences below 0.2 mm and 1.2%, respectively. Simulated central beam positions and widths differed <1 mm to the measurements with films. For both heterogenous phantoms, differences within 1 mm between measured and simulated central beam positions and widths were obtained, confirming a good agreement of the MC model.

Conclusions

A GATE/Geant4 beam model for protons interacting with magnetic fields up to 1 T was developed and benchmarked to experimental data. For the first time, the GATE/Geant4 model was successfully validated not only for single energy beams, but for SOBP, in homogeneous and heterogeneous phantoms. EBT3 film dosimetry demonstrated to be a powerful dosimetric tool, once the film response function is LET corrected, for measurements in‐line and transverse to the beam direction in magnetic fields. The proposed MC beam model is foreseen to support treatment planning and quality assurance (QA) activities toward MRPT.

Carbon ion and proton beam irradiation of a normal human TK6 lymphoblastoid cell line within a magnetic field of 1.0 tesla

Background

Considering the increasing simultaneous application of magnetic resonance imaging (MRI) for more precise photon radiotherapy, it will be likely for particle radiotherapy to adopt MRI for future image guiding. It will then be imperative to evaluate the potential biological effects of a magnetic field (MF) on particle irradiation. This study explores such effects on the highly radiosensitive TK6 lymphoblastoid human cell line.

Methods

The following three parameters were measured after irradiation with either carbon ion or proton beams using spread out Bragg peaks and applying different doses within a perpendicular 1.0 T MF: (1) cell survival fraction (14 days postirradiation), (2) treatment-specific apoptosis, which was determined through the measurement of population in the sub-G1 phase, and (3) cell cycle progression by means of flow cytometry. These were compared to the same parameters measured without an MF.

Results

The clonogenic assay in both treatment groups showed almost identical survival curves with overlapping error bars. The calculated α values with and without an MF were 2.18 (σ=0.245) and 2.17 (σ=0.234) for carbon ions and 1.08 (σ=0.138) and 1.13 (σ=0.0679) for protons, respectively. Similarly, the treatment-specific apoptosis and cell cycle progression showed almost identical curves with overlapping error bars. A two-sample, unpooled t-test analysis was implemented for comparison of all mean values and showed p-values >0.05.

Conclusion

No statistically significant difference in biological response of the TK6 cells was observed when they were irradiated using spreadout Bragg peaks within a perpendicular 1.0 T MF as compared to those, which received the same dose without the MF. This should serve as another supporting piece of evidence toward the implementation of MRI in particle radiotherapy, though further research is necessary.

MRI-guided proton therapy planning: accounting for an inline MRI fringe field

MRI-guided proton therapy is being pursued for its promise to provide a more conformal, accurate proton therapy. However, the presence of the magnetic field imposes a challenge for the beam delivery as protons are deflected due to the Lorenz force. In this study, the impact of realistic inline MRI fringe field on IMPT plan delivery is investigated for a water phantom, liver tumor and prostate cancer differing in target volume, shape, and field configuration using Monte Carlo simulations. A method to correct for the shift of the beam spot positions in the presence of the inline magnetic field is presented. Results show that when not accounting for the effect of the magnetic field on the pencil beam delivery, the spot positions are substantially shifted and the quality of delivered plans is significantly deteriorated leading to dose inhomogeneities and creation of hot and cold spots. However, by correcting the pencil beam delivery, the dose quality of the IMPT plans is restored to a high degree. Nevertheless, adaptation of beam delivery alone is not robust regarding different treatment sites. By fully accounting during plan optimization for the dose distortions caused by the fringe and imaging fields, highly conformal IMPT plans are achieved. These results demonstrate proton pencil beam scanning and treatment planning can be adapted for precise delivery of state-of-the-art IMPT plans in MR-guided proton therapy in the presence of an inline MRI fringe field.

Effect of External Magnetic Fields on Biological Effectiveness of Proton Beams

PURPOSE

The purpose is to verify experimentally whether application of magnetic fields longitudinal and perpendicular to a proton beam alters the biological effectiveness of the radiation.

METHODS AND MATERIALS

Proton beams with linear energy transfer of 1.1 and 3.3 keV/μm irradiated human cancer and normal cells under a longitudinal (perpendicular) magnetic field of B_L (B_P) = 0, 0.3, or 0.6 T. Cell survival curves were constructed to evaluate the effects of the magnetic fields on the biological effectiveness. The ratio of dose that would result in a survival fraction of 10% without the magnetic field D_wo to the dose with the magnetic field D_w, R₁₀ = D_wo/D_w, was determined for each cell line and magnetic field.

RESULTS

For cancer cells exposed to the 1.1- (3.3-) keV/μm proton beams, R₁₀s were increased to 1.10 ± 0.07 (1.11 ± 0.07) and 1.11 ± 0.07 (1.12 ± 0.07) by the longitudinal magnetic fields of B_L = 0.3 and 0.6 T, respectively. For normal cells, R₁₀s were increased to 1.13 ± 0.06 (1.17 ± 0.06) and 1.17 ± 0.06 (1.30 ± 0.06) by the B_Ls. In contrast, R₁₀s were not changed significantly from 1 by the perpendicular magnetic fields of B_P = 0.3 and 0.6 T for both cancer and normal cells exposed to 1.1- and 3.3-keV/μm proton beams.

CONCLUSIONS

The biological effectiveness of proton beams was significantly enhanced by longitudinal magnetic fields of B_L = 0.3 and 0.6 T, whereas the biological effectiveness was not altered by perpendicular magnetic fields of the same strengths. This enhancement effect should be taken into account in magnetic resonance imaging guided proton therapy with a longitudinal magnetic field.

Real-time intrafraction motion monitoring in external beam radiotherapy

Radiotherapy (RT) aims to deliver a spatially conformal dose of radiation to tumours while maximizing the dose sparing to healthy tissues. However, the internal patient anatomy is constantly moving due to respiratory, cardiac, gastrointestinal and urinary activity. The long term goal of the RT community to 'see what we treat, as we treat' and to act on this information instantaneously has resulted in rapid technological innovation. Specialized treatment machines, such as robotic or gimbal-steered linear accelerators (linac) with in-room imaging suites, have been developed specifically for real-time treatment adaptation. Additional equipment, such as stereoscopic kilovoltage (kV) imaging, ultrasound transducers and electromagnetic transponders, has been developed for intrafraction motion monitoring on conventional linacs. Magnetic resonance imaging (MRI) has been integrated with cobalt treatment units and more recently with linacs. In addition to hardware innovation, software development has played a substantial role in the development of motion monitoring methods based on respiratory motion surrogates and planar kV or Megavoltage (MV) imaging that is available on standard equipped linacs. In this paper, we review and compare the different intrafraction motion monitoring methods proposed in the literature and demonstrated in real-time on clinical data as well as their possible future developments. We then discuss general considerations on validation and quality assurance for clinical implementation. Besides photon RT, particle therapy is increasingly used to treat moving targets. However, transferring motion monitoring technologies from linacs to particle beam lines presents substantial challenges. Lessons learned from the implementation of real-time intrafraction monitoring for photon RT will be used as a basis to discuss the implementation of these methods for particle RT.

Characterization of EBT3 radiochromic films for dosimetry of proton beams in the presence of magnetic fields

PURPOSE

Radiochromic film dosimetry is extensively used for quality assurance in photon and proton beam therapy. So far, GafchromicTM EBT3 film appears as a strong candidate to be used in future magnetic resonance (MR) based therapy systems. The response of Gafchromic EBT3 films in the presence of magnetic fields has already been addressed for different MR-linacs systems. However, a detailed evaluation of the influence of external magnetic fields on the film response and calibration curves for proton therapy has not yet been reported. This study aims to determine the dose responses of EBT3 films for clinical proton beams exposed to magnetic field strengths up to 1 T in order to investigate the feasibility of EBT3 film as an accurate dosimetric tool for a future MR particle therapy system (MRPT).

METHODS

The dosimetric characteristics of EBT3 films were studied for a proton beam passing through magnetic field strengths of B = 0, 0.5, and 1 T. Absorbed dose calibration and measurements were performed using clinical proton beams in the nominal energy range of 62.4-252.6 MeV. Irradiations were done using an in-house developed PMMA slab phantom placed in the center of a dipole research magnet. Monte Carlo (MC) simulations using the GATE/Geant4 toolkit were performed to predict the effect of magnetic fields on the energy deposited by proton beams in the phantom. Planned and measured doses from 3D box cube irradiations were compared to assess the accuracy of the dosimetric method using EBT3 films with/without the external magnetic field.

RESULTS

Neither for the mean pixel value nor for the net optical density, any significant deviations were observed due to the presence of an external magnetic field (B ≤ 1T) for doses up to 10 Gy. Dose-response curves for the red channel were fitted by a three-parameter function for the field-free case and for B = 1T, showing for both cases an R-square coefficient of unity and almost identical fitting parameters. Independently of the magnetic field, EBT3 films showed an under-response as high as 8% in the Bragg peak region, similarly to previously reported effects for particle therapy. No noticeable influence of the magnetic field strength was observed on the quenching effect of the EBT3 films.

CONCLUSIONS

For the first time detailed absorbed dose calibrations of EBT3 films for proton beams in magnetic field regions were performed. Results showed that EBT3 films represent an attractive solution for the dosimetry of a future MRPT system. As film response functions for protons are not affected by the magnetic field strenght, they can be used for further investigations to evaluate the dosimetric effects induced due to particle beams bending in magnetic fields regions.

Technical Note: Effect of magnetic fields on the microdosimetry of carbon-ion beams

PURPOSE

To investigate the influence of magnetic fields on the microdosimetry of carbon-ion beams and the scaling effect of tissue equivalent proportional counter (TEPC) defined as the change of energy deposition in a TEPC from that in a microscopic scale region of interest due to the presence of a magnetic field in combination with the TEPC larger physical dimensions.

METHODS

Geant4-based Monte Carlo simulations were conducted to calculate the microdosimetric quantities for carbon-ion beams with different initial energies (10-290 MeV/u) under magnetic fields of various strengths (0.5-3 T). The calculations were performed for a 1 μm spherical volume made of tissue, and for spherical TEPCs of 1 and 10 mm in diameter. Then, values of dose-averaged lineal energy (y_D ) were acquired for the different scenarios to analyze the effect of magnetic fields on the microdosimetry of carbon-ion beams and the scaling effect of TEPC.

RESULTS

The y_D values and lineal energy spectra in the 1 μm spherical tissue volume for the scenarios without magnetic field and with magnetic fields of different strengths and directions remained nearly the same for the various energy carbon-ion beams. However, compared with those of the 1 μm spherical tissue volume, an increase of of y_D values and an obvious shift of the lineal energy spectra for the TEPCs of 1 and 10 mm in diameter under magnetic fields were found.

CONCLUSIONS

The application of magnetic fields under 3 T has no significant influence on the microdosimetric results of carbon-ion beams. However, there is definitely a scaling effect when using TEPC for microdosimetric study, which alters the reading of TEPC in the presence of magnetic fields. Novel methods to correct the reading of TEPC or scaling effect-resistant microdosimetric measurement detectors are urgently needed to perform experimental microdosimetric studies under magnetic fields.

Influence of a perpendicular magnetic field on biological effectiveness of carbon-ion beams

Purpose: Our previous study revealed that the application of a magnetic field longitudinal to a carbon-ion beam of 0.1 ≤ B_//≤ 0.6 T enhances the biological effectiveness of the radiation. The purpose of this study is to experimentally verify whether the application of a magnetic field perpendicular to the beam also alters the biological effectiveness. Methods and materials: Most experimental conditions other than the magnetic field direction were the same as those used in the previous study to allow comparison of their results. Human cancer and normal cells were exposed to low (12 keV/μm) and high (50 keV/μm) linear energy transfer (LET) carbon-ion beams under the perpendicular magnetic fields of B⊥ = 0, 0.15, 0.3, or 0.6 T generated by a dipole magnet. The effects of the magnetic fields on the biological effectiveness were evaluated by clonogenic cell survival. Doses that would result in the survival of 10%, D₁₀s, were determined for the exposures and analyzed using Student's t-tests. Results: For both cancer and normal cells treated by low- and high-LET carbon-ion beams, the D₁₀s measured in the presence of the perpendicular magnetic fields of B⊥ ≥ 0.15 T were not statistically different (p ≫ .05) from the D₁₀s measured in the absence of the magnetic fields, B⊥ = 0 T. Conclusions: Exposure of human cancer and normal cells to the perpendicular magnetic fields of B⊥ ≤ 0.6 T did not alter significantly the biological effectiveness of the carbon-ion beams, unlike the exposure to longitudinal magnetic fields of the same strength. Although the mechanisms underlying the observed results still require further exploration, these findings indicate that the influence of the magnetic field on biological effectiveness of the carbon-ion beam depends on the applied field direction with respect to the beam.

Proton beam electron return effect: Monte Carlo simulations and experimental verification

Proton therapy (PT) is expected to benefit from integration with magnetic resonance (MR) imaging. However, the magnetic field distorts the dose distribution and enhances the dose at tissue-air interfaces by the electron return effect (ERE). The objectives were (a) to provide experimental evidence for the ERE in proton beams and (b) to systematically characterise the dependence of the dose enhancement ratio (DER) on magnetic field strength, orientation, proton energy and voxel size by computer simulations. EBT3 films were irradiated with 200 MeV protons with and without a 0.92 T transverse field of a permanent magnet to determine the DER at effective measurement depths of 0.156 and 0.467 mm from an air interface. High-resolution Monte Carlo simulations were performed to reproduce the irradiation experiments and to calculate the DER for proton energies between 50-200 MeV and magnetic field strengths between 0.35-3 T as function of distance from the air interface. Voxel sizes of 0.05, 0.5 and 1 mm were analysed. DERs of (2.2  ±  0.4)% and (0.5  ±  0.6)% were measured at 0.156 and 0.467 mm from the air interface, respectively. Measurements and simulations agreed within 0.15%. For a 200 MeV proton beam, the maximum DER in 0.05 mm voxels increased with magnetic field strength from 2.6% to 8.2% between 0.35 and 1.5 T, respectively. For a 1.0 T magnetic field, maximum DER increased from 3.2% to 7.6% between 50 and 200 MeV, respectively. Voxel sizes of 0.5 and 1 mm resulted in maximum DER values of 2.6% and 1.4%, respectively. The ERE for proton beams in transverse magnetic fields is measurable. The local dose enhancement is significant, well predictable, decreases rapidly with distance from the air interface, and is negligible beyond 1 mm depth. Its impact on air-filled ionisation chambers and porous tissues (e.g. lung) needs to be considered.

Current state and future applications of radiological image guidance for particle therapy

In this review paper, we first give a short overview of radiological image guidance in photon radiotherapy, placing emphasis on the fact that linac based radiotherapy has outpaced particle therapy in the adoption of volumetric image guidance. While cone beam computed tomography (CBCT) has been an established technique in linac treatment rooms for almost two decades, the widespread adoption of volumetric image guidance in particle therapy, whether by means of CBCT or in-room CT imaging, is recent. This lag may be attributable to the bespoke nature and lower number of particle therapy installations, as well as the differences in geometry between those installations and linac treatment rooms. In addition, for particle therapy the so called shift invariance of the dose distribution rarely applies. An overview of the different volumetric image guidance solutions found at modern particle therapy facilities is provided, covering gantry, nozzle, C-arm, and couch-mounted CBCT as well different in-room CT configurations. A summary of the use of in-room volumetric imaging data beyond anatomy-based positioning is also presented as well as the necessary corrections to CBCT images for accurate water equivalent thickness calculation. Finally, the use of non-ionizing imaging modalities is discussed.

Integrating a low-field open MR scanner with a static proton research beam line: proof of concept

On-line image guidance using magnetic resonance (MR) imaging is expected to improve the targeting accuracy of proton therapy. However, to date no combined system exists. In this study, for the first time a low-field open MR scanner was integrated with a static proton research beam line to test the feasibility of simultaneous irradiation and imaging. The field-of-view of the MR scanner was aligned with the beam by taking into account the Lorentz force induced beam deflection. Various imaging sequences for extremities were performed on a healthy volunteer and on a patient with a soft-tissue sarcoma of the upper arm, both with the proton beam line switched off. T ₁-weighted spin echo images of a tissue-mimicking phantom were acquired without beam, with energised beam line magnets and during proton irradiation. Beam profiles were acquired for the MR scanner's static magnetic field alone and in combination with the dynamic gradient fields during the acquisition of different imaging sequences. It was shown that MR imaging is feasible in the electromagnetically contaminated environment of a proton therapy facility. The observed quality of the anatomical MR images was rated to be sufficient for target volume definition and positioning. The tissue-mimicking phantom showed no visible beam-induced image degradation. The beam profiles depicted no influence due to the dynamic gradient fields of the imaging sequences. This study proves that simultaneous irradiation and in-beam MR imaging is technically feasible with a low-field MR scanner integrated with a static proton research beam line.

Clinical implementations of 4D pencil beam scanned particle therapy: Report on the 4D treatment planning workshop 2016 and 2017

In 2016 and 2017, the 8th and 9th 4D treatment planning workshop took place in Groningen (the Netherlands) and Vienna (Austria), respectively. This annual workshop brings together international experts to discuss research, advances in clinical implementation as well as problems and challenges in 4D treatment planning, mainly in spot scanned proton therapy. In the last two years several aspects like treatment planning, beam delivery, Monte Carlo simulations, motion modeling and monitoring, QA phantoms as well as 4D imaging were thoroughly discussed. This report provides an overview of discussed topics, recent findings and literature review from the last two years. Its main focus is to highlight translation of 4D research into clinical practice and to discuss remaining challenges and pitfalls that still need to be addressed and to be overcome.

Technical Note: Experimental verification of magnetic field-induced beam deflection and Bragg peak displacement for MR-integrated proton therapy

PURPOSE

Given its sensitivity to anatomical variations, proton therapy is expected to benefit greatly from integration with magnetic resonance imaging for online anatomy monitoring during irradiation. Such an integration raises several challenges, as both systems mutually interact. The proton beam will experience quasi-continuous energy loss and energy-dependent electromagnetic deflection at the same time, giving rise to a deflected beam trajectory and an altered dose distribution with a displaced Bragg peak. So far, these effects have only been predicted using Monte Carlo and analytical models, but no clear consensus has been reached and experimental benchmark data are lacking. We measured proton beam trajectories and Bragg peak displacement in a homogeneous phantom placed inside a magnetic field and compared them to simulations.

METHODS

Planar dose distributions of proton pencil beams (80-180 MeV) traversing the field of a 0.95 T NdFeB permanent magnet while depositing energy in a PMMA slab phantom were measured using EBT3 radiochromic films and simulated using the Geant4 toolkit. Deflected beam trajectories and the Bragg peak displacement were extracted from the measured planar dose distributions and compared against the simulations.

RESULTS

The lateral beam deflection was clearly visible on the EBT3 films and ranged from 1 to 10 mm for 80 to 180 MeV, respectively. Simulated and measured beam trajectories and Bragg peak displacement agreed within 0.8 mm for all studied proton energies.

CONCLUSIONS

These results prove that the magnetic field-induced Bragg peak displacement is both measurable and accurately predictable in a homogeneous phantom at 0.95 T, and allows Monte Carlo simulations to be used as gold standard for proton beam trajectory prediction in similar frameworks for MR-integrated proton therapy.

A pencil beam algorithm for magnetic resonance image-guided proton therapy

Purpose

The feasibility of magnetic resonance image (MRI)‐based proton therapy is based, among several other factors, on the implementation of appropriate extensions on current dose calculation methods. This work aims to develop a pencil beam algorithm (PBA) for dose calculation of proton beams within magnetic field regions of up to 3 T.

Methods

Monte Carlo (MC) simulations using the GATE 7.1/GEANT4.9.4p02 toolkit were performed to generate calibration and benchmarking data for the PBA. Dose distributions from proton beams in the clinical required energy range 60–250 MeV impinging on a 400 × 400 × 400 mm³ water phantom and transverse magnetic fields ranging from 0 to 3 T were considered. Energy depositions in homogeneous and heterogeneous phantoms filled with water, adipose, bone, and air were evaluated for proton energies of 80, 150, and 240 MeV, combining a trajectory calculation method and look‐up tables (LUT). A novel parametrization model, using independent tailed Gauss fitting functions, was employed to describe the nonsymmetric shape of lateral beam profiles. Integrated depth‐dose curves (IDD), lateral dose profiles, and two‐dimensional dose distributions calculated with the PBA were compared with results from MC simulations to assess the performance of the algorithm. A gamma index criterion of 2%/2 mm was used for analysis.

Results

A close to perfect agreement was observed for PB‐based dose calculations in water in magnetic fields of 0.5, 1.5, and 3 T. IDD functions showed differences between the PBA and MC of less than 0.1% before the Bragg peak, and deviations of 2–8% in the distal energy falloff region. Gamma index pass rates higher than 99% and mean values lower than 0.1 were encountered for all analyzed configurations. For homogeneous phantoms, only the full bone configuration offered deviations in the Bragg peak position of up to 1.7% and overestimations of the lateral beam spot width for high‐energy protons and magnetic field intensities. An excellent agreement between PBA and MC dose calculation was also achieved using slab‐like and lateral heterogeneous phantoms, with gamma index passing rates above 98% and mean values between 0.1 and 0.2. As expected, agreement reduced for high‐energy protons and high‐intensity magnetic fields, although results remained good enough to be considered for future implementation in clinical practice.

Conclusions

The proposed pencil beam algorithm for protons can accurately account for dose distortion effects induced by external magnetic fields. The application of an analytical model for dose estimation and corrections reduces the calculation times considerably, making the presented PBA a suitable candidate for integration in a treatment planning system.

Analytical investigation of magnetic field effects on Proton lateral deflection and penetrating depth in the water phantom: A relativistic approach

Background

Integrated proton therapy - MRI systems are capable of delivering high doses to the target tissues near sensitive organs and achieve better therapeutic results; however, the applied magnetic field for imaging, influences the protons path, changes the penetration depth and deflects the particles, laterally, leading to dose distribution variations.

Objective

To determine the effects of a magnetic field on the range and the lateral deflection of protons, analytically.

Methods

An analytical survey based on protons energy and range power law relation, without using small angle assumption was done. The penetration depth and lateral deflection of protons with therapeutic energy ranges 60–250 MeV in the presence of uniform magnetic fields of 0–10T intensities, were calculated analytically. Calculations were done for relativistic conditions with Mathematica software version 7.0, and MATLAB 7.0 was applied to plot curves and curve fittings.

Results

In the presence of a magnetic field, the depth of Bragg peak was decreased and it was shifted laterally. A second order polynomial model with power equation for its coefficients and a power model with quadratic polynomial coefficients predicted the maximum lateral deflection (y_max) and maximum penetration depth (z_max) variations with energy and magnetic field intensity, respectively.

Conclusion

The applied correction for deflection angle will give more reliable results in initial energy of 250 MeV and 3T magnetic field intensity. For lower energies and magnetic field intensities the differences are negligible, clinically.