Automating the Treatment Planning Process for Volumetric Modulated Arc Therapy Craniospinal Irradiation.

PURPOSE: The purpose of this work is to develop a method to automate the treatment planning process of craniospinal irradiation (CSI) using volumetric modulated arc therapy. METHODS AND MATERIALS: Two scripts were developed using the Eclipse Scripting Application Programming Interface to perform auto-plan preparation and optimization. Ten patients (age, 5-44 years) previously treated at our institution with low dose volumetric modulated arc therapy CSI (prescription of 12 Gy) before total body irradiation were selected to evaluate the efficacy of the proposed auto-planning process. Paired t tests compared the dosimetric indices of the auto-plans to the manually generated clinical plans. All plans were normalized to 95% of planning target volume (PTV) coverage with the prescription dose. Two physicians and one physicist were asked to evaluate the manual plans and auto-plans of each patient in a blinded retrospective review and to indicate clinical acceptability and which plans were preferred for treatment. RESULTS: Compared with the manual CSI plans, the auto plans obtained significant reductions in Dmean to the parotids, submandibular glands, larynx, thyroid, and significant reduction in the plan PTV Dmax and D0.03 cc. The standard deviation range of the dosimetric parameters was greatly reduced for auto plans (range, 0.1-1.3 Gy) relative to manual plans (range, 0.4-5.9 Gy) indicating better plan consistency. Among the 10 patients, the auto-plans were preferred over the manual plans 90% of the time by the reviewing experts. The required time for auto-planning was approximately 1 hour compared with estimated 4 or more hours for manual planning. CONCLUSIONS: Reductions in planning time without sacrifices in plan quality were obtained using the auto-planning process compared with manual planning. Variation in plan quality was also reduced. The auto-planning scripts will be made freely available to other institutions and clinics.

Tungsten Filled 3-Dimensional Printed Lung Blocks for Total Body Irradiation.

PURPOSE: Lung blocks for total-body irradiation are commonly used to reduce lung dose and prevent radiation pneumonitis. Currently, molten Cerrobend containing toxic materials, specifically lead and cadmium, is poured into molds to construct blocks. We propose a streamlined method to create 3-dimensional (3D)-printed lung block shells and fill them with tungsten ball bearings to remove lead and improve overall accuracy in the block manufacturing workflow. METHODS AND MATERIALS: 3D-printed lung block shells were automatically generated using an inhouse software, printed, and filled with 2 to 3 mm diameter tungsten ball bearings. Clinical Cerrobend blocks were compared with the physician drawn blocks as well as our proposed tungsten filled 3D-printed blocks. Physical and dosimetric comparisons were performed on a linac. Dose transmission through the Cerrobend and 3D-printed blocks were measured using point dosimetry (ion-chamber) and the on-board Electronic-Portal-Imaging-Device (EPID). Dose profiles from the EPID images were used to compute the full-width-half-maximum and to compare with the treatment-planning-system. Additionally, the coefficient-of-variation in the central 80% of full-width-half-maximum was computed and compared between Cerrobend and 3D-printed blocks. RESULTS: The geometric difference between treatment-planning-system and 3D-printed blocks was significantly lower than Cerrobend blocks (3D: -0.88 ± 2.21 mm, Cerrobend: -2.28 ± 2.40 mm, P = .0002). Dosimetrically, transmission measurements through the 3D-printed and Cerrobend blocks for both ion-chamber and EPID dosimetry were between 42% to 48%, compared with the open field. Additionally, coefficient-of-variation was significantly higher in 3D-printed blocks versus Cerrobend blocks (3D: 4.2% ± 0.6%, Cerrobend: 2.6% ± 0.7%, P < .0001). CONCLUSIONS: We designed and implemented a tungsten filled 3D-printed workflow for constructing total-body-irradiation lung blocks, which serves as an alternative to the traditional Cerrobend based workflow currently used in clinics. This workflow has the capacity of producing clinically useful lung blocks with minimal effort to facilitate the removal of toxic materials from the clinic.

Radiation Therapy Without Anesthesia for a 2-Year-Old Child Using Audio-Visual Assisted Therapeutic Ambience in Radiation Therapy (AVATAR).

