Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential theory”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Modified Embedded-Atom Method Potential for a Quaternary Fe-Cr-Si-Mo Solid Solution Alloy

Ferritic-martensitic steels, such as T91, are candidate materials for high-temperature applications, including superheaters, heat exchangers, and advanced nuclear reactors. Considering these alloys’ wide applications, an atomistic understanding of the underlying mechanisms responsible for their excellent mechano-chemical properties is crucial. Here, we developed a modified embedded-atom method (MEAM) potential for the Fe-Cr-Si-Mo quaternary alloy system—i.e., four major elements of T91—using a multi-objective optimization approach to fit thermomechanical properties reported using density functional theory (DFT) calculations and experimental measurements. Elastic constants calculated using the proposed potential for binary interactions agreed well with ab initio calculations. Furthermore, the computed thermal expansion and self-diffusion coefficients employing this potential are in good agreement with other studies. This potential will offer insightful atomistic knowledge to design alloys for use in harsh environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

Freja Studies of the Current-Voltage Relation in Substorm-Related Events

Field-aligned currents and electrostatic potentials play important roles in the coupling between the magnetosphere and the ionosphere. If one assumes that the ionosphere-magnetosphere potential difference is mainly due to the mirror force, one can use the single particle adiabatic kinetic theory to describe the system. From this theory, a linear relationship j(sub II) = KV between field-aligned current density j(sub II) and potential drop V along the same field line can be derived, provided that the potential drop is not too large and not too small. With rare exceptions, observational tests of this relation have mainly concentrated on quiet magnetospheric situations, with acceleration voltages V approx. less than 5 kV. Here we use observations from the Freja satellite of precipitating auroral electrons at 1.700 km altitude to study substorm related events, with acceleration voltages up to 20 keV. The observations are found to be consistent with a linear current-voltage relation even i n these conditions, although with values of the field aligned K lower than previously reported (1-5 x 10(exp 11 S/sq m). This can be explained by lower densities and higher characteristic electron energies in the magnetospheric source region of the precipitating electrons. We analyze the data by three different methods, which are all found to be in general agreement. The results are in agreement with a previous study, where the spectra of precipitating electrons --were indirectly inferred by inversion of data from the EISCAT incoherent scatter radar, thereby validating the use of radar data for studies of auroral electrons. Comparisons with previous studies are made, emphasizing the dependence of the results on the type of auroral structure and magnetospheric conditions.

Olsson, A.↗
DOE OSTI2023

Impact of the Electronic Properties of Chalcogenide Ligands in their Complexation with Uranyl Nitrate Complexes

The idea of covalency being directly correlated with selectivity in f-element separations has gained attention over the years. However, “covalency” as a concept is still under debate as it is not an experimental nor quantum mechanical observable. Regardless, covalency along with ionicity are still the main concepts used to explain the nature of the chemical bond of coordination complexes. Another concept that has gained the attention in the separations community is the hard-soft acid-base (HSAB) theory, which favors soft-donor ligands having the potential to increase selectivity for actinides (An) over the lanthanides (Ln) due to their ability to engage more effectively f-electrons in covalent interactions. However, it has been shown that while using softer donor ligands in fact increased orbital mixing and covalency with the actinides, it compromises the strength of the interaction owing to a decrease in the metal-ligand electrostatic interactions.1 Herein, we study the electronic structure of different types of chalcogenide ligands (LE, E = O, S, Se), their interaction with the uranyl ion, and their complexation free energies. Surprisingly, contrary to the expected order of covalency, i.e. An–LSe > An–LS > An–LO, our results suggest that the degree of covalency and associated complexation free energies depend on the nature of the coordinating moiety rather than solely on the directly coordinated atom. 1) Sadhu, Biswajit, and Michael Dolg. Enhancing actinide (III) over lanthanide (III) selectivity through hard-by-soft donor substitution: exploitation and implication of near-degeneracy-driven covalency. Inorganic Chemistry, 2019, 58, 9738 – 9748.

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

New avenues for |∆ B | = 2 processes beyond neutron-antineutron oscillations

