Search NASA⌕ Search

SEARCH · Search NASA

Results for “liquid liquid equilibrium”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Liquid–Liquid Equilibrium Prediction in Fast Pyrolysis Bio-Oil Systems: A Framework for Incorporating Bio-Oil Complexity

The study of mixtures of bio-oil, water and organic solvents in different proportions can serve as a cost-effective analysis of its content due to the formation of immiscible phases. This manuscript attempts to replicate experimentally determined partition coefficients (K OW ) of relevant species present in fast pyrolysis bio-oil (FPBO). A commercial flowsheeting simulator with surrogate bio-oil model representation is used. Concurrently, pyrolytic lignins in FPBO (‘pyrolignin’) do not have an agreed-upon structural representation, and the literature is ripe with wide variations of said representations. Thus, during the description of FPBO, this pyrolignin fraction was modeled using 20 possible structures (phenolic dimers to tetramers), with the goal of determining the structures for which the experimental data are best described. Two cases were considered: Case 1 normalized the reported experimental mass balance, while Case 2 included the unreported fraction in the mass balance to the total pyroligin. Please, add here a comment on the prediction of the Water oil equilibrium. The best KOW predictions for levoglucosan (LVG) were obtained when the system was modeled with no pyrolignin, presenting an MRE under 10% for both systems WO and BO. Among the possible structures, D2 (dimer), F1 (trimer) and I1, and I3 (tetramers) presented MRE ≤ 13% for both cases.

09 BIOMASS FUELS↗
DOE OSTI2025

Bench-Scale Apparatus for Rapid, Simultaneous, Comprehensive Vapor–Liquid Equilibrium and Kinetic Property Measurements to Advance CO 2 Capture Solvent Technology

We present here, the design, fabrication, and testing of a new "pressure, volume, temperature (PVT)" apparatus that enables rapid standardized testing for viscosity, vapor–liquid equilibria, and kinetics for water-lean carbon capture solvents. This unit is the first of its kind where equilibrium data are collected during operation as a PTx (pressure–temperature–composition) cell while kinetic data are collected simultaneously with an internal mini wetted-wall contactor (WWC) using controlled adjustments of CO 2 injections to allow for measurements of gas flux (in and out of the liquid). Additionally, in situ measurements of viscosity data are also continuously collected while a solvent is in circulation during gas absorption. This cell empowers comprehensive testing of critical CO 2 capture solvent properties in a single measurement, ensuring all data are collected at the same temperature, pressure, and CO 2 loading. This apparatus also expedites screening of materials since less than 50 mL of sample is needed as compared to 2–3 L needed to get similar data from a conventional WWC. In conclusion, we describe here the methodology of data collected on this new PVT cell, nicknamed “Gary” in honor of Professor Gary Rochelle, for multiple water-lean amine solvents, which we compare to data collected from conventional instrumentation for amine testing.

