Fall-off Region Research Articles

Electron beam response corrections for an ultra-high-dose-rate capable diode dosimeter.

Ultra-high-dose-rate (UHDR) electron beams have been commonly utilized in FLASH studies and the translation of FLASH Radiotherapy (RT) to the clinic. The EDGE diode detector has potential use for UHDR dosimetry albeit with a beam energy dependency observed. The purpose is to present the electron beam response for an EDGE detector in dependence on beam energy, to characterize the EDGE detector's response under UHDR conditions, and to validate correction factors derived from the first detailed Monte Carlo model of the EDGE diode against measurements, particularly under UHDR conditions. Percentage depth doses (PDDs) for the UHDR Mobetron were measured with both EDGE detectors and films. A detailed Monte Carlo (MC) model of the EDGE detector has been configured according to the blueprint provided by the manufacturer under an NDA agreement. Water/silicon dose ratios of EDGE detector for a series of mono-energetic electron beams have been calculated. The dependence of the water/silicon dose ratio on depth for a FLASH relevant electron beam was also studied. An analytical approach for the correction of PDD measured with EDGE detectors was established. Water/silicon dose ratio decreased with decreasing electron beam energy. For the Mobetron 9 MeV UHDR electron beam, the ratio decreased from 1.09 to 1.03 in the build-up region, maintained in range of 0.98-1.02 at the fall-off region and raised to a plateau in value of 1.08 at the tail. By applying the corrections, good agreement between the PDDs measured by the EDGE detector and those measured with film was achieved. Electron beam response of an UHDR capable EDGE detector was derived from first principles utilizing a sophisticated MC model. An analytical approach was validated for the PDDs of UHDR electron beams. The results demonstrated the capability of EDGE detector in measuring PDDs of UHDR electron beams.

Calculation of biological effectiveness of SOBP proton beams: a TOPAS Monte Carlo study

Objective. This study aims to investigate the biological effectiveness of Spread-Out Bragg-Peak (SOBP) proton beams with initial kinetic energies 50–250 MeV at different depths in water using TOPAS Monte Carlo code. Approach. The study modelled SOBP proton beams using TOPAS time feature. Various LET-based models and Repair-Misrepair-Fixation model were employed to calculate Relative Biological Effectiveness (RBE) for V79 cell lines at different on-axis depths based on TOPAS. Microdosimetric Kinetic Model and biological weighting function-based models, which utilize microdosimetric distributions, were also used to estimate the RBE. A phase-space-based method was adopted for calculating microdosimetric distributions. Main results. The trend of variation of RBE with depth is similar in all the RBE models, but the absolute RBE values vary based on the calculation models. RBE sharply increases at the distal edge of SOBP proton beams. In the entrance region of all the proton beams, RBE values at 4 Gy i.e. RBE(4 Gy) resulting from different models are in the range of 1.04–1.07, comparable to clinically used generic RBE of 1.1. Moving from the proximal to distal end of the SOBP, RBE(4 Gy) is in the range of 1.15–1.33, 1.13–1.21, 1.11–1.17, 1.13–1.18 and 1.17–1.21, respectively for 50, 100, 150, 200 and 250 MeV SOBP beams, whereas at the distal dose fall-off region, these values are 1.68, 1.53, 1.44, 1.42 and 1.40, respectively. Significance. The study emphasises application of depth-, dose- and energy- dependent RBE values in clinical application of proton beams.

Multicentric characterisation of lateral beam profiles generated by 6FFF beam of three 0.35 T MR-linac systems

