Assessment of the Sun Nuclear ArcCHECK to detect errors in 6MV FFF VMAT delivery of brain SABR using ROC analysis.

Institutions use a range of different detector systems for patient-specific quality assurance (QA) measurements conducted to assure that the dose delivered by a patient's radiotherapy treatment plan matches the calculated dose distribution. However, the ability of different detectors to detect errors from different sources is often unreported. This study contains a systematic evaluation of Sun Nuclear's ArcCHECK in terms of the detectability of potential machine-related treatment errors. The five investigated sources of error were multileaf collimator (MLC) leaf positions, gantry angle, collimator angle, jaw positions, and dose output. The study encompassed the clinical treatment plans of 29 brain cancer patients who received stereotactic ablative radiotherapy (SABR). Six error magnitudes were investigated per source of error. In addition, the Eclipse AAA beam model dosimetric leaf gap (DLG) parameter was varied with four error magnitudes. Error detectability was determined based on the area under the receiver operating characteristic (ROC) curve (AUC). Detectability of DLG errors was good or excellent (AUC >0.8) at an error magnitude of at least ±0.4 mm, while MLC leaf position and gantry angle errors reached good or excellent detectability at error magnitudes of at least 1.0 mm and 0.6°, respectively. Ideal thresholds, that is, gamma passing rates, to maximize sensitivity and specificity ranged from 79.1% to 98.7%. The detectability of collimator angle, jaw position, and dose output errors was poor for all investigated error magnitudes, with an AUC between 0.5 and 0.6. The ArcCHECK device's ability to detect errors from treatment machine-related sources was evaluated, and ideal gamma passing rate thresholds were determined for each source of error. The ArcCHECK was able to detect errors in DLG value, MLC leaf positions, and gantry angle. The ArcCHECK was unable to detect the studied errors in collimator angle, jaw positions, and dose output.

Radiation protection modelling for 2.5 Petawatt-laser production of ultrashort x-ray, proton and ion bunches: Monte Carlo model of the Munich CALA facility.

The 'Centre for Advanced Laser Applications' (CALA) is a new research institute for laser-based acceleration of electron beams for brilliant x-ray generation, laser-driven sub-nanosecond bunches of protons and heavy ions for biomedical applications like imaging and tumour therapy as well as for nuclear and high-field physics.The radiation sources emerging from experiments using the up to 2.5 petawatt laser pulses with 25 femtosecond duration will be mixed particle-species of high intensity, high energy and pulsed, thus posing new challenges compared to conventional radiation protection. Such worldwide pioneering laser experiments result in source characteristics that require careful a-priori radiation safety simulations.The FLUKA Monte-Carlo code was used to model the five CALA experimental caves, including the corridors, halls and air spaces surrounding the caves. Beams of electrons (<5 GeV), protons (<200 MeV),12C (<400MeV/u) and197Au (<10MeV/u) ions were simulated using spectra, divergences and bunch-charges based on expectations from recent scientific progress.Simulated dose rates locally can exceed 1.5 kSv h-1inside beam dumps. Vacuum pipes in the cave walls for laser transport and extraction channels for the generated x-rays result in small dose leakage to neighboring areas. Secondary neutrons contribute to most of the prompt dose rate outside caves into which the beam is delivered. This secondary radiation component causes non-negligible dose rates to occur behind walls to which large fluences of secondary particles are directed.By employing adequate beam dumps matched to beam-divergence, magnets, passive shielding and laser pulse repetition limits, average dose rates in- and outside the experimental building stay below design specifications (<0.5µSv h-1) for unclassified areas,<2.5µSv h-1for supervised areas,<7.5µSv h-1maximum local dose rate) and regulatory limits (<1mSv a-1for unclassified areas).

Evaluation of proton and photon dose distributions recalculated on 2D and 3D Unet-generated pseudoCTs from T1-weighted MR head scans.

