Validation and automation of the DYNJAWS component module of the BEAMnrc Monte Carlo code

Comparison of dosimetric characteristics of physical wedge and enhanced dynamic wedge in inhomogeneous medium using Monte Carlo simulations

Background

Widely used physical wedges in clinical radiotherapy lead to beam intensity attenuation as well as the beam hardening effect, which must be considered. Dynamic wedges devised to overcome the physical wedges (PWs) problems result in dosimetry complications due to jaw movement while the beam is on. This study was aimed to investigate the usability of physical wedge data instead of enhanced dynamic wedge due to the enhanced dynamic wedge (EDW) dosimetry measurement hardships of Varian 2100CD in inhomogeneous phantom by Monte Carlo code as a reliable method in radiation dosimetry.

Materials and methods

A PW and EDW-equipped-linac head was simulated using BEAMnrc code. DOSXYZnrc was used for three-dimensional dosimetry calculation in the CIRS phantom.

Results

Based on the isodose curves, EDW generated a less scattered as well as lower penumbra width compared to the PW. The depth dose variations of PWs and EDWs were more in soft tissue than the lung tissue. Beam profiles of PW and EDW indicated good coincidence in all points, except for the heel area.

Conclusion

Results demonstrated that it is possible to apply PW data instead of EDW due to the dosimetry and commissioning hardships caused by EDW in inhomogeneous media.

Monte Carlo commissioning of radiotherapy LINAC-Introducing an improved methodology

Purpose

Monte Carlo (MC) commissioning of medical linear accelerator (LINAC) is a time-consuming process involving a comparison between measured and simulated cross beam/lateral profiles and percentage depth doses (PDDs) for various field sizes. An agreement between these two data sets is sought by trial and error method while varying the incident electron beam parameters, such as electron beam energy or width, etc. This study aims to improve the efficiency of MC commissioning of a LINAC by assessing the feasibility of using a limited number of simulated PDDs.

Materials and methods

Using EGSnrc codes, a Varian Clinac 2100 unit has been commissioned for 6 MV photon beam, and a methodology has been proposed to identify the incident electron beam parameters in a speedier fashion. Impact of voxel size in 3-dimensions and cost functions used for comparison of the measured and simulated data have been investigated along with the role of interpolation.

Results

A voxel size of 1 × 1×0.5 cm³ has been identified as suitable for accurate and fast commissioning of the LIANC. The optimum number of simulated PDDs (required for further interpolation) has been found to be five.

Conclusion

The present study suggests that PDDs alone at times can be insufficient for an unambiguous commissioning process and should be supported by including the lateral beam profiles in the process.

A comparison of Monte Carlo, anisotropic analytical algorithm (AAA) and Acuros XB algorithms in assessing dosimetric perturbations during enhanced dynamic wedged radiotherapy deliveries in heterogeneous media

Background

A comparison of anisotropic analytical algorithm (AAA) and Acuros XB (AXB) dose calculation algorithms with Electron Gamma Shower (EGSnrc) Monte Carlo (MC) for modelling lung and bone heterogeneities encountered during enhanced dynamic wedged (EDWs) radiotherapy dose deliveries was carried out.

Materials and methods

In three heterogenous slab phantoms: water–bone, lung–bone and bone–lung, wedged percentage depth doses with EGSnrc, AAA and AXB algorithms for 6 MV photons for various field sizes (5×5, 10×10 and 20×20 cm2) and EDW angles (15°, 30°, 45° and 60°) have been scored.

Results

For all the scenarios, AAA and AXB results were within ±1% of the MC in the pre-inhomogeneity region. For water–bone AAA and AXB deviated by 6 and 1%, respectively. For lung–bone an underestimation in lung (AAA: 5%, AXB: 2%) and overestimation in bone was observed (AAA: 13%, AXB: 4%). For bone–lung phantom overestimation in bone (AAA: 7%, AXB: 1%), a lung underdosage (AAA: 8%, AXB: 5%) was found. Post bone up to 12% difference in the AAA and MC results was observed as opposed to 6% in case of AXB.

Conclusion

This study demonstrated the limitation of the AAA (in certain scenarios) and accuracy of AXB for dose estimation inside and around lung and bone inhomogeneities. The dose perturbation effects were found to be slightly dependent on the field size with no obvious EDW dependence.

Soft tissue and water substitutes for megavoltage photon beams: An EGSnrc-based evaluation

