Low-cost flexible thin-film detector for medical dosimetry applications.

The purpose of this study is to characterize dosimetric properties of thin film photovoltaic sensors as a platform for development of prototype dose verification equipment in radiotherapy. Towards this goal, flexible thin-film sensors of dose with embedded data acquisition electronics and wireless data transmission are prototyped and tested in kV and MV photon beams. Fundamental dosimetric properties are determined in view of a specific application to dose verification in multiple planes or curved surfaces inside a phantom. Uniqueness of the new thin-film sensors consists in their mechanical properties, low-power operation, and low-cost. They are thinner and more flexible than dosimetric films. In principle, each thin-film sensor can be fabricated in any size (mm² - cm² areas) and shape. Individual sensors can be put together in an array of sensors spreading over large areas and yet being light. Photovoltaic mode of charge collection (of electrons and holes) does not require external electric field applied to the sensor, and this implies simplicity of data acquisition electronics and low power operation. The prototype device used for testing consists of several thin film dose sensors, each of about 1.5 cm × 5 cm area, connected to simple readout electronics. Sensitivity of the sensors is determined per unit area and compared to EPID sensitivity, as well as other standard photodiodes. Each sensor independently measures dose and is based on commercially available flexible thin-film aSi photodiodes. Readout electronics consists of an ultra low-power microcontroller, radio frequency transmitter, and a low-noise amplification circuit implemented on a flexible printed circuit board. Detector output is digitized and transmitted wirelessly to an external host computer where it is integrated and processed. A megavoltage medical linear accelerator (Varian Tx) equipped with kilovoltage online imaging system and a Cobalt source are used to irradiate different thin-film detector sensors in a Solid Water phantom under various irradiation conditions. Different factors are considered in characterization of the device attributes: energies (80 kVp, 130 kVp, 6 MV, 15 MV), dose rates (different ms × mA, 100-600 MU/min), total doses (0.1 cGy-500 cGy), depths (0.5 cm-20 cm), irradiation angles with respect to the detector surface (0°-180°), and IMRT tests (closed MLC, sweeping gap). The detector response to MV radiation is both linear with total dose (~1-400 cGy) and independent of dose rate (100-600 Mu/min). The sensitivity per unit area of thin-film sensors is lower than for aSi flat-panel detectors, but sufficient to acquire stable and accurate signals during irradiations. The proposed thin-film photodiode system has properties which make it promising for clinical dosimetry. Due to the mechanical flexibility of each sensor and readout electronics, low-cost, and wireless data acquisition, it could be considered for quality assurance (e.g., IMRT, mechanical linac QA), as well as real-time dose monitoring in challenging setup configurations, including large area and 3D detection (multiple planes or curved surfaces).

An efficient inverse radiotherapy planning method for VMAT using quadratic programming optimization.

PURPOSE: The purpose of this study is to investigate the feasibility of an inverse planning optimization approach for the Volumetric Modulated Arc Therapy (VMAT) based on quadratic programming and the projection method. The performance of this method is evaluated against a reference commercial planning system (eclipse(TM) for rapidarc(TM)) for clinically relevant cases. METHODS: The inverse problem is posed in terms of a linear combination of basis functions representing arclet dose contributions and their respective linear coefficients as degrees of freedom. MLC motion is decomposed into basic motion patterns in an intuitive manner leading to a system of equations with a relatively small number of equations and unknowns. These equations are solved using quadratic programming under certain limiting physical conditions for the solution, such as the avoidance of negative dose during optimization and Monitor Unit reduction. The modeling by the projection method assures a unique treatment plan with beneficial properties, such as the explicit relation between organ weightings and the final dose distribution. Clinical cases studied include prostate and spine treatments. The optimized plans are evaluated by comparing isodose lines, DVH profiles for target and normal organs, and Monitor Units to those obtained by the clinical treatment planning system eclipse(TM). RESULTS: The resulting dose distributions for a prostate (with rectum and bladder as organs at risk), and for a spine case (with kidneys, liver, lung and heart as organs at risk) are presented. Overall, the results indicate that similar plan qualities for quadratic programming (QP) and rapidarc(TM) could be achieved at significantly more efficient computational and planning effort using QP. Additionally, results for the quasimodo phantom [Bohsung et al., "IMRT treatment planning: A comparative inter-system and inter-centre planning exercise of the estro quasimodo group," Radiother. Oncol. 76(3), 354-361 (2005)] are presented as an example for an extreme concave case. CONCLUSION: Quadratic programming is an alternative approach for inverse planning which generates clinically satisfying plans in comparison to the clinical system and constitutes an efficient optimization process characterized by uniqueness and reproducibility of the solution.

