Multiple-source models for electron beams of a medical linear accelerator using BEAMDP computer code

AimThe aim of this work was to develop multiple-source models for electron beams of the NEPTUN 10PC medical linear accelerator using the BEAMDP computer code.BackgroundOne of the most accurate techniques of radiotherapy dose calculation is the Monte Carlo (MC) simulation of radiation transport, which requires detailed information of the beam in the form of a phase-space file. The computing time required to simulate the beam data and obtain phase-space files from a clinical accelerator is significant. Calculation of dose distributions using multiple-source models is an alternative method to phase-space data as direct input to the dose calculation system.Materials and methodsMonte Carlo simulation of accelerator head was done in which a record was kept of the particle phase-space regarding the details of the particle history. Multiple-source models were built from the phase-space files of Monte Carlo simulations. These simplified beam models were used to generate Monte Carlo dose calculations and to compare those calculations with phase-space data for electron beams.ResultsComparison of the measured and calculated dose distributions using the phase-space files and multiple-source models for three electron beam energies showed that the measured and calculated values match well each other throughout the curves.ConclusionIt was found that dose distributions calculated using both the multiple-source models and the phase-space data agree within 1.3%, demonstrating that the models can be used for dosimetry research purposes and dose calculations in radiotherapy.     

Establishing the impact of temporary tissue expanders on electron and photon beam dose distributions

PurposeThis study investigates the effects of temporary tissue expanders (TTEs) on the dose distributions in breast cancer radiotherapy treatments under a variety of conditions.MethodsUsing EBT2 radiochromic film, both electron and photon beam dose distribution measurements were made for different phantoms, and beam geometries. This was done to establish a more comprehensive understanding of the implant's perturbation effects under a wider variety of conditions.ResultsThe magnetic disk present in a tissue expander causes a dose reduction of approximately 20% in a photon tangent treatment and 56% in electron boost fields immediately downstream of the implant. The effects of the silicon elastomer are also much more apparent in an electron beam than a photon beam.ConclusionsEvidently, each component of the TTE attenuates the radiation beam to different degrees. This study has demonstrated that the accuracy of photon and electron treatments of post-mastectomy patients is influenced by the presence of a tissue expander for various beam orientations. The impact of TTEs on dose distributions establishes the importance of an accurately modelled high-density implant in the treatment planning system for post-mastectomy patients.     

Implementation of a dose gradient method into optimization of dose distribution in prostate cancer 3D-CRT plans

AimThe aim of this work is to present a method of beam weight and wedge angle optimization for patients with prostate cancer.Background3D-CRT is usually realized with forward planning based on a trial and error method. Several authors have published a few methods of beam weight optimization applicable to the 3D-CRT. Still, none on these methods is in common use.Materials and methodsOptimization is based on the assumption that the best plan is achieved if dose gradient at ICRU point is equal to zero. Our optimization algorithm requires beam quality index, depth of maximum dose, profiles of wedged fields and maximum dose to femoral heads. The method was tested for 10 patients with prostate cancer, treated with the 3-field technique. Optimized plans were compared with plans prepared by 12 experienced planners. Dose standard deviation in target volume, and minimum and maximum doses were analyzed.ResultsThe quality of plans obtained with the proposed optimization algorithms was comparable to that prepared by experienced planners. Mean difference in target dose standard deviation was 0.1% in favor of the plans prepared by planners for optimization of beam weights and wedge angles. Introducing a correction factor for patient body outline for dose gradient at ICRU point improved dose distribution homogeneity. On average, a 0.1% lower standard deviation was achieved with the optimization algorithm. No significant difference in mean dose–volume histogram for the rectum was observed.ConclusionsOptimization shortens very much time planning. The average planning time was 5 min and less than a minute for forward and computer optimization, respectively.     

Characterization of irregular electron beam for boost dose after whole breast irradiation

