Diego M. González-Castaño

Automatic determination of primary electron beam parameters in Monte Carlo simulation

In order to obtain realistic and reliable Monte Carlo simulations of medical linac photon beams, an accurate determination of the parameters that define the primary electron beam that hits the target is a fundamental step. In this work we propose a new methodology to commission photon beams in Monte Carlo simulations that ensures the reproducibility of a wide range of clinically useful fields. For such purpose accelerated Monte Carlo simulations of 2 x 2, 10 x 10, and 20 x 20 cm2 fields at SSD = 100 cm are carried out for several combinations of the primary electron beam mean energy and radial FWHM. Then, by performing a simultaneous comparison with the correspondent measurements for these same fields, the best combination is selected. This methodology has been employed to determine the characteristics of the primary electron beams that best reproduce a Siemens PRIMUS and a Varian 2100 CD machine in the Monte Carlo simulations. Excellent agreements were obtained between simulations and measurements for a wide range of field sizes. Because precalculated profiles are stored in databases, the whole commissioning process can be fully automated, avoiding manual fine-tunings. These databases can also be used to characterize any accelerators of the same model from different sites.

Correction factors for A1SL ionization chamber dosimetry in TomoTherapy: Machine‐specific, plan‐class, and clinical fields

PURPOSE Recently, an international working group on nonstandard fields presented a new formalism for ionization chamber reference dosimetry of small and nonstandard fields [Alfonso et al., Med. Phys. 35, 5179-5186 (2008)] which has been adopted by AAPM TG-148. This work presents an experimental determination of the correction factors for reference dosimetry with an Exradin A1SL thimble ionization chamber in a TomoTherapy unit, focusing on: (i) machine-specific reference field, (ii) plan-class-specific reference field, and (iii) two clinical treatments. METHODS Ionization chamber measurements were performed in the TomoTherapy unit for intermediate (machine-specific and plan-class-specific) calibration fields, based on the reference conditions defined by AAPM TG-148, and two clinical treatments (lung and head-and-neck). Alanine reference dosimetry was employed to determine absorbed dose to water at the point of interest for the fields under investigation. The corresponding chamber correction factors were calculated from alanine to ionization chamber measurements ratios. RESULTS Two different methods of determining the beam quality correction factor k(Q,Q(0) ) for the A1SL ionization chamber in this TomoTherapy unit, where reference conditions for conventional beam quality determination cannot be met, result in consistent values. The observed values of overall correction factors obtained for intermediate and clinical fields are consistently around 0.98 with a typical expanded relative uncertainty of 2% (k = 2), which when considered make such correction factors compatible with unity. However, all of them are systematically lower than unity, which is shown to be significant when a hypothesis test assuming a t-student distribution is performed (p=1.8×10(-2)). Correction factors k(Q(clin),Q(pcsr) ) (f(clin),f(pcsr) ) and k(Q(clin),Q(msr) ) (f(clin),f(msr) ), which are needed for the computation of field factors for relative dosimetry of clinical beams, have been found to be very close to unity for two clinical treatments. CONCLUSIONS The results indicate that the helical field deliveries in this study (including two clinical fields) do not introduce changes on the ion chamber correction factors for dosimetry. For those two specific clinical cases, ratios of chamber readings accurately represent field output factors. The values observed here for intermediate calibration fields are in agreement with previously published data based on alanine dosimetry but differ from values recently reported obtained via radiochromic dosimetry.

The determination of beam quality correction factors: Monte Carlo simulations and measurements

Modern dosimetry protocols are based on the use of ionization chambers provided with a calibration factor in terms of absorbed dose to water. The basic formula to determine the absorbed dose at a users beam contains the well-known beam quality correction factor that is required whenever the quality of radiation used at calibration differs from that of the users radiation. The dosimetry protocols describe the whole ionization chamber calibration procedure and include tabulated beam quality correction factors which refer to 60Co gamma radiation used as calibration quality. They have been calculated for a series of ionization chambers and radiation qualities based on formulae, which are also described in the protocols. In the case of high-energy photon beams, the relative standard uncertainty of the beam quality correction factor is estimated to amount to 1%. In the present work, two alternative methods to determine beam quality correction factors are prescribed-Monte Carlo simulation using the EGSnrc system and an experimental method based on a comparison with a reference chamber. Both Monte Carlo calculations and ratio measurements were carried out for nine chambers at several radiation beams. Four chamber types are not included in the current dosimetry protocols. Beam quality corrections for the reference chamber at two beam qualities were also measured using a calorimeter at a PTB Primary Standards Dosimetry Laboratory. Good agreement between the Monte Carlo calculated (1% uncertainty) and measured (0.5% uncertainty) beam quality correction factors was obtained. Based on these results we propose that beam quality correction factors can be generated both by measurements and by the Monte Carlo simulations with an uncertainty at least comparable to that given in current dosimetry protocols.

