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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
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 OSTI2020

Wannier–Koopmans method calculations for transition metal oxide band gaps

The widely used density functional theory (DFT) has a major drawback of underestimating the band gaps of materials. Wannier–Koopmans method (WKM) was recently developed for band gap calculations with accuracy on a par with more complicated methods. WKM has been tested for main group covalent semiconductors, alkali halides, 2D materials, and organic crystals. Here we apply the WKM to another interesting type of material system: the transition metal (TM) oxides. TM oxides can be classified as either with d 0 or d 10 closed shell occupancy or partially occupied open shell configuration, and the latter is known to be strongly correlated Mott insulators. We found that, while WKM provides adequate band gaps for the d0 and d10 TM oxides, it fails to provide correct band gaps for the group with partially occupied d states. This issue is also found in other mean-field approaches like the GW calculations. We believe that the problem comes from a strong interaction between the occupied and unoccupied d-state Wannier functions in a partially occupied d-state system. We also found that, for pseudopotential calculations including deep core levels, it is necessary to remove the electron densities of these deep core levels in the Hartree and exchange–correlation energy functional when calculating the WKM correction parameters for the d-state Wannier functions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accelerating multiscale electronic stopping power predictions with time-dependent density functional theory and machine learning

Knowing the rate at which particle radiation releases energy in a material, the “stopping power,” is key to designing nuclear reactors, medical treatments, semiconductor and quantum materials, and many other technologies. While the nuclear contribution to stopping power, i.e., elastic scattering between atoms, is well understood in the literature, the route for gathering data on the electronic contribution has for decades remained costly and reliant on many simplifying assumptions, including that materials are isotropic. We establish a method that combines time-dependent density functional theory (TDDFT) and machine learning to reduce the time to assess new materials to hours on a supercomputer and provide valuable data on how atomic details influence electronic stopping. Our approach uses TDDFT to compute the electronic stopping from first principles in several directions and then machine learning to interpolate to other directions at a cost of 10 million times fewer core-hours. We demonstrate the combined approach in a study of proton irradiation in aluminum and employ it to predict how the depth of maximum energy deposition, the “Bragg Peak,” varies depending on the incident angle—a quantity otherwise inaccessible to modelers and far outside the scales of quantum mechanical simulations. The lack of any experimental information requirement makes our method applicable to most materials, and its speed makes it a prime candidate for enabling quantum-to-continuum models of radiation damage. The prospect of reusing valuable TDDFT data for training the model makes our approach appealing for applications in the age of materials data science.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Pursuing Dark Energy with Large Galaxy Redshift Surveys: Baryon Acoustic Oscillations and Beyond (Final Technical Report)

We have developed and applied methods of large-scale structure from galaxy redshift surveys to the study of dark energy. This program supported my group's development of the baryon acoustic method and our participation in the Dark Energy Spectroscopic Instrument (DESI) survey. In addition to my leadership role in the DESI collaboration, my group at Harvard focused on several topics addressing key needs of the DESI science program. 1) We continued our design, validation, and documentation of DESI target selection of the luminous red galaxy and emission-line galaxy samples. 2) We applied simulations from our novel and extremely fast N-body code to make DESI catalogs. With these, we will search for possible systematic shifts in the acoustic scale due to clustering bias with an unprecedented volume of simulations. 3) We extended our method for the fast construction of covariance matrices for correlation functions. 4) We developed new large-scale structure methods, including density-field reconstruction and three-point function analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Cu–Ni Oxidation Mechanism Unveiled: A Machine Learning-Accelerated First-Principles and in Situ TEM Study

Here, the development of accurate methods for determining how alloy surfaces spontaneously restructure under reactive and corrosive environments is a key, long-standing, grand challenge in materials science. Using machine learning-accelerated density functional theory and rare-event methods, in conjunction with in situ environmental transmission electron microscopy (ETEM), we examine the interplay between surface reconstructions and preferential segregation tendencies of CuNi(100) surfaces under oxidation conditions. Our modeling approach predicts that oxygen-induced Ni segregation in CuNi alloys favors Cu(100)-O c(2 × 2) reconstruction and destabilizes the Cu(100)-O (2√2 × √2)R45° missing row reconstruction (MRR). In situ ETEM experiments validate these predictions and show Ni segregation followed by NiO nucleation and growth in regions without MRR, with secondary nucleation and growth of Cu 2 O in MRR regions. Our approach based on combining disparate computational components and in situ ETEM provides a holistic description of the oxidation mechanism in CuNi, which applies to other alloy systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Theoretical Investigation of the Adsorbate and Potential–Induced Stability of Cu Facets During Electrochemical CO 2 and CO Reduction

