A new method for output factor determination in MLC shaped narrow beams

A new method for the measurement of output factors of narrow beams is presented in this work. By combining a new large area parallel plane ionization chamber (PTW model T34070) with a relative film dosimetry the output factors of small square fields of a 6 MV beam shaped by a MLC were measured. Several detectors (three ionization chambers, two solid state detectors and film) and Monte Carlo simulation were also employed to validate this new methodology and also to determine those detectors more suitable for narrow beam output factor determination.The proposed method for output factor measurement has shown to be in a very good agreement with diamond, diode and Monte Carlo results while it is insensitive to position displacements. Several uncertainties associated to the process of narrow beam dosimetry have also been addressed.     

Towards MR-guided electron therapy: Measurement and simulation of clinical electron beams in magnetic fields

PurposeIn the current era of MRI-linac radiotherapy, dose optimization with arbitrary dose distributions is a reality. For the first time, we present new and targeted experiments and modeling to aid in evaluating the potential dose improvements offered with an electron beam mode during MRI-linac radiotherapy.MethodsSmall collimated (1 cm diameter and 1.5 × 1.5 cm² square) electron beams (6, 12 and 20 MeV) from a clinical linear accelerator (Varian Clinac 2100C) are incident perpendicular and parallel to the strong and localized magnetic fields (0–0.7 T) generated by a permanent magnet device. Gafchromic EBT3 film is placed inside a slab phantom to measure two-dimensional dose distributions. A benchmarked and comprehensive Monte Carlo model (Geant4) is established to directly compare with experiments.ResultsWith perpendicular fields a 5% narrowing of the beam FWHM and a 10 mm reduction in the 15% isodose penetration is seen for the 20 MeV beam. In the inline setup the penumbral width is reduced by up to 20%, and a local central dose enhancement of 100% is observed. Monte Carlo simulations are in agreement with the measured dose distributions (2% or 2 mm).ConclusionA new range of experiments have been performed to offer insight into how an electron beam mode could offer additional choices in MRI-linac radiotherapy. The work extends on historic studies to bring a successful unified experimental and Monte Carlo modeling approach for studying small field electron beam dosimetry inside magnetic fields. The results suggest further work, particularly on the inline magnetic field scenario.     

Experimental and theoretical studies of the spatial distribution of the ionization rate in a non-self-sustained discharge controlled by a fast electron beam

The spatial distribution of the current density of fast electrons and the ionization rate in a gap filled with atmospheric-pressure air under the conditions of a non-self-sustained discharge controlled by a fast electron beam were investigated. The experiments were carried out in a gas-discharge chamber with a grid electrode arranged in parallel to the exit window of the ionization source. Spatial variations in the current density of fast electrons resulting from the grid were measured. The propagation of the electron beam through the discharge system was simulated numerically by the Monte Carlo method in the so-called “effective collision” approximation. The calculated results agree well with the experimental data.     

Monte Carlo simulation and measurement of radiation leakage from applicators used in external electron radiotherapy

External electron radiotherapy is performed using a cone or applicator to collimate the beam. However, because of a trade-off between collimation and scattering/bremsstrahlung X-ray production, applicators generate a small amount of secondary radiation (leakage). We investigate the peripheral dose outside the radiation field of a Varian-type applicator. The dose and fluence outside the radiation field were analyzed in a detailed Monte Carlo simulation. The differences between the calculation results and data measured in a water phantom in an ionization chamber were less than ±1% in regions more than 3 mm below the surface of the phantom and at the depth of dose maximum. The calculated fluence was analyzed inside and outside the radiation field on a plane just above the water phantom surface. Changing the electron energy affected the off-axis fluence distribution outside the radiation field; however, the size of the applicator had little effect on this distribution. For each energy, the distributions outside the radiation field were similar to the dose distribution at shallow depths in the water phantom. The effect of secondary electrons generation by photon transmission through the alloy making up the lowest scraper was largest in the region from the field edge to directly below the cutout and at higher beam energies. The results of the Monte Carlo simulation confirm that the peripheral dose outside the field is significantly affected by radiation scattered or transmitted from the applicator, and the effect increases with the electron energy.     