Evaluation of chamber response function influence on IMRT verification using 2D commercial detector arrays

This work is devoted to studying the influence of chamber response functions on the standard IMRT verification for the different detector technologies available on commercial devices. We have tested three of the most used 2D detector arrays for radiotherapy dosimetry verification, based on air-ionization chambers and diode detectors. The response function has been carefully simulated using the Monte Carlo method and measured through slit and pinhole collimators. Although the response function of air-ionization detectors is considerably different with respect to that of standard diodes, the impact on a verification based in the gamma function with tolerances 3 mm and 3% is quite limited. The results show that the standard air-ionization detector arrays perform in a similar way whenever the tolerances for the gamma function are not lowered below 1.5 mm and 1.5%. Additionally, the sensitivity of these devices to fluence perturbations was measured by intentionally modifying some leaf positions in the multileaf collimator. The wider response function of air-ionization chamber arrays made them slightly more sensitive to random fluence perturbations, although silicon diode arrays are more accurate to describe the dose distribution in a point by point basis.

The change of response of ionization chambers in the penumbra and transmission regions: impact for IMRT verification

Significant deviations from the expected dose have been reported in the absolute dosimetry validation of an intensity modulated radiation therapy treatment when individual segments are analyzed. However, when full treatment is considered and all segment doses are added together, these discrepancies fade out, leading to overall dose deviations below a 5% action level. This contradictory behavior may be caused by a partial compensation between detector over-responding and under-responding for measurement conditions far from radiation equilibrium. We consider three treatment verification scenarios that may lead to ionization chamber miss-responding, namely: narrow beam irradiation, field penumbra location and multi-leaf collimator transmission contribution. In this work we have analyzed the response of three different ionization chambers with different active volume under these conditions by means of Monte Carlo (MC) simulation methods. Correction factors needed to convert the detector readout into actual dose to water were calculated by inserting the specific detector geometry (carefully modeled) into the simulations. This procedure required extensive use of parallel computing resources in order to achieve the desired level of uncertainty in the final results. The analysis of the simulations shows the relative contribution of each of the three previously mentioned miss-responding scenarios. Additionally, we provide some evidence on dose deviation compensation in multi-segment radiotherapy treatment verification.

A liquid-filled ionization chamber for high precision relative dosimetry

Radiosurgery and intensity modulated radiation therapy (IMRT) treatments are based on the delivery of narrow and/or irregularly shaped megavoltage photon beams. This kind of beams present both lack of charged particle equilibrium and steep dose gradients. Quality assurance (QA) measurements involved in these techniques must therefore be carried out with a dosimeter featuring high small volume. In order to obtain a good signal to noise ratio, a relatively dense material is needed as active medium. Non-polar organic liquids were proposed as active mediums with both good tissue equivalence and showing high signal to noise ratio. In this work, a liquid-filled ionization chamber is presented. Some results acquired with this detector in relative dosimetry are studied and compared with results obtained with unshielded diode. Medium-term stability measurements were also carried out and its results are shown. The liquid-filled ionization chamber presented here shows its ability to perform profile measurements and penumbrae determination with excellent accuracy. The chamber features a proper signal stability over the period studied.

A convolution model for obtaining the response of an ionization chamber in static non standard fields.

