A fast GPU-based Monte Carlo simulation of proton transport with detailed modeling of nonelastic interactions.

Abstract:

PURPOSE:Very fast Monte Carlo (MC) simulations of proton transport have been implemented recently on graphics processing units (GPUs). However, these MCs usually use simplified models for nonelastic proton-nucleus interactions. Our primary goal is to build a GPU-based proton transport MC with detailed modeling of elastic and nonelastic proton-nucleus collisions. METHODS:Using the cuda framework, the authors implemented GPU kernels for the following tasks: (1) simulation of beam spots from our possible scanning nozzle configurations, (2) proton propagation through CT geometry, taking into account nuclear elastic scattering, multiple scattering, and energy loss straggling, (3) modeling of the intranuclear cascade stage of nonelastic interactions when they occur, (4) simulation of nuclear evaporation, and (5) statistical error estimates on the dose. To validate our MC, the authors performed (1) secondary particle yield calculations in proton collisions with therapeutically relevant nuclei, (2) dose calculations in homogeneous phantoms, (3) recalculations of complex head and neck treatment plans from a commercially available treatment planning system, and compared with (GEANT)4.9.6p2/TOPAS. RESULTS:Yields, energy, and angular distributions of secondaries from nonelastic collisions on various nuclei are in good agreement with the (GEANT)4.9.6p2 Bertini and Binary cascade models. The 3D-gamma pass rate at 2%-2 mm for treatment plan simulations is typically 98%. The net computational time on a NVIDIA GTX680 card, including all CPU-GPU data transfers, is ∼ 20 s for 1 × 10(7) proton histories. CONCLUSIONS:Our GPU-based MC is the first of its kind to include a detailed nuclear model to handle nonelastic interactions of protons with any nucleus. Dosimetric calculations are in very good agreement with (GEANT)4.9.6p2/TOPAS. Our MC is being integrated into a framework to perform fast routine clinical QA of pencil-beam based treatment plans, and is being used as the dose calculation engine in a clinically applicable MC-based IMPT treatment planning system. The detailed nuclear modeling will allow us to perform very fast linear energy transfer and neutron dose estimates on the GPU.

journal_name

Med Phys

journal_title

Medical physics

authors

Wan Chan Tseung H,Ma J,Beltran C

doi

10.1118/1.4921046

subject

Has Abstract

pub_date

2015-06-01 00:00:00

pages

2967-78

issue

6

eissn

0094-2405

issn

2473-4209

journal_volume

42

pub_type

