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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Monitoring Local Electric Fields using Stark Shifts on Napthyl Nitrile-Functionalized Silicon Photoelectrodes

in this paper, we report spectroscopic measurements of the local electric field using vibrational Stark shifts of napthyl nitrile-functionalized silicon under electrochemical working conditions. The C$\equiv$N bond is particularly sensitive to applied electric fields and serves as a good probe for the local electric fields at the silicon–aqueous interface. Here, surface-enhanced Raman spectra (SERS) are collected at a silicon surface using a water immersion lens as a function of the reference potential in a three-terminal potentiostat. In deionized (DI) water and KCl solutions, the nitrile (i.e., C$\equiv$N) stretch downshifts by 4.7 and 8.6 cm –1 , respectively, under an applied potential of –1 V vs Ag/AgCl. Density functional theory (DFT) calculations of the napthyl nitrile complex carried out under various electric fields establish the Stark tuning rate to be 0.5622 cm–1/(MV cm –1 ). Based on this relation, electric fields of –8.4 and –15.2 MV/cm were obtained under negative applied potentials. These measurements report the electric field strength within the double (i.e., Helmholtz) layer, which is responsible for pulling positively charged ions (e.g., H + ) toward the surface in reduction reaction processes.

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

Atomic-resolution imaging as a mechanistic tool for studying single-site heterogeneous catalysis

Heterogeneous catalysts dominate the chemical industry but typically feature diverse, incompletely defined active sites. Thus, describing structure-activity relationships, unlike homogeneous catalysts, remains challenging. In contrast, molecularly defined single-site heterogeneous catalysts (SSHCs), using appropriate tools, are poised to address these challenges and provide new avenues for catalysis research and development. The present study explores eco-friendly H 2 production mediated by discrete MoO 2 sites supported on carbon nanohorns (CNHs) and active for alcohol dehydrogenation. While informative, detailed ensemble EXAFS/XANES/XPS, kinetic measurements, and DFT analysis alone cannot provide a full molecular picture of the reaction pathway. Here, using single-molecule atomic-resolution time-resolved electron microscopy (SMART-EM), we propose the identification of four key catalytic intermediates anchored to CNHs and uncover a new reaction pathway involving alkoxide/hemiacetal equilibration and acetal oligomerization. Furthermore, these intermediates are inferred through a combination of theory and SMART-EM, showcasing the potential of SMART-EM as a complementary tool for exploring mechanistic hypotheses in catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Factors Affecting the Electron Conductivity in Single Crystal Li 7 La 3 Zr 2 O 12 and Li 7 P 3 S 11

One of the serious challenges in all solid-state Li ion batteries is neutral Li intrusion into the solid-state electrolyte that can ultimately cause catastrophic failure. One possibility for this is due to n-type electron conductivity that induces the reaction Li + + e – → Li 0 at sites where the potential is less than the Li + /Li potential. This paper reports hybrid density functional theory calculations of the electronic conductivity in two prototype single crystalline solid-state electrolytes, cubic Li 7 La 3 Zr 2 O 12 (c-LLZO) and Li 7 P 3 S 11 (LPS). The formation energies of important point defects that can affect electron conductivity are determined, and we find that the mechanism of n-type electron conductivity for both solid-state electrolytes is via “small” electron polaron hopping, where the quotes signify that substantial Li ion rearrangement is associated with the polaron formation and its migration. In both electrolytes, the formation energies for the small polarons at the Fermi energy are too high to generate measurable electron conductivity at room temperature. For c-LLZO, the concentration of electron polarons necessary to ensure charge neutrality from positively charged oxygen vacancies formed in synthesis can be significantly higher. Hence, the electron conductivity could be significant when measured with ion-blocking metal electrodes, and we discuss how the synthesis conditions could affect this magnitude. However, in the solid-state battery, these polarons are replaced by negatively charged Li vacancies so that the electron conductivity should remain minimal. For LPS single crystals, the inherent minimal electron conductivity is independent of synthesis conditions. Here we also show that the cost of forming Li 0 in bulk c-LLZO is enormous due to strain effects so that it could only potentially form at voids, grain boundaries, or around vacancy defects which relax the lattice strain.

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