Localization of deformable tumors from short-arc projections using Bayesian estimation.

PURPOSE: The authors present a stochastic framework for radiotherapy patient positioning directly utilizing radiographic projections. This framework is developed to be robust against anatomical nonrigid deformations and to cope with challenging imaging scenarios, involving only a few cone beam CT projections from short arcs. METHODS: Specifically, a Bayesian estimator (BE) is explicitly derived for the given scanning geometry. This estimator is compared to reference methods such as chamfer matching (CM) and the minimization of the median absolute error adapted as tools of robust image processing and statistics. In order to show the performance of the stochastic short-arc patient positioning method, a CIRS IMRT thorax phantom study is presented with movable markers and the utilization of an Elekta Synergy(®) XVI system. Furthermore, a clinical prostate CBCT scan of a Varian(®) On-Board Imager(®) system is utilized to investigate the robustness of the method for large variations of image quality (anterior-posterior vs lateral views). RESULTS: The results show that the BE shifts reduce the initial setup error of up to 3 cm down to 3 mm at maximum for an imaging arc as short as 10° while CM achieves residual errors of 7 mm at maximum only for arcs longer than 40°. Furthermore, the BE can compensate robustly for low image qualities using several low quality projections simultaneously. CONCLUSIONS: In conclusion, an estimation method for marker-based patient positioning for short imaging arcs is presented and shown to be robust and accurate for deformable anatomies.

Strip detector array (SDA) for beam monitoring in radiotherapy: reconstruction of MLC parameters from multiple projections of flux.

Objective.In this paper we propose and investigate a new detector with multiple strip detector arrays (SDA) for monitoring MLC shaped x-ray beams for radiotherapy treatment.Approach.Each SDA measures 1D dose profiles equivalent to dose projections. The goal of such a detector is to determine individual MLC leaf positions as well as the Monitor Units (MU) per MLC segment during radiotherapy. In the present work we investigate an optimal SDA detector configuration and reconstruction algorithm. We determine the accuracy of SDA for different treatment sites (spine, pelvis, retroperitoneum, prostate, brain SRT, SRS, lung and head and neck). We perform a simulation study accounting for different type of MLC leaf positional errors: random MLC leaf, systematic for the whole leaf bank and systematic for an individual leaf. In a similar fashion, we also account for errors in Monitor Units per segment.Main results.We demonstrate that for a broad range of IMRT treatment plans a robust reconstruction of errors is achievable with only 3 projections (3 sets of SDA oriented at at 0°, 45° and 135°). The SDA is capable of capturing both systematic errors in leaf banks and individual leaves as well as random errors sufficient for practical clinical purposes.Significance.These features of the SDA detector makes it suitable for real-time Quality Control of MLC collimated linac output.

Stochastic formulation of patient positioning using linac-mounted cone beam imaging with prior knowledge.

PURPOSE: In this work, a novel stochastic framework for patient positioning based on linac-mounted CB projections is introduced. Based on this formulation, the most probable shifts and rotations of the patient are estimated, incorporating interfractional deformations of patient anatomy and other uncertainties associated with patient setup. METHODS: The target position is assumed to be defined by and is stochastically determined from positions of various features such as anatomical landmarks or markers in CB projections, i.e., radiographs acquired with a CB-CT system. The patient positioning problem of finding the target location from CB projections is posed as an inverse problem with prior knowledge and is solved using a Bayesian maximum a posteriori (MAP) approach. The prior knowledge is three-fold and includes the accuracy of an initial patient setup (such as in-room laser and skin marks), the plasticity of the body (relative shifts between target and features), and the feature detection error in CB projections (which may vary depending on specific detection algorithm and feature type). For this purpose, MAP estimators are derived and a procedure of using them in clinical practice is outlined. Furthermore, a rule of thumb is theoretically derived, relating basic parameters of the prior knowledge (initial setup accuracy, plasticity of the body, and number of features) and the parameters of CB data acquisition (number of projections and accuracy of feature detection) to the expected estimation accuracy. RESULTS: MAP estimation can be applied to arbitrary features and detection algorithms. However, to experimentally demonstrate its applicability and to perform the validation of the algorithm, a water-equivalent, deformable phantom with features represented by six 1 mm chrome balls were utilized. These features were detected in the cone beam projections (XVI, Elekta Synergy) by a local threshold method for demonstration purposes only. The accuracy of estimation (strongly varying for different plasticity parameters of the body) agreed with the rule of thumb formula. Moreover, based on this rule of thumb formula, about 20 projections for 6 detectable features seem to be sufficient for a target estimation accuracy of 0.2 cm, even for relatively large feature detection errors with standard deviation of 0.5 cm and spatial displacements of the features with standard deviation of 0.5 cm. CONCLUSIONS: The authors have introduced a general MAP-based patient setup algorithm accounting for different sources of uncertainties, which are utilized as the prior knowledge in a transparent way. This new framework can be further utilized for different clinical sites, as well as theoretical developments in the field of patient positioning for radiotherapy.

