Improvement in clinical step and shoot intensity modulated radiation therapy delivery accuracy on an integrated linear accelerator control system

Purpose: The dose delivery accuracy of 30 clinical step and shoot intensity modulated radiation therapy plans was investigated using the single integrated multileaf collimator controller of the Varian Truebeam linear accelerator (linac) (Varian Medical Systems, Palo Alto, CA) and compared with the dose delivery accuracy on a previous generation Varian 2100CD C-Series linac.

Methods and Materials: Ten prostate, 10 prostate and pelvic node, and 10 head-and-neck cases were investigated in this study. Dose delivery accuracy on each linac was assessed using Farmer ionization chamber point dose measurements, 2-dimensional planar ionization chamber array measurements, and the corresponding Varian dynamic log files. Absolute point dose measurements, fluence delivery accuracy, leaf position accuracy, and the overshoot effect were assessed for each plan.

Results: Absolute point dose delivery accuracy increased by 1.5% on the Truebeam compared with the 2100CD linac. No improvement in fluence delivery accuracy between the linacs, at a gamma criterion of 3%/3 mm was measured using the 2-dimensional ionization chamber array, with median (interquartile range) gamma passing rates of 98.99% (97.70%-99.72%) and 99.28% (98.26%-99.75%) for the Truebeam and 2100CD linacs, respectively. Varian log files also showed no improvement in fluence delivery between the linacs at 3%/3 mm, with median gamma passing rates of 99.97% (99.93%-99.99%) and 99.98% (99.94%-100%) for the Truebeam and 2100CD linacs, respectively. However, log files revealed improved leaf position accuracy and fluence delivery at 1%/1 mm criterion on the Truebeam (99.87%; 99.78%-99.94%) compared with the 2100CD linac (97.87%; 91.93%-99.49%). The overshoot effect, characterized on the 2100CD linac, was not observed on the Truebeam.

Conclusions: The integrated multileaf collimator controller on the Varian Truebeam improves clinical treatment delivery accuracy of step and shoot intensity modulated radiation therapy fields compared with delivery on a Varian C-series linac. © 2014.

Monitoring daily MLC positional errors using trajectory log files and EPID measurements for IMRT and VMAT deliveries

This work investigated the differences between multileaf collimator (MLC) positioning accuracy determined using either log files or electronic portal imaging devices (EPID) and then assessed the possibility of reducing patient specific quality control (QC) via phantom-less methodologies. In-house software was developed, and validated, to track MLC positional accuracy with the rotational and static gantry picket fence tests using an integrated electronic portal image. This software was used to monitor MLC daily performance over a 1 year period for two Varian TrueBeam linear accelerators, with the results directly compared with MLC positions determined using leaf trajectory log files. This software was validated by introducing known shifts and collimator errors. Skewness of the MLCs was found to be 0.03 ± 0.06° (mean ±1 standard deviation (SD)) and was dependent on whether the collimator was rotated manually or automatically. Trajectory log files, analysed using in-house software, showed average MLC positioning errors with a magnitude of 0.004 ± 0.003 mm (rotational) and 0.004 ± 0.011 mm (static) across two TrueBeam units over 1 year (mean ±1 SD). These ranges, as indicated by the SD, were lower than the related average MLC positioning errors of 0.000 ± 0.025 mm (rotational) and 0.000 ± 0.039 mm (static) that were obtained using the in-house EPID based software. The range of EPID measured MLC positional errors was larger due to the inherent uncertainties of the procedure. Over the duration of the study, multiple MLC positional errors were detected using the EPID based software but these same errors were not detected using the trajectory log files. This work shows the importance of increasing linac specific QC when phantom-less methodologies, such as the use of log files, are used to reduce patient specific QC. Tolerances of 0.25 mm have been created for the MLC positional errors using the EPID-based automated picket fence test. The software allows diagnosis of any specific leaf that needs repair and gives an indication as to the course of action that is required.

An EPID-based method for comprehensive verification of gantry, EPID and the MLC carriage positional accuracy in Varian linacs during arc treatments