The activity and product selectivity of electrocatalysts for reactions like the carbon dioxide reduction reaction (CO 2 RR) are intimately dependent on the catalyst's structure and composition. While engineering catalytic surfaces can improve performance, discovering the key sets of rational design principles remains challenging due to limitations in modeling catalyst stability under operating conditions. Herein, we perform first-principles density functional calculations adopting implicit solvation methods with potential control to study the influence of adsorbates and applied potential on the stability of different facets of model Cu electrocatalysts. Using coverage dependencies extracted from microkinetic models, we describe an approach for calculating potential and adsorbate-dependent contributions to surface energies under reaction conditions, where Wulff constructions are used to understand the morphological evolution of Cu electrocatalysts under CO 2 RR conditions. Here we identify that CO*, a key reaction intermediate, exhibits higher kinetically and thermodynamically accessible coverages on (100) relative to (111) facets, which can translate into an increased relative stabilization of the (100) facet during CO 2 RR. Our results support the known tendency for increased (111) faceting of Cu nanoparticles under more reducing conditions and that the relative increase in (100) faceting observed under CO 2 RR conditions is likely attributed to differences in CO* coverage between these facets.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

Electron transfer calculations between edge sharing octahedra in hematite, goethite, and annite

A key reaction underlying the charge transport in iron containing oxides, clays, micas is the Fe$^{2+}$-Fe$^{3+}$ exchange reaction between edge-sharing iron octahedra. These reactions facilitate conduction in these minerals by the thermally-activated hopping of small polarons across the lattice. Depending on the mineral and local charge state the small polaron can either encase an electron or hole. The probability for conduction of small polarons depends strongly on the height and adiabicity of the reaction barrier, with larger and more diabatic barriers yielding slow conduction associated with either weak coupling or a large prerequisite rearrangement of the lattice during charge transport. To model these reactions, a first principle electron transfer (ET) method was developed to model the small polaron hopping between the edge-sharing octahedra sites in hematite ($e^{-}$ polaron), goethite ($e^{-}$ polaron), and annite ($h^{+}$ polaron) bulk structures. The ET method is based on electronic structure methods (i.e., plane-wave Density Functional Theory) capable of performing calculations with periodic cells and large size systems efficiently while at the same time being accurate enough to be used in the estimation of the electron-transfer coupling matrix element, $V_{AB}$, and the electron transfer transmission factor, $\kappa_{el}$. Additionally, the calculations confirmed the existence of small polarons in all three minerals, and the reactions were predicted to be strongly adiabatic. It was found that transfer of a hole in the octahedral layer of annite had an adiabatic barrier of $0.311$ eV, and the transfer of an extra electron in hematite and goethite had adiabatic barriers of $0.242$ eV and $0.232$ eV respectively. The electronic coupling parameters,$V_{AB}$, were found to be $0.188$ eV, $0.196$ eV, and $0.102$ eV respectively for hematite, goethite, and annite. While similar bonding topologies pertain, the underlying basis for the differences is the subtle differences in local structures.

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

Effects of charged interfaces on electrolyte decomposition at the lithium metal anode

Lithium–Sulfur batteries are promising candidates to substitute conventional Li-ion batteries due to their higher energy density and reduced cost. However, several challenges related to the reactivity of the lithium metal anode have prevented this technology from becoming broadly commercialized. Lithium's high reactive nature leads to the continuous decomposition of the electrolyte and the formation of the solid-electrolyte interphase (SEI) layer. Thus, a comprehensive understanding of how the SEI film is formed is crucial to help elucidate improvements for this battery technology. In this work, we use density functional theory (DFT) based computational methods to investigate the effect of charged interfaces on the electrolyte degradation at the lithium anode. Here, several electrolyte mixtures were considered including 1,2-dimethoxyethane (DME) and 1,3-dioxolane (DOL) as solvents, and lithium bis(fluorosulfonyl)imide (LiFSI) and lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) as salts. It is found that the extent of salt decomposition is higher when the interface is charged. In addition, two types of solvent reduction mechanisms are identified: C–O bond cleavage and radical attacks. Charge transfer and evolution analysis are characterized in detail. Finally, solvent decomposition reactions are found to become more thermodynamically favorable and their activation barriers to be diminished under the effect of constant potentials and electric fields.