PURPOSE: Radiation therapy (RT) is essential to managing many pediatric malignancies but can provoke anxiety, fear, and discomfort for children owing to prolonged treatment time, extended course, and restrictive immobilization. Patients younger than 10 years frequently require daily general anesthesia (GA), which is resource intensive, expensive, potentially toxic, and anxiety and fear provoking. Audio-Visual Assisted Therapeutic Ambience in Radiation Therapy (AVATAR), a video streaming device, has been proposed as an alternative to anesthesia in patients aged 3 to 10 years. A pilot study evaluating the efficacy of this novel innovation is accruing, but patients younger than 3 years are ineligible. METHODS AND MATERIALS: We simulated a 2-year-old with stage IV Wilms tumor for bilateral whole-lung and left-flank irradiation without GA. Using AVATAR, we attempted to deliver RT to this patient without sedation. Patient anxiety at the time of simulation and at the beginning, middle, and end of the treatment course was characterized using the validated Modified Yale Preoperative Anxiety Score (mYPAS) measurement tool. RESULTS: Although the patient tolerated computed tomography simulation without GA or AVATAR use, his mYPAS of 14 out of 18 indicated significant anxiety. Using AVATAR, all treatments were delivered without GA; his mYPASs were 5 and 4 (the lowest possible) and 4 at the first, midcourse, and final treatments, indicating no significant anxiety and a decrease from the pre-AVATAR baseline. Without GA, the time to deliver RT decreased by 66% from 90 to 30 minutes. CONCLUSIONS: We describe an expanded, previously unreported indication for AVATAR by demonstrating the feasibility of this approach to reduce or omit anesthesia in appropriate younger patients currently excluded from ongoing trials. The financial and quality-of-life benefits (including decreased stress, anxiety, toxic effects, cost, and appointment time) of AVATAR use may be extendable to a younger patient population than previously thought. In older children, prospective validation is ongoing, but additional study in patients younger than 3 years is needed.

A robotically assisted 3D printed quality assurance lung phantom for Calypso.

Purpose. Radiation dose delivered to targets located near the upper-abdomen or in the thorax are significantly affected by respiratory-motion. Relatively large-margins are commonly added to compensate for this motion, limiting radiation-dose-escalation. Internal-surrogates of target motion, such as a radiofrequency (RF) tracking system, i.e. Calypso®System, are used to overcome this challenge and improve normal-tissue sparing. RF tracking systems consist of implanting transponders in the vicinity of the tumor to be tracked using radiofrequency-waves. Unfortunately, although the manufacture provides a universal quality-assurance (QA) phantom, QA-phantoms specifically for lung-applications are limited, warranting the development of alternative solutions to fulfil the tests mandated by AAPM's TG142. Accordingly, our objective was to design and develop a motion-phantom to evaluate Calypso for lung-applications that allows the Calypso®Beacons to move in different directions to better simulate truelung-motion.Methods and Materials.A Calypso lung QA-phantom was designed, and 3D-printed. The design consists of three independent arms where the transponders were attached. A pinpoint-chamber with a buildup-cap was also incorporated. A 4-axis robotic arm was programmed to drive the motion-phantom to mimic breathing. After acquiring a four-dimensional-computed-tomography (4DCT) scan of the motion-phantom, treatment-plans were generated and delivered on a Varian TrueBeam®with Calypso capabilities. Stationary and gated-treatment plans were generated and delivered to determine the dosimetric difference between gated and non-gated treatments. Portal cine-images were acquired to determine the temporal-accuracy of delivery by calculating the difference between the observed versus expected transponders locations with the known speed of the transponders' motion.Results.Dosimetric accuracy is better than the TG142 tolerance of 2%. Temporal accuracy is greater than, TG142 tolerance of 100 ms for beam-on, but less than 100 ms for beam-hold.Conclusions.The robotic QA-phantom designed and developed in this study provides an independent phantom for performing Calypso lung-QA for commissioning and acceptance testing of Calypso for lung treatments.

Precision radiotherapy using monochromatic inverse Compton x-ray sources.