Assessment of out-of-field doses in radiotherapy treatments of paediatric patients using Monte Carlo methods and measurements

PurposeTo assess out-of-field doses in radiotherapy treatments of paediatric patients, using Monte Carlo methods to implement a new model of the linear accelerator validated against measurements and developing a voxelized anthropomorphic paediatric phantom.MethodsCT images of a physical anthropomorphic paediatric phantom were acquired and a dosimetric planning using a TPS was obtained. The CT images were used to perform the voxelization of the physical phantom using the ImageJ software and later implemented in MCNP. In order to validate the Monte Carlo model, dose measurements of the 6 MV beam and Linac with 120 MLC were made in a clinical setting, using ionization chambers and a water phantom. Afterwards TLD measurements in the physical anthropomorphic phantom were performed in order to assess the out-of-field doses in the eyes, thyroid, c-spine, heart and lungs.ResultsThe Monte Carlo model was validated for in-field and out-of-field doses with average relative differences below 3%. The average relative differences between TLD measurements and Monte Carlo is 14,3% whilst the average relative differences between TLD and TPS is 55,8%. Moreover, organs up to 22.5 cm from PTV center show TLD and MCNP6 relative differences and TLD and TPS relative differences up to 21.2% and 92.0%, respectively.ConclusionsOur study provides a novel model that could be used in clinical research, namely in dose evaluation outside the treatment fields. This is particularly relevant, especially in pediatric patients, for studying new radiotherapy treatment techniques, since it can be used to estimate the development of secondary tumours.     

Polarity and ion recombination corrections in continuous and pulsed beams for ionization chambers with high Z chamber walls

In this work, the response of Farmer-type ionization chambers fitted with high atomic number (Z) walls is studied, and results of the effects of such walls on polarity and ion recombination correction factors in both continuous and pulsed beams are presented.Measurements were made in a continuous Co-60 beam and a pulsed 6 MV linac beam using an Exradin-A12 ionization chamber fitted with the manufacturer’s C-552 plastic wall, as well as geometrically identical walls made from aluminum, copper and molybdenum. The bias voltage was changed between 10 values (range: +50 to +560 V). Ion recombination was determined from Jaffé plots and by using the “two-voltage technique”. The saturation charge measured with each chamber wall was extrapolated from Jaffé plots. Additionally, the effect of different wall materials on chamber response was studied using MCNP simulations.Results showed that the polarity correction factor is not significantly affected by changes in chamber wall material (within 0.1%). Furthermore, although the saturation charges greatly vary with each chamber wall material, and charge multiplication increases for higher atomic number wall materials, the standard methods of calculating ion recombination yielded results that differed by only 0.2%. Therefore, polarity and ion recombination correction factors are not greatly affected by the chamber wall material. The experimental saturation charges for all the different wall materials agreed well within the uncertainty with MCNP simulations. The breakdown of the linear relationship in Jaffé plots that was previously reported to exist for conventional chamber walls was also observed with the different wall materials.     

Derivation of radiation quality average parameters in neutron-gamma radiation fields with the high-pressure ionization chamber: theory and practice

The established radiation quality parameters in mixed neutron-gamma radiation fields may be measured by applying the initial (columnar) recombination of ions in tissue-equivalent (TE) high-pressure ionization chambers (recombination chambers). The mean quality factor can be determined to within 10-15% for mixed fields with neutrons ranging from thermal to 10 MeV, and the dose mean LET of the proton component can be determined to within 10-15% if the gamma-ray absorbed dose fraction is known. These average parameters are derived by measuring the ratio of the ionization currents collected at two high-field strengths and constant gas pressure applied to the ionization chamber. By utilizing approximate correlations between physical parameters in the neutron energy region from thermal to 10 MeV, the dose mean LET of the heavy ion component, the overall dose mean LET, and the microdosimetric parameter y0,D of the mixed field can also be derived. Experimental verification of the method is presented for various neutron-gamma radiation spectra in air and in water by comparison to theoretical calculations and results from low-pressure proportional counter measurements. Good agreement is shown. The TE high-pressure ionization chamber appears to have wide potential for use as a dose-equivalent meter in radiation protection or as a beam characterization device in radiobiology.     

