skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Spread-out Bragg peak and monitor units calculation with the Monte Carlo Code MCNPX

Abstract

The aim of this work was to study the dosimetric potential of the Monte Carlo code MCNPX applied to the protontherapy field. For series of clinical configurations a comparison between simulated and experimental data was carried out, using the proton beam line of the MEDICYC isochronous cyclotron installed in the Centre Antoine Lacassagne in Nice. The dosimetric quantities tested were depth-dose distributions, output factors, and monitor units. For each parameter, the simulation reproduced accurately the experiment, which attests the quality of the choices made both in the geometrical description and in the physics parameters for beam definition. These encouraging results enable us today to consider a simplification of quality control measurements in the future. Monitor Units calculation is planned to be carried out with preestablished Monte Carlo simulation data. The measurement, which was until now our main patient dose calibration system, will be progressively replaced by computation based on the MCNPX code. This determination of Monitor Units will be controlled by an independent semi-empirical calculation.

Authors:
; ; ;  [1];  [2];  [3]
  1. Centre Antoine Lacassagne, Cyclotron Biomedical, 227 avenue de la Lanterne, 06200 Nice (France)
  2. (Monaco)
  3. (France)
Publication Date:
OSTI Identifier:
20951058
Resource Type:
Journal Article
Resource Relation:
Journal Name: Medical Physics; Journal Volume: 34; Journal Issue: 2; Other Information: DOI: 10.1118/1.2431473; (c) 2007 American Association of Physicists in Medicine; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
61 RADIATION PROTECTION AND DOSIMETRY; BRAGG CURVE; CALIBRATION; COMPUTERIZED SIMULATION; DEPTH DOSE DISTRIBUTIONS; DOSIMETRY; ISOCHRONOUS CYCLOTRONS; MONTE CARLO METHOD; PATIENTS; PROTON BEAMS; QUALITY CONTROL; RADIATION DOSES; RADIOTHERAPY

Citation Formats

