Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional calculation”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effects of substitutional Ga on the phonons of δ -phase Pu based on density functional theory calculations

The delta phase of Pu is stabilized by Ga doping, but the mechanism of this stabilization remains an open question. Density functional theory calculations focused on how Ga doping affects the phonons sheds some light on the phonons' contribution to the stabilization. Here, the calculated phonon modes of Ga-doped delta phase Pu fall into two distinct types: localized, high frequency Ga-dominated phonon modes, and Pu-dominated modes at lower frequencies. Increasing the Ga concentration has an effect on the Pu-dominated phonon modes opposite to that of compression: higher-frequency modes soften, and lower-frequency modes stiffen. The latter provides an indication that the stabilization mechanism is not due to a thermodynamic contribution from the phonons. Furthermore, the stiffened phonon modes include candidate modes that describe possible pathways into low-temperature phases, suggesting that doping with Ga could impede such pathways.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Enhanced Activity of Heterogeneous Pd(II) Catalysts on Acid-Functionalized Metal–Organic Frameworks

Single-site heterogeneous catalysts (SSHCs) play important roles in fundamental science and technology, owing to the molecular level control of structure–support interactions that is possible in these systems. Recently, SSHCs supported by acidic oxides have attracted particular interest because catalytically active metal centers can be formed at the surface sites. Here, we incorporated a palladium SSHC in phosphated and sulfated metal–organic frameworks (MOFs), hafnium-based MOF-808 (Hf-MOF-808-PO 4 and Hf-MOF-808-SO 4 ). The structural and electronic properties of the Pd(II) sites coordinated to the acidic sites in these MOFs were investigated through X-ray photoelectron spectroscopy, vibrational spectroscopy, X-ray crystallographic techniques, catalytic studies, and quantum mechanical electronic structure calculations employing density functional theory. We demonstrated that the presence of node-bound acidic functional groups stabilizes the Pd(II) site in these MOFs, resulting in enhanced catalytic activities (compared to in the nonacid functionalized Hf-MOF-808) in the oxidative Heck reaction where Pd(II) is the active species. The density functional calculations support the interpretation that the acid functionalization of the MOF node can stabilize the Pd(0) intermediate state during the catalytic reactions, thereby suppressing Pd(0) aggregation leading to catalyst deactivation. Furthermore, these findings offer insights and methodology for the catalytic investigation of SSHCs in MOFs.

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

Density functional theory calculation of helium segregation and decohesion effect in W110/Ni111 interphase boundary

Density functional theory (DFT) is employed to study the preferential distribution and decohesion effect of He in a W-NiFe composite consisting of W particles embedded in an Ni-based solid solution matrix. A slab containing {110}<100>W//{111}<110>Ni interface is employed as a surrogate model for the W-NiFe system. Firstly, the fracture energy of the W/Ni interphase boundary (IB) (4.37 J/m 2 ) is higher than the cleavage energy of Ni{111} (3.82 J/m^2) and lower than the cleavage energy of W{110} (6.60 J/m 2 ). The comparison shows that the cohesion of the IB is stronger than the Ni{111} planes that are away from the IB. However, the cohesion between the Ni{111} planes adjacent to the IB is found to be the weakest in this system, with a cleavage energy of 3.11 J/m 2 . Subsequently, formation energy of He is calculated in the Ni slab, W slab, and various interstitial sites in the IB. The DFT calculations show that He is significantly more stable in Ni than in W, by about 1.75 eV. Interestingly, He does not prefer to segregate at the IB as compared to bulk Ni. Nevertheless, it prefers to segregate to the region between the Ni{111} planes adjacent to the IB, and decreases the cohesion of the already weakest region. Based on an estimated amount of He gas production in 5 years under first wall neutron irradiation (neutron flux of 1.04x10 15 n/cm 2 /s), the He will decrease the cleavage energy of the weakest region by 21.2%, 15.4%, and 12.2% at 800, 1000, and 1200 C, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Selective Host-Guest Interaction between Metal Ions and Metal-Organic Frameworks Using Dynamic Nuclear Polarization Enhanced Solid-State NMR Spectroscopy

The host–guest interaction between metal ions (Pt²⁺ and Cu²⁺) and a zirconium metal–organic framework (UiO-66-NH₂) was explored using dynamic nuclear polarization-enhanced ¹⁵N{¹H} CPMAS NMR spectroscopy supported by X-ray absorption spectroscopy and density functional calculations. The combined experimental results conclude that each Pt²⁺ coordinates with two NH₂ groups from the MOF and two Cl⁻ from the metal precursor, whereas Cu²⁺ do not form chemical bonds with the NH₂ groups of the MOF framework. Density functional calculations reveal that Pt²⁺ prefers a square-planar structure with the four ligands and resides in the octahedral cage of the MOF in either cis or trans configurations.

