Experimentally validated pencil beam scanning source model in TOPAS

Liyong Lin, Minglei Kang, Timothy D. Solberg, Christopher G. Ainsley, James E. McDonough

Research output: Contribution to journalArticle

18 Citations (Scopus)

Abstract

The presence of a low-dose envelope, or 'halo', in the fluence profile of a proton spot can increase the output of a pencil beam scanning field by over 10%. This study evaluated whether the Monte Carlo simulation code, TOPAS 1.0-beta 8, based on Geant4.9.6 with its default physics list, can predict the spot halo at depth in phantom by incorporating a halo model within the proton source distribution. Proton sources were modelled using three 2D Gaussian functions, and optimized until simulated spot profiles matched measurements at the phantom surface out to a radius of 100 mm. Simulations were subsequently compared with profiles measured using EBT3 film in Solidwater® phantoms at various depths for 100, 115, 150, 180, 210 and 225 MeV proton beams. Simulations predict measured profiles within a 1 mm distance to agreement for 2D profiles extending to the 0.1% isodose, and within 1 mm/1% Gamma criteria over the integrated curve of spot profile as a function of radius. For isodose lines beyond 0.1% of the central spot dose, the simulated primary spot sigma is smaller than the measurement by up to 15%, and can differ by over 1 mm. The choice of particle interaction algorithm and phantom material were found to cause ∼1 mm range uncertainty, a maximal 5% (0.3 mm) difference in spot sigma, and maximal 1 mm and ∼2 mm distance to agreement in isodoses above and below the 0.1% level, respectively. Based on these observations, therefore, the selection of physics model and the application of Solidwater® as water replacement material in simulation and measurement should be used with caution.

Original languageEnglish (US)
Pages (from-to)6859-6873
Number of pages15
JournalPhysics in Medicine and Biology
Volume59
Issue number22
DOIs
StatePublished - Nov 21 2014

Fingerprint

Protons
Physics
Uncertainty
Water

Keywords

  • Monte Carlo simulation
  • pencil beam scanning
  • proton therapy
  • Solidwater®
  • treatment planning

ASJC Scopus subject areas

  • Radiology Nuclear Medicine and imaging
  • Radiological and Ultrasound Technology

Cite this

Lin, L., Kang, M., Solberg, T. D., Ainsley, C. G., & McDonough, J. E. (2014). Experimentally validated pencil beam scanning source model in TOPAS. Physics in Medicine and Biology, 59(22), 6859-6873. https://doi.org/10.1088/0031-9155/59/22/6859

Experimentally validated pencil beam scanning source model in TOPAS. / Lin, Liyong; Kang, Minglei; Solberg, Timothy D.; Ainsley, Christopher G.; McDonough, James E.

In: Physics in Medicine and Biology, Vol. 59, No. 22, 21.11.2014, p. 6859-6873.

Research output: Contribution to journalArticle

Lin, L, Kang, M, Solberg, TD, Ainsley, CG & McDonough, JE 2014, 'Experimentally validated pencil beam scanning source model in TOPAS', Physics in Medicine and Biology, vol. 59, no. 22, pp. 6859-6873. https://doi.org/10.1088/0031-9155/59/22/6859
Lin, Liyong ; Kang, Minglei ; Solberg, Timothy D. ; Ainsley, Christopher G. ; McDonough, James E. / Experimentally validated pencil beam scanning source model in TOPAS. In: Physics in Medicine and Biology. 2014 ; Vol. 59, No. 22. pp. 6859-6873.
@article{1b17d492f7ad4452b987c6e0ad2cb19c,
title = "Experimentally validated pencil beam scanning source model in TOPAS",
abstract = "The presence of a low-dose envelope, or 'halo', in the fluence profile of a proton spot can increase the output of a pencil beam scanning field by over 10{\%}. This study evaluated whether the Monte Carlo simulation code, TOPAS 1.0-beta 8, based on Geant4.9.6 with its default physics list, can predict the spot halo at depth in phantom by incorporating a halo model within the proton source distribution. Proton sources were modelled using three 2D Gaussian functions, and optimized until simulated spot profiles matched measurements at the phantom surface out to a radius of 100 mm. Simulations were subsequently compared with profiles measured using EBT3 film in Solidwater{\circledR} phantoms at various depths for 100, 115, 150, 180, 210 and 225 MeV proton beams. Simulations predict measured profiles within a 1 mm distance to agreement for 2D profiles extending to the 0.1{\%} isodose, and within 1 mm/1{\%} Gamma criteria over the integrated curve of spot profile as a function of radius. For isodose lines beyond 0.1{\%} of the central spot dose, the simulated primary spot sigma is smaller than the measurement by up to 15{\%}, and can differ by over 1 mm. The choice of particle interaction algorithm and phantom material were found to cause ∼1 mm range uncertainty, a maximal 5{\%} (0.3 mm) difference in spot sigma, and maximal 1 mm and ∼2 mm distance to agreement in isodoses above and below the 0.1{\%} level, respectively. Based on these observations, therefore, the selection of physics model and the application of Solidwater{\circledR} as water replacement material in simulation and measurement should be used with caution.",
keywords = "Monte Carlo simulation, pencil beam scanning, proton therapy, Solidwater{\circledR}, treatment planning",
author = "Liyong Lin and Minglei Kang and Solberg, {Timothy D.} and Ainsley, {Christopher G.} and McDonough, {James E.}",
year = "2014",
month = "11",
day = "21",
doi = "10.1088/0031-9155/59/22/6859",
language = "English (US)",
volume = "59",
pages = "6859--6873",
journal = "Physics in Medicine and Biology",
issn = "0031-9155",
publisher = "IOP Publishing Ltd.",
number = "22",

}