We explore baryon-number-violating (|∆ B | = 2) processes beyond the well-known neutron-antineutron ($n - \bar{n}$) oscillations, focusing on the $Λ - \bar{Λ}$ system. The presence of a strange quark in the Λ baryon introduces a new set of six-quark operators roughly of the form (uds) 2 , which are different from the (udd) 2 operators responsible for oscillations. Using the Standard Model Effective Field Theory (SMEFT), we classify all dimension-9 operators that cause |∆ B | = 2 transitions and study their UV completions mediated by exotic scalar fields with trilinear interactions. We demonstrate that in these models, oscillations can occur at tree level, with $n - \bar{n}$ mixing potentially appearing at higher loop levels. We employ a chiral effective theory to constrain the effective mass mixing δm Λ , deriving bounds from current experimental limits on $n - \bar{n}$ oscillations and dinucleon decays such as pp → K + K + . These bounds indicate that $Λ - \bar{Λ}$ oscillations probe a complementary parameter space, sensitive to baryon-number violation at scales up to 10 2 − 10 3 TeV. We show that the existing indirect bounds make it challenging to provide a competitive bound on δm Λ at BESIII.

Baryon/Lepton Number Violation↗
DOE OSTI2020

Direct Demonstration of Unified Brønsted-Evans-Polanyi Relationships for Proton-Coupled Electron Transfer Reactions on Transition Metal Surfaces

Brønsted-Evans-Polanyi (BEP) relationships, which relate elementary reaction barriers to reaction thermodynamics, have long been discussed in electrochemical science, but demonstration of their existence in heterogeneous electrocatalysis, across many different catalyst surfaces and voltages, is lacking. In this contribution, the BEP model is demonstrated to describe both multiple catalyst surfaces and variable voltages in electrochemical environments. For proton-coupled electron transfer (PCET) to surface nitrogen (N*) and nitric oxide (NO*), reaction energies and activation barriers are calculated using Density Functional Theory (DFT) on a parallel plate capacitor model at three different potentials and for nine different transition metal surfaces. Linear BEP relations that describe all potentials and catalyst surfaces are obtained for these elementary reactions, and the slopes of the correlations are shown to be directly related to the fractional coordinate of the transition states (FCTS) of the reactions. The results, which are explained in terms of Marcus Theory, prove a direct equivalence between unified BEP coefficients, describing both variable catalyst surfaces and voltages, and electrochemical symmetry factors and provide a straightforward means of estimating this quantify for PCET reactions on transition metal surfaces. The resulting relationships, in turn, could lead to predictions of electrocatalytic reactivity trends of enhanced accuracy and efficiency.

Electrochemistry↗
DOE OSTI2024

Enhancing Nitrogen Activation in Electrochemical Reduction: The Role of Rare Earth Oxide Surface Configurations

The pursuit of sustainable ammonia synthesis has prompted the exploration of ambient electrochemical nitrogen reduction reaction (e-NRR) as an alternative to the energy-intensive Haber-Bosch process. Here, this study conducted a theoretical investigation into the use of rare earth oxide materials, specifically dysprosium oxide (Dy 2 O 3 ), as potential electrocatalysts for NRR. Utilizing spin-polarized density functional theory calculations, we explored the interaction between Dy 2 O 3 surfaces and nitrogen (N 2 ) molecules, examining the capability of Dy 2 O 3 to adsorb and activate N 2 under ambient conditions. The results indicate that Dy 2 O 3 surfaces exhibit diverse configurations and bonding environments, providing a variety of reactive sites that display different behaviors in N 2 adsorption and activation. The distinctive electronic structure and surface chemistry of a particular Dy 2 O 3 surface configuration were found to significantly enhance the activation of N 2 by promoting charge transfer, which facilitates the NRR process. This research provides deep insights into the mechanistic pathways of N 2 reduction over Dy 2 O 3 , highlighting the surface properties as pivotal in catalysis. These theoretical insights serve as a foundation for the development of novel rare earth-based electrocatalytic materials for efficient ambient e-NRR, potentially transforming ammonia production into a greener and more energy-efficient process.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1994

Discrete-Layer Piezoelectric Plate and Shell Models for Active Tip-Clearance Control

