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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

First-Principles Study of Tritium Trapping by Point Defects in Fe-Al Aluminide Coating Phases

Density functional theory simulations have been carried out to investigate the potential for tritium trapping by metal vacancies in five different Fe-Al aluminide coating phases. It was found that tritiation of Fe and Ni vacancies is generally less favorable than the tritiation of Al vacancies. However, for the first tritiation, a trend in the defect formation energy can be obtained such that metal defects in the Fe 2 Al x family of materials (i.e., Fe 2 Al 4 , Fe 2 Al 5 , and Fe 2 Al 6 ) trap tritium species more favorably than metal vacancies in Fe 4 Al 13 and FeNiAl 5 . Further investigations using ab initio thermodynamics calculations confirmed that trend for a range of tritium partial pressure at a temperature of 700 K. Especially, it was shown that the energy difference between tritiated and non-tritiated metal vacancies is smaller and more favorable for the Fe 2 Al x family, followed by Fe 4 Al 13 , and FeNiAl 5 . While this study shows that tritium interacts differently in the various Fe-Al aluminide phases, it also suggests that tritium trapping and retention could be more efficient if metal defects are present in some Fe-Al phases.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thermodynamics of Tritium Trapping by Point Defects in Intermetallic Al 12 (TM) 2.35 Aluminide Coating Phase

Density functional theory simulations have been carried out to investigate the potential for tritium trapping by metal vacancies in intermetallic Al 12 (TM) 2.35 phase (TM = Fe, Cr, and Ni) as function of temperature and tritium partial pressure. It was found that tritium could be favorably trapped by Fe and Ni vacancies and not favorably trapped by Al and Cr vacancies. However, due to the presence of partially occupied Al sites in bulk Al 12 (TM) 2.35 , leading to the approximate number of ~255 Al atoms in the unit cell, 86 sites were found energetically favorable to the creation of an Al vacancy. While adding a tritium atom in an Al vacancy is not energetically favorable, the tritiated defect still has a negative Gibbs free energy because the energy gain for creating an Al vacancy overcome the energy cost of adding the tritium species. Based on the calculated Gibbs free energy, the first tritiation of a metal vacancy, at conditions relevant to in-reactor operations, should be more favorable for Al, followed Fe, Ni, and Cr vacancies. By comparing the behavior of tritium in Al 12 (TM) 2.35 with previously studied Fe-Al coating phases (i.e., FeNiAl 5 , Fe 4 Al 13 , and Fe 2 Al 5.6 ), we found that there is a correlation between interstitial tritium solubility and the potential for vacancy trapping. The current trend suggests that if the insertion of an interstitial tritium cost more than 0.3 eV, then trapping by metal vacancies should be preferred. By combining the simulations results obtained to date, we noticed different trapping mechanisms of tritium in the Al coating. Tritium is mostly trapped by Fe and Ni vacancies in the outer Fe-Al coating phase Al 12 (TM) 2.35 while tritium should be preferentially trapped by Al and Fe vacancies for the inner Fe-Al coating phases (FeNiAl 5 , Fe 4 Al 13 , Fe 2 Al 5.6 ). Altogether, these studies show that tritium interacts differently with the various Fe-Al aluminide phases, they also suggest that tritium trapping and retention could be more efficient if metal defects are present and if the solubility of interstitial tritium in the different phases is low.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Characterization of dissipative regions of a N-doped superconducting radio-frequency cavity

We report radio-frequency measurements of quality factors and temperature mapping of a nitrogen doped Nb superconducting RF cavity. Cavity cutouts of hot and cold spots were studied with low temperature scanning tunneling microscopy and spectroscopy, X-ray photoelectron spectroscopy and secondary electron microscopy. Temperature mapping revealed a substantial reduction of the residual resistance upon cooling the cavity with a greater temperature gradient and hysteretic losses at the quench location, pointing to trapped vortices as the dominant source of residual surface resistance. Analysis of the tunneling spectra in the framework of a proximity effect theory shows that hot spots have a reduced pair potential and a wider distribution of the contact resistance between the Nb and the top Nb oxide. Alone, these degraded superconducting properties account for a much weaker excess dissipation as compared with the vortex contribution. Based on the correlation between the quasiparticle density of states and temperature mapping, we suggest that degraded superconducting properties may facilitate vortex nucleation or settling of trapped flux during cooling the cavity through the critical temperature.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