Zheng, Richard [STARS Technology Corporation, Rich↗
DOE OSTI2026

Data for "Plasmon-driven exciton formation in a non-equilibrium Fermi liquid"

This repository contains source data for key plots presented in the manuscript "Plasmon-driven exciton formation in a non-equilibrium Fermi liquid." Experimental data that was analyzed in Igor Pro 8 are presented as the .pxp files used to generate individual sub-plots. Electronic spectral function calculations are provided as .txt files, in which consecutive rows refer to the meshgrid x coordinate, y coordinate, spectral function (and, where relevant, axis-projected local angular momentum). We additionally include the Wannier model and DFT-obtained bulk band structure on which the Wannier model was based. Files are named as the number of the figure in the manuscript to which they correspond, with additional details included where necessary. Details of file names: 2a_DOS_Lxz_Ek_KGM_40layer_xnum_800kpt_tot.txt: Density of states, xz-axis projected local orbital angular momentum, for 800 points along the K-Gamma-M path, for a 40-layer model. 2c_composite_y.pxp: ARPES (angle-resolved photoemission spectroscopy) spectra along the ky axis, including both a scan near the Fermi level and a scan at high kinetic energies. 2d_LCP_RCP_diff_Sect_20K.pxp: difference between ARPES constant energy cuts at T=20 K at E0 + 0.23 eV taken with left- and right-circularly polarized photons. The polarization-integrated intensity at the constant energy cut is also included. 2e_DOS_L45_E11pt79_m0pt25to0pt25_xnum_800kpt_tot.txt: Density of states, xz-projected local orbital angular momentum, and corresponding k-points in two dimensions from ab-initio electronic structure calculations for a constant-energy cut. 3a_[x]_[y]ps: ARPES cut under excitation at a fluence of x uJ/cm2, measured y ps after photoexcitation. Measurements were performed at 9 K. 3b_[x]: Energy distribution curves under excitation at a fluence x uJ/cm2 at selected delay times after photoexcitation. 4a_ImSigma_vs_temperature.pxp: Imaginary self energy (extracted from ARPES linewidths) at different energies above E0 for selected lattice temperatures. 4b_EELS_lowE.pxp: Electron energy loss spectrum over a low energy range 5b_diff_55m15.pxp: Difference between momentum-integrated Tr-ARPES traces at 55 uJ/cm2 and 15 uJ/cm2 photoexcitation. Time-dependent intensity at each energy level has been normalized to a maximum of 1 for each individual fluence prior to subtraction. 5d_invtau_at_EX_vs_fluence.pxp: decay rate at a specified energy EX for different excitation fluences, from single exponential fits. NOTE: Analyses based on the Wannier model presented here should cite both the associated Article and this dataset. For all other files in the repository, citing the dataset alone is sufficient.

Acharya, Rishi [University of Illinois] (ORCID:000↗
DOE OSTI2025

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗
DOE OSTI2024

A structural underpinning of the lower critical solution temperature (LCST) behavior behind temperature-switchable liquids

In this work, we use state-of-the-art X-ray scattering and molecular dynamics to analyze amine-water mixtures that show the unusual lower critical solution temperature (LCST) behavior. The goal is to provide direct experimental evidence for the entropy-lowering molecular cluster formation hypothesized as necessary for LCST behavior. Differential wide-angle X-ray scattering and pair distribution analysis and small-angle X-ray scattering measurements were combined with molecular modeling and liquid-liquid equilibrium measurements, revealing direct experimental evidence for the hypothesis. However, the response of the amine phase to accommodating water is even more subtle than the simple hypothesis suggests, with the formation of robust nanoscale reverse micelles. The techniques developed in this paper can be expected to yield insights in the use of temperature-switchable liquids in solvent extraction and other separations, and the stabilization of organelles in living cells that do not have physical membranes but do require compositional gradients to operate.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Deconvoluting sources of variability in aerosol jet printing using light scattering measurements

Aerosol jet printing (AJP) is a digital additive manufacturing technique for hybrid and conformal electronics, where its contactless deposition readily enables patterning over 3D surface topography. However, the complexity of aerosol transport physics makes deposition rate sensitive to process variability, hindering widespread industry adoption. Light scattering measurements have recently been demonstrated as a capable tool for measuring deposition rate in real time, enabling closed-loop control and in-situ qualification frameworks. These inline optical measurements present further opportunity as a diagnostic tool to understand the effects of secondary process parameters affecting vapor–liquid equilibrium and heat and mass transport during deposition. Downstream of the optical measurement cell, changing the sheath gas flow rate, increasing temperature via an inline heater, and adding solvent vapor via a sheath gas bubbler were observed to alter drying physics during sheath-collimation, and thereby impaction efficiency. Further upstream, the introduction of solvent vapor via a carrier gas bubbler and liquid build-up in the printhead were seen to affect the quantitative relationship between the light scattering data and the true deposition rate. Using this information, tightened controls of meaningful secondary parameters were implemented to improve the batch-to-batch consistency for printing a dielectric ink. In addition to advancing AJP process reliability for production, this work demonstrates the capability of inline optical measurements to provide insight into process physics and deconvolute competing mechanisms that underly variability.

additive manufacturing↗
DOE OSTI2025

Thermodynamic modeling of aqueous acetic acid, butyric acid, and lactic acid solutions

Based on the activity coefficient – fugacity coefficient approach, a rigorous thermodynamic modeling study is presented for accurate correlation of vapor-liquid equilibrium data of aqueous solutions of acetic acid (293 to 391 K), butyric acid (325 to 436 K), lactic acid (378 to 409 K), and acetic acid + butyric acid binary mixture (358 to 421 K). In addition, the pH data of the three aqueous, single carboxylic acid solutions were measured at 298 to 328 K and successfully correlated. Given that these aqueous carboxylic acid solutions exhibit various degrees of association behavior in both vapor and liquid phases, the thermodynamic models considered for this study include the Redlich-Kwong equation of state (RK-EoS) and the Hayden-O’Connell equation of state (HOC-EoS) for the vapor phase fugacity coefficients and the electrolyte non-random two-liquid model (eNRTL) and the association electrolyte non-random two-liquid model (AeNRTL) for the liquid phase activity coefficients. The combination of the HOC-EoS for the vapor phase and the AeNRTL model for the liquid phase is found to provide the best correlation results, consistent with the fact that the HOC-EoS and the AeNRTL model explicitly account for association behaviors in the vapor phase and the liquid phase, respectively.

09 BIOMASS FUELS↗
DOE OSTI2024

Identifying Green Solvent Mixtures for Bioproduct Separation Using Bayesian Experimental Design

Liquid–liquid extraction (LLE) is a widely used technique for the separation and purification of liquid-phase products with applications in various industries, including pharmaceuticals, petrochemicals, and renewable chemistry. A critical step in the design of an LLE process is the selection of appropriate solvents. This study presents a new methodology for identifying solvent mixtures for bioproduct separation using Bayesian experimental design (BED). Motivated by the need for environmentally friendly and effective separation methods, we address the challenge of selecting solvent systems that balance separation efficiency, selectivity, and environmental impact while also tackling the difficulty of separating multiple bioproducts using complex solvent systems. Our approach specifically seeks to predict product partition coefficients (log10 Kp values) as thermodynamic parameters underlying solvent selection. The iterative approach integrates Bayesian optimization with experimental measurements to guide solvent selection and leverages COSMO-RS simulations to enhance high-throughput experimentation. Using the design of solvent systems for the separation of lignin-derived aromatic products via centrifugal partition chromatography (CPC) as a case study, we show that within seven iterations/cycles of the methodology, we can identify new mixtures of green solvents that align with CPC design principles. Furthermore, these results demonstrate the efficacy of the BED framework in optimizing green solvent systems for complex separations, highlighting the potential of this method to advance the field of green chemistry and contribute to the development of sustainable industrial processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Optimizing CO 2 -Loaded Aqueous Amine Solutions for Higher Electrocatalytic CO 2 Reduction Activity

The activity of aqueous-based carbon dioxide reduction (CO 2 R) reactions is often limited by the solubility of CO 2 . The addition of amines can increase the total dissolved carbon in water through the formation of bicarbonate and carbamate species, which has been used to a great effect to capture CO 2 from dilute streams. Here, in this study, we explore the effect of 12 primary and secondary amines of varying Brønsted basicity, steric profile, and hydrogen-bonding capabilities on the aqueous CO 2 R to CO activity of a molecular Ni(cyclam)Cl 2 catalyst with a Hg electrode. Addition of some of the amines results in greater activity and selectivity for CO production compared to equivalent aqueous solutions without added amines. Under optimal conditions (0.4 M 3-amino-propionitrile), there is an over sevenfold increase in partial current density and greater selectivity for CO compared to equivalent conditions with no amine. Interestingly, the increase in activity did not correlate to any single property across the 12 amines. To elucidate the effect of the amine additives on catalysis, we used vapor–liquid equilibrium modeling (VLE), 13 C NMR spectroscopy, and computational analysis to determine the carbon speciation of the solutions. These results indicate that for amines without ethylalcohol functionalities, CO 2 R activity correlates with carbamate concentration, which is in turn governed by amine basicity and steric effects. However, this correlation does not persist for amines with ethylalcohol functionalities, which can form more stable carbamates through intramolecular-hydrogen bonding. These studies demonstrate that amine additives can enhance aqueous CO 2 R activity and selectivity and describe amine properties that lead to these higher performance metrics.

Amines↗
DOE OSTI2025

Efficient machine learning interatomic potentials robust for liquid and multiple solid polymorphs of NaF and KF

Achieving atomic-level understanding of crystallization of molten salts is of importance to a wide range of technological applications. Recent work [Fan et al., Proc. Natl. Acad. Sci. USA 122, e2425702122 (2025)] revealed that crystal nucleation in molten LiF salt is a multistage process according to the molecular-dynamics (MD) simulations based on an atomic cluster expansion (ACE) machine-learning interatomic potential (MLIP). In order to understand the influence of increasing cation size on nucleation pathways and nucleation rates of molten fluoride salts, here we develop two new ACE MLIPs for NaF and KF. The two ACE MLIPs feature DFT-SCAN-level accuracy for liquid and multiple solid polymorphs over a wide temperature (0–2000 K) and pressure (0–100 GPa) range, and also reproduce well a number of experimental data for solid and liquid equilibrium properties. The efficiency of the two ACE MLIPs enable million-atom-scale or microsecond-scale MD simulations. The two general-purpose ACE MLIPs are expected to be useful for atomistic simulations for different purposes, in addition to studying crystallization of molten NaF and KF salts.

Crystal melting↗
DOE OSTI2025

Liquid–Vapor Phase Equilibrium in Molten Aluminum Chloride (AlCl 3 ) Enabled by Machine Learning Interatomic Potentials

Molten salts are promising candidates in numerous clean energy applications, where knowledge of thermophysical properties and vapor pressure across their operating temperature ranges is critical for safe operations. Due to challenges in evaluating these properties using experimental methods, fast and scalable molecular simulations are essential to complement the experimental data. In this study, we developed machine learning interatomic potentials (MLIP) to study the AlCl 3 molten salt across varied thermodynamic conditions (T = 473–613 K and P = 2.7–23.4 bar), which allowed us to predict temperature-surface tension correlations and liquid–vapor phase diagram from direct simulations of two-phase coexistence in this molten salt. Two MLIP architectures, a Kernel-based potential and neural network interatomic potential (NNIP), were considered to benchmark their performance for AlCl 3 molten salt using experimental structure and density values. The NNIP potential employed in two-phase equilibrium simulations yields the critical temperature and critical density of AlCl 3 that are within 10 K (∼3%) and 0.03 g/cm 3 (∼7%) of the reported experimental values. An accurate correlation between temperature and viscosities is obtained as well. In doing so, we report that the inclusion of low-density configurations in their training is critical to more accurately represent the AlCl 3 system across a wide phase-space. The MLIP trained using PBE-D3 functional in the ab initio molecular dynamics (AIMD) simulations (120 atoms) also showed close agreement with experimentally determined molten salt structure comprising Al 2 Cl 6 dimers, as validated using Raman spectra and neutron structure factor. Furthermore, the PBE-D3 as well as its trained MLIP showed better liquid density and temperature correlation for AlCl 3 system when compared to several other density functionals explored in this work. Overall, the demonstrated approach to predict temperature correlations for liquid and vapor densities in this study can be employed to screen nuclear reactors-relevant compositions, helping to mitigate safety concerns.

Ab initio molecular dynamics↗
DOE OSTI2025

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Selection of solvents for integrated CO 2 absorption and electrochemical reduction systems

Abstract Solvent‐based electrochemical CO 2 reduction (CO 2 R) enables the production of chemicals or fuels using CO 2 from a preceding absorption process. Employing previously tested CO 2 capture solvents does not ensure their suitability for either CO 2 R or integrated CO 2 absorption‐reduction. We propose solvent selection criteria that include the CO 2 solubility, kinetic constant, ionic conductivity, concentration of the bicarbonate, carbamate, and solvent cation in the CO 2 ‐loaded solution, and sustainability indicators. They are implemented for solvent selection (a) from novel, aqueous mixtures of N ‐methylcyclohexylamine (MCA) with piperazine (PZ), 2‐amino‐2‐methyl‐1‐propanol (AMP), potassium hydroxide (KOH), and potassium chloride (KCl) and (b) from aqueous monoethanolamine (MEA), AMP, KOH, MCA, and PZ solutions. Versions of a modified Kent‐Eisenberg model for strong bases, carbamate, and non‐carbamate‐forming amine solutions are developed and parameterized through experimental equilibrium measurements. CO 2 R experimental results are presented for solutions of KOH and MCA + KOH, as these indicate desired trade‐offs for CO 2 absorption and reduction.

Amines↗
DOE OSTI2025

Machine Learning-Driven Solvent Screening for Biobased 2,3-Butanediol Extraction

Biobased 2,3-butanediol (2,3-BDO) is a valuable biomass-derived chemical due to its versatility in being transformed into a wide variety of products. However, the separation and purification of 2,3-BDO from fermentation broth remain a significant challenge owing to its high boiling point and hydrophilic nature. Herein, we developed a machine learning (ML)-based screening workflow that uses molecular calculations as training data and requires only a small number of experimental measurements for validation to identify alternative solvent candidates for the liquid–liquid extraction (LLE) of 2,3-BDO from aqueous solution. In particular, 130 density functional theory (DFT) calculations with the implicit solvation method not only built a correlation between the computational partition coefficient and the experimental distribution coefficient of 2,3-BDO but also parameterized an Extra-Trees ML model to screen the distribution coefficient for a wider range of 6717 organic solvents. The experimental measurements of only 24 solvents were needed to validate the computational results. A list of 50 prioritized solvents was proposed for 2,3-BDO LLE, and seven additional experimental measurements were conducted to further verify our selected solvents. The impact of the extraction temperature and solvent-to-feed ratio was also investigated for selected solvents in experiments. Furthermore, this work suggested alternative solvents for 2,3-BDO LLE and proposed a versatile workflow that requires fewer experiments and can be applied to a broader range of LLE studies.

Extraction↗
DOE OSTI2025

Multiscale modeling and laser diagnostics to reveal non-equilibrium reaction chemistry at a plasma-liquid interface

Low-temperature, atmospheric-pressure plasmas in contact with liquid are at the core of a wide range of applications including water treatment, medicine, materials synthesis, and chemical transformation. In general, plasma-liquid processes are scientifically compelling because reactivity can be produced without a catalyst, in relatively inert molecules, such as air or nitrogen and liquid water, in their native states at ambient conditions. However, reactions at a plasma-liquid interface are extremely complex, occurring at a multiphase, gas-liquid boundary where excited or dissociate gaseous species dissolve and react with solution-phase species, and unique reaction pathways are induced by non-equilibrium chemistry. In particular, detailed knowledge of the physical and chemical processes, including what species are produced in the gas phase and how these species are subsequently transported across and react near the interfacial region, remains largely unanswered. In this project, modeling of the non-equilibrium reaction chemistry at the interface of low-temperature, atmospheric-pressure plasmas and liquid water is developed, supported by advanced laser diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Liquid Sorption-Enhanced Haber–Bosch Process

The use of a liquid sorbent in a traditional Haber-Bosch process enables significant improvements in energy efficiency and potential cost savings for arguably the most important chemical process on the planet. The approach presented in this report employs an incompressible liquid sorbent that absorbs and releases ammonia (NH 3 ) under specific conditions. To achieve this, we investigate reactions of ammonia and pure phosphoric acid (H 3 PO 4 , PA), which rapidly neutralize to form an equilibrated solution of monoammonium phosphate (MAP) and diammonium phosphate (DAP) that functions as a reversible and regenerable sorbent. Through intimate contact of the gas-phase Haber-Bosch reaction mixture with this liquid absorbent, complete equilibrium uptake may be achieved in an appropriately sized separator, and facile separation occurs through the use of independent liquid and gas phases. Following depressurization and release of the ammonia product, only the incompressible fluid needs to be repressurized and returned to the reactor. This study documents proof-of-concept absorption and desorption experiments carried out in 75 mL batch reactors, predominantly charged with precise MAP and DAP mixtures that equilibrate at process-relevant temperatures and pressures. We then assemble the first thermodynamic relationships that underlie this advantaged separation strategy, validated by reactive force field (ReaxFF) interatomic potential simulations, and benchmarked with traditional separation routes via process modeling and technoeconomic analysis. The scale of energy consumption in the century-old Haber-Bosch process is massive, and the elegant liquid sorption approach reported here offers opportunities to enhance its energy efficiency for the next frontier of ammonia synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Structure of molten ytterbium aluminum garnet

Rare earth aluminum garnets are important materials in optical, dielectric, and thermal barrier applications. To advance the understanding of their melt processing and glass forming ability, we report the atomic structure of molten Yb 3 Al 5 O 12 over 1770–2630 K, which spans the equilibrium and supercooled liquid regimes. The melt density at T m = 2283 K is 5.50 g cm –3 , measured via silhouette imaging of electrostatically levitated drops over 1010–2420 K. Four separate structure measurements were made with aerodynamically levitated melts using x-ray and neutron diffraction with isotope substitution of Yb ( 172 Yb, 174 Yb, or nat Yb). Empirical potential structure refinement models were developed, which are in excellent agreement with the experiments. Coordination environments for Al–O are predominantly 4- and 5-coordinate, with a mean coordination of n AlO = 4.43(8), while Yb–O environments mostly range from 5- to 8-coordinate, with n YbO = 6.26(8). The cation–oxygen polyhedra are connected primarily by corner-sharing, with edge-sharing constituting up to ~1/3 of the connectivity among polyhedra with Yb or higher-coordinated Al–O. Structurally, the –Al–O– network in molten Yb 3 Al 5 O 12 appears conducive to glass formation: n OAl = 1.85(3), there are 1.86 AlO x –AlO x connections per Al atom (e.g., a mixture of Q 3 and Q 4 units), and the modal ring size is six cations. These characterize a network that is somewhat less constrained compared to SiO 2 glass, yet Yb 3 Al 5 O 12 cannot be quenched into crystal-free glass. Here, aluminum garnet compositions with larger rare earth cations do form glass, so these characterizations help reveal the structural characteristics corresponding to the limit of glass forming ability in rare earth aluminates.

36 MATERIALS SCIENCE↗