In this work, soft‐tissue equivalence of water, polystyrene, PMMA and water equivalence of polystyrene, and PMMA has been assessed for multiple megavoltage photon beams and field sizes. EGSnrc based Monte Carlo (MC) codes, BEAMnrc and DOSXYZnrc are used for the linac head modeling and the phantom dose calculations, respectively. Percentage depth doses (PDDs) are scored for two field sizes (5×5 cm2, 10×10 cm2) and photon energies (6 MV and 10 MV) in water, polystyrene, PMMA, and soft tissue. The comparisons of PDDs show that soft‐tissue equivalence of various materials varies with the depth in the phantom, field size, and photon energy. Water and PMMA are found to be the closest soft‐tissue and water substitutes, respectively. Soft‐tissue and water equivalence of dosimetry materials need to be evaluated for a range of photon energies and field sizes before their application in complex radiation beams.

PACS numbers: 87.55.Gh, 87.55.K‐

Jaw position uncertainty and adjacent fields in breast cancer radiotherapy

Locoregional treatment of breast cancer involves adjacent, half blocked fields matched at isocenter. The objective of this work is to study the dosimetric effects of the uncertainties in jaw positioning for such a case, and how a treatment planning protocol including adjacent field overlap of 1 mm affects the dose distribution. A representative treatment plan, involving 6 and 15 photon beams, for a patient treated at our hospital is chosen. Monte Carlo method (EGSnrc/BEAMnrc) is used to simulate the treatment. Uncertainties in jaw positioning of ± 1 mm are addressed, which implies extremes in reality of 2 mm field gap/overlap when planning adjacent fields without overlap and 1 mm gap or 3 mm overlap for a planning protocol with 1 mm overlap. Dosimetric parameters for PTV, lung and body are analyzed. Treatment planning protocol with 1 mm overlap of the adjacent fields does not considerably counteract possible underdosage of the target in the case studied. PTV-V95% is for example reduced from 95% for perfectly aligned fields to 90% and 91% for 2 mm and 1 mm gap, respectively. However, the risk of overdosage in PTV and in healthy soft tissue is increased when following the protocol with 1 mm overlap. A 3 mm overlap compared to 2 mm overlap results in an increase in maximum dose to PTV, PTV-D2%, from 113% to 121%. V120% for 'Body-PTV' is also increased from 5 cm(3) to 14 cm(3). A treatment planning protocol with 1 mm overlap does not considerably improve the coverage of PTV in the case of erroneous jaw positions causing gap between fields, but increases the overdosage in PTV and doses to healthy tissue, in the case of overlapping fields, for the case investigated.

Monte Carlo simulation of linac irradiation with dynamic wedges

This study aims to simulate the dose distributions of LINAC with dynamic wedges (DWs) under various field sizes and wedge angles by the BEAMnrc code with DYNJAWS component module. These were compared with those calculated by the treatment planning system (TPS) and the measured data. All percentage depth doses (PDDs) were found to be in good agreement between TPS, Monte Carlo (MC) and measurements made in open fields and fields with DWs. For dose profiles, compared with the MC and the measurements, TPS gives reliable results for large field sizes (>10 × 10 cm(2)) but results in significant errors in small field sizes (5 × 5 cm(2)). The entrance surface doses calculated by TPS were found to be significantly overestimated. For depths deeper than 0.5 cm, TPS yields PDDs in agreement with MC simulations.

A comparison of physical and dosimetric properties of lung substitute materials

PURPOSE

The need for an accurate estimate of absorbed doses within and around irradiated thorax tissues necessitates the use of carefully selected materials from which phantoms are constructed. A lung substitute is more difficult to establish mostly due to its low physical density. Although many researchers have used cork as a lung substitute, very little research data address cork's characteristics to determine which type of cork is optimal as a substitute for lung tissue.

METHODS

Natural cork, composition cork, rubber cork, ATOM, RANDO, and a reference lung material (ICRU-44 lung tissue) were investigated to establish comparisons of physical properties. Following the determination of the respective physical properties, the dose distributions from 6 MV photon beams in water/lung substitute/water phantoms were assessed using the Monte Carlo method. Physical and electron densities affecting the dose distributions through lung tissues in different field size conditions were investigated.

RESULTS

The physical properties (physical density, electronic density, and effective atomic number) of the composition cork are the most similar to those of the ICRU-44 lung, and the CT number of the composition cork is very similar to that of humans aged 30-60. PDD of the composition cork and the RANDO phantom are the most comparable to that of ICRU-44 lung in 1 × 1 cm(2) field size due to the combined properties of physical density (PD) and electron density per gram (EDG) of the studied lung materials. PD and EDG affect the lung dose primarily in small field size. The effects of PD are minimal in large fields, having a more rapid lateral electron equilibrium. EDG dominates PDD pattern in lung material when large fields are applied. Combined effects of PD and EDG are nonlinear for all field sizes.

CONCLUSIONS

The composition cork is the preferred lung substitute based on physical and dosimetric properties.