Background

]In modern radiotherapy, it is crucial to monitor the performance of all linac components including gantry, collimation system and electronic portal imaging device (EPID) during arc deliveries. In this study, a simple EPID-based measurement method has been introduced in conjunction with an algorithm to investigate the stability of these systems during arc treatments with the aim of ensuring the accuracy of linac mechanical performance.

The Varian EPID sag, gantry sag, changes in source-to-detector distance (SDD), EPID and collimator skewness, EPID tilt, and the sag in MLC carriages as a result of linac rotation were separately investigated by acquisition of EPID images of a simple phantom comprised of 5 ball-bearings during arc delivery. A fast and robust software package was developed for automated analysis of image data. Twelve Varian linacs of different models were investigated.

The average EPID sag was within 1 mm for all tested linacs. All machines showed less than 1 mm gantry sag. Changes in SDD values were within 1.7 mm except for three linacs of one centre which were within 9 mm. Values of EPID skewness and tilt were negligible in all tested linacs. The maximum sag in MLC leaf bank assemblies was around 1 mm. The EPID sag showed a considerable improvement in TrueBeam linacs.

The methodology and software developed in this study provide a simple tool for effective investigation of the behaviour of linac components with gantry rotation. It is reproducible and accurate and can be easily performed as a routine test in clinics.

IRIDE:Interdisciplinary research infrastructure based on dual electron linacs and lasers

This paper describes the scientific aims and potentials as well as the preliminary technical design of IRIDE, an innovative tool for multi-disciplinary investigations in a wide field of scientific, technological and industrial applications. IRIDE will be a high intensity "particles factory", based on a combination of high duty cycle radio-frequency superconducting electron linacs and of high energy lasers. Conceived to provide unique research possibilities for particle physics, for condensed matter physics, chemistry and material science, for structural biology and industrial applications, IRIDE will open completely new research possibilities and advance our knowledge in many branches of science and technology. IRIDE is also supposed to be realized in subsequent stages of development depending on the assigned priorities. © 2013 Elsevier B.V.

ELM-induced transient tungsten melting in the JET divertor

The original goals of the JET ITER-like wall included the study of the impact of an all W divertor on plasma operation (Coenen et al 2013 Nucl. Fusion 53 073043) and fuel retention (Brezinsek et al 2013 Nucl. Fusion 53 083023). ITER has recently decided to install a full-tungsten (W) divertor from the start of operations. One of the key inputs required in support of this decision was the study of the possibility of W melting and melt splashing during transients. Damage of this type can lead to modifications of surface topology which could lead to higher disruption frequency or compromise subsequent plasma operation. Although every effort will be made to avoid leading edges, ITER plasma stored energies are sufficient that transients can drive shallow melting on the top surfaces of components. JET is able to produce ELMs large enough to allow access to transient melting in a regime of relevance to ITER.

Transient W melt experiments were performed in JET using a dedicated divertor module and a sequence of I-P = 3.0 MA/B-T = 2.9 T H-mode pulses with an input power of P-IN = 23 MW, a stored energy of similar to 6 MJ and regular type I ELMs at Delta W-ELM = 0.3 MJ and f(ELM) similar to 30 Hz. By moving the outer strike point onto a dedicated leading edge in the W divertor the base temperature was raised within similar to 1 s to a level allowing transient, ELM-driven melting during the subsequent 0.5 s. Such ELMs (delta W similar to 300 kJ per ELM) are comparable to mitigated ELMs expected in ITER (Pitts et al 2011 J. Nucl. Mater. 415 (Suppl.) S957-64).

Although significant material losses in terms of ejections into the plasma were not observed, there is indirect evidence that some small droplets (similar to 80 mu m) were released. Almost 1 mm (similar to 6 mm(3)) of W was moved by similar to 150 ELMs within 7 subsequent discharges. The impact on the main plasma parameters was minor and no disruptions occurred. The W-melt gradually moved along the leading edge towards the high-field side, driven by j x B forces. The evaporation rate determined from spectroscopy is 100 times less than expected from steady state melting and is thus consistent only with transient melting during the individual ELMs. Analysis of IR data and spectroscopy together with modelling using the MEMOS code Bazylev et al 2009 J. Nucl. Mater. 390-391 810-13 point to transient melting as the main process. 3D MEMOS simulations on the consequences of multiple ELMs on damage of tungsten castellated armour have been performed.

