Modelling of an Orthovoltage X-ray Therapy Unit with the EGSnrc Monte Carlo Package

Clinical treatment planning for kilovoltage radiotherapy using EGSnrc and Python

Kilovoltage radiotherapy dose calculations are generally performed with manual point dose calculations based on water dosimetry. Tissue heterogeneities, irregular surfaces, and introduction of lead cutouts for treatment are either not taken into account or crudely approximated in manual calculations. Full Monte Carlo (MC) simulations can account for these limitations but require a validated treatment unit model, accurately segmented patient tissues and a treatment planning interface (TPI) to facilitate the simulation setup and result analysis. EGSnrc was used in this work to create a model of Xstrahl kilovoltage unit extending the range of energies, applicators, and validation parameters previously published. The novel functionality of the Python-based framework developed in this work allowed beam modification using custom lead cutouts and shields, commonly present in kilovoltage treatments, as well as absolute dose normalization using the output of the unit. 3D user-friendly planning interface of the developed framework facilitated non-co-planar beam setups for CT phantom MC simulations in DOSXYZnrc. The MC models of 49 clinical beams showed good agreement with measured and reference data, to within 2% for percentage depth dose curves, 4% for beam profiles at various depths, 2% for backscatter factors, 0.5 mm of absorber material for half-value layers, and 3% for output factors. End-to-end testing of the framework using custom lead cutouts resulted in good agreement to within 3% of absolute dose distribution between simulations and EBT3 GafChromic film measurements. Gamma analysis demonstrated poor agreement at the field edges which was attributed to the limitations of simulating smooth cutout shapes. Dose simulated in a heterogeneous phantom agreed to within 7% with measured values converted using the ratio of mass energy absorption coefficients of appropriate tissues and air.

Bone and mucosal dosimetry in skin radiation therapy: a Monte Carlo study using kilovoltage photon and megavoltage electron beams

This study examines variations of bone and mucosal doses with variable soft tissue and bone thicknesses, mimicking the oral or nasal cavity in skin radiation therapy. Monte Carlo simulations (EGSnrc-based codes) using the clinical kilovoltage (kVp) photon and megavoltage (MeV) electron beams, and the pencil-beam algorithm (Pinnacle(3) treatment planning system) using the MeV electron beams were performed in dose calculations. Phase-space files for the 105 and 220 kVp beams (Gulmay D3225 x-ray machine), and the 4 and 6 MeV electron beams (Varian 21 EX linear accelerator) with a field size of 5 cm diameter were generated using the BEAMnrc code, and verified using measurements. Inhomogeneous phantoms containing uniform water, bone and air layers were irradiated by the kVp photon and MeV electron beams. Relative depth, bone and mucosal doses were calculated for the uniform water and bone layers which were varied in thickness in the ranges of 0.5-2 cm and 0.2-1 cm. A uniform water layer of bolus with thickness equal to the depth of maximum dose (d(max)) of the electron beams (0.7 cm for 4 MeV and 1.5 cm for 6 MeV) was added on top of the phantom to ensure that the maximum dose was at the phantom surface. From our Monte Carlo results, the 4 and 6 MeV electron beams were found to produce insignificant bone and mucosal dose (<1%), when the uniform water layer at the phantom surface was thicker than 1.5 cm. When considering the 0.5 cm thin uniform water and bone layers, the 4 MeV electron beam deposited less bone and mucosal dose than the 6 MeV beam. Moreover, it was found that the 105 kVp beam produced more than twice the dose to bone than the 220 kVp beam when the uniform water thickness at the phantom surface was small (0.5 cm). However, the difference in bone dose enhancement between the 105 and 220 kVp beams became smaller when the thicknesses of the uniform water and bone layers in the phantom increased. Dose in the second bone layer interfacing with air was found to be higher for the 220 kVp beam than that of the 105 kVp beam, when the bone thickness was 1 cm. In this study, dose deviations of bone and mucosal layers of 18% and 17% were found between our results from Monte Carlo simulation and the pencil-beam algorithm, which overestimated the doses. Relative depth, bone and mucosal doses were studied by varying the beam nature, beam energy and thicknesses of the bone and uniform water using an inhomogeneous phantom to model the oral or nasal cavity. While the dose distribution in the pharynx region is unavailable due to the lack of a commercial treatment planning system commissioned for kVp beam planning in skin radiation therapy, our study provided an essential insight into the radiation staff to justify and estimate bone and mucosal dose.

Kilovoltage x-ray dosimetry--an experimental comparison between different dosimetry protocols

Kilovoltage dosimetry protocols by the IAEA (TRS-277 and TRS-398), DIN (6809), IPEMB (with addendum), AAPM (TG-61) and NCS (report 10) were compared experimentally in four clinical beams. The beams had acceleration potentials of 30, 80, 120 and 200 kV, with half-value layers ranging from 0.6 mm Al to 1 mm Cu. Dosimetric measurements were performed and data were collected under reference conditions as stipulated within each separate protocol under investigation. The Monte Carlo method was used to derive backscatter factors for the actual x-ray machine. In general, the agreement of the dosimetric data at the surface of a full-scatter water phantom obtained using the guidelines of the various protocols was fairly good, i.e. within 1-2%. However, the in-air calibration method using the IPEMB and AAPM TG-61 protocols yielded an absorbed dose about 7% lower than the IAEA TRS-398 protocol in the 120 kV beam. By replacing the backscatter factors given in the protocols with Monte Carlo calculated backscatter factors, the convergence between the protocols improved (within 4%). The internal consistency obtained for protocols supporting more than one geometry for dosimetry under reference conditions was better than 0.2% for the DIN protocol (120 kV beam), 2-3% for the AAPM TG-61 (120 and 200 kV beams) and about 2% for the IPEMB protocol (200 kV beam). The present study shows that the current-supported dosimetry protocols in the kilovoltage range were in fairly good agreement, and there were only a few exceptions of clinical significance.