A Comparison of Three Ways to Measure Time-Dependent Densities With Quantum Simulators

Quantum algorithms are touted as a way around some classically intractable problems such as the simulation of quantum mechanics. At the end of all quantum algorithms is a quantum measurement whereby classical data is extracted and utilized. In fact, many of the modern hybrid-classical approaches are essentially quantum measurements of states with short quantum circuit descriptions. Here, we compare and examine three methods of extracting the time-dependent one-particle probability density from a quantum simulation: direct Z -measurement, Bayesian phase estimation, and harmonic inversion. We have tested these methods in the context of the potential inversion problem of time-dependent density functional theory. Our test results suggest that direct measurement is the preferable method. We also highlight areas where the other two methods may be useful and report on tests using Rigetti's quantum virtual device. This study provides a starting point for imminent applications of quantum computing.

Yang, Jun↗
DOE OSTI2026

The Molecular Structures of Liquid and Glassy Nifedipine and Felodipine and Their Incorporation into PVP

Background: Amorphous drug formulations are commonly used to improve the solubility and bioavailability of poorly soluble molecular pharmaceuticals, yet less is known about their molecular conformations and local bonding interactions than their crystalline phases. Methods: High-energy X-ray diffraction structure factor measurements have been made on liquid and glassy nifedipine (NIF), felodipine (FEL), NIF 1:3 polyvinylpyrrolidone (PVP), and FEL 1:3 PVP wt.% mixtures. The corresponding X-ray pair distribution functions have been interpreted using empirical potential structure refinement using different models and density functional theory conformer calculations. Results: In both NIF and FEL, the NH···O inter-molecular hydrogen bonds between the pyridyl nitrogen and ester carbonyls are found to be considerably weaker than those observed in the crystalline polymorphs. For nifedipine, it is proposed that either inter-molecular NH…ON nitro bonds are present and/or a fraction (<20%) of conformational changes, with the aryl ring flipped, occur in the liquid state. For felodipine, the models indicate significant disorder associated with the methyl and ethyl side chains in the liquid state, with the main peak intensity at 3.0 Å arising from intra-molecular Cl-Cl atom pairs. When nifedipine molecules are incorporated into PVP, our models show they possess stronger NH···O bonds to the PVP polymer than felodipine molecules, which have stronger affinity for bonding to the polymer than to other felodipine molecules. Conclusions: The amorphous forms of both NIF and FEL show much weaker hydrogen bonding than found in their crystalline phases. Liquid NIF also exhibits configurations which are not observed in the crystal phases.

Amorphous Solid Dispersion↗
DOE OSTI2025

Bottomonium suppression in pNRQCD and open quantum system approach

By employing the potential non-relativistic quantum chromodynamics (pNRQCD) effective field theory within an open quantum system framework, we derive a Lindblad equation governing the evolution of the heavy-quarkonium reduced density matrix, accurate to next-to-leading order (NLO) in the ratio of the state's binding energy to the medium's temperature [1]. The derived NLO Lindblad equation provides a more reliable description of heavy-quarkonium evolution in the quark-gluon plasma at low temperatures compared to the leading-order truncation. For phenomenological applications, we numerically solve this equation using the quantum trajectories algorithm. By averaging over Monte Carlo-sampled quantum jumps, we obtain solutions without truncation in the angular momentum quantum number of the considered states. Our analysis highlights the importance of quantum jumps in the nonequilibrium evolution of bottomonium states within the quark-gluon plasma [2]. Additionally, we demonstrate that the quantum regeneration of singlet states from octet configurations is essential to explain experimental observations of bottomonium suppression. The heavy-quarkonium transport coefficients used in our study align with recent lattice QCD determinations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Acid anion electrolyte effects on platinum for oxygen and hydrogen electrocatalysis

