Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PHASE DIAGRAMS”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Using Neural Network Force Fields to Ascertain the Quality of Ab Initio Simulations of Liquid Water

Accurately simulating the properties of bulk water, despite the apparent simplicity of the molecule, is still a challenge. In order to fully understand and reproduce its complex phase diagram, it is necessary to perform simulations at the ab initio level, including quantum mechanical effects both for electrons and nuclei. This comes at a high computational cost, given that the structural and dynamical properties tend to require long timescales and large simulation cells. In this work, we evaluate the errors that density functional theory (DFT)-based simulations routinely incur into due time- and size-scale limitations. These errors are evaluated using neural-network-trained force fields that are accurate at the level of DFT methods. We compare different exchange and correlation potentials for properties of bulk water that require large timescales. Here we show that structural properties are less dependent on the system size and that dynamical properties such as the diffusion coefficient have a strong dependence on the simulation size and timescale. Our results facilitate comparisons of DFT-based simulation results with experiments and offer a path to discriminate between model and convergence errors in these simulations.

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

First-Principles Insights into the Thermocatalytic Cracking of Ammonia-Hydrogen Blends on Fe(110): 1. Thermodynamics

Ammonia (NH 3 ) is being considered as a practical means to transport hydrogen (H 2 ) because of its higher volumetric energy density for the same temperature and pressure. Thermodynamics suggest high temperature is needed to decompose NH 3 to nitrogen (N 2 ) and H 2 . Furthermore, overcoming decomposition kinetic barriers requires a catalyst. Via density functional theory, we study this reaction on a model catalyst: the close-packed (110) facet of α-Fe. Specifically, we predict detailed in-operando temperature- and pressure-dependent surface phase diagrams on this benchmark catalyst that offer insights for the design of optimal NH 3 decomposition catalysts. Here, we explore the equilibrium composition(s) of the Fe(110) surface when exposed to NH 3 –H 2 mixtures. Further, we predict that both N and NH partially cover the Fe(110) surface at 300–400 °C (far above the NH 3 decomposition-formation coexistence temperature at standard partial pressures of 1 bar: ~ 180 °C) and 2–4 bar of total reactor pressure. At the equilibrium N/NH coverage, these species inhibit coadsorption of H, indicating that direct H 2 production may occur. However, from thermodynamics alone, removal of N/NH as N 2 (g) is extremely unfavorable even at these elevated temperatures–effectively deactivating the surface toward further NH 3 decomposition. Thus, catalysts that remain active toward N–H bond breaking but with weakened N binding relative to Fe are needed to enable high-turnover catalytic NH 3 decomposition to release H 2 .

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

Determining Catalytically Relevant Surfaces through Coverage-Dependent Lattice Gas Models: Carbon Adsorption on Fe(100)

Here, we have quantified the C–C lateral interactions on Fe(100) using a density functional theory (DFT)-parameterized lattice gas cluster expansion (LG CE) model trained using 265 unique configurations spanning a C coverage from 0 to 1 monolayer (ML). Our LG CE model shows high predictive accuracy with a leave-multiple-out cross-validation score of 10.2 and 16.6 meV/site for systems with and without the top two layers of Fe atoms fixed, respectively. Electronic ground-state structures identified from the lattice gas model (including the structures at 0 and 1 monolayers) were further used to generate ab initio phase diagrams under a range of temperatures and pressures. At low temperatures (<400 K), we found that the 1.0 monolayer structure is dominant, whereas at higher temperatures (>500 K), the 0.88 ML structure is most likely to form on the Fe surface. Interestingly, our model identified a c (2 × 2) ordered structure at 1/2 ML, which correlates well with previous DFT studies for carbon adsorption on iron surfaces and matches with the experimentally observed low-energy electron diffraction structure. Overall, the DFT-parameterized energies for the C/Fe system including effects of coverage and configurational space can further help in developing multiscale models for various heterogeneous reactions involving C–C and C–Fe interactions.

08 HYDROGEN↗
DOE OSTI2023

Ab Initio Molecular Dynamics Spectra for Characterization of Hydrated Al 2 O 3 Supported MoO x

The structure of supported MoO x /Al 2 O 3 catalysts is investigated using ab initio molecular dynamics (AIMD) and density-functional theory (DFT). Phase diagrams were computed to understand the hydroxylation coverage as a function of the temperature, H 2 O partial pressure, and MoO x loading. We relate the shifts in experimental Raman vibrational frequencies under hydrated conditions to hydrogen bonding interactions and the MoO x anchoring location. We showcase that the use of AIMD as a benchmark for DFT can provide insight into how, under certain hydrated conditions, DFT excels as an inexpensive method for computing vibrational frequencies, while, under dehydrated conditions, it is susceptible to the largest errors. Additionally, to facilitate the analysis of the hydroxyl region of experimental infrared spectra, we compute the power spectra of individual hydroxyls and see a strong relationship between OH stretching frequencies and the surrounding coordination environment, which are also impacted by anchored MoO x . Here, we compare computational and experimental IR and Raman spectra of catalysts synthesized herein under the same conditions.

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

