A particle track-repeating algorithm for proton beam dose calculation

A particle track-repeating algorithm has been developed for proton beam dose calculation for radiotherapy. Monoenergetic protons with 250 MeV kinetic energy were simulated in an infinite water phantom using the GEANT3 Monte Carlo code. The changes in location, angle and energy for every transport step and the energy deposition along the track were recorded for the primary protons and all secondary particles. When calculating dose for a patient with a realistic proton beam, the pre-generated particle tracks were repeated in the patient geometry consisting of air, soft tissue and bone. The medium and density for each dose scoring voxel in the patient geometry were derived from patient CT data. The starting point, at which a proton track was repeated, was determined according to the incident proton energy. Thus, any protons with kinetic energy less than 250 MeV can be simulated. Based on the direction of the incident proton, the tracks were first rotated and for the subsequent steps, the scattering angles were simply repeated for air and soft tissue but adjusted properly based on the scattering power for bone. The particle step lengths were adjusted based on the density for air and soft tissue and also on the stopping powers for bone while keeping the energy deposition unchanged in each step. The difference in nuclear interactions and secondary particle generation between water and these materials was ignored. The algorithm has been validated by comparing the dose distributions in uniform water and layered heterogeneous phantoms with those calculated using the GEANT3 code for 120, 150, 180 and 250 MeV proton beams. The differences between them were within 2%. The new algorithm was about 13 times faster than the GEANT3 Monte Carlo code for a uniform phantom geometry and over 700 times faster for a heterogeneous phantom geometry.     

4D Monte Carlo simulation of proton beam scanning: modelling of variations in time and space to study the interplay between scanning pattern and time-dependent patient geometry

When dosimetric effects in time-dependent geometries are studied, usually either the results of individual three-dimensional (3D) calculations are combined or probability-based approaches are applied. These methods may become cumbersome and time-consuming if high time resolution is required or if the geometry is complex. Furthermore, it is difficult to study double-dynamic systems, e.g., to investigate the influence of time-dependent beam delivery (i.e., magnetically moving beam spots in proton beam scanning) on the dose deposition in a moving target. We recently introduced the technique of 4D Monte Carlo dose calculation to model continuously changing geometries. In intensity modulated proton therapy, dose is delivered by individual pristine Bragg curves. Dose spots are positioned in the patient by varying magnetic field and beam energy. If the movement of these dose spots occurs during significant respiratory motion, interplay effects can take place. Because of the inhomogeneity of individual subfields, the consequences of motion can be more severe than in conventional proton therapy. We demonstrate how the technique of 4D Monte Carlo can be used to study interplay effects in proton beam scanning. Time-dependent beam delivery to a changing patient geometry is simulated in a single 4D dose calculation. Interplay effects between respiratory motion and beam scanning speed are demonstrated.     

Beam characteristics of upper and lower physical wedge systems of Varian accelerators

The beam characteristics of a dual physical wedge system, upper and lower, for Varian accelerators are studied over the energy range 6-18 MV. Wedge factors for both systems are measured in a water phantom as a function of field size, depth and source-to-wedge (SWD) distance. Our results indicate that apart from their physical differences, dosimetrically, the two wedge systems have <2% difference in central axis percentage depth dose beyond the build-up region. The lower wedge central axis percentage depth dose is consistently lower than that of the corresponding upper wedge, with the effect more pronounced for large field sizes. The wedge profiles are identical within 2% for all field sizes, depths and energies. The wedge factors for both wedge systems are also within 2% for all field sizes and depths for both 6 and 15 MV photons and slightly higher for the 18 MV beam and 45 degrees-60 degrees wedge angle. The wedge factor variation with SWD reveals an interesting fact that thinner wedges (15 degrees and 30 degrees) result in a higher surface dose in the central axis region than thicker wedges. As the SWD increases beyond 80 cm, the reverse is true, i.e. thicker wedges produce higher surface dose than thinner wedges. It is also verified that the wedge factor at any depth and for any field size can be calculated from the wedged and open field central axis percentage depth dose, and the wedge factor at dmax, resulting in nearly 44% reduction in water phantom scanning and 80% reduction in point measurements during commissioning.     

Nuclear interactions in proton therapy: dose and relative biological effect distributions originating from primary and secondary particles