Melting point of iron at high pressure: An assessment of uncertainties and effect of electronic temperature

An accurate calculation of the melting point of iron at various pressures in the Earth's core is important for understanding the core structure, geodynamo, and the Earth's history. Previous studies have assessed the melt line of iron at these extreme conditions using various experimental measurement techniques as well as both ab initio and classic molecular dynamics simulations. However, experimental measurements have uncertainties up to several hundred Kelvin, and inconsistencies remain among simulation results. Here in this work, we propose an iterative framework that couples density functional theory (DFT) calculations and molecular dynamics simulations performed using an ensemble of interatomic potentials to assess the effect of electronic temperature on the melting point. We systematically validate the potentials by comparing lattice constants and phonon dispersion curves at 0 K and enthalpy differences between liquid and HCP, FCC, BCC phases of iron close to the melt line at 300 GPa with DFT. Our results show that HCP iron melts at 6144 K (at 300 GPa), BCC phase is thermodynamically unstable, and FCC is metastable at this temperature. The melting points of FCC and BCC phases at 300 GPa are 5858 and 5647 K, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Demonstration of an x-ray Raman spectroscopy setup to study warm dense carbon at the high energy density instrument of European XFEL

We present a proof-of-principle study demonstrating x-ray Raman Spectroscopy (XRS) from carbon samples at ambient conditions in conjunction with other common diagnostics to study warm dense matter, performed at the high energy density scientific instrument of the European x-ray Free Electron Laser (European XFEL). We obtain sufficient spectral resolution to identify the local structure and chemical bonding of diamond and graphite samples, using highly annealed pyrolytic graphite spectrometers. Due to the high crystal reflectivity and XFEL brightness, we obtain signal strengths that will enable accurate XRS measurements in upcoming pump–probe experiments with a high repetition-rate, where the samples will be pumped with high-power lasers. Molecular dynamics simulations based on density functional theory together with XRS simulations demonstrate the potential of this technique and show predictions for high-energy-density conditions. Our setup allows simultaneous implementation of several different diagnostic methods to reduce ambiguities in the analysis of the experimental results, which, for warm dense matter, often relies on simplifying model assumptions. The promising capabilities demonstrated here provide unprecedented insights into chemical and structural dynamics in warm dense matter states of light elements, including conditions similar to the interiors of planets, low-mass stars, and other celestial bodies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A Transferable Force Field for Simulating Adsorption in Metal–Organic Frameworks with Open Metal Sites Based on the 12–6–4 Lennard-Jones Potential

Metal−organic frameworks (MOFs) that contain coordinatively unsaturated open metal sites (OMSs) provide strong host− guest interactions, making them promising sorbents for low-concentration gas adsorption applications such as direct air capture and atmospheric water harvesting. However, accurately modeling host−guest interactions involving OMSs remains challenging for classical force fields (FFs) based on the 12−6 Lennard−Jones (LJ) potential, as the polarization effect of the guest molecule induced by the positively charged OMS is not considered. Here, we introduce an FF based on the 12−6−4 LJ potential, which incorporates charge−induced dipole interactions and is parametrized against a diverse set of host−guest potential energy surfaces (PESs) obtained from density functional theory (DFT). The resulting FF, trained on a generic trimetallic cluster, performs well in both host−guest binding energetics and gas adsorption isotherms across different OMS-containing MOFs, including MOF-74 series and Cu-BTC. These results highlight the excellent transferability of our approach and its potential to enhance the accuracy and robustness of high-throughput MOF discovery workflows, particularly for gas adsorption and separation in large and diverse MOF databases.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Adsorbate Free Energies from DFT-Derived Translational Energy Landscapes