These experiments provide the first experimental evidence for the absence of significant melt splashing at transient events resembling mitigated ELMs on ITER and establish a key experimental benchmark for the MEMOS code.

X-ray laser-induced ablation of lead compounds

The recent commissioning of a X-ray free-electron laser triggered an extensive research in the area of X-ray ablation of high-Z, high-density materials. Such compounds should be used to shorten an effective attenuation length for obtaining clean ablation imprints required for the focused beam analysis. Compounds of lead (Z=82) represent the materials of first choice. In this contribution, single-shot ablation thresholds are reported for PbWO4 and PbI2 exposed to ultra-short pulses of extreme ultraviolet radiation and X-rays at FLASH and LCLS facilities, respectively. Interestingly, the threshold reaches only 0.11 J/cm(2) at 1.55 nm in lead tungstate although a value of 0.4 J/cm(2) is expected according to the wavelength dependence of an attenuation length and the threshold value determined in the XUV spectral region, i.e., 79 mJ/cm(2) at a FEL wavelength of 13.5 nm. Mechanisms of ablation processes are discussed to explain this discrepancy. Lead iodide shows at 1.55 nm significantly lower ablation threshold than tungstate although an attenuation length of the radiation is in both materials quite the same. Lower thermal and radiation stability of PbI2 is responsible for this finding.

The complex ion structure of warm dense carbon measured by spectrally resolved x-ray scattering

We present measurements of the complex ion structure of warm dense carbon close to the melting line at pressures around 100 GPa. High-pressure samples were created by laser-driven shock compression of graphite and probed by intense laser-generated x-ray sources with photon energies of 4.75 keV and 4.95 keV. High-efficiency crystal spectrometers allow for spectrally resolving the scattered radiation. Comparing the ratio of elastically and inelastically scattered radiation, we find evidence for a complex bonded liquid that is predicted by ab-initio quantum simulations showing the influence of chemical bonds under these conditions. Using graphite samples of different initial densities we demonstrate the capability of spectrally resolved x-ray scattering to monitor the carbon solid-liquid transition at relatively constant pressure of 150 GPa. Showing first single-pulse scattering spectra from cold graphite of unprecedented quality recorded at the Linac Coherent Light Source, we demonstrate the outstanding possibilities for future high-precision measurements at 4th Generation Light Sources.

Relativistic radiatively damped R-matrix calculation of the electron-impact excitation of W46+

The current design plans for the International Thermonuclear\nExperimental Reactor ( ITER) call for tungsten to be employed for\ncertain plasma facing components in the divertor region. Thus, accurate\natomic collision data are needed for emission modelling of tungsten.\nElectron-impact excitation and radiative rates are of particular\nimportance for Ni-like W, since this ion emits some of the most intense\nspectral lines of all ionization stages. We report on a fully\nrelativistic 115-level R-matrix calculations of W46+, which includes the\neffects of radiation damping. Although radiation damping is very\nimportant in most highly ionized species, its effects are reduced in\nthis case because of the closed-shell Ni-like ground state. The rates\nfrom these relativistic atomic calculations will be employed for\ncollisional-radiative modelling of this ion.

Non-equilibrium modeling of the FE XVII 3C/3D Line Ratio in an Intense X-ray Free-electron Laser Excited Plasma

Recent measurements using an X-ray Free Electron Laser (XFEL) and an Electron Beam Ion Trap at the Linac Coherent Light Source facility highlighted large discrepancies between the observed and theoretical values for the Fe XVII 3C/3D line intensity ratio. This result raised the question of whether the theoretical oscillator strengths may be significantly in error, due to insufficiencies in the atomic structure calculations. We present time-dependent spectral modeling of this experiment and show that non-equilibrium effects can dramatically reduce the predicted 3C/3D line intensity ratio, compared with that obtained by simply taking the ratio of oscillator strengths. Once these non-equilibrium effects are accounted for, the measured line intensity ratio can be used to determine a revised value for the 3C/3D oscillator strength ratio, giving a range from 3.0 to 3.5. We also provide a framework to narrow this range further, if more precise information about the pulse parameters can be determined. We discuss the implications of the new results for the use of Fe XVII spectral features as astrophysical diagnostics and investigate the importance of time-dependent effects in interpreting XFEL-excited plasmas.