Monte Carlo determination of a nanoDot OSLD response using quality index for diagnostic kilovoltage X-ray beams

PurposeThis study aims to investigate the energy response of an optically stimulated luminescent dosimeter known as nanoDot for diagnostic kilovoltage X-ray beams via Monte Carlo calculations.MethodsThe nanoDot response is calculated as a function of X-ray beam quality in free air and on a water phantom surface using Monte Carlo simulations. The X-ray fluence spectra are classified using the quality index (QI), which is defined as the ratio of the effective energy to the maximum energy of the photons. The response is calculated for X-ray fluence spectra with QIs of 0.4, 0.5, and 0.6 with tube voltages of 50–137.6 kVp and monoenergetic photon beams. The surface dose estimated using the calculated response is verified by comparing it with that measured using an ionization chamber.ResultsThe nanoDot response in free air for monoenergetic photon beams (QI = 1.0) varies significantly at photon energies below 100 keV and reaches a factor of 3.6 at 25–30 keV. The response differs by up to approximately 6% between QIs of 0.4 and 0.6 for the same half-value layer (HVL). The response at the phantom surface decreases slightly owing to the backscatter effect, and it is almost independent of the field size. The agreement between the surface dose estimated using the nanoDot and that measured using the ionization chamber for assessing X-ray beam qualities is less than 2%.ConclusionsThe nanoDot response is indicated as a function of HVL for the specified QIs, and it enables the direct surface dose measurement.     

Assessment of ionization chamber correction factors in photon beams using a time saving strategy with PENELOPE code

The purpose of this study was to investigate Monte Carlo-based perturbation and beam quality correction factors for ionization chambers in photon beams using a saving time strategy with PENELOPE code. Simulations for calculating absorbed doses to water using full spectra of photon beams impinging the whole water phantom and those using a phase-space file previously stored around the point of interest were performed and compared. The widely used NE2571 ionization chamber was modeled with PENELOPE using data from the literature in order to calculate absorbed doses to the air cavity of the chamber. Absorbed doses to water at reference depth were also calculated for providing the perturbation and beam quality correction factors for that chamber in high energy photon beams. Results obtained in this study show that simulations with phase-space files appropriately stored can be up to ten times shorter than using a full spectrum of photon beams in the input-file. Values of k_Q and its components for the NE2571 ionization chamber showed good agreement with published values in the literature and are provided with typical statistical uncertainties of 0.2%. Comparisons to k_Q values published in current dosimetry protocols such as the AAPM TG-51 and IAEA TRS-398 showed maximum percentage differences of 0.1% and 0.6% respectively. The proposed strategy presented a significant efficiency gain and can be applied for a variety of ionization chambers and clinical photon beams.     

Measurement of the restricted dose mean LET outside the primary beam of a 60Co radiotherapy unit

A recently developed method for the direct measurement of the restricted dose mean LET (cutoff energy = 500 eV) of an unspecified photon or electron radiation field with the high-pressure ionization chamber has been utilized to investigate the variation of this radiobiologically important parameter outside the primary beam of a clinical 60Co unit. A small high-pressure tissue-equivalent ionization chamber was used, and its characteristics and experimental considerations for the present investigation are reported. Measurement of the restricted dose mean LET at the examined points outside the primary 60Co beam showed an increase of 50% with respect to the restricted dose mean LET of the uncollimated 60Co beam. No significant variation was noted with off-central axis distance, field size, wedge filter, or depth below Perspex slabs. Dose rates at the points of measurement outside the primary 60Co beam were 1-5% of the dose rate in the primary beam.     

Influence of the Presence of Tissue Expanders on Energy Deposition for Post-Mastectomy Radiotherapy