Introduction: The recent developments of magnetic resonance (MR) based adaptive strategies for photon and, potentially for proton therapy, require a fast and reliable conversion of MR images to X-ray computed tomography (CT) values. CT values are needed for photon and proton dose calculation. The improvement of conversion results employing a 3D deep learning approach is evaluated. Material and methods: A database of 89 T1-weighted MR head scans with about 100 slices each, including rigidly registered CTs, was created. Twenty-eight validation patients were randomly sampled, and four patients were selected for application. The remaining patients were used to train a 2D and a 3D U-shaped convolutional neural network (Unet). A stack size of 32 slices was used for 3D training. For all application cases, volumetric modulated arc therapy photon and single-field uniform dose pencil-beam scanning proton plans at four different gantry angles were optimized for a generic target on the CT and recalculated on 2D and 3D Unet-based pseudoCTs. Mean (absolute) error (MAE/ME) and a gradient sharpness estimate were used to quantify the image quality. Three-dimensional gamma and dose difference analyses were performed for photon (gamma criteria: 1%, 1 mm) and proton dose distributions (gamma criteria: 2%, 2 mm). Range (80% fall off) differences for beam's eye view profiles were evaluated for protons. Results: Training 36 h for 1000 epochs in 3D (6 h for 200 epochs in 2D) yielded a maximum MAE of 147 HU (135 HU) for the application patients. Except for one patient gamma pass rates for photon and proton dose distributions were above 96% for both Unets. Slice discontinuities were reduced for 3D training at the cost of sharpness. Conclusions: Image analysis revealed a slight advantage of 2D Unets compared to 3D Unets. Similar dose calculation performance was reached for the 2D and 3D network.

Towards a novel small animal proton irradiation platform: the SIRMIO project.

Background: Precision small animal radiotherapy research is a young emerging field aiming to provide new experimental insights into tumor and normal tissue models in different microenvironments, to unravel complex mechanisms of radiation damage in target and non-target tissues and assess efficacy of novel therapeutic strategies. For photon therapy, modern small animal radiotherapy research platforms have been developed over the last years and are meanwhile commercially available. Conversely, for proton therapy, which holds potential for an even superior outcome than photon therapy, no commercial system exists yet. Material and methods: The project SIRMIO (Small Animal Proton Irradiator for Research in Molecular Image-guided Radiation-Oncology) aims at realizing and demonstrating an innovative portable prototype system for precision image-guided small animal proton irradiation, suitable for installation at existing clinical treatment facilities. The proposed design combines precise dose application with in situ multi-modal anatomical image guidance and in vivo verification of the actual treatment delivery. Results and conclusions: This manuscript describes the status of the different components under development, featuring a dedicated beamline for degradation and focusing of clinical proton beams, along with novel detector systems for in situimaging and range verification. The foreseen workflow includes pre-treatment proton transmission imaging, complemented by ultrasonic tumor localization, for treatment planning and position verification, followed by image-guided delivery with on site range verification by means of ionoacoustics (for pulsed beams) and positron-emission-tomography (PET, for continuous beams). The proposed compact and cost-effective system promises to open a new era in small animal proton therapy research, contributing to the basic understanding of in vivo radiation action to identify areas of potential breakthroughs for future translation into innovative clinical strategies.

Sensitivity of post treatment positron emission tomography/computed tomography to detect inter-fractional range variations in scanned ion beam therapy.

BACKGROUND: Ion therapy, especially with modern scanning beam delivery, offers very sharp dose gradients for highly conformal cancer treatment. However, it is very sensitive to uncertainties of tissue stopping properties as well as to anatomical changes and setup errors, making range verification highly desirable. To this end, positron emission tomography (PET) can be used to measure decay products of ß+-emitters created in interactions inside the patient. This work investigates the sensitivity of post treatment PET/CT (computed tomography) to detect inter-fractional range variations. MATERIAL AND METHODS: Fourteen patients of different indication underwent PET/CT monitoring after selected treatment fractions with scanned proton or carbon ion beams. In addition to PET/CT measurements, PET and dose distributions were simulated on different co-registered CT data. Pairs of PET data were then analyzed in terms of longitudinal shifts along the beam path, as surrogate of inter-fractional range deviations. These findings were compared to changes of dose-volume-histogram indexes and corresponding dose as well as CT shifts to disentangle the origin of possible PET shifts. RESULTS: Biological washout modeling (PET simulations) and low (<55 Bq/ml) activity concentrations (offline PET measurements, especially for 12C ions) were the main limitations for clinical treatment verification. For two selected cases, the benefit of improved washout modeling based on organ segmentation could be demonstrated. Overall, inter-fractional range shifts up to ±3 mm could be deduced from both PET measurements and simulations, and found well correlated (typically within 1.8 mm) to anatomical changes derived from CT scans, in agreement with dose data. CONCLUSIONS: Despite known limitations of post treatment PET/CT imaging, this work indicates its potential for assessing inter-fractional changes and points to future developments for improved PET-based treatment verification.

