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 37 records · Page 2

Density-functional calculations of the surface tension of liquid Al and Na

Calculations of the surface tensions of liquid Al and Na are described using the full ionic density functional formalism of Wood and Stroud (1983). Surface tensions are in good agreement with experiment in both cases, with results substantially better for Al than those found previously in the gradient approximation. Preliminary minimization with respect to surface profile leads to an oscillatory profile superimposed on a nearly steplike ionic density disribution; the oscillations have a wavellength of about a hardsphere diameter.

Stroud, D.↗

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↗

All-electron density functional calculations for electron and nuclear spin interactions in molecules and solids

The interaction between electronic and nuclear spins in the presence of external magnetic fields can be described by a spin Hamiltonian, with parameters obtained from first principles, electronic structure calculations. We describe an approach to compute these parameters, applicable to both molecules and solids, which is based on density functional theory (DFT) and real-space, all-electron calculations using finite elements (FE). We report results for hyperfine tensors, zero field splitting tensors (spin-spin component), and nuclear quadrupole tensors of a series of molecules and of the nitrogen-vacancy center in diamond. We compare our results with calculations using Gaussian orbitals and plane-wave basis sets, and we discuss their numerical accuracy. Furthermore, we show that calculations based on FE can be systematically converged with respect to the basis set, thus allowing one to establish reference values for the spin Hamiltonian parameters, at a given level of DFT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Density-functional calculations of transport properties in the nondegenerate limit and the role of electron-electron scattering

We compute electrical and thermal conductivities of hydrogen plasmas in the non-degenerate regime using Kohn-Sham Density Functional Theory (DFT) and an application of the Kubo- Greenwood response formula, and demonstrate that for thermal conductivity, the mean-field treatment of the electron-electron (e-e) interaction therein is insufficient to reproduce the weak-coupling limit obtained by plasma kinetic theories. An explicit e-e scattering correction to the DFT is posited by appealing to Matthiessen's Rule and the results of our computations of conductivities with the quantum Lenard-Balescu (QLB) equation. Further motivation of our correction is provided by an argument arising from the Zubarev quantum kinetic theory approach. Significant emphasis is placed on our efforts to produce properly converged results for plasma transport using Kohn-Sham DFT, so that an accurate assessment of the importance and efficacy of our e-e scattering corrections to the thermal conductivity can be made.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

How well do self-interaction corrections repair the overestimation of static polarizabilities in density functional calculations?

Herein, we examine the effect of removing self-interaction error (SIE) on the calculation of molecular polarizabilities in the local spin density (LSDA) and generalized gradient approximations (GGA). To this end, we utilize a database of 132 molecules taken from a recent benchmark study [Hait and Head-Gordon, Phys. Chem. Chem. Phys. , 2018, 20 , 19800] to assess the influence of SIE on polarizabilities by comparing results with accurate reference data. Our results confirm that the general overestimation of molecular polarizabilities by these density functional approximations can be attributed to SIE. However, removing SIE using the Perdew–Zunger self-interaction-correction (PZ-SIC) method, implemented using the Fermi–Löwdin Orbital SIC approach, leads to an underestimation of molecular polarizabilities, showing that PZ-SIC overcorrects when combined with LSDA or GGA. Application of a recently proposed locally scaled SIC [Zope, et al., J. Chem. Phys. , 2019, 151 , 214108] is found to provide more accurate polarizabilities. We attribute this to the ability of the local scaling scheme to selectively correct for SIE in the regions of space where the correction is needed most.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reliable density functional calculations for the electronic structure of thermoelectric material ZnSb

In this paper, we present the results of systematic test calculations for the electronic structure of thermoelectric material ZnSb using a first-principles full-potential all electron computational method. We used a linear combination of atomic orbitals (LACO) formalism, based on density functional theory (DFT). The exchange-correlation interaction potential of the many electron system was described by using a generalized gradient approximation (GGA). We compared the calculated indirect and direct band gaps as well as the effective masses of holes and electrons in ZnSb with experimental measurement results. The calculated indirect band gap of ZnSb is 0.56 eV, which agrees very well with the experimentally measured values of 0.50 eV ~ 0.61 eV. The calculated direct band gap at X point is 0.89 eV. The calculated effective masses of electrons and holes in ZnSb also agree with experimental data. The systematical test calculations as well as the comparisons of the calculated results with experimental measurements show that the obtained electronic structure of ZnSb would be reliable. We did not observe a major deficiency of the first-principles DFT calculation for the electronic structure of ZnSb, using full-potential all electron LACO method. The reported electronic structure of single crystal ZnSb from this work may provide a fundamental knowledge base for further research and applications for this important thermoelectric material.

36 MATERIALS SCIENCE↗

Ground State Structure of BaFeO{sub 3}: Density Functional Theory Calculations

Using density functional theory calculations, the ground state structure of BaFeO{sub 3} (BFO) is investigated with local spin density approximation (LSDA). Cubic, tetragonal, orthorhombic, and rhombohedral types BFO are considered to calculate the formation enthalpy. The formation enthalpies reveal that cubic is the most stable structure of BFO. Small energy difference between the cubic and tetragonal suggests a possible tetragonal BFO. Ferromagnetic (FM) and anitiferromagnetic (AFM) coupling between the Fe atoms show that all the striochmetric BFO are FM. The energy difference between FM and AFM shows room temperature ferromagnetism in cubic BFO in agreement with the experimental work. The LSDA-calculated electronic structures are metallic in all studied crystallographic phases of BFO. Calculations including the Hubbard potential U, i.e. LSDA + U, show that all phases of BFO are half-metallic consistent with the integer magnetic moments. The presence of half-metallicity is discussed in terms of electronic band structures of BFO.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Study of Self-Interaction Errors in Density Functional Calculations of Magnetic Exchange Coupling Constants Using Three Self-Interaction Correction Methods