PURPOSE: The dosimetric properties of inverse Compton (IC) x-ray sources were investigated to determine their utility for stereotactic radiation therapy. METHODS: Monte Carlo simulations were performed using the egs brachy user code of EGSnrc. Nominal IC source x-ray energies of 80 and 150 keV were considered in this work. Depth-dose and lateral dose profiles in water were calculated, as was dose enhancement in the bone. Further simulations were performed for brain and spine treatment sites. The impact of gold nanoparticle doping was also investigated for the brain treatment site. Analogous dose calculations were performed in a clinical treatment planning system using a clinical 6 MV photon beam model and were compared to the Monte Carlo simulations. RESULTS: Both 80 and 150 keV IC beams were observed to have sharp 80-20 penumbra (i.e., < 0.1 mm) with broad low-dose tails in water. For reference, the calculated penumbra for the 6 MV clinical beam was 3 mm. Maximum dose enhancement factors in bone of 3.1, 1.4, and 1.1 were observed for the 80, 150 keV, and clinical 6 MV beams, respectively. The plan quality for the single brain metastasis case was similar between the IC beams and the 6 MV beam without gold nanoparticles. As the concentration of gold within the target increased, the V12 Gy to the normal brain tissue and D max within the target volume significantly decreased and the conformity significantly improved, which resulted in superior plan quality over the clinical 6 MV beam plan. In the spine cases, the sharp penumbra and enhanced dose to bone of the IC beams produced superior plan quality (i.e., better conformity, normal tissue sparing, and spinal cord sparing) as compared to the clinical 6 MV beam plans. CONCLUSIONS: The findings from this work indicate that inverse Compton x-ray sources are well suited for stereotactic radiotherapy treatments due to their sharp penumbra and dose enhancement around high atomic number materials. Future work includes investigating the properties of intensity-modulated inverse Compton x-ray sources to improve the homogeneity within the target tissue.

Technical Note: Evaluation of audiovisual biofeedback smartphone application for respiratory monitoring in radiation oncology.

PURPOSE: Radiation dose delivered to targets located near the upper abdomen or thorax are significantly affected by respiratory motion, necessitating large margins, limiting dose escalation. Surrogate motion management devices, such as the Real-time Position Management (RPM™) system (Varian Medical Systems, Palo Alto, CA), are commonly used to improve normal tissue sparing. Alternative to current solutions, we have developed and evaluated the feasibility of a real-time position management system that leverages the motion data from the onboard hardware of Apple iOS devices to provide patients with visual coaching with the potential to improve the reproducibility of breathing as well as improve patient compliance and reduce treatment delivery time. METHODS AND MATERIALS: The iOS application, coined the Instant Respiratory Feedback (IRF) system, was developed in Swift (Apple Inc., Cupertino, CA) using the Core-Motion library and implemented on an Apple iPhone® devices. Operation requires an iPhone®, a three-dimensional printed arm, and a radiolucent projector screen system for feedback. Direct comparison between IRF, which leverages sensor fusion data from the iPhone®, and RPM™, an optical-based system, was performed on multiple respiratory motion phantoms and volunteers. The IRF system and RPM™ camera tracking marker were placed on the same location allowing for simultaneous data acquisition. The IRF surrogate measurement of displacement was compared to the signal trace acquired using RPM™ with univariate linear regressions and Bland-Altman analysis. RESULTS: Periodic motion shows excellent agreement between both systems, and subject motion shows good agreement during regular and irregular breathing motion. Comparison of IRF and RPM™ show very similar signal traces that were significantly related across all phantoms, including those motion with different amplitude and frequency, and subjects' waveforms (all r > 0.9, P < 0.0001). We demonstrate the feasibility of performing four-dimensional cone beam computed tomography using IRF which provided similar image quality as RPM™ when reconstructing dynamic motion phantom images. CONCLUSIONS: Feasibility of an iOS application to provide real-time respiratory motion is demonstrated. This system generated comparable signal traces to a commercially available system and offers an alternative method to monitor respiratory motion.

Intensity modulated Ir-192 brachytherapy using high-Z 3D printed applicators.

