Temporal characterization and in vitro comparison of cell survival following the delivery of 3D-conformal, intensity-modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT)

A phantom was designed and implemented for the delivery of treatment plans to cells in vitro. Single beam, 3D-conformal radiotherapy (3D-CRT) plans, inverse planned five-field intensity-modulated radiation therapy (IMRT), nine-field IMRT, single-arc volumetric modulated arc therapy (VMAT) and dual-arc VMAT plans were created on a CT scan of the phantom to deliver 3 Gy to the cell layer and verified using a Farmer chamber, 2D ionization chamber array and gafchromic film. Each plan was delivered to a 2D ionization chamber array to assess the temporal characteristics of the plan including delivery time and 'cell's eye view' for the central ionization chamber. The effective fraction time, defined as the percentage of the fraction time where any dose is delivered to each point examined, was also assessed across 120 ionization chambers. Each plan was delivered to human prostate cancer DU-145 cells and normal primary AGO-1522b fibroblast cells. Uniform beams were delivered to each cell line with the delivery time varying from 0.5 to 20.54 min. Effective fraction time was found to increase with a decreasing number of beams or arcs. For a uniform beam delivery, AGO-1552b cells exhibited a statistically significant trend towards increased survival with increased delivery time. This trend was not repeated when the different modulated clinical delivery methods were used. Less sensitive DU-145 cells did not exhibit a significant trend towards increased survival with increased delivery time for either the uniform or clinical deliveries. These results confirm that dose rate effects are most prevalent in more radiosensitive cells. Cell survival data generated from uniform beam deliveries over a range of dose rates and delivery times may not always be accurate in predicting response to more complex delivery techniques, such as IMRT and VMAT.

In-vitro investigation of out-of-field cell survival following the delivery of conformal, intensity-modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT) plans

The aim of this work is to determine the out-of-field survival of cells irradiated with either the primary field or scattered radiation in the presence and absence of intercellular communication following delivery of conformal, IMRT and VMAT treatment plans. Single beam, conformal, IMRT and VMAT plans were created to deliver 3 Gy to half the area of a T80 flask containing either DU-145 or AGO-1522 cells allowing intercellular communication between the in-and out-of-field cell populations. The same plans were delivered to a similar custom made phantom used to hold two T25 culture flasks, one flask in-field and one out-of-field to allow comparison of cell survival responses when intercellular communication is physically inhibited. Plans were created for the delivery of 8 Gy to the more radio-resistant DU-145 cells only in the presence and absence of intercellular communication. Cell survival was determined by clonogenic assay. In both cell lines, the out-of-field survival was not statistically different between delivery techniques for either cell line or dose. There was however, a statistically significant difference between survival out-of-field when intercellular communication was intact (single T80 culture flask) or inhibited (multiple T25 culture flasks) to in-field for all plans. No statistically significant difference was observed in-field with or without cellular communication to out-of-field for all plans. These data demonstrate out-of-field effects as important determinants of cell survival following exposure to modulated irradiation fields when cellular communication between differentially irradiated cell populations is present. This data is further evidence that refinement of existing radiobiological models to include indirect cell killing effects is required.

Use of arc-length method for capturing mode jumping in postbuckling aerostructures

The arc-length method has become a widely established solution technique for studying nonlinear structural behavior. By augmenting the set of nonlinear equilibrium equations with a constraint equation, which is a function of both the displacements and load increment, it is capable of traversing limit points. Numerous investigations have shown that highly nonlinear behavior such as sharp "snap-backs" can still lead to numerical difficulties. Two practical examples are presented to assess the effectiveness of this solution technique in capturing secondary instabilities in postbuckling structures, which present themselves as abrupt mode jumps. Although the first example poses no special difficulties, in the second case the nonlinear procedure fails to converge. An improvement to the method's formulation is suggested, which accounts for the residual forces that are usually neglected, when proceeding to the next increment once convergence is reached on the current increment. The choice of a correct load increment at the first iteration, within a predictor-corrector scheme, is central to the method's effectiveness. Current strategies for a choice of this load increment are discussed and are shown to be no longer consistent with the modified formulation; therefore, a new approach is proposed.

On thin ice: surface order and disorder during pre-melting

The effect of temperature on the structure of the ice Ih (0001) surface is considered through a series of molecular dynamics simulations on an ice slab. At relatively low temperatures (200K) a small fraction of surface self-interstitials (i.e. admolecules) appear that are formed exclusively from molecules leaving the outermost bilayer. At higher temperatures (ca. 250 K), vacancies start to appear in the inner part of the outermost bilayer exposing the underlying bilayer and providing sites with a high concentration of the dangling hydrogen bonds. Around 250-260 K aggregates of molecules formed on top of the outermost bilayer from self-interstitials become more mobile and have diffusivities approaching that of liquid water. At similar to 270-280 K the inner bilayer of one surface noticeably destructures and it appears that at above 285 K both surfaces are melting. The observed disparity in the onset of melting between the two sides of the slab is rationalised by considering the relationship between surface energy and the spatial distribution of protons at the surface; thermodynamic stability is conferred on the surface by maximising separations between dangling protons at the crystal exterior. Local hotspots associated with a high dangling proton density are suggested to be susceptible to pre-melting and may be more efficient at trapping species at the external surface than regions with low concentrations of protons thus potentially helping ice particles to catalyse reactions. A preliminary conclusion of this work is that only about 10-20 K below the melting temperature of the particular water potential employed is major disruption of the crystalline lattice noted which could be interpreted as being "liquid", the thickness of this film being about a nanometre.

