A hybrid semiconductor detector consisting of a 65 mu m thick 4H-SiC sensor bump-bonded to a Timepix3 ASIC was used to characterize carbon-and helium-ion beams (up to 402 MeV/u). Compared to a 500 mu m thick Si sensor, the SiC detector shows an 87 +/- 3% reduction in pixel occupancy for 290 MeV/u carbon, preventing ToT saturation while maintaining tracking properties; this corresponds to approximate to 7.7x higher tolerable particle flux at matched settings. For helium, mean cluster sizes of 1.0 px (SiC) versus 5.0 px (Si) imply 80% lower occupancy and approximate to 5x higher flux headroom. For ultra-high dose-rate (UHDR/FLASH) applications, out-of-field monitoring of 226 MeV pencil-beam scanning deliveries (169 spots, of similar to 2 ms duration) achieved 100% spots detection efficiency and a linear dose-rate response across 100-500 nA, with sensitivity down to 10 mGy/s. Overall, SiC-based MiniPIX-Timepix3 enables compact, non-invasive, particle-resolved beam monitoring with improved dynamic-range utilization and scalable flux tolerance, supporting advanced in-room QA across clinical and UHDR conditions.

BackgroundCarbon ion beams are well accepted as densely ionizing radiation with a high linear energy transfer (LET). However, the current clinical practice does not fully exploit the highest possible dose-averaged LET (LETd) and, consequently, the biological potential in the target. This aspect becomes worse in larger tumors for which inferior clinical outcomes and corresponding lower LETd was reported.PurposeThe vicinity to critical organs in general and the inferior overall survival reported for larger sacral chordomas treated with carbon ion radiotherapy (CIRT), makes the treatment of such tumors challenging. In this work it was aimed to increase the LETd in large volume tumors while maintaining the relative biological effectiveness (RBE)-weighted dose, utilizing the LETd optimization functions of a commercial treatment planning system (TPS).MethodsTen reference sequential boost carbon ion treatment plans, designed to mimic clinical plans for large sacral chordoma tumors, were generated. High dose clinical target volumes (CTV-HD) larger than 250cm3$250 \,{\rm cm}<^>{3}$ were considered as large targets. The total RBE-weighted median dose prescription with the local effect model (LEM) was DRBE,50%=73.6Gy$\textrm {D}_{\rm RBE, 50\%}=73.6 \,{\rm Gy}$ in 16 fractions (nine to low dose and seven to high dose planning target volume). No LETd optimization was performed in the reference plans, while LETd optimized plans used the minimum LETd (Lmin) optimization function in RayStation 2023B. Three different Lmin values were investigated and specified for the seven boost fractions: Lmin=60keV/mu m$\textrm {L}_{\rm min}=60 \,{\rm keV}/{\umu }{\rm m}$, Lmin=80keV/mu m$\textrm {L}_{\rm min}=80 \,{\rm keV}/{\umu }{\rm m}$ and Lmin=100keV/mu m$\textrm {L}_{\rm min}=100 \,{\rm keV}/{\umu }{\rm m}$. To compare the LETd optimized against reference plans, LETd and RBE-weighted dose based goals similar to and less strict than clinical ones were specified for the target. The goals for the organs at risk (OAR) remained unchanged. Robustness evaluation was studied for eight scenarios (+/- 3.5%$\pm 3.5\%$ range uncertainty and +/- 3mm$\pm 3 \,{\rm mm}$ setup uncertainty along the main three axes).ResultsThe optimization method with Lmin=60keV/mu m$\textrm {L}_{\rm min}=60 \,{\rm keV}/{\umu }{\rm m}$ resulted in an optimal LETd distribution with an average increase of LETd,98%${\rm {LET}}_{{\rm {d,}}98\%}$ (and LETd,50%${\rm {LET}}_{{\rm {d,}}50\%}$) in the CTV-HD by 8.9 +/- 1.5keV/mu m$8.9\pm 1.5 \,{\rm keV}/{\umu }{\rm m}$ (27%$27\%$) (and 6.9 +/- 1.3keV/mu m$6.9\pm 1.3 \,{\rm keV}/{\umu }{\rm m}$ (17%$17\%$)), without significant difference in the RBE-weighted dose. By allowing +/- 5%$\pm 5\%$ over- and under-dosage in the target, the LETd,98%${\rm {LET}}_{{\rm {d,}}98\%}$ (and LETd,50%${\rm {LET}}_{{\rm {d,}}50\%}$) can be increased by 11.3 +/- 1.2keV/mu m$11.3\pm 1.2 \,{\rm keV}/{\umu }{\rm m}$ (34%$34\%$) (and 11.7 +/- 3.4keV/mu m$11.7\pm 3.4 \,{\rm keV}/{\umu }{\rm m}$ (29%$29\%$)), using the optimization parameters Lmin=80keV/mu m$\textrm {L}_{\rm min}=80 \,{\rm keV}/{\umu }{\rm m}$. The pass rate for the OAR goals in the LETd optimized plans was in the same level as the reference plans. LETd optimization lead to less robust plans compared to reference plans.ConclusionsCompared to conventionally optimized treatment plans, the LETd in the target was increased while maintaining the RBE-weighted dose using TPS LETd optimization functionalities. Regularly assessing RBE-weighted dose robustness and acquiring more in-room images remain crucial and inevitable aspects during treatment.