Gynecologic cancers are often asymmetric, yet current Ir-192 brachytherapy techniques provide only limited radial modulation of the dose. The shielded solutions investigated here solve this by providing the ability to modulate between highly asymmetric and radially symmetric dose distributions at a given location. To find applicator designs that can modulate between full dose and less than 50% dose, at the dimensions of the urethra, a 2D calculation algorithm was developed to narrow down the search space. Two shielding design types were then further investigated using Monte Carlo and Boltzmann-solver dose calculation algorithms. 3D printing techniques using ISO 10993 certified biocompatible plastics and 3D printable tungsten-loaded plastics were tested. It was also found that shadowing effects set by the shape of the shielding cannot be easily modulated out, hence careful design is required. The shielded applicator designs investigated here, allow for reduction of the dose by over 50% at 5 mm from the applicator surface in desired regions, while also allowing radially symmetric dose with isodose line deviations less than 0.5 mm from circular. The shielding designs were also chosen with treatment delivery time in mind. Treatment times for these shielded designs were found to be less than 1.4 times longer than a 6-channel unshielded cylinder for the equivalent fully symmetric dose distribution. The 2D calculation methods developed here provide a simple way to rapidly evaluate shielding designs, while the 3D printing techniques also allow for devices with novel shapes to be rapidly prototyped. Both TOPAS Monte Carlo and Acuros BV calculations show that significant dose shaping and organ at risk sparing can be achieved without significantly compromising the plan in regions that require the full dose.

An integrated quality assurance phantom for frameless single-isocenter multitarget stereotactic radiosurgery.

Brain stereotactic radiosurgery (SRS) treatments require multiple quality assurance (QA) procedures to ensure accurate and precise treatment delivery. As single-isocenter multitarget SRS treatments become more popular, the quantification of off-axis accuracy of the linear accelerator is crucial. In this study, a novel brain SRS integrated phantom was developed and validated to enable SRS QA with a single phantom to facilitate implementation of a frameless single-isocenter, multitarget SRS program. This phantom combines the independent verification of each positioning system, the Winston-Lutz, off-axis accuracy evaluation (i.e. off-axis Winston-Lutz), and the dosimetric accuracy utilizing both point dose measurements as well as film measurement, without moving the phantom. A novel 3D printed phantom, coined OneIso, was designed with a movable insert which can switch between the Winston-Lutz test target and dose measurement without moving the phantom itself. For dose verification, ten brain SRS clinical treatment plans with 10 MV flattening-filter-free beams were delivered on a Varian TrueBeam with a high-definition multileaf collimator (HD-MLC). Radiochromic film and pinpoint ion chamber comparison measurements were made between the OneIso and solid water (SW) phantom setups. For the off-axis Winston-Lutz measurements, a row of off-axis ball bearings (BBs) was integrated into the OneIso. To quantify the spatial accuracy versus distance from the isocenter, two-dimensional displacements were calculated between the planned and delivered BB locations relative to their respective MLC defined field border. OneIso and the SW phantoms agree within 1%, for both film and point dose measurements. OneIso identified a reduction in spatial accuracy further away from the isocenter. Differences increased as distance from the isocenter increased, exceeding recommended SRS accuracy tolerances at 7 cm away from the isocenter. OneIso provides a streamlined, single-setup workflow for single-isocenter multitarget frameless linac-based SRS QA. Additionally, with the ability to quantify off-axis spatial discrepancies, we can determine limitations on the maximum distance between targets to ensure a single-isocenter multitarget SRS program meets recommended guidelines.

Patient motion tracking for non-isocentric and non-coplanar treatments via fixed frame-of-reference 3D camera.

PURPOSE: As C-arm linac radiation therapy evolves toward faster, more efficient delivery, and more conformal dosimetry, treatments with increasingly complex couch motions are emerging. Monitoring the patient motion independently of the couch motion during non-coplanar, non-isocentric, or dynamic couch treatments is a key bottleneck to their clinical implementation. The goal of this study is to develop a prototype real-time monitoring system for unconventional beam trajectories to ensure a safe and accurate treatment delivery. METHODS: An in-house algorithm was developed for tracking using a couch-mounted three-dimensional (3D) depth camera. The accuracy of patient motion detection on the couch was tested on a 3D printed phantom created from the body surface contour exported from the treatment planning system. The technique was evaluated against a commercial optical surface monitoring system with known phantom displacements of 3, 5, and 7 mm in lateral, longitudinal, and vertical directions by placing a head phantom on a dynamic platform on the treatment couch. The stability of the monitoring system was evaluated during dynamic couch trajectories, at speeds between 10.6 and 65 cm/min. RESULTS: The proposed monitoring system agreed with the ceiling mounted optical surface monitoring system in longitudinal, lateral, and vertical directions within 0.5 mm. The uncertainty caused by couch vibration increased with couch speed but remained sub-millimeter for speeds up to 32 cm/min. For couch speeds of 10.6, 32.2, and 65 cm/min, the uncertainty ranges were 0.27- 0.73 mm, 0.15-0.87 mm, and 0.28-1.29 mm, respectively. CONCLUSION: By mounting a 3D camera in the same frame-of-reference as the patient and eliminating dead spots, this proof of concept demonstrates real-time patient monitoring during couch motion. For treatments with non-coplanar beams, multiple isocenters, or dynamic couch motion, this provides additional safety without additional radiation dose and avoids some of the complexity and limitations of room mounted systems.