Irradiating a tumor bed with boost dose after whole breast irradiation helps reducing the probability of local recurrence. However, the success of electron beam treatment with a small area aiming to cover a superficial lesion is a dual challenge as it requires an adequate dosimetry beside a double check for dose coverage with an estimation of various combined uncertainty of tumor location and losing lateral electron equilibrium within small field dimensions.Aim of workthis work aims to measure the electron beam fluence within different field dimensions and the deviation from measurement performed in standard square electron applicator beam flatness and symmetry, then to calculate the average range of the correction factor required to overcome the loss of lateral electron equilibrium.Material and methodthe electron beam used in this work generated from the linear accelerator model ELEKTA Precise and dosimetry system used were a pair of PTW Pin Point ion chambers for electron beam dosimetry at standard conditions and assessment of beam quality at a reference depth of measurement, with an automatic water phantom, then a Roos ion chamber was used for absolute dose measurement, and PTW 2Darray to investigate the beam fluence of four applicators 6, 10, 14 and 20 cm² and 4 rectangular cutouts 6 × 14, 8 × 14, 6 × 17 and 8 × 17 cm², the second part was clinical application which was performed in a precise treatment planning system and examined boost dose after whole breast irradiation.Resultsrevealed that lower energy (6MeV and 8MeV) showed the loss of lateral electron equilibrium and deviation from measurements of a standard applicator more than the high energy (15 MeV) which indicated that the treatment of superficial dose with 6MeV required higher monitor unit to allow for the loss of lateral electron equilibrium and higher margin as well.     

A patient immobilization device for prone breast radiotherapy: Dosimetric effects and inclusion in the treatment planning system

PurposeTo assess the dosimetric impact of a patient positioning device for prone breast radiotherapy and assess the accuracy of a treatment planning system (TPS) in predicting this impact.MethodsBeam attenuation and build-up dose perturbations, quantified by ionization chamber and radiochromic film dosimetry, were evaluated for 3 components of the patient positioning device: the carbon fiber baseplate, the support cushions and the support wedge for the contralateral breast. Dose calculations were performed using the XVMC dose engine implemented in the Monaco TPS. All components were included during planning CT acquisition.ResultsBeam attenuation amounted to 7.57% (6 MV) and 5.33% (15 MV) for beams obliquely intersecting the couchtop–baseplate combination. Beams traversing large sections of the support wedge were attenuated by 12.28% (6 MV) and 9.37% (15 MV). For the support cushion foam, beam attenuation remained limited to 0.11% (6 MV) and 0.08% (15 MV) per centimeter thickness. A substantial loss of dose build-up was detected when irradiating through any of the investigated components. TPS dose calculations accurately predicted beam attenuation by the baseplate and support wedge. A manual density overwrite was needed to model attenuation by the support cushion foam. TPS dose calculations in build-up regions differed considerably from measurements for both open beams and beams traversing the device components.ConclusionsIrradiating through the components of the positioning device resulted in a considerable degradation of skin sparing. Inclusion of the device components in the treatment planning CT allowed to accurately model the most important attenuation effect, but failed to accurately predict build-up doses.     

Precision IORT – Image guided intraoperative radiation therapy (igIORT) using online treatment planning including tissue heterogeneity correction

BackgroundTo the present date, IORT has been eye and hand guided without treatment planning and tissue heterogeneity correction. This limits the precision of the application and the precise documentation of the location and the deposited dose in the tissue. Here we present a set-up where we use image guidance by intraoperative cone beam computed tomography (CBCT) for precise online Monte Carlo treatment planning including tissue heterogeneity correction.Materials and methodsAn IORT was performed during balloon kyphoplasty using a dedicated Needle Applicator. An intraoperative CBCT was registered with a pre-op CT. Treatment planning was performed in Radiance using a hybrid Monte Carlo algorithm simulating dose in homogeneous (MCwater) and heterogeneous medium (MChet). Dose distributions on CBCT and pre-op CT were compared with each other. Spinal cord and the metastasis doses were evaluated.ResultsThe MCwater calculations showed a spherical dose distribution as expected. The minimum target dose for the MChet simulations on pre-op CT was increased by 40% while the maximum spinal cord dose was decreased by 35%. Due to the artefacts on the CBCT the comparison between MChet simulations on CBCT and pre-op CT showed differences up to 50% in dose.ConclusionsigIORT and online treatment planning improves the accuracy of IORT. However, the current set-up is limited by CT artefacts. Fusing an intraoperative CBCT with a pre-op CT allows the combination of an accurate dose calculation with the knowledge of the correct source/applicator position. This method can be also used for pre-operative treatment planning followed by image guided surgery.     

Dose distribution comparison in volumetric-modulated arc therapy plans for head and neck cancers with and without an external body contour extended technique