Inverse planned constant dose rate volumetric modulated arc therapy (VMAT) as an efficient alternative to five-field intensity modulated radiation therapy (IMRT) for prostate

Purpose: The aim of this work was to determine if volumetric modulated arc therapy (VMAT) plans, created for constant dose-rate (cdrVMAT) delivery are a viable alternative to step and shoot five-field intensity modulated radiation therapy (IMRT). Materials and methods: The cdrVMAT plans, inverse planned on a treatment planning system with no solution to account for couch top or rails, were created for delivery on a linear accelerator with no variable dose rate control system. A series of five-field IMRT and cdrVMAT plans were created using dual partial arcs (gantry rotating between 260° and 100°) with 4° control points for ten prostate patients with the average rectal constraint incrementally increased. Pareto fronts were compared for the planning target volume homogeneity and average rectal dose between the two techniques for each patient. Also investigated were tumour control probability and normal tissue complication probability values for each technique. The delivery parameters [monitor units (MU) and time] and delivery accuracy of the IMRT and VMAT plans were also compared. Results: Pareto fronts showed that the dual partial arc plans were superior to the five-field IMRT plans, particularly for the clinically acceptable plans where average rectal doses were less for rotational plans (p = 0·009) with no statistical difference in target homogeneity. The cdrVMAT plans had significantly more MU (p = 0·005) but the average delivery time was significantly less than the IMRT plans by 42%. All clinically acceptable cdrVMAT plans were accurate in their delivery (gamma 99·2 ± 1·1%, 3%3 mm criteria). Conclusions Accurate delivery of dual partial arc cdrVMAT avoiding the couch top and rails has been demonstrated.

Magma Evolution in the Primitive, Intra-oceanic Tonga Arc: Petrogenesis of Basaltic Andesites at Tofua Volcano

Tofua volcano is situated midway along the Tonga oceanic arc and has undergone two phases of ignimbrite-forming activity. The eruptive products are almost entirely basaltic andesites (52 center dot 5-57 wt % SiO2) with the exception of a volumetrically minor pre-caldera dacite. The suite displays a strong tholeiitic trend with K2O <1 wt %. Phenocryst assemblages typically comprise plagioclase + clinopyroxene +/- orthopyroxene with microlites of Ti-magnetite. Olivine (Fo(83-88)) is rare and believed to be dominantly antecrystic. An increase in the extent and frequency of reverse zoning in phenocrysts, sieve-textured plagioclase and the occurrence of antecrystic phases in post-caldera lavas record a shift to dynamic conditions, allowing the interaction of magma batches that were previously distinct. Pyroxene thermobarometry suggests crystallization at 950-1200 degrees C and 0 center dot 8-1 center dot 8 kbar. Volatile measurements of glassy melt inclusions indicate a maximum H2O content of 4 center dot 16 wt % H2O, and CO2-H2O saturation curves indicate that crystallization occurred at two levels, at depths of 4-5 center dot 5 km and 1 center dot 5-2 center dot 5 km. Major and trace element models suggest that the compositions of the majority of the samples represent a differentiation trend whereby the dacite was produced by 65% fractional crystallization of the most primitive basaltic andesite. Trace element models suggest that the sub-arc mantle source is the residuum of depleted Indian mid-ocean ridge basalt mantle (IDMM-1% melt), whereas radiogenic isotope data imply addition of 0 center dot 2% average Tongan sediment melt and a fluid component derived from the subducted altered Pacific oceanic crust. A horizontal array on the U-Th equiline diagram and Ra excesses of up to 500% suggest fluid addition to the mantle wedge within the last few thousand years. Time-integrated (Ra-226/Th-230) vs Sr/Th and Ba/Th fractionation models imply differentiation timescales of up to 4500 years for the dacitic magma compositions at Tofua.

Recent contribution of sediments and fluids to the mantle's volatile budget

Subduction modifies the cycling of Earth's volatile elements. Fluid-rich sediments and hydrated oceanic lithosphere enter the convecting mantle at subduction zones. Some of the sediments and volatile components are released from the subducting slab, promote mantle melting and are returned to the surface by volcanism. The remainder continue into the deeper mantle. Quantification of the fate of these volatiles requires an understanding of both the nature and timing of fluid release and mantle melting(1). Here we analyse the trace element and isotopic geochemistry of fragments of upper mantle rocks that were transported to the surface by volcanic eruptions above the Batan Island subduction zone, Philippines. We find that the mantle fragments exhibit extreme disequilibrium between their U-Th-Ra isotopic ratios, which we interpret to result from the interaction of wet sediment melts and slab-derived fluids with rocks in the overlying mantle wedge. We infer that wet sediments were delivered from the slab to the mantle wedge between 8,000 and 10,000 years ago, whereas aqueous fluids were delivered separately much later. We estimate that about 625 ppm of water is retained in the wedge. A significant volume of water could therefore be delivered to the mantle transition zone at the base of the upper mantle, or even to the deeper mantle.

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.