Search

Background and purpose: Monte Carlo simulations as well as analytical computations of proton transport in material media require accurate values of multiple Coulomb scattering (MCS) angles. High-quality experimental data on MCS angles in the energy range for proton therapy are, however, sparse. In this work, MCS modeling in proton transport was evaluated employing an experimental method to measure these angles on a medical proton beamline in clinically relevant materials. Results are compared to Monte Carlo simulations and analytical models. Materials and methods: Aluminum, brass, and lucite (PMMA) scatterers of clinically relevant thicknesses were irradiated with protons at 100, 160, and 220 MeV. Resulting spatial distributions of individual pencil beams were measured with a scintillating screen. The MCS angles were determined by deconvolution and a virtual point source approach. Results were compared to those obtained with the Monte Carlo codes PENH, TOPAS, and RayStation Monte Carlo, as well as the analytical models RayStation Pencil Beam Algorithm and the Molière/Fano/Hanson variant of the Molière theory. Results: Experimental data obtained with the presented methodology agree with previously published results within 6%, with an average deviation of 3%. The combined average uncertainty of the experimental data yielded 1.8%, while the combined maximum uncertainty was below 4%. The obtained Monte Carlo results for PENH, TOPAS, and RayStation deviate on average for all considered energies, materials and thicknesses, by 2.5%, 3.4%, and 2.8% from the experimental data, respectively. For the analytical models, the average deviations amount to 4.5% and 2.9% for the RayStation Pencil Beam Algorithm and the Molière/Fano/Hanson model, respectively. Conclusion: The experimental method developed for the present work allowed to measure MCS angles in clinical proton facilities with good accuracy. The presented method permits to extend the database on experimental MCS angles which is rather limited. This work further provides benchmark data for lucite in thicknesses relevant for clinical applications. The data may serve to validate dose engines of treatment planning systems and secondary dose check software. The Monte Carlo and analytical algorithms studied are capable of reproducing MCS data within the required accuracy for clinical applications.

Poster-Abstract

Background and purpose: PENH is a recently coded module for simulation of proton transport in conjunction with the Monte Carlo code PENELOPE. PENELOPE applies class II simulation to all type of interactions, in particular, to elastic collisions. PENH uses calculated differential cross sections for proton elastic collisions that include electron screening effects as well as nuclear structure effects. Proton-induced nuclear reactions are simulated from information in the ENDF-6 database or from alternative nuclear databases in ENDF format. The purpose of this work is to benchmark this module by simulating absorbed dose distributions from a single finite spot size proton pencil beam in water. Materials and methods: Monte Carlo simulations with PENH are compared with simulation results from TOPAS Monte Carlo (v3.1p2) and RayStation Monte Carlo (v6). Different beam models are examined in terms of mean energy and energy spread to match the measured profiles. The phase-space file is derived from experimental measurements. Simulated absorbed dose distributions are compared to experimental data obtained with the ionization chamber array MatriXX 2D detector (IBA Dosimetry) in a water tank. The experiments were conducted with a clinical IBA pencil beam scanning dedicated nozzle. In all simulations a Fermi-Eyges phase-space representation of a single finite spot size proton pencil beam is used. Results: In general, there is a good agreement between simulated results and experimental data up to a distance of 3 cm from the central axis. In the core region (region where the dose is more than 10% of the maximum dose) PENH shows, overall, the smallest deviations from experimental data, with the largest radial rms (root mean square) smaller than 0.2. The results achieved by TOPAS and RayStation in that region are very close to those of PENH. For the halo region, that is the area of the dose distribution outside the core region reaching down to 0.01% of the maximum intensity, the largest rms achieved by TOPAS is always smaller than 0.5, yielding better results than the rest of the codes.

Background The availability of photon and electron spectra in digital form from current accelerators and Monte Carlo (MC) systems is scarce, and one of the packages widely used refers to linacs with a reduced clinical use nowadays. Such spectra are mainly intended for the MC calculation of detector-related quantities in conventional broad beams, where the use of detailed phase-space files (PSFs) is less critical than for MC-based treatment planning applications, but unlike PSFs, spectra can easily be transferred to other computer systems and users. Methods A set of spectra for a range of Varian linacs has been calculated using the PENELOPE/PRIMO MC system. They have been extracted from PSFs tallied for field sizes of 10 cm × 10 cm and 15 cm × 15 cm for photon and electron beams, respectively. The influence of the spectral bin width and of the beam central axis region used to extract the spectra have been analyzed. Results Spectra have been compared to those by other authors showing good agreement with those obtained using the, now superseded, EGS4/BEAM MC code, but significant differences with the most widely used photon data set. Other spectra, particularly for electron beams, have not been published previously for the machines simulated in this work. The influence of the bin width on the spectrum mean energy for 6 and 10 MV beams has been found to be negligible. The size of the region used to extract the spectra yields differences of up to 40% for the mean energies in 10 MV beams, but the maximum difference for TPR 20,10 values derived from depth-dose distributions does not exceed 2% relative to those obtained using the PSFs. This corresponds to kQ differences below 0.2% for a typical Farmer-type chamber, considered to be negligible for reference dosimetry. Different configurations for using electron spectra have been compared for 6 MeV beams, concluding that the geometry used for tallying the PSFs used to extract the spectra must be accounted for in subsequent calculations using the spectra as a source. Conclusions An up-to-date set of consistent spectra for Varian accelerators suitable for the calculation of detector-related quantities in conventional broad beams has been developed and made available in digital form. Electronic supplementary material The online version of this article (10.1186/s13014-018-1186-8) contains supplementary material, which is available to authorized users.

Background/Aims: Ruthenium plaques are used for the treatment of ocular tumors. The aim of this work is the comparison between simulated absorbed dose distributions tallied in an anthropomorphic phantom, obtained from ideal homogeneous plaques, and real eye plaques in which the actual heterogeneous distribution of 106Ru was measured. The placement of the plaques with respect to the tumor location was taken into consideration to optimize the effectiveness of the treatment. Methods: The generic CCA and CCB, and the specific CCA1364 and CCB1256 106Ru eye plaques were modeled with the Monte Carlo code PENELOPE. To compare the suitability of each treatment for an anterior, equatorial and posterior tumor location, cumulative dose-volume histograms for the tumors and structures at risk were calculated. Results: Eccentric placements of the plaques, taking into account the inhomogeneities of the emitter map, can substantially reduce the dose delivered to structures at risk while maintaining the prescribed dose at the tumor apex. Conclusions: The emitter map distribution of the plaque and the computerized tomography of the patient used in a Monte Carlo simulation allow an accurate determination of the plaque position with respect to the tumor with the potential to reduce the dose to sensitive structures.

Poster Abstract