Practical implications for the quality assurance of modulated radiation therapy techniques using point detector arrays.

PURPOSE: Linac parameters potentially influencing the delivery quality of IMRT and VMAT plans are investigated with respect to threshold ranges, consequently to be considered in a linac based quality assurance procedure. Three commercially available 2D arrays are used to further investigate the influence of the measurement device. METHODS: Using three commercially available 2D arrays (Mx: MatriXXevolution , Oc: Octavius1500 , Mc: MapCHECK2), simple static measurements, measurements for MLC characterization and dynamic interplay of gantry movement, MLC movement and variable dose rate were performed. The results were evaluated with respect to each single array as well as among each other. RESULTS: Simple static measurements showed different array responses to dose, dose rate and profile homogeneity and revealed instabilities in dose delivery and profile shape during linac ramp up. Using the sweeping gap test, all arrays were able to detect small leaf misalignments down to ±0.1 mm, but this test also demonstrated up to 15% dose deviation due to profile instabilities and fast accelerating leaves during linac ramp up. Tests including gantry rotation showed different stability of gantry mounts for each array. Including gantry movement and dose rate variability, differences compared to static delivery were smaller compared to dose differences when simultaneously controling interplay of gantry movement, leaf movement and dose rate variability. CONCLUSION: Linac based QA is feasible with the tested commercially available 2D arrays. Limitations of each array and the linac ramp up characteristics should be carefully considered during individual plan generation and regularly checked in linac QA. Especially the dose and dose profile during linac ramp up should be checked regularly, as well as MLC positioning accuracy using a sweeping gap test. Additionally, dynamic interplay tests including various gantry rotation speeds and angles, various leaf speeds and various dose rates should be included.

Linearity of patient positioning detection : a phantom study of skin markers, cone beam computed tomography, and 3D ultrasound.

BACKGROUND: Three-dimensional ultrasound (3D-US) is a modality complementary to kilovoltage cone beam computed tomography (kV-CBCT) and skin markers for patient positioning detection. This study compares the linearity of evaluations based on measurements using a modern 3D-US system (Elekta Clarity®; Elekta, Stockholm, Sweden), a kV-CBCT system (Elekta iView®), and skin markers. MATERIALS AND METHODS: An investigator deliberately displaced a multimodal phantom by up to ± 30 mm along different axes. The following data points were acquired: 27 along the lateral axis, 29 along the longitudinal axis, 27 along the vertical axis, and 27 along the space diagonal. At each of these 110 positions, the displacements according to skin' markers were recorded and scans were performed using both 3D-US and kV-CBCT. Shifts were detected by matching bony anatomy or soft tissue density to a reference planning CT in the case of kV-CBCT and for 3D-US, by matching ultrasound volume data to a reference planning volume. A consensus value was calculated from the average of the four modalities. With respect to this consensus value, the linearity (offset and regression coefficient, i.e., slope), average offset, systematic error, and random error of all four modalities were calculated for each axis. RESULTS: Linearity was similar for all four modalities, with regression coefficients between 0.994 and 1.012, and all offsets below 1 mm. The systematic errors of skin markers and 3D-US were higher than for kV-CBCT, but random errors were similar. In particular, 3D-US demonstrated an average offset of 0.36 mm to the right, 0.08 mm inferiorly, and 0.15 mm anteriorly; the systematic error was 0.36 mm laterally, 0.35 mm longitudinally, and 0.22 mm vertically; the random error was 0.15 mm laterally, 0.30 mm longitudinally, and 0.12 mm vertically. A total of 109 out of 110 (99 %) 3D-US measurements were within 1 mm of the consensus value on either axis. CONCLUSION: The linearity of 3D-US was no worse than that of skin markers or kV-CBCT. Average offsets, systematic errors, and random errors were all below 1 mm. Optimal margins in the order of 1 mm could be achieved in the controlled laboratory setting of this phantom study.

Comparing cone-beam CT intensity correction methods for dose recalculation in adaptive intensity-modulated photon and proton therapy for head and neck cancer.

