title: Evaluation of uncertainty-based stopping criteria for monte carlo calculations of intensity-modulated radiotherapy and arc therapy patient dose distributions.
abstract: OBJECTIVE
  To formulate uncertainty-based stopping criteria for Monte Carlo (MC) calculations of intensity-modulated radiotherapy and intensity-modulated arc therapy patient dose distributions and evaluate their influence on MC simulation times and dose characteristics.
  
  
  METHODS
  For each structure of interest, stopping criteria were formulated as follows: sigma(rel) <or=sigma(rel,tol) or Dsigma(rel) <or=D(lim)sigma(rel,tol) within >or=95% of the voxels, where sigma(rel) represents the relative statistical uncertainty on the estimated dose, D. The tolerated uncertainty (sigma(rel,tol)) was 2%. The dose limit (D(lim)) equaled the planning target volume (PTV) prescription dose or a dose value related to the organ at risk (OAR) planning constraints. An intensity-modulated radiotherapy-lung, intensity-modulated radiotherapy-ethmoid sinus, and intensity-modulated arc therapy-rectum patient case were studied. The PTV-stopping criteria-based calculations were compared with the PTV+OAR-stopping criteria-based calculations.
  
  
  RESULTS
  The MC dose distributions complied with the PTV-stopping criteria after 14% (lung), 21% (ethmoid), and 12% (rectum) of the simulation times of a 100 million histories reference calculation, and increased to 29%, 44%, and 51%, respectively, by the addition of the OAR-stopping criteria. Dose-volume histograms corresponding to the PTV-stopping criteria, PTV+OAR-stopping criteria, and reference dose calculations were indiscernible. The median local dose differences between the PTV-stopping criteria and the reference calculations amounted to 1.4% (lung), 2.1% (ethmoid), and 2.5% (rectum).
  
  
  CONCLUSIONS
  For the patient cases studied, the MC calculations using PTV-stopping criteria only allowed accurate treatment plan evaluation. The proposed stopping criteria provided a flexible tool to assist MC patient dose calculations. The structures of interest and appropriate values of sigma(rel,tol) and D(lim) should be selected for each patient individually according to the clinical treatment planning goals.
authors.affiliation: Department of Medical Physics, Ghent University, Gent, Belgium. Barbara.Vanderstraeten@UGent.be
authors.forename: Barbara
authors.initials: B
authors.lastname: Vanderstraeten
authors.forename: Ana Maria Luiza
authors.initials: AM
authors.lastname: Olteanu
authors.forename: Nick
authors.initials: N
authors.lastname: Reynaert
authors.forename: Antonio
authors.initials: A
authors.lastname: Leal
authors.forename: Wilfried
authors.initials: W
authors.lastname: De Neve
authors.forename: Hubert
authors.initials: H
authors.lastname: Thierens
mesh_terms: D005005:Ethmoid Sinus; D006801:Humans; D008175:Lung Neoplasms / Q000473:pathology / Q000532:radiotherapy; D009010:Monte Carlo Method; D009369:Neoplasms / Q000473:pathology / Q000532:radiotherapy*; D010255:Paranasal Sinus Neoplasms / Q000473:pathology / Q000532:radiotherapy; D011832:Radiation Injuries / Q000517:prevention & control; D011879:Radiotherapy Dosage; D050397:Radiotherapy, Intensity-Modulated / Q000379:methods*; D012004:Rectal Neoplasms / Q000473:pathology / Q000532:radiotherapy; D047368:Tumor Burden; D035501:Uncertainty