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 667 records · Page 37

Contribution of three-body potentials to the binding energy of heavy atoms

The conversion of quantum electrodynamics to a configuration-space Hamiltonian formalism introduces three-electron potentials of relativistic origin. For heavy atoms, it is found that the contribution of these potentials to the inner-shell binding energy is no more than 0.21 eV. This is too small to explain the discrepancy between current theory and experiment. The uniqueness of the potentials obtained in the configuration-space Hamiltonian is also discussed.

Zygelman, B.↗

All-electron molecular Dirac-Hartree-Fock calculations - Properties of the XH4 and XH2 molecules and the reaction energy XH4 yields XH2 + H2, X = Si, Ge, Sn, Pb

Relativistic corrections to a number of properties of the Group IV hydrides are calculated using the Dirac-Hartree-Fock method. The use of first-order perturbation theory is sufficient to obtain relativistic corrections for Ge, but the effects of spin-orbit interaction and other higher-order effects begin to show for Sn and become important for Pb. The energy of the reaction XH4 yields XH2 + H2 (X = Si, Ge, Sn, and Pb) is also calculated. The results are compared with relativistic effective core potential calculations, first-order perturbation theory calculations, and limited experimental data.

Dyall, Kenneth G.↗

Fluorination Enables Simultaneous Improvements of a Dialkoxybenzene-Based Redoxmer for Nonaqueous Redox Flow Batteries

Redoxmers or redox-active organic materials, arc one critical component for nonaqueous redox flow batteries (RFBs), which hold high promise in enabling the time domain of the grid. While tuning redox potentials of redoxmers is a very effective way to enhance energy densities of NRFBs, those improvements often accompany accelerated kinetics of the charged species, undermining stability and cycling performance. In this paper, a strategy for designing redoxmers with simultaneous improvements in redox potential and stability is proposed. Specifically, the redoxmer 1,4-di-tert-butyl-2,5-bis(2,2,2-trifluoroethoxy)benzene (ANL-C46) is developed by incorporating fluorinated substitutions into the dialkoxybenzene-based platform. Compared to the non-fluorinated analogue, ANL-C46 demonstrates not only an increased (similar to 0.41 V) redox potential but also much enhanced stability (1.6 times) and cyclability (4 times) evidenced by electron paramagnetic resonance kinetic study, H-cell and flow cell cycling. In fact, the cycling performance of ANL-C46 is among the best of high potential (>1.0 V vs Ag/Ag + ) redoxmers ever reported. Density functional theory calculations suggest that while the introduced fluorine substitutions elevate the redox potentials, they also help to depress the decomposition reactions of the charged redoxmers, affording excellent stability. The findings represent an interesting strategy for simultaneously improving energy density and stability, which could further prompt the development of high-performance redoxmers.

25 ENERGY STORAGE↗

A Van der Waals-like theory of plasma double layers

A theory describing plasma double layers in terms of multiple roots of the charge density expression is presented. The theory presented uses the fact that equilibrium plasmas shield small potential perturbations linearly; for high potentials, the shielding decreases. The approach is analogous to Van der Waals' theory of simple fluids in which inclusion of approximate expressions for both excluded volume and long range attractive forces sufficiently describes the first-order liquid-gas phase transition.

Katz, Ira↗

Resolving Discrepancies Between Theory and Experiment for the NCN + H Reaction

Flame modeling studies have highlighted the role of branching in the kinetics for the prompt NO switch reaction, NCN + H. In this reaction, there is strong competition between the CH + N2 (R1) and HCN + N (R2) product channels. Increased branching towards reaction R2 promotes the subsequent formation of Fenimore or Prompt NO. Recent direct shock tube studies on this complex reaction conclude that branching to reaction R2 is much more favorable than the predicted branching from theoretical studies. The experimentally predicted prompt NO switch temperature (TS), i.e., the temperature at which k2/(k1+k2) = 0.5, is at 1670 K, in contrast to theory that predicts Ts > 3200 K. The present theory/modeling work was initiated to highlight potential causes for this significant discrepancy between experiment and theory. Our analysis of the simulations reported in the experimental studies indicate that including a fall-off representation for C2H5 radical dissociation (with C2H5 produced from C2H5I dissociation and acting as an H-atom source) has a noticeable influence on the simulations for NCN decay and HCN formation in the shock tube studies. With the inclusion of a theory-based pressure fall-off representation for C2H5 dissociation, we show that this radical persists for much longer timescales even at the high-temperatures in the shock tube studies. This persistence then facilitates the rapid conversion of reactive H-atoms to CH3 radicals at the shock tube conditions via the H + C2H5 → CH3 + CH3 reaction. An updated theoretical analysis for this classic addition-elimination reaction is also provided in this work. The rapid formation of CH3 radicals necessitates the inclusion of additional reactions to adequately simulate the NCN and HCN data. Our simulations conclude that, with the inclusion of these additional reactions, the NCN and HCN profiles can be reasonably well simulated by the most recent theoretical predictions.