The objectives of this work were to develop computational tools for the analysis of active-sensory composite structures with added or embedded piezoelectric layers. The targeted application for this class of smart composite laminates and the analytical development is the accomplishment of active tip-clearance control in turbomachinery components. Two distinct theories and analytical models were developed and explored under this contract: (1) a discrete-layer plate theory and corresponding computational models, and (2) a three dimensional general discrete-layer element generated in curvilinear coordinates for modeling laminated composite piezoelectric shells. Both models were developed from the complete electromechanical constitutive relations of piezoelectric materials, and incorporate both displacements and potentials as state variables. This report describes the development and results of these models. The discrete-layer theories imply that the displacement field and electrostatic potential through-the-thickness of the laminate are described over an individual layer rather than as a smeared function over the thickness of the entire plate or shell thickness. This is especially crucial for composites with embedded piezoelectric layers, as the actuating and sensing elements within these layers are poorly represented by effective or smeared properties. Linear Lagrange interpolation polynomials were used to describe the through-thickness laminate behavior. Both analytic and finite element approximations were used in the plane or surface of the structure. In this context, theoretical developments are presented for the discrete-layer plate theory, the discrete-layer shell theory, and the formulation of an exact solution for simply-supported piezoelectric plates. Finally, evaluations and results from a number of separate examples are presented for the static and dynamic analysis of the plate geometry. Comparisons between the different approaches are provided when possible, and initial conclusions regarding the accuracy and limitations of these models are given.

Heyliger, P. R.↗
DOE OSTI2023

Preparations for quantum simulations of quantum chromodynamics in 1+1 dimensions. II. Single-baryon β-decay in real time

Here, a framework for quantum simulations of real-time weak decays of hadrons and nuclei in a two-flavor lattice theory in one spatial dimension is presented. A single generation of the Standard Model is found to require 16 qubits per spatial lattice site after mapping to spin operators via the Jordan-Wigner transformation. Both quantum chromodynamics and flavor-changing weak interactions are included in the dynamics, the latter through four-Fermi effective operators. Quantum circuits that implement time evolution in this lattice theory are developed and run on Quantinuum’s H1-1 20-qubit trapped ion system to simulate the β-decay of a single baryon on one lattice site. These simulations include the initial state preparation and are performed for both one and two Trotter time steps. The potential intrinsic error-correction properties of this type of lattice theory are discussed, and the leading lattice Hamiltonian required to simulate 0νββ-decay of nuclei induced by a neutrino Majorana mass term is provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1987

Theoretical approach to oxygen atom degradation of silver

Based on available Rutherford backscattering spectrometry (RBS), proton induced X-ray emission (PIXE) and ellipsometry data obtained on silver specimens subjected to atomic oxygen attack in low Earth orbit STS flight 41-G, a theory was developed to model the oxygen atom degradation of silver. The diffusion of atomic oxygen in a microscopically nonuniform medium is an essential constituent of the theory. The driving force for diffusion is the macroscopic electrochemical potential gradient developed between the specimen surface exposed to the ambient and the bulk of the silver specimen. The longitudinal electric effect developed parallel to the gradient is modified by space charge of the diffusing charged species. Lateral electric fields and concentration differences also exist due to the nonuniform nature of the medium. The lateral concentration differences are found to be more important than the lateral electric fields in modifying the diffusion rate. The model was evaluated numerically. Qualitative agreement exists between the kinetics predicted by the theory and kinetic data taken in ground-based experiments utilizing a plasma asher.

Fromhold, Albert T., Jr.↗
NASA NTRS1990

Application of symbolic computations to the constitutive modeling of structural materials

In applications involving elevated temperatures, the derivation of mathematical expressions (constitutive equations) describing the material behavior can be quite time consuming, involved and error-prone. Therefore intelligent application of symbolic systems to faciliate this tedious process can be of significant benefit. Presented here is a problem oriented, self contained symbolic expert system, named SDICE, which is capable of efficiently deriving potential based constitutive models in analytical form. This package, running under DOE MACSYMA, has the following features: (1) potential differentiation (chain rule), (2) tensor computations (utilizing index notation) including both algebraic and calculus; (3) efficient solution of sparse systems of equations; (4) automatic expression substitution and simplification; (5) back substitution of invariant and tensorial relations; (6) the ability to form the Jacobian and Hessian matrix; and (7) a relational data base. Limited aspects of invariant theory were also incorporated into SDICE due to the utilization of potentials as a starting point and the desire for these potentials to be frame invariant (objective). The uniqueness of SDICE resides in its ability to manipulate expressions in a general yet pre-defined order and simplify expressions so as to limit expression growth. Results are displayed, when applicable, utilizing index notation. SDICE was designed to aid and complement the human constitutive model developer. A number of examples are utilized to illustrate the various features contained within SDICE. It is expected that this symbolic package can and will provide a significant incentive to the development of new constitutive theories.

Arnold, Steven M.↗
DOE OSTI2020

Needle-Like Ferroelastic Domains in Individual Ferroelectric Nanoparticles