An increasing number of studies have shown that post-mastectomy radiotherapy presents benefits associated with the patients survival and a significant fraction of the treated patients makes use of tissue expanders for breast reconstruction. Some models of tissue expanders have a magnetic disk on their surface that constitutes heterogeneity in the radiation field, which can affect the dose distribution during the radiotherapy treatment. In this study, the influence of a metallic heterogeneity positioned in a breast tissue expander was evaluated by means of Monte Carlo simulations using the MCNPX code and using Eclipse treatment planning system. Deposited energy values were calculated in structures which have clinical importance for the treatment. Additionally, the effect in the absorbed energy due to backscattering and attenuation of the incident beam caused by the heterogeneity, as well as due to the expansion of the prosthesis, was evaluated in target structures for a 6 MV photon beam by simulations. The dose distributions for a breast treatment were calculated using a convolution/superposition algorithm from the Eclipse treatment planning system. When compared with the smallest breast expander volume, underdosage of 7% was found for the largest volume of breast implant, in the case of frontal irradiation of the chest wall, by Monte Carlo simulations. No significant changes were found in dose distributions for the presence of the heterogeneity during the treatment planning of irradiation with an opposed pair of beams. Even considering the limitation of the treatment planning system, the results obtained with its use confirm those ones found by Monte Carlo simulations for a tangent beam irradiation. The presence of a heterogeneity didńt alters the dose distributions on treatment structures. The underdosage of 7% observed with Monte Carlo simulations were found for irradiation at 0°, not used frequently in a clinical routine.     

Deterministic model for the transport of energetic particles: Application in the electron radiotherapy

A new deterministic method for calculating the dose distribution in the electron radiotherapy field is presented. The aim of this work was to validate our model by comparing it with the Monte Carlo simulation toolkit, GEANT4. A comparison of the longitudinal and transverse dose deposition profiles and electron distributions in homogeneous water phantoms showed a good accuracy of our model for electron transport, while reducing the calculation time by a factor of 50. Although the Bremsstrahlung effect is not yet implemented in our model, we propose here a method that solves the Boltzmann kinetic equation and provides a viable and efficient alternative to the expensive Monte Carlo modeling.     

Myosin regulatory domain orientation in skeletal muscle fibers: application of novel electron paramagnetic resonance spectral decomposition and molecular modeling methods

Reorientation of the regulatory domain of the myosin head is a feature of all current models of force generation in muscle. We have determined the orientation of the myosin regulatory light chain (RLC) using a spin-label bound rigidly and stereospecifically to the single Cys-154 of a mutant skeletal isoform. Labeled RLC was reconstituted into skeletal muscle fibers using a modified method that results in near-stoichiometric levels of RLC and fully functional muscle. Complex electron paramagnetic resonance spectra obtained in rigor necessitated the development of a novel decomposition technique. The strength of this method is that no specific model for a complex orientational distribution was presumed. The global analysis of a series of spectra, from fibers tilted with respect to the magnetic field, revealed two populations: one well-ordered (+/-15 degrees ) with the spin-label z axis parallel to actin, and a second population with a large distribution (+/-60 degrees ). A lack of order in relaxed or nonoverlap fibers demonstrated that regulatory domain ordering was defined by interaction with actin rather than the thick filament surface. No order was observed in the regulatory domain during isometric contraction, consistent with the substantial reorientation that occurs during force generation. For the first time, spin-label orientation has been interpreted in terms of the orientation of a labeled domain. A Monte Carlo conformational search technique was used to determine the orientation of the spin-label with respect to the protein. This in turn allows determination of the absolute orientation of the regulatory domain with respect to the actin axis. The comparison with the electron microscopy reconstructions verified the accuracy of the method; the electron paramagnetic resonance determined that axial orientation was within 10 degrees of the electron microscopy model.     

Characterization of Recombination Effects in a Liquid Ionization Chamber Used for the Dosimetry of a Radiosurgical Accelerator

