First Monte Carlo beam model for ultra-high dose rate radiotherapy with a compact electron LINAC

Background: FLASH radiotherapy based on ultra-high dose rate (UHDR) is actively being studied by the radiotherapy community. Dedicated UHDR electron devices are currently a mainstay for FLASH studies. Purpose: To present the first Monte Carlo (MC) electron beam model for the UHDR capable Mobetron (FLASH-IQ) as a dose calculation and treatment planning platform for preclinical research and FLASH-radiotherapy (RT) clinical trials. Methods: The initial beamline geometry of the Mobetron was provided by the manufacturer, with the first-principal implementation realized in the Geant4-based GAMOS MC toolkit. The geometry and electron source characteristics, such as energy spectrum and beamline parameters, were tuned to match the central-axis percentage depth dose (PDD) and lateral profiles for the pristine beam measured during machine commissioning. The thickness of the small foil in secondary scatter affected the beam model dominantly and was fine tuned to achieve the best agreement with commissioning data. Validation of the MC beam modeling was performed by comparing the calculated PDDs and profiles with EBT-XD radiochromic film measurements for various combinations of applicators and inserts. Results: The nominal 9 MeV electron FLASH beams were best represented by a Gaussian energy spectrum with mean energy of 9.9 MeV and variance (sigma) of 0.2 MeV. Good agreement between the MC beam model and commissioning data were demonstrated with maximal discrepancy < 3% for PDDs and profiles. Hundred percent gamma pass rate was achieved for all PDDs and profiles with the criteria of 2 mm/3%. With the criteria of 2 mm/2%, maximum, minimum and mean gamma pass rates were (100.0%, 93.8%, 98.7%) for PDDs and (100.0%, 96.7%, 99.4%) for profiles, respectively. Conclusions: A validated MC beam model for the UHDR capable Mobetron is presented for the first time. The MC model can be utilized for direct dose calculation or to generate beam modeling input required for treatment planning systems for FLASH-RT planning. The beam model presented in this work should facilitate translational and clinical FLASH-RT for trials conducted on the Mobetron FLASH-IQ platform.