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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1987

4U 1820-30 as a potential test of the nonsymmetric gravitational theory of Moffat

Recent observations of the X-ray burst source 4U 1820-30 have revealed a 685 s modulation of the luminosity. How this system could provide a stringent test of the nonsymmetric gravitational theory (NGT) of Moffat (1979), provided the observed periodicity is due to orbital motion of a binary system, is discussed. The possible orbital period change predicted by general relativity may be detectable in this system.

Krisher, Timothy P.↗
NASA NTRS1979

General Theory of Aerodynamic Instability and the Mechanism of Flutter

The aerodynamic forces on an oscillating airfoil or airfoil-aileron combination of three independent degrees of freedom were determined. The problem resolves itself into the solution of certain definite integrals, which were identified as Bessel functions of the first and second kind, and of zero and first order. The theory, based on potential flow and the Kutta condition, is fundamentally equivalent to the conventional wing section theory relating to the steady case. The air forces being known, the mechanism of aerodynamic instability was analyzed. An exact solution, involving potential flow and the adoption of the Kutta condition, was derived. The solution is of a simple form and is expressed by means of an auxiliary parameter k. The flutter velocity, treated as the unknown quantity, was determined as a function of a certain ratio of the frequencies in the separate degrees of freedom for any magnitudes and combinations of the airfoil-aileron parameters.

Theodorsen, Theodore↗
NASA NTRS2018

Clockwork Axions in Cosmology. Is Chromonatural Inflation Chrononatural?

Many cosmological models rely on large couplings of axions to gauge fields. Examples include theories of magnetogenesis, inflation on a steep potential, chiral gravitational waves, and chromonatural inflation. Such theories require a mismatch between the axion field range and the mass scale appearing in the aFF~ coupling. This mismatch suggests an underlying monodromy, with the axion winding around its fundamental period a large number of times. We investigate the extent to which this integer can be explained as a product of smaller integers in a UV completion: in the parlance of our times, can the theory be “clockworked”? We argue that a clockwork construction producing a potential μ4cos(ajFa) for an axion of fundamental period Fa will obey the constraint μ < Fa. For a some applications, including chromonatural inflation with sub-Planckian field range, this constraint obstructs a clockwork UV completion. Alternative routes to a large coupling include fields of large charge (an approach limited by strong coupling) or kinetic mixing (requiring a lighter axion). Our results suggest that completions of axion cosmologies that explain the large parameter in the theory potentially alter the phenomenological predictions of the model.

Cosmology of Theories beyond the SM↗
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↗
NASA NTRS1977

Unsteady subsonic and supersonic potential aerodynamics for complex configurations

A recently developed general theory for unsteady compressible potential fluid dynamics for complex-configuration aircraft is reviewed. The method is based on a combination of the following techniques: Green's function method (to transform the differential equation into an integral differential-delay equation), finite element method (to transform the equation into a set of differential-delay equations in time), and the Laplace transform method (to transform the differential-delay equations into algebraic equations).

Morino, L.↗
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,↗
NASA NTRS1986

The theory of Langmuir probes in strong electrostatic potential structures

The operation of collecting and emitting Langmuir probes and double probes within time-stationary strong electrostatic potential structures is analyzed. The cross sections of spherical and cylindrical probes to charged particles within the structures are presented and used to obtain the current-voltage characteristics of idealized probes. The acquisition of plasma parameters from these characteristics is outlined, and the operation of idealized floating double-probe systems is analyzed. Probe surface effects are added to the idealized theory, and some surface effects pertinent to spacecraft probes are quantified. Magnetic field effects on idealized probes are examined, and the time required for floating probes to change their potentials by collecting charge and by emitting photoelectrons is discussed. Calculations on the space-charge effects of probe-perturbed beams and on the space-charge limiting of electron emission are given in an appendix.

Borovsky, J. E.↗
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↗
NASA NTRS1996

Polymer-Induced Depletion Interaction and Its Effect on Colloidal Sedimentation in Colloid-Polymer Mixtures

In this paper we focus on the polymer-induced depletion attraction and its effect on colloidal sedimentation in colloid-polymer mixtures. We first report a small angle neutron scattering (SANS) study of the depletion effect in a mixture of hard-sphere-like colloid and non-adsorbing polymer. Then we present results of our recent sedimentation measurements in the same colloid-polymer mixture. A key parameter in controlling the sedimentation of heavy colloidal particles is the interparticle potential U(tau), which is the work required to bring two colloidal particles from infinity to a distance tau under a give solvent condition. This potential is known to affect the average settling velocity of the particles and experimentally one needs to have a way to continuously vary U(tau) in order to test the theory. The interaction potential U(tau) can be altered by adding polymer molecules into the colloidal suspension. In a mixture of colloid and non-adsorbing polymer, the potential U(tau) can develop an attractive well because of the depletion effect, in that the polymer chains are expelled from the region between two colloidal particles when their surface separation becomes smaller than the size of the polymer chains. The exclusion of polymer molecules from the space between the colloidal particles leads to an unbalanced osmotic pressure difference pushing the colloidal particles together, which results in an effective attraction between the two colloidal particles. The polymer-induced depletion attraction controls the phase stability of many colloid-polymer mixtures, which are directly of interest to industry.

Tong, Penger↗
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↗