Evaluation of the new electron-transport algorithm in MCNP6.1 for the simulation of dose point kernel in water

被引:9
作者
Antoni, Rodolphe [1 ]
Bourgois, Laurent
机构
[1] CEA, DEN, CAD, F-13108 Cadarache, St Paul Lez Dur, France
关键词
MCNP6; Electrons transport algorithm; Condensed history method; Single event algorithm; Dose Point Kernel (DPK); MONTE-CARLO CODES; DISTRIBUTIONS; ENERGIES; EGSNRC;
D O I
10.1016/j.nimb.2017.09.026
中图分类号
TH7 [仪器、仪表];
学科分类号
0804 ; 080401 ; 081102 ;
摘要
In this work, the calculation of specific dose distribution in water is evaluated in MCNP6.1 with the regular condensed history algorithm the "detailed electron energy-loss straggling logic" and the new electrons transport algorithm proposed the "single event algorithm". Dose Point Kernel (DPK) is calculated with monoenergetic electrons of 50, 100, 500, 1000 and 3000 keV for different scoring cells dimensions. A comparison between MCNP6 results and well-validated codes for electron-dosimetry, i.e., EGSnrc or Penelope, is performed. When the detailed electron energy-loss straggling logic is used with default setting (down to the cut-off energy 1 keV), we infer that the depth of the dose peak increases with decreasing thickness of the scoring cell, largely due to combined step-size and boundary crossing artifacts. This finding is less prominent for 500 key, 1 MeV and 3 MeV dose profile. With an appropriate number of sub-steps (ESTEP value in MCNP6), the dose-peak shift is almost complete absent to 50 keV and 100 keV electrons. However, the dose-peak is more prominent compared to EGSnrc and the absorbed dose tends to be underestimated at greater depths, meaning that boundaries crossing artifact are still occurring while step-size artifacts are greatly reduced. When the single-event mode is used for the whole transport, we observe the good agreement of reference and calculated profile for 50 and 100 keV electrons. Remaining artifacts are fully vanished, showing a possible transport treatment for energies less than a hundred of keV and accordance with reference for whatever scoring cell dimension, even if the single event method initially intended to support electron transport at energies below 1 keV. Conversely, results for 500 keNT, 1 MeV and 3 MeV undergo a dramatic discrepancy with reference curves. These poor results and so the current unreliability of the method is for a part due to inappropriate elastic cross section treatment from the ENDF/B-VI.8 library in those energy ranges. Accordingly, special care has to be taken in setting choice for calculating electron dose distribution with MCNP6, in particular with regards to dosimetry or nuclear medicine applications.
引用
收藏
页码:102 / 108
页数:7
相关论文
共 26 条
  • [21] Evaluation of computational models and cross sections used by MCNP6 for simulation of electron backscattering
    Poskus, Andrius
    [J]. NUCLEAR INSTRUMENTS & METHODS IN PHYSICS RESEARCH SECTION B-BEAM INTERACTIONS WITH MATERIALS AND ATOMS, 2016, 368 : 15 - 27
  • [22] Parameter dependence of the MCNP electron transport in determining dose distributions
    Reynaert, N
    Palmans, H
    Thierens, H
    Jeraj, R
    [J]. MEDICAL PHYSICS, 2002, 29 (10) : 2446 - 2454
  • [23] Salvat F., 2015, PENELOPE2014 CODE SY
  • [24] A comparison of MCNP4C electron transport with ITS 3.0 and experiment at incident energies between 100 keV and 20 MeV: influence of voxel size, substeps and energy indexing algorithm
    Schaart, DR
    Jansen, JTM
    Zoetelief, J
    de Leege, PFA
    [J]. PHYSICS IN MEDICINE AND BIOLOGY, 2002, 47 (09) : 1459 - 1484
  • [25] DOSE POINT KERNEL SIMULATION FOR MONOENERGETIC ELECTRONS AND RADIONUCLIDES USING MONTE CARLO TECHNIQUES
    Wu, J.
    Liu, Y. L.
    Chang, S. J.
    Chao, M. M.
    Tsai, S. Y.
    Huang, D. E.
    [J]. RADIATION PROTECTION DOSIMETRY, 2012, 152 (1-3) : 119 - 124
  • [26] X-5 Monte Carlo Team, MCNP 2003 GEN PURP M