The dose distribution delivered in charged particle therapy is due to both primary and secondary particles. The secondaries, originating from non-elastic nuclear interactions, are of interest for three reasons. First, if fast Monte Carlo treatment planning is envisaged, the question arises whether all nuclear interaction products deliver a significant contribution to the total dose and, hence, need to be tracked. Second, there could be an enhanced relative biological effectiveness (RBE) due to low energy and/or heavy secondaries. Third, neutrons originating from nuclear interactions may deliver dose outside the target volume. The particle yield from different nuclear interaction channels as a function of proton penetration depth was studied theoretically for different proton beam energies. Three-dimensional dose distributions from primary and secondary particles were simulated for an unmodulated 160 MeV proton beam with and without including a slice of bone material and for a spread-out Bragg peak (SOBP) of 3 x 3 x 3 cm3 in water. Secondary protons deliver up to 10% of the total dose proximal to the Bragg peak of an unmodulated proton beam and they affect the flatness of the SOBP. Furthermore, they cause a dose build-up due to forward emission of secondary particles from nuclear interactions. The dose deposited by d, t, 3He and alpha-particles was found to contribute less than 0.1% of the total dose. The dose distal to the target volume caused by liberated neutrons was studied for four proton beam energies in the range of 160-250 MeV and found to be below 0.05% (2 cm distal to SOBP) of the prescribed target dose for a 3 x 3 x 3 cm3 target. RBE values relative to 60Co were calculated proximal to and within the SOBP. The RBE proximal to the Bragg peak (100% dose) is influenced by secondary particles (mainly protons and a-particles) with a strong dose dependency resulting in RBE values up to 1.2 (2 Gy; inactivation of V79). Depending on the endpoint considered, secondary particles cause a shift in RBE by up to 8% at 2 Gy. In contrast, the RBE in the Bragg peak is almost entirely determined by primary protons due to a decreasing secondary particle fluence with depth. RBE values up to 1.3 (2 Gy; inactivation of V79) at 1 cm distal to the Bragg peak maximum were found. The inactivations of human skin fibroblasts and mouse lymphoma cells were also analysed and reveal a substantial tissue dependency of the total RBE. The outcome of this study shows that elevated RBE values occur not only at the distal edge of the SOBP. Although the variations are modest, and in most cases might have no observable clinical effect, they might have to be considered in certain treatment situations. The biological effect downstream of the target caused by neutrons was analysed using a radiation quality factor of 10. The biological dose was found to be below 0.5% of the prescribed target dose (for a 3 x 3 x 3 cm3 SOBP) but depends on the size of the SOBP. This dose should not be significant with respect to late effects, e.g. cancer induction.     

光子辐射输运中次级效应对辐射剂量深度分布的影响

目的:模拟光子输运的过程,记录各相互作用和次级粒子对剂量计算的贡献,总结分析其对剂量贡献的大小.方法:PENELOPE程序包提供了模拟光子和电子输运的基本MC模块.基于所关心的物理问题本文对PENELOPE程序包进行二次编程,以在模拟过程中追踪光子输运详细过程,记录各相互作用及次级粒子对剂量的贡献.结果:首先研究在相同照射条件下,4种能量(10 keV,100keV,1 MeV,10MeV)的光子产生的中心轴剂量分布,次级粒子的软碰撞和硬碰撞产生的中心轴剂量分布,以及各级次级粒子的中心轴剂量分布;然后研究在相同照射条件下,4种能量(30keV,40keV,50keV,60keV)的光子产生的次级康普顿效应和次级光电效应对中心轴剂量分布.结论:不同能量下,次级电子软碰撞对于中心轴剂量的贡献起主要作用,次级光电效应对中心轴剂量的贡献随能量的增加而减小,而第一代次级粒子对于中心轴剂量的贡献大于其它代粒子的贡献.     

Simple beam models for Monte Carlo photon beam dose calculations in radiotherapy

