Skip to main content

Lund University Publications

LUND UNIVERSITY LIBRARIES

Monte Carlo calculations of thermal neutron capture in gadolinium : A comparison of GEANT4 and MCNP with measurements

Enger, Shirin A. ; Af Rosenschöld, Per Munck LU orcid ; Rezaei, Arash and Lundqvist, Hans (2006) In Medical Physics 33(2). p.337-341
Abstract

GEANT4 is a Monte Carlo code originally implemented for high-energy physics applications and is well known for particle transport at high energies. The capacity of GEANT4 to simulate neutron transport in the thermal energy region is not equally well known. The aim of this article is to compare MCNP, a code commonly used in low energy neutron transport calculations and GEANT4 with experimental results and select the suitable code for gadolinium neutron capture applications. To account for the thermal neutron scattering from chemically bound atoms [S(α,β)] in biological materials a comparison of thermal neutron fluence in tissue-like poly(methylmethacrylate) phantom is made with MCNP4B, GEANT4 6.0 patch1, and measurements from the neutron... (More)

GEANT4 is a Monte Carlo code originally implemented for high-energy physics applications and is well known for particle transport at high energies. The capacity of GEANT4 to simulate neutron transport in the thermal energy region is not equally well known. The aim of this article is to compare MCNP, a code commonly used in low energy neutron transport calculations and GEANT4 with experimental results and select the suitable code for gadolinium neutron capture applications. To account for the thermal neutron scattering from chemically bound atoms [S(α,β)] in biological materials a comparison of thermal neutron fluence in tissue-like poly(methylmethacrylate) phantom is made with MCNP4B, GEANT4 6.0 patch1, and measurements from the neutron capture therapy (NCT) facility at the Studsvik, Sweden. The fluence measurements agreed with MCNP calculated results considering S(α,β). The location of the thermal neutron peak calculated with MCNP without S(α,β) and GEANT4 is shifted by about 0.5 cm towards a shallower depth and is 25%-30% lower in amplitude. Dose distribution from the gadolinium neutron capture reaction is then simulated by MCNP and compared with measured data. The simulations made by MCNP agree well with experimental results. As long as thermal neutron scattering from chemically bound atoms are not included in GEANT4 it is not suitable for NCT applications.

(Less)
Please use this url to cite or link to this publication:
author
; ; and
publishing date
type
Contribution to journal
publication status
published
subject
keywords
MCNP, GEANT4, Monte Carlo calculations, neutron capture therapy
in
Medical Physics
volume
33
issue
2
pages
337 - 341
publisher
American Association of Physicists in Medicine
external identifiers
  • pmid:16532938
  • wos:000235682600010
  • scopus:31844438710
  • pmid:16532938
ISSN
0094-2405
DOI
10.1118/1.2150787
language
English
LU publication?
no
id
0243a915-c369-4321-ae1e-272e04a1b36f (old id 416692)
date added to LUP
2016-04-01 15:19:03
date last changed
2024-01-10 13:38:38
@article{0243a915-c369-4321-ae1e-272e04a1b36f,
  abstract     = {{<p>GEANT4 is a Monte Carlo code originally implemented for high-energy physics applications and is well known for particle transport at high energies. The capacity of GEANT4 to simulate neutron transport in the thermal energy region is not equally well known. The aim of this article is to compare MCNP, a code commonly used in low energy neutron transport calculations and GEANT4 with experimental results and select the suitable code for gadolinium neutron capture applications. To account for the thermal neutron scattering from chemically bound atoms [S(α,β)] in biological materials a comparison of thermal neutron fluence in tissue-like poly(methylmethacrylate) phantom is made with MCNP4B, GEANT4 6.0 patch1, and measurements from the neutron capture therapy (NCT) facility at the Studsvik, Sweden. The fluence measurements agreed with MCNP calculated results considering S(α,β). The location of the thermal neutron peak calculated with MCNP without S(α,β) and GEANT4 is shifted by about 0.5 cm towards a shallower depth and is 25%-30% lower in amplitude. Dose distribution from the gadolinium neutron capture reaction is then simulated by MCNP and compared with measured data. The simulations made by MCNP agree well with experimental results. As long as thermal neutron scattering from chemically bound atoms are not included in GEANT4 it is not suitable for NCT applications.</p>}},
  author       = {{Enger, Shirin A. and Af Rosenschöld, Per Munck and Rezaei, Arash and Lundqvist, Hans}},
  issn         = {{0094-2405}},
  keywords     = {{MCNP; GEANT4; Monte Carlo calculations; neutron capture therapy}},
  language     = {{eng}},
  number       = {{2}},
  pages        = {{337--341}},
  publisher    = {{American Association of Physicists in Medicine}},
  series       = {{Medical Physics}},
  title        = {{Monte Carlo calculations of thermal neutron capture in gadolinium : A comparison of GEANT4 and MCNP with measurements}},
  url          = {{http://dx.doi.org/10.1118/1.2150787}},
  doi          = {{10.1118/1.2150787}},
  volume       = {{33}},
  year         = {{2006}},
}