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

Microstate counting via Bethe Ansätze in the 4d $ \mathcal{N} $ = 1 superconformal index

We study the superconfomal index of four-dimensional toric quiver gauge theories using a Bethe Ansatz approach recently applied by Benini and Milan. Relying on a particular set of solutions to the corresponding Bethe Ansatz equations we evaluate the superconformal index in the large N limit, thus avoiding to take any Cardy-like limit. We present explicit results for theories arising as a stack of N D 3 branes at the tip of toric Calabi-Yau cones: the conifold theory, the suspended pinch point gauge theory, the first del Pezzo theory and Y p,q quiver gauge theories. For a suitable choice of the chemical potentials of the theory we find agreement with predictions made for the same theories in the Cardy-like limit. However, for other regions of the domain of chemical potentials the superconformal index is modified and consequently the associated black hole entropy receives corrections. We work out explicitly the simple case of the conifold theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Molecular interactions from the density functional theory for chemical reactivity: Interaction chemical potential, hardness, and reactivity principles

In the first paper of this series, the authors derived an expression for the interaction energy between two reagents in terms of the chemical reactivity indicators that can be derived from density functional perturbation theory. While negative interaction energies can explain reactivity, reactivity is often more simply explained using the “|dμ| big is good” rule or the maximum hardness principle. Expressions for the change in chemical potential (μ) and hardness when two reagents interact are derived. A partial justification for the maximum hardness principle is that the terms that appear in the interaction energy expression often reappear in the expression for the interaction hardness, but with opposite sign.

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

Collision integrals within the Chapman–Enskog theory for a generalized Lennard-Jones potential

Here we report the values of the collision integrals, needed for the calculation of the macroscopic transport properties such as viscosity (η) and diffusion coefficient (D) of gases within the Chapman–Enskog kinetic gas theory, for a generalized Lennard-Jones potential (gLJ), a more general potential with an adjustable long range 1/r dependence that can describe a wide range of intermolecular interactions.

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

Deep learning-based quasi-continuum theory for structure of confined fluids

Predicting the structural properties of water and simple fluids confined in nanometer scale pores and channels is essential in, for example, energy storage and biomolecular systems. Classical continuum theories fail to accurately capture the interfacial structure of fluids. In this work, we develop a deep learning-based quasi-continuum theory (DL-QT) to predict the concentration and potential profiles of a Lennard-Jones (LJ) fluid and water confined in a nanochannel. The deep learning model is built based on a convolutional encoder–decoder network (CED) and is applied for high-dimensional surrogate modeling to relate the fluid properties to the fluid–fluid potential. The CED model is then combined with the interatomic potential-based continuum theory to determine the concentration profiles of a confined LJ fluid and confined water. Further, we show that the DL-QT model exhibits robust predictive performance for a confined LJ fluid under various thermodynamic states and for water confined in a nanochannel of different widths. The DL-QT model seamlessly connects molecular physics at the nanoscale with continuum theory by using a deep learning model.

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

Low-temperature breakdown of many-body perturbation theory for thermodynamics

It is shown analytically and numerically that the finite-temperature many-body perturbation theory in the grand canonical ensemble has zero radius of convergence at zero temperature when the energy ordering or degree of degeneracy for the ground state changes with the perturbation strength. When the degeneracy of the reference state is partially or fully lifted at the first-order Hirschfelder-Certain degenerate perturbation theory, the grand potential and internal energy diverge as T → 0. Contrary to earlier suggestions of renormalizability by the chemical potential μ, this nonconvergence, first suspected by W. Kohn and J. M. Luttinger, is caused by the nonanalytic nature of the Boltzmann factor e –E/k e T at T = 0, also plaguing the canonical ensemble, which does not involve chemical potential. The finding reveals a fundamental flaw in perturbation theory, which is deeply rooted in the mathematical limitation of power-series expansions and is unlikely to be removed within its framework.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Ca Cobaltites as Potential Cathode Materials for Rechargeable Ca-Ion Batteries: Theory and Experiment

Rechargeable Ca-based batteries can potentially achieve higher energy capacity than Li-ion batteries. The development of Ca-ion batteries, however, remains in its infancy, especially due to the challenge of finding cathode materials with high reversible capacity. In this work, we investigate properties of Ca cobaltites as Ca-ion intercalation cathodes, by means of density functional theory calculations validated by synthesis and electrochemical measurements. Computationally, we accessed thermodynamic, diffusion properties, energy capacity, and the voltage profiles for four Ca cobaltite compounds of different stoichiometry and Co oxidation states: Ca 3 Co 2 O 6 , CaCo 2 O 4 , Ca 2 Co 2 O 5 , and [Ca 2 CoO 3 ][CoO 2 ] 1.62 . We found good stability and relatively low migration barriers of some Ca cobaltites during cycling with the layered CaCo 2 O 4 having the lowest Ca migration barrier of 0.7 eV and the highest theoretical capacity. To validate our calculations, we synthesized Ca 3 Co 2 O 6 , CaCo 2 O 4 , and [Ca 2 CoO 3 ][CoO 2 ] 1.62 with the Pechini method and, subsequently, to test the electrochemical extraction of Ca.