BackgroundThe physical characterisation of FFF-beam profiles in the presence of a magnetic field requires a new standardization procedure and formulation. PurposeThe aim of this multicentric experience is to propose new normalisation factors to allow for the calculation of standard parameters typical of flattened beams, such as dosimetric field size and penumbra, for a 6 MV FFF beam from an MR-linac. MethodsThe measurements were carried out on three ViewRay-MRIdiansystems. An equal set of measurements was acquired using the same equipment. Transverse beam profiles were acquired at seven different depthsand for six beam dimensions.The inflection point was estimated as the position of the maximum of a Gaussian fit obtained from the first derivative of the profiles. The position of the minimum and maximum points of the second derivative of the above Gaussian described the fall-off region, and the external peaks of the third derivative were considered as the in-field and out-field points. The profile normalisation was determined by imposing a 55% dose level at the inflection point and the renormalisation factors were calculated. ResultsThe position of the inflection point, and the second and third derivatives peaks were analysed,and the renormalisation factors as a function of field size and depth were determined. The values of the unflatness and the slope have been calculated for different depths and field sizes. ConclusionThis study represents the first multi-centric evaluation of the profiles on different low-field MR-Linac systems and theset of renormalisation parameters to analyse the FFF-beam on that system was effectively proposed.

Sensitivity of a mini-TEPC to radiation quality variations in clinical proton beams

PurposeThis work aims at studying the sensitivity of a miniaturized Tissue-Equivalent Proportional Counter to variations of beam quality in clinical radiation fields, to further investigate its performances as radiation quality monitor. MethodsMeasurements were taken at the CATANA facility (INFN-LNS, Catania, Italy), in a monoenergetic and an energy-modulated proton beam with the same initial energy of 62 MeV. PMMA layers were placed in front of the detector to measure at different depths along the depth-dose profile. The frequency- and dose-mean lineal energy were compared to the track- and dose-averaged LET calculated by Monte Carlo simulations. A microdosimetric evaluation of the Relative Biological Effectiveness (RBE) was performed and compared with cell survival experiments. ResultsMicrodosimetric distributions measured at identical depths in the two beams show spectral differences reflecting their different radiation quality. Discrepancies are most evident at depths corresponding to the Spread-Out Bragg Peak, while spectra at the entrance and in the dose fall-off regions are similar. This can be explained by the different energy components that compose the pristine and spread-out peaks at each depth. The trend of microdosimetric mean values matches that of calculated LET averages along the entire penetration depth, and the microdosimetric estimation of RBE is consistent with radiobiological data not only at 2 Gy but also at lower dose levels, such as those absorbed by healthy tissues. ConclusionsThe mini-TEPC is sensitive to differences in radiation quality resulting from different modulations of the proton beam, confirming its potential for beam quality monitoring in proton therapy.

Unveiling the oxygen addition kinetics of 1-Hydroxy-2-Hexyl and 2-Hydroxy-1-Hexyl at low temperature

Hydroxyalkyl radicals formed from addition of OH to alkenes play a crucial role in the oxidation of alkenes. The oxygen addition kinetics of 1‑hydroxy-2-hexyl and 2‑hydroxy-1-hexyl at low temperature has been investigated by using quantum chemical calculations coupled with the RRKM/master equation simulations. The main reaction pathways that occur on the potential energy surfaces of HOC6H12 + O2 were studied at the CCSD(T)/CBS(cc-pV[T,Q]Z)//M06–2X/jun-cc-pVTZ level. Rate coefficients of key reactions have been determined at the high-pressure limit and the falloff region. The isomerizations of hydroxyhexylperoxy to hydroperoxy-hydroxyalkyl are found to be of great importance. The formed hydroperoxy-hydroxyalkyl radicals can either dissociate rapidly into aldehydes and ketones, or be accumulated to participate in the second oxygen addition. The formation of formaldehyde, n-pentanal and OH via the Waddington mechanism is among the most important product channels. H2O elimination pathways are identified for the most stable HOQOOH intermediates, with the radical site at α-carbon to the hydroxy group. The present study recognizes those hydroperoxy-hydroxyalkyl radicals that are able to undergo the second oxygen addition. Kinetic analysis reveals that the facile formation via 1,5 or 1,6 H-transfer and lack of fast decomposition pathways combine to enable the accumulation of particular hydroperoxy-hydroxyalkyl radicals, facilitating their involvement in subsequent chain-branching chemistry. This study further illuminates the intricate understanding of the low temperature chemistry of larger alkenes.

Beam-wise dose composition learning for head and neck cancer dose prediction in radiotherapy