25 ENERGY STORAGE↗
DOE OSTI2020

A Generalized Variational Principle with Applications to Excited State Mean Field Theory

We present a generalization of the variational principle that is compatible with any Hamiltonian eigenstate that can be specified uniquely by a list of properties. This variational principle appears to be compatible with a wide range of electronic structure methods, including mean field theory, density functional theory, multireference theory, and quantum Monte Carlo. Like the standard variational principle, this generalized variational principle amounts to the optimization of a nonlinear function that, in the limit of an arbitrarily flexible wave function, has the desired Hamiltonian eigenstate as its global minimum. Unlike the standard variational principle, it can target excited states and select individual states in cases of degeneracy or near-degeneracy. As an initial demonstration of how this approach can be useful in practice, we employ it to improve the optimization efficiency of excited state mean field theory by an order of magnitude. With this improved optimization, we are able to demonstrate that the accuracy of the corresponding second-order perturbation theory rivals that of singles-and-doubles equation-of-motion coupled cluster in a substantially broader set of molecules than could be explored by our previous optimization methodology.

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

2D-Covalent Organic Frameworks with Interlayer Hydrogen Bonding Oriented through Designed Nonplanarity

We report the synthesis and characterization of a new class of 2D-covalent organic frameworks, called COFamides, whose layers are held together by amide hydrogen bonds. To accomplish this, we have designed monomers with a nonplanar structure that arises from steric crowding, forcing the amide side groups out of plane with the COF sheets orienting the hydrogen bonds between the layers. Here, the presence of these hydrogen bonds provides significant structural stabilization as demonstrated by comparison to control structures that lack hydrogen bonding capability, resulting in lower surface area and crystallinity. We have characterized both azine and imine-linked versions of these COFs, named COFamide-1 and -2, respectively, for their surface areas, pore sizes, and crystallinity. In addition to these more conventional characterization methods, we also used variable temperature infrared spectroscopy methods and van der Waals density functional calculations to directly observe the presence of hydrogen bonding.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electromagnetic properties of indium isotopes illuminate the doubly magic character of 100 Sn

Understanding the nuclear properties in the vicinity of 100 Sn, which has been suggested to be the heaviest doubly magic nucleus with proton number Z equal to neutron number N, has been a long-standing challenge for experimental and theoretical nuclear physics. In particular, contradictory experimental evidence exists regarding the role of nuclear collectivity in this region of the nuclear chart. Here, we provide further evidence for the doubly magic character of 100 Sn by measuring the ground-state electromagnetic moments and nuclear charge radii of indium (Z = 49) isotopes as N approaches 50 from above using precision laser spectroscopy. Our results span almost the complete range between the two major closed neutron shells at N = 50 and N = 82 and reveal parabolic trends as a function of the neutron number, with a clear reduction towards these two closed neutron shells. Furthermore, a detailed comparison between our experimental results and numerical results from two complementary nuclear many-body frameworks (density functional theory and ab initio methods) exposes deficiencies in nuclear models and establishes a benchmark for future theoretical developments.

100Sn↗
DOE OSTI2025

Mechanistic insights into rubidium ion adsorption at the quartz (101) surface from quantum chemical metadynamics

Ion sorption extent and mechanism depend in part on the mineral surface termination, which can be highly complex. Variations in surface functional groups, particularly with defect density and surface roughness, influence mineral reactivity towards solutes. In this work, we investigate the adsorption of a rubidium cation (Rb + ) at pristine and defect quartz (101) surface sites using well-tempered metadynamics, based on simulations with the quantum chemical density-functional tight-binding (DFTB) method. We compare the relative energetics of Rb + adsorption across selected sites for each surface, with nanosecond-level sampling, highlighting similarities between vicinal and geminal silanol sites. We find that the positive Rb + partial atomic charge can increase by as much as 0.5e as it approaches the surface, with implications for modulation of ion adsorption strength and extent at the quartz (101) surface with surface vacancies, silanol coverage, and charge.

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

Efficacy of the radial pair potential approximation for molecular dynamics simulations of dense plasmas