Platinum is an important material with applications in oxygen and hydrogen electrocatalysis. To better understand how its activity can be modulated through electrolyte effects in the double layer microenvironment, herein we investigate the effects of different acid anions on platinum for the oxygen reduction/evolution reaction (ORR/OER) and hydrogen evolution/oxidation reaction (HER/HOR) in pH 1 electrolytes. Experimentally, we see the ORR activity trend of HClO 4 > HNO 3 > H 2 SO 4 , and the OER activity trend of HClO 4 > HNO 3 ~ H 2 SO 4 . HER/HOR performance is similar across all three electrolytes. Notably, we demonstrate that ORR performance can be improved 4-fold in nitric acid compared to in sulfuric acid. Assessing the potential-dependent role of relative anion competitive adsorption with density functional theory, we calculate unfavorable adsorption on Pt(111) for all the anions at HER/HOR conditions while under ORR/OER conditions ClO$^{-}_{4}$ binds the weakest followed by NO$^{-}_{3}$ and SO$^{2-}_{4}$. Our combined experimental-theoretical work highlights the importance of understanding the role of anions across a large potential range and reveals nitrate-like electrolyte microenvironments as interesting possible sulfonate alternatives to mitigate the catalyst poisoning effects of polymer membranes/ionomers in electrochemical systems. These findings help inform rational design approaches to further enhance catalyst activity via microenvironment engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Real-Time Time-Dependent Electronic Structure Theory

Real-time electronic structure methods provide an unprecedented view of electron dynamics and ultrafast spectroscopy on the atto- and femtosecond time scale with vast potential to yield new insights into the electronic behavior of molecules and materials. Here, we discuss the fundamental theory underlying various real-time electronic structure methods as well as advantages and disadvantages of each. We give an overview of the numerical techniques that are widely used for real-time propagation of the quantum electron dynamics with an emphasis on Gaussian basis set methods. We also showcase many of the chemical applications and scientific advances made by using real-time electronic structure calculations and provide an outlook of possible new directions.

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

Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces

Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.

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

Solving the one-dimensional penetration problem for the fission channel in the statistical Hauser-Feshbach theory

Here we solve the Schrödinger equation for an arbitrary one-dimensional potential energy to calculate the transmission coefficient in the fission channel of compound nucleus reactions. We incorporate the calculated transmission coefficients into the statistical Hauser-Feshbach model calculation for neutron-induced reactions on 235,238 U and 239 Pu. The one-dimensional model reproduces the evaluated fission cross section data reasonably well considering the limited number of model parameters involved. A resonance-like structure appears in the transmission coefficient for a double-humped fission barrier shape that includes an intermediate well, which is understood to be a quantum mechanical effect in the fission channel. The calculated fission cross sections for the neutron-induced reactions on 235,238 U and 239 Pu all exhibit a similar structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Learning local and semi-local density functionals from exact exchange-correlation potentials and energies

Finding accurate exchange-correlation (XC) functionals remains the defining challenge in density functional theory (DFT). Despite 40 years of active development, attaining general purpose chemical accuracy is still elusive with existing functionals. We present a data-driven pathway to learn the XC functional by using the exact density, XC energy, and XC potential. While the exact densities are obtained from accurate configuration interaction (CI), the exact XC energies and XC potentials are obtained via inverse DFT calculations on the CI densities. We demonstrate how simple neural network (NN)–based local density approximation (LDA) and generalized gradient approximation (GGA), trained on just five atoms and two molecules, provide remarkable improvement in total energies and densities. Particularly, the NN-based GGA functional attains similar accuracy as the higher rung SCAN meta-GGA on various thermochemistry datasets. These results underscore the promise of using the XC potential in modeling XC functionals and can pave the way for systematic learning of increasingly accurate XC functionals.

Science & Technology - Other Topics↗
DOE OSTI2022

Neutron Star Radii, Deformabilities, and Moments of Inertia from Experimental and Ab Initio Theory Constraints of the 208Pb Neutron Skin Thickness

Recent experimental and ab initio theory investigations of the 208Pb neutron skin thickness have the potential to inform the neutron star equation of state. In particular, the strong correlation between the 208Pb neutron skin thickness and the pressure of neutron matter at normal nuclear densities leads to modified predictions for the radii, tidal deformabilities, and moments of inertia of typical 1.4M⊙ neutron stars. In the present work, we study the relative impact of these recent analyses of the 208Pb neutron skin thickness on bulk properties of neutron stars within a Bayesian statistical analysis. Two models for the equation of state prior are employed in order to highlight the role of the highly uncertain high-density equation of state. From our combined Bayesian analysis of nuclear theory, nuclear experiment, and observational constraints on the dense matter equation of state, we find at the 90% credibility level R1.4=12.36−0.73+0.38 km for the radius of a 1.4M⊙ neutron star, R2.0=11.96−0.71+0.94 km for the radius of a 2.0M⊙ neutron star, Λ1.4=440−144+103 for the tidal deformability of a 1.4M⊙ neutron star, and I1.338=1.425−0.146+0.074×1045gcm2 for the moment of inertia of PSR J0737-3039A whose mass is 1.338M⊙.