Automatic and accurate dose distribution prediction plays an important role in radiotherapy plan. Although previous methods can provide promising performance, most methods did not consider beam-shaped radiation of treatment delivery in clinical practice. This leads to inaccurate prediction, especially on beam paths. To solve this problem, we propose a beam-wise dose composition learning (BDCL) method for dose prediction in the context of head and neck (H&N) radiotherapy plan. Specifically, a global dose network is first utilized to predict coarse dose values in the whole-image space. Then, we propose to generate individual beam masks to decompose the coarse dose distribution into multiple field doses, called beam voters, which are further refined by a subsequent beam dose network and reassembled to form the final dose distribution. In particular, we design an overlap consistency module to keep the similarity of high-level features in overlapping regions between different beam voters. To make the predicted dose distribution more consistent with the real radiotherapy plan, we also propose a dose-volume histogram (DVH) calibration process to facilitate feature learning in some clinically concerned regions. We further apply an edge enhancement procedure to enhance the learning of the extracted feature from the dose falloff regions. Experimental results on a public H&N cancer dataset from the AAPM OpenKBP challenge show that our method achieves superior performance over other state-of-the-art approaches by significant margins. Source code is released at https://github.com/TL9792/BDCLDosePrediction.

Extension of matRad with a modified microdosimetric kinetic model for carbon ion treatment planning: Comparison with Monte Carlo calculation.

Treatment planning is essential for in silico particle therapy studies. matRad is an open-source research treatment planning system (TPS) based on the local effect model, which is a type of relative biological effectiveness (RBE) model. This study aims to implement a microdosimetric kinetic model (MKM) in matRad and develop an automation algorithm for Monte Carlo (MC) dose recalculation using the TOPAS code. In addition, we provide the developed MKM extension as open-source tool for users. Carbon beam data were generated using TOPAS MC pencil beam irradiation. We parameterized the TOPAS MC beam data with a double-Gaussian fit and modeled the integral depth doses and lateral spot profiles in the range of 100-430 MeV/u. To implement the MKM, the specific energy data table for Z=1-6 and integrated depth-specific energy data were acquired based on the Kiefer-Chatterjee track structure and TOPAS MC simulation, respectively. Generic data were integrated into matRad, and treatment planning was performed based on these data. The optimized plan parameters were automatically converted into MC simulation input. Finally, the matRad TPS and TOPAS MC simulations were compared using the RBE-weighted dose calculation results. A comparison was made for three geometries: homogeneous water phantom, inhomogeneous phantom, and patient. The RBE-weighted dose (DRBE ) distribution agreed with TOPAS MC within 1.8% for all target sizes for the homogeneous phantom. For the inhomogeneous phantom, the relative difference in the range of 80% of the prescription dose in the distal fall-off region (R80) between the matRad TPS and TOPAS MC was 0.6% (1.1mm). DRBE between the TPS and the MC was within 4.0%. In the patient case, the difference in the dose-volume histogram parameters for the target volume between the TPS and the MC was less than 2.7%. The relative difference in R80 was 0.7% (1.2mm). The MKM was successfully implemented in matRad TPS, and the RBE-weighted dose was comparable to that of TOPAS MC. The MKM-implemented matRad was released as an open-source tool. Further investigations with MC simulations can be conducted using this tool, providing a good option for carbon ion research.

Sensitivity Analysis of a 6 MeV Photon Beam Monte Carlo Model.

This study aimed to optimize efficiency in Monte Carlo (MC) simulation using sensitivity analysis of a beam model. The BEAMnrc-based model of 6 MV beam of a Siemens Primus linac was developed. For sensitivity analysis, the effect of the electron source, treatment head, and virtual phantom specifications on calculated percent depth dose (PDD) and lateral dose profiles was evaluated. The optimum mean energy (E) and the full width at half maximum (FWHM) of the intensity distribution of the electron beam were calculated as 6.7 MeV and 3 mm, respectively. Increasing E from 6.1 to 6.7 MeV, increased the PDD in the fall-off region by 4.70% and decreased the lateral profile by 8.76%. Changing the FWHM had a significant effect on the buildup region of PDD and the horns and out-of-field regions of the lateral profile. Increasing the collimators opening by 0.5 mm, PDD increased by 2.13% and the central and penumbra regions of profiles decreased by 1.98% and 11.40% respectively. Collimator properties such as thickness and density were effective in changing the penumbra (11.32% for 0.25 cm increment) and the out-of-field (22.82% for 3 g/cm3) regions of the lateral profiles. Analysis of a 6 MV model showed that PDD profiles were more sensitive to changes in energy than to FWHM of the electron source. The lateral profiles were sensitive to E, FWHM, and collimator opening. The density of the collimator affected only the out-of-field region of lateral profiles. The findings of this study may be used to make benchmarking of an MC beam model more efficient.