A modification of flattening filter free linac for IMRT.

PURPOSE: This study investigates the benefits of a modified flattening filter free (FFF) linac over the standard (STD) linac equipped with the flattening filter. Energy and angular spread of the electron beam of the FFF linac were modified. Modification of FFF beam parameters is explored to maximize the monitor unit efficiency and to minimize the head scatter in IMRT delivery for large target volumes or targets lying away from the central axis. METHODS: The EGSnrc code is used to model FFF and STD linacs and study basic beam properties for both linac types in various beam configurations. Increasing energy of FFF linac results in similar beam attenuation properties and maximized dose rate compared to STD linac. Matching beam attenuation properties allows a more direct exploration of beam flatness of FFF linac in regard to IMRT delivery, especially away from the central axis where the effective dose rate is considerably smaller than the one at the central axis. Flatness of open beam dose profile of FFF linac is improved by increasing the angular spread of the electron beam. The resulting dose rate within the treatment field and outside of the field (peripheral dose) are characterized and compared to the unmodified FFF and STD linacs, RESULTS: In order to match beam penetration properties, the energy of FFF is adjusted from 6.5 to 8.0 MeV for small to medium field sizes and from 6.5 to 8.5 MeV for larger ones. Dose rate of FFF vs STD linac increased by a factor of 1.9 (6.5 MeV) and 3.4-4.1 (8.0-8.5 MeV). Adjusting the mean angular spread of the electron beam from 0 degrees to 5 degrees-10 degrees resulted in complete flattening of photon beam for field sizes between 10 x 10 cm2 and 15 x 15 cm2 and partial flattening for field sizes from 15 x 15 cm2 to 30 x 30 cm2. Values of angular spread > or =14 degrees are not recommended as they exceed the opening of the primary collimator, affecting the area at the edges of the field. FFF fields of sizes smaller than 6 x 6 cm2 are already flat and beam flattening is not necessary. Overall, the angular spread of 5 degrees-10 degrees is sufficient and can satisfactorily flatten open beam dose profiles even for larger field sizes. Increasing the electron beam angular spread amounts to a slight decrease of dose rate of FFF linac. However, for angular spread, 5 degrees-10 degrees dose rate factor of FFF vs STD is still about 1.6-2.6, depending on the field size (and the adjusted energy). Similarly, in case of peripheral dose, a moderate increase in dose can be observed for angular spread of 5 degrees-10 degrees and for field sizes 10 x 10 cm2 to 30 x 30 cm2. Lastly, beam flatness of not modified FFF linac can be conveniently described by an analytical function representing a ratio of STD vs FFF doses: 1 + b|r|(n). CONCLUSIONS: A modified FFF beamline with increased energy and electron beam angular spread results in satisfactory flattened beam and high dose rate within the field. Peripheral dose remaining at similar (or smaller) level than that of STD linac for the same delivered dose within the treatment field.

Gold nanoparticle detection and quantification in therapeutic MV beams via pair production.