BACKGROUND:Large tumor size has been reported as a predicting factor for inferior clinical outcome in carbon ion radiotherapy (CIRT). Besides the clinical factors accompanied with such tumors, larger tumors receive typically more low linear energy transfer (LET) contributions than small ones which may be the underlying physical cause. Although dose averaged LET is often used as a single parameter descriptor to quantify the beam quality, there is no evidence that this parameter is the optimal clinical predictor for the complex mixed radiation fields in CIRT. PURPOSE:Purpose of this study was to investigate on a novel dosimetric quantity, namely high-LET-dose ( D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ , the physical dose filtered based on an LET threshold) as a single parameter estimator to differentiate between carbon ion treatment plans (cTP) with a small and large tumor volume. METHODS:Ten cTPs with a planning target volume, PTV ≥ 500 cm 3 $\mathrm{PTV}\ge {500}\,{{\rm cm}^{3}}$ (large) and nine with a PTV < 500 cm 3 $\mathrm{PTV}<{500}\,{{\rm cm}^{3}}$ (small) were selected for this study. To find a reasonable LET threshold ( L thr $\textrm {L}_{\textrm {thr}}$ ) that results in a significant difference in terms of D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ , the voxel based normalized high-LET-dose ( D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ ) distribution in the clinical target volume (CTV) was studied on a subset (12 out of 19 cTPs) for 18 LET thresholds, using standard distribution descriptors (mean, variance and skewness). The classical dose volume histogram concept was used to evaluate the D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ and D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ distributions within the target of all 19 cTPs at the before determined L thr $\textrm {L}_{\textrm {thr}}$ . Statistical significance of the difference between the two groups in terms of mean D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ and D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ volume histogram parameters was evaluated by means of (two-sided) t-test or Mann-Whitney-U-test. In addition, the minimum target coverage at the above determined L thr $\textrm {L}_{\textrm {thr}}$ was compared and validated against three other thresholds to verify its potential in differentiation between small and large volume tumors. RESULTS:An L thr $\textrm {L}_{\textrm {thr}}$ of approximately 30 keV / μ m ${30}\,{\rm keV/}\umu {\rm m}$ was found to be a reasonable threshold to classify the two groups. At this threshold, the D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ and D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ were significantly larger ( p < 0.05 $p<0.05$ ) in small CTVs. For the small tumor group, the near-minimum and median D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ (and D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ ) in the CTV were in average 9.3 ± 1.5 Gy $9.3\pm {1.5}\,{\rm Gy}$ (0.31 ± 0.08) and 13.6 ± 1.6 Gy $13.6\pm {1.6}\,{\rm Gy}$ (0.46 ± 0.06), respectively. For the large tumors, these parameters were 6.6 ± 0.2 Gy $6.6\pm {0.2}\,{\rm Gy}$ (0.20 ± 0.01) and 8.6 ± 0.4 Gy $8.6\pm {0.4}\,{\rm Gy}$ (0.28 ± 0.02). The difference between the two groups in terms of mean near-minimum and median D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ ( D ̂ > L thr $\hat{\textrm {D}}_{>\textrm {L}_{\textrm {thr}}}$ ) was 2.7 Gy (11%) and 5.0 Gy (18%), respectively. CONCLUSIONS:The feasibility of high-LET-dose based evaluation was shown in this study where a lower D > L thr $\textrm {D}_{>\textrm {L}_{\textrm {thr}}}$ was found in cTPs with a large tumor size. Further investigation is needed to draw clinical conclusions. The proposed methodology in this work can be utilized for future high-LET-dose based studies.