Ethyl Radicals↗

Fluorine spillover for ceria- vs silica-supported palladium nanoparticles: A MD study using machine learning potentials

Supported metallic nanoparticles play a central role in catalysis. However, predictive modeling is particularly challenging due to the structural and dynamic complexity of the nanoparticle and its interface with the support, given that the sizes of interest are often well beyond those accessible via traditional ab initio methods. With recent advances in machine learning, it is now feasible to perform MD simulations with potentials retaining near-density-functional theory (DFT) accuracy, which can elucidate the growth and relaxation of supported metal nanoparticles, as well as reactions on those catalysts, at temperatures and time scales approaching those relevant to experiments. Furthermore, the surfaces of the support materials can also be modeled realistically through simulated annealing to include effects such as defects and amorphous structures. We study the adsorption of fluorine atoms on ceria and silica supported palladium nanoparticles using machine learning potential trained by DFT data using the DeePMD framework. We show defects on ceria and Pd/ceria interfaces are crucial for the initial adsorption of fluorine, while the interplay between Pd and ceria and the reverse oxygen migration from ceria to Pd control spillover of fluorine from Pd to ceria at later stages. In contrast, silica supports do not induce fluorine spillover from Pd particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Weighted active space protocol for multireference machine-learned potentials

Multireference methods such as multiconfiguration pair-density functional theory accurately capture electronic correlation in systems with strong multiconfigurational character, but their cost precludes direct use in molecular dynamics. Combining these methods with machine-learned interatomic potentials (MLPs) can extend their reach. However, the sensitivity of multireference calculations to the choice of the active space complicates the consistent evaluation of energies and gradients across structurally diverse nuclear configurations. To overcome this limitation, we introduce the weighted active space protocol (WASP), a systematic approach to assign a consistent active space for a given system across uncorrelated configurations. By integrating WASP with MLPs and enhanced sampling techniques, we propose a data-efficient active learning cycle that enables the training of an MLP on multireference data. We demonstrated the approach on the TiC + -catalyzed C–H activation of methane, a reaction that poses challenges for Kohn–Sham density functional theory due to its significant multireference character. This framework enables accurate and efficient modeling of catalytic dynamics, establishing a paradigm for simulating complex reactive processes beyond the limits of conventional electronic-structure methods.

enhanced sampling↗

An Atomistic Study of Reactivity in Solid-State Electrolyte Interphase Formation for Li/Li7P3S11

Lithium metal batteries offer superior volumetric and gravimetric specific capacities compared to those based on traditional graphite anodes. Although advancements in solid-state electrolytes address safety concerns, challenges remain, particularly regarding interphase formation in lithium metal anodes. This work presents a computational framework based on high-throughput first-principles density functional theory and machine-learning interatomic potentials (MLIPs) including automated iterative, active learning to enable robust computational exploration of interphase formation between lithium metal anodes and an inorganic solid-state electrolyte. As a demonstration, we apply the framework to a Li/Li7P3S11 interface and find that it accurately identifies the experimentally observed, thermodynamically stable interphase products as well as their overall spatial arrangement within a heterogeneous, amorphous layered structure, with Li2S domains of nanocrystallinity. Our simulations show two stages, a fast and slow diffusion reaction regime, that corroborate the relative phase formation rate of Li x P, Li2S, and Li3P. Using the Onsager transport theory, we capture time-dependent ionic diffusion within the reacting interface, including cross-correlation effects. We found that cross-correlation effects between Li-P and P-S ionic motion significantly influence P-ion diffusion, making it highly sensitive to the local environment and potentially leading to "kinetic trapping" of Li-P phases. The passivation of the interface is shown as the ionic fluxes all approach zero, effectively halting interphase growth.

Diffusion↗

Unified Aeroacoustics Analysis for High Speed Turboprop Aerodynamics and Noise: Development of Theory for Blade Loading, Wakes, and Noise - Volume 1

A unified theory for the aerodynamics and noise of advanced turboprops are presented. Aerodynamic topics include calculation of performance, blade load distribution, and non-uniform wake flow fields. Blade loading can be steady or unsteady due to fixed distortion, counter-rotating wakes, or blade vibration. The aerodynamic theory is based on the pressure potential method and is therefore basically linear. However, nonlinear effects associated with finite axial induction and blade vortex flow are included via approximate methods. Acoustic topics include radiation of noise caused by blade thickness, steady loading (including vortex lift), and unsteady loading. Shielding of the fuselage by its boundary layer and the wing are treated in separate analyses that are compatible but not integrated with the aeroacoustic theory for rotating blades.