Characterisation and Quenching Correction for an Al2O3:C Optical Fibre Real Time System in Therapeutic Proton, Helium, and Carbon-Charged Beams.

Real time radioluminescence fibre-based detectors were investigated for application in proton, helium, and carbon therapy dosimetry. The Al2O3:C probes are made of one single crystal (1 mm) and two droplets of micro powder in two sizes (38 μm and 4 μm) mixed with a water-equivalent binder. The fibres were irradiated behind different thicknesses of solid slabs, and the Bragg curves presented a quenching effect attributed to the nonlinear response of the radioluminescence (RL) signal as a function of linear energy transfer (LET). Experimental data and Monte Carlo simulations were utilised to acquire a quenching correction method, adapted from Birks' formulation, to restore the linear dose-response for particle therapy beams. The method for quenching correction was applied and yielded the best results for the '4 μm' optical fibre probe, with an agreement at the Bragg peak of 1.4% (160 MeV), and 1.5% (230 MeV) for proton-charged particles; 2.4% (150 MeV/u) for helium-charged particles and of 4.8% (290 MeV/u) and 2.9% (400 MeV/u) for the carbon-charged particles. The most substantial deviations for the '4 μm' optical fibre probe were found at the falloff regions, with ~3% (protons), ~5% (helium) and 6% (carbon).

Microdosimetry with a mini-TEPC in the spread-out Bragg peak of 148 MeV protons

A new miniaturized Tissue Equivalent Proportional Counter (mini‒TEPC) with cylindrical sensitive volume of 1.0 mm in diameter and height, and with external diameter of 1.53 cm, has been developed to work without gas flow. With this mini counter we have measured the radiation quality of the 148 MeV energy-modulated proton beam at the radiobiological research line of the Trento Proton Therapy Centre. The irradiation field had a diameter of about 6 cm, a modulation width of 2.5 cm and a range of about 13.5 cm. At four specific depths, the response function of the new microdosimeter has also been compared to that of a commercial spherical microdosimeter, the LET-1/2” Spherical TEPC produced by Far West Technology. Afterwards, measurements were performed with the mini‒TEPC at nine precise positions across the Spread-Out Bragg Peak (SOBP).From the microdosimetric spectra the frequency and the dose‒mean lineal energy values, y‾F and y‾D, were derived for both detectors. Afterwards, microdosimetric assessment of Relative Biological Effectiveness (RBEμ) has been performed by weighting the dose distribution of the lineal energy, d(y), with the Loncol’s biological weighting function r(y).RBEμ has been used to calculate the RBEμD at different dose levels D. In the distal fall-off region the RBEμD(d) increases significantly, up to a factor 3. Finally, the RBEμD-weighted dose was calculated and compared to that obtained with a fixed RBE = 1.1. The use of a variable RBE to weight the physical dose leads to a biological range extension evaluated at R80 of about 2.0 mm, with respect to the R80 obtained with a fixed RBE = 1.1.The use of a mini‒TEPC in proton therapy has been so far limited to 60 MeV proton beams, used for the treatment of ocular melanoma. This is the first time that a mini-TEPC was used with an higher energy proton beam of 148 MeV, and the spectra measured with the mini-TEPC were compared to those measured with the commercial FWT LET1/2 counter for the first time.These new results corroborate the usability of mini-TEPCs in proton therapy.

3D high resolution clonogenic survival measurement of xrs-5 cells in low-dose region of carbon ion plans