Adsorption free energies are fundamental to surface chemistry and catalysis. Standard models combine some assumed analytical form of the translational potential energy surface, often parameterized against density functional theory (DFT) calculations, with an analytical expression for the resultant translational densities of states (DOS), free energy, and entropy. Here we compare the performance of such models against numerical evaluations of the DOS and thermodynamic functions derived from solutions to the translational Schrodinger equation. We compare results for a translational potential energy surface (PES) derived from nudged eleastic band calculations with those obtained from adsorbate rastering across a series of monatomic (O, S, C, N, H) and polyatomic (NH x ) adsorbates on (100) Pt and Au facets. We find that analytical models as commonly parameterized have mixed performance for describing the translational PES and that the consequences for computed free energies are modest but potentially significant in microkinetic models. Here, numerical solutions are possible for modest to no additional computational cost over analytical models and thus should be considered when reliable free energy estimates are needed or translational potential energy surfaces are available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Neural network potential from bispectrum components: A case study on crystalline silicon

In this article, we present a systematic study on developing machine learning force fields (MLFFs) for crystalline silicon. While the main-stream approach of fitting a MLFF is to use a small and localized training set from molecular dynamics simulations, it is unlikely to cover the global features of the potential energy surface. Additionally, to remedy this issue, we used randomly generated symmetrical crystal structures to train a more general Si-MLFF. Furthermore, we performed substantial benchmarks among different choices of material descriptors and regression techniques on two different sets of silicon data. Our results show that neural network potential fitting with bispectrum coefficients as descriptors is a feasible method for obtaining accurate and transferable MLFFs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multi-fidelity learning for interatomic potentials: low-level forces and high-level energies are all you need

The promise of machine learning interatomic potentials (MLIPs) has led to an abundance of public quantum mechanical (QM) training datasets. The quality of an MLIP is directly limited by the accuracy of the energies and atomic forces in the training dataset. Unfortunately, most of these datasets are computed with relatively low-accuracy QM methods, e.g. density functional theory with a moderate basis set. Due to the increased computational cost of more accurate QM methods, e.g. coupled-cluster theory with a complete basis set (CBS) extrapolation, most high-accuracy datasets are much smaller and often do not contain atomic forces. The lack of high-accuracy atomic forces is quite troubling, as training with force data greatly improves the stability and quality of the MLIP compared to training to energy alone. Because most datasets are computed with a unique level of theory, traditional single-fidelity (SF) learning is not capable of leveraging the vast amounts of published QM data. In this study, we apply multi-fidelity learning (MFL) to train an MLIP to multiple QM datasets of different levels of accuracy, i.e. levels of fidelity. Specifically, we perform three test cases to demonstrate that MFL with both low-level forces and high-level energies yields an extremely accurate MLIP—far more accurate than a SF MLIP trained solely to high-level energies and almost as accurate as a SF MLIP trained directly to high-level energies and forces. Therefore, MFL greatly alleviates the need for generating large and expensive datasets containing high-accuracy atomic forces and allows for more effective training to existing high-accuracy energy-only datasets. Indeed, low-accuracy atomic forces and high-accuracy energies are all that are needed to achieve a high-accuracy MLIP with MFL.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mathematical Foundations for Nonlocal Interface Problems: Multiscale Simulations of Heterogeneous Materials (Final LDRD Report)

Nonlocal models provide a much-needed predictive capability for important Sandia mission applications, ranging from fracture mechanics for nuclear components to subsurface flow for nuclear waste disposal, where traditional partial differential equations (PDEs) models fail to capture effects due to long-range forces at the microscale and mesoscale. However, utilization of this capability is seriously compromised by the lack of a rigorous nonlocal interface theory, required for both application and efficient solution of nonlocal models. To unlock the full potential of nonlocal modeling we developed a mathematically rigorous and physically consistent interface theory and demonstrate its scope in mission-relevant exemplar problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Conserving asphalt resources: Rethinking rejuvenator performance evaluation through peptizing efficiency

