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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Spin crossover transition driven by pressure: Barocaloric applications

This article describes a mean-field theoretical model for Spin-Crossover (SCO) materials and explores its implications. It is based on a simple Hamiltonian that yields the high spin molar fraction as a function of temperature and pressure, as well as a temperature–pressure phase diagram for the SCO transition. In order to test the model, it was compared with the giant Barocaloric Effect (BCE) of the SCO material [FeL 2 ][B 4 ] 2 . Furthermore, we found that optical phonons are responsible for 92% of the total barocaloric entropy change. DFT calculations successively indicates that, as expected, the majority of this effect can be traced to low frequencies modes of vibration (400 cm -1 ), associated to the Fe coordination.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Ordinary and exotic mesons in the extended Linear Sigma Model

The extended Linear Sigma Model (eLSM) is a hadronic model based on the global symmetries of QCD and the corresponding explicit, anomalous, and spontaneous breaking patterns. In its basic three-flavor form, its mesonic part contains the dilaton/glueball as well as the nonets of pseudoscalar, scalar, vector, and axial–vector mesons, thus chiral symmetry is linearly realized. In the chiral limit and neglecting the chiral anomaly, only one term – within the dilaton potential – breaks dilatation invariance, and all terms are chirally symmetric. Spontaneous symmetry breaking is implemented by a generalization of the Mexican-hat potential, with explicit symmetry breaking responsible for its tilting. The overall mesonic phenomenology up to ~2 GeV is in agreement with the PDG compilation of masses and partial and total decay widths. The eLSM was enlarged in a straightforward way to include other conventional quark–antiquark nonets (pseudovector and orbitally excited vector mesons, tensor and axial-tensor mesons, radially excited (pseudo)scalar mesons, etc.), as well as two nonets of hybrid mesons, the lightest one with exotic quantum numbers J $\mathscr{PC}$ = 1 -+ not allowed for $\bar{q}q$ objects, such as the resonance π 1 (1600) and the recently discovered η 1 (1855) . In doing so, different types of chiral multiplets are introduced: heterochiral and homochiral multiplets, which differ in the way they transform under chiral transformations. Moreover, besides the scalar glueball that is present from the beginning as dilaton, other glueballs, the tensor, the pseudoscalar and the vector glueballs were coupled to the eLSM: the scalar resonance ƒ 0 (1710) turns out to be mostly gluonic, the tensor glueball couples strongly to vector mesons, and the pseudoscalar glueball couples sizably to $ππη'$ and can be assigned to X (2370) or X (2600). In all cases above, masses and decays can be analyzed allowing for a better understanding of both conventional and non-conventional mesons: whenever data are available, a comparison is performed and, when this is not the case, predictions of decay widths and decay ratios are outlined. The eLSM contains chiral partners on an equal footing and is therefore well suited for studies of chiral symmetry restoration at nonzero temperature and densities: this is done by coupling it to the Polyakov loop. Finally, the QCD phase diagram and the location of the critical endpoint were investigated within this framework.

Chiral symmetry↗
DOE OSTI2024

Type-II superconductivity at 9K in Pb–Bi alloy