AimThis study compared volumetric-modulated arc therapy (VMAT) plans for head and neck cancers with and without an external body contour extended technique (EBCT).BackgroundDose calculation algorisms for VMAT have limitations in the buildup region.Materials and methodsThree VMAT plans were enrolled, with one case having a metal artifact from an artificial tooth. The proper dose was calculated using Eclipse version 11.0. The body contours were extended 2 cm outward from the skin surface in three-dimensional space, and the dose was recalculated with an anisotropic analytical algorithm (AAA) and Acuros XB (AXB). Monitor units (MUs) were set, and the dose distributions in the planning target volume (PTV), clinical target volume, and organ at risk (OAR) and conformity index (CI) with and without an EBCT were compared. The influence of a metal artifact outside of the thermoplastic head mask was also compared.ResultsThe coverage of PTV by the 95% dose line near the patient’s skin was increased drastically by using an EBCT. Plan renormalization had a negligible impact on MUs and doses delivered to OARs. CI of PTV with a 6-MV photon beam was closer to 1 than that with a 10-MV photon beam when both AAA and AXB were used in all cases. Metal artifacts outside the head mask had no effect on dose distribution.ConclusionsAn EBCT is needed to estimate the proper dose at object volumes near the patient’s skin and can improve the accuracy of the calculated dose at target volumes.     

Application of Monte Carlo calculations for validation of a treatment planning system in high dose rate brachytherapy

AimThe accuracy of treatment planning systems is of vital importance in treatment outcomes in brachytherapy. In the current study the accuracy of dose calculations of a high dose rate (HDR) brachytherapy treatment planning system (TPS) was validated using the Monte Carlo method.Materials and methodsThree ⁶⁰Co sources of the GZP6 afterloading brachytherapy system were modelled using MCNP4C Monte Carlo (MC) code. The dose distribution around all the sources was calculated by MC and a dedicated treatment planning system. The results of both methods were compared.ResultsThere was good agreement (<2%) between TPS and MC calculated dose distributions except at a point near the sources (<1 cm) and beyond the tip of the sources.ConclusionsOur study confirmed the accuracy of TPS calculated dose distributions for clinical use in HDR brachytherapy.     

Towards MR-guided electron therapy: Measurement and simulation of clinical electron beams in magnetic fields

PurposeIn the current era of MRI-linac radiotherapy, dose optimization with arbitrary dose distributions is a reality. For the first time, we present new and targeted experiments and modeling to aid in evaluating the potential dose improvements offered with an electron beam mode during MRI-linac radiotherapy.MethodsSmall collimated (1 cm diameter and 1.5 × 1.5 cm² square) electron beams (6, 12 and 20 MeV) from a clinical linear accelerator (Varian Clinac 2100C) are incident perpendicular and parallel to the strong and localized magnetic fields (0–0.7 T) generated by a permanent magnet device. Gafchromic EBT3 film is placed inside a slab phantom to measure two-dimensional dose distributions. A benchmarked and comprehensive Monte Carlo model (Geant4) is established to directly compare with experiments.ResultsWith perpendicular fields a 5% narrowing of the beam FWHM and a 10 mm reduction in the 15% isodose penetration is seen for the 20 MeV beam. In the inline setup the penumbral width is reduced by up to 20%, and a local central dose enhancement of 100% is observed. Monte Carlo simulations are in agreement with the measured dose distributions (2% or 2 mm).ConclusionA new range of experiments have been performed to offer insight into how an electron beam mode could offer additional choices in MRI-linac radiotherapy. The work extends on historic studies to bring a successful unified experimental and Monte Carlo modeling approach for studying small field electron beam dosimetry inside magnetic fields. The results suggest further work, particularly on the inline magnetic field scenario.     

Integration of the M6 Cyberknife in the Moderato Monte Carlo platform and prediction of beam parameters using machine learning

PurposeThis work describes the integration of the M6 Cyberknife in the Moderato Monte Carlo platform, and introduces a machine learning method to accelerate the modelling of a linac.MethodsThe MLC-equipped M6 Cyberknife was modelled and integrated in Moderato, our in-house platform offering independent verification of radiotherapy dose distributions. The model was validated by comparing TPS dose distributions with Moderato and by film measurements. Using this model, a machine learning algorithm was trained to find electron beam parameters for other M6 devices, by simulating dose curves with varying spot size and energy. The algorithm was optimized using cross-validation and tested with measurements from other institutions equipped with a M6 Cyberknife.ResultsOptimal agreement in the Monte Carlo model was reached for a monoenergetic electron beam of 6.75 MeV with Gaussian spatial distribution of 2.4 mm FWHM. Clinical plan dose distributions from Moderato agreed within 2% with the TPS, and film measurements confirmed the accuracy of the model. Cross-validation of the prediction algorithm produced mean absolute errors of 0.1 MeV and 0.3 mm for beam energy and spot size respectively. Prediction-based simulated dose curves for other centres agreed within 3% with measurements, except for one device where differences up to 6% were detected.ConclusionsThe M6 Cyberknife was integrated in Moderato and validated through dose re-calculations and film measurements. The prediction algorithm was successfully applied to obtain electron beam parameters for other M6 devices. This method would prove useful to speed up modelling of new machines in Monte Carlo systems.     