Superior structural, physical and electronic properties make ferroelectric nanocrystals essential in enabling a range of next-generation devices. Ferroelectric responses are determined by crystal structure and domain morphology. The ability to reversibly displace, create and annihilate elastic domains is critical to device applications. Whereas electric-field control has been demonstrated for ferroelectric 180° surface domain walls and vortices, similar control of ferroelastic domains and domain boundaries within individual nanocrystals remains challenging. Using controlled external compressive and tensile axial stress, deterministic and reversible control of highly mobile ferroelastic domains and axial polarization in three dimensions is demonstrated in this study. While many studies exist on ferroelastic domains in thin films and bulk, little is known about ferroelastic interactions at the single nanoparticle level, especially involving domain boundaries. Through combining Bragg coherent X-ray diffractive imaging and Landau theory, strain gradients in individual BaTiO3 nanocrystals are shown to stabilize needle-like ferroelastic twin domains. These domains are highly labile under applied axial stress, producing a locally enhanced electric polarization mediated by a ferroelectric phase transition. The efficacy of Bragg coherent X-ray diffractive imaging in studying in operando domains in three-dimensions is demonstrated, while synergy with theory provides a paradigm for domain boundary engineering and potential for nanoscale functional devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ambipolar Electrochemistry of Pre‐Intercalated Ti 3 C 2 T x MXene in Ionic Liquid Electrolyte

Abstract Cation intercalation with or without redox remains the dominant charge storage mechanism for two‐dimensional (2D) Ti 3 C 2 T x MXene. Anion‐based charge storage remains unexplored due to intrinsic negative surface charge of MXenes preventing spontaneous intercalation of anions and irreversible oxidation of Ti at anodic potentials in aqueous electrolytes. In this work, we report on the ambipolar electrochemical behavior of the Ti 3 C 2 T x in ionic liquid electrolyte over a 2.5 V electrochemically stable window. The experiments are conducted on a thin Ti 3 C 2 T x film current collector coated with an electroactive layer of small flakes (∼150 nm) of Ti 3 C 2 T x pre‐intercalated with 1‐ethyl‐3‐methylimidazolium bis‐(trifluoromethylsulfonyl)‐imide (EMIM‐TFSI) ionic liquid. Couples of redox peaks with a very small potential separation during the voltage sweep are observed at high negative (−0.75 V vs. Ag wire) and high positive (+0.75 V vs. Ag wire) potentials. Our experimental electrochemical data combined with density functional theory (DFT) calculations suggest feasibility of pseudo‐intercalation of TFSI anions between Ti 3 C 2 T x flakes. This study provides a pathway for elucidating anion intercalation for different MXene chemistries in solvent‐free electrolytes, which can lead to development of MXene based energy storage devices with improved performance.

Electrochemistry↗
DOE OSTI2021

The role of precursor decomposition in the formation of samarium doped ceria nanoparticles via solid-state microwave synthesis