25 ENERGY STORAGE↗
DOE OSTI2021

Renormalization group improvement of the effective potential: an EFT approach

We apply effective field theory (EFT) methods to compute the renormalization group improved effective potential for theories with a large mass hierarchy. Our method allows one to compute the effective potential in a systematic expansion in powers of the mass ratio, as well as to sum large logarithms of mass ratios using renormalization group evolution. The effective potential is the sum of one-particle irreducible diagrams (1PI) but information about which diagrams are 1PI is lost after matching to the EFT, since heavy lines get shrunk to a point. We therefore introduce a tadpole condition in place of the 1PI condition, and use the renormalization group improved value of the tadpole in computing the effective potential. We explain why the effective potential computed using an EFT is not the same as the effective potential of the EFT. We illustrate our method using the O(N) model, a theory of two scalars in the unbroken and broken phases, and the Higgs-Yukawa model. Our leading-log result, obtained by integrating the one-loop β-functions, correctly reproduces the log-squared term in explicit two-loop calculations. Our method does not have a Goldstone boson infrared divergence problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Shadow Molecular Dynamics and Atomic Cluster Expansions for Flexible Charge Models

Here, a shadow molecular dynamics scheme for flexible charge models is presented where the shadow Born–Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, is modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born–Oppenheimer molecular dynamics (BOMD). XL-BOMD provides stable dynamics while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO 2 ) and a molecular system of liquid water. The combined ACE+XL-QEq molecular dynamics simulations are stable over a wide range of temperatures both for the oxide and for the molecular systems and provide a precise sampling of the Born–Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

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

Shadow energy functionals and potentials in Born–Oppenheimer molecular dynamics

In Born–Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born–Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born–Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born–Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born–Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original “zeroth-level” shadow energy functionals and potentials used in extended Lagrangian BOMD. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree–Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

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

Shadow molecular dynamics and atomic cluster expansions for flexible charge models

A shadow molecular dynamics scheme for flexible charge models is presented, where the shadow Born-Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, are modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born-Oppenheimer molecular dynamics (BOMD) [Eur. Phys. J. B 94, 164 (2021)]. XL-BOMD provides a stable dynamics, while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using the atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO2) and a molecular system of liquid water. The combined ACE + XL-QEq dynamics are stable over a wide range of temperatures both for the oxide and the molecular systems, and provide a precise sampling of the Born-Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Medium-range order and compositional correlation in metallic glasses

The compositional atomic ordering in metallic glasses was studied by simulation focusing on the medium-range order (MRO). Many metallic alloy liquids and glasses show MRO characterized by the oscillations in the atomic pair-distribution function (PDF) beyond the first peak, which decay exponentially with distance. To study the effects of the local chemical order on MRO, we examine the compositionally resolved PDF and its MRO for models of various binary metallic alloy glasses. We show that compositional ordering is limited mostly to the nearest-neighbor atoms and the MRO is largely independent of the compositional order. For some elements that strongly repel each other in the alloy, a second MRO periodicity is observed owing to the distinct correlations among them. These results are discussed in light of the idea that the MRO oscillations in the PDF describe the correlations in the atomic density fluctuations, rather than the detailed local atomic structure.

Atomic structure↗
DOE OSTI2022

Simulations of x-ray absorption spectra for CO desorbing from Ru(0001) with transition-potential and time-dependent density functional theory approaches

The desorption of a carbon monoxide molecule from a Ru(0001) surface was studied by means of X-ray Absorption Spectra (XAS) computed with Transition Potential (TP-DFT) and Time Dependent (TD-DFT) DFT methods. By unraveling the evolution of the CO electronic structure upon desorption, we observed that at 2.3 Å from the surface, the CO molecule has already predominantly gas-phase character. While C 1s XAS is quite insensitive to changes in the C–O bond length, the O 1s excitation is very sensitive with the π* coming down in energy upon CO bond stretching, which competes with the increase in orbital energy due to the repulsive interaction with the metallic surface. We show in a systematic way that the TP-DFT method can describe the XAS rather well at the endpoints (chemisorbed and gas phase) but is affected by artificial charge transfer and/or incorrect spin treatment in the transition region in cases like CO, where there are low-lying π* orbitals and large exchange interactions between the core 1s and valence-acceptor π* orbitals. As an alternative, we demonstrate by comparing with experimental data that a linear response approach using TD-DFT employing common exchange-correlation functionals and finite-size clusters can yield a good description of the spectral evolution of the 1s → π* transition with correct spin and gas-to-chemisorbed chemical shifts in good agreement with experiment.

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

