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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Native defects in Tl{sub 6}SI{sub 4}: Density functional calculations

Tl{sub 6}SI{sub 4} is a promising room-temperature semiconductor radiation detection material. Here, we report density functional calculations of native defects and dielectric properties of Tl{sub 6}SI{sub 4}. Formation energies and defect levels of native point defects and defect complexes are calculated. Donor-acceptor defect complexes are shown to be abundant in Tl{sub 6}SI{sub 4}. High resistivity can be obtained by Fermi level pinning by native donor and acceptor defects. Deep donors that are detrimental to electron transport are identified and methods to mitigate such problem are discussed. Furthermore, we show that mixed ionic-covalent character of Tl{sub 6}SI{sub 4} gives rise to enhanced Born effective charges and large static dielectric constant, which provides effective screening of charged defects and impurities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Charge compensation and electrostatic transferability in three entropy-stabilized oxides: Results from density functional theory calculations

Density functional theory calculations were carried out for three entropic rocksalt oxides, (Mg{sub 0.1}Co{sub 0.1}Ni{sub 0.1}Cu{sub 0.1}Zn{sub 0.1})O{sub 0.5}, termed J14, and J14 + Li and J14 + Sc, to understand the role of charge neutrality and electronic states on their properties, and to probe whether simple expressions may exist that predict stability. The calculations predict that the average lattice constants of the ternary structures provide good approximations to that of the random structures. For J14, Bader charges are transferable between the binary, ternary, and random structures. For J14 + Sc and J14 + Li, average Bader charges in the entropic structures can be estimated from the ternary compositions. Addition of Sc to J14 reduces the majority of Cu, which show large displacements from ideal lattice sites, along with reduction of a few Co and Ni cations. Addition of Li to J14 reduces the lattice constant, consistent with experiment, and oxidizes some of Co as well as some of Ni and Cu. The Bader charges and spin-resolved density of states (DOS) for Co{sup +3} in J14 + Li are very different from Co{sup +2}, while for Cu and Ni the Bader charges form continuous distributions and the two DOS are similar for the two oxidation states. Experimental detection of different oxidation states may therefore be challenging for Cu and Ni compared to Co. Based on these results, empirical stability parameters for these entropic oxides may be more complicated than those for non-oxide entropic solids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Electronic Structure, Phonon Dynamical Properties, and CO 2 Capture Capability of Na 2 - x M x Zr O 3 ( M = Li ,K): Density-Functional Calculations and Experimental Validations

The electronic structural and phonon properties of Na 2-α M α ZrO 3 (M ¼ Li,K, α = ¼ 0.0,0.5,1.0,1.5,2.0) are investigated by first-principles density-functional theory and phonon dynamics. The thermodynamic properties of CO 2 absorption and desorption in these materials are also analyzed. With increasing doping level α, the binding energies of Na 2-α Li α ZrO 3 are increased while the binding energies of Na 2-α K α ZrO 3 are decreased to destabilize the structures. The calculated band structures and density of states also show that, at the same doping level, the doping sites play a significant role in the electronic properties. The phonon dispersion results show that few soft modes are found in several doped configurations, which indicates that these structures are less stable than other configurations with different doping levels. From the calculated relationships among the chemical-potential change, the CO 2 pressure, and the temperature of the CO 2 capture reactions by Na 2-α M α ZrO 3 , and from thermogravimetric-analysis experimental measurements, the Li- and K-doped mixtures Na 2-α M α ZrO 3 have lower turnover temperatures (T t ) and higher CO 2 capture capacities, compared to pure Na 2 ZrO 3 . The Li-doped systems have a larger T t decrease than the K-doped systems. When increasing the Li-doping level α, the T t of the corresponding mixture Na 2-α Li α ZrO 3 decreases further to a low-temperature range. However, in the case of K-doped systems Na 2-α K α ZrO 3 , although doping K into Na 2 ZrO 3 initially shifts its T t to lower temperatures, further increases of the K-doping level α causes T t to increase. Therefore, doping Li into Na 2 ZrO 3 has a larger influence on its CO 2 capture performance than the K-doped Na 2 ZrO 3 . Compared with pure solidsM 2 ZrO 3 , after doping with other elements, these doped systems’ CO 2 capture performances are improved.

Duan, Yuhua↗
DOE OSTI2015

Large-scale real-space density-functional calculations: Moiré-induced electron localization in graphene

We show that our real-space finite-difference scheme allows us to perform density-functional calculations for nanometer-scale targets containing more than 100 000 atoms. This real-space scheme is applied to twisted bilayer graphene, clarifying that Moiré pattern induced in the slightly twisted bilayer graphene drastically modifies the atomic and electronic structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Adaptively Compressed Exchange Operator for Large-Scale Hybrid Density Functional Calculations with Applications to the Adsorption of Water on Silicene