In the present work, we report the synthesis of Pb–Bi alloy with enhanced Tc of up to 9K, which is higher than that of Pb. The alloy is synthesized via a solid-state reaction route in the vacuum-encapsulated quartz tube at 700°C in an automated furnace. The synthesized sample is characterized by X-ray Diffraction(XRD) and Energy dispersive X-ray analysis(EDAX) for its phase purity and elemental composition. Rietveld refinement of XRD reveals that the end product is a majority hexagonal Pb 7 Bi 3 , with minor rhombohedral Bi. The electronic transport measurement shows metallic behavior with the Debye temperature of 108K and a superconductivity transition temperature (T c ) below 9K, which is the maximum to date for any reported Pb–Bi alloy, Pb or Bi at ambient pressure. Partial substitution of Bi at the Pb site may modify the free density of electronic states within the BCS model to attain the optimum T c , which is higher by around 2K from the reported T c of Pb. The superconductor phase diagram derived from magneto-transport measurements reveals that the synthesized alloy is a conventional superconductor with an upper critical field (H c2 ) of 3.9 T, which lies well within the Pauli paramagnetic limit. The magnetization measurements carried out following ZFC(Zero Field Cool) protocols infer that the synthesized alloy is a bulk superconductor below 9K. The isothermal M-H(Magnetization vs. Field) measurements performed below T c establish it as a type-II superconductor. Furthermore, the specific heat capacity measurements show that the Pb–Bi alloy is a strongly coupled bulk superconductor below around 9K with possibly two superconducting gaps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Edge-Propagation Discharge Mechanism in CFx Batteries—A First-Principles and Experimental Study

Graphite fluoride (CF x ) cathodes coupled with lithium anodes yield one of the highest theoretical specific capacities (>860 mAh/g) among primary batteries. In practice, the observed discharge voltage (~2.5 V) is significantly lower than thermodynamic limits (>4.5 V), the discharge rate is low, and so far Li/CF x has only been used in primary batteries. Understanding the discharge mechanism at atomic length scales will improve practical CFx energy density, rate capability, and rechargeability. So far, purely experimental techniques have not identified the correct discharge mechanism or explained the discharge voltage. In this work, we apply density functional theory calculations to demonstrate that a CF x -edge propagation discharge mechanism based on lithium insertion at the CF/C boundary in partially discharged CF x exhibits a voltage range of 2.5 to 2.9 V—depending on whether solvent molecules are involved. The voltages and solvent dependence agree with our discharge and galvanostatic intermittent titration technique measurements. The predicted discharge kinetics are consistent with CF x operations. Finally, we predict some Li/CF x rechargeability under the application of high potentials, along a charging pathway different from that of discharge. Our work represents a general, quasi-kinetic framework to understand the discharge of conversion cathodes, circumventing the widely used phase diagram approach which most likely does not apply to Li/CF x because equilibrium conditions are not attained in this system.

25 ENERGY STORAGE↗
DOE OSTI2021

Understanding the Reactions Between Fe and Se Binary Diffusion Couples

Spurred by recent discoveries of high-temperature superconductivity in Fe-Se based materials, the magnetic, electronic, and catalytic properties of iron-chalcogenides have drawn significant attention. Furthermore, much remains to be understood about the sequence of phase formation in these systems. In this work, we shed light on this issue by preparing a series of binary Fe-Se ultrathin diffusion couples via designed thin film precursors and investigating their structural evolution as a function of composition and annealing temperature. Two previously unreported Fe-Se phases crystallized during the deposition process on a nominally room-temperature Si substrate in the 27-33% and 37-47% Fe (atomic percent) composition regimes. Both phases completely decompose after annealing to 200°C in a nitrogen glovebox. At higher temperatures, the sequence of phase formation is governed by Se loss in the annealing process, consistent with what would be expected from the phase diagram. Films rich in Fe (53-59% Fe) crystalized during deposition as β-FeSe (P4/nmm) with preferred c-axis orientation to the amorphous SiO 2 substrate surface, providing a means to non-epitaxial self-assembly of crystallographically aligned, iron-rich β-FeSe for future research. Our findings suggest the crystallization of binary Fe-Se compounds at room temperature via near diffusionless transformations should be a significant consideration in future attempts to prepare metastable ternary and higher order compounds containing Fe and Se.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Weak-Bonding Elements Lead to High Thermoelectric Performance in BaSnS 3 and SrSnS 3 : A First-Principles Study