BACKGROUND: Adaptive intensity-modulated photon and proton radiotherapy (IMRT and IMPT) of head and neck (H&N) cancer requires frequent three-dimensional (3D) dose calculation. We compared two approaches for dose recalculation on the basis of intensity-corrected cone-beam (CB) x-ray computed tomography (CT) images. MATERIAL AND METHODS: For nine H&N tumor patients, virtual CTs (vCT) were generated by deformable image registration of the planning CT (pCT) to the CBCT. The second intensity correction approach used population-based lookup tables for scaling CBCT intensities to the pCT HU range (CBCTLUT). IMRT and IMPT plans were generated with a commercial treatment planning system. Dose recalculations on vCT and CBCTLUT were analyzed using a (3%, 3 mm) gamma-index analysis and comparison of normal tissue and tumor dose/volume parameters. A replanning CT (rpCT) acquired within three days of the CBCT served as reference. Single field uniform dose (SFUD) proton plans were created and recalculated on vCT and CBCTLUT for proton range comparison. RESULTS: Dose/volume parameters showed minor differences between rpCT, vCT and CBCTLUT in IMRT, but clinically relevant deviations between CBCTLUT and rpCT in the spinal cord for IMPT. Gamma-index pass-rates were found increased for vCT with respect to CBCTLUT in IMPT (by up to 21 percentage points) and IMRT (by up to 9 percentage points) for most cases. The SFUD-based proton range assessment showed improved agreement of vCT and rpCT, with 88-99% of the depth dose profiles in beam's eye view agreeing within 3 mm. For CBCTLUT, only 80-94% of the profiles fulfilled this criterion. CONCLUSION: vCT and CBCTLUT are suitable options for dose recalculation in adaptive IMRT. In the scope of IMPT, the vCT approach is preferable.

A quantitative assessment of Geant4 for predicting the yield and distribution of positron-emitting fragments in ion beam therapy.

Objective.To compare the accuracy with which different hadronic inelastic physics models across ten Geant4 Monte Carlo simulation toolkit versions can predict positron-emitting fragments produced along the beam path during carbon and oxygen ion therapy.Approach.Phantoms of polyethylene, gelatin, or poly(methyl methacrylate) were irradiated with monoenergetic carbon and oxygen ion beams. Post-irradiation, 4D PET images were acquired and parent11C,10C and15O radionuclides contributions in each voxel were determined from the extracted time activity curves. Next, the experimental configurations were simulated in Geant4 Monte Carlo versions 10.0 to 11.1, with three different fragmentation models-binary ion cascade (BIC), quantum molecular dynamics (QMD) and the Liege intranuclear cascade (INCL++) - 30 model-version combinations. Total positron annihilation and parent isotope production yields predicted by each simulation were compared between simulations and experiments using normalised mean squared error and Pearson cross-correlation coefficient. Finally, we compared the depth of the maximum positron annihilation yield and the distal point at which the positron yield decreases to 50% of peak between each model and the experimental results.Main results.Performance varied considerably across versions and models, with no one version/model combination providing the best prediction of all positron-emitting fragments in all evaluated target materials and irradiation conditions. BIC in Geant4 10.2 provided the best overall agreement with experimental results in the largest number of test cases. QMD consistently provided the best estimates of both the depth of peak positron yield (10.4 and 10.6) and the distal 50%-of-peak point (10.2), while BIC also performed well and INCL generally performed the worst across most Geant4 versions.Significance.The best predictions of the spatial distribution of positron annihilations and positron-emitting fragment production along the beam path during carbon and oxygen ion therapy was obtained using Geant4 10.2.p03 with BIC or QMD. These version/model combinations are recommended for future heavy ion therapy research.

On the potential biological impact of radiation-induced acoustic emissions during ultra-high dose rate electron radiotherapy: a preliminary study.

Ionizing radiation pulses delivered at ultra-high dose rates in emerging FLASH radiotherapy can result in high-intensity low-frequency thermoacoustic emissions that may have a biological impact. This study aims at providing insights into the thermoacoustic emissions expected during FLASH radiotherapy and their likelihood of inducing acoustic cavitation. The characteristics of acoustic waves induced by the energy deposition of a pulsed electron beam similar to previous pre-clinical FLASH radiotherapy studies and their propagation in murine head-like phantoms are investigated in-silico. The results show that the generated pressures are sufficient to produce acoustic cavitation due to resonance in the irradiated object. It suggests that thermoacoustics may, in some irradiation scenarios, contribute to the widely misunderstood FLASH effect or cause adverse effects if not taken into account at the treatment planning stage.

