Fast and accurate proton therapy dose calculations
Simulating particle beam interactions includes identifying a set of n functions F 1 , F 2 , . . . , F n corresponding to a plurality of different physical aspects of a particle beam, performing simulations of each F i using a full physics model, selecting for each F i a distribution function f i that models relevant behavior and reducing computation of the full physics model for each F i by replacing F i with a distribution function f i . The computation reduction includes comparing a set of simulations wherein each f i replaces its respective F i to determine if relevant behavior is accurately modeled and selecting one of f i or F i for each n, for a Monte Carlo simulation based on a runtime and accuracy criteria.
1. A method for simulating particle beam interactions for reducing computation time of particle dosage applied to a lesion, comprising:
identifying a set of n functions F 1 , F 2 , . . . , F n corresponding to a plurality of different physical aspects of a particle beam;
performing simulations of each F i using a full physics model;
selecting for each F a distribution function f i that models relevant behavior of the particle beam;
reducing computation of the full physics model corresponding to each F i by replacing F i with the distribution function f i , wherein reducing computation includes:
computing the particle dosage using at least one f;
computing energy loss fluctuations in the particle beam;
comparing a set of simulations wherein each f i replaces Its respective F i to determine if relevant behavior is accurately modeled; and
selecting one of f i or F i for each n, for a Monte Carlo simulation based on a runtime and accuracy criteria;
treating the lesion based on the computed particle dosage and the computed energy loss fluctuations.
2. The method as recited in claim 1 , wherein the set of n functions includes one or more of particle energy levels, scattering angles, absorption material or particle density.
3. The method as recited in claim 1 , wherein the particle beam includes protons.
4. The method as recited in claim 3 , wherein the energy loss fluctuations are based on multipliers of cumulative distribution functions.
5. The method as recited in claim 1 , wherein reducing computation includes eliminating distributions outside of an operating range.
6. The method as recited in claim 1 , wherein reducing computations includes reducing complexity by generating a prediction model or an estimating model that maintains accuracy.
7. The method as recited in claim 1 , wherein the distribution function includes one or more of probability distribution functions and/or cumulative distribution functions.
8. The method as recited in claim 1 , wherein the distribution function f i includes one or more of an analytical function or table.
9. The method as recited in claim 1 , wherein selecting for each F i a distribution function f i includes selecting f i based on a computation system being used.
10. The method as recited in claim 1 , further comprising controlling a radiation source in accordance with a reduced physics model including at least one f i .
11. A non-transitory computer readable storage medium comprising a computer readable program for simulating particle beam interactions for reducing computation time of particle dosage applied to a lesion, wherein the computer readable program when executed on a computer causes the computer to perform steps of:
identifying a set of n functions F i , F 2 , . . . , F n corresponding to a plurality of different physical aspects of a particle beam;
performing simulations of each F i using a full physics model;
selecting for each F i a distribution function f that models relevant behavior of the particle beam;
reducing computation of the full physics model corresponding to each F i by replacing F i with the distribution function f i , wherein reducing computation includes:
computing the particle dosage using at least one f i ;
computing energy loss fluctuations in the particle beam; and
comparing a set of simulations wherein each f i replaces its respective F i to determine if relevant behavior is accurately modeled; and
selecting one of f i or F i for each n, for a Monte Carlo simulation based on a runtime and accuracy criteria; and
treating the lesion based on the computed particle dosage and the computed energy loss fluctuations.
12. The non-transitory computer readable storage medium as recited in claim 11 , wherein the set of n functions includes one or more of particle energy levels, scattering angles, absorption material or particle density.
13. The non-transitory computer readable storage medium as recited in claim 11 , wherein the particle beam includes protons.
14. The non-transitory computer readable storage medium as recited in claim 13 , wherein the energy loss fluctuations are based on multipliers of cumulative distribution functions.
15. The non-transitory computer readable storage medium as recited in claim 11 , wherein reducing computation includes eliminating distributions outside of an operating range.
16. The non-transitory computer readable storage medium as recited in claim 11 , wherein reducing computations includes reducing complexity by generating a prediction model or an estimating model that maintains accuracy.
17. The non-transitory computer readable storage medium as recited in claim 11 , wherein the distribution function includes one or more of probability distribution functions and/or cumulative distribution functions.
18. The non-transitory computer readable storage medium as recited in claim 11 , wherein the distribution function f i includes one or more of an analytical function or table.
19. The non-transitory computer readable storage medium as recited in claim 11 , wherein selecting for each F i a distribution function f i includes selecting f i based on a computation system being used.
20. The non-transitory computer readable storage medium as recited in claim 11 , further comprising controlling a radiation source in accordance with a reduced physics model including at least one f i .
21. A system for simulating particle beam interactions for reducing computation time of therapy-doses particle dosage applied to a lesion, comprising:
a dose computation module configured in non-transitory stored media to reduce computation of a full physics model corresponding to a set of n functions F 1 , F 2 , . . . , F n for a plurality of different physical aspects of a particle beam, wherein at least one F i is replaced with a distribution function f i to reduce computation by computing the particle dosage using, at least one f i by computing energy loss fluctuations in the particle beam, where f i models relevant behavior of the particle beam;
a simulation and optimization module configured to compare a set of simulations wherein each f i replaces its respective F i to determine if relevant behavior is accurately modeled and to selecting one of f i or F i for each n, for a Monte Carlo simulation based on a runtime and accuracy criteria; and
treat the lesion based on the computed particle dosage and the computed energy loss fluctuations from the dose computation module.
22. The system as recited in claim 21 , wherein the set of n functions includes one or more of particle energy levels, scattering angles, absorption material or particle density.
23. The system as recited in claim 21 , wherein the particle beam includes protons.
24. The system as recited in claim 21 , wherein the energy loss fluctuations are based on multipliers of cumulative distribution functions.
25. The system as recited in claim 21 , wherein the simulation and optimization module selects f i based on a computation system being used.