Density functional theory (DFT) calculations using hybrid exchange-correlation functionals have been shown to provide an accurate description of the electronic structures of nanosystems. However, such calculations are often limited to small system sizes due to the high computational cost associated with the construction and application of the Hartree-Fock (HF) exchange operator. Here, we demonstrate that the recently developed adaptively compressed exchange (ACE) operator formulation [J. Chem. Theory Comput. 2016, 12, 2242-2249] can enable hybrid functional DFT calculations for nanosystems with thousands of atoms. The cost of constructing the ACE operator is the same as that of applying the exchange operator to the occupied orbitals once, while the cost of applying the Hamiltonian operator with a hybrid functional (after construction of the ACE operator) is only marginally higher than that associated with applying a Hamiltonian constructed from local and semilocal exchange-correlation functionals. Therefore, this new development significantly lowers the computational barrier for using hybrid functionals in large-scale DFT calculations. We demonstrate that a parallel planewave implementation of this method can be used to compute the ground-state electronic structure of a 1000-atom bulk silicon system in less than 30 wall clock minutes and that this method scales beyond 8000 computational cores for a bulk silicon system containing about 4000 atoms. The efficiency of the present methodology in treating large systems enables us to investigate adsorption properties of water molecules on Ag-supported two-dimensional silicene. Our computational results show that water monomer, dimer, and trimer configurations exhibit distinct adsorption behaviors on silicene. In particular, the presence of additional water molecules in the dimer and trimer configurations induces a transition from physisorption to chemisorption, followed by dissociation on Ag-supported silicene. This is caused by the enhanced effect of hydrogen bonds on charge transfer and proton transfer processes. Lastly, such a hydrogen bond autocatalytic effect is expected to have broad applications for silicene as an efficient surface catalyst for oxygen reduction reactions and water dissociation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Accurate Hellmann–Feynman forces from density functional calculations with augmented Gaussian basis sets

The Hellmann–Feynman (HF) theorem provides a way to compute forces directly from the electron density, enabling efficient force calculations for large systems through machine learning (ML) models for the electron density. The main issue holding back the general acceptance of the HF approach for atom-centered basis sets is the well-known Pulay force which, if naively discarded, typically constitutes an error upward of 10 eV/Å in forces. In this work, we demonstrate that if a suitably augmented Gaussian basis set is used for density functional calculations, the Pulay force can be suppressed, and HF forces can be computed as accurately as analytical forces with state-of-the-art basis sets, allowing geometry optimization and molecular dynamics to be reliably performed with HF forces. Our results pave a clear path forward for the accurate and efficient simulation of large systems using ML densities and the HF theorem.

Pathak, Shivesh (ORCID:0000000172798195)↗
DOE OSTI2018

Accurate double excitations from ensemble density functional calculations

The recent use of a new ensemble in density functional theory (DFT) to produce direct corrections to the Kohn-Sham transitions yields the elusive double excitations that are missed by time-dependent DFT (TDDFT) with the standard adiabatic approximation. But accuracies are lower than for single excitations, and formal arguments about TDDFT suggest that a correction kernel is needed. In principle, ensemble DFT with direct corrections at the exchange level must yield accurate doubles in the weakly correlated limit. We illustrate with exact calculations and analytic results on the Hubbard dimer. Lastly, we also explain the error in formal arguments in TDDFT.

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

Density Functional Calculations of Native Defects in CH 3 NH 3 PbI 3 : Effects of Spin–Orbit Coupling and Self-Interaction Error

We know that native point defects play an important role in carrier transport properties of CH3NH3PbI3. However, the nature of many important defects remains controversial due partly to the conflicting results reported by recent density functional theory (DFT) calculations. In this Letter, we show that self-interaction error and the neglect of spin–orbit coupling (SOC) in many previous DFT calculations resulted in incorrect positions of valence and conduction band edges, although their difference, which is the band gap, is in good agreement with the experimental value. Moreover, this problem has led to incorrect predictions of defect-level positions. Hybrid density functional calculations, which partially correct the self-interaction error and include the SOC, show that, among native point defects (including vacancies, interstitials, and antisites), only the iodine vacancy and its complexes induce deep electron and hole trapping levels inside of the band gap, acting as nonradiative recombination centers.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular NMR shieldings, J -couplings, and magnetizabilities from numeric atom-centered orbital based density-functional calculations