Timely rejuvenation and restoration of asphalt are essential conservation practices that help preserve and extend the service life of roads, bridges, driveways, and parking lots. The performance of asphalt rejuvenators is often assessed based on their diffusion rates and softening power, yet these metrics alone fail to capture true rejuvenation potential. Here, this study integrates density functional theory (DFT) modeling with experimental analysis to show that rejuvenation effectiveness is primarily governed by molecular interactions with oxidized asphaltene nanoaggregates. DFT results revealed that amide- and unsaturated-chain compounds, such as hexadecanamide, oleic acid, and 9,17-octadecadienal, act cooperatively to disrupt π–π stacking and exfoliate asphaltene layers, reducing binding strength and enhancing dispersion. Unlike bulky, rigid molecules that remain trapped, these components remain mobile and repeatedly interact with multiple aggregation sites. Experimental validation using rheometry and FTIR confirmed that such molecularly compatible rejuvenators restore the binder flexibility and polydispersity, even when diffusion is relatively slow. Building on this mechanistic foundation, six rejuvenators (A2, A5, A7, A8, A9, and A10) were evaluated using a three-metric performance framework encompassing cracking resistance (Glover–Rowe parameter), UV stability, and surface hydrophobicity. Among the six rejuvenators evaluated, A9 exhibited the highest overall performance, reducing the Glover–Rowe cracking parameter by 85% (from 351 kPa to 53 kPa), demonstrating the greatest resistance to UV-induced aging with a stability index of 4.17 h·kPa⁻¹, and increasing surface hydrophobicity to a contact angle of 103.6°. Its superior performance is primarily attributed to its amide- and unsaturated-chain components, which act cooperatively to disrupt π–π stacking interactions and exfoliate asphaltene layers, thereby promoting molecular deagglomeration and enhancing long-term durability. These results shift the criteria for selecting rejuvenators: effective candidates must pair the electronic capability to unlock aged asphaltenes with sufficient structural stability to resist secondary aging and restore hydrophobicity. Restoring hydrophobicity is critical, as aging reduces the asphalt’s water repellency and increases water diffusion, which in turn accelerates moisture-related damage; consequently, a high-performing rejuvenator must effectively restore the binder’s hydrophobic characteristic. Collectively, these findings provide a framework for the rational design of next-generation of bio-based rejuvenators that enhance pavement longevity and promote long-term sustainability.

Aging↗
DOE OSTI2025

Encapsulation of Ba in InSb Framework Introduces Chirality in Clathrate-Like BaIn 4 Sb 4

The discovery of Zintl compounds remains a powerful strategy for identifying materials with tunable electronic and thermal transport properties. During a concerted search for new inorganic clathrates with In–Sb frameworks, we discovered BaIn 4 Sb 4 . The composition of this phase deviates from that expected for a type-I clathrate with tetrahedral coordination of all In and Sb atoms (Ba 8 In 31 Sb 15 ). Instead, in the chiral structure of BaIn 4 Sb 4 (space group P 3 1 21, No. 152), a part of the In atoms have a trigonal planar coordination of 1In + 2Sb, forming Sb 2 –In–In–Sb 2 nonplanar fragments isostructural to diborane(4) B 2 H 4 with D 2 d symmetry. Ba atoms are located inside 16-vertex In 8 Sb 8 polyhedra, which share vertices and edges to form a chiral framework around the 3 1 screw axes. The title compound is electron-balanced, [Ba 2+ ][In 2+ ] 2 [In 3+ ] 2 [Sb3–] 4 , which was confirmed by characterization of the charge and heat transport properties. BaIn 4 Sb 4 exhibits a low thermal conductivity and high Seebeck coefficient, suggesting its untapped potential for thermoelectric applications. Finally, density functional theory (DFT) calculations indicate that chemical doping may enhance carrier concentration and improve the originally low electrical conductivity, thus enhancing thermoelectric performance.