Dielectronic recombination of W 35 +

Accurate data for dielectronic recombination (DR) of the ions of tungsten are of significant interest in the modelling of tungsten impurity transport and radiative power loss in current tokamaks and in ITER. However, the complexity of the atomic structure for many of these ions makes level-resolved DR calculations untenable on currently available computers, especially for open d- and f-subshell ions. The majority of DR data presently available for ITER modelling are based on an average-atom approximation. To improve upon these baseline calculations, we investigate the use of the configuration-average distorted-wave (CADW) method to calculate DR rate coefficients for complex open d-shell systems. The aim is to produce rate coefficients that are sufficiently accurate in terms of modelling, yet greatly reduced in term of computational complexity compared to level-resolved calculations. In this paper, we consider the DR of W 35 + . Initially, we carry out several large-scale level-resolved calculations for the DR associated with the 4d → 4f and 4p → 4d excitations in this ion, using both the level-resolved distorted-wave and Dirac R -matrix methods. These calculations allow us to test the validity of the CADW approach on these same excitations by comparing cross sections and rate coefficients. These comparisons demonstrate that the CADW method is relatively accurate in relation to these level-resolved methods for the temperature range for which W 35 + should exist in a collisionally ionized plasma. We then present results for CADW rate coefficients for both Δ n = 0 and Δ n = 1 excitations for this ion. This study indicates that it is now feasible to generate a much improved comprehensive set of DR data for the entire tungsten isonuclear sequence.

The effects of radiative cascades on the x-ray diagnostic lines of Fe16+

We present complete collisional-radiative modelling results for the soft x-ray emission lines of Fe16+ in the 15 Å–17 Å range. These lines have been the subject of much controversy in the astrophysical and laboratory plasma community. Radiative transition rates are generated from fully relativistic atomic structure calculations. Electron-impact excitation cross sections are determined using a fully relativistic R-matrix method employing 139 coupled atomic levels through n = 5. We find that, in all cases, using a simple ratio of the collisional rate coefficient times a radiative branching factor is not sufficient to model the widely used diagnostic line ratios. One has to include the effects of collisional-radiative cascades in a population model to achieve accurate line ratios. Our line ratio results agree well with several previous calculations and reasonably well with tokamak experimental measurements, assuming a Maxwellian electron-energy distribution. Our modelling results for four EBIT line ratios, assuming a narrow Gaussian electron-energy distribution, are in generally poor agreement with all four NIST measurements but are in better agreement with the two LLNL measurements. These results suggest the need for an investigation of the theoretical polarization calculations that are required to interpret the EBIT line ratio measurements.

Electron-impact ionization of the boron atom

Accurate knowledge of the electron-impact ionization of the B atom is urgently needed in current fusion plasma experiments to help design ITER wall components. Since no atomic measurements exist, nonperturba- tive time-dependent close-coupling (TDCC) calculations are carried out to accurately determine the direct ionization cross sections of the outer two subshells of B. Perturbative distorted-wave and semiempirical binary encounter calculations are found to yield cross sections from 26% lower to an order of magnitude higher than the current TDCC results. Unlike almost all neutral atoms, large excitation-autoionization contributions are found for the B atom. Nonperturbative R matrix with pseudostates (RMPS) calculations are also carried out to accurately determine the total ionization cross section of B. Previous 60 LS-term RMPS calculations are found to yield cross sections up to 40% higher than the current more extensive 476 LS-term RMPS results

The application of atomic physics within impurity diagnostics for fusion plasmas

With the focus of ITER on the transport and emission properties of tungsten, generating atomic data for complex species has received much interest. Focusing on impurity influx diagnostics, we discuss recent work on heavy species. Perturbative approaches do not work well for near neutral systems so non-perturbative data are required, presenting a particular challenge for these influx diagnostics. Recent results on Mo+ are given as an illustration of how the diagnostic applications can guide the theoretical calculations for such systems.