Experience and new prospects of PET imaging for ion beam therapy monitoring.

Pioneering investigations on the usage of positron-emission-tomography (PET) for the monitoring of ion beam therapy with light (protons, helium) and heavier (stable and radioactive neon, carbon and oxygen) ions started shortly after the first realization of planar and tomographic imaging systems, which were able to visualize the annihilation of positrons resulting from irradiation induced or implanted positron emitting nuclei. And while the first clinical experience was challenged by the utilization of instrumentation directly adapted from nuclear medicine applications, new detectors optimized for this unconventional application of PET imaging are currently entering the phase of (pre)clinical testing for more reliable monitoring of treatment delivery during irradiation. Moreover, recent advances in detector technologies and beam production open several new exciting opportunities which will not only improve the performance of PET imaging under the challenging conditions of in-beam applications in ion beam therapy, but will also likely expand its field of application. In particular, the combination of PET and Compton imaging can enable the most efficient utilization of all possible radiative emissions for both stable and radioactive ion beams, while positronium lifetime imaging may enable probing new features of the underlying tumour and normal tissue environment. Thereby, PET imaging will not only provide means for volumetric reconstruction of the delivered treatment and in-vivo verification of the beam range, but can also shed new insights for biological optimization of the treatment or treatment response assessment.

Translation of AI into oncology clinical practice.

Artificial intelligence (AI) is a transformative technology that is capturing popular imagination and can revolutionize biomedicine. AI and machine learning (ML) algorithms have the potential to break through existing barriers in oncology research and practice such as automating workflow processes, personalizing care, and reducing healthcare disparities. Emerging applications of AI/ML in the literature include screening and early detection of cancer, disease diagnosis, response prediction, prognosis, and accelerated drug discovery. Despite this excitement, only few AI/ML models have been properly validated and fewer have become regulated products for routine clinical use. In this review, we highlight the main challenges impeding AI/ML clinical translation. We present different clinical use cases from the domains of radiology, radiation oncology, immunotherapy, and drug discovery in oncology. We dissect the unique challenges and opportunities associated with each of these cases. Finally, we summarize the general requirements for successful AI/ML implementation in the clinic, highlighting specific examples and points of emphasis including the importance of multidisciplinary collaboration of stakeholders, role of domain experts in AI augmentation, transparency of AI/ML models, and the establishment of a comprehensive quality assurance program to mitigate risks of training bias and data drifts, all culminating toward safer and more beneficial AI/ML applications in oncology labs and clinics.

Comparison of reconstructed prompt gamma emissions using maximum likelihood estimation and origin ensemble algorithms for a Compton camera system tailored to proton range monitoring.

Compton-based prompt gamma (PG) imaging is being investigated by several groups as a potential solution for in vivo range monitoring in proton therapy. The performance of this technique depends on the detector system as well as the ability of the reconstruction method to obtain good spatial resolution to establish a quantitative correlation between the PG emission and the proton beam range in the patient. To evaluate the feasibility of PG imaging for range monitoring, we quantitatively evaluated the emission distributions reconstructed by a Maximum Likelihood Expectation Maximization (MLEM) and a Stochastic Origin Ensemble (SOE) algorithm. To this end, we exploit experimental and Monte Carlo (MC) simulation data acquired with the Polaris-J Compton Camera (CC) prototype. The differences between the proton beam range (RD) defined as the 80% distal dose fall-off and the PG range (RPG), obtained by fitting the distal end of the reconstructed profile with a sigmoid function, were quantified. A comparable performance of both reconstruction algorithms was found. For both experimental and simulated irradiation scenarios, the correlation between RD and RPG was within 5 mm. These values were consistent with the ground truth distance (RD-RPGg≈ 3 mm) calculated by using the expected PG emission available from MC simulation. Furthermore, shifts of 3 mm in the proton beam range were resolved with the MLEM algorithm by calculating the relative difference between the RPG for each reconstructed profile. In non-homogeneous targets, the spatial changes in the PG emission due to the different materials could not be fully resolved from the reconstructed profiles; however, the fall-off region still resembled the ground truth emission. For this scenario, the PG correlation (RD-RPG) varied from 0.1 mm to 4 mm, which is close to the ground truth correlation (3 mm). This work provides a framework for the evaluation of the range monitoring capabilities of a CC device for PG imaging. The two investigated image reconstruction algorithms showed a comparable and consistent performance for homogeneous and heterogeneous targets.