Purpose In this study, we performed biological verification measurements of cell survival of a 12C ion irradiation plan employing a high-resolution 3D culture setup. This allowed, in particular, to access the cell inactivation in the low-dose regions close to the target area. Materials and methods We established the protocol for a 3D culture setup where xrs-5 cells were grown inside a layered matrigel structure in 384-well plates. Their radiosensitivity to conventional and 12C ion radiation was evaluated by irradiating them either with 250 kV X-rays at GSI or with monoenergetic 12C beams of 110 MeV/u at MIT, and compared with those of monolayers. A treatment plan for a rectangular target was prepared using the GSI research treatment planning system TRiP98. xrs-5 cells were seeded in the matrigel-based setup and irradiated in dose fall-off regions using active scanning 12C ion beams. In addition, film dosimetry utilizing radiochromic EBT3 film has been performed to assess the field homogeneity downstream of 384-well V-bottom plates with or without additional agarose coating of the well plate bottom. Results Dose response curves following X-ray and 12C ion irradiation had linear shape and showed a significant decrease in survival fraction at even moderate doses. Survival measurements in the low-dose regions of the plan for the extended target showed good agreement to the predicted survival fraction. The irradiated film profiles yielded a flat dose distribution without apparent artifacts or inhomogeneities for well plates both with and without agarose coating, confirming the suitability of the experimental setup. Conclusions We conclude that the V-bottom 384-well plates in combination with the radiation-sensitive xrs-5 cell line constitute a suitable radiobiological verification tool which can be used especially for low doses. Furthermore, the measured survival of xrs-5 cells show a good agreement with the expected survival in the low-dose out-of-field regions, both laterally and downstream of the target.

Accelerating and improving radiochromic film calibration by utilizing the dose ratio in photon and proton beams.

PurposeRadiochromic films are versatile 2D dosimeters with high‐resolution and near tissue equivalence. To assure high precision and accuracy, a time‐consuming calibration process is required. To improve the time efficiency, a novel calibration method utilizing the ratio of the same dose profile measured at different monitor units (MUs) is introduced and tested in a proton and photon beam.MethodsThe calibration procedure employs the dose ratio of film measurements of the same relative profile for different absolute dose values. Hence, the ratio of the dose is constant at any point of the profile, but the ratio of the net optical densities is not constant. The key idea of the method is to optimize the calibration function until the ratio of the calculated doses is constant. The proposed method was tested in the dose range between 0.25–12 and 1–6 Gy in a proton and photon beam, respectively. A radial symmetric profile and a rectangular profile were created, both having a central plateau region of about 3 cm diameter and a dose falloff of about 1.5 cm at larger distances. The dose falloff region was used as input for the optimization method and the central plateau region served as dose reference points. Only the plateau region of the highest dose entered the optimization as an additional objective. The measured data were randomly split into differently sized training and test sets. The optimization was repeated 1000 times with random start value initialization using the same start values for the standard and the gradient method. Finally, a proton plan with four dose levels was created, which were separated spatially, to test the possibility of a full calibration within a single measurement.ResultsParameter estimation was possible with as low as one dose ratio used for optimization in both the photon and the proton case, yet exhibiting a high sensitivity on the dose level. The root mean squared deviation (RMSD) of the dose was less than 1% when the dose ratio was in the order of 20, whereas the median RMSD of all optimizations was 1.7%. Using four dose levels for optimization resulted in a median RMSD of 1% when randomly selecting the dose levels. Having at least one dose ratio of about 20 included in the optimization considerably improved the RMSD of the calibration function. Using six or eight dose levels reduced the sensitivity on the dose level selection and the median RMSD was 0.8%. A full calibration was possible in a single measurement having four dose levels in one plan but spatially separated.ConclusionsThe number of measurements required to obtain an EBT3 film calibration function could be reduced using the proposed dose ratio method while maintaining the same accuracy as with the standard method.

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.

Evaluation of calibration methods of Exradin W2 plastic scintillation detector for CyberKnife small-field dosimetry