Most modern radiation therapy devices allow the use of very small fields, either through beamlets in Intensity-Modulated Radiation Therapy (IMRT) or via stereotactic radiotherapy where positioning accuracy allows delivering very high doses per fraction in a small volume of the patient. Dosimetric measurements on medical accelerators are conventionally realized using air-filled ionization chambers. However, in small beams these are subject to nonnegligible perturbation effects. This study focuses on liquid ionization chambers, which offer advantages in terms of spatial resolution and low fluence perturbation. Ion recombination effects are investigated for the microLion detector (PTW) used with the Cyberknife system (Accuray). The method consists of performing a series of water tank measurements at different source-surface distances, and applying corrections to the liquid detector readings based on simultaneous gaseous detector measurements. This approach facilitates isolating the recombination effects arising from the high density of the liquid sensitive medium and obtaining correction factors to apply to the detector readings. The main difficulty resides in achieving a sufficient level of accuracy in the setup to be able to detect small changes in the chamber response.     

A comparison of symmetry and flatness measurements in small electron fields by different dosimeters in electron beam radiotherapy

BackgroundSymmetry and flatness are two quantities which should be evaluated in the commissioning and quality control of an electron beam in electron beam radiotherapy. The aim of this study is to compare symmetry and flatness obtained using three different dosimeters for various small and large fields in electron beam radiotherapy with linac.Materials and methodsBeam profile measurements were performed in a PTW water phantom for 10, 15 and 18 MeV electron beams of an Elekta Precise linac for small and large beams (1.5 × 1.5 cm² to 20 × 20 cm² field sizes). A Diode E detector and Semiflex-3D and Advanced Markus ionization chambers were used for dosimetry.ResultsBased on the obtained results, there are minor differences between the responses from different dosimeters (Diode E detector and Semiflex-3D and Advanced Markus ionization chambers) in measurement of symmetry and flatness for the electron beams. The symmetry and flatness values increase with increasing field size and electron beam energy for small and large field sizes, while the increases are minor in some cases.ConclusionsThe results indicate that the differences between the symmetry and flatness values obtained from the three dosimeter types are not practically important.     

Measurement of depth-dose of linear accelerator and simulation by use of Geant4 computer code

Radiation therapy is an established method of cancer treatment. New technologies in cancer radiotherapy need a more accurate computation of the dose delivered in the radiotherapy treatment plan. This study presents some results of a Geant4-based application for simulation of the absorbed dose distribution given by a medical linear accelerator (LINAC). The LINAC geometry is accurately described in the Monte Carlo code with use of the accelerator manufacturer''s specifications. The capability of the software for evaluating the dose distribution has been verified by comparisons with measurements in a water phantom; the comparisons were performed for percentage depth dose (PDD) and profiles for various field sizes and depths, for a 6-MV electron beam. Experimental and calculated dose values were in good agreement both in PDD and in transverse sections of the water phantom.     

Evaluation of ionization chamber stability checks using various sources

PurposeIt is important to check stability of ionization chambers in between regular calibration cycles. Stability checks can include individual ⁶⁰Co irradiations, use of a beta-emitting check source, or redundant measurements in megavoltage photon beams. While ⁶⁰Co irradiators are considered stable, they are rarely found in the clinical setting. Thus, this study seeks to compare the precision and efficiency in monitoring chamber stability using ⁹⁰Sr check sources and linear accelerator beams which are both commonly found in the clinical setting, and compare these sources to ⁶⁰Co.MethodsMeasurements were made with a ⁹⁰Sr beta-emitting check source and a 6 MV photon beam using a Constancy Check Phantom with three custom inserts to hold the ionization chambers. A comparison of both methods was performed with an Exradin A28 scanning chamber, Wellhofer IC69 Farmer-type chamber, and Exradin A12 Farmer-type chamber. Chamber stability was evaluated with individual charge readings and charge ratios among the three chambers. Results were compared to measurements taken in ⁶⁰Co with three Farmer-type chambers: the NEL 2571, PTW N30001G, and Exradin A12.ResultsStability of individual charge reading was found to be within ±1.0% for ⁹⁰Sr source measurements and ±0.5% for external beam measurements, including the ⁶⁰Co comparison. Additionally, the standard deviation of the mean charge ratios ranged from 0.15% to 0.40% for ⁹⁰Sr measurements and from 0.10% to 0.30% for the external beam measurements.ConclusionsThis work provides a comparison of techniques used to assess stability of ionization chambers in order to better inform the clinical physicist.     