We examine the role of self-interaction errors (SIE) removal on the evaluation of magnetic exchange coupling constants. In particular we analyze the effect of scaling down the self-interaction-correction (SIC) for three non-empirical density functional approximations (DFAs) namely, the local spin density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation, and recent SCAN family of meta-GGA functionals. To this end, we employ three one-electron SIC methods: Perdew-Zunger [Perdew, J. P.; Zunger, A. Phys. Rev. B, 1981, 23, 5048] SIC, the orbitalwise scaled SIC method [Vydrov, O. A. et al., J. Chem. Phys. 2006, 124, 094108], and the recent local scaling method [Zope, R. R. et al., J. Chem. Phys. 2019, 151, 214108]. We compute the magnetic exchange coupling constants using the spin projection and non projection approaches for sets of molecules composed of dinuclear and polynuclear H-He models, organic radical molecules, and chlorocuprate, and compare these results against accurate theories and experiment. Our results show that for the systems that mainly consist of single electron regions, PZSIC performs well but for more complex organic systems and the chlorcuprates, an overcorrecting tendency of PZSIC combined with the DFAs utilized in this work is more pronounced, and in such cases LSIC with kinetic energy density ratio performs better than PZSIC. Furthermore, analysis of the results in terms of SIC corrections to the density and to the total energy shows that both density and energy correction are required to obtain an improved prediction of magnetic exchange couplings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synthesis, Spectra, and Theoretical Investigations of 1,3,5-Triazines Compounds as Ultraviolet Rays Absorber Based on Time-Dependent Density Functional Calculations and three-Dimensional Quantitative Structure-Property Relationship

A series of 1,3,5-triazines were synthesized and their UV absorption properties were tested. The computational chemistry methods were used to construct quantitative structure-property relationship (QSPR), which was used to computer aided design of new 1,3,5-triazines ultraviolet rays absorber compounds. The experimental UV absorption data are in good agreement with those predicted data using the Time-dependent density functional theory (TD-DFT) [B3LYP/6–311 + G(d,p)]. A suitable forecasting model (R > 0.8, P < 0.0001) was revealed. Predictive three-dimensional quantitative structure-property relationship (3D-QSPR) model was established using multifit molecular alignment rule of Sybyl program, which conclusion is consistent with the TD-DFT calculation. The exceptional photostability mechanism of such ultraviolet rays absorber compounds was studied and confirmed as principally banked upon their ability to undergo excited-state deactivation via an ultrafast excited-state proton transfer (ESIPT). The intramolecular hydrogen bond (IMHB) of 1,3,5-triazines compounds is the basis for the excited state proton transfer, which was explored by IR spectroscopy, UV spectra, structural and energetic aspects of different conformers and frontier molecular orbitals analysis.

74 ATOMIC AND MOLECULAR PHYSICS↗

Self-interaction corrected density functional calculations of Rydberg states of molecular clusters: N,N-dimethylisopropylamine

Theoretical calculations of Rydberg excited states of molecular clusters consisting of N,N-dimethylisopropylamine molecules using a Perdew-Zunger self-interaction corrected energy functional are presented and compared with results of resonant multiphoton ionization measurements. The binding energy of the Rydberg electron in the monomer is calculated to be 2.79 eV and 2.27 eV in the 3s and 3p state, respectively, which compares well with measured values of 2.88 eV and 2.21 eV. Three different stable configurations of the dimer in the ground state were found using an energy functional that includes van der Waals interaction. The lowest ground state energy conformation has the two N-atoms widely separated, by 6.2 Å, while the Rydberg state energy is lowest for a configuration where the N-atoms of the two molecules come close together, separated by 3.7 Å. This conformational change is found to lower the Rydberg electron binding energy by 0.2 eV. The self-interaction corrected functional gives a highly localized hole on one of the two molecules, unlike results obtained using the PBE functional or the hybrid B3LYP functional which give a delocalized hole. For the trimer, the self-interaction corrected calculation gives a Rydberg electron binding energy lowered further by 0.13 eV as compared with the dimer. The calculated results compare well with trends observed in experimental measurements. The reduction of the Rydberg electron binding energy with cluster size can be ascribed to an effective delocalization of the positive charge of the hole by the induced and permanent dipole moments of the neighboring molecules. A further decrease observed to occur on a time scale of tens of ps can be ascribed to a structural rearrangement of the clusters in the Rydberg state where molecules rotate to orient their dipoles in response to the formation of the localized hole.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fragment approach to constrained density functional theory calculations using Daubechies wavelets

In a recent paper, we presented a linear scaling Kohn-Sham density functional theory (DFT) code based on Daubechies wavelets, where a minimal set of localized support functions are optimized in situ and therefore adapted to the chemical properties of the molecular system. Thanks to the systematically controllable accuracy of the underlying basis set, this approach is able to provide an optimal contracted basis for a given system: accuracies for ground state energies and atomic forces are of the same quality as an uncontracted, cubic scaling approach. This basis set offers, by construction, a natural subset where the density matrix of the system can be projected. In this paper, we demonstrate the flexibility of this minimal basis formalism in providing a basis set that can be reused as-is, i.e., without reoptimization, for charge-constrained DFT calculations within a fragment approach. Support functions, represented in the underlying wavelet grid, of the template fragments are roto-translated with high numerical precision to the required positions and used as projectors for the charge weight function. We demonstrate the interest of this approach to express highly precise and efficient calculations for preparing diabatic states and for the computational setup of systems in complex environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