Guo, Zhiyong↗
DOE OSTI2016

Difficulty of carrier generation in orthorhombic PbO

Polycrystalline β-PbO films were grown by pulsed laser deposition in atmospheres ranging from oxygen-poor (the oxygen pressure of 0.01 Pa) to oxygen-rich (13 Pa) conditions, and the oxygen chemical potential was further enhanced by ozone annealing to examine hole doping. It was found that each of the as-grown β-PbO films showed poor electrical conductivity, σ < 1.4 × 10{sup −7} S cm{sup −1}, regardless of the oxygen pressure. The density functional calculations revealed that native defects including Pb and O vacancies have deep transition levels and extremely high formation enthalpies, which indicates difficulty of carrier generation in β-PbO and explains the experimentally observed poor electrical conductivity. The analysis of the electronic structures showed that the interaction between Pb 6s and O 2p orbitals is weak due to the deep energy level of Pb 6s and does not raise the valence band maximum (VBM) level unlike that observed in SnO, which is also supported by ultraviolet photoemission spectroscopy measurements. The deep acceptor transition levels of the native defects are attributed to the deep VBM of β-PbO. On the other hand, annealing β-PbO films in reactive oxygen-containing atmospheres (i.e., O{sub 3}) led to a significantly enhanced electrical conductivity (i.e., σ > 7.1 × 10{sup 2} S cm{sup −1}) but it is the result of the formation of an n-type PbO{sub 2} phase because oxygen chemical potential exceeded the phase boundary limit. The striking difference in carrier generation between PbO and SnO is discussed based on the electronic structures calculated by density functional theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Spectral scheme for atomic structure calculations in density functional theory

In this study, we present a spectral scheme for atomic structure calculations in pseudopotential Kohn-Sham density functional theory. In particular, after applying an exponential transformation of the radial coordinates, we employ global polynomial interpolation on a Chebyshev grid, with derivative operators approximated using the Chebyshev differentiation matrix, and integrations using Clenshaw-Curtis quadrature. We demonstrate the accuracy and efficiency of the scheme through spin-polarized and unpolarized calculations for representative atoms, while considering local, semilocal, and hybrid exchange-correlation functionals. In particular, we find that $\mathcal{O}$(200) grid points are sufficient to achieve an accuracy of 1 microhartree in the eigenvalues for optimized norm conserving Vanderbilt pseudopotentials spanning the periodic table from atomic number Ζ = 1 to 83.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

DGDFT: A massively parallel method for large scale density functional theory calculations

We describe a massively parallel implementation of the recently developed discontinuous Galerkin density functional theory (DGDFT) method, for efficient large-scale Kohn-Sham DFT based electronic structure calculations. The DGDFT method uses adaptive local basis (ALB) functions generated on-the-fly during the self-consistent field iteration to represent the solution to the Kohn-Sham equations. The use of the ALB set provides a systematic way to improve the accuracy of the approximation. By using the pole expansion and selected inversion technique to compute electron density, energy, and atomic forces, we can make the computational complexity of DGDFT scale at most quadratically with respect to the number of electrons for both insulating and metallic systems. We show that for the two-dimensional (2D) phosphorene systems studied here, using 37 basis functions per atom allows us to reach an accuracy level of 1.3 × 10{sup −4} Hartree/atom in terms of the error of energy and 6.2 × 10{sup −4} Hartree/bohr in terms of the error of atomic force, respectively. DGDFT can achieve 80% parallel efficiency on 128,000 high performance computing cores when it is used to study the electronic structure of 2D phosphorene systems with 3500-14 000 atoms. This high parallel efficiency results from a two-level parallelization scheme that we will describe in detail.

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

Finding the Pareto front for high-entropy-alloy catalysts

Finding catalysts that have both high activity and high stability presents a long-standing challenge. Since optimizing activity and stability are conflicting objectives, the best one can do is find the Pareto front that yields optimal tradeoffs between these features. On the Pareto front, there is a trade-off where a portion of catalytic activity must be sacrificed to gain further stability and vice versa . Here, we provide a method to optimize the front by designing a multi-objective genetic algorithm that combines machine learning, graph neural network calculations, and density functional calculations. The application considered is the oxygen evolution reaction catalyzed by high-entropy alloys. We find that the Pareto front generally contains alloys with diverse elements, but that enhancing stability inevitably inflicts a toll on activity. We compare the general conclusions of our work to a survey of 545 experiments.

