Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

LaMnO 3 Dopants for Efficient Thermochemical Water Splitting Identified by Density Functional Theory Calculations

While ceria is the standard material for two-step water splitting, perovskites are emerging as viable alternatives. In this work, based on the orthorhombic LaMnO 3 supercell, we substitute Li Na K Rb Mg Ca Sr Ba on the A-sites (La sites) and Al Ga In Mg Zn on the B-sites (Mn sites) at a concentration of 37.5%. The range of temperature and oxygen partial pressure at which each composition is stable is predicted. For compositions that are stable in relevant temperature and pressure ranges, the oxygen vacancy formation energies are determined for all of the oxygen vacancy site positions available in the computational supercell. Mg, Ca, Sr, and Ba A-site-substituted LaMnO 3 and Al and In B-site-substituted LaMnO 3 meet these two criteria for candidates in solar-thermal water splitting applications. Finally, oxygen vacancy formation energy can also be controlled by adjusting the doping strategy.

08 HYDROGEN↗
DOE OSTI2023

Feasibility of a Novel Density Functional Method Outside the Kohn-Sham Framework for Modeling Global Potential Energy Surfaces of Molecular Chemical Reactions (Final Technical Report)

The project aimed towards the construction of computational methods capable of modeling the global potential energy surface (PES) of small molecules, molecular ions, and radicals — including the parts of the PES which correspond to chemical reactions, and the reaction paths connecting reactants, intermediates, and reaction products. This goal may seem humble at first glance, but for chemical systems with more than about six atoms in total, at the time the project was proposed, there were no established theoretical methods capable of simulating such systems reliably, not even for small molecules in the gas phase which are electronically benign. This restriction severely hampers our ability to model and control chemical processes under harsh conditions. With the goal of constructing a method capable of modeling such chemical systems, we proposed to pursue a novel approach towards a Multi-Configuration (MC) DFT outside the traditional frameworks of Kohn-Sham theory and other methods of coupling wave function theory with DFT. Rather than being a complete active space (CAS) method, the proposed DOCI-DFT would employ a special restricted form of the active space wave function, called Doubly-Occupied CI (DOCI)—this wave function form is sufficient to describe not only heterolytic, but also homolytic bond dissociation processes at the zeroth order (i.e., as active space wave function). There are no other standard mean field methods which can do so. In the original proposal, we also outlined strong formal and practical arguments speaking for this method. The project proposal was accepted by the Department of Energy and provided two years of seed funding for one graduate student, as well as two years of two weeks PI summer-salary for the PI. However, ultimately the project could not be effectively pursued due to severe interferences outside the context scientific problems and was cancelled. Apart from two weeks of PI salary (and associated fringe & overhead costs issues by the performing institution), no costs were charged against the award, with the entire rest of funding returned to the Department of Energy. This report reiterates the primary background information regarding the project, its original goals, and outlines the preliminary work performed during the two weeks of DOE funding charged.

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

Global properties of nuclei at finite-temperature within the covariant energy density functional theory

In stellar environments nuclei appear at finite temperatures, becoming extremely hot in core-collapse supernovae and neutron-star mergers. However, due to theoretical and computational complexity, most model calculations of nuclear properties are performed at zero temperature, while those existing at finite temperatures are limited only to selected regions of the nuclide chart. In this study we perform the global calculation of nuclear properties for even-even 8 ≤ Z ≤ 104 nuclei at temperatures in range 0 ≤ T ≤ 2 MeV. Calculations are based on the finite-temperature relativistic Hartree-Bogoliubov model supplemented by the Bonche-Levit-Vautherin vapor subtraction procedure. We find that near the neutron-drip line the continuum states have significant contribution already at moderate temperature T ≈ 1 MeV, thus emphasizing the necessity of the vapor subtraction procedure. Results include neutron emission lifetimes, quadrupole deformations, neutron-skin thickness, proton and neutron pairing gaps, entropy and excitation energy. Up to the temperature T ≈ 1 MeV, the nuclear landscape is influenced only moderately by the finite-temperature effects, mainly by reducing the pairing correlations. Here, as the temperature increases further, the effects on nuclear structures become pronounced, reducing both the deformations and the shell effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Density Functional Theory Study of the Initial Stages of Cl-Induced Degradation of α -Cr 2 O 3 Passive Film