PURPOSE: We propose a new detection method of gold nanoparticles (AuNP) in therapeutic megavoltage (MV) x-ray beams by means of coincidence counting of annihilation photons following pair production in gold. METHODS: The proposed MV x-ray induced positron emission (MVIPE) imaging technique is studied by radiation transport computations using MCNP6 (3D) and CEPXS/ONEDANT (1D) codes for two water phantoms: a 35 cm slab and a similarly sized cylinder, both having a 5 cm AuNP filled region in the center. MVIPE is compared to the standard x-ray fluorescence computed tomography (XFCT). MVIPE adopts MV x-ray sources (Co-60, 2 MV, 6 MV, 6 MV with closed MLC and 15 MV) and relies on the detection of 511 keV photon-pairs. XFCT uses kilovoltage sources (100 kVp, 120 kVp and 150 kVp) and imaging is characterized by analysis of k α1,2 Au characteristic lines. Three levels of AuNP concentration were studied: 0.1%, 1% and 10% by weight. RESULTS: Annihilation photons in the MVIPE technique originate both in the AuNP and in water along the x-ray beam path with significantly larger production in the AuNP-loaded region. MVIPE signal from AuNP is linearly increasing with AuNP concentration up to 10%wt, while XFCT signal reaches saturation due to self-absorption within AuNP. The production of annihilation photons is proportional to the MV source energy. MVIPE technique using a 15 MV pencil beam and 10 wt% AuNP detects about 4.5 × 103 511 keV-photons cm-2 at 90° w/r to the incident beam per 109 source photons cm-2; 500 of these come from AuNP. In contrast, the XFCT technique using 150 kVp detects only about 100 k α1-photons cm-2 per 109 source photons cm-2. CONCLUSIONS: In MVIPE, the number of annihilation photons produced for different MV-beam energies and AuNP concentrations is significantly greater than the k α1 photons generated in XFCT. Coincidence counting in MVIPE allows to avoid collimation, which is a major limiting factor in XFCT. MVIPE challenges include the filtering of Compton scatter and annihilation photons originating in water.

Balancing dose and image registration accuracy for cone beam tomosynthesis (CBTS) for breast patient setup.

PURPOSE: To balance dose reduction and image registration accuracy in breast setup imaging. In particular, the authors demonstrate the relationship between scan angle and dose delivery for cone beam tomosynthesis (CBTS) when employed for setup verification of breast cancer patients with surgical clips. METHODS: The dose measurements were performed in a female torso phantom for varying scan angles of CBTS. Setup accuracy was measured using three registration methods: Clip centroid localization accuracy and the accuracy of two semiautomatic registration algorithms. The dose to the organs outside of the ipsilateral breast and registration accuracy information were compared to determine the optimal scan angle for CBTS for breast patient setup verification. Isocenter positions at the center of the patient and at the breast-chest wall interface were considered. RESULTS: Image registration accuracy was within 1 mm for the CBTS scan angles theta above 20 degrees for some scenarios and as large as 80 degrees for the worst case, depending on the imaged breast and registration algorithm. Registration accuracy was highest based on clip centroid localization. For left and right breast imaging with the isocenter at the chest wall, the dose to the contralateral side of the patient was very low (<0.5 cGy) for all scan angles considered. For central isocenter location, the optimal scan angles were 30 degrees - 50 degrees for the left breast imaging and 40 degrees - 50 degrees for the right breast imaging, with the difference due to the geometric asymmetry of the current clinical imaging system. CONCLUSIONS: The optimal scan angles for CBTS imaging were found to be between 10 degrees and 50 degrees, depending on the isocenter location and ipsilateral breast. Use of the isocenter at the breast-chest wall locations always resulted in greater accuracy of image registration (<1 mm) at smaller angles (10 degrees - 20 degrees) and at lower doses (<0.1 cGy) to the contralateral organs. For chest wall isocenters, doses delivered to organs outside of the target breast were much smaller than the scattered and leakage doses of the treatment beams. The complete volumetric information of all clips in the region of interest, combined with the small dose to the contralateral organs and the small scan angle, could result in an advantage for small angle CBTS with off center isocenters over simple orthogonal pairs.

An alternative VMAT with prior knowledge about the type of leaf motion utilizing projection method for concave targets.