Dosimetric impact of failing to apply correction factors to ion recombination in percentage depth dose measurements and the volume-averaging effect in flattening filter-free beams

PurposeTo examine whether it is essential to apply correction factors for ion recombination (k_S) to percentage depth dose (PDD) measurements and to the volume-averaging effect (k_vol) to ensure accurate absolute dose calibration for flattening filter-free (FFF) beams for the most commonly used ionization chambers.MethodsWe surveyed medical physicists worldwide (n = 159) to identify the five most common ionization chamber combinations used for absolute and relative reference dosimetry of FFF beams. We then assessed the overall absolute dose calibration error for FFF beams of the Artiste Siemens and TrueBeam Varian linear accelerators resulting from failing to apply correction factors k_S in the PDD(10) and the volume-averaging effect (k_vol) to such chamber combinations.ResultsAll the chamber combinations examined—the Farmer PTW 30013 ionization chamber used for absolute dosimetry, and the PTW 31010, PTW 30013, IBA CC04, IBA CC13, and PTW 31021 ionization chambers used for PDD curves measurements—showed non-negligible errors (≥0.5%). The largest error (1.6%) was found for the combination of the Farmer PTW 30013 chamber with the IBA CC13 chamber, which was the most widely used chamber combination in our survey.ConclusionsBased on our findings, we strongly recommend assessing the impact of failing to apply correction factors k_S in the PDD(10) and k_vol prior to using any chamber type for FFF beam reference dosimetry purposes.     

The track structure of alpha-particles from the point of view of ionization-cluster formation in "nanometric" volumes of nitrogen

Probability distributions of the size of ion clusters created in "nanometric" volumes of nitrogen by single alpha-particles of a gold-plated 241Am source, were measured and compared with those calculated by Monte Carlo methods in the same geometry. The diameter of the sensitive volumes had a mass per area of between 0.015 microgram/cm2 and 1.3 micrograms/cm2 which, for a material at unit density, corresponds to a nanometric target volume 0.15-13 nm in diameter. These nanometre sizes were simulated experimentally in a device called the Jet Counter. This consists of a pulse-operated valve which injects into an interaction chamber an expansion jet of molecular nitrogen gas, which is crossed by a narrow beam of alpha-particles. The resulting ions are counted and analyzed from the point of view of ionization cluster formation. The measured or calculated cluster size probabilities prove that the formation of ionization clusters along a "nanometre" track is governed by Poisson's law only in the case of very small target volumes, due to the contributions by secondary electrons. The present ionization cluster probabilities produced in "nanometric" volumes 0.15-13 nm in diameter, are the first ever determined experimentally and confirmed by Monte Carlo simulation.     

A semi-analytical model for calculating the penetration depth of a high energy electron beam in a water phantom with a magnetic field

PurposeAs an electron beam is incident on a uniform water phantom in the presence of a lateral magnetic field, the depth-dose distribution of the electron beam changes significantly and forms the well-known ‘Bragg peak’, with a depth-dose distribution similar to that of heavy ions. This phenomenon has pioneered a new field in the clinical application of electron beams. For such clinical applications, evaluating the penetration depth of electron beams quickly and accurately is the critical problem.MethodsThis paper describes a model for calculating the penetration depth of an electron beam rapidly and correctly in a water phantom under the influence of a magnetic field. The model was used to calculate the penetration depths under different conditions: the energies of electron beams of 6, 8, 12 and 15 MeV and the magnetic induction intensities of 0.75, 1.0, 1.5, 2.0 and 3.0 T. In addition, the calculation results were compared with the results of a Monte Carlo simulation.ResultsThe comparison results indicate that the difference between the two calculation methods was less than 0.5 cm. Moreover, the computing time of the calculation model was less than a second.ConclusionsThe semi-analytical model proposed in the present study enables the penetration depth of the electron beam in the presence of a magnetic field to be obtained with a computational efficiency higher than that of the Monte Carlo approach; thus, the proposed model has high potential for application.