Silicon Carbide is a suitable semiconductor sensor for radiation measurement and nuclear applications. It has been recently implemented as a position-sensitive and radiation imaging device coupled to the Timepix3 ASIC chip in the form of a miniaturized radiation camera MiniPIX-Timepix3 SiC. In this work we systematically evaluate the detection resolving power to different radiation species: electrons, fast neutrons, ions and protons. Experimental calibrations were made at well-defined reference fields in terms of particle type, energy and direction. The spectral-sensitive tracking response and characteristic morphology of the particle tracks are analyzed by pattern recognition algorithms. The low thickness of the radiation sensitive volume (65 mu m) of the SiC sensor limits the directional tracking response and angular resolution. Three broad classes of particle-type events are resolved. The detector together with suitable data processing can be used for radiation dosimetry and particle tracking tasks in space, particle therapy, nuclear physics and nuclear reactors and particle accelerator environments.

Introduction: Particle therapy using pencil beam scanning (PBS) faces large uncertainties related to ranges and target motion. One possibility to improve existing mitigation strategies is a 2D range modulator (2DRM). A 2DRM offers faster irradiation times by reducing the number of layers and spots needed to create a spread-out Bragg peak. We have investigated the impact of 2DRM on microdosimetric spectra measured in proton and carbon ion beams. Materials and Methods: Two 2DRMs were designed and 3D printed, one for. 124.7 MeV protons and one for 238.6 MeV/u carbon ions. Their dosimetric validation was performed using Roos and PinPoint ionization chamber and EBT3 films. Monte Carlo simulations were done using GATE. A silicon-based solid-state microdosimeter was used to collect pulse-height spectra along three depths for two irradiation modalities, PBS and a single central beam. Results: For both particle types, the original pin design had to be optimized via GATE simulations. The difference between the R80 of the simulated and measured depth dose curve was 0.1 mm. The microdosimetric spectra collected with the two irradiation modalities overlap well. Their mean lineal energy values differ over all positions by 5.2 % for the proton 2DRM and 2.1 % for the carbon ion 2DRM. Conclusion: Radiation quality in terms of lineal energy was independent of the irradiation method. This supports the current approach in reference dosimetry, where the residual range is chosen as a beam quality index to select stopping power ratios. (C) 2023 The Author(s). Published by Elsevier B.V. Radiotherapy and Oncology 182 (2023) 1-6 This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