This paper reports and benchmarks a new implementation of nuclear magnetic resonance shieldings, magnetizabilities, and J-couplings for molecules within semilocal density functional theory, based on numeric atom-centered orbital (NAO) basis sets. NAO basis sets are attractive for the calculation of these nuclear magnetic resonance (NMR) parameters because NAOs provide accurate atomic orbital representations especially near the nucleus, enabling high-quality results at modest computational cost. Moreover, NAOs are readily adaptable for linear scaling methods, enabling efficient calculations of large systems. Here, the paper has five main parts: (1) It reviews the formalism of density functional calculations of NMR parameters in one comprehensive text to make the mathematical background available in a self-contained way. (2) The paper quantifies the attainable precision of NAO basis sets for shieldings in comparison to specialized Gaussian basis sets, showing similar performance for similar basis set size. (3) The paper quantifies the precision of calculated magnetizabilities, where the NAO basis sets appear to outperform several established Gaussian basis sets of similar size. (4) The paper quantifies the precision of computed J-couplings, for which a group of customized NAO basis sets achieves precision of ~Hz for smaller basis set sizes than some established Gaussian basis sets. (5) The paper demonstrates that the implementation is applicable to systems beyond 1000 atoms in size.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

GPU acceleration of local and semilocal density functional calculations in the SPARC electronic structure code

We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing Kohn–Sham density functional theory calculations within the local density and generalized gradient approximations. In particular, we develop a modular math-kernel based implementation for NVIDIA architectures wherein the computationally expensive operations are carried out on the GPUs, with the remainder of the workload retained on the central processing units (CPUs). Here, using representative bulk and slab examples, we show that relative to CPU-only execution, GPUs enable speedups of up to 6× and 60× in node and core hours, respectively, bringing time to solution down to less than 30 s for a metallic system with over 14 000 electrons and enabling significant reductions in computational resources required for a given wall time.

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

Interpolative Separable Density Fitting Decomposition for Accelerating Hybrid Density Functional Calculations with Applications to Defects in Silicon

We present a new efficient way to perform hybrid density functional theory (DFT)-based electronic structure calculations. The new method uses an interpolative separable density fitting (ISDF) procedure to construct a set of numerical auxiliary basis vectors and a compact approximation of the matrix consisting of products of occupied orbitals represented in a large basis set such as the planewave basis. Such an approximation allows us to reduce the number of Poisson solves from O ($N$ $^{2}_{ϵ}$) to O ( N ϵ ) when we apply the exchange operator to occupied orbitals in an iterative method for solving the Kohn–Sham equations, where Ne is the number of electrons in the system to be studied. We show that the ISDF procedure can be carried out in O ($N$ $^{3}_{ϵ}$) operations, with a much smaller preconstant compared to methods used in existing approaches. When combined with the recently developed adaptively compressed exchange (ACE) operator formalism, which reduces the number of times the exchange operator needs to be updated, the resulting ACE-ISDF method significantly reduces the computational cost associated with the exchange operator by nearly 2 orders of magnitude compared to existing approaches for a large silicon system with 1000 atoms. We demonstrate that the ACE-ISDF method can produce accurate energies and forces for insulating and metallic systems and that it is possible to obtain converged hybrid functional calculation results for a 1000-atom bulk silicon within 10 min on 2000 computational cores. We also show that ACE-ISDF can scale to 8192 computational cores for a 4096-atom bulk silicon system. We use the ACE-ISDF method to geometrically optimize a 1000-atom silicon system with a vacancy defect using the HSE06 functional and computes its electronic structure. Here, we find that that the computed energy gap from the HSE06 functional is much closer to the experimental value compared to that produced by semilocal functionals in the DFT calculations.

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

Instability of Rock-Salt Cubic NbN in Density Functional Calculations

All-nitride semiconductor/superconductor heterojunctions utilizing cubic niobium nitride (NbN) are a promising approach to superconducting quantum circuits for next-generation quantum-information systems. However, there are fundamental open questions on the atomic structure of NbN. In calculations with several levels of density functional theory (DFT), we find that the cubic rocksalt structure NbN (Fm-3m, 225) is energetically very unstable against the ground state hexagonal NbN in a tungsten carbide (P-6m2, 187) type lattice. To better understand the appearance of a cubic phase in numerous experiments, we perform a DFT study on possible NbN structures, determining the energy ordering between different polymorphs from databases and structure prediction. We perform supercell calculations of disordered NbN, finding that the rocksalt structure is dynamically unstable and relaxes to a lower energy monoclinic phase (C2/m, 12), which retains an approximate average cubic symmetry. However, the associated energy gain is not substantial enough for a plausible explanation for cubic NbN. We further investigate the role of external factors such as in-plane strain during epitaxial growth as well as the presence of off-stoichiometry and impurity doping on the energy ordering.

CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SU↗
DOE OSTI2015

Density functional calculation of the structural and electronic properties of germanium quantum dots

We apply first principles density functional computational methods to study the structures, densities of states (DOS), and higher occupied molecular orbital (HOMO) – lowest unoccupied molecular orbital (LUMO) gaps of selected free-standing Ge semiconductor quantum dots up to 1.8nm. Our calculations are performed using numerical atomic orbital approach where linear combination of atomic orbital was applied. The surfaces of the quantum dots was passivized by hydrogen atoms. We find that surface passivation does affect the electronic properties associated with the changes of surface state, electron localization, and the energy gaps of germanium nanocrystals as well as the confinement of electrons inside the quantum dots (QDs). Our study shows that the energy gaps of germanium quantum dots decreases with the increasing dot diameter. The size-dependent variations of the computed HOMO-LUMO gaps in our quantum dots model were found to be consistent with the effects of quantum confinement reported in others theoretical and experimental calculation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Diffusion quantum Monte Carlo and density functional calculations of the structural stability of bilayer arsenene

We have studied the structural stability of monolayer and bilayer arsenene (As) in the buckled (b) and washboard (w) phases with diffusion quantum Monte Carlo (DMC) and density functional theory (DFT) calculations. DMC yields cohesive energies of 2.826(2) eV/atom for monolayer b-As and 2.792(3) eV/atom for w-As. In the case of bilayer As, DMC and DFT predict that AA-stacking is the more stable form of b-As, while AB is the most stable form of w-As. The DMC layer-layer binding energies for b-As-AA and w-As-AB are 30(1) and 53(1) meV/atom, respectively. The interlayer separations were estimated with DMC at 3.521(1) Å for b-As-AA and 3.145(1) Å for w-As-AB. A comparison of DMC and DFT results shows that the van der Waals density functional method yields energetic properties of arsenene close to DMC, while the DFT + D3 method closely reproduced the geometric properties from DMC. The electronic properties of monolayer and bilayer arsenene were explored with various DFT methods. The bandgap values vary significantly with the DFT method, but the results are generally qualitatively consistent. Finally, we expect the present work to be useful for future experiments attempting to prepare multilayer arsenene and for further development of DFT methods for weakly bonded systems.

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↗
DOE OSTI2017

Ab Initio Density Functional Calculations and Infra-Red Study of CO Interaction with Pd Atoms on θ-Al 2 O 3 (010) Surface

The ab initio density functional theoretical studies show that energetics favor CO oxidation on single Pd atoms supported on θ-alumina. The diffuse reflectance infra-red spectroscopy (DRIFTS) results show that carbonates are formed as intermediates when single supported Pd atoms are exposed to a gaseous mixture of CO + O 2 . The rapid agglomeration of Pd atoms under CO oxidation conditions even at 6°C leads to the presence of Pd particles along with single atoms during CO oxidation experiments. Thus, the observed CO oxidation has contributions from both single Pd atoms and Pd particles.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Importance of self-interaction-error removal in density functional calculations on water cluster anions

Accurate description of the excess charge in water cluster anions is challenging for standard semi-local and (global) hybrid density functional approximations (DFAs). Using the recent unitary invariant implementation of the Perdew–Zunger self-interaction correction (SIC) method using Fermi-Löwdin orbitals, we assess the effect of self-interaction error on the vertical detachment energies of water cluster anions with the local spin density approximation (LSDA), Perdew–Burke–Ernzerhof (PBE) generalized gradient approximation, and the strongly constrained and appropriately normed (SCAN) meta-GGA functionals. Our results show that for the relative energies of isomers with respect to reference CCSD(T) values, the uncorrected SCAN functional has the smallest deviation of 21 meV, better than that for the MP2 method. The performance of SIC-SCAN is comparable to that of MP2 and is better than SIC-LSDA and SIC-PBE, but it reverses the ordering of the two lowest isomers for water hexamer anions. Removing self interaction error (SIE) corrects the tendency of LSDA, PBE, and SCAN to over-bind the extra electron. The vertical detachment energies (VDEs) of water cluster anions, obtained from the total energy differences of corresponding anion and neutral clusters, are significantly improved by removing self-interaction and are better than the hybrid B3LYP functional, but fall short of MP2 accuracy. Removing SIE results in substantial improvement in the position of the eigenvalue of the extra electron. The negative of the highest occupied eigenvalue after SIC provides an excellent approximation to the VDE, especially for SIC-PBE where the mean absolute error with respect to CCSD(T) is only 17 meV, the best among all approximations compared in this work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