Coverage-Dependent Adsorption of Hydrogen on Fe(100): Determining Catalytically Relevant Surface Structures via Lattice Gas Models

Hydrogen adatoms are a critical surface species for several reactions catalyzed by Fe surfaces such as Fischer–Tropsch synthesis, ammonia synthesis, and the hydrodeoxygenation of biomass-derived oxygenates. Parameterizing the energetics for H/Fe in terms of both coverage and configuration space can significantly aid in the development of multiscale models as well as provide atomic level insight into the dominant surface structures present under realistic reaction conditions. Here, we construct a lattice gas model for H/Fe(100), where the lateral interactions are determined from first-principles using density functional theory. Using 950 symmetrically unique H/Fe(100) configurations, we generate a cluster expansion with a predictive accuracy in terms of surface energy of 3.8 meV/site over a coverage range from 0 to 3 monolayers. Ten electronic ground state structures are identified from this thorough scan (including the structures at 0 and 3 monolayers), which were subsequently used to generate ab initio phase diagrams under a range of temperatures and pressures. Under reaction conditions typical of Fischer–Tropsch synthesis, ammonia synthesis, and biomass oxygenate hydrodeoxygenation, we find that the 1.0 monolayer structure is dominant. Furthermore, examination of the total H–H lateral interactions for the H/Fe(100) electronic ground state structures shows that H/Fe(100) can be accurately modeled via a mean-field ideal lattice gas model for coverages less than 1.0 monolayers. Altogether, this work enables the incorporation of H–H lateral interactions on Fe(100) into multiscale models, via either mean-field or site-dependent techniques, and provides atomic insight into the catalytically relevant H/Fe(100) structures for a range of heterogeneous reactions.

09 BIOMASS FUELS↗
DOE OSTI2020

Acetonitrile Transition Metal Interfaces from First Principles

Acetonitrile is among the most commonly used nonaqueous solvents in catalysis and electrochemistry. Here, we study its interfaces with multiple facets of the metals Ag, Cu, Pt, and Rh using density functional theory calculations; the structures reported shed new light on experimental observations and underscore the importance of solvent–solvent interactions at high coverage. We investigate the relationship of potential of zero charge (PZC) to metal work function, reporting results in agreement with experimental measurements. We develop a model to explain the effects of solvent chemisorption and orientation on the PZC to within a mean absolute deviation of 0.08–0.12 V for all facets studied. Our electrostatic field dependent phase diagram agrees with spectroscopic observations and sheds new light on electrostatic field effects. This work provides new insight into experimental observations on acetonitrile metal interfaces and provides guidance for future studies of acetonitrile and other nonaqueous solvent interfaces with transition metals.

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

Anomalies and Local Structure of Liquid Water from Boiling to the Supercooled Regime as Predicted by the Many-Body MB-pol Model

For the past 50 years, researchers have sought molecular models that can accurately reproduce water’s microscopic structure and thermophysical properties across broad ranges of its complex phase diagram. Herein, molecular dynamics simulations with the many-body MB-pol model are performed to monitor the thermodynamic response functions and local structure of liquid water from the boiling point down to deeply supercooled temperatures at ambient pressure. The isothermal compressibility and isobaric heat capacity show maxima near 223 K, in excellent agreement with recent experiments, and the liquid density exhibits a minimum at ~208 K. A local tetrahedral arrangement, where each water molecule accepts and donates two hydrogen bonds, is found to be the most probable hydrogen-bonding topology at all temperatures. Further, this work suggests that MB-pol may provide predictive capability for studies of liquid water’s physical properties across broad ranges of thermodynamic states, including the so-called water’s “no man’s land” which is difficult to probe experimentally.

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

Topological Frustration as a New Parameter to Tune Morphology Revealed through Exploring the Continuum between A-B-C 3-Arm Star and Linear Triblock Polymers