Herault, J., Iborra, N., Serrano, B., Chauvel, P., Centre Hospitalier Princesse Grace, Service de Medecine Nucleaire, avenue Pasteur, 98000 Monaco, and Centre Antoine Lacassagne, Cyclotron Biomedical, 227 avenue de la Lanterne, 06200 Nice. Spread-out Bragg peak and monitor units calculation with the Monte Carlo Code MCNPX. United States: N. p., 2007. Web. doi:10.1118/1.2431473.
Herault, J., Iborra, N., Serrano, B., Chauvel, P., Centre Hospitalier Princesse Grace, Service de Medecine Nucleaire, avenue Pasteur, 98000 Monaco, & Centre Antoine Lacassagne, Cyclotron Biomedical, 227 avenue de la Lanterne, 06200 Nice. Spread-out Bragg peak and monitor units calculation with the Monte Carlo Code MCNPX. United States. doi:10.1118/1.2431473.
Herault, J., Iborra, N., Serrano, B., Chauvel, P., Centre Hospitalier Princesse Grace, Service de Medecine Nucleaire, avenue Pasteur, 98000 Monaco, and Centre Antoine Lacassagne, Cyclotron Biomedical, 227 avenue de la Lanterne, 06200 Nice. Thu . "Spread-out Bragg peak and monitor units calculation with the Monte Carlo Code MCNPX". United States. doi:10.1118/1.2431473.
@article{osti_20951058,
title = {Spread-out Bragg peak and monitor units calculation with the Monte Carlo Code MCNPX},
author = {Herault, J. and Iborra, N. and Serrano, B. and Chauvel, P. and Centre Hospitalier Princesse Grace, Service de Medecine Nucleaire, avenue Pasteur, 98000 Monaco and Centre Antoine Lacassagne, Cyclotron Biomedical, 227 avenue de la Lanterne, 06200 Nice},
abstractNote = {The aim of this work was to study the dosimetric potential of the Monte Carlo code MCNPX applied to the protontherapy field. For series of clinical configurations a comparison between simulated and experimental data was carried out, using the proton beam line of the MEDICYC isochronous cyclotron installed in the Centre Antoine Lacassagne in Nice. The dosimetric quantities tested were depth-dose distributions, output factors, and monitor units. For each parameter, the simulation reproduced accurately the experiment, which attests the quality of the choices made both in the geometrical description and in the physics parameters for beam definition. These encouraging results enable us today to consider a simplification of quality control measurements in the future. Monitor Units calculation is planned to be carried out with preestablished Monte Carlo simulation data. The measurement, which was until now our main patient dose calibration system, will be progressively replaced by computation based on the MCNPX code. This determination of Monitor Units will be controlled by an independent semi-empirical calculation.},
doi = {10.1118/1.2431473},
journal = {Medical Physics},
number = 2,
volume = 34,
place = {United States},
year = {Thu Feb 15 00:00:00 EST 2007},
month = {Thu Feb 15 00:00:00 EST 2007}
}
  • Purpose: In broad-beam carbon-ion radiotherapy performed using the heavy-ion medical accelerator in Chiba, the number of monitor units is determined by measuring the physical dose at the center of the spread-out Bragg peak (SOBP) for the treatment beam. The total measurement time increases as the number of treatment beams increases, which hinders the treatment of an increased number of patients. Hence, Kusano et al.[Jpn. J. Med. Phys. 23(Suppl. 2), 65-68 (2003)] proposed a method to calculate the physical dose at the center of the SOBP for a treatment beam. Based on a recent study, the authors here propose a moremore » accurate calculation method.Methods: The authors measured the physical dose at the center of the SOBP while varying the circular field size and range-shifter thickness. The authors obtained the physical dose at the center of the SOBP for an irregularly shaped beam using Clarkson integration based on these measurements.Results: The difference between the calculated and measured physical doses at the center of the SOBP varied with a change in the central angle of the sector segment. The differences between the calculated and measured physical doses at the center of the SOBP were within {+-}1% for all irregularly shaped beams that were used to validate the calculation method.Conclusions: The accuracy of the proposed method depends on both the number of angular intervals used for Clarkson integration and the fineness of the basic data used for calculations: sampling numbers for the field size and thickness of the range shifter. If those parameters are properly chosen, the authors can obtain a calculated monitor unit number with high accuracy sufficient for clinical applications.« less
  • Purpose: Aim of this study was to analyze the modulating, broadening effect on the Bragg Peak due to heterogeneous geometries like multi-wire chambers in the beam path of a particle therapy beam line. The effect was described by a mathematical model which was implemented in the Monte-Carlo code FLUKA via user-routines, in order to reduce the computation time for the simulations. Methods: The depth dose curve of 80 MeV/u C12-ions in a water phantom was calculated using the Monte-Carlo code FLUKA (reference curve). The modulating effect on this dose distribution behind eleven mesh-like foils (periodicity ∼80 microns) occurring in amore » typical set of multi-wire and dose chambers was mathematically described by optimizing a normal distribution so that the reverence curve convoluted with this distribution equals the modulated dose curve. This distribution describes a displacement in water and was transferred in a probability distribution of the thickness of the eleven foils using the water equivalent thickness of the foil’s material. From this distribution the distribution of the thickness of one foil was determined inversely. In FLUKA the heterogeneous foils were replaced by homogeneous foils and a user-routine was programmed that varies the thickness of the homogeneous foils for each simulated particle using this distribution. Results: Using the mathematical model and user-routine in FLUKA the broadening effect could be reproduced exactly when replacing the heterogeneous foils by homogeneous ones. The computation time was reduced by 90 percent. Conclusion: In this study the broadening effect on the Bragg Peak due to heterogeneous structures was analyzed, described by a mathematical model and implemented in FLUKA via user-routines. Applying these routines the computing time was reduced by 90 percent. The developed tool can be used for any heterogeneous structure in the dimensions of microns to millimeters, in principle even for organic materials like lung tissue.« less
  • Purpose: The application of optically stimulated luminescence dosimeters (OSLDs) may be extended to clinical investigations verifying irradiated doses in small animal models. In proton beams, the accurate positioning of the Bragg peak is essential for tumor targeting. The purpose of this study was to estimate the displacement of a pristine Bragg peak when an Al2O3:C nanodot (Landauer, Inc.) is placed on the surface of a water phantom and to evaluate corresponding changes in dose. Methods: Clinical proton pencil beam simulations were carried out with using TOPAS, a Monte Carlo platform layered on top of GEANT4. Point-shaped beams with no energymore » spread were modeled for energies 100MV, 150MV, 200MV, and 250MV. Dose scoring for 100,000 particle histories was conducted within a water phantom (20cm × 20cm irradiated area, 40cm depth) with its surface placed 214.5cm away from the source. The modeled nanodot had a 4mm radius and 0.2mm thickness. Results: A comparative analysis of Monte Carlo depth dose profiles modeled for these proton pencil beams did not demonstrate an energy dependent in the Bragg peak shift. The shifts in Bragg Peak depth for water phantoms modeled with a nanodot on the phantom surface ranged between 2.7 to 3.2 mm. In all cases, the Bragg Peaks were shifted closer to the irradiation source. The peak dose in phantoms with an OSLD remained unchanged with percent dose differences less than 0.55% when compared to phantom doses without the nanodot. Conclusion: Monte Carlo calculations show that the presence of OSLD nanodots in proton beam therapy will not change the position of a pristine Bragg Peak by more than 3 mm. Although the 3.0 mm shift will not have a detrimental effect in patients receiving proton therapy, this effect may not be negligible in dose verification measurements for mouse models at lower proton beam energies.« less
  • Proton treatments with spread-out Bragg peak (SOBP) fields often use a rotating modulation wheel of varying thickness to modulate the pristine Bragg peak in depth and intensity. The technique of modulating also the beam current independently over the wheel rotation provides an additional control over the intensities of the pulled-back Bragg peaks. As a result, a single wheel can be used over a large range of energies and SOBP parameters and field-specific wheels are no longer necessary. An essential task in commissioning a particular treatment depth is the determination of this current modulation function. We have developed a method formore » the optimization of the current modulation function. The basic idea is to treat the entire beam nozzle, housing the various beam scattering and modulating components, as a whole and to characterize its effect as a transformation from a modulating beam current to a depth-dose distribution. While this transformation is difficult to calculate theoretically due to the complex scattering paths in the nozzle and the phantom, it can, however, be determined by time-resolved dose measurements. Using this transformation, we can calculate SOBP depth-dose distributions for any current modulation function and optimize it by a simple numerical optimization. We have applied the new method to a number of proton beams with satisfactory results.« less
  • A spread-out Bragg peak (SOBP) is used in proton beam therapy to create a longitudinal conformality of the required dose to the target. In order to create this effect in a passive beam scattering system, a variety of components must operate in conjunction to produce the desired beam parameters. We will describe how the SOBP is generated and will explore the tolerances of the various components and their subsequent effect on the dose distribution. A specific aspect of this investigation includes a case study involving the use of a beam current modulated system. In such a system, the intensity ofmore » the beam current can be varied in synchronization with the revolution of the range-modulator wheel. As a result, the weights of the pulled-back Bragg peaks can be individually controlled to produce uniform dose plateaus for a large range of treatment depths using only a small number of modulator wheels.« less