electrical properties↗
DOE OSTI2021

Teaching a neural network to attach and detach electrons from molecules

Abstract Interatomic potentials derived with Machine Learning algorithms such as Deep-Neural Networks (DNNs), achieve the accuracy of high-fidelity quantum mechanical (QM) methods in areas traditionally dominated by empirical force fields and allow performing massive simulations. Most DNN potentials were parametrized for neutral molecules or closed-shell ions due to architectural limitations. In this work, we propose an improved machine learning framework for simulating open-shell anions and cations. We introduce the AIMNet-NSE (Neural Spin Equilibration) architecture, which can predict molecular energies for an arbitrary combination of molecular charge and spin multiplicity with errors of about 2–3 kcal/mol and spin-charges with error errors ~0.01e for small and medium-sized organic molecules, compared to the reference QM simulations. The AIMNet-NSE model allows to fully bypass QM calculations and derive the ionization potential, electron affinity, and conceptual Density Functional Theory quantities like electronegativity, hardness, and condensed Fukui functions. We show that these descriptors, along with learned atomic representations, could be used to model chemical reactivity through an example of regioselectivity in electrophilic aromatic substitution reactions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Efficient simulations of charge density waves in the transition metal Dichalcogenide TiSe 2

Charge density waves (CDWs) in transition metal dichalcogenides are the subject of growing scientific interest due to their rich interplay with exotic phases of matter and their potential technological applications. Here, using density functional theory with advanced meta-generalized gradient approximations (meta-GGAs) and linear response time-dependent density functional theory (TDDFT) with state-of-the-art exchange-correlation kernels, we investigate the electronic, vibrational, and optical properties in 1T-TiSe 2 with and without CDW. In both bulk and monolayer TiSe 2 , the electronic bands and phonon dispersions in either normal or CDW (semiconducting) phase are described well via meta-GGAs, which separate the valence and conduction bands just as HSE06 does but with significantly more computational feasibility. The experimentally observed humps of electron energy loss spectroscopy are successfully reproduced in TDDFT. Our work opens the door to simulating these complexities in CDW compounds from first principles by revealing meta-GGAs as an accurate low-cost alternative to HSE06.

Electronic properties and materials↗
DOE OSTI2021

Chemisorbed oxygen or surface oxides steer the selectivity in Pd electrocatalytic propene oxidation observed by operando Pd L-edge X-ray absorption spectroscopy

Controlled electrochemical oxidation of hydrocarbons to desired products is an attractive approach in catalysis. Here we study the electrochemical propene oxidation under operando conditions using Pd L-edge X-ray absorption spectroscopy (XAS) as a sensitive probe to elucidate surface processes occurring during catalysis. Together with ab initio multiple-scattering calculations, our XAS results enable assignment of characteristic changes of the Pd L-edge intensity and energy position in terms of a mechanistic understanding of the selective oxidation of propene. The results, supported by electrochemical density functional theory DFT simulations, show that in the potential range of 0.8–1.0 V vs. the reversible hydrogen electrode (RHE), selective oxidation of propene to acrolein and acrylic acid occurs on the metallic Pd surface. Here, these reactions are proposed to proceed via the Langmuir–Hinshelwood mechanism. In contrast, for the potential range of 1.1–1.3 V vs. RHE, selective oxidation of propene to propylene glycol takes place on a Pd oxide surface.

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

Impact of secondary salts, temperature, and pH on the colloidal stability of graphene oxide in water