杂志文章
  • Portal dose image (PDI) prediction for dosimetric treatment verification in radiotherapy. I. An algorithm for open beams.

    abstract::A method is presented for calculation of transmission functions for high energy photon beams through patients. These functions are being used in our clinic for prediction of portal dose images (PDIs) which are compared with PDIs measured with an electronic portal imaging device (EPID). The calculations are based on th...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.598293

    authors: Pasma KL,Heijmen BJ,Kroonwijk M,Visser AG

    更新日期:1998-06-01 00:00:00

  • Maximum likelihood estimation of proton irradiated field and deposited dose distribution.

    abstract::In proton therapy, it is important to evaluate the field irradiated with protons and the deposited dose distribution in a patient's body. Positron emitters generated through fragmentation reactions of target nuclei can be used for this purpose. By detecting the annihilation gamma rays from the positron emitters, the a...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.2712572

    authors: Inaniwa T,Kohno T,Yamagata F,Tomitani T,Sato S,Kanazawa M,Kanai T,Urakabe E

    更新日期:2007-05-01 00:00:00

  • Systematic calibration of an integrated x-ray and optical tomography system for preclinical radiation research.

    abstract:PURPOSE:The cone beam computed tomography (CBCT) guided small animal radiation research platform (SARRP) has been developed for focal tumor irradiation, allowing laboratory researchers to test basic biological hypotheses that can modify radiotherapy outcomes in ways that were not feasible previously. CBCT provides exce...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4914860

    authors: Yang Y,Wang KK,Eslami S,Iordachita II,Patterson MS,Wong JW

    更新日期:2015-04-01 00:00:00

  • Evaluating the impact of extended field-of-view CT reconstructions on CT values and dosimetric accuracy for radiation therapy.

    abstract:PURPOSE:Wide bore CT scanners use extended field-of-view (eFOV) reconstruction algorithms to attempt to recreate tissue truncated due to large patient habitus. Radiation therapy planning systems rely on accurate CT numbers in order to correctly plan and calculate radiation dose. This study looks at the impact of eFOV r...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1002/mp.13299

    authors: Cheung JP,Shugard E,Mistry N,Pouliot J,Chen J

    更新日期:2019-02-01 00:00:00

  • Temporal stability assessment in shear wave elasticity images validated by deep learning neural network for chronic liver disease fibrosis stage assessment.

    abstract:PURPOSE:To automatically detect and isolate areas of low and high stiffness temporal stability in shear wave elastography (SWE) image sequences and define their impact in chronic liver disease (CLD) diagnosis improvement by means of clinical examination study and deep learning algorithm employing convolutional neural n...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1002/mp.13521

    authors: Gatos I,Tsantis S,Spiliopoulos S,Karnabatidis D,Theotokas I,Zoumpoulis P,Loupas T,Hazle JD,Kagadis GC

    更新日期:2019-05-01 00:00:00

  • State-of-the-art in two-dimensional ultrasonic transducer array technology.

    abstract::Attention is given to means of sensing ultrasonic energy distributions over an area. Under the restraints appropiate to real-time imaging of deep body organs, piezoelectric arrays offer the most promising method. Adaptation of integrated circuit techniques to array assembly permits very large arrays of small elements ...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.594247

    authors: Maginness MG,Plummer JD,Beaver WL,Meindl JD

    更新日期:1976-09-01 00:00:00

  • A residual correction method for high-resolution PET reconstruction with application to on-the-fly Monte Carlo based model of positron range.

    abstract:PURPOSE:The quality of tomographic images is directly affected by the system model being used in image reconstruction. An accurate system matrix is desirable for high-resolution image reconstruction, but it often leads to high computation cost. In this work the authors present a maximum a posteriori reconstruction algo...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.3284980

    authors: Fu L,Qi J

    更新日期:2010-02-01 00:00:00

  • Dynamic field-of-view imaging to increase temporal resolution in the early phase of contrast media uptake in breast DCE-MRI: A feasibility study.

    abstract:PURPOSE:To increase diagnostic accuracy of breast MRI by increasing temporal resolution and more accurately sampling the early kinetics of contrast media uptake. We tested the feasibility of accelerating bilateral breast DCE-MRI by reducing the FOV, allowing aliasing, and unfolding the resulting images. METHODS:Previo...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1002/mp.12747

    authors: Pineda FD,Easley TO,Karczmar GS

    更新日期:2018-03-01 00:00:00

  • Design and experimental testing of air slab caps which convert commercial electron diodes into dual purpose, correction-free diodes for small field dosimetry.

    abstract:PURPOSE:Two diodes which do not require correction factors for small field relative output measurements are designed and validated using experimental methodology. This was achieved by adding an air layer above the active volume of the diode detectors, which canceled out the increase in response of the diodes in small f...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4894728

    authors: Charles PH,Cranmer-Sargison G,Thwaites DI,Kairn T,Crowe SB,Pedrazzini G,Aland T,Kenny J,Langton CM,Trapp JV

    更新日期:2014-10-01 00:00:00

  • SU-E-T-404: A Prototype Program for Analyzing 4D Image Guidance Shifts for Lung SBRT.

    abstract:PURPOSE:A prototype in-house program was developed to assist deriving image guidance shifts clinically using our 4DCT method for lung SBRT. METHODS:Our in-house prototype program was implemented in MATLAB. The 4D DICOM CT dataset of interest is imported by the program from AW and subsequently accepts input from the us...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4735493

    authors: Wang B,Sarkar V,Busselberg C,Rassiah-Szegedi P,Zhao H,Huang Y,Szegedi M,Salter B

    更新日期:2012-06-01 00:00:00

  • Practical cone-beam lambda tomography.

    abstract::As a potentially important technology for medical x-ray computed tomography (CT), lambda tomography (LT) is to reconstruct a gradient-like image only from local projection data. Based on our recently derived exact fan-beam LT formula, [H. Y. Gu and G. Wang, Int. J. Biomed. Imaging 2006(1), 1-9 (2006)] here we propose ...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.2348767

    authors: Yu H,Ye Y,Wang G

    更新日期:2006-10-01 00:00:00

  • Four-dimensional inverse treatment planning with inclusion of implanted fiducials in IMRT segmented fields.

    abstract::The purpose of this study is to develop a 4D inverse planning strategy capable of controlling the appearance of the implanted fiducial(s) in segmented IMRT fields for cine MV or combined MV/kV image-guided IMRT. This work is focused on enhancing the visibility of the implanted fiducial(s) in 4D IMRT inverse planning, ...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.3121425

    authors: Ma Y,Lee L,Keshet O,Keall P,Xing L

    更新日期:2009-06-01 00:00:00

  • Interleaf leakage for 5 and 10 mm dynamic multileaf collimation systems incorporating patient motion.

    abstract::Use of dynamic multileaf collimation (DMLC) for intensity modulated radiation therapy (IMRT) is accelerating. Delivery systems have the ailment of interleaf leakage (IL). This is compounded by the inefficiency of IMRT delivery, estimated to be a factor of 5 for DMLC. With IL on the order of 4%, it is possible to deliv...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.1388216

    authors: Klein EE,Low DA

    更新日期:2001-08-01 00:00:00

  • Measurement of primary bremsstrahlung spectrum from an 8-MeV linear accelerator.

    abstract::The bremsstrahlung spectrum from an 8-MeV linear accelerator has been measured using a NaI(T1) spectrometer system. The spectrum shows a low-energy cutoff at 0.4 MeV and the maximum photon energy to be approximately 6% greater than the nominal energy. The maximum emission of energy fluence was 1.6 and 1.8 MeV for meas...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.594221

    authors: Levy LB,Waggener RG,Wright AE

    更新日期:1976-05-01 00:00:00

  • Eye tracking and gating system for proton therapy of orbital tumors.

    abstract:PURPOSE:A new motion-based gated proton therapy for the treatment of orbital tumors using real-time eye-tracking system was designed and evaluated. METHODS:We developed our system by image-pattern matching, using a normalized cross-correlation technique with LabVIEW 8.6 and Vision Assistant 8.6 (National Instruments, ...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4729708

    authors: Shin D,Yoo SH,Moon SH,Yoon M,Lee SB,Park SY

    更新日期:2012-07-01 00:00:00

  • Monte Carlo study of ionization chamber magnetic field correction factors as a function of angle and beam quality.

    abstract:PURPOSE:To use EGSnrc Monte Carlo simulations for magnetic field dosimetry to determine optimal measurement orientations, calculate beam quality conversion factors for 32 cylindrical and three parallel-plate (PP) ion chambers, evaluate the beam quality and angular dependence of these factors, and examine the magnetic f...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1002/mp.12716

    authors: Malkov VN,Rogers DWO

    更新日期:2018-02-01 00:00:00

  • Approximate 3D iterative reconstruction for SPECT.

    abstract::Compared with slice-by-slice approaches for SPECT reconstruction, three-dimensional iterative methods provide a more accurate physical model and an improved SPECT image. Clinical application of these methods, however, is limited primarily to their computational demands. This paper investigates the methods for approxim...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.598030

    authors: Gilland DR,Jaszczak RJ,Riauka TA,Coleman RE

    更新日期:1997-09-01 00:00:00

  • An accurate method for direct dual-energy calibration and decomposition.

    abstract::We propose the use of conic and cubic surface equations (surfaces of second and third order) to directly approximate the dual-energy equations (the integral equations for the dual-energy log-signal functions, i.e., the negative logarithms of the relative detector signals, considered as functions of the basis-material ...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.596512

    authors: Cardinal HN,Fenster A

    更新日期:1990-05-01 00:00:00

  • Effects of protocol and obesity on dose conversion factors in adult body CT.

    abstract:PURPOSE:In computed tomography (CT), organ dose, effective dose, and risk index can be estimated from volume-weighted CT dose index (CTDI(vol)) or dose-length product (DLP) using conversion coefficients. Studies have investigated how these coefficients vary across scanner models, scan parameters, and patient size. Howe...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4754584

    authors: Li X,Samei E,Williams CH,Segars WP,Tward DJ,Miller MI,Ratnanather JT,Paulson EK,Frush DP

    更新日期:2012-11-01 00:00:00

  • Performance evaluation of a diode array for enhanced dynamic wedge dosimetry.

    abstract::The performance of a diode array (Profiler) was evaluated by comparing its enhanced dynamic wedge (EDW) profiles measured at various depths with point measurements using a 0.03 cm3 ionization chamber on a commercial linear accelerator. The Profiler, which covers a 22.5 cm width, was used to measure larger field widths...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.598019

    authors: Zhu TC,Ding L,Liu CR,Palta JR,Simon WE,Shi J

    更新日期:1997-07-01 00:00:00

  • Assessment of image registration accuracy in three-dimensional transrectal ultrasound guided prostate biopsy.

    abstract:PURPOSE:Prostate biopsy, performed using two-dimensional (2D) transrectal ultrasound (TRUS) guidance, is the clinical standard for a definitive diagnosis of prostate cancer. Histological analysis of the biopsies can reveal cancerous, noncancerous, or suspicious, possibly precancerous, tissue. During subsequent biopsy s...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.3298010

    authors: Karnik VV,Fenster A,Bax J,Cool DW,Gardi L,Gyacskov I,Romagnoli C,Ward AD

    更新日期:2010-02-01 00:00:00

  • High resolution 2D dose measurement device based on a few long scintillating fibers and tomographic reconstruction.

    abstract:PURPOSE:Patient-specific QA of highly conformal radiotherapy treatments are usually conducted using 2D or 3D dosimetry of the incident dose distribution in a water-equivalent phantom. However, dosimeters typically used for this task usually lack in either spatial resolution or dose accuracy. The purpose of this work is...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.4736526

    authors: Goulet M,Archambault L,Beaulieu L,Gingras L

    更新日期:2012-08-01 00:00:00

  • Robust optimization in IMPT using quadratic objective functions to account for the minimum MU constraint.

    abstract:PURPOSE:Currently, in clinical practice of intensity-modulated proton therapy (IMPT), the influence of the minimum monitor unit (MU) constraint is taken into account through postprocessing after the optimization is completed. This may degrade the plan quality and plan robustness. This study aims to mitigate the impact ...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1002/mp.12677

    authors: Shan J,An Y,Bues M,Schild SE,Liu W

    更新日期:2018-01-01 00:00:00

  • Microdosimetric study for secondary neutrons in phantom produced by a 290 MeV/nucleon carbon beam.

    abstract::Absorbed doses from main charged-particle beams and charged-particle fragments have been measured with high accuracy for particle therapy, but there are few reports for doses from neutron components produced as fragments. This study describes the measurements on neutron doses produced by carbon beams; microdosimetric ...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.2767933

    authors: Endo S,Tanaka K,Takada M,Onizuka Y,Miyahara N,Sato T,Ishikawa M,Maeda N,Hayabuchi N,Shizuma K,Hoshi M

    更新日期:2007-09-01 00:00:00

  • Minimax optimization for handling range and setup uncertainties in proton therapy.

    abstract:PURPOSE:Intensity modulated proton therapy (IMPT) is sensitive to errors, mainly due to high stopping power dependency and steep beam dose gradients. Conventional margins are often insufficient to ensure robustness of treatment plans. In this article, a method is developed that takes the uncertainties into account duri...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.3556559

    authors: Fredriksson A,Forsgren A,Hårdemark B

    更新日期:2011-03-01 00:00:00

  • Nine-degrees-of-freedom flexmap for a cone-beam computed tomography imaging device with independently movable source and detector.

    abstract:PURPOSE:Couch-mounted cone-beam computed tomography (CBCT) imaging devices with independently rotatable x-ray source and flat-panel detector arms for acquisitions of arbitrary regions of interest (ROI) have recently been introduced in image-guided radiotherapy (IGRT). This work analyzes mechanical limitations and gravi...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1002/mp.12033

    authors: Keuschnigg P,Kellner D,Fritscher K,Zechner A,Mayer U,Huber P,Sedlmayer F,Deutschmann H,Steininger P

    更新日期:2017-01-01 00:00:00

  • A generic TG-186 shielded applicator for commissioning model-based dose calculation algorithms for high-dose-rate 192 Ir brachytherapy.

    abstract:PURPOSE:A joint working group was created by the American Association of Physicists in Medicine (AAPM), the European Society for Radiotherapy and Oncology (ESTRO), and the Australasian Brachytherapy Group (ABG) with the charge, among others, to develop a set of well-defined test case plans and perform calculations and ...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1002/mp.12459

    authors: Ma Y,Vijande J,Ballester F,Tedgren ÅC,Granero D,Haworth A,Mourtada F,Fonseca GP,Zourari K,Papagiannis P,Rivard MJ,Siebert FA,Sloboda RS,Smith R,Chamberland MJP,Thomson RM,Verhaegen F,Beaulieu L

    更新日期:2017-11-01 00:00:00

  • Comparison of planned dose on different CT image sets to four-dimensional Monte Carlo dose recalculation using the patient's actual breathing trace for lung stereotactic body radiation therapy.

    abstract:PURPOSE:The need for four-dimensional (4D) treatment planning becomes indispensable when it comes to radiation therapy for moving tumors in the thoracic and abdominal regions. The primary purpose of this study is to combine the actual breathing trace during each individual treatment fraction with the Linac's log file i...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1002/mp.13579

    authors: Freislederer P,von Münchow A,Kamp F,Heinz C,Gerum S,Corradini S,Söhn M,Reiner M,Roeder F,Floca R,Alber M,Belka C,Parodi K

    更新日期:2019-07-01 00:00:00

  • Reduction of motion blurring artifacts using respiratory gated CT in sinogram space: a quantitative evaluation.

    abstract::Techniques have been developed for reducing motion blurring artifacts by using respiratory gated computed tomography (CT) in sinogram space and quantitatively evaluating the artifact reduction. A synthetic sinogram was built from multiple scans intercepting a respiratory gating window. A gated CT image was then recons...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1118/1.2074187

    authors: Lu W,Parikh PJ,Hubenschmidt JP,Politte DG,Whiting BR,Bradley JD,Mutic S,Low DA

    更新日期:2005-11-01 00:00:00

  • A unified material decomposition framework for quantitative dual- and triple-energy CT imaging.

    abstract:PURPOSE:Many clinical applications depend critically on the accurate differentiation and classification of different types of materials in patient anatomy. This work introduces a unified framework for accurate nonlinear material decomposition and applies it, for the first time, in the concept of triple-energy CT (TECT)...

    journal_title:Medical physics

    pub_type: 杂志文章

    doi:10.1002/mp.12933

    authors: Zhao W,Vernekohl D,Han F,Han B,Peng H,Yang Y,Xing L,Min JK

    更新日期:2018-07-01 00:00:00