On the use of a novel Ferrous Xylenol-orange gelatin dosimeter for HDR brachytherapy commissioning and quality assurance testing

PurposeTo evaluate a commercially available Ferrous-Xylenol Orange-Gel (FXG) dosimeter (TrueView™) coupled with Optical-Computed Tomography (OCT) read out, for 3D dose verification in an Ir-192 superficial brachytherapy application.MethodsTwo identical polyethylene containers filled with gel from the same batch were used. One was irradiated with an 18 MeV electron field to examine the dose-response linearity and obtain a calibration curve. A flap surface applicator was attached to the other to simulate treatment of a skin lesion. The dose distribution in the experimental set up was calculated with the TG-43 and the model based dose calculation (MBCA) algorithms of a commercial treatment planning system (TPS), as well as Monte Carlo (MC) simulation using the MCNP code. Measured and calculated dose distributions were spatially registered and compared.ResultsApart from a region close to the container’s neck, where gel measurements exhibited an over-response relative to MC calculations (probably due to stray light perturbation), an excellent agreement was observed between measurements and simulations. More than 97% of points within the 10% isodose line (80 cGy) met the gamma index criteria established from uncertainty analysis (5%/2 mm). The corresponding passing rates for the comparison of experiment to calculations using the TG-43 and MBDCA options of the TPS were 57% and 92%, respectively.ConclusionTrueView™ is suitable for the quality assurance of demanding radiotherapy applications. Experimental results of this work confirm the advantage of the studied MBDCA over TG-43, expected from the improved account of scatter radiation in the treatment geometry.     

Comparison of the RayStation photon Monte Carlo dose calculation algorithm against measured data under homogeneous and heterogeneous irradiation geometries

PurposeThis work compares Monte Carlo dose calculations performed using the RayStation treatment planning system against data measured on a Varian Truebeam linear accelerator with 6 MV and 10 MV FFF photon beams.MethodsThe dosimetric performance of the RayStation Monte Carlo calculations was evaluated in a variety of irradiation geometries employing homogeneous and heterogeneous phantoms. Profile and depth dose comparisons against measurement were carried out in relative mode using the gamma index as a quantitative measure of similarity within the central high dose regions.ResultsThe results demonstrate that the treatment planning system dose calculation engine agrees with measurement to within 2%/1 mm for more than 95% of the data points in the high dose regions for all test cases. A systematic underestimation was observed at the tail of the profile penumbra and out of field, with mean differences generally <0.5 mm or 1% of curve dose maximum respectively. Out of field agreement varied between evaluated beam models.ConclusionsThe RayStation implementation of photon Monte Carlo dose calculations show good agreement with measured data for the range of scenarios considered in this work and is deemed sufficiently accurate for introduction into clinical use.     

A machine learning-based framework for delivery error prediction in proton pencil beam scanning using irradiation log-files

PurposeThis study aims to investigate the use of machine learning models for delivery error prediction in proton pencil beam scanning (PBS) delivery.MethodsA dataset of planned and delivered PBS spot parameters was generated from a set of 20 prostate patient treatments. Planned spot parameters (spot position, MU and energy) were extracted from the treatment planning system (TPS) for each beam. Delivered spot parameters were extracted from irradiation log-files for each beam delivery following treatment. The dataset was used as a training dataset for three machine learning models which were trained to predict delivered spot parameters based on planned parameters. K-fold cross validation was employed for hyper-parameter tuning and model selection where the mean absolute error (MAE) was used as the model evaluation metric. The model with lowest MAE was then selected to generate a predicted dose distribution for a test prostate patient within a commercial TPS.ResultsAnalysis of the spot position delivery error between planned and delivered values resulted in standard deviations of 0.39 mm and 0.44 mm for x and y spot positions respectively. Prediction error standard deviation values of spot positions using the selected model were 0.22 mm and 0.11 mm for x and y spot positions respectively. Finally, a three-way comparison of dose distributions and DVH values for select OARs indicates that the random-forest-predicted dose distribution within the test prostate patient was in closer agreement to the delivered dose distribution than the planned distribution.ConclusionsPBS delivery error can be accurately predicted using machine learning techniques.     