Cost Analysis of Audiovisual-Assisted Therapeutic Ambiance in Radiation Therapy (AVATAR)-Aided Omission of Anesthesia in Radiation for Pediatric Malignancies.

PURPOSE: Because children cannot reliably remain immobile during radiation therapy (RT) for cancer anatomy targeting requiring millimeter precision, daily anesthesia plays a large role in each RT session. Unfortunately, anesthesia is a source of financial burden for patients' families and is invasive and traumatic. This study attempts to assess the cost-savings benefit of audiovisual-assisted therapeutic ambiance in radiation therapy (AVATAR)-aided omission of pediatric anesthesia in RT. METHODS AND MATERIALS: The baseline time of anesthesia during RT was derived from documented anesthesia billing time during RT simulation at our institution and from the published literature. Current Procedural Terminology and relative value unit codes encompassing anesthesia-related charges from radiation oncology and anesthesia were analyzed in concert with this value to calculate the total cost of pediatric anesthesia per RT session. RESULTS: The mean number of RT fractions administered per patient with AVATAR-directed anesthesia omission at our institution was 19.0, similar to the 17.6 previously reported. At a mean anesthesia time exceeding 30 minutes (with mean RT duration of 4 weeks), the cost of pediatric anesthesia per RT fraction in non-AVATAR sessions was $1,904.35, yielding a total RT treatment anesthesia cost of $38,087.00 per patient (including simulation). Patients at our institution were not billed for AVATAR-assisted RT. CONCLUSIONS: The ability of AVATAR to obviate the need for daily anesthesia in pediatric RT provides substantial cost-savings. These findings argue for increased utilization of AVATAR and for analyses of RT targeting the accuracy of AVATAR versus conventional anesthesia-guided treatment of pediatric malignancies.

Factor 10 Expedience of Monthly Linac Quality Assurance via an Ion Chamber Array and Automation Scripts.

PURPOSE: While critical for safe and accurate radiotherapy, monthly quality assurance of medical linear accelerators is time-consuming and takes physics resources away from other valuable tasks. The previous methods at our institution required 5 hours to perform the mechanical and dosimetric monthly linear accelerator quality assurance tests. An improved workflow was developed to perform these tests with higher accuracy, with fewer error pathways, in significantly less time. METHODS: A commercial ion chamber array (IC profiler, Sun Nuclear, Melbourne, Florida) is combined with automation scripts to consolidate monthly linear accelerator QA. The array was used to measure output, flatness, symmetry, jaw positions, gated dose constancy, energy constancy, collimator walkout, crosshair centering, and dosimetric leaf gap constancy. Treatment plans were combined with automation scripts that interface with Sun Nuclear's graphical user interface. This workflow was implemented on a standard Varian clinac, with no special adaptations, and can be easily applied to other C-arm linear accelerators. RESULTS: These methods enable, in 30 minutes, measurement and analysis of 20 of the 26 dosimetric and mechanical monthly tests recommended by TG-142. This method also reduces uncertainties in the measured beam profile constancy, beam energy constancy, field size, and jaw position tests, compared to our previous methods. One drawback is the increased uncertainty associated with output constancy. Output differences between IC profiler and farmer chamber in plastic water measurements over a 6-month period, across 4 machines, were found to have a 0.3% standard deviation for photons and a 0.5% standard deviation for electrons, which is sufficient for verifying output accuracy according to TG-142 guidelines. To minimize error pathways, automation scripts which apply the required settings, as well as check the exported data file integrity were employed. CONCLUSIONS: The equipment, procedure, and scripts used here reduce the time burden of routine quality assurance tests and in most instances improve precision over our previous methods.

The Structure of Liquid and Amorphous Hafnia.

