Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY FUNCTIONAL METHOD”

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

Ab initio Study of the Mechanism for the Initiation of Localized Corrosion on Stainless Steel

Stainless steels are widely used as primary piping material in PWR’s because of their high corrosion resistance largely due to the formation of a passive film of metal oxides, such as chromium (e.g., Cr2O3). However, upon exposure to aggressive ions like chlorides (Cl-), stainless steel becomes susceptible to localized corrosion due to the breakdown of the passive film. It is therefore essential to gain fundamental insights into the mechanism of the passive film breakdown (depassivation) to develop a tailored design to minimize localized corrosion. This project investigates the initial stages of the depassivation of the Cr passive film using density functional theory (DFT), a quantum based computational method.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phonon-Assisted Ballistic Current from First-Principles Calculations

The bulk photovoltaic effect (BPVE) refers to current generation due to illumination by light in a homogeneous bulk material lacking inversion symmetry. In addition to the intensively studied shift current, the ballistic current, which originates from asymmetric carrier generation due to scattering processes, also constitutes an important contribution to the overall kinetic model of the BPVE. In this Letter, we use a perturbative approach to derive a formula for the ballistic current resulting from the intrinsic electron-phonon scattering in a form amenable to first-principles calculation. We then implement the theory and calculate the ballistic current of the prototypical BPVE material BaTiO 3 using quantum-mechanical density functional theory. The magnitude of the ballistic current is comparable to that of the shift current, and the total spectrum (shift plus ballistic) agrees well with the experimentally measured photocurrents. Furthermore, we show that the ballistic current is sensitive to structural change, which could benefit future photovoltaic materials design.

14 SOLAR ENERGY↗
DOE OSTI2025

Advances and Challenges of SCAN and r 2 SCAN Density Functionals in Transition‐Metal Compounds

ABSTRACT Transition‐metal compounds (TMCs) with open‐shell d ‐electrons are characterized by a complex interplay of lattice, charge, orbital, and spin degrees of freedom, giving rise to various fascinating applications. Often exhibiting exotic properties, these compounds are commonly classified as correlated systems due to strong inter‐electronic interactions called Hubbard U . This inherent complexity presents significant challenges to Kohn‐Sham density functional theory (KS‐DFT), the most widely used electronic structure method in condensed matter physics and materials science. While KS‐DFT is, in principle, exact for the ground‐state total energy, its exchange‐correlation energy must be approximated in practice. The mean‐field nature of KS implementations, combined with the limitations of current exchange‐correlation density functional approximations, has led to the perception that DFT is inadequate for correlated systems, particularly TMCs. Consequently, a common workaround involves augmenting DFT with an on‐site Hubbard‐like U correction. In recent years, the strongly constrained and appropriately normed (SCAN) density functional, along with its refined variant r 2 SCAN, has achieved remarkable progress in accurately describing the structural, energetic, electronic, magnetic, and vibrational properties of TMCs, challenging the traditional perception of DFT's limitations. This review explores the design principles of SCAN and r 2 SCAN, highlights their key advancements in studying TMCs, explains the mechanisms driving these improvements, and addresses the remaining challenges in this evolving field.