PURPOSE This work contains an alternative methodology for obtaining correction factors for ionization chamber (IC) dosimetry of small fields and composite fields such as IMRT. The method is based on the convolution/superposition (C/S) of an IC response function (RF) with the dose distribution in a certain plane which includes chamber position. This method is an alternative to the full Monte Carlo (MC) approach that has been used previously by many authors for the same objective. METHODS The readout of an IC at a point inside a phantom irradiated by a certain beam can be obtained as the convolution of the dose spatial distribution caused by the beam and the IC two-dimensional RF. The proposed methodology has been applied successfully to predict the response of a PTW 30013 IC when measuring different nonreference fields, namely: output factors of 6 MV small fields, beam profiles of cobalt 60 narrow fields and 6 MV radiosurgery segments. The two-dimensional RF of a PTW 30013 IC was obtained by MC simulation of the absorbed dose to cavity air when the IC was scanned by a 0.6 × 0.6 mm(2) cross section parallel pencil beam at low depth in a water phantom. For each of the cases studied, the results of the IC direct measurement were compared with the corresponding obtained by the C/S method. RESULTS For all of the cases studied, the agreement between the IC direct measurement and the IC calculated response was excellent (better than 1.5%). CONCLUSIONS This method could be implemented in TPS in order to calculate dosimetry correction factors when an experimental IMRT treatment verification with in-phantom ionization chamber is performed. The miss-response of the IC due to the nonreference conditions could be quickly corrected by this method rather than employing MC derived correction factors. This method can be considered as an alternative to the plan-class associated correction factors proposed recently as part of an IAEA work group on nonstandard field dosimetry.

A portable device for small animal SPECT imaging in clinical gamma-cameras

Molecular imaging is reshaping clinical practice in the last decades, providing practitioners with non-invasive ways to obtain functional in-vivo information on a diversity of relevant biological processes. The use of molecular imaging techniques in preclinical research is equally beneficial, but spreads more slowly due to the difficulties to justify a costly investment dedicated only to animal scanning. An alternative for lowering the costs is to repurpose parts of old clinical scanners to build new preclinical ones. Following this trend, we have designed, built, and characterized the performance of a portable system that can be attached to a clinical gamma-camera to make a preclinical single photon emission computed tomography scanner. Our system offers an image quality comparable to commercial systems at a fraction of their cost, and can be used with any existing gamma-camera with just an adaptation of the reconstruction software.

eIMRT: a web platform for the verification and optimization of radiation treatment plans

The eIMRT platform is a remote distributed computing tool that provides users with Internet access to three different services: Monte Carlo optimization of treatment plans, CRT & IMRT treatment optimization, and a database of relevant radiation treatments/clinical cases. These services are accessible through a user‐friendly and platform independent web page. Its flexible and scalable design focuses on providing the final users with services rather than a collection of software pieces. All input and output data (CT, contours, treatment plans and dose distributions) are handled using the DICOM format. The design, implementation, and support of the verification and optimization algorithms are hidden to the user. This allows a unified, robust handling of the software and hardware that enables these computation‐intensive services. The eIMRT platform is currently hosted by the Galician Supercomputing Center (CESGA) and may be accessible upon request (there is a demo version at http://eimrt.cesga.es:8080/eIMRT2/demo; request access in http://eimrt.cesga.es/signup.html). This paper describes all aspects of the eIMRT algorithms in depth its user interface, and its services. Due to the flexible design of the platform, it has numerous applications including the intercenter comparison of treatment planning, the quality assurance of radiation treatments, the design and implementation of new approaches to certain types of treatments, and the sharing of information on radiation treatment techniques. In addition, the web platform and software tools developed for treatment verification and optimization have a modular design that allows the user to extend them with new algorithms. This software is not a commercial product. It is the result of the collaborative effort of different public research institutions and is planned to be distributed as an open source project. In this way, it will be available to any user; new releases will be generated with the new implemented codes or upgrades. PACS number: 87.55.kh

A subthreshold, low-power, RHBD reference circuit, for earth observation and communication satellites

A low-power, wide temperature range, radiation tolerant CMOS voltage reference is presented. The proposed reference circuit exhibits a voltage deviation of 0.8mV for 3-MeV protons total ionization dose of 2Mrad and a voltage deviation of 3.8mV for 10-keV X-rays total ionization dose of 4Mrad while being biased at the nominal supply voltage of 0.75V during X-ray irradiation. In addition, the circuit consumes only 4μW and exhibits a measured Temperature Drift of 15ppm/°C for a temperature range of 190°C (-60°C to 130°C) at the supply voltage of 0.75V. It utilizes only CMOS transistors, operating in the subthreshold regime, and poly-silicon resistors without using any diodes or external components such as compensating capacitors. The circuit is radiation hardened by design (RHBD), it was fabricated using TowerJazz Semiconductors 0.18μm standard CMOS technology and occupies a silicon area of 0.039mm2. The proposed voltage reference is suitable for high-precision and low-power space applications.