Dose gradient based algorithm for beam weights selection in 3D-CRT plans

AimIn this work we test the usage of dose gradient based algorithm for the selection of beam weights in 3D-CRT plans for different cancer locations. Our algorithm is easy to implement for three fields technique with wedges defined by planner.Background3D-CRT is usually realized with forward planning which is quite time consuming. Several authors published a few methods of beams weights optimization applicable to the 3D-CRT.Materials and methodsOptimization is based on an assumption that the best plan is achieved if dose gradient at ICRU point is equal to zero. Method was tested for 120 patients, treated in our clinic in 2011-2012, with different cancer locations. For each patient, three fields conformal plan (6 MV and 15 MV X-ray) with the same geometry as proposed by experienced planners was prepared. We compared dose distributions achieved with the proposed method and those prepared by experienced planners. The homogeneity of dose distributions was compared in terms of STD and near minimum and near maximum doses in the PTV.ResultsMean difference of STD obtained by the proposed algorithm and by planners was 0.1%: 0.1% for prostate cancer, 0.3% for lung cancer, −0.1% for esophagus cancer, 0.1% for rectum cancer, −0.1% for gynecology cancer, −0.1% for stomach cancer.ConclusionsApplying the proposed algorithm leads to obtain the similar dose distribution homogeneity in the PTV as these achieved by planners and therefore can serve as a support in creating 3D-CRT plans. It is also simple to use and can significantly speed up the treatment planning process.     

Novel Monte Carlo dose calculation algorithm for robotic radiosurgery with multi leaf collimator: Dosimetric evaluation

PurposeAt introduction in 2014, dose calculation for the first MLC on a robotic SRS/SBRT platform was limited to a correction-based Finite-Size Pencil Beam (FSPB) algorithm. We report on the dosimetric accuracy of a novel Monte Carlo (MC) dose calculation algorithm for this MLC, included in the Precision™ treatment planning system.MethodsA phantom was built of one slab (5.0 cm) of lung-equivalent material (0.09…0.29 g/cc) enclosed by 3.5 cm (above) and 5 cm (below) slabs of solid water (1.045 g/cc). This was irradiated using rectangular (15.4 × 15.4 mm² to 53.8 × 53.7 mm²) and two irregular MLC-fields. Radiochromic film (EBT3) was positioned perpendicular to the slabs and parallel to the beam. Calculated dose distributions were compared to film measurements using line scans and 2D gamma analysis.ResultsMeasured and MC calculated percent depth dose curves showed a characteristic dose drop within the low-density region, which was not correctly reproduced by FSPB. Superior average gamma pass rates (2%/1 mm) were found for MC (91.2 ± 1.5%) compared to FSPB (55.4 ± 2.7%). However, MC calculations exhibited localized anomalies at mass density transitions around 0.15 g/cc, which were traced to a simplification in electron transport. Absence of these anomalies was confirmed in a modified build of the MC engine, which increased gamma pass rates to 96.6 ± 1.2%.ConclusionsThe novel MC algorithm greatly improves dosimetric accuracy in heterogeneous tissue, potentially expanding the clinical use of robotic radiosurgery with MLC. In-depth, independent validation is paramount to identify and reduce the residual uncertainties in any software solution.     

On pathlength and energy straggling of megavoltage electrons slowing down

Purposeto elucidate the effects of multiple scattering and energy-loss straggling on electron beams slowing down in materials.MethodsEGSnrc Monte Carlo simulations are done using a purpose-written user-code.ResultsPlots are presented of the primary electron’s energy as a function of pathlength for 20 MeV electrons incident on water and tantalum as are plots of the overall distribution of pathlengths as the 20 MeV electrons slow down under various Monte Carlo scenarios in water and tantalum. The distributions range from 1 % to 135 % of the CSDA range in water and from 1 % to 186 % in tantalum. The effects of energy-loss straggling on energy spectra at depth and electron fluence at depth are also presented.ConclusionsThe role of energy-loss straggling and multiple scattering are shown to play a significant role in the range straggling which determines the dose fall-off region in electron beam dose vs depth curves and a significant role in the energy distributions as a function of depth.     

Combined PET/CT for IMRT treatment planning of NSCLC: Contrast-enhanced CT images for Monte Carlo dose calculation