The stability of graphene oxide (GO) in water is extremely relevant because of its application as an adsorbent material, as well as for its fate and behavior in the environment. Zeta potential was used to study the effect of secondary salts (carbonate, sulfate, and phosphate), temperature (20 to 60 °C), and pH (5 to 9) on the stability of six different GOs produced from natural, synthetic, and amorphous graphite—with and without the use of attrition milling. Generally, GOs produced with attrition-milled graphites had lower ζ-potentials than their unmilled counterparts because of their smaller particle sizes and higher concentration of oxygen-containing functional groups. It was observed that GO produced from graphite and synthetic graphite had ζ-potential values lower than –30 mV, even at 30 °C. However, it was observed that all the GOs studied were unstable in the presence of carbonate and sulfate salts at concentrations between 170 and 1695 mg L –1 , as they reached a ζ-potential of –4.1 mV. Density-functional theory electronic structure calculations suggested that the instability of GO in the presence of carbonate and sulfate was caused by the abstraction of a proton resulting in interaction energies E int of 28.3 and 168.9 kJ mol –1 , respectively. Our results suggest that temperatures above 30 °C, as well as carbonate and sulfate salts at concentrations relevant to arid and semi-arid regions, could promote the formation of agglomerates of GO, thus limiting its use and mobility in water.

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

Ab initio generalized Langevin equation

We introduce a machine learning–based approach called ab initio generalized Langevin equation (AIGLE) to model the dynamics of slow collective variables (CVs) in materials and molecules. In this scheme, the parameters are learned from atomistic simulations based on ab initio quantum mechanical models. Force field, memory kernel, and noise generator are constructed in the context of the Mori–Zwanzig formalism, under the constraint of the fluctuation–dissipation theorem. Combined with deep potential molecular dynamics and electronic density functional theory, this approach opens the way to multiscale modeling in a variety of situations. Here, we demonstrate this capability with a study of two mesoscale processes in crystalline lead titanate, namely the field-driven dynamics of a planar ferroelectric domain wall, and the dynamics of an extensive lattice of coarse-grained electric dipoles. In the first case, AIGLE extends the reach of ab initio simulations to a regime of noise-driven motions not accessible to molecular dynamics. In the second case, AIGLE deals with an extensive set of CVs by adopting a local approximation for the memory kernel and retaining only short-range noise correlations. The scheme is computationally more efficient than molecular dynamics by several orders of magnitude and mimics the microscopic dynamics at low frequencies where it reproduces accurately the dominant far-infrared absorption frequency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Bound states of Ω baryons in light nuclei

Here, we investigate bound states of light Ω 3⁢𝑥 clusters (𝑥=𝑠,𝑐), motivated by the Ω 3⁢𝑠 ⁢𝑁 potential recently developed by the HAL QCD collaboration. To regularize this potential, we remove the deeply attractive core at 𝑟 < 0.4 fm and parametrize the long-range component (𝑟 > 0.4 fm) using a two-range Gaussian form. This procedure preserves the relevant two-body bound-state energy while having a negligible effect on the Ω 3⁢𝑠⁢ 𝑁⁢𝑁 and Ω 3⁢𝑠⁢ Ω 3⁢𝑠 ⁢𝑁 systems. An effective Ω 3⁢𝑠 ⁢𝛼 potential is then constructed by fitting a two-range Gaussian function to the long-range component of the folding potential, enabling calculations of the bound-state energies of the Ω 3⁢𝑠⁢ 𝛼, Ω 3⁢𝑠⁢ 𝛼⁢𝛼, and Ω 3⁢𝑠 ⁢Ω 3⁢𝑠 ⁢𝛼 systems. The regularization procedure leads to a substantial reduction in bound-state energies compared to those obtained with the original potential. We further extend the analysis to Ω 3⁢𝑐 -cluster systems by introducing an Ω 3⁢𝑐 ⁢𝑁 interaction, derived by comparing the existing Ω 3⁢𝑠⁢ Ω 3⁢𝑠 and Ω 3⁢𝑐 ⁢Ω 3⁢𝑐 potentials. Our results suggest that several parametrizations predict bound states in Ω 3⁢𝑐 -containing clusters. Finally, the Ω 3⁢𝑠 ⁢Ω 3⁢𝑠 interaction is described using a contactlike potential approach, motivated by the effective field theory.

binding energy & masses↗