Zhang, Chengyi [Univ. of Auckland (New Zealand)]↗
DOE OSTI2021

Optical response and band structure of LiCoO 2 including electron-hole interaction effects

The optical response functions and band structures of LiCoO 2 are studied at different levels of approximation, from density functional theory (DFT) in the generalized gradient approximation (GGA) to quasiparticle self-consistent QSGW (with G for Green's function and W for screened Coulomb interaction) without and with ladder diagrams ($QSG\hat{W}$) and the Bethe Salpeter Equation (BSE) approach. The QSGW method is found to strongly overestimate the band gap and electron-hole or excitonic effects are found to be important. They lower the quasiparticle gap by only about 11~\% but the lowest energy peaks in absorption are found to be excitonic in nature. The contributions from different band to band transitions and the relation of excitons to band-to-band transitions are analyzed. Finally, the excitons are found to be strongly localized. Finally, a comparison to experimental data is presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Electronically induced defect creation at semiconductor/oxide interface revealed by time-dependent density functional theory

Carrier induced defect creation at the semiconductor-oxide interface has been known as the origin of electronic device degradation for a long time, but how exactly the interface lattice can be damaged by carriers (especially low-energy ones) remains unclear. Here we carry out real-time time-dependent density functional theory simulations on concrete Si/SiO 2 interfaces to study the interaction between excited electrons and interface bonds. We show that the normal interface Si-H bonds are generally resistant to electrons due to the delocalized nature and high energy level of the Si-H antibonding states, and due to the high-energy barrier to break the Si-H bond. However, if an additional hydrogen atom exists by attaching to a nearby oxygen atom (forming a “Si-H···H-O” complex), the Si-H bond will be greatly weakened, including the reduction of energy barrier for bond breaking, and the lowering of the antibonding state energy level which favors electron injection. Together with the multiple vibrational excitation process, the corresponding Si-H bond can be broken much more easily. Thus we propose that the Si-H···H-O complex will be the center for defect creation and device degradation. Finally, we also explain why such a center might be relatively easy to form during the hydrogen annealing process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Highly accurate and constrained density functional obtained with differentiable programming

Using an end-to-end differentiable implementation of the Kohn-Sham self-consistent field equations, we obtain a highly accurate neural network–based exchange and correlation (XC) functional of the electronic density. The functional is optimized using information on both energy and density while exact constraints are enforced through an appropriate neural network architecture. Here we evaluate our model against different families of XC approximations and show that at the meta-GGA level our functional exhibits unprecedented accuracy for both energy and density predictions. For nonempirical functionals, there is a strong linear correlation between energy and density errors. We use this correlation to define an XC functional quality metric that includes both energy and density errors, leading to an improved way to rank different approximations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Pure quantum extension of the semiclassical Boltzmann-Uehling-Uhlenbeck equation

The Boltzmann equation is the traditional framework in which one extends the time-dependent mean field classical description of a many-body system to include the effect of particle-particle collisions in an approximate manner. A semiclassical extension of this approach to quantum many-body systems was suggested by Uehling and Uhlenbeck in 1933 for both Fermi and Bose statistics, and many further developments of this approach are known as the Boltzmann-Uehling-Uhlenbeck (BUU) equations. In this paper I introduce a pure quantum version of the BUU type of equations, which is mathematically equivalent to a generalized time-dependent density functional theory extended to superfluid systems. As expected, during nonequilibrium processes the quantum Boltzmann one-body entropy increases during evolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Fast Real-Time Time-Dependent Density Functional Theory Calculations with the Parallel Transport Gauge

Real-time time-dependent density functional theory (RT-TDDFT) is known to be hindered by the very small time step (attosecond or smaller) needed in the numerical simulation, because of the fast oscillation of electron wave functions, which significantly limits its range of applicability for the study of ultrafast dynamics. In this paper, we demonstrate that such oscillation can be considerably reduced by optimizing the gauge choice using the parallel transport formalism. RT-TDDFT calculations can thus be significantly accelerated using a combination of the parallel transport gauge and implicit integrators, and the resulting scheme can be used to accelerate any electronic structure software that uses a Schrödinger representation. Here, using absorption spectrum, ultrashort laser pulse, and Ehrenfest dynamics calculations for example, we show that the new method can utilize a time step that is on the order of 10–100 attoseconds using a planewave basis set. Thanks to the significant increase of the size of the time step, we also demonstrate that the new method is more than 10 times faster, in terms of the wall clock time, when compared to the standard explicit fourth-order Runge–Kutta time integrator for silicon systems ranging from 32 to 1024 atoms.

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

Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements

Here, we present a computationally efficient approach to perform systematically convergent real-space all-electron Kohn-Sham density functional theory calculations for solids using an enriched finite element (FE) basis. The enriched FE basis is constructed by augmenting the classical FE basis with atom-centered numerical basis functions, comprising of atomic solutions to the Kohn-Sham problem. Notably, to improve the conditioning, we orthogonalize the enrichment functions with respect to the classical FE basis, without sacrificing the locality of the resultant basis. In addition to improved conditioning, this orthogonalization procedure also renders the overlap matrix block diagonal, greatly simplifying its inversion. Subsequently, we use a Chebyshev polynomial based filtering technique to efficiently compute the occupied eigenspace in each self-consistent field iteration. We demonstrate the accuracy and efficiency of the proposed approach on periodic unit cells and supercells. The benchmark studies show a staggering 130× speedup of the orthogonalized enriched FE basis over the classical FE basis. We also present a comparison of the orthogonalized enriched FE basis with the linearized augmented plane-wave + local orbitals basis, both in terms of accuracy and efficiency. Notably, we demonstrate that the orthogonalized enriched FE basis can handle large system sizes of ~10 000 electrons. Finally, we observe good parallel scalability of our implementation with 92% efficiency at 22× speedup for a system with 620 electrons

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2002

Time Dependent Density Functional Theory Calculations of Large Compact PAH Cations: Implications for the Diffuse Interstellar Bands

We investigate the electronic absorption spectra of several maximally pericondensed polycyclic aromatic hydrocarbon radical cations with time dependent density functional theory calculations. We find interesting trends in the vertical excitation energies and oscillator strengths for this series containing pyrene through circumcoronene, the largest species containing more than 50 carbon atoms. We discuss the implications of these new results for the size and structure distribution of the diffuse interstellar band carriers.

Weisman, Jennifer L.↗
DOE OSTI2014

Electronic and structural properties of ultrathin tungsten nanowires and nanotubes by density functional theory calculation

The simulated annealing basin-hopping method incorporating the penalty function was used to predict the lowest-energy structures for ultrathin tungsten nanowires and nanotubes of different sizes. These predicted structures indicate that tungsten one-dimensional structures at this small scale do not possess B.C.C. configuration as in bulk tungsten material. In order to analyze the relationship between multi-shell geometries and electronic transfer, the electronic and structural properties of tungsten wires and tubes including partial density of state and band structures which were determined and analyzed by quantum chemistry calculations. In addition, in order to understand the application feasibility of these nanowires and tubes on nano-devices such as field emitters or chemical catalysts, the electronic stability of these ultrathin tungsten nanowires was also investigated by density functional theory calculations.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2015

Hybrid functional band gap calculation of SnO{sub 6} containing perovskites and their derived structures

We have studied the properties of SnO{sub 6} octahedra-containing perovskites and their derived structures using ab initio calculations with different density functionals. In order to predict the correct band gap of the materials, we have used B3LYP hybrid density functional, and the results of B3LYP were compared with those obtained using the local density approximation and generalized gradient approximation data. The calculations have been conducted for the orthorhombic ground state of the SnO{sub 6} containing perovskites. We also have expended the hybrid density functional calculation to the ASnO{sub 3}/A'SnO{sub 3} system with different cation orderings. We propose an empirical relationship between the tolerance factor and the band gap of SnO{sub 6} containing oxide materials based on first principles calculation. - Graphical abstract: (a) Structure of ASnO{sub 3} for orthorhombic ground state. The green ball is A (Ba, Sr, Ca) cation and the small (red) ball on edge is oxygen. SnO{sub 6} octahedrons are plotted as polyhedron. (b) Band gap of ASnO{sub 3} as a function of the tolerance factor for different density functionals. The experimental values of the band gap are marked as green pentagons. (c) ASnO{sub 3}/A'SnO{sub 3} superlattices with two types cation arrangement: [001] layered structure and [111] rocksalt structure, respectively. (d) B3LYP hybrid functional band gaps of ASnO{sub 3}, [001] ordered superlattices, and [111] ordered superlattices of ASnO{sub 3}/A'SnO{sub 3} as a function of the effective tolerance factor. Note the empirical linear relationship between the band gap and effective tolerance factor. - Highlights: • We report the hybrid functional band gap calculation of ASnO{sub 3} and ASnO{sub 3}/A'SnO{sub 3}. • The band gap of ASnO{sub 3} using B3LYP functional reproduces the experimental value. • We propose the linear relationship between the tolerance factor and the band gap.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