Hanson, D. B.↗

Exchange correlation potentials from full configuration interaction in a Slater orbital basis

Ryabinkin–Kohut–Staroverov (RKS) theory builds a bridge between wave function theory and density functional theory by using quantities from the former to produce accurate exchange-correlation potentials needed by the latter. In this work, the RKS method is developed and tested alongside Slater atomic orbital basis functions for the first time. To evaluate this approach, full configuration interaction computations in the Slater orbital basis are employed to give quality input to RKS, allowing full correlation to be present along with correct nuclei cusps and asymptotic decay of the wavefunction. SlaterRKS is shown to be an efficient algorithm to arrive at exchange-correlation potentials without unphysical artifacts in moderately-sized basis sets. Furthermore, enforcement of the nuclear cusp conditions will be shown to be vital for the success of the Slater-basis RKS method. As a result, examples of weakly and strongly correlated molecular systems will demonstrate the main features of SlaterRKS.

74 ATOMIC AND MOLECULAR PHYSICS↗

A Theory-Guided X-ray Absorption Spectroscopy Approach for Identifying Active Sites in Atomically Dispersed Transition-Metal Catalysts

Atomically dispersed supported metal catalysts offer new properties and the benefits of maximized metal accessibility and utilization. The characterization of these materials, however, remains challenging. Using atomically dispersed platinum supported on crystalline MgO (chosen for its well-defined bonding sites) as a prototypical example, we demonstrate in this work how systematic density functional theory calculations for assessing all the potentially stable platinum sites, combined with automated analysis of extended X-ray absorption fine structure (EXAFS) spectra, leads to unbiased identification of isolated, surface-enveloped platinum cations as the catalytic species for CO oxidation. The catalyst has been characterized by atomic-resolution imaging and EXAFS and high-energy resolution fluorescence detection X-ray absorption near edge spectroscopy. The proposed platinum sites are in agreement with experiment. This theory-guided workflow leads to rigorously determined structural models and provides a more detailed picture of the structure of the catalytically active site than what is currently possible with conventional EXAFS analyses. As this approach is efficient and agnostic to the metal, support, and catalytic reaction, we posit that it will be of broad interest to the materials characterization and catalysis communities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Band gap and electronic structure of defects in the ternary nitride BP 3 N 6 : experiment and theory

Recent advances in methods to access nitride systems by a high-pressure high-temperature approach have made possible the one-step synthesis of mixed ternary non-metal nitrides. As a prerequisite to use in a practical device, it is important to understand important bulk electronic properties, such as the band gap, as well as characterizing the presence and effect of defects that are present. In this work, the novel ternary nitride BP 3 N 6 is studied using techniques sensitive to the partial electronic density of states, specifically X-ray absorption spectroscopy and X-ray emission spectroscopy. Complementary full-potential all-electron density functional theory (DFT) calculations allow important bulk electronic parameters, such as the band gap, to be elucidated. The band gap of BP 3 N 6 has been determined to be 3.9 ± 0.2 eV and 4.1 ± 0.4 eV at the B K- and N K-edges, respectively. This is close to a theoretical value of 4.3 eV predicted by the PBEsol exchange–correlation functional and considerably less than a value of 5.8 eV predicted by the modified Becke–Johnson exchange–correlation functional. X-Ray excited optical luminescence (XEOL) measurements are performed to interrogate the presence of point defects in this system. Together with DFT calculations, these measurements reveal the presence of nitrogen vacancies which lead to multiple mid-gap trap states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gold–Thiolate Nanocluster Dynamics and Intercluster Reactions Enabled by a Machine Learned Interatomic Potential

Monolayer protected metal clusters comprise a rich class of molecular systems and are promising candidate materials for a variety of applications. While a growing number of protected nanoclusters have been synthesized and characterized in crystalline forms, their dynamical behavior in solution, including prenucleation cluster formation, is not well understood due to limitations both in characterization and first-principles modeling techniques. Recent advancements in machine-learned interatomic potentials are rapidly enabling the study of complex interactions such as dynamical behavior and reactivity on the nanoscale. Here, we develop an Au-S-C-H atomic cluster expansion (ACE) interatomic potential for efficient and accurate molecular dynamics simulations of thiolate-protected gold nanoclusters (Au n (SCH 3 ) m ). Trained on more than 30,000 density functional theory calculations of gold nanoclusters, the interatomic potential exhibits ab initio level accuracy in energies and forces and replicates nanocluster dynamics including thermal vibration and chiral inversion. Long dynamics simulations (up to 0.1 μs time scale) reveal a mechanism explaining the thermal instability of neutral Au 25 (SR) 18 clusters. Specifically, we observe multiple stages of isomerization of the Au 25 (SR) 18 cluster, including a chiral isomer. Additionally, we simulate coalescence of two Au 25 (SR) 18 clusters and observe series of clusters where the formation mechanisms are critically mediated by ligand exchange in the form of [Au-S] n rings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Probing the self-ionization of liquid water with ab initio deep potential molecular dynamics