PurposePlastic scintillation detectors (PSD) require corrections of Cerenkov light for accurate radiation dosimetry. The vendor-recommended calibration method of the Exradin W1 requires irradiation with a 40 × 40 cm2 field size, although the CyberKnife robotic radiosurgery system cannot generate such large field. Recently, an Exradin W2 PSD, succeeding Exradin W1, was characterized. The W2 PSD can be calibrated in water with small curvature using a new plastic holder, which has a short pole on the bottom to mount it on the arm of the water phantom. The purpose of this study was to evaluate the calibration method of Exradin W2 using the new plastic holder for calibration with small fields and compare the impact of the calibration methods on the beam data. MethodsAn Exradin W2-1x1 detector with active volume length 1 mm and diameter 1 mm was used. The Cerenkov light ratio (CLR) values were determined by two calibration methods: measurements in a 30 × 30 cm2 solid water and in water using a small plastic holder. Percent-depth-dose (PDD), beam profiles, and detector output factor (OFdet) of TrueBeam 6 MV photon beams with flattening filter (WFF) and flattening filter-free (FFF) beams were measured to compare these two calibration methods. For the CyberKnife photon beams, PDD, beam profiles, and OFdet were measured with the calibration in water, and the data were compared with those measured using a Diode SRS (PTW model 60018) and microDiamond (PTW model 60019). ResultsThe difference of the CLR between two calibration methods were within 3%. The effect of the calibration methods on the PDD, beam profiles, and OFdet of the TrueBeam were within 0.06%, 0.07%, and 0.8%, respectively. For the CyberKnife beams, the maximum difference of PDD in the dose fall-off region, full-width half-maximum (FWHM) of beam profiles, and penumbra among the three detectors were 1.3%, 0.25 mm, and 0.35 mm, respectively. The relative difference of the W2-measured OFdet from the microDiamond-measured field output factor were within 1.7%. ConclusionsThe calibration in water with a small-field size provided equivalent data compared to those obtained with a conventional calibration method. The Exradin W2 provided appropriate data for CyberKnife small-field dosimetry.

Comprehensive Commissioning and Clinical Implementation of GammaTiles STaRT for Intracranial Brain Cancer

PurposeTo validate the dose calculation accuracy and dose distribution of GammaTiles for brain tumors, and to suggest a surgically targeted radiation therapy (STaRT) workflow for planning, delivery, radiation safety documentation, and posttreatment validation. Methods and MaterialsNovel surgically targeted radiation therapy, GammaTiles, uses Cs-131 radiation isotopes embedded in collagen-based tiles that can be resorbed after surgery. GammaTile target delineation and dose calculation were performed on MIM Symphony software. Point-based and complex seed distribution calculations in MIM Symphony were verified with hand calculations and BrachyVision calculations. Vendor-provided 2-dimensional dose distribution calculation accuracy was validated using gafchromic EBT3 film measurements at various depths. A workflow was established for safe and effective GammaTile implants. ResultsGood agreement was observed between different calculations. Calculation accuracy of less than 0.5% was achieved for all points except one, which had rounding issues for very low doses and resulted in just below 5% difference. Differences in anisotropy and geometry positioning were noticed in the delineation of Cs-131 IsoRay seeds in the compared systems, resulting in minor discrepancies in the calculated dosimetry distributions. Film measurements showed profiles with relatively good agreement of 0% to 5% in nongradient regions with higher differences between 5% to 10% in the sharp dose fall-off regions. ConclusionsA comprehensive evaluation of GammaTile geometry, dose distribution, and clinical workflow was conducted. Safe intro-operative implantation of GammaTiles requires extensive preplanning and interdisciplinary collaboration. A STaRT workflow was outlined to provide a guideline for an accurate treatment planning and safe implant process at other institutions.

New reactions of diazene and related species for modelling combustion of amine fuels

Potential energy surfaces for reactions involving NH isomers of diazene (diimide) have been explored using density functional theory, with energies based on coupled-cluster theory. A focus is on processes that create or consume these species, and isomerisation between the E (trans) and Z (cis) forms of HNNH. These include isomerisation and dissociation pathways for HNNH, addition of H atoms to form NH, abstraction by H atoms yielding short-lived NNH, and abstraction reactions of H with NH. Transition state and capture theories are applied for high-pressure-limiting behaviour, while low-pressure and falloff regions are characterised via the methods of Troe and coworkers. Rate constants and thermochemistry are provided to improve models of diamine chemistry, relevant to the combustion of NH especially at high concentrations, high pressures or under reducing conditions. Results indicate that amine radical recombination mainly yields the E HNNH isomer, while H-abstraction from NH results in E HNNH and HNN. However, at elevated temperature E Z isomerisation becomes competitive, and Z HNNH, being more reactive, acts to enhance the diazene consumption rate.

