SEARCH · Search NASA
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.
A user's guide for the RSOPRNT, a TRASYS Master Restart Output Tape (RSO) reader is presented. Background information and sample runstreams, as well as, references, input requirements and options, are included.
Abstract Equilibrium phase transitions are influenced by fluctuations and often discussed within the framework of the Gibbs free energy, wherein the exchange of energy between system and thermal bath is stationary and all regions of the sample exhibit the same phase. Presence of spatial heterogeneity in the magnetic structures such as pinning centers, domain walls, topological defects, etc. may cause temporal heterogeneity that modifies the nature of the magnetic phase transition. This study reports that interplay of nanoscale thermodynamics with spatio‐temporal heterogeneity gives rise to complex phase transition pathways in amorphous Fe x Ge 1‐x thin films with temperature and Fe‐concentration ( x ). Coherent resonant soft X‐ray scattering experiments that have simultaneous spatial, temporal, and spectral sensitivity show that the origin of helical to paramagnetic phase transition in amorphous Fe‐Ge thin films lies in the appearance of enhanced‐fluctuation spots deep inside the ordered state. The fluctuations are heterogeneous, starting over a small fraction of the domains that increases and becomes isotropic over the entire film as the temperature increases or the Fe‐concentration decreases. The fluctuating‐fraction, when normalized to magnetization for different Fe‐concentrations, follows a single power law behavior, suggesting that the nature of the transition can be described in terms of the underlying spatio‐temporal fluctuations.
The MCCCS‒MN (Monte Carlo for Complex Chemical Systems‒Minnesota) software is developed by the Siepmann research group at the University of Minnesota. MCCCS‒MN allows for the simulation of multi-component molecular systems in the canonical, isobaric-isothermal (including constant stress for solids), grand-canonical, semi-grand, and Gibbs (NVT, NPT, and more than two simulation boxes) ensembles. It uses the configurational-bias Monte Carlo method to efficiently sample phase space for linear, branched and cyclic chain molecules, the adiabatic nuclear and electronic sampling Monte Carlo method to treat many-body polarization effects, and the aggregation-volume-bias Monte Carlo algorithm to efficiently sample the spatial distribution of associating molecules. MCCCS-MN employs a molecular representation of the system where force fields contain bonded and non-bonded terms. Funding for the development of MCCCS-MN through grants from the National Science Foundation (simulation of fluid phase equilibria and chromatography) and the Department of Energy (simulation of adsorption equilibria) is gratefully acknowledged.
The high-resolution measurement capability of Fourier-transform mass spectrometry (FT-MS) has made it a necessity for exploring the molecular composition of complex organic mixtures, like soil, plant, aquatic, and petroleum samples. This demand has driven a need for informatics tools to explore and analyze FT-MS data in a robust and reproducible manner. FREDA is an interactive web application developed to enable spectrometrists to format, process, and explore their FT-MS data without the need for statistical programming expertise. FREDA was built to explore outputs from a molecular identification tool, like CoreMS, and provide a suite of methods to filter data, compute chemical properties of peaks, statistically compare samples and groups of samples, conduct exploratory data analysis, and download the results with a report detailing all steps conducted. To demonstrate the utility of FREDA, an example analysis was conducted using FT-MS data from a soil microbiology study of samples collected in two different soil depths at the Sphagnum bog forest north of Grand Rapids, Minnesota. Differences between the two depths are observed using Kendrick, Gibbs free energy, and van Krevelen plots. G-tests are used to quantify a significant difference between the groups. All analyses and plotting are conducted using only the FREDA application. FREDA is an open-source and readily available web application that allows users to explore and make statistically valid conclusions about their FT-MS data. The application is available online (https://map.emsl.pnnl.gov/app/freda) with a tutorial web series (https://youtu.be/k5HLE2kNSBY?si=yB6sGoyvzxrFf5MP) and freely accessible code on Github (https://github.com/EMSL-Computing/FREDA).
For many years we have been developing the ability, through high-resolution spectroscopy targeting ro-vibrational emission in the approximately 3 - 5 micrometer region, to quantify a suite of (approximately 10) parent volatiles in comets using quantum mechanical fluorescence models. Our efforts are ongoing and our latest includes methanol (CH3OH). This is unique among traditionally targeted species in having lacked sufficiently robust models for its symmetric (v(sub 3) band) and asymmetric (v(sub 2) and v(sub 9) bands) C-H3 stretching modes, required to provide accurate predicted intensities for individual spectral lines and hence rotational temperatures and production rates. This has provided the driver for undertaking a detailed empirical study of line intensities, and has led to substantial progress regarding our ability to interpret CH3OH in comets. The present study concentrates on the spectral region from approximately 2970 - 3010 per centimeter (3.367 - 3.322 micrometer), which is dominated by emission in the (v(sub 7) band of C2H6 and the v(sub 2) band of CH3OH, with minor contributions from CH3OH (v(sub 9) band), CH4 (v(sub 3)), and OH prompt emissions (v(sub 1) and v(sub 2)- v(sub 1)). Based on laboratory jet-cooled spectra (at a rotational temperature near 20 K)[1], we incorporated approximately 100 lines of the CH3OH v(sub 2) band, having known frequencies and lower state rotational energies, into our model. Line intensities were determined through comparison with several comets we observed with NIRSPEC at Keck 2, after removal of continuum and additional molecular emissions and correcting for atmospheric extinction. In addition to the above spectral region, NIRSPEC allows simultaneous sampling of the CH3OH v(sub 3) band (centered at 2844 per centimeter, or 3.516 micrometers and several hot bands of H2O in the approximately 2.85 - 2.9 micrometer region, at a nominal spectral resolving power of approximately 25,000 [2]. Empirical g-factors for v(sub 2) lines were based on the production rate as determined from the v(sub 3) Q-branch intensity; application to comets spanning a range of rotational temperatures (approximately 50 - 90 K) will be reported. This work represents an extension of that presented for comet 21P/Giacobini-Zinner at the 2010 Division for Planetary Sciences meeting [3]. Our empirical study also allows for quantifying CH3OH in comets using IR spectrometers for which the v(sub 3) and v(sub 2) bands are not sampled simultaneously, for example CSHELL/NASA IRTF or CRIRES/VLT.
Abstract Ternary compounds with an immiscible pair of elements are relatively unexplored but promising for novel quantum materials discovery. Exploring what third element and its ratio that can be added to make stable ternary compounds out of an immiscible pair of elements remains a great challenge. In this work, we combine a machine learning (ML) method with ab initio calculations to efficiently search for the energetically favorable ternary La-Co-Pb compounds containing immiscible elements Co and Pb. Three previously reported structures are correctly captured by our approach. Moreover, we predict a ground state La 3 CoPb compound and 57 low-energy La-Co-Pb ternary compounds. Attempts to synthesize La 3 CoPb via multiple techniques produce mixed or multi-phases samples with, at best, ambiguous signals of the predicted lowest-energy La 3 CoPb and the second lowest-energy La 18 Co 28 Pb 3 phases. The calculated results of Gibbs free energy are consistent with experiments, and will provide very useful guidance for further experimental synthesis.
Here we introduce a new training algorithm for deep neural networks that utilize random complex exponential activation functions. Our approach employs a Markov chain Monte Carlo sampling procedure to iteratively train network layers, avoiding global and gradient-based optimization while maintaining error control. It consistently attains the theoretical approximation rate for residual networks with complex exponential activation functions, determined by network complexity. Additionally, it enables efficient learning of multiscale and high-frequency features, producing interpretable parameter distributions. Despite using sinusoidal basis functions, we do not observe Gibbs phenomena in approximating discontinuous target functions.
On 2014 May 22 and 24 we characterized the volatile composition of the dynamically new Oort cloud comet C2012 K1 (PanSTARRS) using the long-slit, high resolution ( lambda/delta lambda is approximately or equal to 25,000) near-infrared echelle spectrograph (NIRSPEC) at the 10 m Keck II telescope on Maunakea, Hawaii. We detected fluorescent emission from six primary volatiles (H2O, HCN, CH4, C2H6, CH3OH, and CO). Upper limits were derived for C2H2, NH3, and H2CO. We report rotational temperatures, production rates, and mixing ratios (relative to water). Compared with median abundance ratios for primary volatiles in other sampled Oort cloud comets, trace gas abundance ratios in C2012 K1 (PanSTARRS) for CO and HCN are consistent, but CH3OH and C2H6 are enriched while H2CO, CH4, and possibly C2H2 are depleted. When placed in context with comets observed in the near- infrared to date, the data suggest a continuous distribution of abundances of some organic volatiles (HCN, C2H6, CH3OH, CH4) among the comet population. The level of enrichment or depletion in a given comet does not necessarily correlate across all molecules sampled, suggesting that chemical diversity among comets may be more complex than the simple organics-enriched, organics-normal, and organics-depleted framework.
The first-reported competitive CO 2 –H 2 O sorption isotherms on a naturally occurring expandable clay mineral (saponite) exchanged with Na, K, and Cs at near ambient CO 2 concentrations (1% CO 2 ) are reported by combining ATR-FTIR spectroscopy with gravimetric sorption methods. Little CO 2 sorption was observed at high relative humidity (RH) but increased significantly at low RH. Unlike prior work where CO 2 sorption was reported on samples that were previously subjected to heating and/or evacuation treatment to remove sorbed H 2 O, this study manipulated CO 2 sorption by controlling RH. CO 2 sorption increased in the order Cs- > K- > Na-exchanged saponite and was anticorrelated with the Gibbs energy of hydration of the exchangeable cation. As expected, H 2 O sorption followed the opposite trend, with increased H 2 O sorption occurring in the order Na > K > Cs. The amount of CO 2 sorbed using 1 atm CO 2 (p/p o = 1.0) ranged from 50 to 30 mg CO 2 /g clay and from 1.25 to 0.8 mg CO2 /g clay for 1% CO 2 (p/p o = 0.01) in N 2 . The sensitivity of the ATR-FTIR measurements is demonstrated by detecting the CO 2 stretch of adsorbed CO 2 from air containing 500 ppmv CO 2 . The spectral features of sorbed CO 2 were minimally affected by the nature of the exchangeable cation or variations in H 2 O content. However, the position of the HOH bending band shifted significantly from 1638 to 1610 cm –1 as H 2 O content decreased, reflecting changes in intermolecular hydrogen bonding between H 2 O molecules in the saponite interlayer. The shift in position of the HOH bending band to lower energy was coincident with increased CO 2 sorption. CO 2 sorption in clay mineral interlayers is consistent with sorption on weakly hydrated partially hydrophobic sites found on low charge density smectites and is enhanced by lower RH and by the presence of weakly hydrated exchangeable cations like Cs + or K + . In conclusion, the greatest CO 2 sorption occurs on Cs-saponite because this is the least hydrated cation of the three used in the study and indicates that CO2 sorption occurs on the neutral portion of the siloxane surface.
This data package is associated with the publication “Moisture content modulates DOM thermodynamic regulation of oxygen consumption in drying streambed sediments” published in Scientific Reports (Garayburu-Caruso et al., 2026). The package contains processed data products and scripts used to quantify how drying and re-inundation of riverbed sediments influence dissolved organic matter (DOM) thermodynamic properties and their relationship with sediment oxygen (O₂) consumption across 33 stream sites in the contiguous United States. The data package contains DOM thermodynamic metrics (e.g., Gibbs free energy of carbon oxidation and thermodynamic efficiency), and O₂ consumption along with watershed-scale climate and land-cover metrics used as explanatory variables in the analyses. Underlying unprocessed and processed ultrahigh-resolution mass spectrometry data, oxygen consumption rates from laboratory moisture-manipulation experiments, within-sample environmental properties, sediment moisture content and contextual field measurements are archived separately at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/2428003 (Laan et al., 2024) and https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1923689 (Forbes et al.,2023). A preliminary version of this data package was published in February 2026 at the time of manuscript submission. It was updated in June 2026, at the time of manuscript acceptance, to include the finalized data and additional metadata (readme, data dictionary, and file level metadata). For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. At the top level, the data package is organized into five main folders: (1) Data, (2)Figures, (3) Map, (4) GAM_Reulsts, and (5) src. The Data folder contains analysis-ready tabular files with oxygen consumption rates, DOM thermodynamic properties by site and treatment, site-level environmental variables, watershed-scale metrics, and other derived variables referenced in the manuscript. The Figures folder contains static image files associated with the main text and supplemental figures, while the Map folder includes spatial data and map-layer files used to create the sampling-location map. The GAM results folder contains the results for each of the general additive model (GAM).The src folder contains R scripts used to perform data processing, statistical analyses (including clustering, generalized additive models, and threshold analysis), and figure generation. This data package is associated with a GitHub repository found at https://github.com/WHONDRS-Hub/ECA_DOM_Thermodynamics.
Here, we present a data-centric deep learning (DL) approach using neural networks (NNs) to predict the thermodynamics of ternary solid solutions. We explore how NNs can be trained with a dataset of Gibbs free energies computed from a CALPHAD database to predict ternary systems as a function of composition and temperature. We have chosen the energetics of the FCC solid solution phase in 226 binaries consisting of 23 elements at 11 different temperatures to demonstrate the feasibility. The number of binary data points included in the present study is 102,000. We select six ternaries to augment the binary dataset to investigate their influence on the NN prediction accuracy. We examine the sensitivity of data sampling on the prediction accuracy of NNs over selected ternary systems. It is anticipated that the current DL workflow can be further elevated by integrating advanced descriptors beyond the elemental composition and more curated training datasets to improve prediction accuracy and applicability.
Introduction: HED meteorites are thought to rep-resent igneous rocks from Vesta’s basaltic crust and preserve evidence of early crustal metamorphism. Determining the temperatures and timescales of thermal metamorphism is important for reconstructing crustal evolution in the early solar system. Here, we study basaltic eucrite Graves Nunataks (GRA) 98098, which has been identified as a highly metamorphosed eucrite [1]. We present new estimates on metamorphic temperatures determined via thermodynamic modeling as well as the initial results from diffusion models constraining timescales of thermal metamorphism. Sample Description: GRA 98098 is an unbrecciated eucrite with a granoblastic plagioclase and pyroxene mineralogy. Millimeter to cm-long lathes of tridymite cross-cut and poikilitically enclose plagioclase and pyroxene [1,this work]. Pyroxene grains have exsolved into Ca-rich (~Wo38En29Fs33) and Ca-poor (~Wo4.5En36Fs59.5) lamellae. Both unzoned and zoned plagioclase grains are observed. Unzoned plagioclase grains are found solely with tridymite laths. These grains have ~An92 compositions. The cores of the zoned plagioclase grains have the same composition and thin, relatively sodic rims (~An67), (Fig. 1). The bulk sample is unusually enriched in highly in-compatible elements and has one of the most fractionated REE patterns reported [1]. Maximum metamorphic temperatures of 985±78°C have been estimated using two-pyroxene thermometry [2]. Methods: Thermodynamic modeling. Thermodynamic models were constructed using the software Perple_X, which employs a Gibbs free energy minimization in order to determine the most stable phase assemblage for a given bulk rock composition [3]. Bulk composition was calculated using mineral com-positions acquired via EMPA (this study) and the observed abundancies present in the thin section. Two bulk compositions were estimated; 1) includes all phases present in the thin section, (assumes that all phases are present during metamorphism), 2) excludes tridymite from the bulk calculation (assumes that tridymite was not present during metamorphism). In order to determine whether metamorphic equilibria was achieved and estimate temperatures of metamorphism, we compared measured pyroxene compositions with thermodynamically predicted compositions [4, Fig. 2]. Diffusion Modeling. Several time-temperature de-pendent diffusion profiles were calculated in order to determine the best match for XAn chemical profiles observed at the edges of the zoned plagioclase (Fig. 3). We assumed that the start condition was a stepwise gradient at the plagioclase/pyroxene interface. We also assumed an average diffusion coefficient (D) and a constant temperature using the equation in [5]. D was determined for two temperatures (T = 1060ºC; near eucrite solidus [6] and T = 985ºC; metamorphism reported in [4]) and then XAn was calculated as a function of distance from plagioclase core to rim using an error function solution to Fick’s second law. Results: Thermodynamic model results are summarized in Fig. 2. For a bulk composition that includes all phases in the thin section, pyroxene endmember compositions plot in the following temperature ranges: Fs ~660-860ºC, En~1000ºC & 1150ºC, and Wo~760-900ºC (Fig. 2a). For a bulk composition that excludes tridymite from the peak metamorphic assemblage (i.e., the bulk composition minus the contribution from tridymite), a temperature range could not be determined for the Fs component of pyroxene. For Wo, T~760-900ºC and En, T~1000ºC & 1150ºC (Fig. 2b). Fig. 3 summarizes the diffusion model results. For T = 1060°C & 985°C, the most appropriate time interval was estimated based on which diffusion curve most matched (solid lines, Fig. 3) the EMPA data. For T = 1060°C, the best looking match was t = 500 ka. For T = 985°C, the best match was t = 7 Ma. Discussion and future work: Temperature estimates from thermodynamic models are not conclusive because the temperature ranges determined for pyroxene endmember stability do not overlap (colored fields in Fig. 2), thus implying that there is disequilibrium between pyroxene crystals and the bulk composition considered [4]. Thus, additional exploration is needed to define a metamorphically equilibrated do-main that accurately records peak temperature. The utility of defining metamorphically equilibrated do-mains to improve the accuracy and level of detail elucidated regarding the petrogenetic history of metamorphose samples has been demonstrated previously [4,7]. We suggest that in the case of Fig. 2a, the thin section composition is not representative of the length scales over which metamorphic equilibrium was achieve and in the case of Fig. 2b, the assumption that tridymite was not present during metamorphism was incorrect. However, results from thermodynamic models can provide insight into the relative timing of mineral and compositional textures. For example from texture alone, it is unclear whether tridymite was igneous in origin and represents the last bits of melt in a crystallizing magma chamber, or if it formed during (and possibly initiated) open system thermal metamorphism. The latter could be consistent with a partial melt hypothesis [8,9] while the former implies that simple fractional crystallization can yield the textures present in GRA 98098. The lack of coincidence be-tween pyroxene endmember compositions in Fig. 2b suggest that the bulk composition minus tridymite was not the assemblage in equilibrium with the pyroxene, suggesting that tridymite was present during metamorphism and formed during igneous crystallization. We conclude that the development of the Na-rich plagioclase rims likely occurred during or immediately after peak thermal metamorphism, because eucrites of similar metamorphic grade and texture have unzoned plagioclase (~An92) [2,4,8], and Na zoning is only observed in the plagioclase not included in the tridymite. This suggests that the zoning formed after tridymite formation, and therefore after igneous crystallization. Thus, the timescales calculated via diffusion modeling possibly represent the time interval over which thermal metamorphism occurred. Cooling rates approximated for the Vestan crust predict that the crust cooled below 300°C around 35-40 Ma after formation[10]. This is consistent with our modeling results that predict formation of the Na rich plagioclase rims occurring at higher temperatures over a period of 0.5 to 7 Ma years. Future work. Additional thermodynamic modeling work will focus on selecting an equilibrated bulk rock domain in which to elucidate metamorphic conditions. Diffusion models currently provide a minimum time-scale, since diffusion slows down as the system cools. Future work with will focus on integrating cooling into the diffusion models and constraining the depth at which thermal metamorphism occurs because it could be used to determine whether the range of time-scales calculated for thermal metamorphism are consistent with the geologic environment.
Metabolic potential and activity associated with hydrogen and carbon monoxide were characterized in fluids sampled from the the Coast Range Ophiolite Microbial Observatory (CROMO). CROMO consists of two clusters of science-dedicated wells drilled to varying depths up to 35m in the actively serpentinizing, Jurassic-age Coast Range Ophiolite of Northern California, along with a suite of pre-existing monitoring wells at the same site. Consistent with the fluid chemistry observed in other serpentinizing systems, CROMO fluids are highly alkaline, with pH up to 12.5, high in methane, with concentrations up 1600 micromolar, and low in dissolved inorganic carbon (DIC), with concentrations of 10's to 100's of micromolar. CROMO is conspicuous for fluid H2 concentrations that are consistently sub-micromolar, orders of magnitude lower than is typical of other systems. However, higher H2 concentrations (10's -100's of micromolar) at an earlier stage of fluid chemical evolution are predicted by, or consistent with: thermodynamic models for fluid chemistry based on parent rock composition equivalent to local peridotite and with water:rock ratio constrained by observed pH; the presence of magnetite at several wt% in CROMO drill cores; and concentrations of formate and carbon monoxide that would require elevated H2 if formed in equilibrium with H2 and DIC. Calculated Gibbs energy changes for reaction of H2 and CO in each of several metabolisms, across the range of fluid composition encompassed by the CROMO wells, range from bioenergetically feasible (capable of driving ATP synthesis) to thermodynamically unfavorable. Active consumption relative to killed controls was observed for both CO and H2 during incubation of fluids from the pre-existing monitoring wells; in incubations of freshly cored solids, consumption was only observed in one sample set (corresponding to the lowest pH) out of three. The specific metabolisms by which H2 and CO are consumed remain to be determined.
In this work, the heat capacities of nanocrystalline Ni 3 Fe and control materials with larger crystallites were measured from 0.4–300 K. The heat capacities were integrated to obtain the enthalpy, entropy, and Gibbs free energy and to quantify how these thermodynamic functions are altered by nanocrystallinity. From the phonon density of states (DOS) measured by inelastic neutron scattering, we find that the Gibbs free energy is dominated by phonons and that the larger heat capacity of the nanomaterial below 100 K is attributable to its enhanced phonon DOS at low energies. Besides electronic and magnetic contributions, the nanocrystalline material has an additional contribution at higher temperatures, consistent with phonon anharmonicity. The nanocrystalline material shows a stronger increase with temperature of both the enthalpy and entropy compared to the bulk sample. Its entropy exceeds that of the bulk material by 0.4 $k_B$/atom at 300 K. This is insufficient to overcome the enthalpy of grain boundaries and defects in the nanocrystalline material, making it thermodynamically unstable with respect to the bulk control material.
Abstract. Ultrafine aerosol particles with sizes smaller than 50 nm have been shown in recent studies to serve as a large source of cloud condensation nuclei that can promote additional cloud droplet formation under supersaturation conditions. Knowledge of the microphysics of liquid water in these droplets remains limited, particularly in the sub-10 nm particle size range, due to experimental and theoretical challenges associated with the complexity of aerosol components and the small length scales of interest (e.g., difficulty of precisely sampling the liquid–air interface, questionable validity of mean-field theoretical representations). Here, we carried out molecular dynamics (MD) simulations of aerosol particles with diameters between 1 and 10 nm and characterized atomistic-level structure and water dynamics in well-mixed and phase-separated systems with different particle sizes, NaCl salinities, and pimelic acid (PML) organic surface loadings as a function of distance from the time-averaged Gibbs dividing interface or instantaneous water–air interface. We define a sphericity factor (ϕ) that can shed light on the phase-mixing state of nanodroplets, and we reveal an unexpected dependence of mixing state on droplet size. Our results also evidence an ion concentration enhancement in ultrafine aerosols, which should modulate salt nucleation kinetics in sub-10 nm droplets, and provide detailed characterization of the influence of droplet size on surface tension and on water self-diffusivity near the interface. Analysis of water evaporation free energy and water activity demonstrates the validity of the Kelvin equation and Köhler theory at droplet sizes larger than 4 nm under moderate salinities and organic loadings and the need for further extension to account for ion concentration enhancement in sub-10 nm aerosols, droplet-size-dependent phase separation effects, and a sharp decrease in the cohesiveness of liquid water in sub-4 nm droplets. Finally, we show that an idealized fractional surface coating factor (fs) can be used to categorize and reconcile water accommodation coefficients (α*) observed in MD simulations and experimental results in the presence of organic coatings, and we resolve the droplet size dependence of α*.
We report high entropy alloys (HEAs) have gained interest for structural applications in extreme environments. With a potentially vast chemical and phase space, there are significant opportunities to discover superior performing alloys. Crucial for most high-temperature applications is understanding and mitigating the oxidation behavior of these chemically complex alloys. Most experimental and computational HEA studies have focused on a limited set of compositions and only a fraction these compositions have been characterized for oxidation. We present a high-throughput framework that utilizes density-functional theory (DFT) in concert with a combined machine-learning model and grand-canonical linear programming for assessing phase stability, phase-fraction, chemical activity and high-temperature survivability of arbitrary HEAs. This framework considers temperature dependent contributions to the Gibbs energy of the competing phases arising from short-range order and vibrational entropy. We demonstrate the effectiveness of the framework by assessing the thermodynamic stability, oxidation behavior, chemical activity, and phase decomposition of body-centered cubic Mo-W-Ta-Ti-Zr refractory HEAs. A total of 51 compositions were analyzed and ranked in order of their survivability based on the Pareto-front analysis. Oxidation was performed at 1373 K on fours samples in air showing the difference in oxidation behavior determined experimentally through scale thickness and their mass changes. The insights on oxidation behavior presented in this work will enable the fast assessment of technologically useful HEAs needed for future structural application in extreme conditions
Basaltic eucrites, which are possibly parts of the Vestan crust, formed as either lava flows or intrusions, and can offer insights into early crust formation. However most eucrites have experienced thermal metamorphism, which exacerbates the challenges of understanding these samples. A variety of petrogenetic models such as, partial melting of a primitive source, fractional crystallization of magmas emplaced in a primitive crust, and/or partial melting of a eucritic source coupled with melt mixing and assimilation have been proposed in order to explain observed petrologic and chemical characteristics of eucrites. EET 90020 is an unbrecciated eucrite that has experienced significant thermal metamorphism, with temperatures of metamorphic equilibration calculated to be ~840-1040°C. Its petrologic history remains contentious, in part due to differences in trace element analyses [10,17, this work], which has led to multiple metamorphic interpretations [8,10-11]. Here, we combined petrologic observations, chemical analyses, and thermodynamic modeling, to interpret micro-domain textures identified in EET 90020 and better constrain its petrologic history. Ultimately, the development of such textures is a direct result of the geologic processes operating on a young Vesta or similar asteroid. Sample description: We identified three textural domains in EET 90020 (Figure 1); (I) A coarse grain domain with granoblastic plagioclase and pyroxene, (II) a medium grain domain with curved grain boundaries between plagioclase and pyroxene, and (III) a fine grain domain dominated by tridymite that fills in interstitial space between spherical plagioclase. Major element chemistry in phases is consistent across domains. Pyroxene have distinctive high Ca lamellae (Wo40.92En28.1Fs31.0) forming from the low-Ca pigeonite host (Wo3.5En33.7Fs62.8). There are small amounts of fayalitic olivine along pyroxene/oxide grain boundaries, where pyroxenes have reacted with ilmenite/chromite clasts. Plagioclase is anorthitic (Avg ~An88), with little variation (An86-92.5). Trace element bulk rock analyses were collected via ICP-MS and individual mineral analyses with LA-ICP-MS (Figure 2). The bulk sample has an observable enrichment in light rare earth elements, with a notable depletion in Eu. Plagioclase is enriched in LREEs, depleted in HREEs, and has a positive Eu anomaly while pyroxene is enriched in HREEs with a negative Eu anomaly. Variations in trace element analyses is likely influenced by the presence of phosphates because are highly enriched in REEs and have a dramatic effect on the bulk domain compositions. Thermodynamic modeling: The software package, Perple_X, was used to calculate mineral phase equilibria over a range of conditions using a Gibbs free energy minimization approach [12]. Models used all major/minor elements except P2O5. Thermodynamic properties from [13] were used to constrain endmember phase stabilities. Activity models were used to describe mixing in phases with solid solution. Isochemical P-T phase diagrams were constructed for the bulk thin section, and the coarse, medium and fine grain domain compositions. Domain compositions were determined using microprobe analyses of phases identified in each domain and observed modal mineral abundancies. Thermodynamically stable pyroxene compositions were extracted from model results and compared to pyroxene compositions collected via EMPA, allowing us to calculate temperatures for metamorphic equilibrium. Results: Results that use the coarse domain composition (Fig. 3a) demonstrate overlap at T ~ 1020°C indicating high T equilibration. There was no overlap for the medium and fine grain domains, or the bulk composition (e.g., Fig. 3b) indicating that these are not equilibrium assemblages. These results were replicated using data from [8,20]. Discussion: The disequilibrium results from the bulk rock, fine grain and medium grain models imply that EET 90020 was modified during and/or after peak metamorphism (open system behavior). However, this is in contrast to the coarse domain where it appears as if equilibrium was maintained. We suggest that this conflict is explained by highly localized equilibrium occurring at millimeter length-scales, which has been observed in terrestrial metamorphism where aqueous fluid infiltrate [14] or melt loss occurs [15]. It is unlikely that aqueous alteration was significant in modifying the bulk composition of EET 90020, because it lacks hydrated minerals. However, the presence of melt is supported by the presence of an abundance of curved grain boundaries and spherical grains in the medium and fine grain domains This contrasts the coarse domain, where grain boundaries are angular. Given the evidence for disequilibrium in the medium and fine grain domains, it would be inappropriate to draw additional conclusions about textural evolution using thermodynamic models of those domains. However, modeling results of the coarse grain domain, where equilibrium was maintained, can provide useful insights into the geologic evolution of EET 90020. We approximate maximum metamorphic temperatures, using the Ca component in pyroxene, at ~1020°C (Fig. 4), which is consistent with previous two-pyroxene thermometry [11]. At this temperature, up to 10 vol % melt can be produced in the coarse grain domain, which is at the boundary of minimum melt needed in order to segregate and form a melt network [16]. Tridymite and plagioclase would be the first minerals to melt out during heating starting ~990°C. This is consistent with the observed melt textures formed by tridymite and plagioclase in the fine grain domain. We suggest that a partial melting model best explains the textural development of EET 90020. Upon heating, melt generated was either trapped, resulting in no net change in the localized rock composition, or migrated, resulting in disequilibrium between pyroxene and the surrounding matrix. A positive Eu anomaly in plagioclase indicates that all of the plagioclase was melted previously. This is consistent with previous studies that have suggested that EET 90020 represents a residual eucrite that experienced partial melting and subsequent melt loss [8]. We speculate that partial melting could have occurred when EET 90020 was heated in the asteroid’s lower crust or adjacent to an intrusion [7]. Future work will study other, texturally heterogeneous samples to see if EET 90020 presents a unique circumstance or whether crustal partial melting was ubiquitous on the eucrite asteroid.
Here, a set of self-consistent Gibbs energy functions has been developed to describe the thermochemical behavior of the major reactive fission product ZrCl 4 with key components in chloride molten salt reactors (MSRs): LiCl–ZrCl 4 , NaCl–ZrCl 4 , KCl–ZrCl 4 , CsCl–ZrCl 4 , MgCl 2 –ZrCl 4 and UCl 3 –ZrCl 4 . Low ZrCl 4 concentration phase equilibria most relevant to MSR applications have been confirmed via differential scanning calorimetry for NaCl–ZrCl 4 , and X-ray diffraction analysis performed on equilibrated samples of NaCl–ZrCl 4 , KCl–ZrCl 4 , MgCl 2 –ZrCl 4 , and UCl 3 –ZrCl 4 . Within the framework of the modified quasi-chemical model in the quadruplet approximation, extrapolations of pseudo-binary models were generated to represent KCl–MgCl 2 –ZrCl 4 , KCl–NaCl–ZrCl 4 , and MgCl 2 –NaCl–ZrCl 4 , which show agreement with phase equilibria data where available without the use of ternary interaction parameters. The optimized thermodynamic descriptions for these systems and others are available in the open-source compendium Molten Salt Thermal Properties Database – Thermochemical (MSTDB–TC).