Zhang, Yubo [Minjiang Collaborative Center for The↗
DOE OSTI2026

Enhancing the efficiency of time-dependent density functional theory calculations of dynamic response properties

X-ray Thomson scattering (XRTS) constitutes an essential technique for diagnosing material properties under extreme conditions, such as high pressures and intense laser heating. Time-dependent density functional theory (TDDFT) is one of the most accurate available ab initio methods for modeling XRTS spectra, as well as a host of other dynamic material properties. However, strong thermal excitations, along with the need to account for variations in temperature and density as well as the finite size of the detector significantly increase the computational cost of TDDFT simulations compared to ambient conditions. In this work, we present a broadly applicable method for optimizing and enhancing the efficiency of TDDFT calculations. Our approach is based on a one-to-one mapping between the dynamic structure factor and the imaginary time density–density correlation function, which naturally emerges in Feynman’s path integral formulation of quantum many-body theory. Specifically, we combine rigorous convergence tests in the imaginary time domain with a constraints-based attenuation of narrow-band fluctuations to improve the efficiency of TDDFT modeling without the introduction of any significant bias. As a result, we can report a speed-up by up to an order of magnitude, thus substantially reducing the burden of computational cost required for XRTS analysis.

Moldabekov, Zhandos A. [Helmholtz-Zentrum Dresden-↗
DOE OSTI2023

Hydrogen irradiation-driven computational surface chemistry of lithium oxide and hydroxide

We have investigated, using molecular dynamics, the surface chemistry of hydrogen incident on the amorphous and crystalline lithium oxide and lithium hydroxide surfaces upon being slowed down by a collision cascade and retained in the amorphous surface of either Li 2 O or LiOH. We looked for the bonding of H to the resident structures in the surface to understand a possible chain of chemical reactions that can lead to surface transformation upon H atom impact. Our findings, using Density-Functional Theory (DFT) trained ReaxFF force field/electronegativity equalization method potentials, stress the importance of inclusion of polarization in the dynamics of a Li–O–H system, which is also illustrated by DFT energy minimization and quantum–classical molecular dynamics using tight binding DFT. The resulting polar-covalent chemistry of the studied systems is complex and very sensitive to the instantaneous positions of all atoms as well as the ratio of concentrations of various resident atoms in the surface.

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

Quantum Monte Carlo benchmarking of large noncovalent complexes in the L7 benchmark set

We have used diffusion Monte Carlo (DMC) to perform calculations on the L7 benchmark set. DMC is a stochastic numerical integration scheme in real-space and part of a larger set of quantum Monte Carlo methods. The L7 set was designed to test the ability of electronic structure methods to include dispersive interactions. While the agreement between DMC and quantum-chemical state-of-the-art methods is excellent for some of the structures, there are significant differences in others. In contrast to wavefunction-based quantum chemical methods, DMC is a first-principle many-body method with the many-body wavefunction evolving in real space. It includes explicitly all electron–electron interactions and is relatively insensitive to the size of the basis set.

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

Water Catalysis of the Reaction of Methanol with OH Radical in the Atmosphere is Negligible

Abstract Faced with the contradictory results of two recent experimental studies [Jara‐Toro et al., Angew. Chem. Int. Ed . 2017 , 56 , 2166 and Chao et al., Angew. Chem. Int. Ed . 2019 , 58 , 5013] of the possible catalytic effect of water vapor on CH 3 OH + OH reaction, we report calculations that corroborate the conclusion made by Chao et al. and extend the rate constant evaluation down to 200 K. The rate constants of the CH 3 OH + OH reaction catalyzed by a water molecule are computed as functions of temperature and relative humidity using high‐level electronic structure and kinetics calculations. The Wuhan–Minnesota Scaling (WMS) method is used to provide accurate energetics to benchmark a density functional for direct dynamics. Both high‐frequency and low‐frequency anharmonicities are included. Variational and tunneling effects are treated by canonical variational transition state theory with multidimensional small‐curvature tunneling. And, most significantly, we include multistructural effects in the rate constant calculations. Our calculations show that the catalytic effect of water vapor is not observable at 200–400 K.

Wu, Junjun↗
DOE OSTI2020

Predictive Model of Charge Mobilities in Organic Semiconductor Small Molecules with Force-Matched Potentials

Charge mobility of crystalline organic semiconductors (OSC) is limited by local dynamic disorder. Recently, the charge mobility for several high mobility OSCs, including TIPS-pentacene, were accurately predicted from a density functional theory (DFT) simulation constrained by the crystal structure and the inelastic neutron scattering spectrum, which provide direct measures of the structure and the dynamic disorder in the length scale and energy range of interest. However, the computational expense required for calculating all of the atomic and molecular forces is prohibitive. In this work, we demonstrate the use of density functional tight binding (DFTB), a semiempirical quantum mechanical method that is 2 to 3 orders of magnitude more efficient than DFT. We show that force matching a many-body interaction potential to DFT derived forces yields highly accurate DFTB models capable of reproducing the low-frequency intricacies of experimental inelastic neutron scattering (INS) spectra and accurately predicting charge mobility. We subsequently predicted charge mobilities from our DFTB model of a number of previously unstudied structural analogues to TIPS-pentacene using dynamic disorder from DFTB and transient localization theory. The approach we establish here could provide a truly rapid simulation pathway for accurate materials properties prediction, in our vision applied to new OSCs with tailored properties.

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

Molecular-based analysis of nanoparticle solvation: Classical density functional approach

Proper statistical mechanics understanding of nanoparticle solvation processes requires an accurate description of the molecular structure of the solvent. Achieving this goal with standard molecular dynamics (MD) simulation methods is challenging due to large length scales. An alternative approach to this problem can be formulated using classical density functional theory (cDFT), where a full configurational description of the positions of all the atoms is replaced by collective atomic site densities in the molecule. Using an example of the negatively charged silica-like system in an aqueous polar environment represented by a two-site water model, we demonstrate here that cDFT can reproduce MD data at a fraction of the computational cost. An important implication of this result is the ability to understand how the solvent molecular features may affect the system’s properties at the macroscopic scale. A concrete example highlighted in this work is the analysis of nanoparticle interactions with sizes of up to 100 nm in diameter.

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

Multiscale computational understanding and growth of 2D materials: a review

Abstract The successful discovery and isolation of graphene in 2004, and the subsequent synthesis of layered semiconductors and heterostructures beyond graphene have led to the exploding field of two-dimensional (2D) materials that explore their growth, new atomic-scale physics, and potential device applications. This review aims to provide an overview of theoretical, computational, and machine learning methods and tools at multiple length and time scales, and discuss how they can be utilized to assist/guide the design and synthesis of 2D materials beyond graphene. We focus on three methods at different length and time scales as follows: (i) nanoscale atomistic simulations including density functional theory (DFT) calculations and molecular dynamics simulations employing empirical and reactive interatomic potentials; (ii) mesoscale methods such as phase-field method; and (iii) macroscale continuum approaches by coupling thermal and chemical transport equations. We discuss how machine learning can be combined with computation and experiments to understand the correlations between structures and properties of 2D materials, and to guide the discovery of new 2D materials. We will also provide an outlook for the applications of computational approaches to 2D materials synthesis and growth in general.

Momeni, Kasra↗
DOE OSTI2020

Final Technical Report for project entitled Highly Active, Durable, and Ultra-Low PGM NSTF Thin Film ORR Catalysts and Support

In this project, the objective was to develop new oxygen reduction reaction (ORR) electrocatalysts for proton exchange membrane fuel cells (PEMFCs) which could exceed all of the Department of Energy (DOE) 2020 targets listed in DE-FOA-0001224, Subtopic 1b, Table I. The expected outcome was development of one or more electrocatalysts which are substantially improved over the current state-of-the-art in terms of overall activity, durability, and cost, and are suitable for automotive traction and stationary fuel cell applications. This report summarizes this project’s progress towards meeting the stated objectives and includes a summary of key findings and conclusions. The development activities towards new highly active and durable thin film electrocatalysts has led to the discovery of several electrocatalysts which approached or exceeded several DOE 2020 targets for activity, durability, and performance in proton exchange membrane fuel cell (membrane electrode assemblies (MEAs). In this project, electrocatalyst development focused on systematic physical and electrochemical characterization of electrocatalyst activity, durability, and performance in MEAs as a function of electrocatalyst fabrication, compositional and structural variables. The project approach was to develop relationships between the catalyst functional responses (activity, durability and performance) with catalyst physical properties and catalyst fabrication methods. Additionally, the thin film catalysts were integrated onto the unique and durable 3M Nanostructured Thin Film (NSTF) support, consisting of arrays of self-assembled organic crystalline whiskers. The catalysts were evaluated against the project targets via extensive electrochemical characterization in MEA format at 3M, advanced structural and compositional microscopy at Oak Ridge National Laboratory (ORNL), and atomic structure analysis via XAFS at Argonne National Laboratory (ANL). The extensive electrochemical and physical characterization resulted in development of several trends which correlate the catalysts’ electrochemical properties to their physical properties, including composition, structure, and method of fabrication. The development was guided by density functional theory (DFT) modeling at Purdue University and kinetic Monte Carlo (kMC) modeling at Johns Hopkins University. The catalyst simulations provided key insights into the observed experimental activity and durability trends, and additionally were utilized to assess new electrocatalyst concepts prior to or coincident with physical catalyst development. This development has led to several electrocatalyst candidates with activity and durability which approach or exceed DOE targets. One class of catalysts, based on nm-scale thin layers of Pt on Ir, met or exceeded 6 of the 6 DOE targets the project addressed.

08 HYDROGEN↗
DOE OSTI2020

Beyond the proton drip line: Bayesian analysis of proton-emitting nuclei

The limits of the nuclear landscape are determined by nuclear binding energies. Beyond the proton drip lines, where the separation energy becomes negative, there is not enough binding energy to prevent protons from escaping the nucleus. Predicting properties of unstable nuclear states in the vast territory of proton emitters poses an appreciable challenge for nuclear theory as it often involves far extrapolations. In addition, significant discrepancies between nuclear models in the proton-rich territory call for quantified predictions. With the help of Bayesian methodology, we mix a family of nuclear mass models corrected with statistical emulators trained on the experimental mass measurements, in the proton-rich region of the nuclear chart. Separation energies were computed within nuclear density functional theory using several Skyrme and Gogny energy density functionals. We also considered mass predictions based on two models used in astrophysical studies. Quantified predictions were obtained for each model using Bayesian Gaussian processes trained on separation-energy residuals and combined via Bayesian model averaging. Here, we obtained a good agreement between averaged predictions of statistically corrected models and experiment. In particular, we quantified model results for one- and two-proton separation energies and derived probabilities of proton emission. This information enabled us to produce a quantified landscape of proton-rich nuclei. The most promising candidates for two-proton decay studies have been identified. Finally, the methodology used in this work has broad applications to model-based extrapolations of various nuclear observables. It also provides a reliable uncertainty quantification of theoretical predictions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Primary $\gamma$-ray intensities and $\gamma$-strength functions from discrete two-step $\gamma$-ray cascades in radiative proton-capture experiments

Background: Reaction rates of radiative capture reactions can play a crucial role in the nucleosynthesis of heavy nuclei in explosive stellar environments. These reaction rates depend strongly on $\gamma$-ray decay widths in the reaction products, which are, for nonresonant capture reactions at high excitation energies, derived from the $\gamma$-ray strength function and the nuclear level density. Recently, the ratio method was applied to primary $\gamma$ rays observed from ($\textit{d, p}$) reactions and nuclear resonance fluorescence measurements to extract the dipole strength in atomic nuclei and to test the generalized Brink-Axel hypothesis. Purpose: The purpose of this work is to apply the ratio method to primary $\gamma$-ray intensities of the 63,65 Cu($\textit{p},\gamma$) reactions to extract $\gamma$-ray strength information on the nuclei 64,66 Zn. Here, the impact of spin distribution, total $\gamma$-ray decay widths, level densities, and width fluctuations on the application of the ratio method will be discussed. Additionally, by comparing the relative $\gamma$-ray strength at different excitation energies, conclusions on the validity of the generalized Brink-Axel hypothesis can be made. Method: The radiative proton capture reaction measurements have been performed at the HORUS $\gamma$-ray spectrometer of the University of Cologne at one excitation energy for each reaction. Primary $\gamma$-ray intensities have been determined by normalizing secondary $\gamma$-ray transitions in two-step cascades using their absolute branching ratio. The ratio method was applied to the measured primary $\gamma$-ray intensities as well as to previous measurements by Erlandsson et al. at different excitation energies. Results: The relative strength function curve for 64 Zn from our measurement shows no significant deviation from the previous measurement at a different excitation energy. The same is true for 66 Zn where both measurements were at almost the same excitation energy. Absolute $\gamma$-strength function values have been obtained by normalizing the relative curves to quasiparticle random phase approximation calculations because of the absence of experimental data in the respective energy region. Conclusion: The generalized Brink-Axel hypothesis, i.e., the independence of the strength function on the excitation energy, seems to hold in the studied energy region and nuclei. The method to obtain primary $\gamma$-ray intensities from two-step cascade spectra was shown to be a valuable and sensitive tool although its uncertainties are connected to the knowledge of the low-energy level scheme of the investigated nucleus. The scaling in the ratio method should be taken with care, because the relative strength is not a simple sum of $f_{E1}$ and $f_{M1}$ but a somewhat complex linear combination dependent on the excitation energy of the nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Observation of halogen-like behavior of gold in fluorinated bimetallic CoAuF 1-2 − and CuAuF 1-2 − clusters: Anion photoelectron spectroscopy and density functional theory

Using size-selected anion photoelectron spectroscopy and density functional theory, we investigated the structures and properties of fluorinated bimetallic clusters CoAuF 1-2 − and CuAuF 1-2 − and their neutrals. Both experimental and theoretical results show that in these cluster anions, Au behaves like a halogen atom. For example, the measured vertical detachment energies (VDEs) of CoAuF − (2.00 ± 0.08 eV) and CuAuF − (3.8 ± 0.1 eV) are close to those of CoF 2 − (2.12 ± 0.08 eV) and CuF 2 − (3.58 ± 0.08 eV), respectively. The theoretical results show that the geometries and electronic structures of CoAuF − and CuAuF − are similar to those of CoF 2 − and CuF 2 − . The natural population analysis and natural electron configuration analyses further confirm that the electronic properties of Au in MAuF − (M = Co, Cu) mimic those of MF 2 − . In addition, the electron localization function analyses show that the M-Au chemical bonds are similar to the corresponding M-F chemical bonds, providing evidence for the ionic nature of the interactions. When a second F atom is attached to the CoAuF − and CuAuF − clusters, the VDEs of the resulting CoAuF 2 − and CuAuF 2 − are 4.38 ± 0.08 eV and 3.71 ± 0.08 eV, respectively, indicating their superhalogen character as these values are higher than those of halogen anions. The results may be useful for understanding the properties of gold at the nanoscale that play an important role in catalysis and nanotechnology.

Ab-initio methods↗
DOE OSTI2022

Phonons of metallic hydrogen with quantum Monte Carlo

Here, we describe a simple scheme to perform phonon calculations with quantum Monte Carlo (QMC) methods and demonstrate it on metallic hydrogen. Because of the energy and length scales of metallic hydrogen and the statistical noise inherent to QMC methods, the conventional manner of calculating force constants is prohibitively expensive. We show that our alternate approach is nearly 100 times more efficient in resolving the force constants needed to calculate the phonon spectrum in the harmonic approximation. This requires only the calculation of atomic forces, as in the conventional approach, and otherwise little or no programmatic modification.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