The ion exchange and point defect models are two prominent models describing the role of anions, such as chlorides, in the degradation of passive oxide films. Here the thermodynamic feasibility of critical steps of Cl-induced degradation of a hydroxylated α-Cr 2 O 3 (0001) surface, as proposed by these two models, are studied. Both models begin with Cl substitution of surface OH and H 2 O, which becomes less favorable with increasing Cl coverage. The initial stages of Cl-induced breakdown of the α-Cr 2 O 3 depend on Cl coverage and the presence of O vacancy near the surface as follows: (1) neither Cl insertion (supporting the ion exchange model) nor Cr vacancy formation (supporting the point defect model) is feasible at low Cl coverages except in the presence of O vacancies near the surface, where Cl insertion is thermodynamically feasible even at low coverages, (2) in the absence of O vacancies, Cr vacancy formation becomes feasible from 10/12 ML onwards whereas Cl insertion by exchange with subsurface OH only becomes feasible at full coverage. This implies that at higher coverages Cl-induced degradation first initiates through a vacancy formation mechanism, but both insertion and vacancy formation would be feasible at full coverage.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Elucidation of Structure and Physical Properties of Pyrolytic Sugar Oligomers Derived from Cellulose Depolymerization/Dehydration Reactions: A Density Functional Theory Study

Fast pyrolysis of lignocellulosic materials is a promising research area to produce renewable fuels and chemicals. Dehydration is known to be among the most important reaction families during cellulose pyrolysis; water is the most important product. Together with water, dehydration reactions also form a range of poorly known oligomer species of varying molecular sizes, often collected as part of bio-oil water-soluble (WS) fraction. In this work, we used electronic structure calculations to evaluate the relative thermodynamic stabilities of several oligomer species resulting from up to three consecutive dehydration events from cellulose depolymerization intermediates. A library of the thermodynamically favored candidate molecular structures was compiled. Results revealed that most of the water molecules are eliminated from the non-reducing end, forming thermodynamically more stable conjugated compounds. This is consistent with results reported by other researchers in literature where dehydration reactions occur preferably at the non-reducing ends of oligomers. The physical-chemical properties of the proposed structures were estimated using quantitative structure-property relationships (QSPRs) and quantitative property-property relationships (QPPRs). The anhydro-sugars derived from cellulose are often blamed for coke formation during bio-oil hydrotreatment. Understanding their chemical structure could help to develop rational strategies to mitigate coke formation. Furthermore, the thermo-physical properties reported (boiling point, melting point, Gibb’s free energy of formation, enthalpy of formation, and solubility parameters among others) are also fundamental to conducting first principle engineering calculations to design and analyze new pyrolysis reactors and bio-oil up-grading units.

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

A density-functional theory study of the Al/AlOx/Al tunnel junction

The aluminum oxide tunnel junction is a key component of the majority of superconducting quantum devices. For high-quality, reproducible, and scalably manufacturable qubits, the ability to fabricate Josephson junctions (JJs) with a targeted critical current and high uniformity is essential. In this study, we use first-principles modeling to assess fundamental aspects of the atomic structure of both amorphous and crystalline aluminum oxide tunnel junctions and relate the structure to predicted performance metrics. We use modified ab initio molecular dynamics to develop realistic models of the tunnel junction, from which interface roughness and local thickness fluctuations are analyzed in an unbiased manner by training a neural network to identify the boundary between metal and oxide. We show that the effective thickness of the insulating part of the junction can be different from the apparent physical thickness. We calculate the rate of Cooper pair tunneling for the atomically resolved electrostatic potential using direct numerical solution in 3D, which shows a channeling effect that impacts the junction critical current. The predicted critical current is a useful JJ design parameter that can be accessed from the ab initio calculations without fitting parameters. To assess the limits of uniformity and fabrication choices (e.g., oxidation vs epitaxy), we compare the amorphous junctions to crystalline models, which show order of magnitude more efficient tunneling compared to the amorphous case, underlining the connection between atomistic structure and Cooper pair tunneling efficiency. Further, this work provides a foundation for ab initio materials design and evaluation to help accelerate future development of improved tunnel junctions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1983