BACKGROUND The increasing number of studies dealing with linear energy transfer (LET)-based evaluation and optimization in the field of carbon ion radiotherapy (CIRT) indicates the rising demand for LET implementation in commercial treatment planning systems (TPS). Benchmarking studies could play a key role in detecting (and thus preventing) computation errors prior implementing such functionalities in a TPS. PURPOSE This in silico study was conducted to benchmark the following two LET-related functionalities in a commercial TPS against Monte Carlo simulations: (1) dose averaged LET (LETd ) scoring and (2) physical dose filtration based on LET for future LET-based treatment plan evaluation and optimization studies. METHODS The LETd scoring and LET-based dose filtering (in which the deposited dose can be separated into the dose below and above the user specified LET threshold) functionalities for carbon ions in the research version RayStation (RS) 9A-IonPG TPS (RaySearch Laboratories, Sweden) were benchmarked against GATE/Geant4 simulations. Pristine Bragg peaks (BPs) and cuboid targets, positioned at different depths in a homogeneous water phantom and a setup with heterogeneity were used for this study. RESULTS For all setups (homogeneous and heterogeneous), the mean absolute (and relative) LETd difference was less than 1 keV/ μ $\umu$ m (3.5%) in the plateau and target and less than 2 keV/ μ $\umu$ m (8.3%) in the fragmentation tail. The maximum local differences were 4 and 6 keV/ μ $\umu$ m, respectively. The mean absolute (and relative) physical dose differences for both low- and high-LET doses were less than 1 cGy (1.5%) in the plateau, target and tail with a maximum absolute difference of 2 cGy. CONCLUSIONS No computation error was found in the tested functionalities except for LETd in lateral direction outside the target, showing the limitation of the implemented monochrome model in the tested TPS version.

The aim of this work is to present the first microdosimetric spectra measured with a miniaturised tissue-equivalent proportional counter in the clinical environment of the MedAustron ion-beam therapy facility. These spectra were gathered with a 62.4-MeV proton beam and have been compared with microdosimetric spectra measured in the 62-MeV clinical proton beam of the CATANA beam line. Monte Carlo simulations were performed using the Geant4 toolkit GATE and a fully commissioned clinical beam line model. Finally, similarities and discrepancies of the measured data to simulations based on a simple and complex detector geometry are discussed.

Introduction: A growing interest in using proton pencil beam scanning in combination with collimators for the treatment of small, shallow targets, such as ocular melanoma or pre-clinical research emerged recently. This study aims at demonstrating that the dose of a synchrotron-based PBS system with a dedicated small, shallow field nozzle can be accurately predicted by a commercial treatment planning system (TPS) following appropriate tuning of both, nozzle and TPS.Materials: A removable extension to the clinical nozzle was developed to modify the beam shape passively. Five circular apertures with diameters between 5 to 34 mm, mounted 72 cm downstream of a range shifter were used. For each collimator treatment plans with spread-out Bragg peaks (SOBP) with a modulation of 3 to 30 mm were measured and calculated with GATE/Geant4 and the research TPS RayStation (RS11B-R). The dose grid, multiple coulomb scattering and block discretization resolution were varied to find the optimal balance between accuracy and performance.Results: For SOBPs deeper than 10 mm, the dose in the target agreed within 1% between RS11B-R, GATE/Geant4 and measurements for aperture diameters between 8 to 34 mm, but deviated up to 5% for smaller apertures. A plastic taper was introduced reducing scatter contributions to the patient (from the pipe) and improving the dose calculation accuracy of the TPS to a 5% level in the entrance region for large apertures.Conclusion: The commercial TPS and GATE/Geant4 can accurately calculate the dose for shallow, small proton fields using a collimator and pencil beam scanning.

Purpose: To evaluate the dosimetric accuracy for small field proton irradiation relevant for pre-clinical in vivo studies using clinical infrastructure and technology. In this context additional beam collimation and range reduction was implemented. Methods and materials: The clinical proton beam line employing pencil beam scanning (PBS) was adapted for the irradiation of small fields at shallow depths. Cylindrical collimators with apertures of 15, 12, 7 and 5 mm as well as two different range shifter types, placed at different distances relative to the target, were tested: a bolus range shifter (BRS) attached to the collimator and a clinical nozzle mounted range shifter (CRS) placed at a distance of 72 cm from the collimator. The Monte Carlo (MC) based dose calculation engine implemented in the clinical treatment planning system (TPS) was commissioned for these two additional hardware components. The study was conducted with a phantom and cylindrical target sizes between 2 and 25 mm in diameter following a dosimetric end-to-end test concept. Results: The setup with the CRS provided a uniform dose distribution across the target. An agreement of better than 5% between the planned dose and the measurements was obtained for a target with 3 mm diameter (collimator 5 mm). A 2 mm difference between the collimator and the target diameter (target being 2 mm smaller than the collimator) sufficed to cover the whole target with the planned dose in the setup with CRS. Using the BRS setup (target 8 mm, collimator 12 mm) resulted in non-homogeneous dose distributions, with a dose discrepancy of up to 10% between the planned and measured doses. Conclusion: The clinical proton infrastructure with adequate beam line adaptations and a state-of-the-art TPS based on MC dose calculations enables small animal irradiations with a high dosimetric precision and accuracy for target sizes down to 3 mm.