PurposeCombined PET/CT imaging has been proposed as an integral part of radiotherapy treatment planning (TP). Contrast-enhanced CT (ceCT) images are frequently acquired as part of the PET/CT examination to support target delineation. The aim of this dosimetric planning study was to investigate the error introduced by using a ceCT for intensity modulated radiotherapy (IMRT) TP with Monte Carlo dose calculation for non-small cell lung cancer (NSCLC).Material and methodsNine patients with NSCLC prior to chemo-RT were included in this retrospective study. For each patient non-enhanced, low-dose CT (neCT), ceCT and [¹⁸F]-FDG-PET emission data were acquired within a single examination. Manual contouring and TP were performed on the ceCT. An additional set of independent target volumes was auto-segmented in PET images. Dose distributions were recalculated on the neCT. Differences in dosimetric parameters were evaluated.ResultsDose differences in PTV and lungs were small for all patients. The maximum difference in all PTVs when using ceCT images for dose calculation was ?2.1%, whereas the mean difference was less than ?1.7%. Maximum differences in the lungs ranged from ?1.8% to 2.1% (mean: ?0.1%). In four patients an underestimation of the maximum spinal cord dose between 2% and 3.2% was observed, but treatment plans remained clinically acceptable.ConclusionsMonte Carlo based IMRT planning for NSCLC patients using ceCT allows for correct dose calculation. A direct comparison to neCT-based treatment plans revealed only small dose differences. Therefore, ceCT-based TP is clinically safe as long as the maximum acceptable dose to organs at risk is not approached.     

eXaSkin: A novel high-density bolus for 6MV X-rays radiotherapy

PurposeTo evaluate eXaSkin, a novel high-density bolus alternative to commercial tissue-equivalent Superflab, for 6MV photon-beam radiotherapy.Materials and methodsWe delivered a 10 × 10 cm² open field at 90° and head-and-neck clinical plan, generated with the volumetric modulated arc therapy (VMAT) technique, to an anthropomorphic phantom in three scenarios: with no bolus on the phantom’s surface, with Superflab, and with eXaSkin. In each scenario, we measured dose to a central planning target volume (PTV) in the nasopharynx region with an ionization chamber, and we measured dose to the skin, at three different positions within the vicinity of a neck lymph node PTV, with MOSkin™, a semiconductor dosimeter. Measurements were compared against calculations with the treatment planning system (TPS).ResultsFor the static field, MOSkin results underneath the eXaSkin were in agreement with calculations to within 1.22%; for VMAT, to within 5.68%. Underneath Superflab, those values were 3.36% and 11.66%. The inferior agreement can be explained by suboptimal adherence of Superflab to the phantom’s surface as well as difficulties in accurately reproducing its placement between imaging and treatment session. In all scenarios, dose measured at the central target agreed to within 1% with calculations.ConclusionseXaSkin was shown to have superior adaptation to the phantom’s surface, producing minimal air gaps between the skin surface and bolus, allowing for accurate positioning and reproducibility of set-up conditions. eXaSkin with its high density material provides sufficient build-up to achieve full skin dose with less material thickness than Superflab.     

Effects of shielding on pelvic and abdominal IORT dose distributions

PurposeTo study the impact of shielding elements in the proximity of Intra-Operative Radiation Therapy (IORT) irradiation fields, and to generate graphical and quantitative information to assist radiation oncologists in the design of optimal shielding during pelvic and abdominal IORT.MethodAn IORT system was modeled with BEAMnrc and EGS++ Monte Carlo codes. The model was validated in reference conditions by gamma index analysis against an experimental data set of different beam energies, applicator diameters, and bevel angles. The reliability of the IORT model was further tested considering shielding layers inserted in the radiation beam. Further simulations were performed introducing a bone-like layer embedded in the water phantom. The dose distributions were calculated as 3D dose maps.ResultsThe analysis of the resulting 2D dose maps parallel to the clinical axis shows that the bevel angle of the applicator and its position relative to the shielding have a major influence on the dose distribution. When insufficient shielding is used, a hotspot nearby the shield appears near the surface. At greater depths, lateral scatter limits the dose reduction attainable with shielding, although the presence of bone-like structures in the phantom reduces the impact of this effect.ConclusionsDose distributions in shielded IORT procedures are affected by distinct contributions when considering the regions near the shielding and deeper in tissue: insufficient shielding may lead to residual dose and hotspots, and the scattering effects may enlarge the beam in depth. These effects must be carefully considered when planning an IORT treatment with shielding.