The impact on the final morphology of ceria (CeO 2 ) nanoparticles made from different precursors (commercial: cerium acetate/nitrate) and in house: cerium tri(methylsilyl)amide (Ce-TMSA)) via a microwave solid state reaction has been determined. In all instances, powder X-ray diffraction indicated that the cubic fluorite CeO 2 phase (PDF# 04–004-9150, with the space group Fm-3 m) had formed. Scanning electron microscopy (SEM) images revealed spherical nanoparticles were produced from the Ce-TMSA precursor. The commercial acetate and nitrate precursors produced particles with irregular morphology. The roles of the precursor decomposition and binding energy in the synthesis of the nanocrystals with various morphologies, as well as a possible growth mechanism, were evaluated based on experimental and computational data. The formation of spherical shaped nanoparticles was determined to be due to the preferential single-step decomposition of the Ce-TMSA as well as the low activation energy to overcome decomposition. Due to the complicated decomposition of the commercial precursors and high activation energy the resulting particles adopted an irregular morphology. Highly uniform samarium doped ceria (Sm x Ce 1-x O 2-δ ) nanospheres were also synthesized from Ce-TMSA and samarium tri(methylsilyl)amide (Sm-TMSA). The effects of reaction time and temperature, on the final morphology were observed through SEM. The rapid single-step decomposition of TMSA-based precursors as observed through thermogravimetric analysis (TGA) and confirmed through the calculation of potential energy surfaces and binding energies from density functional theory (DFT) calculations, indicated that nanoparticle formation follows LaMer’s classical nucleation theory.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Bayesian optimization of entropy-stabilized metal fluoride conversion cathodes and their synthesis

Fluoride-based conversion cathodes are promising for next-generation Li ion batteries because of their high voltage and energy densities. High entropy earth abundant fluoride cathodes are attractive for sustainable battery because of the potentially higher cycling stability. While Ni and Co containing equimolar alloys stabilized by high entropy were considered previously, Ni and Co free non-equimolar metal alloy compositions might offer new opportunities to increase materials stability and average voltage. This opportunity was pursued using Bayesian optimization of average voltage for the computational prediction of (Fe x1 Cu x2 Zn x3 Mn x4 Mg x5 )F 2 cathode chemistries. Theoretical voltage, room temperature free energy, and electronic structure of potentially high voltage compositions were computed using density functional theory calculations. Here, the experimental realization of non-equimolar quinary Co and Ni free fluorides with stable rutile single phase and improved conversion reaction voltage by up to 28 % compared to the equimolar alloy demonstrated the power of Bayesian optimization methods for the computationally guided experimental synthesis of fluoride cathodes. Finally, non-equimolar cathode chemistries were found with cycling stability superior to known multi-element fluorides.

25 ENERGY STORAGE↗
DOE OSTI2024

Minimum-Energy Conical Intersections by Compressed Multistate Pair-Density Functional Theory

Compressed multistate pair-density functional theory (CMS-PDFT) is a multistate version of multiconfiguration pair-density functional theory that can capture the correct topology of coupled potential energy surfaces (PESs) around conical intersections. Here, in this work, we develop interstate coupling vectors (ISCs) for CMS-PDFT in the OpenMolcas and PySCF/mrh electronic structure packages. Yet, the main focus of this work is using ISCs to calculate minimum-energy conical intersections (MECIs) by CMS-PDFT. This is performed using the projected constrained optimization method in OpenMolcas, which uses ISCs to restrain the iterations to the conical intersection seam. We optimize the S 1 /S 0 MECIs for ethylene, butadiene, and benzene and show that CMS-PDFT gives smooth PESs in the vicinities of the MECIs. Furthermore, the CMS-PDFT MECIs are in good agreement with the MECI calculated by the more expensive XMS-CASPT2 method.

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

Collisional Excitation of HCN by CO to Refine the Modeling of Cometary Comae

Here, we present the first dataset of collisional (de)-excitation rate coefficients of HCN induced by CO, one of the main perturbing gases in cometary atmospheres. The dataset spans the temperature range of 5–50 K. It includes both state-to-state rate coefficients involving the lowest ten and nine rotational levels of HCN and CO, respectively, and the so-called “thermalized” rate coefficients over the rotational population of CO at each kinetic temperature. The derivation of these coefficients exploited the good performance of the statistical adiabatic channel model (SACM) on top of an accurate interaction potential computed at the CCSD(T)-F12b/CBS level of theory. The reliability of the SACM approach was validated by comparison with full quantum calculations restricted at the lowest total angular momentum of the system. These results provide essential input to accurately model the distribution among the rotational energy levels and the abundance of HCN in cometary atmospheres, accounting for deviations from local thermodynamic equilibrium that typically occurs in such environments.

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

Machine-Learned Force Field for Molecular Dynamics Simulations of Nonequilibrium Ammonia Synthesis on Iron Catalysts

Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.

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

Relationship Between Hydrogen-Bonding Motifs and the 1b1 Splitting in the X-ray Emission Spectrum of Liquid Water

The split of the 1??! peak observed in the X-ray emission (XE) spectrum of liquid water has been the focus of intense research over the last two decades. Although several hypotheses have been proposed to explain the origin of the 1??! splitting, a general consensus has not yet been reached. In this study, we introduce a novel theoretical/computational approach which, combining pathintegral molecular dynamics (PIMD) simulations carried out with the MB-pol potential energy function and time-dependent density functional theory (TD-DFT) calculations, correctly predicts the split of the 1??! peak in liquid water and not in crystalline ice. A systematic analysis in terms of the underlying local structure of liquid water at ambient conditions indicates that several different hydrogen-bonding motifs contribute to the overall XE lineshape in the energy range corresponding to emissions from the 1??! orbitals, which suggests that it is not possible to unambiguously attribute the split of the 1??! peak to only two specific structural arrangements of the underlying hydrogen-bonding network.

Cruzeiro, Vinicius↗