Block polymers assemble into a variety of phase-separated morphologies based on volume fraction (φ) and interactions (χ) of the respective blocks. The arrangement of three different polymer blocks could either be a 3-arm star, with each block having one terminus attached to a common junction point or a linear A-B-C architecture. A versatile strategy is reported to synthesize a series of well-defined graft polymers that lie along the unexplored continuum between a 3-arm star and an A-B-C linear triblock polymer architecture. Using the technique of single-molecule insertion, precise control over the position of graft arm C along the B chain was achieved. A series of discrete graft polymers (PMMA-b-PS-g-PEO) with fixed φ and prescribed ω values that lie on the continuum between a 3-arm star (ω = 0) and linear triblock polymer (ω = 1) were synthesized. Morphological studies using small-angle X-ray scattering and conventional and energy-filtered transmission electron microscopy reveal the transition between lamellae, perforated lamellae, and cylindrical morphologies with systematic variation in the ω values, a trend attributed to the topological frustration and the associated χ values between the three blocks. Molecular dynamic simulations of coarse-grained models were found to predict phase diagrams that are consistent with the experimentally observed morphologies. Furthermore, our results suggest that changes in ω lead to topological frustration which is an important additional new design parameter that can be used to tune the morphology of multiblock polymers in addition to φ and χ.

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

Correction to “Frank–Kasper Phases of Diblock Copolymer Melts Studied with the DPD Model: SCF Results”