Lim, Yeunhwan↗
DOE OSTI2023

Double-Atom Catalysts Featuring Inverse Sandwich Structure for CO 2 Reduction Reaction: A Synergetic First-Principles and Machine Learning Investigation

Electrocatalytic CO 2 reduction reactions (CO 2 RR) based on scalable and highly efficient catalysis provide an attractive strategy for reducing CO 2 emissions. Here in this work, we combined first-principles density functional theory (DFT) and machine learning (ML) to comprehensively explore the potential of double-atom catalysts (DACs) featuring an inverse sandwich structure anchored on defective graphene (gra) to catalyze CO 2 RR to generate C 1 products. We started with five homonuclear M 2 ⊥gra (M = Co, Ni, Rh, Ir, and Pt), followed by 127 heteronuclear MM'⊥gra (M = Co, Ni, Rh, Ir, and Pt, M' = Sc–Au). Stable DACs were screened by evaluating their binding energy, formation energy, and dissolution potential of metal atoms, as well as conducting first-principles molecular dynamics simulations with and without solvent water molecules. Based on DFT calculations, Rh 2 ⊥gra DAC was found to outperform the other four homonuclear DACs and the Rh-based single- and double-atom catalysts of noninverse sandwich structures. Out of the 127 heteronuclear DACs, 14 were found to be stable and have good catalytic performance. An ML approach was adopted to correlate key factors with the activity and stability of the DACs, including the sum of radii of metal and ligand atoms (d M–M' , d M–C , and d M'–C ), the sum and difference of electronegativity of two metal atoms (P M + P M' , P M – P M '), the sum and difference of first ionization energy of two metal atoms (I M + I M' , I M – I M '), the sum and difference of electron affinity of two metal atoms (A M + A M' , A M – A M '), and the number of d-electrons of the two metal atoms (Nd). The obtained ML models were further used to predict 154 potential electrocatalysts out of 784 possible DACs featuring the same inverse sandwich configuration. Overall, this work not only identified promising CO 2 RR DACs featuring the reported inverse sandwich structure but also provided insights into key atomic characteristics associated with high CO 2 RR activity.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Effect of applied potential on metal surfaces: Surface energy, Wulff shape and charge distribution

Here we use grand canonical density functional theory to predict the surface energies, Wulff shapes, charge distributions and catalytically active sites of different metal surfaces under electrochemical conditions. We propose a method for computing surface energies from grand canonical density functional theory (GC-DFT) calculations of periodic slab models and use it to compute the surface energies of the facets of Pt, Cu, and Ag crystals to predict their Wulff shapes under electrochemical conditions. GC-DFT predicts that, for the pure metals studied, solvation only slightly affects the Wulff shape while applied potentials considerably affect the surface energies and corresponding Wulff shapes. We used Bader charge analysis of GC-DFT computed electron densities to investigate the effect of applied potential on the distribution of electron density over the atoms of the surfaces of Pt, Cu, Ag, and the 75–25 Ag-Pt and Au-Ni alloys. This analysis shows that, under an applied potential, the electron density is unevenly distributed over the surface atoms and that the charges of atoms more exposed to solvent are more sensitive to bias. Our results show that the most sensitive atom to bias can be used to identify the most favorable adsorption site and thus, the active sites of electrochemical reactions, which is computationally less demanding than calculating the adsorption energies on all possible adsorption sites.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Single-gap isotropic s- wave superconductivity in single crystals AuSn 4