v-representability and density functional theory

It is shown that if n(r) is the discrete density on a lattice (enclosed in a finite box) associated with a nondegenerate ground state in an external potential v(r) (i.e., is 'v-representable'), then the density n(r) + mu(r), with m(r) arbitrary (apart from trivial constraints) and mu small enough, is also associated with a nondegenerate ground state in an external potential v'(r) near v(r); i.e., n(r) + m(r) is also v-representable. Implications for the Hohenberg-Kohn variational principle and the Kohn-Sham equations are discussed.

Kohn, W.↗
DOE OSTI2024

Putting error bars on density functional theory dataset

This dataset contains submission files and raw output files from high-throughput DFT simulations to analyze the systemic errors in lattice constant, bulk moduli and formation energy predictions for a range of binary and ternary oxides using four exchange correlation functionals (LDA, PBE, PBEsol and vdW-DF-C09). This data was then used as the basis for employing materials informatics methods to predict the expected errors in the lattice constants of the studied compounds. Predicted errors were also used to better the DFT-predicted lattice parameters. Our results emphasize the link between the computed errors and the electron density and hybridization errors of a functional. In essence, these results provide “error bars” for choosing a functional for the creation of high-accuracy, high-throughput datasets as well as avenues for the development of XC functionals with enhanced performance, thereby enabling the accelerated discovery and design of new materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Exact restoration of Galilei invariance in density functional calculations with quantum Monte Carlo

We investigate some aspects of the connection between mean-field and ab-initio calculations for any-nucleon systems employing the Gogny effective interaction. In particular, we present variational Monte Carlo calculations of 4 He and 16 O, for which we extended the formalism of the linear method to the case of a density-dependent Hamiltonian. Monte Carlo calculations using an uncorrelated basis reduce to Hartree-Fock ones, but Galilean invariance can be preserved, contrary to the usual treatment. An analysis of the results obtained with the Gogny force when used in conjunction with correlated wave functions has also been performed. Here, the calculations show the collapse of the deuteron, leading to a set of unphysical predictions for the structure of heavier nuclei. This behavior remarks the different role of the Gogny effective interaction as compared to bare interactions typically used in fully microscopic calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2018

Random Variables with Moment-Matching Staircase Density Functions

This paper proposes a family of random variables for uncertainty modeling. The variables of interest have a bounded support set, and prescribed values for the first four moments. We present the feasibility conditions for the existence of any of such variables, and propose a class of variables that conforms to such constraints. This class is called staircase because the density of its members is a piecewise constant function. Convex optimization is used to calculate their distributions according to several optimality criteria, including maximal entropy and maximal log-likelihood. The flexibility and efficiency of staircases enable modeling phenomena having a possibly skewed and/or multimodal response at a low computational cost. Furthermore, we provide a means to account for the uncertainty in the distribution caused by estimating staircases from data. These ideas are illustrated by generating empirical staircase predictor models. We consider the case in which the predictor matches the sample moments exactly (a setting applicable to large datasets), as well as the case in which the predictor accounts for the sampling error in such moments (a setting applicable to sparse datasets). A predictor model for the dynamics of an aeroelastic airfoil subject to flutter instability is used as an example. The resulting predictor not only describes the system's response accurately, but also enables carrying out a risk analysis for safe flight.

Luis G. Crespo↗
NASA NTRS1992