In a recent study, (1) we performed the self-consistent field (SCF) calculations of a dissipative particle dynamics (DPD) model of conformationally asymmetric diblock copolymer A-B melts to study the stability of various Frank–Kasper (FK) phases formed by such A-B melts. While all the equations in ref (1) are correct, we did not calculate the dimensionless (mean-field) Helmholtz free energies per chain βf c = βu c,χ + βu c,κ – s c /k B and the dimensionless stresses dβf c /dθ of ordered phases accurately. For example, the dimensionless (nonbonded) internal energy per chain due to the system compressibility βu c,κ = (N/2κ)∑ q βû 0 (|q|)(ϕ̂ A (q) + ϕ̂ B (q) – δ q,0 )(ϕ̂ A (−q) + ϕ̂ B (−q) – δ– q,0 ) should be calculated in the reciprocal space as given (note that βû 0 (|q|) is obtained analytically, and that since ϕ A (r) is real 𝜙̂ A (−𝐪)= $\overline{𝜙̂ A}$(𝐪), i.e., the complex conjugate of ϕ̂ A (q)); this corresponds to numerically evaluating βu c,κ = (N/2κ)∫dr(ϕ A (r) + ϕ B (r) – 1)ψ(r) with ψ(r) ≡ (∫dr′βu 0 (|r–r′|)(ϕ A (r′) + ϕ B (r′) – 1)/V in the real space via the (composite) trapezoidal rule. While the trapezoidal rule is less accurate than the Romberg integration (2) in general, the opposite occurs when the integration domain is periodic; this was pointed out by Matsen in ref (3) and explained in detail in refs (4and5). With our ignorance of this exceptional convergence property of the trapezoidal rule for periodic functions, we calculated all convolutions having the form of ψ(r) in the reciprocal space via the fast Fourier transforms (6) but those having the form of βu c,κ in the real space via Romberg integration; the latter leads to inaccurate results and unfortunately incorrect phase diagrams reported in ref (1), which are corrected here.

He, Juntong↗
DOE OSTI2025

Data-Efficient Methods for Determining Flory–Huggins χ Parameters in Multicomponent Polymer Formulations

Polymer formulations are essential in diverse applications including personal care products, coatings, paints, adhesives, and plastic materials. Designing these formulations requires navigating large, complex design spaces, where phase and self-assembly behavior critically impact performance. The Flory–Huggins χ parameter, which quantifies segmental miscibility, is widely used to parametrize the excess free energy of mixing in formulation models. In this work, we introduce two data-efficient, top-down methods for estimating χ parameters using the Random Phase Approximation (RPA): (i) Boundary Nonlinear Regression (Boundary-NLR), which fits theoretical spinodal boundaries to experimental phase boundaries, and (ii) Surrogate Model Inverse Parameter Estimation (SMIPE), which uses a Gaussian Process Classifier to fit sparse phase maps via a surrogate model. Both methods allow rapid parametrization of polymer field-theoretic models without the need for additional experiments. We evaluate these approaches on data sets involving polymer–solvent–nonsolvent ternary mixtures and block copolymer–solvent systems, demonstrating their robustness to experimental noise and their relevance for real-world formulation design.

copolymers↗
DOE OSTI2021

Hybrid Symmetry Epitaxy of the Superconducting Fe(Te,Se) Film on a Topological Insulator

It is challenging to grow an epitaxial 4-fold compound superconductor (SC) on a 6-fold topological insulator (TI) platform due to the stringent lattice-matching requirement. Here, we demonstrate that Fe(Te,Se) can grow epitaxially on a TI (Bi 2 Te 3 ) layer due to accidental, uniaxial lattice match, which is dubbed as “hybrid symmetry epitaxy”. This new growth mode is critical to stabilizing robust superconductivity with TC as high as 13 K. Furthermore, the superconductivity in this FeTe 1–x Se x /Bi 2 Te 3 system survives in the Te-rich phase with Se content as low as x = 0.03 but vanishes at Se content above x = 0.56, exhibiting a phase diagram that is quite different from that of the conventional Fe(Te,Se) systems. This unique heterostructure platform that can be formed in both TI-on-SC and SC-on-TI sequences opens a route to unprecedented topological heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

LASSO for CALPHAD Model Selection Enables Data-Efficient Thermodynamic Modeling: An Application in Thermochemical Hydrogen Production Materials

Phenomenological CALPHAD (CALculation of PHAse Diagrams) models, widely used for multicomponent materials, often contain a considerable number of parameters and require fitting using data from a relatively small number of experimental measurements or theoretical calculations. Sometimes these parameters are introduced for the purpose of improving model fits but without clear physical justification, which leads to overparametrized models with poor generalization performance. Automated approaches for optimal model selection based on the available data therefore become critical. Here, in this work, a least absolute shrinkage and selection operator (LASSO)-based approach is developed for model selection by leveraging the linearity of the CALPHAD model with respect to its parameters to convert the model selection and fitting to a LASSO minimization problem. We demonstrate its utility for thermodynamic modeling of thermochemical hydrogen (TCH) production materials using lanthanum strontium manganite (LSM) as an example. Various TCH-relevant properties, including oxygen stoichiometry as a function of oxygen partial pressure, enthalpy of reduction, and entropy of reduction, are successfully predicted with reasonable accuracy using a minimal set of model parameters. Importantly, the model selection and fitting involve minimal human decision; it can therefore be applied to high-throughput DFT defect calculations and yield efficient workflows for TCH material modeling and optimization.

CALPHAD↗
DOE OSTI2020

Structural Rearrangements of Subnanometer Cu Oxide Clusters Govern Catalytic Oxidation

Sub-nanometer metal oxide clusters are very important materials that are widely used, for example, in catalysis or electronic devices such as sensors. Hence, it is critical to understand the atomic structures and properties of sub-nanometer metal oxide clusters under a reactive gas environment, such as O 2 . We consider here experimentally accessible precise-size Cu clusters (Cu 4 ) supported on partially hydroxylated amorphous alumina and show that such clusters can access, in catalytic conditions at high temperature under a pressure of O 2 , a large ensemble of oxidized structures, representing a large variety of oxygen content, of geometries in link with the support, and of catalytic activities for oxidation reactions, as seen from their reducibilities. A grand canonical basin hopping method based on first-principles energy reveals an ensemble of 24 configurations for the Cu 4 O x cluster of low free energy, less than 0.8 eV above the global minimum. The low free energy ensemble consists of clusters of different stoichiometries, which are mainly Cu 4 O 3 and Cu 4 O 4 in the temperature range of 200–400 °C and under a pressure of 0.5 bar of O 2 . The presence of several competitive isomers at each composition implies that cluster fluxionality impacts the phase diagram, which should be ensemble-averaged. In terms of catalytic oxidation activity, Cu 4 O 3 isomers present highly variable O abstraction energy: the most stable isomers are inactive for alkane oxidative dehydrogenation, but isomerization to metastable isomers, that proceed with low barrier, enable to create active configurations with low O abstraction energy. O atoms with the lowest anionic character, and thus of more electrophilic nature, present the best oxidation capability. In contrast, all Cu 4 O 4 isomers show a low O abstraction energy and a high potential catalytic activity. Here, this manuscript demonstrates the unique structural and electronic properties of sub-nano Cu oxide clusters and illustrates the critical roles of configuration ensembles and rearrangement to highly reactive metastable cluster isomers in nanocatalysis.

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

Hydrogen Evolution on Restructured B-Rich WB: Metastable Surface States and Isolated Active Sites

Borides have been recently identified to be a class of earth-abundant and low-cost materials that are surprisingly active toward oxidative dehydrogenation and hydrogen evolution reaction (HER) catalysis. Here, we explain from first-principles calculations the HER activity of WB, an industrial material known for its superior physical properties and chemical inertness. We find that, compared to bulk termination, a major surface reconstruction takes place, leading to the off-stoichiometric B-rich surface termination that contains the active sites. The hydrogen adsorbate configurations are further investigated under various adsorbate coverages. Many competing configurations appear to be accessible in reaction conditions, and thus, a grand canonical ensemble is established to describe the catalytic system. A phase diagram of adsorbate coverages is constructed as a function of pH and the applied potential. A complex reaction network is presented based on the ensemble-averaged reaction rates, and the electrocatalytic activities are evaluated under different overpotentials. Here, the major contributors to the activity are found to be a few metastable surface states with a distinct electronic structure that are only accessible at high adsorbate coverages in reaction conditions. In addition, while the activity of the dominant active site is nearly the same as that on the unreconstructed WB, the B-rich formations play an important role of isolating the active sites and preventing the passivation of the surface with H 2 bubble formation.

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

Formic Acid Electrooxidation on Pt or Pd Monolayer on Transition-Metal Single Crystals: A First-Principles Structure Sensitivity Analysis

Here, we present a density functional theory analysis of trends for the electrooxidation of formic acid (FAO) on a single Pt or Pd monolayer supported on the close-packed (111) facet of transition metals (Pt*/M or Pd*/M): Au, Ag, Cu, Pt, Pd, Ir, and Rh, the close-packed (0001) facet of Os, Ru, and Re, and the open (100) facet of Au, Ag, Pt, Pd, Ir, and Rh. We show that the deposition of Pt or Pd pseudomorphic monolayers on these single crystals modifies the electronic structure of the Pt or Pd monolayer. Specifically, we found a direct correlation between the d-band center of the Pt and Pd monolayer and the free energy of adsorbed CO*, the latter being a reactivity descriptor for FAO. Together with the free energy of adsorbed OH* as a second reactivity descriptor, we depict the thermochemistry of the reaction network as phase diagrams showing calculated free energies across regions of rate-determining steps. We found that FAO is structure-sensitive on most surfaces studied. Pt*/Au(111) is predicted to be the most active among all Pt*/M(111/0001) surfaces studied, despite binding CO*, the strongest among the close-packed facets. This is the case because of its superiority in activating water to OH*, thus removing CO* at lower potentials than other surfaces. On similar grounds, Pt*/Pd(100), Pd*/Re(0001), Pd*/Au(111), Pd*/Ag(111), and Pd*/Pt(111) are predicted to show higher FAO activity than the corresponding monometallic Pt and Pd surfaces.

(100) facets↗
DOE OSTI2021

Elucidating Mass Transport Regimes in Gas Diffusion Electrodes for CO 2 Electroreduction

Gas diffusion electrodes (GDEs) have shown promising performance for the electrochemical reduction of CO 2 (CO 2 R). In this study, a resolved, pore scale model of electrochemical reduction of CO 2 within a liquid-filled catalyst layer is developed. Three CO 2 mass transport regimes are identified in which the CO 2 penetration depth is controlled by CO 2 consumption in the electrolyte, CO 2 conversion along the solid-electrolyte double-phase boundaries (DPBs), and CO 2 conversion concentrated around the gas–solid–electrolyte triple-phase boundaries (TPBs). While it is possible for CO 2 R to be localized around the TPBs, in systems with submicron pore radii operating at <1 A cm –2 CO 2 R will be distributed across the DPBs within the catalyst layer. This validates the assumption of pore-scale uniformity implicit in popular, volume-averaged GDE models. Overall, the CO 2 conversion efficiency depends strongly on the governing mass transport regime, and operational-phase diagrams are constructed to guide the catalyst layer design.

10 SYNTHETIC FUELS↗
DOE OSTI2021

Influence of Charge Fraction on the Phase Behavior of Symmetric Single-Ion Conducting Diblock Copolymers

A series of symmetric poly[(oligo(ethylene glycol) methyl ether methacrylate-co-oligo(ethylene glycol) propyl sodium sulfonate methacrylate)]-block-polystyrene (PsOEGMA-PS) diblock copolymers were synthesized as a model system to probe the effect of charge fraction on the phase behavior of charged-neutral single-ion conducting diblock copolymers. Small-angle X-ray scattering (SAXS) experiments showed that increasing the charge fraction does not alter the ordered phase morphology (lamellar) but increases the order–disorder transition temperature (T ODT ) significantly. Additionally, the effective Flory–Huggins interaction parameter (χ eff ) was found to increase linearly with the charge fraction, similar to the case of conventional salt-doped diblock copolymers. This indicates that the effect of counterion solvation, attributed to the significant mismatch between the dielectric constant of each block, provides the dominant effect in tuning the phase behavior of this charged diblock copolymer. Lastly, we therefore infer that electrostatic cohesion (local charge ordering induced by Coulombic interactions), which is predicted to suppress microphase separation and lead to asymmetric phase diagrams, only plays a minor role in this model system.

36 MATERIALS SCIENCE↗