Scatter correction of 4D cone beam computed tomography to detect dosimetric effects due to anatomical changes in proton therapy for lung cancer.

BACKGROUND: The treatment of moving tumor entities is expected to have superior clinical outcomes, using image-guided adaptive intensity-modulated proton therapy (IMPT). PURPOSE: For 21 lung cancer patients, IMPT dose calculations were performed on scatter-corrected 4D cone beam CTs (4DCBCTcor ) to evaluate their potential for triggering treatment adaptation. Additional dose calculations were performed on corresponding planning 4DCTs and day-of-treatment 4D virtual CTs (4DvCTs). METHODS: A 4DCBCT correction workflow, previously validated on a phantom, generates 4DvCT (CT-to-CBCT deformable registration) and 4DCBCTcor images (projection-based correction using 4DvCT as a prior) with 10 phase bins, using day-of-treatment free-breathing CBCT projections and planning 4DCT images as input. Using a research planning system, robust IMPT plans administering eight fractions of 7.5 Gy were created on a free-breathing planning CT (pCT) contoured by a physician. The internal target volume (ITV) was overridden with muscle tissue. Robustness settings for range and setup uncertainties were 3% and 6 mm, and a Monte Carlo dose engine was used. On every phase of planning 4DCT, day-of-treatment 4DvCT, and 4DCBCTcor , the dose was recalculated. For evaluation, image analysis as well as dose analysis were performed using mean error (ME) and mean absolute error (MAE) analysis, dose-volume histogram (DVH) parameters, and 2%/2-mm gamma pass rate analysis. Action levels (1.6% ITV D98 and 90% gamma pass rate) based on our previous phantom validation study were set to determine which patients had a loss of dosimetric coverage. RESULTS: Quality enhancements of 4DvCT and 4DCBCTcor over 4DCBCT were observed. ITV D98% and bronchi D2% had its largest agreement for 4DCBCTcor -4DvCT, and the largest gamma pass rates (>94%, median 98%) were found for 4DCBCTcor -4DvCT. Deviations were larger and gamma pass rates were smaller for 4DvCT-4DCT and 4DCBCTcor -4DCT. For five patients, deviations were larger than the action levels, suggesting substantial anatomical changes between pCT and CBCT projections acquisition. CONCLUSIONS: This retrospective study shows the feasibility of daily proton dose calculation on 4DCBCTcor for lung tumor patients. The applied method is of clinical interest as it generates up-to-date in-room images, accounting for breathing motion and anatomical changes. This information could be used to trigger replanning.

Impact of secondary particles on the magnetic field generated by a proton pencil beam: a finite-element analysis based on Geant4-DNA simulations.

PURPOSE: To investigate the static magnetic field generated by a proton pencil beam as a candidate for range verification by means of Monte Carlo simulations, thereby improving upon existing analytical calculations. We focus on the impact of statistical current fluctuations and secondary protons and electrons. METHODS: We considered a pulsed beam (10 µ ${\umu}$ s pulse duration) during the duty cycle with a peak beam current of 0.2 µ $\umu$ A and an initial energy of 100 MeV. We ran Geant4-DNA Monte Carlo simulations of a proton pencil beam in water and extracted independent particle phase spaces. We calculated longitudinal and radial current density of protons and electrons, serving as an input for a magnetic field estimation based on a finite element analysis in a cylindrical geometry. We made sure to allow for non-solenoidal current densities as is the case of a stopping proton beam. RESULTS: The rising proton charge density toward the range is not perturbed by energy straggling and only lowered through nuclear reactions by up to 15%, leading to an approximately constant longitudinal current. Their relative low density however (at most 0.37 protons/mm3 for the 0.2  µ ${\umu}$ A current and a beam cross-section of 2.5 mm), gives rise to considerable current density fluctuations. The radial proton current resulting from lateral scattering and being two orders of magnitude weaker than the longitudinal current is subject to even stronger fluctuations. Secondary electrons with energies above 10 eV, that far outnumber the primary protons, reduce the primary proton current by only 10% due to their largely isotropic flow. A small fraction of electrons (<1%), undergoing head-on collisions, constitutes the relevant electron current. In the far-field, both contributions to the magnetic field strength (longitudinal and lateral) are independent of the beam spot size. We also find that the nuclear reaction-related losses cause a shift of 1.3 mm to the magnetic field profile relative to the actual range, which is further enlarged to 2.4 mm by the electron current (at a distance of ρ = 50 $\rho =50$  mm away from the central beam axis). For ρ > 45 $\rho >45$  mm, the shift increases linearly. While the current density variations cause significant magnetic field uncertainty close to the central beam axis with a relative standard deviation (RSD) close to 100%, they average out at a distance of 10 cm, where the RSD of the total magnetic field drops below 2%. CONCLUSIONS: With the small influence of the secondary electrons together with the low RSD, our analysis encourages an experimental detection of the magnetic field through sensitive instrumentation, such as optical magnetometry or SQUIDs.

