Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dipolar molecule emulator of lattice gauge theories (Final Report)

The quantum many-body problem is a great unsolved, cross-cutting challenge in physics that is of fundamental importance. Our understanding of phenomena related to dense quark matter, in particular, is challenged by this practical intractability of classical simulations. Because of the sheer cost and challenge of performing experiments that probe the length and energy scales relevant to such physics, there are practical motivations for finding theoretical methods to address the many-body problem. One promising approach is based on the use of a programmable and controllable analog quantum systems to emulate the physics of many-body problems of interest. While there are currently broad efforts to develop mid- to large-scale quantum computers, we are still likely many years away from such devices outperforming classical supercomputers for useful calculations. Even though digital quantum computers are still at a stage too premature for such tasks, approaches based on analog quantum simulators have advanced rapidly over the past two decades and can now treat many-body problems of interest on small- to medium-scale systems. Tackling the challenging many-body problems of relevance to high energy physics represents one of the great new opportunities and challenges for the field of analog quantum simulation. To address this challenge, this award brought together an interdisciplinary team of physicists with expertise in high energy theory (El-Khadra and Draper), computational and condensed matter physics (Clark), and experimental quantum and atomic, molecular, optical (AMO) science (DeMarco and Gadway), with the goal of developing realistic and optimized strategies for the analog quantum simulation of lattice gauge theories using cold atoms and molecules. This track 1 project included both theoretical and experimental goals, stated as follows. The theory goals were to "develop and verify a novel approach to quantum emulation of (1+1)d quantum link models, based on arrays of trapped quantum spins with long-range interactions." This final project report details the progress made under this award, as well as new research directions developed under this award.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Size-dependent attraction of Cu solutes to clusters formed at Ag grain boundaries

We report a size-dependent solute clustering mechanism at grain boundaries in a sputtered ultrafine-grained Ag-Cu alloy, where large Cu clusters form despite weak individual solute-solute interactions. X-ray diffraction confirms limited Cu solubility in the Ag matrix, while scanning transmission electron microscopy reveals Cu clustering at both ordinary GBs and GB junctions. Density functional theory calculations show that 12-atom Cu clusters are energetically preferred, while smaller three-atom clusters are significantly less stable. Additional calculations demonstrate a marked increase in solute-cluster attraction energy with cluster size. As a result, these findings point to a previously unrecognized pathway for grain-boundary solute clustering in immiscible systems, driven by collective solute-cluster interactions, with implications for segregation behavior and stability in nanocrystalline and ultrafine-grained alloys.

Alloys↗
NASA NTRS2019

AladynPi – Adaptive Neural Network Molecular Dynamics Simulation Code with Physically Informed Potential: Computational Materials Mini-Application