TY - JOUR

T1 - Experimentally validated pencil beam scanning source model in TOPAS

AU - Lin, Liyong

AU - Kang, Minglei

AU - Solberg, Timothy D.

AU - Ainsley, Christopher G.

AU - McDonough, James E.

PY - 2014/11/21

Y1 - 2014/11/21

N2 - The presence of a low-dose envelope, or 'halo', in the fluence profile of a proton spot can increase the output of a pencil beam scanning field by over 10%. This study evaluated whether the Monte Carlo simulation code, TOPAS 1.0-beta 8, based on Geant4.9.6 with its default physics list, can predict the spot halo at depth in phantom by incorporating a halo model within the proton source distribution. Proton sources were modelled using three 2D Gaussian functions, and optimized until simulated spot profiles matched measurements at the phantom surface out to a radius of 100 mm. Simulations were subsequently compared with profiles measured using EBT3 film in Solidwater® phantoms at various depths for 100, 115, 150, 180, 210 and 225 MeV proton beams. Simulations predict measured profiles within a 1 mm distance to agreement for 2D profiles extending to the 0.1% isodose, and within 1 mm/1% Gamma criteria over the integrated curve of spot profile as a function of radius. For isodose lines beyond 0.1% of the central spot dose, the simulated primary spot sigma is smaller than the measurement by up to 15%, and can differ by over 1 mm. The choice of particle interaction algorithm and phantom material were found to cause ∼1 mm range uncertainty, a maximal 5% (0.3 mm) difference in spot sigma, and maximal 1 mm and ∼2 mm distance to agreement in isodoses above and below the 0.1% level, respectively. Based on these observations, therefore, the selection of physics model and the application of Solidwater® as water replacement material in simulation and measurement should be used with caution.

AB - The presence of a low-dose envelope, or 'halo', in the fluence profile of a proton spot can increase the output of a pencil beam scanning field by over 10%. This study evaluated whether the Monte Carlo simulation code, TOPAS 1.0-beta 8, based on Geant4.9.6 with its default physics list, can predict the spot halo at depth in phantom by incorporating a halo model within the proton source distribution. Proton sources were modelled using three 2D Gaussian functions, and optimized until simulated spot profiles matched measurements at the phantom surface out to a radius of 100 mm. Simulations were subsequently compared with profiles measured using EBT3 film in Solidwater® phantoms at various depths for 100, 115, 150, 180, 210 and 225 MeV proton beams. Simulations predict measured profiles within a 1 mm distance to agreement for 2D profiles extending to the 0.1% isodose, and within 1 mm/1% Gamma criteria over the integrated curve of spot profile as a function of radius. For isodose lines beyond 0.1% of the central spot dose, the simulated primary spot sigma is smaller than the measurement by up to 15%, and can differ by over 1 mm. The choice of particle interaction algorithm and phantom material were found to cause ∼1 mm range uncertainty, a maximal 5% (0.3 mm) difference in spot sigma, and maximal 1 mm and ∼2 mm distance to agreement in isodoses above and below the 0.1% level, respectively. Based on these observations, therefore, the selection of physics model and the application of Solidwater® as water replacement material in simulation and measurement should be used with caution.

KW - Monte Carlo simulation

KW - pencil beam scanning

KW - proton therapy

KW - Solidwater®

KW - treatment planning

UR - http://www.scopus.com/inward/record.url?scp=84908568354&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84908568354&partnerID=8YFLogxK

U2 - 10.1088/0031-9155/59/22/6859

DO - 10.1088/0031-9155/59/22/6859

M3 - Article

C2 - 25349982

AN - SCOPUS:84908568354

VL - 59

SP - 6859

EP - 6873

JO - Physics in Medicine and Biology

JF - Physics in Medicine and Biology

SN - 0031-9155

IS - 22

ER -