Purpose Radiochromic 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. Methods The 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. Results Parameter 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. Conclusions The 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.

PURPOSE:The Local Effect Model version one (LEM I) is applied clinically across Europe to quantify the relative biological effectiveness (RBE) of carbon ion beams. It requires the full particle fluence spectrum differential in energy in each voxel as input parameter. Treatment planning systems (TPSs) use beamline-specific look-up tables generated with Monte Carlo (MC) codes. In this study, the changes in RBE weighted dose were quantified using different levels of details in the simulation or different MC codes.METHODS:The particle fluence differential in energy was simulated with FLUKA and Geant4 at 500 depths in water in 1-mm steps for 58 initial carbon ion energies (between 120.0 and 402.8 MeV/u). A dedicated beam model was applied, including the full description of the Nozzle using GATE-RTionV1.0 (Geant4.10.03p03). In addition, two tables generated with FLUKA were compared. The starting points of the FLUKA simulations were phase space (PhS) files from, firstly, the Geant4 nozzle simulations, and secondly, a clinical beam model where an analytic approach was used to mimic the beamline. Treatment plans (TPs) were generated with RayStation 8B (RaySearch Laboratories AB, Sweden) for cubic targets in water and 10 clinical patient cases using the clinical beam model. Subsequently, the RBE weighted dose was re-computed using the two other fluence tables (FLUKA PhS or Geant4).RESULTS:The fluence spectra of the primary and secondary particles simulated with Geant4 and FLUKA generally agreed well for the primary particles. Differences were mainly observed for the secondary particles. Interchanging the two energy spectra (FLUKA vs. GEANT4) to calculate the RBE weighted dose distributions resulted in average deviations of less than 1% in the entrance up to the end of the target region, with a maximum local deviation at the distal edge of the target. In the fragment tail, larger discrepancies of up to 5% on average were found for deep-seated targets. The patient and water phantom cases demonstrated similar results.CONCLUSION:RBE weighted doses agreed well within all tested setups, confirming the clinical beam model provided by the TPS vendor. Furthermore, the results showed that the open source and generally available MC code Geant4 (in particular using GATE or GATE-RTion) can also be used to generate basic beam data required for RBE calculation in carbon ion therapy.

The central Gaussian shaped high dose region of a pencil beam (PB) in light ion beam therapy (LIBT) is enveloped by a low dose region causing non-negligible field size effects and impairs the dose calculation accuracy considerably if the low dose envelope is not well modeled. The purpose of this study was to calculate the practical radius, Rc, at which a PB does not influence a field more than a certain accuracy level. Lateral dose profiles of proton beams in water were simulated using GATE/Geant4. Those lateral dose profiles were integrated numerically and used to calculate field size factors (FSFs). The Rc was then determined such, that the lateral dose at radii exceeding Rc can be neglected without compromising the FSF of a 20cm×20cm field more than a desired accuracy level c. The practical radius Rc yielding c=0.5% was compared to the frequently applied concept of full width at a ratio x of the maximum (FWxM). The sensitivity to variations of the beam width was tested by increasing the initial beam width σC of the clinical beam model by 0.5 and 1mm, respectively. Neglecting the dose at radii exceeding Rc resulted in the desired FSF accuracy, whereas using the FW0.01%M cut resulted in varying accuracy. In order to yield a constant FSF accuracy, the ratio x in FWxM ranged from 0.003% to 0.065% of the maximum. In contrast to Rc, FWxM was sensitive to variations of the initial beam width. The maximum Rc over all depths was less than 7cm for the low(62.4MeV) and medium(148.2MeV) proton energy beam, which suggests that a plane parallel ionization chamber exceeding that radius is sufficient to acquire laterally integrated depth dose distributions for those energies. However, this holds not true for the highest energy (252.7MeV) or when including a range shifter (RaShi). The values of Rc are specific to our beam line configuration as the maximum Rc was depending on both, the scattering material in the Nozzle as well as the distance of the air-gap between Nozzle and phantom.