Range verification of a clinical proton beam in an abdominal phantom by co-registration of ionoacoustics and ultrasound.

Objective.The range uncertainty in proton radiotherapy is a limiting factor to achieve optimum dose conformity to the tumour volume. Ionoacoustics is a promising approach forin siturange verification, which would allow to reduce the size of the irradiated volume relative to the tumour volume. The energy deposition of a pulsed proton beam leads to an acoustic pressure wave (ionoacoustics), the detection of which allows conclusion about the distance between the Bragg peak and the acoustic detector. This information can be transferred into a co-registered ultrasound image, marking the Bragg peak position relative to the surrounding anatomy.Approach.A CIRS 3D abdominal phantom was irradiated with 126 MeV protons at a clinical proton therapy centre. Acoustic signals were recorded on the beam axis distal to the Bragg peak with a Cetacean C305X hydrophone. The ionoacoustic measurements were processed with a correlation filter using simulated filter templates. The hydrophone was rigidly attached to an ultrasound device (Interson GP-C01) recording ultrasound images of the irradiated region.Main results.The time of flight obtained from ionoacoustic measurements were transferred to an ultrasound image by means of an optoacoustic calibration measurement. The Bragg peak position was marked in the ultrasound image with a statistical uncertainty ofσ= 0.5 mm of 24 individual measurements depositing 1.2 Gy at the Bragg peak. The difference between the evaluated Bragg peak position and the one obtained from irradiation planning (1.0 mm) is smaller than the typical range uncertainty (≈4 mm) at the given penetration depth (10 cm).Significance.The measurements show that it is possible to determine the Bragg peak position of a clinical proton beam with submillimetre precision and transfer the information to an ultrasound image of the irradiated region. The dose required for this is smaller than that used for a typical irradiation fraction.

Multi-stage image registration based on list-mode proton radiographies for small animal proton irradiation: A simulation study.

We present a multi-stage and multi-resolution deformable image registration framework for image-guidance at a small animal proton irradiation platform. The framework is based on list-mode proton radiographies acquired at different angles, which are used to deform a 3D treatment planning CT relying on normalized mutual information (NMI) or root mean square error (RMSE) in the projection domain. We utilized a mouse X-ray micro-CT expressed in relative stopping power (RSP), and obtained Monte Carlo simulations of proton images in list-mode for three different treatment sites (brain, head and neck, lung). Rigid transformations and controlled artificial deformation were applied to mimic position misalignments, weight loss and breathing changes. Results were evaluated based on the residual RMSE of RSP in the image domain including the comparison of extracted local features, i.e. between the reference micro-CT and the one transformed taking into account the calculated deformation. The residual RMSE of the RSP showed that the accuracy of the registration framework is promising for compensating rigid (>97% accuracy) and non-rigid (∼95% accuracy) transformations with respect to a conventional 3D-3D registration. Results showed that the registration accuracy is degraded when considering the realistic detector performance and NMI as a metric, whereas the RMSE in projection domain is rather insensitive. This work demonstrates the pre-clinical feasibility of the registration framework on different treatment sites and its use for small animal imaging with a realistic detector. Further computational optimization of the framework is required to enable the use of this tool for online estimation of the deformation.

On the robustness of multilateration of ionoacoustic signals for localization of the Bragg peak at pre-clinical proton beam energies in water.

