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

Title: Application of adjoint Monte Carlo to accelerate simulations of mono-directional beams in treatment planning for Boron Neutron Capture Therapy

Abstract

This paper deals with the application of the adjoint transport theory in order to optimize Monte Carlo based radiotherapy treatment planning. The technique is applied to Boron Neutron Capture Therapy where most often mixed beams of neutrons and gammas are involved. In normal forward Monte Carlo simulations the particles start at a source and lose energy as they travel towards the region of interest, i.e., the designated point of detection. Conversely, with adjoint Monte Carlo simulations, the so-called adjoint particles start at the region of interest and gain energy as they travel towards the source where they are detected. In this respect, the particles travel backwards and the real source and real detector become the adjoint detector and adjoint source, respectively. At the adjoint detector, an adjoint function is obtained with which numerically the same result, e.g., dose or flux in the tumor, can be derived as with forward Monte Carlo. In many cases, the adjoint method is more efficient and by that is much quicker when, for example, the response in the tumor or organ at risk for many locations and orientations of the treatment beam around the patient is required. However, a problem occurs when the treatment beammore » is mono-directional as the probability of detecting adjoint Monte Carlo particles traversing the beam exit (detector plane in adjoint mode) in the negative direction of the incident beam is zero. This problem is addressed here and solved first with the use of next event estimators and second with the application of a Legendre expansion technique of the angular adjoint function. In the first approach, adjoint particles are tracked deterministically through a tube to a (adjoint) point detector far away from the geometric model. The adjoint particles will traverse the disk shaped entrance of this tube (the beam exit in the actual geometry) perpendicularly. This method is slow whenever many events are involved that are not contributing to the point detector, e.g., neutrons in a scattering medium. In the second approach, adjoint particles that traverse an adjoint shaped detector plane are used to estimate the Legendre coefficients for expansion of the angular adjoint function. This provides an estimate of the adjoint function for the direction normal to the detector plane. In a realistic head model, as described in this paper, which is surrounded by 1020 mono-directional neutron/gamma beams and from which the best ones are to be selected, the example calculates the neutron and gamma fluxes in ten tumors and ten organs at risk. For small diameter beams (5 cm), and with comparable relative errors, forward Monte Carlo is seen to be 1.5 times faster than the adjoint Monte Carlo techniques. For larger diameter neutron beams (10 and 15 cm), the Legendre technique is found to be 6 and 20 times faster, respectively. In the case of gammas alone, for the 10 and 15 cm diam beams, both adjoint Monte Carlo Legendre and point detector techniques are respectively 2 and 3 times faster than forward Monte Carlo.« less

Authors:
; ; ; ; ;  [1];  [2];  [2];  [2];  [2];  [2]
  1. Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft (Netherlands) and Joint Research Centre of the European Commission, P.O. Box 2, 1755ZG Petten (Netherlands)
  2. (Netherlands)
Publication Date:
OSTI Identifier:
20951150
Resource Type:
Journal Article
Resource Relation:
Journal Name: Medical Physics; Journal Volume: 34; Journal Issue: 4; Other Information: DOI: 10.1118/1.2712573; (c) 2007 American Association of Physicists in Medicine; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
62 RADIOLOGY AND NUCLEAR MEDICINE; COMPUTERIZED SIMULATION; ERRORS; MONTE CARLO METHOD; NEUTRON BEAMS; NEUTRON CAPTURE THERAPY; OPTIMIZATION; PHANTOMS; RADIATION DOSES

Citation Formats