The authors present an alternative approach to inverse planning optimization and apply it to volumetric modulated are therapy (VMAT) in one rotation with a prior knowledge about the type of leaf motions. The optimization is based on the projection theorem in inner product spaces. MLC motion is directly considered in the optimization, thus avoiding leaf segmentation characteristic of IMRT optimization. In this work they realize the method for concave irregular targets encompassing an organ at risk leading to a repetitive MLC motion pattern. Applying the projection theorem leads to a noniterative optimization method and reduces to solving few systems of linear equations with small numbers of dimensions. The solution of the inverse problem is unique, and false minima are naturally excluded. They divided the full rotation into about 50 short arc segments and for each segment decomposed dose into separate contributions related to stages of MLC motion. This results generally in an inverse problem with just four free parameters per arc segment. Practically three degrees of freedom will be used for the purpose of a constant angular speed of the gantry. Therefore the total number of degrees of freedom for a 3D problem is about 3 x 50 x number of collimator leaf pairs for irradiating the whole target volume in one rotation. Two 2D and one 3D concave target volumes are applied for a slice by slice optimization. A 6 MV photon beam model is used, including realistic scattering and attenuation, and a maximal leaf velocity of 3 cm/s is regarded. The resulting dose distributions cover the PTVs very well and have maxima at about 108% of dose in the PTVs. The OAR is spared very strong in all cases. As a result of optimization, the MLC apertures are repetitively opening and closing and can be interpreted in an intuitive way. Applying the projection method for this knowledge-based VMAT delivery scheme for concave target volumes is an alternative technique for dose optimization. There are several properties, such as uniqueness of MLC motions and their continuous dependence on geometry and prescribed dose, that make this approach interesting to inverse planning. This method is still in an investigational stage, but promising results are presented. In future work it will be extended directly (without conceptual changes) in several directions to be more clinically applicable.

Comment on: Polarity effects and apparent ion recombination in microionization chambers Med. Phys. 43(5), 2141-2152 (2016).