Objectives.The energy deposited in a medium by a pulsed proton beam results in the emission of thermoacoustic waves, also called ionoacoustics (IA). The proton beam stopping position (Bragg peak) can be retrieved from a time-of-flight analysis (ToF) of IA signals acquired at different sensor locations (multilateration). This work aimed to assess the robustness of multilateration methods in proton beams at pre-clinical energies for the development of a small animal irradiator.Approach.The accuracy of multilateration performed using different algorithms; namely, time of arrival and time difference of arrival, was investigatedin-silicofor ideal point sources in the presence of realistic uncertainties on the ToF estimation and ionoacoustic signals generated by a 20 MeV pulsed proton beam stopped in a homogeneous water phantom. The localisation accuracy was further investigated experimentally based on two different measurements with pulsed monoenergetic proton beams at energies of 20 and 22 MeV.Main results.It was found that the localisation accuracy mainly depends on the position of the acoustic detectors relative to the proton beam due to spatial variation of the error on the ToF estimation. By optimally positioning the sensors to reduce the ToF error, the Bragg peak could be locatedin-silicowith an accuracy better than 90µm (2% error). Localisation errors going up to 1 mm were observed experimentally due to inaccurate knowledge of the sensor positions and noisy ionoacoustic signals.Significance.This study gives a first overview of the implementation of different multilateration methods for ionoacoustics-based Bragg peak localisation in two- and three-dimensions at pre-clinical energies. Different sources of uncertainty were investigated, and their impact on the localisation accuracy was quantifiedin-silicoand experimentally.

Design study of a novel geometrical arrangement for an in-beam small animal positron emission tomography scanner.

Objective.We designed a geometrical solution for a small animal in-beam positron emission tomography (PET) scanner to be used in the project SIRMIO (Small animal proton irradiator for research in molecular image-guided radiation-oncology). The system is based on 56 scintillator blocks of pixelated LYSO crystals. The crystals are arranged providing a pyramidal-step shape to optimize the geometrical coverage in a spherical configuration.Approach.Different arrangements have been simulated and compared in terms of spatial resolution and sensitivity. The chosen setup enables us to reach a good trade-off between a solid angle coverage and sufficient available space for the integration of additional components of the first design prototype of the SIRMIO platform. The possibility of moving the mouse holder inside the PET scanner furthermore allows for achieving the optimum placement of the irradiation area for all the possible tumor positions in the body of the mouse. The work also includes a study of the scintillator material where LYSO and GAGG are compared with a focus on the random coincidence noise due to the natural radioactivity of Lutetium in LYSO, justifying the choice of LYSO for the development of the final system.Main results.The best imaging performance can be achieved with a sub-millimeter spatial resolution and sensitivity of 10% in the center of the scanner, as verified in thorough simulations of point sources. The simulation of realistic irradiation scenarios of proton beams in PMMA targets with/without air gaps indicates the ability of the proposed PET system to detect range shifts down to 0.2 mm.Significance.The presented results support the choice of the identified optimal design for a novel spherical in-beam PET scanner which is currently under commissioning for application to small animal proton and light ion irradiation, and which might find also application, e.g. for biological image-guidance in x-ray irradiation.

Towards real-time EPID-based 3D in vivo dosimetry for IMRT with Deep Neural Networks: A feasibility study.

We investigate the potential of the Deep Dose Estimate (DDE) neural network to predict 3D dose distributions inside patients with Monte Carlo (MC) accuracy, based on transmitted EPID signals and patient CTs. The network was trained using as input patient CTs and first-order dose approximations (FOD). Accurate dose distributions (ADD) simulated with MC were given as training targets. 83 pelvic CTs were used to simulate ADDs and respective EPID signals for subfields of prostate IMRT plans (gantry at 0∘). FODs were produced as backprojections from the EPID signals. 581 ADD-FOD sets were produced and divided into training and test sets. An additional dataset simulated with gantry at 90∘ (lateral set) was used for evaluating the performance of the DDE at different beam directions. The quality of the FODs and DDE-predicted dose distributions (DDEP) with respect to ADDs, from the test and lateral sets, was evaluated with gamma analysis (3%,2 mm). The passing rates between FODs and ADDs were as low as 46%, while for DDEPs the passing rates were above 97% for the test set. Meaningful improvements were also observed for the lateral set. The high passing rates for DDEPs indicate that the DDE is able to convert FODs into ADDs. Moreover, the trained DDE predicts the dose inside a patient CT within 0.6 s/subfield (GPU), in contrast to 14 h needed for MC (CPU-cluster). 3D in vivo dose distributions due to clinical patient irradiation can be obtained within seconds, with MC-like accuracy, potentially paving the way towards real-time EPID-based in vivo dosimetry.