Monte Carlo (code GEANT) produced 6 and 15 MV phase space (PS) data were used to define several simple photon beam models. For creating the PS data the energy of starting electrons hitting the target was tuned to get correct depth dose data compared to measurements. The modeling process used the full PS information within the geometrical boundaries of the beam including all scattered radiation of the accelerator head. Scattered radiation outside the boundaries was neglected. Photons and electrons were assumed to be radiated from point sources. Four different models were investigated which involved different ways to determine the energies and locations of beam particles in the output plane. Depth dose curves, profiles, and relative output factors were calculated with these models for six field sizes from 5x5 to 40x40cm2 and compared to measurements. Model 1 uses a photon energy spectrum independent of location in the PS plane and a constant photon fluence in this plane. Model 2 takes into account the spatial particle fluence distribution in the PS plane. A constant fluence is used again in model 3, but the photon energy spectrum depends upon the off axis position. Model 4, finally uses the spatial particle fluence distribution and off axis dependent photon energy spectra in the PS plane. Depth dose curves and profiles for field sizes up to 10x10cm2 were not model sensitive. Good agreement between measured and calculated depth dose curves and profiles for all field sizes was reached for model 4. However, increasing deviations were found for increasing field sizes for models 1-3. Large deviations resulted for the profiles of models 2 and 3. This is due to the fact that these models overestimate and underestimate the energy fluence at large off axis distances. Relative output factors consistent with measurements resulted only for model 4.     

A beam source model for scanned proton beams

A beam source model, i.e. a model for the initial phase space of the beam, for scanned proton beams has been developed. The beam source model is based on parameterized particle sources with characteristics found by fitting towards measured data per individual beam line. A specific aim for this beam source model is to make it applicable to the majority of the various proton beam systems currently available or under development, with the overall purpose to drive dose calculations in proton beam treatment planning. The proton beam phase space is characterized by an energy spectrum, radial and angular distributions and deflections for the non-modulated elementary pencil beam. The beam propagation through the scanning magnets is modelled by applying experimentally determined focal points for each scanning dimension. The radial and angular distribution parameters are deduced from measured two-dimensional fluence distributions of the elementary beam in air. The energy spectrum is extracted from a depth dose distribution for a fixed broad beam scan pattern measured in water. The impact of a multi-slab range shifter for energy modulation is calculated with an own Monte Carlo code taking multiple scattering, energy loss and straggling, non-elastic and elastic nuclear interactions in the slab assembly into account. Measurements for characterization and verification have been performed with the scanning proton beam system at The Svedberg Laboratory in Uppsala. Both in-air fluence patterns and dose points located in a water phantom were used. For verification, dose-in-water was calculated with the Monte Carlo code GEANT 3.21 instead of using a clinical dose engine with approximations of its own. For a set of four individual pencil beams, both with the full energy and range shifted, 96.5% (99.8%) of the tested dose points satisfied the 1%/1 mm (2%/2 mm) gamma criterion.     

On the field-size dependence of relative output from a linear accelerator.

The radiation output in air on the central axis of a linac photon beam has been modeled as the sum of two components. These are a point source representing radiation direct from the target and a distributed source representing scatter in the flattening filter and primary collimator. By fitting only two parameters, the ratio of the two components for a 20 x 20 field and a width parameter for the distributed source this semi-empirical model describes the relative outputs of square, symmetric rectangular, and asymmetric rectangular fields with an average error of 0.25% for the field sizes studied.     

Using a photon phase-space source for convolution/superposition dose calculations in radiation therapy

For a given linac design, the dosimetric characteristics of a photon beam are determined uniquely by the energy and radial distributions of the electron beam striking the x-ray target. However, in the usual commissioning of a beam from measured data, a large number of variables can be independently tuned, making it difficult to derive a unique and self-consistent beam model. For example, the measured dosimetric penumbra in water may be attributed in various proportions to the lateral secondary electron range, the focal spot size and the transmission through the tips of a non-divergent collimator; the head-scatter component in the tails of the transverse profiles may not be easy to resolve from phantom scatter and head leakage; and the head-scatter tails corresponding to a certain extra-focal source model may not agree self-consistently with in-air output factors measured on the central axis. To reduce the number of adjustable variables in beam modelling, we replace the focal and extra-focal sources with a single phase-space plane scored just above the highest adjustable collimator in a EGS/BEAM simulation of the linac. The phase-space plane is then used as photon source in a stochastic convolution/superposition dose engine. A photon sampled from the uncollimated phase-space plane is first propagated through an arbitrary collimator arrangement and then interacted in the simulation phantom. Energy deposition kernel rays are then randomly issued from the interaction points and dose is deposited along these rays. The electrons in the phase-space file are used to account for electron contamination. 6 MV and 18 MV photon beams from an Elekta SL linac are used as representative examples. Except for small corrections for monitor backscatter and collimator forward scatter for large field sizes (<0.5% with <20 x 20 cm2 field size), we found that the use of a single phase-space photon source provides accurate and self-consistent results for both relative and absolute dose calculations.     