Reachability Deficits in Quantum Approximate Optimization of Graph Problems

The quantum approximate optimization algorithm (QAOA) has become a cornerstone of contemporary quantum applications development. Here we show that the density of problem constraints versus problem variables acts as a performance indicator. Density is found to correlate strongly with approximation inefficiency for fixed depth QAOA applied to random graph minimization problem instances. Further, the required depth for accurate QAOA solution to graph problem instances scales critically with density. Motivated by Google's recent experimental realization of QAOA, we preform a reanalysis of the reported data reproduced in an ideal noiseless setting. We found that the reported capabilities of instances addressed experimentally by Google, approach a rapid fall-off region in approximation quality experienced beyond intermediate-density. Our findings offer new insight into performance analysis of contemporary quantum optimization algorithms and contradict recent speculation regarding low-depth QAOA performance benefits.

Extension and validation of a GPU-Monte Carlo dose engine gDPM for 1.5 T MR-LINAC online independent dose verification.

To extend and validate the accuracy and efficiency of a graphics processing unit (GPU)-Monte Carlo dose engine for Elekta Unity 1.5 T Magnetic Resonance-Linear Accelerator (MR-LINAC) online independent dose verification. Electron/positron propagation physics in a uniform magnetic field was implemented in a previously developed GPU-Monte Carlo dose engine-gDPM. The dose calculation accuracy in the magnetic field was first evaluated in heterogeneous phantom with EGSnrc. The dose engine was then commissioned to a Unity machine with a virtual two photon-source model and compared with the Monaco treatment planning system. Fifteen patient plans from five tumor sites were included for the quantification of online dose verification accuracy and efficiency. The extended gDPM accurately calculated the dose in a 1.5 T external magnetic field and was well matched with EGSnrc. The relative dose difference along central beam axis was less than 0.5% for the homogeneous region in water-lung phantom. The maximum difference was found at the build-up regions and heterogeneous interfaces, reaching 1.9% and 2.4% for 2 and 6MeV mono-energy photon beams, respectively. The root mean square errors for depth-dose fall-off region were less than 0.2% for all field sizes and presented a good match between gDPM and Monaco GPUMCD. For in-field profiles, the dose differences were within 1% for cross-plane and in-plane directions for all calculated depths except dmax. For penumbra regions, the distance-to-agreements between two dose profiles were less than 0.1cm. For patient plan verification, the maximum relative average dose difference was 1.3%. The gamma passing rates with criteria 3% (2mm) for dose regions above 20% were between 93% and 98%. gDPM can complete the dose calculation for less than 40 s with 5 × 108 photons on a single NVIDIA GTX-1080Ti GPU and achieve a statistical uncertainty of 0.5%-1.1% for all evaluated cases. A GPU-Monte Carlo package-gDPM was extended and validated for Elekta Unity online plan verification. Its calculation accuracy and efficiency make it suitable for online independent dose verification for MR-LINAC.

Characterization of the ZAP-X® Peripheral Dose Fall-Off.