Local thermodynamic mapping for effective liquid density-functional theory

The structural-mapping approximation introduced by Lutsko and Baus (1990) in the generalized effective-liquid approximation is extended to include a local thermodynamic mapping based on a spatially dependent effective density for approximating the solid phase in terms of the uniform liquid. This latter approximation, called the local generalized effective-liquid approximation (LGELA) yields excellent predictions for the free energy of hard-sphere solids and for the conditions of coexistence of a hard-sphere fcc solid with a liquid. Moreover, the predicted free energy remains single valued for calculations with more loosely packed crystalline structures, such as the diamond lattice. The spatial dependence of the weighted density makes the LGELA useful in the study of inhomogeneous solids.

Kyrlidis, Agathagelos↗
DOE OSTI2020

Real-time time-dependent density functional theory in Quantum Espresso and related codes

A set of subroutines that interact with an existing computer program, Quantum Espresso (QE), that takes a set of electronic orbitals from QE and propagates them forward in time in response to a number of perturbations. These perturbations include the motion of ions and a wide spectrum of applied electromagnetic fields (e.g., static magnetic fields, slowly-varying electric fields, x-rays, etc.). From the time evolution of the electronic orbitals comes time-varying values of physical observables like the electronic charge and current densities and atomic forces. These are further post-processed to yield properties of interest like stopping powers, conductivities, and the dynamic structure factor. The electronic orbitals are described in a plane wave basis and the equations of motion are carried out using a unitary or approximately unitary time propagation algorithm. Related work will include subroutines for interfacing this capability, as well as others that already exist within QE, with a fork of QE that is more extensible and modular.SAND2019-12611 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Baczewski, Andrew↗
DOE OSTI2023

Distributed memory, GPU accelerated Fock construction for hybrid, Gaussian basis density functional theory

With the growing reliance of modern supercomputers on accelerator-based architecture such a graphics processing units (GPUs), the development and optimization of electronic structure methods to exploit these massively parallel resources has become a recent priority. While significant strides have been made in the development GPU accelerated, distributed memory algorithms for many modern electronic structure methods, the primary focus of GPU development for Gaussian basis atomic orbital methods has been for shared memory systems with only a handful of examples pursing massive parallelism. Here in this work, we present a set of distributed memory algorithms for the evaluation of the Coulomb and exact exchange matrices for hybrid Kohn–Sham DFT with Gaussian basis sets via direct density-fitted (DF-J-Engine) and seminumerical (sn-K) methods, respectively. The absolute performance and strong scalability of the developed methods are demonstrated on systems ranging from a few hundred to over one thousand atoms using up to 128 NVIDIA A100 GPUs on the Perlmutter supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Novel Ultrabright and Air-Stable Photocathodes Discovered from Machine Learning and Density Functional Theory Driven Screening

The high brightness, low emittance electron beams achieved in modern X-ray free-electron lasers (XFELs) have enabled powerful X-ray imaging tools, allowing molecular systems to be imaged at picosecond time scales and sub-nanometer length scales. One of the most promising directions for increasing the brightness of XFELs is through the development of novel photocathode materials. Whereas past efforts aimed at discovering photocathode materials have typically employed trial-and-error-based iterative approaches, this work represents the first data-driven screening for high brightness photocathode materials. Through screening over 74 000 semiconducting materials, a vast photocathode dataset is generated, resulting in statistically meaningful insights into the nature of high brightness photocathode materials. This screening results in a diverse list of photocathode materials that exhibit intrinsic emittances that are up to 4x lower than currently used photocathodes. In a second effort, multiobjective screening is employed to identify the family of M 2 O (M = Na, K, Rb) that exhibits photoemission properties that are comparable to the current state-of-the-art photocathode materials, but with superior air stability. This family represents perhaps the first intrinsically bright, visible light photocathode materials that are resistant to reactions with oxygen, allowing for their transport and storage in dry air environments.

36 MATERIALS SCIENCE↗