SnS2, an earth-abundant and ecofriendly material, is limited as a thermoelectric material because of the high lattice thermal conductivity κ L and low carrier mobility μ. By introducing weak-bonding elements Ba or Sr into the SnS 2 framework, we discovered two SnS 2 -based materials BaSnS3 and SrSnS3 with the calculated low κL values of 0.15 and 0.17 W m -1 K -1 , respectively, along the a-axis. The low group velocity and high lattice anharmonicity originating from the weakened and distorted Sn–S bonding network are found in both systems. Moreover, the vibrations of Ba and Sr induce low-lying optical phonons, which strongly couple with the acoustic phonons and strengthen the phonon scattering rates. Compared to SnS 2 , both compounds present lower single-band effective masses, smaller deformation potential constants, and better band convergence, which enhance μ with an insignificantly reduced effective mass. By solving the linearized Boltzmann transport equation with a nonempirical carrier lifetime, we predict excellent ZT values of 2.89 and 2.77 along the a-axis at 900 K in BaSnS 3 and SrSnS 3 , respectively. Further phase diagram calculations of Ba 1–x Sr x SnS 3 solid solutions propose a new compound, Ba 0.5 Sr 0.5 SnS 3 , with an even higher ZT of 3.0. Our work analyzes explicitly how weak-bonding elements enhance μ and suppress κL simultaneously in SnS 2 -analogous systems with a series of compounds nominated as potential high-performance thermoelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Tuning the Morphology and Electronic Properties of Single-Crystal LiNi 0.5 Mn 1.5 O 4-δ : Exploring the Influence of LiCl–KCl Molten Salt Flux Composition and Synthesis Temperature

Single-crystal materials have played a unique role in the development of high-performance cathode materials for Li batteries due to their favorable chemomechanical stability. The molten salt synthesis method has become one of the most prominent techniques used to synthesize single-crystal layered and spinel materials. In this work, the molten salt synthesis method is used as a technique to tune both the morphology and Mn 3+ content of high-voltage LiNi 0.5 Mn 1.5 O 4 (LNMO) cathodes. The resulting materials are thoroughly characterized by a suite of analytical techniques, including synchrotron X-ray core-level spectroscopy, which are sensitive to the material properties on multiple length scales. Furthermore, the multidimensional characterization allows us to build a materials library according to the molten salt phase diagram as well as to establish the relationship among synthesis, material properties, and battery performance. The results of this work show that the Mn 3+ content is primarily dependent on the synthesis temperature and increases as the temperature is increased. Additionally, the particle morphology is mostly dependent on the composition of the molten salt flux, which can be tailored to obtain well-defined octahedrons enclosed by (111) facets, plates with predominant ($11\tilde{2}$) facets, irregularly shaped particles, or mixtures of these. The electrochemical measurements indicate that the Mn 3+ content has a larger contribution to the battery performance of LNMO than do morphological characteristics and that a significant amount of Mn 3+ could become detrimental to the battery performance. However, with similar Mn 3+ contents, morphology still plays a role in influencing the battery cycle life and rate performance. The insights of molten salt synthesis parameters on the formation of LNMO, with deconvolution of the roles of Mn 3+ and morphology, are crucial to continuing studies in the rational design of LNMO cathode materials for high-energy Li batteries.

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

Confusion-Driven Machine Learning of Structural Phases of a Flexible, Magnetic Stockmayer Polymer

We use a semisupervised, neural-network-based machine learning technique, the confusion method, to investigate structural transitions in magnetic polymers, which we model as chains of magnetic colloidal nanoparticles characterized by dipole–dipole and Lennard-Jones interactions. As input for the neural network, we use the particle positions and magnetic dipole moments of equilibrium polymer configurations, which we generate via replica-exchange Wang–Landau simulations. We demonstrate that by measuring the classification accuracy of neural networks, we can effectively identify transition points between multiple structural phases without any prior knowledge of their existence or location. We corroborate our findings by investigating relevant conventional order parameters. Our study furthermore examines previously unexplored low-temperature regions of the phase diagram, where we find new structural transitions between highly ordered helicoidal polymer configurations.

36 MATERIALS SCIENCE↗
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↗