Initial beam size study for passive scatter proton therapy. II. Changes in delivered depth dose profiles

In passively scattered proton radiotherapy, a clinically useful treatment beam is produced by spreading a small proton "pencil beam" extracted from the accelerator to create both a uniform dose profile laterally and a uniform spread-out Bragg peak (SOBP) in depth. Lateral spreading and range modulation of the beam are accomplished using specially designed components within the treatment delivery nozzle. The purpose of this study was to determine how changes in the size of the initial proton pencil beam affect the delivery of dose with a passive scatter treatment nozzle. Monte Carlo calculations were used to study changes of the beam's in-air energy distribution at the exit of the nozzle and the central axis depth dose profiles in water resulting from changes in the incident beam size. Our results indicate that the width of the delivered SOBP decreases as the size of the initial beam increases.     

Commissioning 6 MV photon beams of a stereotactic radiosurgery system for Monte Carlo treatment planning

The goal of this work is to implement a beam commissioning procedure to generate a multiple source model using a set of standard measurement data for possible Monte Carlo treatment planning in the clinic for a Cyberknife stereotactic radiosurgery system. The required measurement data include the central axis depth dose curve (PDD), the dose profile at dmax(= 1.5 cm) of 60 mm cone at 80 cm source-to-surface distance (SSD), and the cone output factors for cones of 5 mm to 60 mm at 80 cm source-to-axis distance (SAD). The employed dual source model has the same structure as the one that has been studied in our previous work while most of the parameters of each source are extracted from the measurement data rather than the beam phase space. The energy spectra will be extracted from the central axis PDD, the fluence distributions will be deconvoluted from the dose profile at dmax, and the source distributions will be determined from the measured cone output factors. Monte Carlo dose calculations in various water phantoms have been performed to verify the beam commissioning procedure. The agreement between the measurements and the commissioning results was within 2%/1 mm for the central axis PDDs and the dose profiles at various depths when an IC-3 chamber was used and within 2% for the cone output factors for various collimator sizes of 5 to 60 mm. Largest difference (9.5%) was observed for the 7.5 mm cone when an IC-10 chamber was used. The large differences can be attributed to the volumetric averaging effect of the IC-10 chamber, whose dimension is comparable to the field of the small cones. The overall agreement between the measurements and the commissioning results is clinically acceptable, which implies that our commissioning tool is adequate for clinical applications of Monte Carlo dose calculations for the Cyberknife stereotactic radiosurgery system.     

A Monte Carlo pencil beam scanning model for proton treatment plan simulation using GATE/GEANT4

This work proposes a generic method for modeling scanned ion beam delivery systems, without simulation of the treatment nozzle and based exclusively on beam data library (BDL) measurements required for treatment planning systems (TPS). To this aim, new tools dedicated to treatment plan simulation were implemented in the Gate Monte Carlo platform. The method was applied to a dedicated nozzle from IBA for proton pencil beam scanning delivery. Optical and energy parameters of the system were modeled using a set of proton depth-dose profiles and spot sizes measured at 27 therapeutic energies. For further validation of the beam model, specific 2D and 3D plans were produced and then measured with appropriate dosimetric tools. Dose contributions from secondary particles produced by nuclear interactions were also investigated using field size factor experiments. Pristine Bragg peaks were reproduced with 0.7 mm range and 0.2 mm spot size accuracy. A 32 cm range spread-out Bragg peak with 10 cm modulation was reproduced with 0.8 mm range accuracy and a maximum point-to-point dose difference of less than 2%. A 2D test pattern consisting of a combination of homogeneous and high-gradient dose regions passed a 2%/2 mm gamma index comparison for 97% of the points. In conclusion, the generic modeling method proposed for scanned ion beam delivery systems was applicable to an IBA proton therapy system. The key advantage of the method is that it only requires BDL measurements of the system. The validation tests performed so far demonstrated that the beam model achieves clinical performance, paving the way for further studies toward TPS benchmarking. The method involves new sources that are available in the new Gate release V6.1 and could be further applied to other particle therapy systems delivering protons or other types of ions like carbon.     