Ab-initio simulation of spin-vibronic spectra of methoxy radical

Despite the fact that experimental and theoretical work on the spectrum of methoxy has stretched from the microwave to the ultraviolet and proceeded for nearly 50 years, parts of the spectrum have remained a challenge to simulate theoretically and make reliable line-by-line assignments. The spectral complexity arises because the radical has a non-zero electron spin and significant vibronic coupling between the two elec- tronic components of the ground state due to the presence of a conical intersection. This work describes a completely ab initio effort to understand and assign the spin- vibronic levels of the X 2E state from 0 to above 3000 cm−1, a region that includes the fundamental transitions of the C-H symmetric and asymmetric stretches that have not previously been identified uniquely. A potential energy surface for methoxy was calculated at the EOM-CCSDT/ANO1 level of theory. Subsequently this potential energy surface was fit to a quartic power series expansion of all nine vibrational nor- mal coordinates (as determined at the minimum of the conical intersection) by the use of a machine-learning-based algorithm. After the addition of spin-orbit coupling, the spin-vibronic problem was solved using both the Krylov-Schur and Lanczos algorithms with the SOCJT3 software to converge eigenvalues up to 3500 cm−1 and their eigen- vectors. The latter were used, in conjunction with the calculated dipole moment and its derivatives (calculated using finite differences at EOM-CCSDT/ANO1 level), to determine spectral intensities for the spin-vibronic spectra. The calculated transition frequencies and intensities were used to simulate and assign the observed transitions of the spin-vibronic spectra of the radical. The credibility of the assignments and their significance is discussed in detail.