The chemical equilibrium between self-ionized and molecular water dictates the acid–base chemistry in aqueous solutions, yet understanding the microscopic mechanisms of water self-ionization remains experimentally and computationally challenging. Herein, Density Functional Theory (DFT)–based deep neural network (DNN) potentials are combined with enhanced sampling techniques and a global acid–base collective variable to perform extensive atomistic simulations of water self-ionization for model systems of increasing size. The explicit inclusion of long-range electrostatic interactions in the DNN potential is found to be crucial to accurately reproduce the DFT free energy profile of solvated water ion pairs in small (64 and 128 H 2 O) cells. The reversible work to separate the hydroxide and hydronium to a distance S is found to converge for simulation cells containing more than 500 H 2 O, and a distance of ∼ 8 Å is the threshold beyond which the work to further separate the two ions becomes approximately zero. The slow convergence of the potential of mean force with system size is related to a restructuring of water and an increase of the local order around the water ions. Calculation of the dissociation equilibrium constant illustrates the key role of long-range electrostatics and entropic effects in the water autoionization process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessment of empirical interatomic potential to predict thermal conductivity in ThO 2 and UO 2

Computing vibrational properties of crystals in the presence of complex defects often necessitates the use of (semi-)empirical potentials, which are typically not well characterized for perfect crystals. In this study we explore the efficacy of a commonly used embedded-atomempirical interatomic potential for the U x Th 1- x O 2 system, to compute phonon dispersion, lifetime, and branch specific thermal conductivity. Our approach for ThO 2 involves using lattice dynamics and the linearized Boltzmann transport equation to calculate phonon transport properties based on second and third order force constants derived from the empirical potential and from first-principles calculations. For UO 2 , to circumvent the accuracy issues associated with first-principles treatments of strong electronic correlations, we compare results derived from the empirical interatomic potential to previous experimental results. It is found that the empirical potential can reasonably capture the dispersion of acoustic branches, but exhibits significant discrepancies for the optical branches, leading to overestimation of phonon lifetime and thermal conductivity. The branch specific conductivity also differs significantly with either first-principles based results (ThO 2 ) or experimental measurements (UO 2 ). These findings suggest that the empirical potential needs to be further optimized for robust prediction of thermal conductivity both in perfect crystals and in the presence of complex defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Machine-learned quantum molecular dynamics calculations of warm dense equation of state and ionic transport coefficients of deuterated water

White dwarf models require accurate equations of state and ionic transport coefficients in the warm dense matter regime, where kinetic theory models and tabulated equations of state are often inaccurate. In this work, spectral-partitioned density functional theory and machine-learned interatomic potentials are combined to perform large-scale, first-principles quantum molecular dynamics simulations of deuterated water (D 2 O) near the principal Hugoniot. This approach retains Kohn-Sham accuracy while achieving orders-of-magnitude speedup, yielding converged equation of state and transport properties over a broad pressure and temperature range. The results reveal the thermodynamic conditions under which ionic transport models for interdiffusivity and shear viscosity converge and identify those in closest agreement with density functional theory benchmarks at temperatures in the warm dense matter regime. The present framework extends first-principles transport calculations to higher temperatures than previously achieved, and provides an efficient, scalable, and general approach for studying transport properties in complex multicomponent mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗

Experimental tests of relativistic gravitation theories

Experimental tests were studied for determining the potential uses of future deep space missions in studies of relativistic gravity. The extensions to the parametrized post-Newtonian framework to take explicit account of the solar system's center of mass relative to the mean rest frame of the Universe is reported. Discoveries reported include the Machian effects of motion relative to the universal rest frame. Summaries of the JPL research are included.

Anderson, J. D.↗

A new approach for solving the three-dimensional steady Euler equations. I - General theory

The present iterative procedure combines the Clebsch potentials and the Munk-Prim (1947) substitution principle with an extension of a semidirect Cauchy-Riemann solver to three dimensions, in order to solve steady, inviscid three-dimensional rotational flow problems in either subsonic or incompressible flow regimes. This solution procedure can be used, upon discretization, to obtain inviscid subsonic flow solutions in a 180-deg turning channel. In addition to accurately predicting the behavior of weak secondary flows, the algorithm can generate solutions for strong secondary flows and will yield acceptable flow solutions after only 10-20 outer loop iterations.

Chang, S.-C.↗