Photonuclear dose calculations for high-energy photon beams from Siemens and Varian linacs

The dose from photon-induced nuclear particles (neutrons, protons, and alpha particles) generated by high-energy photon beams from medical linacs is investigated. Monte Carlo calculations using the MCNPX code are performed for three different photon beams from two different machines: Siemens 18 MV, Varian 15 MV, and Varian 18 MV. The linac head components are simulated in detail. The dose distributions from photons, neutrons, protons, and alpha particles are calculated in a tissue-equivalent phantom. Neutrons are generated in both the linac head and the phantom. This study includes (a) field size effects, (b) off-axis dose profiles, (c) neutron contribution from the linac head, (d) dose contribution from capture gamma rays, (e) phantom heterogeneity effects, and (f) effects of primary electron energy shift. Results are presented in terms of absolute dose distributions and also in terms of DER (dose equivalent ratio). The DER is the maximum dose from the particle (neutron, proton, or alpha) divided by the maximum photon dose, multiplied by the particle quality factor and the modulation scaling factor. The total DER including neutrons, protons, and alphas is about 0.66 cSv/Gy for the Siemens 18 MV beam (10 cm x 10 cm). The neutron DER decreases with decreasing field size while the proton (or alpha) DER does not vary significantly except for the 1 cm x 1 cm field. Both Varian beams (15 and 18 MV) produce more neutrons, protons, and alphas particles than the Siemens 18 MV beam. This is mainly due to their higher primary electron energies: 15 and 18.3 MeV, respectively, vs 14 MeV for the Siemens 18 MV beam. For all beams, neutrons contribute more than 75% of the total DER, except for the 1 cm x 1 cm field (approximately 50%). The total DER is 1.52 and 2.86 cSv/Gy for the 15 and 18 MV Varian beams (10 cm x 10 cm), respectively. Media with relatively high-Z elements like bone may increase the dose from heavy charged particles by a factor 4. The total DER is sensitive to primary electron energy shift. A Siemens 18 MV beam with 15 MeV (instead of 14 MeV) primary electrons would increase by 40% the neutron DER and by 210% the proton + alpha DER. Comparisons with measurements (neutron yields from different materials and neutron dose equivalent) are also presented. Using the NCRP risk assessment method, we found that the dose equivalent from leakage neutrons (at 50-cm off-axis distance) represent 1.1, 1.1, and 2.0% likelihood of fatal secondary cancer for a 70 Gy treatment delivered by the Siemens 18 MV, Varian 15 MV, and Varian 18 MV beams, respectively.     

Particle selection and beam collimation system for laser-accelerated proton beam therapy

In a laser-accelerated proton therapy system, the initial protons have broad energy and angular distributions, which are not suitable for direct therapeutic applications. A compact particle selection and collimation device is needed to deliver small pencil beams of protons with desired energy spectra. In this work, we characterize a superconducting magnet system that produces a desired magnetic field configuration to spread the protons with different energies and emitting angles for particle selection. Four magnets are set side by side along the beam axis; each is made of NbTi wires which carry a current density of approximately 10(5) A/cm2 at 4.2 K, and produces a magnetic field of approximately 4.4 T in the corresponding region. Collimation is applied to both the entrance and the exit of the particle selection system to generate a desired proton pencil beam. In the middle of the magnet system, where the magnetic field is close to zero, a particle selection collimator allows only the protons with desired energies to pass through for therapy. Simulations of proton transport in the presence of the magnetic field show that the selected protons have successfully refocused on the beam axis after passing through the magnetic field with the optimal magnet system. The energy spread for any given characteristic proton energy has been obtained. It is shown that the energy spread is a function of the magnetic field strength and collimator size and reaches the full width at half maximum of 25 MeV for 230 MeV protons. Dose distributions have also been calculated with the GEANT3 Monte Carlo code to study the dosimetric properties of the laser-accelerated proton beams for radiation therapy applications.     