This report provides an overview and description of commands used in the Computational Materials mini-application, AladynPi. AladynPi is an extension of a previously released mini-application, Aladyn (https://github.com/nasa/aladyn; Yamakov, V.I., and Glaessgen, E.H., NASA/TM-2018-220104). Aladyn and AladynPi are basic molecular dynamics codes written in FORTRAN 2003, which are designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory method. An input for the ANN is a set of structure coefficients, characterizing the local atomic environment of each atom, for which the atomic energy is obtained in the ANN inference process. In Aladyn, the ANN gives directly the energy of interatomic interactions. In AladynPi, the ANN gives optimized parameters for a predefined empirical function, known as bond-order-potential (BOP). The parameterized BOP function is then used to calculate the energy. AladynPi code is being released to serve as a training testbed for students and professors in academia to explore possible optimization algorithms for parallel computing on multicore central processing unit (CPU) computers or computers utilizing manycore architectures based on graphic processing units (GPUs). The effort is supported by the High Performance Computing incubator (HPCi) project at NASA Langley Research Center.

Yamakov, Vesselin I.↗
DOE OSTI2022

Improving Predictive Capability in REHEDS Simulations with Fast, Accurate, and Consistent Non-Equilibrium Material Properties

Predictive design of REHEDS experiments with radiation-hydrodynamic simulations requires knowledge of material properties (e.g. equations of state (EOS), transport coefficients, and radiation physics). Interpreting experimental results requires accurate models of diagnostic observables (e.g. detailed emission, absorption, and scattering spectra). In conditions of Local Thermodynamic Equilibrium (LTE), these material properties and observables can be pre-computed with relatively high accuracy and subsequently tabulated on simple temperature-density grids for fast look-up by simulations. When radiation and electron temperatures fall out of equilibrium, however, non-LTE effects can profoundly change material properties and diagnostic signatures. Accurately and efficiently incorporating these non-LTE effects has been a longstanding challenge for simulations. At present, most simulations include non-LTE effects by invoking highly simplified inline models. These inline non-LTE models are both much slower than table look-up and significantly less accurate than the detailed models used to populate LTE tables and diagnose experimental data through post-processing or inversion. Because inline non-LTE models are slow, designers avoid them whenever possible, which leads to known inaccuracies from using tabular LTE. Because inline models are simple, they are inconsistent with tabular data from detailed models, leading to ill-known inaccuracies, and they cannot generate detailed synthetic diagnostics suitable for direct comparisons with experimental data. This project addresses the challenge of generating and utilizing efficient, accurate, and consistent non-equilibrium material data along three complementary but relatively independent research lines. First, we have developed a relatively fast and accurate non-LTE average-atom model based on density functional theory (DFT) that provides a complete set of EOS, transport, and radiative data, and have rigorously tested it against more sophisticated first-principles multi-atom DFT models, including time-dependent DFT. Next, we have developed a tabular scheme and interpolation methods that compactly capture non-LTE effects for use in simulations and have implemented these tables in the GORGON magneto-hydrodynamic (MHD) code. Finally, we have developed post-processing tools that use detailed tabulated non-LTE data to directly predict experimental observables from simulation output.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Univariate Prediction of Hammett Parameters and Select Relative Reaction Rates Using Loewdin Atomic Charges

Loewdin charges from density functional theory calculations were used here to obtain general, univariate linear correlations for the prediction of experimental Hammett parameters and relative reaction rates. While previous studies have established that Hirshfeld and CM5 charges perform strongly as univariate predictors, the near-ubiquitous Loewdin charges have not yet been evaluated. To this end, we assess the predictive capability of Loewdin charges for three chemical systems. First, we show that Loewdin charges outperform Hirshfeld and CM5 charges for Hammett parameter prediction. Second, we see that Loewdin charges generally perform comparably to Hirshfeld charges for predicting the relative rates of olefin cleavage by photoexcited nitroarenes. The single case of poor correlation, between relative rates and the Loewdin charges on nitrogen sites, is ameliorated when considering the net charge on the NO 2 group. Third, we show that Loewdin, Hirshfeld, and CM5 charges all perform very well for generating correlations for relative reaction rates for C–H activation of 9-(4-X-phenyl)-9H-fluorene substrates by a transition metal catalyst. The equations generated throughout the study enable the prediction of Hammett parameters and relative reaction rates. Finally, these tools can accelerate synthetic and experimental studies by enabling the in silico prediction of uncharacterized chemical properties.

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

Optimization of Ag Electrocatalyst Performance for CO2 to CO Conversion: Pairing Atomic Simulations with Experiments

Density-functional theory- based calculations that complement our series of experimental efforts on optimizing Ag electrocatalyst performance for CO2 to CO conversion were presented. Three key findings are drawn out from the combined UHV surface science, STM and electrochemical measurements: (1) Using a series of Ag nanoparticles with 2-6 nm average diameters, CO2 reduction reaction (CO2RR) activity increases, with particle between 2 nm and ∼4 nm demonstrating the highest combination of activity and selectivity; (2) Electronic metal−support interactions (EMSIs) between Ag and C dramatically improve CO2RR performance as evidenced by a scaling relationship between particle size and the relative Ag−C EMSIs strength, which improves the CO2-to-CO Faradaic efficiency of sub-2 nm Ag particles from 2 to ∼100% and increases the CO turnover frequency ∼15-fold compared with similarly sized bare Ag particles; (3) The performance Ag electrocatalysts is improved when supported on S-doped C materials. Computational modeling of 1−10 nm Ag particles predicts a nearly identical size-dependent trend with maximum CO2RR activity predicted for 3.7nm particles. The calculations support the promotional effect of C materials, showing a large charge transfer of 1.02 e from Ag clusters to defective C and a 0.41 eV less endergonic step of forming COOH intermediate on Ag/defective-C compared to the Ag/C system. Calculations indicate a more favorable energetic pathway of CO2-to-CO at the C-S-Ag interface, consistent with experiments.

CO2 utilization↗
NASA NTRS1969

Radiative effects in semiclassical theory.

Unquantized field calculations extended to include atomic field effect on atom, predicting spontaneous decay rate from excited state and light frequency time dependence

Crisp, M. D.↗
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 OSTI2020

Catalytic consequences of hydrogen addition events and solvent-adsorbate interactions during guaiacol-H 2 reactions at the H 2 O-Ru(0 0 0 1) interface

Catalytic reactions of biomass-derived phenolics and H 2 occur on transition metal surfaces via competitive C–O cleavage and ring saturation pathways, with both requiring multiple hydrogen addition events before forming their respective rate limiting transition states. These events are markedly affected by solvent chemical identity, with polar protic solvents ionizing hydrogen adatoms (H*) to interfacial protons (H + ) and opening up new catalytic routes. Here, we establish the reaction coordinate space for guaiacol-H 2 reactions on Ru(0 0 0 1) using density functional theory and describe the atomic-scale effect of a polar protic solvent, H 2 O. Coupled H + and H* attack leads to quasi-equilibrated enol and keto intermediates as the precursors for C–O cleavage and ring saturation, respectively. For C–O cleavage, H 2 O solvent enables a lower energy pathway via concomitant transfer of the hydroxyl H + to the methoxy oxygen during C–OCH 3 cleavage, forming a charge separated [Ru(s)–(C 6 H 5 O – )…(H + )…OCH 3 ] transition state and reducing the barrier by up to 0.8 eV as compared to unassisted C–OCH 3 cleavage. For ring saturation, H* attack onto an unsaturated meta carbon is rate limiting with no direct solvent participation, suggesting that protic polar solvents selectively promote the C–O cleavage pathway. Taken together, we show that activating guaiacol for either C–O bond cleavage or ring saturation product formation depends on the reactive hydrogen identity (H* or H + ), enol/keto isomerization equilibrium, and accessibility of the proton assisted Car–OCH 3 cleavage transition state. Here, all such factors are tunable via changes to the solvent or metal identity.

09 BIOMASS FUELS↗
DOE OSTI2025

Effect of dilute Rh on oxygen dissociation, spillover, and the oxidation of Cu across many orders of magnitude pressure

Knowledge of how trace amounts of more reactive metals influence the oxidation rate and mechanism of Cu surfaces is essential for developing strategies to optimize the performance of Cu-based catalysts. We find that the addition of 1% Rh to Cu(111) increases the initial O 2 dissociation rate by approximately 9-fold. CO poisoning experiments reveal that single Rh atoms activate O 2 and facilitate the spillover of atomic oxygen to Cu sites. Scanning tunneling microscopy (STM) and in situ X-ray photoelectron spectroscopy (XPS) support this mechanism, showing enhanced surface oxygen near Rh atoms. Here, a density functional theory (DFT)-based model demonstrates that Rh binds the O 2 precursor 0.15 eV more strongly than Cu(111) and lowers the O 2 dissociation barrier by 0.02 eV. Both single-crystal and nanoparticle experiments show that at low oxygen pressures, Rh enhances Cu oxidation, whereas at higher pressures, it inhibits deeper oxidation, as evidenced by in situ ultraviolet-visible (UV-vis) spectra.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The role of chalcogen vacancies for atomic defect emission in MoS 2

For two-dimensional (2D) layered semiconductors, control over atomic defects and understanding of their electronic and optical functionality represent major challenges towards developing a mature semiconductor technology using such materials. Here, we correlate generation, optical spectroscopy, atomic resolution imaging, and ab initio theory of chalcogen vacancies in monolayer MoS 2 . Chalcogen vacancies are selectively generated by in-vacuo annealing, but also focused ion beam exposure. The defect generation rate, atomic imaging and the optical signatures support this claim. We discriminate the narrow linewidth photoluminescence signatures of vacancies, resulting predominantly from localized defect orbitals, from broad luminescence features in the same spectral range, resulting from adsorbates. Vacancies can be patterned with a precision below 10 nm by ion beams, show single photon emission, and open the possibility for advanced defect engineering of 2D semiconductors at the ultimate scale.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Polarizabilities of neutral atoms and atomic ions with a noble gas electron configuration

Atomic polarizabilities play an important role in the development of force fields for molecular simulations, as well as for the development of qualitative concepts of atomic and molecular behavior. Coupled cluster theory at the coupled cluster singles doubles triples level with very large correlation-consistent basis sets with extended diffuse functions has been used to predict the polarizabilities of the atomic neutrals, mono-cations and mono-anions with a noble gas configuration. Additional corrections for scalar relativistic and spin–orbit effects were also included for the electron configurations of Kr, Xe, and Rn. The results are in excellent agreement with experiment or with other high level calculations where available. The current results for most of these species represent the best available values for the polarizabilities. The results show that the polarizability of H - is very difficult to calculate without extremely diffuse functions. The polarizability of H - is the largest value, 34.05 Å 3 , calculated for all species in the current study. The polarizabilities of the remaining halogen anions are also the best available values. The polarizabilities of the halogen anions (excluding F - ) and H - have a linear correlation with the electron affinity of the neutral atom. Spin– orbit effects, even for closed shell species, cannot be ignored for quantitative accuracy, and the inclusion of spin–orbit effects for Fr + , Rn, and At - increases the polarizability by 4%, 6%, and 15%, respectively.

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

Density-functional theory of material design: fundamentals and applications-I

This article is part-I of a review of density-functional theory (DFT) that is the most widely used method for calculating electronic structure of materials. The accuracy and ease of numerical implementation of DFT methods has resulted in its extensive use for materials design and discovery and has thus ushered in the new field of computational material science. In this article, we start with an introduction to Schrödinger equation and methods of its solutions. After presenting exact results for some well-known systems, difficulties encountered in solving the equation for interacting electrons are described. How these difficulties are handled using the variational principle for the energy to obtain approximate solutions of the Schrödinger equation is discussed. The resulting Hartree and Hartree–Fock theories are presented along with results they give for atomic and solid-state systems. We then describe Thomas–Fermi theory and its extensions which were the initial attempts to formulate many-electron problem in terms of electronic density of a system. Having described these theories, we introduce modern DFT by discussing Hohenberg–Kohn theorems that form its foundations. We then go on to discuss Kohn–Sham (KS) formulation of DFT in its exact form. Next, local density approximation (LDA) is introduced and solutions of KS equation for some representative systems, obtained using the LDA, are presented. We end part-I of the review describing the contents of part-II.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Defect-Promoted Ni-Based Layer Double Hydroxides with Enhanced Deprotonation Capability for Efficient Biomass Electrooxidation

Ni-based hydroxides are promising electrocatalysts for biomass oxidation reactions, supplanting the oxygen evolution reaction (OER) due to lower overpotentials while producing value-added chemicals. The identification and subsequent engineering of their catalytically active sites are essential to facilitate these anodic reactions. Herein, the proportional relationship between catalysts’ deprotonation propensity and Faradic efficiency of 5-hydroxymethylfurfural (5-HMF)-to-2,5 furandicarboxylic acid (FDCA, FEFDCA) is revealed by thorough density functional theory (DFT) simulations and atomic-scale characterizations, including in situ synchrotron diffraction and spectroscopy methods. The deprotonation capability of ultrathin layer-double hydroxides (UT-LDHs) is regulated by tuning the covalency of metal (M)-oxygen (O) motifs through defect site engineering and selection of M 3+ co-chemistry. NiMn UT-LDHs show an ultrahigh FE FDCA of 99% at 1.37 V versus reversible hydrogen electrode (RHE) and retain a high FE FDCA of 92.7% in the OER-operating window at 1.52 V, about 2× that of NiFe UT-LDHs (49.5%) at 1.52 V. Ni–O and Mn–O motifs function as dual active sites for HMF electrooxidation, where the continuous deprotonation of Mn–OH sites plays a dominant role in achieving high selectivity while suppressing OER at high potentials. The results showcase a universal concept of modulating competing anodic reactions in aqueous biomass electrolysis by electronically engineering the deprotonation behavior of metal hydroxides, anticipated to be translatable across various biomass substrates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Local p‐ and n‐Type Doping of an Oxide Semiconductor via Electric‐Field‐Driven Defect Migration

Layered oxides exhibit high ionic mobility and chemical flexibility, attracting interest as cathode materials for lithium-ion batteries and the pairing of hydrogen production and carbon capture. Recently, layered oxides emerged as highly tunable semiconductors. For example, by introducing anti-Frenkel defects, the electronic hopping conductance in hexagonal manganites is increased locally by orders of magnitude. Here, local acceptor and donor doping in Er(Mn,Ti)O3 is demonstrated, facilitated by the controlled splitting of anti-Frenkel defects under applied d.c. voltage. By combining density functional theory calculations, scanning probe microscopy, atom probe tomography, and scanning transmission electron microscopy, it is shown that the oxygen defects can readily be moved through the layered crystal structure, leading to nano-sized interstitial-rich (p-type) and vacancy-rich (n-type) regions. The resulting pattern is comparable to dipolar npn-junctions and stable on the timescale of days. These findings reveal the possibility of temporarily functionalizing oxide semiconductors at the nanoscale, giving additional opportunities for the field of oxide electronics and the development of transient electronics in general.

He, Jiali↗