Understanding the atomic structure of amorphous solids is important in predicting and tuning their macroscopic behavior. Here, we use a combination of high-energy X-ray diffraction, neutron diffraction, and molecular dynamics simulations to benchmark the atomic interactions in the high temperature stable liquid and low-density amorphous solid states of hafnia. The diffraction results reveal an average Hf-O coordination number of ~7 exists in both the liquid and amorphous nanoparticle forms studied. The measured pair distribution functions are compared to those generated from several simulation models in the literature. We have also performed ab initio and classical molecular dynamics simulations that show density has a strong effect on the polyhedral connectivity. The liquid shows a broad distribution of Hf-Hf interactions, while the formation of low-density amorphous nanoclusters can reproduce the sharp split peak in the Hf-Hf partial pair distribution function observed in experiment. The agglomeration of amorphous nanoparticles condensed from the gas phase is associated with the formation of both edge-sharing and corner-sharing HfO6,7 polyhedra resembling that observed in the monoclinic phase.

Beyond sixfold coordinated Si in SiO2 glass at ultrahigh pressures.

We investigated the structure of SiO2 glass up to 172 GPa using high-energy X-ray diffraction. The combination of a multichannel collimator with diamond anvil cells enabled the measurement of structural changes in silica glass with total X-ray diffraction to previously unachievable pressures. We show that SiO2 first undergoes a change in Si-O coordination number from fourfold to sixfold between 15 and 50 GPa, in agreement with previous investigations. Above 50 GPa, the estimated coordination number continuously increases from 6 to 6.8 at 172 GPa. Si-O bond length shows first an increase due to the fourfold to sixfold coordination change and then a smaller linear decrease up to 172 GPa. We reconcile the changes in relation to the oxygen-packing fraction, showing that oxygen packing decreases at ultrahigh pressures to accommodate the higher than sixfold Si-O coordination. These results give experimental insight into the structural changes of silicate glasses as analogue materials for silicate melts at ultrahigh pressures.

Probing the triplet correlation function in liquid water by experiments and molecular simulations.

Despite very significant developments in scattering experiments like X-ray and neutron diffraction, it has been challenging to elucidate the nature of tetrahedral molecular configurations in liquid water. A key question is whether the pair correlation functions, which can be obtained from scattering experiments, are sufficient to describe the tetrahedral ordering of water molecules. In our previous study (Dhabal et al., J. Chem. Phys., 2014, 141, 174504), using data-sets generated from reverse Monte Carlo and molecular dynamics simulations, we showed that the triplet correlation functions contain important information on the tetrahedrality of water in the liquid state. In the present study, X-ray scattering experiments and molecular dynamics (MD) simulations are used to link the isothermal pressure derivative of the structure factor with the triplet correlation functions for water. Triplet functions are determined for water up to 3.3 kbar at 298 K to display the effect of pressure on the water structure. The results suggest that triplet functions (H[combining tilde](q)) obtained using a rigid-body TIP4P/2005 water model are consistent with the experimental results. The triplet functions obtained in experiment as well as in simulations evince that in the case of tetrahedral liquids, exertion of higher pressure leads to a better agreement with the Kirkwood superposition approximation (KSA). We further validate this observation using the triplet correlation functions (g(3)(r,s,t)) calculated directly from simulation trajectory, revealing that both H[combining tilde](q) in q-space and g(3)(r,s,t) in real-space contain similar information on the tetrahedrality of liquids. This study demonstrates that the structure factor, even though it has only pair correlation information of the liquid structure, can shed light on three-body correlations in liquid water through its isothermal pressure derivative term.

The temperature dependence of intermediate range oxygen-oxygen correlations in liquid water.

We analyze the recent temperature dependent oxygen-oxygen pair-distribution functions from experimental high-precision x-ray diffraction data of bulk water by Skinner et al. [J. Chem. Phys. 141, 214507 (2014)] with particular focus on the intermediate range where small, but significant, correlations are found out to 17 Å. The second peak in the pair-distribution function at 4.5 Å is connected to tetrahedral coordination and was shown by Skinner et al. to change behavior with temperature below the temperature of minimum isothermal compressibility. Here we show that this is associated also with a peak growing at 11 Å which strongly indicates a collective character of fluctuations leading to the enhanced compressibility at lower temperatures. We note that the peak at ∼13.2 Å exhibits a temperature dependence similar to that of the density with a maximum close to 277 K or 4 °C. We analyze simulations of the TIP4P/2005 water model in the same manner and find excellent agreement between simulations and experiment albeit with a temperature shift of ∼20 K.