Characterisation of a two-dimensional liquid-filled ion chamber detector array using flattened and unflattened beams for small fields, small MUs and high dose-rates

In this study, the PTW 1000SRS array with Octavius 4D phantom was characterised for FF and FFF beams. MU linearity, field size, dose rate, dose per pulse (DPP) response and dynamic conformal arc treatment accuracy of the 1000SRS array were assessed for 6MV, 6FFF and 10FFF beams using a Varian TrueBeam STx linac. The measurements were compared with a pinpoint IC, microdiamond IC and EBT3 Gafchromic film. Measured dose profiles and FWHMs were compared with film measurements. Verification of FFF volumetric modulated arc therapy (VMAT) clinical plans were assessed using gamma analysis with 3%/3 mm and 2%/2 mm tolerances (10% threshold). To assess the effect of cross calibration dose rate, clinical plans with different dose rates were delivered and analysed. Output factors agreed with film measurements to within 4.5% for fields between 0.5 and 1 cm and within 2.7% for field sizes between 1.5 and 10 cm and were highly correlated with the microdiamond IC detector. Field sizes measured with the 1000SRS array were within 0.5 mm of film measurements. A drop in response of up to 1.8%, 2.4% and 5.2% for 6MV, 6FFF and 10FFF beams respectively was observed with increasing nominal dose rate. With an increase in DPP, a drop of up to 1.7%, 2.4% and 4.2% was observed in 6MV, 6FFF and 10FFF respectively. The differences in dose following dynamic conformal arc deliveries were less than 1% (all energies) from calculated. Delivered VMAT plans showed an average pass percentage of 99.5(±0.8)% and 98.4(±3.4)% with 2%/2 mm criteria for 6FFF and 10FFF respectively. A drop to 97.7(±2.2)% and 88.4(±9.6)% were observed for 6FFF and 10FFF respectively when plans were delivered at the minimum dose rate and calibrated at the maximum dose rate. Calibration using a beam with the average dose rate of the plan may be an efficient method to overcome the dose rate effects observed by the 1000SRS array.

Time-resolved dosimetric verification of respiratory-gated radiotherapy exposures using a high-resolution 2D ionisation chamber array

The aim of this work was to track and verify the delivery of respiratory-gated irradiations, performed with three versions of TrueBeam linac, using a novel phantom arrangement that combined the OCTAVIUS® SRS 1000 array with a moving platform. The platform was programmed to generate sinusoidal motion of the array. This motion was tracked using the real-time position management (RPM) system and four amplitude gating options were employed to interrupt MV beam delivery when the platform was not located within set limits. Time-resolved spatial information extracted from analysis of x-ray fluences measured by the array was compared to the programmed motion of the platform and to the trace recorded by the RPM system during the delivery of the x-ray field. Temporal data recorded by the phantom and the RPM system were validated against trajectory log files, recorded by the linac during the irradiation, as well as oscilloscope waveforms recorded from the linac target signal. Gamma analysis was employed to compare time-integrated 2D x-ray dose fluences with theoretical fluences derived from the probability density function for each of the gating settings applied, where gamma criteria of 2%/2 mm, 1%/1 mm and 0.5%/0.5 mm were used to evaluate the limitations of the RPM system. Excellent agreement was observed in the analysis of spatial information extracted from the SRS 1000 array measurements. Comparisons of the average platform position with the expected position indicated absolute deviations of  <0.5 mm for all four gating settings. Differences were observed when comparing time-resolved beam-on data stored in the RPM files and trajectory logs to the true target signal waveforms. Trajectory log files underestimated the cycle time between consecutive beam-on windows by 10.0  ±  0.8 ms. All measured fluences achieved 100% pass-rates using gamma criteria of 2%/2 mm and 50% of the fluences achieved pass-rates  >90% when criteria of 0.5%/0.5 mm were used. Results using this novel phantom arrangement indicate that the RPM system is capable of accurately gating x-ray exposure during the delivery of a fixed-field treatment beam.