Nievaart, V. A., Legrady, D., Moss, R. L., Kloosterman, J. L., Hagen, T. H. J. J. van der, Dam, H. van, Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft, Joint Research Centre of the European Commission, P.O. Box 2, 1755ZG Petten, Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft, Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft, The Netherlands and Reactor Institute Delft, Delft University Technology, Mekelweg 15, 2629JB Delft, and Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft. Application of adjoint Monte Carlo to accelerate simulations of mono-directional beams in treatment planning for Boron Neutron Capture Therapy. United States: N. p., 2007. Web. doi:10.1118/1.2712573.
Nievaart, V. A., Legrady, D., Moss, R. L., Kloosterman, J. L., Hagen, T. H. J. J. van der, Dam, H. van, Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft, Joint Research Centre of the European Commission, P.O. Box 2, 1755ZG Petten, Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft, Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft, The Netherlands and Reactor Institute Delft, Delft University Technology, Mekelweg 15, 2629JB Delft, & Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft. Application of adjoint Monte Carlo to accelerate simulations of mono-directional beams in treatment planning for Boron Neutron Capture Therapy. United States. doi:10.1118/1.2712573.
Nievaart, V. A., Legrady, D., Moss, R. L., Kloosterman, J. L., Hagen, T. H. J. J. van der, Dam, H. van, Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft, Joint Research Centre of the European Commission, P.O. Box 2, 1755ZG Petten, Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft, Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft, The Netherlands and Reactor Institute Delft, Delft University Technology, Mekelweg 15, 2629JB Delft, and Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft. Sun . "Application of adjoint Monte Carlo to accelerate simulations of mono-directional beams in treatment planning for Boron Neutron Capture Therapy". United States. doi:10.1118/1.2712573.
@article{osti_20951150,
title = {Application of adjoint Monte Carlo to accelerate simulations of mono-directional beams in treatment planning for Boron Neutron Capture Therapy},
author = {Nievaart, V. A. and Legrady, D. and Moss, R. L. and Kloosterman, J. L. and Hagen, T. H. J. J. van der and Dam, H. van and Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft and Joint Research Centre of the European Commission, P.O. Box 2, 1755ZG Petten and Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft and Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft, The Netherlands and Reactor Institute Delft, Delft University Technology, Mekelweg 15, 2629JB Delft and Department of Applied Sciences, Delft University of Technology, Lorentzweg 1, 2628CJ Delft},
abstractNote = {This paper deals with the application of the adjoint transport theory in order to optimize Monte Carlo based radiotherapy treatment planning. The technique is applied to Boron Neutron Capture Therapy where most often mixed beams of neutrons and gammas are involved. In normal forward Monte Carlo simulations the particles start at a source and lose energy as they travel towards the region of interest, i.e., the designated point of detection. Conversely, with adjoint Monte Carlo simulations, the so-called adjoint particles start at the region of interest and gain energy as they travel towards the source where they are detected. In this respect, the particles travel backwards and the real source and real detector become the adjoint detector and adjoint source, respectively. At the adjoint detector, an adjoint function is obtained with which numerically the same result, e.g., dose or flux in the tumor, can be derived as with forward Monte Carlo. In many cases, the adjoint method is more efficient and by that is much quicker when, for example, the response in the tumor or organ at risk for many locations and orientations of the treatment beam around the patient is required. However, a problem occurs when the treatment beam is mono-directional as the probability of detecting adjoint Monte Carlo particles traversing the beam exit (detector plane in adjoint mode) in the negative direction of the incident beam is zero. This problem is addressed here and solved first with the use of next event estimators and second with the application of a Legendre expansion technique of the angular adjoint function. In the first approach, adjoint particles are tracked deterministically through a tube to a (adjoint) point detector far away from the geometric model. The adjoint particles will traverse the disk shaped entrance of this tube (the beam exit in the actual geometry) perpendicularly. This method is slow whenever many events are involved that are not contributing to the point detector, e.g., neutrons in a scattering medium. In the second approach, adjoint particles that traverse an adjoint shaped detector plane are used to estimate the Legendre coefficients for expansion of the angular adjoint function. This provides an estimate of the adjoint function for the direction normal to the detector plane. In a realistic head model, as described in this paper, which is surrounded by 1020 mono-directional neutron/gamma beams and from which the best ones are to be selected, the example calculates the neutron and gamma fluxes in ten tumors and ten organs at risk. For small diameter beams (5 cm), and with comparable relative errors, forward Monte Carlo is seen to be 1.5 times faster than the adjoint Monte Carlo techniques. For larger diameter neutron beams (10 and 15 cm), the Legendre technique is found to be 6 and 20 times faster, respectively. In the case of gammas alone, for the 10 and 15 cm diam beams, both adjoint Monte Carlo Legendre and point detector techniques are respectively 2 and 3 times faster than forward Monte Carlo.},
doi = {10.1118/1.2712573},
journal = {Medical Physics},
number = 4,
volume = 34,
place = {United States},
year = {Sun Apr 15 00:00:00 EDT 2007},
month = {Sun Apr 15 00:00:00 EDT 2007}
}