Particle selection for laser-accelerated proton therapy feasibility study

In this paper we present calculations for the design of a particle selection system for laser-accelerated proton therapy. Laser-accelerated protons coming from a thin high-density foil have broad energy and angular spectra leading to dose distributions that cannot be directly used for therapeutic applications. Our solution to this problem is a compact particle selection and collimation device that delivers small pencil beams of protons with desired energy spectra. We propose a spectrometer-like particle selection and beam modulation system in which the magnetic field will be used to spread the protons spatially according to their energies and emitting angles. Subsequently, an aperture will be used to select the protons within a therapeutic window of energy (energy modulation). It will be shown that for the effective proton spatial differentiation, the primary collimation device should be used, which will collimate protons to the desired angular distribution and limit the spatial mixing of different energy protons once they have traveled through the magnetic system. Due to the angular proton distribution, the spatial mixing of protons of different energies will always be present and it will result in a proton energy spread with the width depending on the energy. For 250 MeV protons, the width (from the maximum to the minimum energy) is found to be 50 MeV for the magnetic field configuration used in our calculations. As the proton energy decreases, its energy width decreases as well, and for 80 MeV protons it equals 9 MeV. The presence of the energy width in the proton energy distribution will modify the depth dose curves needed for the energy modulation calculation. The matching magnetic field setup will ensure the refocusing of the selected protons and the final beam will be collimated by the secondary collimator. The calculations presented in this article show that the dose rate that the selection system can yield is on the order of D=260 Gy/min for a field size of 1 x 1 cm2.     

Photon dose calculation based on electron multiple-scattering theory: primary dose deposition kernels.

The transport of the secondary electrons resulting from high-energy photon interactions is essential to energy redistribution and deposition. In order to develop an accurate dose-calculation algorithm for high-energy photons, which can predict the dose distribution in inhomogeneous media and at the beam edges, we have investigated the feasibility of applying electron transport theory [Jette, Med. Phys. 15, 123 (1988)] to photon dose calculation. In particular, the transport of and energy deposition by Compton electron and electrons and positrons resulting from pair production were studied. The primary photons are treated as the source of the secondary electrons and positrons, which are transported through the irradiated medium using Gaussian multiple-scattering theory [Jette, Med. Phys. 15, 123 (1988)]. The initial angular and kinetic energy distribution(s) of the secondary electrons (and positrons) emanating from the photon interactions are incorporated into the transport. Due to different mechanisms of creation and cross-section functions, the transport of and the energy deposition by the electrons released in these two processes are studied and modeled separately based on first principles. In this article, we focus on determining the dose distribution for an individual interaction site. We define the Compton dose deposition kernel (CDK) or the pair-production dose deposition kernel (PDK) as the dose distribution relative to the point of interaction, per unit interaction density, for a monoenergetic photon beam in an infinite homogeneous medium of unit density. The validity of this analytic modeling of dose deposition was evaluated through EGS4 Monte Carlo simulation. Quantitative agreement between these two calculations of the dose distribution and the average energy deposited per interaction was achieved. Our results demonstrate the applicability of the electron dose-calculation method to photon dose calculation.     

Monte Carlo simulations of neutron spectral fluence, radiation weighting factor and ambient dose equivalent for a passively scattered proton therapy unit

Stray neutron exposures pose a potential risk for the development of secondary cancer in patients receiving proton therapy. However, the behavior of the ambient dose equivalent is not fully understood, including dependences on neutron spectral fluence, radiation weighting factor and proton treatment beam characteristics. The objective of this work, therefore, was to estimate neutron exposures resulting from the use of a passively scattered proton treatment unit. In particular, we studied the characteristics of the neutron spectral fluence, radiation weighting factor and ambient dose equivalent with Monte Carlo simulations. The neutron spectral fluence contained two pronounced peaks, one a low-energy peak with a mode around 1 MeV and one a high-energy peak that ranged from about 10 MeV up to the proton energy. The mean radiation weighting factors varied only slightly, from 8.8 to 10.3, with proton energy and location for a closed-aperture configuration. For unmodulated proton beams stopped in a closed aperture, the ambient dose equivalent from neutrons per therapeutic absorbed dose (H*(10)/D) calculated free-in-air ranged from about 0.3 mSv/Gy for a small scattered field of 100 MeV proton energy to 19 mSv/Gy for a large scattered field of 250 MeV proton energy, revealing strong dependences on proton energy and field size. Comparisons of in-air calculations with in-phantom calculations indicated that the in-air method yielded a conservative estimation of stray neutron radiation exposure for a prostate cancer patient.     