In this article, London, λ L (T), and Campbell, λ C (T), penetration depths were measured in single crystals of a topological superconductor candidate AuSn 4 . At low temperatures, λ L (T) is exponentially attenuated and, if fitted with the power law, λ(T) ~ T n , gives exponents n > 4, indistinguishable from the isotropic single s- wave gap Bardeen-Cooper-Schrieffer (BCS) asymptotic. The superfluid density fits perfectly in the entire temperature range to the BCS theory. The superconducting transition temperature, T c = 2.40 ± 0.05 K, does not change after 2.5 MeV electron irradiation, indicating the validity of the Anderson theorem for isotropic s- wave superconductors. Campbell penetration depth before and after electron irradiation shows no hysteresis between the zero-field cooling (ZFC) and field cooling (FC) protocols, consistent with the parabolic pinning potential. Interestingly, the critical current density estimated from the original Campbell theory decreases after irradiation, implying that a more sophisticated theory involving collective effects is needed to describe vortex pinning in this system. In general, our thermodynamic measurements strongly suggest that the bulk response of the AuSn 4 crystals is fully consistent with the isotropic s- wave weak-coupling BCS superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Concordant Mode Approach (CMA): Vibrational Analysis of New and Upgraded Intermolecular Benchmarks for Noncovalent Bonding

The Concordant Mode Approach (CMA) is a novel method that offers tremendous potential for increasing the system size and the level of theory attainable in quantum chemical computations of molecular vibrational frequencies. To investigate the extension of CMA to intermolecular vibrations, computations with coupled cluster singles and doubles with perturbative triples theory [CCSD(T)] using two augmented correlation-consistent polarized-valence triple-ζ basis sets (aug-cc-pVTZ or h-aug-ccpVTZ) were performed on 17 prototypical loosely bound complexes of hydrogen-bonded, dispersion, and mixed character. These Level A results provide new and upgraded benchmarks for noncovalent bonding and a severe test for CMA vibrational analyses. The Level A target frequencies were recovered remarkably well using second-order Møller−Plesset perturbation theory (MP2) with h-aug-cc-pVTZ for generating the underlying (Level B) normal modes of the CMA scheme. Employing this Level B within the lowest-rung CMA-0A method reproduces the 435 benchmark frequencies with a mean absolute error (MAE) of 0.23 cm −1 and a corresponding standard deviation (σ) of 0.84 cm −1 ; strikingly, the corresponding subset of 106 interfragment frequencies exhibits MAE = 0.34 cm −1 and σ = 0.90 cm −1 . Subsequent application of the higher-rung CMA-2A scheme eliminates all outliers and reduces the overall MAE to a minuscule 0.08 cm −1 with the inclusion of only 3.0% of the off-diagonal couplings not accounted for by CMA-0A. Accordingly, the highly efficient CMA methodology proves to be robust even for vibrations on flat potential energy surfaces.

Aromatic compounds↗
DOE OSTI2024

Transferable Water Potentials Using Equivariant Neural Networks

Machine learning interatomic potentials (MLIPs) have emerged as a technique that promises quantum theory accuracy for reduced cost. It has been proposed [J. Chem. Phys. 2023, 158, 084111] that MLIPs trained on solely liquid water data cannot accurately transfer to the vapor–liquid equilibrium while recovering the many-body decomposition (MBD) analysis of gas-phase water clusters. This suggests that MLIPs do not directly learn the physically correct interactions of water molecules, limiting transferability. In this work, we show that MLIPs using equivariant architecture and trained on 3200 liquid water structures reproduces liquid-phase water properties (e.g., density within 0.003 g/cm 3 between 230 and 365 K), vapor–liquid equilibrium properties up to 550 K, the MBD analysis of gas-phase water cluster up to six-body interactions, and the relative energy and the vibrational density of states of ice phases. We show that potentials developed using equivariant MLIPs allow transferability for arbitrary phases of water that remain stable in nanosecond long simulations.

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

Isospin composition of fission barriers

We employ a microscopic method to study how isospin affects the fission potential of 240 Pu. Our approach uses constrained Hartree-Fock theory (CHF) which allows us to separately investigate the isoscalar and isovector properties of the nuclear energy density functional (EDF). By analyzing the isoscalar and isovector components of the EDF along the fission path we can assess the isovector contribution to fission barriers. Here, we study this effect for the fully adiabatic path to scission. The isovector component of the fission potential is found to increase in magnitude as the nucleus evolves towards scission, exemplifying the importance of stringent constraints on the isovector sector of the nuclear EDF for reliable predictions of fission properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