Macroscopic simulations of dense plasmas rely on detailed microscopic information that can be computationally expensive and is difficult to verify experimentally. In this work, we delineate the accuracy boundary between microscale simulation methods by comparing Kohn–Sham density functional theory molecular dynamics (KS-MD) and radial pair potential molecular dynamics (RPP-MD) for a range of elements, temperature, and density. By extracting the optimal RPP from KS-MD data using force matching, we constrain its functional form and dismiss classes of potentials that assume a constant power law for small interparticle distances. Our results show excellent agreement between RPP-MD and KS-MD for multiple metrics of accuracy at temperatures of only a few electron volts. The use of RPPs offers orders of magnitude decrease in computational cost and indicates that three-body potentials are not required beyond temperatures of a few eV. Due to its efficiency, the validated RPP-MD provides an avenue for reducing errors due to finite-size effects that can be on the order of ∼20%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Multiple shape coexistence near {sup 118}Sn: First 0+ 3 lifetime measurement

The intruder bands in Sn isotopes, built on the 2p-2h excitation across the Z = 50 proton shell gap, are well-known examples of shape coexistence near the neutron mid-shell region. Spectroscopic signatures for shape coexistence include enhanced E0 transitions between the 0+ band heads. However, the underlying shape coexistence and mixing has been unclear because lifetime information for the excited 0+ states was incomplete in 118Sn. We thus present here the first measurement of the 0+ 3 lifetime in 118Sn using the fast-timing technique following thermal-neutron capture. The observed enhancement in ρ2(E0;0+ 3 → 0+ 2 ) of 150(30) milliunits provides compelling indications for multiple shape coexistence in 118Sn. Additionally, three distinct shapes in 116,118,120Sn naturally emerged in theoretical calculations based on the quantum-number-projected generator coordinate method employing a relativistic energy density functional.

Wu, F.↗
DOE OSTI2020

Analysis of diagonal G and subspace W approximations within fully self-consistent GW calculations for bulk semiconducting systems

Fully self-consistent GW (sc-GW) methods are now available to evaluate quasiparticle and spectral properties of various molecular and bulk systems. However, such techniques based on the full matrix of G and W are computationally demanding. Additionally, the routinely used single-shot GW approximation (G 0 W 0 ) has an undesirable dependency on the choice of initial exchange-correlation functional. In the literature, many so-called self-consistent GW methods are based on diagonal approximation of G and low-ranking approximation of W. It is thus worth checking how good such approximations are in comparison with the full matrix method. In this work, we consider AlAs, AlP, GaP, and ZnS as the prototype systems to perform sc-GW calculations by expressing the full G matrix using a plane-wave basis set. We compared our sc-GW results with the diagonal G and subspace W approximated sc-GW results (sc-GW-diagG and sc-GW-subW methods). In the sc-GW-diagG method, interacting G is expanded in the eigenvectors of noninteracting G such that only diagonal elements are retained, whereas the number of eigenmodes is truncated in sc-GW-subW calculations. A systematic analysis of the results obtained from the above techniques is presented. The differences in the quasiparticle band gap between the approximated and the full matrix sc-GW approaches are mostly less than 1.7%, which validates such widely adopted approximations, and also shows how such low-ranking approximation can be used to include higher-order terms such as the vertex correction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Phonon transport governed by intrinsic scattering in short-period AlN/GaN superlattices

Here we employ density functional theory based phonon transport methods to provide a rigorous understanding of the nature of thermal transport in coherent short-period AlN/GaN superlattices (SLs), with period lengths up to three unit cells of each, and compare these with properties of their bulk constituents. Increasing the period length leads to phonon band folding with frequency gaps and thus smaller phonon velocities and more phonon scattering than in bulk, both of which reduce lattice thermal conductivity (κ). Contrary to expectations, we find that velocity variations among larger-period AlN/GaN SLs play only a minor role in cross-plane κ reductions, while variations in intrinsic phonon scattering are strongly correlated with their transport behaviors. This work provides insights into the microscopic behaviors of technologically relevant nanostructured materials, which are likely applicable to a wider range of SL systems and other nanostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Time-dependent density functional theory of charge, energy, and spin transport, and quantum continuum mechanics (Final Report)

The overarching goal of this project, which unfolded over a period of seventeen years, including two no-cost extensions in the last two years, was to extend and deepen our ability to use time-dependent density functional theory and related theoretical methods for the description of the transport of charge, spin, or energy in quantum electronic systems at the nanoscale. The output of our efforts is contained in 52 papers, listed below, which we have published since 2006. There have also been many presentations at international conferences, workshops, colloquia and seminars. This body of work has been carried out with the essential input of 8 postdoctoral students (who were supported by the grant) and several external senior collaborators and their postdocs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