Neutron flux-density and secondary-particle energy spectra at the 184-inch synchrocyclotron medical facility

We have identified the sources of neutron production in the beam transport system of the 720-MeV helium beam used for radiation therapy at the 184-in synchrocyclotron of the Lawrence Berkeley Laboratory, and determined their magnitude. Measurements with activation detectors of differing energy response were used to unfold secondary particle spectra at various locations on the patient table. The effect of charged particles was estimated using a calculation of neutron-flux densities derived from published cross sections. The absorbed dose, as a function of distance from the beam axis, was calculated using the unfolded spectra and evaluated fluence-to-dose conversion factors. The values of absorbed dose obtained from the unfolding of experimental data agree with the values obtained from the calculated spectra within the estimated uncertainty of +/- 25%. These values are approximately 5 X 10(-3) rad on the beam axis and approximately 1 X 10(-3) rad at distances greater than 20 cm, perpendicular to the beam axis, per rad deposited by the incident alpha-particle beam in the plateau. Estimates of upper limits of dose to two critical organs, the lens of the eye and red bone marrow, are approximately 25 rad and approximately 5 rad, respectively, for a typical treatment plan.     

Radiobiological significance of beamline dependent proton energy distributions in a spread-out Bragg peak

Similar target doses can be achieved with different mixed radiation fields, i.e., particle energy distributions, produced by a practical proton beam and a range modulator. The dose delivered in particle therapy can be described as the integral of fluence times the total mass stopping power over the particle energy distributions. We employed Monte Carlo simulations to explore the influence on the relative biological effectiveness (RBE) of the energy and the energy spread of the proton beam incident on a range modulator system. Using different beams, the conditions of beam delivery were adjusted so that similar spread out Bragg peak (SOBP) doses were delivered to a simulated water phantom. We calculated the RBE for inactivation of three different cell lines using the track structure model. The RBE depends on the details of the dose deposition and the biological characteristics of the irradiated tissue. Our calculations show that, for differing beam conditions, the corresponding differences in the total mass stopping power distributions are reflected in differences in the RBE. However, these differences are remarkable only at the very distal edge of the SOBP, for low doses, and/or for large differences in beam setup.     

Parametrization and application of scatter kernels for modelling scanned proton beam collimator scatter dose

Collimators are routinely used in proton radiotherapy to laterally confine the field and improve the penumbra. Collimator scatter contributes up to 15% of the local dose and is therefore important to include in treatment planning dose calculation. We present a method for reconstruction of the collimator scatter phase space based on the parametrization of pre-calculated scatter kernels. Collimator scatter distributions, generated by the Monte Carlo (MC) package GEANT4.8.2, were scored differential in direction and energy. The distributions were then parametrized so as to enable a fast reconstruction by sampling. MC calculated dose distributions in water based on the parametrized phase space were compared to full MC simulations that included the collimator in the simulation geometry, as well as to experimental data. The experiments were performed at the scanned proton beam line at the The Svedberg Laboratory (TSL) in Uppsala, Sweden. Dose calculations using the parametrization of this work and the full MC for isolated typical cases of collimator scatter were compared by means of the gamma index. The result showed that in total 96.7% (99.3%) of the voxels fulfilled the gamma 2.0%/2.0 mm (3.0%/3.0 mm) criterion. The dose distribution for a collimated field was calculated based on the phase space created by the collimator scatter model incorporated into the generation of the phase space of a scanned proton beam. Comparing these dose distributions to full MC simulations, including particle transport in the MLC, yielded that in total for 18 different collimated fields, 99.1% of the voxels satisfied the gamma 1.0%/1.0 mm criterion and no voxel exceeded the gamma 2.6%/2.6 mm criterion. The dose contribution of collimator scatter along the central axis as predicted by the model showed good agreement with experimental data.