Atomistic insight into viscosity and density of silicate melts under pressure.

A defining characteristic of silicate melts is the degree of polymerization (tetrahedral connectivity), which dictates viscosity and affects compressibility. While viscosity of depolymerized silicate melts increases with pressure consistent with the free-volume theory, isothermal viscosity of polymerized melts decreases with pressure up to ~3-5 GPa, above which it turns over to normal (positive) pressure dependence. Here we show that the viscosity turnover in polymerized liquids corresponds to the tetrahedral packing limit, below which the structure is compressed through tightening of the inter-tetrahedral bond angle, resulting in high compressibility, continual breakup of tetrahedral connectivity and viscosity decrease with increasing pressure. Above the turnover pressure, silicon and aluminium coordination increases to allow further packing, with increasing viscosity and density. These structural responses prescribe the distribution of melt viscosity and density with depth and play an important role in magma transport in terrestrial planetary interiors.

Benchmark oxygen-oxygen pair-distribution function of ambient water from x-ray diffraction measurements with a wide Q-range.

Four recent x-ray diffraction measurements of ambient liquid water are reviewed here. Each of these measurements represents a significant development of the x-ray diffraction technique applied to the study of liquid water. Sources of uncertainty from statistical noise, Q-range, Compton scattering, and self-scattering are discussed. The oxygen-hydrogen contribution to the measured x-ray scattering pattern was subtracted using literature data to yield an experimental determination, with error bars, of the oxygen-oxygen pair-distribution function, g(OO)(r), which essentially describes the distribution of molecular centers. The extended Q-range and low statistical noise of these measurements has significantly reduced truncation effects and related errors in the g(OO)(r) functions obtained. From these measurements and error analysis, the position and height of the nearest neighbor maximum in g(OO)(r) were found to be 2.80(1) Å and 2.57(5) respectively. Numerical data for the coherent differential x-ray scattering cross-section I(X)(Q), the oxygen-oxygen structure factor S(OO)(Q), and the derived g(OO)(r) are provided as benchmarks for calibrating force-fields for water.

Structure of the floating water bridge and water in an electric field.

The floating water bridge phenomenon is a freestanding rope-shaped connection of pure liquid water, formed under the influence of a high potential difference (approximately 15 kV). Several recent spectroscopic, optical, and neutron scattering studies have suggested that the origin of the bridge is associated with the formation of anisotropic chains of water molecules in the liquid. In this work, high energy X-ray diffraction experiments have been performed on a series of floating water bridges as a function of applied voltage, bridge length, and position within the bridge. The two-dimensional X-ray scattering data showed no direction-dependence, indicating that the bulk water molecules do not exhibit any significant preferred orientation along the electric field. The only structural changes observed were those due to heating, and these effects were found to be the same as for bulk water. These X-ray scattering measurements are supported by molecular dynamics (MD) simulations which were performed under electric fields of 10(6) V/m and 10(9) V/m. Directional structure factor calculations were made from these simulations parallel and perpendicular to the E-field. The 10(6) V/m model showed no significant directional-dependence (anisotropy) in the structure factors. The 10(9) V/m model however, contained molecules aligned by the E-field, and had significant structural anisotropy.

Comment on 'Molecular arrangement in water: random but not quite'.

Accurate high energy x-ray diffraction data are presented on liquid water measured at room temperature. Sources of both systematic and statistical errors within the experiment are considered and data consistency checks are discussed. It is found that the resulting x-ray pair distribution function is smoothly varying in real space and shows no evidence of small peaks in the 3-5 Å region. Our results are in contrast to the recent findings reported in Petkov et al 2012 J. Phys.: Condens. Matter 24 155102.

Liquid-liquid phase transition in supercooled yttria-alumina.

The structure and thermal characteristics of aerodynamically levitated samples of yttria-alumina in the liquid, supercooled liquid and solid phases were explored in an extensive series of high energy x-ray diffraction, small angle neutron scattering, and pyrometric cooling measurements. Particular focus was placed on the compound (Y2O3)(x)(Al2O3)(1-x) with x = 0.2 for which a liquid-liquid phase transition at a temperature of 1788 K has recently been reported. No structural or thermal signature in support of this metastable phase transition could be found.