In a recent paper of Miller et al. (Med. Phys. 43(5), 2141-2152 (2016), the voltage dependence polarity effect of microionization chambers was studied. It was identified that polarity effect arises from the difference in electric potential between the collecting electrode and the guard electrode, which in turn deforms the electric field and affects the charge collection. Nevertheless, the cause of such potential difference has not been identified. In this letter, we explain that the electric potential perturbation arises from the work function difference of the disparate materials electrodes (collecting vs guard electrodes in the particular case).

A Monte Carlo study of I-125 prostate brachytherapy with gold nanoparticles: dose enhancement with simultaneous rectal dose sparing via radiation shielding.

We investigate via Monte Carlo simulations a new 125I brachytherapy treatment technique for high-risk prostate cancer patients via injection of Au nanoparticle (AuNP) directly into the prostate. The purpose of using the nanoparticles is to increase the therapeutic index via two synergistic effects: enhanced energy deposition within the prostate and simultaneous shielding of organs at risk from radiation escaping from the prostate. Both uniform and non-uniform concentrations of AuNP are studied. The latter are modeled considering the possibility of AuNP diffusion after the injection using brachy needles. We study two extreme cases of coaxial AuNP concentrations: centered on brachy needles and centered half-way between them. Assuming uniform distribution of 30 mg g-1 of AuNP within the prostate, we obtain a dose enhancement larger than a factor of 2 to the prostate. Non-uniform concentration of AuNP ranging from 10 mg g-1 and 66 mg g-1 were studied. The higher the concentration in a given region of the prostate the greater is the enhancement therein. We obtain the highest dose enhancement when the brachytherapy needles are coincident with AuNP injection needles but, at the same time, the regions in the tail are colder (average dose ratio of 0.7). The best enhancement uniformity is obtained with the seeds in the tail of the AuNP distribution. In both uniform and non-uniform cases the urethra and rectum receive less than 1/3 dose compared to an analog treatment without AuNP. Remarkably, employing AuNP not only significantly increases dose to the target but also decreases dose to the neighboring rectum and even urethra, which is embedded within the prostate. These are mutually interdependent effects as more enhancement leads to more shielding and vice-versa. Caution must be paid since cold spot or hot spots may be created if the AuNP concentration versus seed position is not properly distributed respect to the seed locations.

Kilovoltage radiosurgery with gold nanoparticles for neovascular age-related macular degeneration (AMD): a Monte Carlo evaluation.

This work uses Monte Carlo radiation transport simulation to assess the potential benefits of gold nanoparticles (AuNP) in the treatment of neovascular age-related macular degeneration with stereotactic radiosurgery. Clinically, a 100 kVp x-ray beam of 4 mm diameter is aimed at the macula to deliver an ablative dose in a single fraction. In the transport model, AuNP accumulated at the bottom of the macula are targeted with a source representative of the clinical beam in order to provide enhanced dose to the diseased macular endothelial cells. It is observed that, because of the AuNP, the dose to the endothelial cells can be significantly enhanced, allowing for greater sparing of optic nerve, retina and other neighboring healthy tissue. For 20 nm diameter AuNP concentration of 32 mg g(-1), which has been shown to be achievable in vivo, a dose enhancement ratio (DER) of 1.97 was found to be possible, which could potentially be increased through appropriate optimization of beam quality and/or AuNP targeting. A significant enhancement in dose is seen in the vicinity of the AuNP layer within 30 µm, peaked at the AuNP-tissue interface. Different angular tilting of the 4 mm beam results in a similar enhancement. The DER inside and in the penumbra of the 4 mm irradiation-field are almost the same while the actual delivered dose is more than one order of magnitude lower outside the field leading to normal tissue sparing. The prescribed dose to macular endothelial cells can be delivered using almost half of the radiation allowing reduction of dose to the neighboring organs such as retina/optic nerve by 49% when compared to a treatment without AuNP.

Comprehensive quality assurance phantom for the small animal radiation research platform (SARRP).

PURPOSE: To develop and test the suitability and performance of a comprehensive quality assurance (QA) phantom for the Small Animal Radiation Research Platform (SARRP). METHODS AND MATERIALS: A QA phantom was developed for carrying out daily, monthly and annual QA tasks including: imaging, dosimetry and treatment planning system (TPS) performance evaluation of the SARRP. The QA phantom consists of 15 (60 × 60 × 5 mm(3)) kV-energy tissue equivalent solid water slabs. The phantom can incorporate optically stimulated luminescence dosimeters (OSLD), Mosfet or film. One slab, with inserts and another slab with hole patterns are particularly designed for image QA. RESULTS: Output constancy measurement results showed daily variations within 3%. Using the Mosfet in phantom as target, results showed that the difference between TPS calculations and measurements was within 5%. Annual QA results for the Percentage depth dose (PDD) curves, lateral beam profiles, beam flatness and beam profile symmetry were found consistent with results obtained at commissioning. PDD curves obtained using film and OSLDs showed good agreement. Image QA was performed monthly, with image-quality parameters assessed in terms of CBCT image geometric accuracy, CT number accuracy, image spatial resolution, noise and image uniformity. CONCLUSIONS: The results show that the developed QA phantom can be employed as a tool for comprehensive performance evaluation of the SARRP. The study provides a useful reference for development of a comprehensive quality assurance program for the SARRP and other similar small animal irradiators, with proposed tolerances and frequency of required tests.

Method of identifying dynamic multileaf collimator irradiation that is highly sensitive to a systematic MLC calibration error.

In intensity modulated radiotherapy (IMRT), radiation is delivered in a multiple of multileaf collimator (MLC) subfields. A subfield with a small leaf-to-leaf opening is highly sensitive to a leaf-positional error. We introduce a method of identifying and rejecting IMRT plans that are highly sensitive to a systematic MLC gap error (sensitivity to possible random leaf-positional errors is not addressed here). There are two sources of a systematic MLC gap error: centerline mechanical offset (CMO) and, in the case of a rounded end MLC, radiation field offset (RFO). In IMRT planning system, using an incorrect value of RFO introduces a systematic error ARFO that results in all leaf-to-leaf gaps that are either too large or too small by (2*DeltaRFO), whereas assuming that CMO is zero introduces systematic error DeltaCMO that results in all gaps that are too large by DeltaCMO=CMO. We introduce a concept of the average leaf pair Opening (ALPO) that can be calculated from a dynamic MLC delivery file. We derive an analytic formula for a fractional average fluence error resulting from a systematic gap error of Deltax and show that it is inversely proportional to ALPO; explicitly it is equal to Deltax/(ALPO+ 2 * RFO+ epsilon), in which epsilon is generally of the order of 1 mm and Deltax =2 * Delta RFO + CMO. This analytic relationship is verified with independent numerical calculations.

A method of calculating a lung clinical target volume DVH for IMRT with intrafractional motion.

The motion of lung tumors from respiration has been reported in the literature to be as large as 1-2 cm. This motion requires an additional margin between the Clinical Target Volume (CTV) and the Planning Target Volume (PTV). In Intensity Modulated Radiotherapy (IMRT), while such a margin is necessary, the margin may not be sufficient to avoid unintended high and low dose regions to the interior on moving CTV. Gated treatment has been proposed to improve normal tissues sparing as well as to ensure accurate dose coverage of the tumor volume. The following questions have not been addressed in the literature: (a) what is the dose error to a target volume without a gated IMRT treatment? (b) What is an acceptable gating window for such a treatment. In this study, we address these questions by proposing a novel technique for calculating the three-dimensional (3-D) dose error that would result if a lung IMRT plan were delivered without a gated linac beam. The method is also generalized for gated treatment with an arbitrary triggering window. IMRT plans for three patients with lung tumors were studied. The treatment plans were generated with HELIOS for delivery with 6 MV on a CL2100 Varian linear accelerator with a 26 pair MLC. A CTV to PTV margin of 1 cm was used. An IMRT planning system searches for an optimized fluence map phi(x,y) for each port, which is then converted into a dynamic MLC file (DMLC). The DMLC file contains information about MLC subfield shapes and the fractional Monitor Units (MUs) to be delivered for each subfield. With a lung tumor, a CTV that executes a quasiperiodic motion z(t) does not receive phi(x,y), but rather an Effective Incident Fluence EIF(x,y). We numerically evaluate the EIF(x,y) from a given DMLC file by a coordinate transformation to the Target's Eye View (TEV). In the TEV coordinate system, the CTV itself is stationary, and the MLC is seen to execute a motion -z(t) that is superimposed on the DMLC motion. The resulting EIF(x,y) is input back into the dose calculation engine to estimate the 3-D dose to a moving CTV. In this study, we model respiratory motion as a sinusoidal function with an amplitude of 10 mm in the superior-inferior direction, a period of 5 s, and an initial phase of zero.

The measurement of proton stopping power using proton-cone-beam computed tomography.

A cone-beam computed tomography (CT) system utilizing a proton beam has been developed and tested. The cone beam is produced by scattering a 160 MeV proton beam with a modifier that results in a signal in the detector system, which decreases monotonically with depth in the medium. The detector system consists of a Gd2O2S:Tb intensifying screen viewed by a cooled CCD camera. The Feldkamp-Davis-Kress cone-beam reconstruction algorithm is applied to the projection data to obtain the CT voxel data representing proton stopping power. The system described is capable of reconstructing data over a 16 x 16 x 16 cm3 volume into 512 x 512 x 512 voxels. A spatial and contrast resolution phantom was scanned to determine the performance of the system. Spatial resolution is significantly degraded by multiple Coulomb scattering effects. Comparison of the reconstructed proton CT values with x-ray CT derived proton stopping powers shows that there may be some advantage to obtaining stopping powers directly with proton CT. The system described suggests a possible practical method of obtaining this measurement in vivo.

Beam quality and dose perturbation of 6 MV flattening-filter-free linac.

The aim of this study is twofold: (a) determination of the spectral differences for flattening-filter-free (FFF) versus standard (STD) linac under various clinical conditions, (b) based on an extensive list of clinically important beam configurations, identification of clinical scenarios that lead to higher macroscopic dose perturbations due to the presence of high-Z material. The focus is on dose enhancement due to contrast agents including high-Z elements such as gold or gadolinium. EGSnrc was used to simulate clinical beams under various irradiation conditions: open/IMRT/spit-IMRT fields, in/out-off-field areas, different depths and field sizes. Spectra were calculated and analyzed for about 80 beams and for a total of 480 regions. Quantitative differential effects in beam quality were characterized using energy-dependent and cumulative dose perturbation metrics. Analysis of the spectral database showed that even though the general trends for both linacs (FFF/STD) were the same, there were crucial differences. In general, the relative changes between different conditions were smaller for FFF spectra. This was because of the higher component of low-energy photons of the FFF linac, which already lead to higher dose enhancement than for the STD linac (photon energies were more "uniformly" distributed for FFF spectra and henceforth their perturbation resulted in lesser relative changes). For out-of-field FFF spectra and split-IMRT fields the strongest enhancement were observed (∼25 and ∼5 respectively). Different spectral scenarios lead to different dose enhancements, however, they scale with the higher effective-Z of the materials and were directly related to the lower range of the spectra (<200 keV).

Low-Z linac targets for low-MV gold nanoparticle radiation therapy.

PURPOSE: To investigate the potential of low-Z/low-MV (low-Z) linac targets for gold nanoparticle radiotherapy (GNPT) and to determine the microscopic dose enhancement ratio (DER) due to GNP for the alternative beamlines. In addition, to evaluate the degradation of dose enhancement arising from the increased attenuation of x rays and larger skin dose in water for the low-MV beams compared to the standard linac. METHODS: Monte Carlo simulations were used to compute dose and DER for various flattening-filter-free beams (2.5, 4, 6.5 MV). Target materials were beryllium, diamond, and tungsten-copper high-Z target. Target thicknesses were selected based on 20%, 60%, 70%, and 80% of the continuous slowing down approximation electron ranges for a given target material and energy. Evaluation of the microscopic DER was carried out for 100 nm GNP including the degradation factors due to beam attenuation. RESULTS: The greatest increase in DER compared to the standard 6.5 MV linac was for a 2.5 MV Be-target (factor of ∼ 2). Skin dose ranged from ∼ 10% (Be, 6.5 MV-80%) to ∼ 85% (Be, 2.5 MV-20%) depending on the target case. Attenuation of 2.5 MV beams at 22 cm was higher by ∼ 75% compared with the standard beam. Taking into account the attenuation at 22 cm depth, the effective dose enhancement was up to ∼ 60% above the DER of the high-Z target. For these cases the effective DER ranged between ∼ 1.6 and 6 compared with the standard linac. CONCLUSIONS: Low-Z (2.5 MV) GNPT is possible even after accounting for greater beam attenuation for deep-seated tumors (22 cm) and the increased skin dose. Further, it can lead to significant sparing of normal tissue while simultaneously escalating the dose in the tumor cells.

High throughput film dosimetry in homogeneous and heterogeneous media for a small animal irradiator.

PURPOSE: We have established a high-throughput Gafchromic film dosimetry protocol for narrow kilovoltage beams in homogeneous and heterogeneous media for small-animal radiotherapy applications. The kV beam characterization is based on extensive Gafchromic film dosimetry data acquired in homogeneous and heterogeneous media. An empirical model is used for parameterization of depth and off-axis dependence of measured data. METHODS: We have modified previously published methods of film dosimetry to suit the specific tasks of the study. Unlike film protocols used in previous studies, our protocol employs simultaneous multi-channel scanning and analysis of up to nine Gafchromic films per scan. A scanner and background correction were implemented to improve accuracy of the measurements. Measurements were taken in homogeneous and inhomogeneous phantoms at 220 kVp and a field size of 5 × 5 mm(2). The results were compared against Monte Carlo simulations. RESULTS: Dose differences caused by variations in background signal were effectively removed by the corrections applied. Measurements in homogeneous phantoms were used to empirically characterize beam data in homogeneous and heterogeneous media. Film measurements in inhomogeneous phantoms and their empirical parameterization differed by about 2%-3%. The model differed from MC by about 1% (water, lung) to 7% (bone). Good agreement was found for measured and modelled off-axis ratios. CONCLUSIONS: EBT2 films are a valuable tool for characterization of narrow kV beams, though care must be taken to eliminate disturbances caused by varying background signals. The usefulness of the empirical beam model in interpretation and parameterization of film data was demonstrated.

A fast analytic dose calculation method for arc treatments for kilovoltage small animal irradiators.

Arc treatments require calculation of dose for collections of discrete gantry angles. The sampling of angles must balance between short computation time of small angle sets and the better calculation reliability of large sets. In this paper, an analytical formula is presented that allows calculation of dose delivered during continuous rotation of the gantry. The formula holds valid for continuous short arcs of up to about 30° and is derived by integrating a dose formula over gantry angles within a small angle approximation. Doses for longer arcs may be obtained in terms of doses for shorter arcs. The formula is derived with an empirical beam model in water and extended to inhomogeneous media. It is validated with experimental data obtained by applying arc treatment using kV small animal irradiator to a phantom of solid water and lung-equivalent material. The results are a promising step towards efficient 3D dose calculation and inverse planning purposes. In principle, this method also applies to VMAT dose calculation and optimization but requires extensions.