Various small-field radiation dose detectors were systematically compared and their impact on measured beam performance of the ZAP-X® dedicated stereotactic radiosurgery system (ZAP Surgical Systems, Inc., San Carlos, CA, USA) was determined. Three Physikalische Technische Werkstaetten (PTW) diodes, i.e., the microSilicon, the microDiamond, and the Stereotactic Radiosurgery (SRS) diode detectors of (PTW-Freiburg, Freiburg, Germany), as well as Gafchromic™ External Beam Therapy 3 (EBT) film (Ashland, Inc., Wilmington, DE, USA), were used and compared to arrive at a recommended standard for this critical component of small-field beam measurements. Beam profiles, including the dose fall-off region near the edge of the beam, were measured with the PTW diodes and EBT3 film and subsequently contrasted. The impact of detector physical and dosimetric characteristics on the results of the measurements was investigated and compared with film measurements. The beam penumbra was used to quantify the dose fall-off. The measurement acquired with the diodes and film showed the most significant differences in the fall-off region near the field edge. The film-based measurements clearly showed the steepest dose gradient verified by the penumbra value of 1.21 mm, followed by the SRS diode with 1.60 mm, the microSilicon diode with 1.67 mm, and the microDiamond diode with 1.83 mm. A clear correlation of each detector’s sensitive area with the penumbra was found, with the microDiamond detector at 2.2 mm diameter sensitive area having the largest penumbra, followed by the microSilicon and SRS diodes. Beam measurements for the purposes of system characterization or treatment planning system beam data acquisition depend, to a large extent, on detector characteristics. This is especially true for small-field dosimetry performed during stereotactic radiosurgery beam measurements. Careful consideration should be practiced which allows for the measurements to represent true beam characteristics and minimize the impact of the detector on the measurements. We conclude that film should be considered the reference method for such measurements with the ZAP-X due to its smallest physical measurement resolution of 23.1 µm. Potential drawbacks to this methodology are the need to calibrate the film relative to the dose and possible problems with saturation and non-linear film response for very high and very low optical densities.

Application of radiochromic gel dosimetry to commissioning of a megavoltage research linear accelerator for small-field animal irradiation studies.

To develop and implement an efficient and accurate commissioning procedure for small-field static beam animal irradiation studies on an MV research linear accelerator (Linatron-M9) using radiochromic gel dosimetry. The research linear accelerator (Linatron-M9) is a 9MV linac with a static fixed collimator opening of 5.08cm diameter. Lead collimators were manually placed to create smaller fields of 2×2cm2 , 1×1cm2 , and 0.5×0.5cm2 . Relative dosimetry measurements were performed, including profiles, percent depth dose (PDD) curves, beam divergence, and relative output factors using various dosimetry tools, including a small volume ionization chamber (A14), GAFCHROMIC™ EBT3 film, and Clearview gel dosimeters. The gel dosimeter was used to provide a 3D volumetric reference of the irradiated fields. The Linatron profiles and relative output factors were extracted at a reference depth of 2cm with the output factor measured relative to the 2×2cm2 reference field. Absolute dosimetry was performed using A14 ionization chamber measurements, which were verified using a national standards laboratory remote dosimetry service. Absolute dosimetry measurements were confirmed within 1.4% (k=2, 95% confidence=5%). The relative output factor of the small fields measured with films and gels agreed with a maximum relative percent error difference between the two methods of 1.1 % for the 1×1cm2 field and 4.3 % for the 0.5×0.5cm2 field. These relative errors were primarily due to the variability in the collimator positioning. The measured beam profiles demonstrated excellent agreement for beam size (measured as FWHM), within approximately 0.8mm (or less). Film measurements were more accurate in the penumbra region due to the film's finer resolution compared with the gel dosimeter. Following the van Dyk criteria, the PDD values of the film and gel measurements agree within 11% in the buildup region starting from 0.5cm depth and within 2.6 % beyond maximum dose and into the fall-off region for depths up to 5cm. The 2D beam profile isodose lines agree within 0.5mm in all regions for the 0.5×0.5cm2 and the 1×1cm2 fields and within 1mm for the larger field of 2×2cm2 . The 2D PDD curves agree within approximately 2% of the maximum in the typical therapy region (1-4cm) for the 1×1cm2 and 2×2cm2 and within 5% for the 0.5×0.5cm2 field. This work provides a commissioning process to measure the beam characteristics of a fixed beam MV accelerator with detailed dosimetric evaluation for its implementation in megavoltage small animal irradiation studies. Radiochromic gel dosimeters are efficient small-field relative dosimetry tools providing 3D dose measurements allowing for full representation of dose, dosimeter misalignment corrections and high reproducibility with low inter-dosimeter variability. Overall, radiochromic gels are valuable for fast, full relative dosimetry commissioning in comparison to films for application in high-energy small-field animal irradiation studies.