Purpose The accurate knowledge of the effective point of measurement (P-eff) is particularly important for measurements in proximity to high dose gradients such as in the distal fall-off of particle beams. For plane-parallel ionization chambers (ICs), P-eff is well known and located at the center of the inner surface of the entrance window. For cylindrical ICs, P-eff is shifted from the chamber's center toward the beam source. According to IAEA TRS-398, this shift can be calculated as 0.75 center dot r(cyl) for light ions with r(cyl) being the radius of the cavity. For proton beams and in absence of a dose gradient, no shift is recommended. We have experimentally determined P-eff for the 0.125 cc Semiflex IC in both proton and carbon ion beams. Methods The first method consisted of simultaneous irradiation of a plane-parallel IC and the Semiflex in a 4-cm wide spread-out Bragg peak. In the second method, a single-energy beam was used, and both ICs were positioned successively at the same measurement depths. For both approaches, the shift of the distal edge of the depth ionization distributions recorded by the two chambers at different reference points was used to calculate P-eff of the Semiflex. Both methods were applied in carbon ion beams, and only the latter was applied in proton beams. Results Both methods yielded a similar P-eff for carbon ions, 0.88 center dot r(cyl), and 0.84 center dot r(cyl), which results in a difference of only 0.1 mm. The difference to the recommended value of 0.75 center dot r(cyl) is 0.4 and 0.3 mm, respectively, which is larger than the positioning uncertainty. In the proton beam, a P-eff of 0.92 center dot r(cyl) was obtained. Conclusions The P-eff for the 0.125 cc Semiflex IC is shifted further from the cavity center as recommended by IAEA TRS-398 for light ions, with the shift for proton beams being even larger than for carbon ion beams.

Purpose To develop a computer‐driven and thus less user‐dependent method, allowing for a simple and straightforward generation of a Monte Carlo (MC) beam model of a scanned proton and carbon ion beam delivery system. Methods In a first step, experimental measurements were performed for proton and carbon ion energies in the available energy ranges. Data included depth dose profiles measured in water and spot sizes in air at various isocenter distances. Using an automated regularization‐based optimization process (AUTO‐BEAM), GATE/Geant4 beam models of the respective beam lines were generated. These were obtained sequentially by using least square weighting functions with and without regularization, to iteratively tune the beam parameters energy, energy spread, beam sigma, divergence, and emittance until a user‐defined agreement was reached. Based on the parameter tuning for a set of energies, a beam model was semi‐automatically generated. The resulting beam models were validated for all centers comparing to independent measurements of laterally integrated depth dose curves and spot sizes in air. For one representative center, three‐dimensional dose cubes were measured and compared to simulations. The method was applied on one research as well as four different clinical beam lines for proton and carbon ions of three different particle therapy centers using synchrotron or cyclotron accelerator systems: (a) MedAustron ion therapy center, (b) University Proton Therapy Dresden, and (c) Center Antoine Lacassagne Nice. Results Particle beam ranges in the MC beam models agreed on average within 0.2 mm compared to measurements for all energies and beam lines. Spot sizes in air (full‐width at half maximum) at all positions differed by less than 0.4% from the measurements. Dose calculation with the beam model for the clinical beam line at MedAustron agreed better than 1.7% in absolute dose for a representative clinical case treated with protons. For protons, beam model generation, including geometry creation, data conversion, and validation, was possible within three working days. The number of iterations required for the optimization process to converge, was found to be similar for all beam line geometries and particle types. Conclusion The presented method was demonstrated to work independently of the beam optics behavior of the different beam lines, particle types, and geometries. Furthermore, it is suitable for non‐expert users and requires only limited user interaction. Beam model validation for different beam lines based on different beam delivery systems, showed good agreement.