Sharma, Ketan [University of Florida, Gainesville,↗
DOE OSTI2022

Multilevel Computational Studies Reveal the Importance of Axial Ligand for Oxygen Reduction Reaction on Fe–N–C Materials

The systematic improvement of Fe–N–C materials for fuel cell applications has proven challenging, due in part to an incomplete atomistic understanding of the oxygen reduction reaction (ORR) under electrochemical conditions. Herein, a multilevel computational approach, which combines ab initio molecular dynamics simulations and constant potential density functional theory calculations, is used to assess proton-coupled electron transfer (PCET) processes and adsorption thermodynamics of key ORR intermediates. These calculations indicate that the potential-limiting step for ORR on Fe–N–C materials is the formation of the Fe III –OOH intermediate. They also show that an active site model with a water molecule axially ligated to the iron center throughout the catalytic cycle produces results that are consistent with the experimental measurements. In particular, reliable prediction of the ORR onset potential and the Fe(III/II) redox potential associated with the conversion of Fe III –OH to Fe II and desorbed H 2 O requires an axial H 2 O co-adsorbed to the iron center. The observation of a five-coordinate rather than four-coordinate active site has significant implications for the thermodynamics and mechanism of ORR. Importantly, these findings highlight the importance of solvent–substrate interactions and surface charge effects for understanding the PCET reaction mechanisms and transition-metal redox couples under realistic electrochemical conditions.

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

Self-interaction-corrected Kohn–Sham effective potentials using the density-consistent effective potential method

Density functional theory (DFT) and beyond-DFT methods are often used in combination with photoelectron spectroscopy to obtain physical insights into the electronic structure of molecules and solids. The Kohn–Sham eigenvalues are not electron removal energies except for the highest occupied orbital. The eigenvalues of the highest occupied molecular orbitals often underestimate the electron removal or ionization energies due to the self-interaction (SI) errors in approximate density functionals. Here, we adapt and implement the density-consistent effective potential method of Kohut, Ryabinkin, and Staroverov (2014) to obtain SI-corrected local effective potentials from the SI-corrected Fermi–Löwdin orbitals and density in the Fermi–Löwdin orbital self-interaction correction scheme. The implementation is used to obtain the density of states (photoelectron spectra) and HOMO–LUMO gaps for a set of molecules and polyacenes. Good agreement with experimental values is obtained compared to a range of SI uncorrected density functional approximations.

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

Developing Machine Learning Interatomic Potential for Fe-Cr-Ni Alloys

Accurate prediction of creep and fatigue behavior of stainless steel at elevated temperatures in hydrogen environment requires fundamental understanding of alloy-hydrogen interaction at cross-scale including bulk lattice and key defects such as vacancies, grain boundaries, surfaces, stacking faults, dislocations, and precipitates. This project aims to predict creep behavior of 347H stainless steel with H using machine learning interatomic potentials based on first-principles density functional theory simulations. The Moment Tensor Potentials platform is adopted for this work since it demonstrates a fine balance between model accuracy and computational efficiency. The potential is well trained based on large amount of high-fidelity density functional theory calculations. The validation is carried out by comparing various important properties including short range order, coefficient of thermal expansion, elastic properties, stacking fault energy, grain boundary energy, and surface energy. This work lays the foundation for reliable atomistic simulation of high temperature hydrogen attack of stainless steel.

density functional theory (DFT)↗
DOE OSTI2020

Charge radii of exotic neon and magnesium isotopes

We compute the charge radii and ground-state energies of even-mass neon and magnesium isotopes from neutron number N = 8 to the dripline. Our calculations are based on nucleon-nucleon and three-nucleon potentials from chiral effective field theory that include Δ isobars. These potentials yield an accurate saturation point and symmetry energy of nuclear matter. We use the coupled-cluster method and start from an axially symmetric reference state. Binding energies and two-neutron separation energies largely agree with data, and the dripline in neon is accurate. The computed charge radii are accurate for many isotopes where data exist. Finer details, such as isotope shifts, however, are not accurately reproduced. These chiral potentials indicate a subshell closure at N = 14 for the radii (but not for two-neutron separation energies) and a decrease in charge radii at N = 8 (observed in neon and predicted for magnesium). Furthermore, they yield a continued increase of charge radii as neutrons are added beyond N = 14 yet underestimate the large increase at N = 20 in magnesium.

20 ≤ A ≤ 38↗
DOE OSTI2022

Precise $\textit{Q}$-value measurements of 112,113 Ag and 115 Cd with the Canadian Penning trap for evaluation of potential ultralow $\textit{Q}$-value $\textit{β}$ decays

Background: An ultralow $\textit{Q}$-value $\textit{β}$ decay can occur from a parent nuclide to an excited nuclear state in the daughter such that $Q_{\text{UL}}$ 1 keV. These decay processes are of interest for nuclear $\textit{β}$-decay theory and as potential candidates in neutrino mass determination experiments. To date, only one ultralow $\textit{Q}$-value $\textit{β}$ decay has been observed—that of 115 In with $Q_β$ = 147(10) eV. A number of other potential candidates exist, but improved mass measurements are necessary to determine if these decay channels are energetically allowed and, in fact, ultralow. Purpose: To perform precise $\textit{β}$-decay $\textit{Q}$-value measurements of 112,113 Ag and 115 Cd and to use them in combination with nuclear energy level data for the daughter isotopes 112,113 Cd and 115 In to determine if the potential ultralow $\textit{Q}$-value $\textit{β}$-decay branches of 112,113 Ag and 115 Cd are energetically allowed and 1 keV. Method: The Canadian Penning Trap at Argonne National Laboratory was used to measure the cyclotron frequency ratios of singly charged 112,113 Ag and 115 Cd ions with respect to their daughters 112,113 Cd and 115 In. From these measurements, the ground-state to ground-state $\textit{β}$-decay $\textit{Q}$ values were obtained. Results: The 112 Ag → 112 Cd, 113 Ag → 113 Cd, and 115 Cd → 115 In $\textit{β}$-decay $\textit{Q}$ values were measured to be $Q_β$( 112 Ag) = 3990.16(22) keV, $Q_β$( 113 Ag) = 2085.7(4.6) keV, and $Q_β$( 115 Cd) = 1451.36(34) keV. These results were compared to energies of excited states in 112 Cd at 3997.75(14) keV, 113 Cd at 2015.6(2.5) and 2080(10) keV, and 115 In at 1448.787(9) keV, resulting in precise $Q_{\text{UL}}$ values for the potential decay channels of –7.59(26) keV, 6(11) keV, and 2.57(34) keV, respectively. Conclusion: The potential ultralow $\textit{Q}$-value decays of 112 Ag and 115 Cd have been ruled out. 113 Ag is still a possible candidate until a more precise measurement of the 2080(10) keV, 1/2+ state of 113 Cd is available. In the course of this work we have found the ground state mass of 113 Ag reported in the 2020 Atomic Mass Evaluation [Wang et al., Chin. Phys. C 45, 030003 (2021)] to be lower than our measurement by 69(17) keV (a 4σ discrepancy).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