Monte Carlo modeling of an electronic brachytherapy source using MCNP5 and EGSnrc Stephen D. Davis and Larry A. DeWerd University of Wisconsin, Madison, WI 13 th UK Monte Carlo User Group meeting March 29, 2007
Xoft Axxent TM electronic brachytherapy source 2 mm 2 mm Anode 10 mm Cathode Photos courtesy of Xoft, Inc. 2/26
Xoft Axxent TM electronic brachytherapy source X-Ray Probe Tip Detail Photos courtesy of Xoft, Inc. 3/26
Motivation Brachytherapy uses air-kerma strength, S K, as the measurement of source strength Air kerma rates can be measured using free-air ionization chambers, but S K is defined in vacuo Correction is required for attenuation in air from the source to the point of measurement (100 cm in our case) Air cross sections change rapidly at photon energies < 20 kev Photon spectrum needs to be accurately known for a low uncertainty on the air attenuation correction 4/26
Photon spectra Measurements with high purity germanium spectrometer Requires correction for the energy response of the detector Monte Carlo simulations Assumes Monte Carlo code has the low energy physics included to model x-ray spectra accurately 5/26
Photon spectra at 178 cm in air 120 Counts (arb. units) 100 80 60 40 NIST HPGe MCNP5 20 0 0 10 20 30 40 50 Photon energy (kev) 6/26
Photon spectra at 178 cm in air 30 25 NIST HPGe MCNP5 Counts (arb. units) 20 15 10 W-L β W-L γ 5 W-L α 0 5 6 7 8 9 10 11 12 13 14 15 Photon energy (kev) 7/26
Photon spectra at 0 cm in air 180 160 Counts (arb. units) 140 120 100 80 60 NIST HPGe MCNP5 40 20 0 0 10 20 30 40 50 Photon energy (kev) 8/26
EGSnrc January 2005 release included model by Iwan Kawrakow for electron impact ionization for K- and L- shells with binding energies above 1 kev October 2005 release included C++ class library with general purpose geometry modeling and user code cavity.cpp February 2007 release allowed for scoring of photon spectra in a circular plane using cavity.cpp 9/26
EGSnrc vs. MCNP5 Two codes have differences in: Electron impact ionization Binding effects and Doppler broadening for Compton interactions Atomic relaxation Electron transport Variance reduction techniques Scoring options 10/26
Doppler broadening Broomstick problem from MCNP5 LANL document 100 kev photons incident on 10 cm graphite rod with 10-6 cm radius 11/26
No Doppler broadening 10-1 10-2 Fluence per unit incident fluence 10-3 10-4 10-5 10-6 10-7 10-8 10-9 10-10 10-11 10-12 10-13 EGSnrc MCNP5 10-14 0 20 40 60 80 100 Photon energy (kev) 12/26
With Doppler broadening 10-1 10-2 Fluence per unit incident fluence 10-3 10-4 10-5 10-6 10-7 10-8 10-9 10-10 10-11 10-12 10-13 EGSnrc MCNP5 10-14 0 20 40 60 80 100 Photon energy (kev) 13/26
Atomic relaxation 30 kev photons incident on 10-5 cm thick yttrium slab 14/26
Atomic relaxation (cont.) 10 2 Fluence per unit incident fluence 10 1 10 0 10-1 10-2 10-3 10-4 EGSnrc MCNP5 10-5 0 5 10 15 20 25 30 Photon energy (kev) 15/26
Atomic relaxation (cont.) 16 14 Fluence per unit incident fluence 12 10 8 6 4 Y-K α EGSnrc MCNP5 (Mode p,e) Y-K β 2 0 13 14 15 16 17 18 19 Photon energy (kev) 16/26
Atomic relaxation (cont.) 16 14 Fluence per unit incident fluence 12 10 8 6 4 Y-K α EGSnrc MCNP5 (Mode p) Y-K β 2 0 13 14 15 16 17 18 19 Photon energy (kev) 17/26
cavity.cpp model of Axxent TM source 18/26
Dose rates at 1 cm in water MC code P (r 0,θ 0 ) [cgy h -1 µa -1 ] MCNP5 345 ± 0.1% EGSnrc 346 ± 2.0% EGSnrc/MCNP5 1.003 ± 2.0% 19/26
Air kerma rates at 100 cm in air MC code P (r 0,θ 0 ) [cgy h -1 µa -1 ] MCNP5 4.81 10-2 ± 0.1% EGSnrc 4.72 10-2 ± 3.1% EGSnrc/MCNP5 0.981 ± 3.1% 20/26
Calculated photon spectra at 1 cm in vacuum 9x10-4 8x10-4 Fluence per unit incident fluence 7x10-4 6x10-4 5x10-4 4x10-4 3x10-4 2x10-4 1x10-4 EGSnrc MCNP5 3.2% higher air kerma for EGSnrc 0 0 5 10 15 20 25 30 35 40 45 50 55 Photon energy (kev) 21/26
Calculated photon spectra at 1 cm in vacuum 1.2x10-4 Fluence per unit incident fluence 1.0x10-4 8.0x10-5 6.0x10-5 4.0x10-5 2.0x10-5 EGSnrc MCNP5 <L>-<M> L III -<M> L II -<M> L II -<N> L I -<M> + L III -<N> EII contributes 2.6% of air kerma L I -<N> 0.0 5 10 15 Photon energy (kev) 22/26
Calculated photon spectra at 1 cm in vacuum 9x10-4 Fluence per unit incident fluence 8x10-4 7x10-4 6x10-4 5x10-4 4x10-4 3x10-4 2x10-4 1x10-4 Y-K α EGSnrc MCNP5 Differences in Y fluorescent peaks led to 2.3% difference in air kerma Y-K β 0 14 15 16 17 18 Photon energy (kev) 23/26
Calculated photon spectra at 1 cm in vacuum 2x10-4 Fluence per unit incident fluence 1x10-4 5x10-5 Removing EII and Y fluorescence differences, EGSnrc produces 2.4% lower air kerma than MCNP5 EGSnrc MCNP5 0 20 25 30 35 40 45 50 55 Photon energy (kev) 24/26
Conclusions Both codes produce similar water dose rates, air kerma rates, and photon spectra, and the small differences can be explained by different treatment of the low energy physics Further work will be necessary to compare these results to measurements, and to establish the required air attenuation corrections with low uncertainties 25/26
Acknowledgements Xoft, Inc. Iwan Kawrakow NIST UWMRRC staff and students UW ADCL customers Grid Laboratory of Wisconsin (GLOW) 26/26