Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional approximations”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

ABINIT: Overview and focus on selected capabilities

ABINIT is probably the first electronic-structure package to have been released under an open-source license about 20 years ago. It implements density functional theory, density-functional perturbation theory (DFPT), many-body perturbation theory (GW approximation and Bethe–Salpeter equation), and more specific or advanced formalisms, such as dynamical mean-field theory (DMFT) and the “temperature-dependent effective potential” approach for anharmonic effects. Relying on planewaves for the representation of wavefunctions, density, and other space-dependent quantities, with pseudopotentials or projector-augmented waves (PAWs), it is well suited for the study of periodic materials, although nanostructures and molecules can be treated with the supercell technique. The present article starts with a brief description of the project, a summary of the theories upon which ABINIT relies, and a list of the associated capabilities. It then focuses on selected capabilities that might not be present in the majority of electronic structure packages either among planewave codes or, in general, treatment of strongly correlated materials using DMFT; materials under finite electric fields; properties at nuclei (electric field gradient, Mössbauer shifts, and orbital magnetization); positron annihilation; Raman intensities and electro-optic effect; and DFPT calculations of response to strain perturbation (elastic constants and piezoelectricity), spatial dispersion (flexoelectricity), electronic mobility, temperature dependence of the gap, and spin-magnetic-field perturbation. The ABINIT DFPT implementation is very general, including systems with van der Waals interaction or with noncollinear magnetism. Community projects are also described: generation of pseudopotential and PAW datasets, high-throughput calculations (databases of phonon band structure, second-harmonic generation, and GW computations of bandgaps), and the library libpaw. ABINIT has strong links with many other software projects that are briefly mentioned.

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

A self-interaction-free local hybrid functional: Accurate binding energies vis-à-vis accurate ionization potentials from Kohn-Sham eigenvalues

We present and test a new approximation for the exchange-correlation (xc) energy of Kohn-Sham density functional theory. It combines exact exchange with a compatible non-local correlation functional. The functional is by construction free of one-electron self-interaction, respects constraints derived from uniform coordinate scaling, and has the correct asymptotic behavior of the xc energy density. It contains one parameter that is not determined ab initio. We investigate whether it is possible to construct a functional that yields accurate binding energies and affords other advantages, specifically Kohn-Sham eigenvalues that reliably reflect ionization potentials. Tests for a set of atoms and small molecules show that within our local-hybrid form accurate binding energies can be achieved by proper optimization of the free parameter in our functional, along with an improvement in dissociation energy curves and in Kohn-Sham eigenvalues. However, the correspondence of the latter to experimental ionization potentials is not yet satisfactory, and if we choose to optimize their prediction, a rather different value of the functional's parameter is obtained. We put this finding in a larger context by discussing similar observations for other functionals and possible directions for further functional development that our findings suggest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Technical Note: Statistical dependences between channels in radiochromic film readings. Implications in multichannel dosimetry

Purpose: This note studies the statistical relationships between color channels in radiochromic film readings with flatbed scanners. The same relationships are studied for noise. Finally, their implications for multichannel film dosimetry are discussed. Methods: Radiochromic films exposed to wedged fields of 6 MV energy were read in a flatbed scanner. The joint histograms of pairs of color channels were used to obtain the joint and conditional probability density functions between channels. Then, the conditional expectations and variances of one channel given another channel were obtained. Noise was extracted from film readings by means of a multiresolution analysis. Two different dose ranges were analyzed, the first one ranging from 112 to 473 cGy and the second one from 52 to 1290 cGy. Results: For the smallest dose range, the conditional expectations of one channel given another channel can be approximated by linear functions, while the conditional variances are fairly constant. The slopes of the linear relationships between channels can be used to simplify the expression that estimates the dose by means of the multichannel method. The slopes of the linear relationships between each channel and the red one can also be interpreted as weights in the final contribution to dose estimation. However, for the largest dose range, the conditional expectations of one channel given another channel are no longer linear functions. Finally, noises in different channels were found to correlate weakly. Conclusions: Signals present in different channels of radiochromic film readings show a strong statistical dependence. By contrast, noise correlates weakly between channels. For the smallest dose range analyzed, the linear behavior between the conditional expectation of one channel given another channel can be used to simplify calculations in multichannel film dosimetry.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Kohn–Sham exchange-correlation potentials from second-order reduced density matrices

We describe a practical algorithm for constructing the Kohn–Sham exchange-correlation potential corresponding to a given second-order reduced density matrix. Unlike conventional Kohn–Sham inversion methods in which such potentials are extracted from ground-state electron densities, the proposed technique delivers unambiguous results in finite basis sets. The approach can also be used to separate approximately the exchange and correlation potentials for a many-electron system for which the reduced density matrix is known. The algorithm is implemented for configuration-interaction wave functions and its performance is illustrated with numerical examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Exploring interlayer coupling in the twisted bilayer PtTe 2

We have investigated interlayer interactions in the bilayer PtTe 2 system, which influence the electronic energy bands near the Fermi level. Our diffusion Monte Carlo (DMC) calculations for the high-symmetry bilayer stackings (AA, AB, AC) manifest distinct interlayer binding characteristics among the stacking modes by revealing significantly different interlayer separations depending on the stacking, which is critical to understanding the interlayer coupling of the twisted bilayers consisting of various local stacking arrangements. Furthermore, a comparison between the interlayer separations obtained from DMC and density functional theory (DFT) shows that meta-generalized gradient approximation (GGA)-based van der Waals–DFT results agree with DMC for different layer stackings, including twisted bilayers, but only the ground-state AA stacking matches well with GGA-based DFT predictions. This underscores the importance of accurate exchange-correlation potentials even for capturing the stacking-dependent interlayer binding properties. We further show that the variability in DFT-predicted interlayer separations is responsible for the large discrepancy of band structures in the 21 . 79 ∘ twisted bilayer PtTe 2 , affecting its classification as metallic or insulating. These results demonstrate the importance of obtaining a correct description of stacking-dependent interlayer coupling in modeling delicate bilayer systems at finite twists. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2013

Galaxy correlations and the BAO in a void universe: structure formation as a test of the Copernican Principle

A suggested solution to the dark energy problem is the void model, where accelerated expansion is replaced by Hubble-scale inhomogeneity. In these models, density perturbations grow on a radially inhomogeneous background. This large scale inhomogeneity distorts the spherical Baryon Acoustic Oscillation feature into an ellipsoid which implies that the bump in the galaxy correlation function occurs at different scales in the radial and transverse correlation functions. We compute these for the first time, under the approximation that curvature gradients do not couple the scalar modes to vector and tensor modes. The radial and transverse correlation functions are very different from those of the concordance model, even when the models have the same average BAO scale. This implies that if void models are fine-tuned to satisfy average BAO data, there is enough extra information in the correlation functions to distinguish a void model from the concordance model. We expect these new features to remain when the full perturbation equations are solved, which means that the radial and transverse galaxy correlation functions can be used as a powerful test of the Copernican Principle.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Comparison of Exact and Model Exchange–Correlation Potentials for Molecules

Accurate exchange–correlation (XC) potentials for three-dimensional systems–via solution of the inverse density functional theory (DFT) problem–are now available to test the quality of DFT approximations. Herein, the exact XC potential for seven molecules–dihydrogen at four different bond-lengths, lithium hydride, water, and ortho-benzyne–are computed from full configuration interaction reference densities. These are compared to model XC potentials from nonlocal (B3LYP, HSE06, SCAN0, and M08-HX) and semilocal/local (SCAN, PBE, and PW92) XC functionals. Whereas for most systems, relative errors in the ground-state densities are O(10 –3 –10 –2 ) , the model XC potentials have much higher errors of O(10 –1 –10 0 ) . Among the model XC functionals, SCAN0 offers the best agreement with the exact XC potential, underlining the significance of satisfying exact conditions as well as including nonlocal effects in XC functionals. This work indicates that tests against the exact XC potential will provide a promising new direction for building more accurate XC functionals for DFT.

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

The influence of defects on magnetic properties of fcc-Pu

The influence of vacancies and interstitial atoms on magnetism in Pu is considered in the framework of the density functional theory. The crystal structure relaxation arising due to various types of defects is calculated using the molecular dynamics method with a modified embedded atom model. The local density approximation with explicit inclusion of Coulomb and spin-orbit interactions is applied in matrix invariant form to describe correlation effects in Pu with these types of defects. The calculations show that both vacancies and interstitials give rise to local moments in the f-shell of Pu in good agreement with experimental data for aged Pu. Magnetism appears due to the destruction of a delicate balance between spin-orbit and exchange interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

The structural, electronic and dynamic properties of the L1{sub 2}- type Co{sub 3}Ti alloy

The structural, electronic and dynamic properties of the cubic Co{sub 3}Ti alloy in L1{sub 2} structure have been investigated using a pseudopotential plane wave (PP-PW) method within the generalized gradient approximation proposed by Perdew–Burke–Ernzerhof (GGA-PBE). The structural properties, including the lattice constant, the bulk modulus and its pressure derivative agree reasonably with the previous results. The density of state (DOS), projected density of state (PDOS) and electronic band structure are also reported. The DOS shows that Co{sub 3}Ti alloy has a metallic character since the energy bands cross the Fermi level. The density of states at Fermi level mainly comes from the Co-3d states. Phonon dispersion curves and their corresponding total densities of states were obtained using a linear response in the framework of the density functional perturbation theory. All computed phonon frequencies are no imaginer and thus, Co{sub 3}Ti alloy is dynamically stable. The zone center phonon modes have been founded to be 9.307, 9.626 and 13.891 THz for Co{sub 3}Ti.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Theoretical study of the decomposition pathways and products of C5- perfluorinated ketone (C5 PFK)

Due to the high global warming potential (GWP) and increasing environmental concerns, efforts on searching the alternative gases to SF{sub 6}, which is predominantly used as insulating and interrupting medium in high-voltage equipment, have become a hot topic in recent decades. Overcoming the drawbacks of the existing candidate gases, C5- perfluorinated ketone (C5 PFK) was reported as a promising gas with remarkable insulation capacity and the low GWP of approximately 1. Experimental measurements of the dielectric strength of this novel gas and its mixtures have been carried out, but the chemical decomposition pathways and products of C5 PFK during breakdown are still unknown, which are the essential factors in evaluating the electric strength of this gas in high-voltage equipment. Therefore, this paper is devoted to exploring all the possible decomposition pathways and species of C5 PFK by density functional theory (DFT). The structural optimizations, vibrational frequency calculations and energy calculations of the species involved in a considered pathway were carried out with DFT-(U)B3LYP/6-311G(d,p) method. Detailed potential energy surface was then investigated thoroughly by the same method. Lastly, six decomposition pathways of C5 PFK decomposition involving fission reactions and the reactions with a transition states were obtained. Important intermediate products were also determined. Among all the pathways studied, the favorable decomposition reactions of C5 PFK were found, involving C-C bond ruptures producing Ia and Ib in pathway I, followed by subsequent C-C bond ruptures and internal F atom transfers in the decomposition of Ia and Ib presented in pathways II + III and IV + V, respectively. Possible routes were pointed out in pathway III and lead to the decomposition of IIa, which is the main intermediate product found in pathway II of Ia decomposition. We also investigated the decomposition of Ib, which can undergo unimolecular reactions to give the formation of IV a, IV b and products of CF{sub 3} + CF-CF{sub 3} in pathway IV. Although IV a is dominant to a lesser extent due to its relative high energy barrier, its complicated decomposition pathway V was also studied and CF{sub 3}, C = CF{sub 2} as well as C-CF{sub 3} species were found as the ultimate products. To complete the decomposition of C5 PFK, pathway V I of Ic decomposition was fully explored and the final products were obtained. Therefore, the integrate decomposition scheme of C5 PFK was proposed, which contains six pathways and forty-eight species (including all the reactants, products and transition states). This work is hopeful to lay a theoretical basis for the insulating properties of C5 PFK.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

ωB97M-V: A combinatorially optimized, range-separated hybrid, meta-GGA density functional with VV10 nonlocal correlation

A combinatorially optimized, range-separated hybrid, meta-GGA density functional with VV10 nonlocal correlation is presented in this paper. The final 12-parameter functional form is selected from approximately 10 × 10 9 candidate fits that are trained on a training set of 870 data points and tested on a primary test set of 2964 data points. The resulting density functional, ωB97M-V, is further tested for transferability on a secondary test set of 1152 data points. For comparison, ωB97M-V is benchmarked against 11 leading density functionals including M06-2X, ωB97X-D, M08-HX, M11, ωM05-D, ωB97X-V, and MN15. Encouragingly, the overall performance of ωB97M-V on nearly 5000 data points clearly surpasses that of all of the tested density functionals. Finally, in order to facilitate the use of ωB97M-V, its basis set dependence and integration grid sensitivity are thoroughly assessed, and recommendations that take into account both efficiency and accuracy are provided.

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

Reflection-asymmetric nuclear deformations within the Density Functional Theory

Within the nuclear density functional theory (DFT) we study the effect of reflection-asymmetric shapes on ground-state binding energies and binding energy differences. To this end, we developed the new DFT solver AXIALHFB that uses an approximate second-order gradient to solve the Hartree-Fock-Bogoliubov equations of superconducting DFT with the quasi-local Skyrme energy density functionals. Illustrative calculations are carried out for even-even isotopes of radium and thorium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Strain induced modification in phonon dispersion curves of monolayer boron pnictides

In the frame work of density functional theory, the biaxial strain induced phonon dispersion curves of monolayer boron pnictides (BX, X = N, P, As, and Sb) have been investigated. The electron-ion interactions have been modelled using ultrasoft pseudopotentials while exchange-correlation energies have been approximated by the method of local density approximation in the parameterization of Perdew-Zunger. The longitudinal and transverse acoustic phonon modes of boron pnictide sheets show linear dependency on wave vector k{sup →} while out of plane mode varies as k{sup 2}. The in-plane longitudinal and out of plane transverse optical modes in boron nitride displaying significant dispersion similar to graphene. We have analyzed the biaxial strain dependent behaviour of out of plane acoustic phonon mode which is linked to ripple for four BX sheets using a model equation with shell elasticity theory. The strain induces the hardening of this mode with tendency to become more linear with increase in strain percentage. The strain induced hardening of out of plane acoustic phonon mode indicates the absence of rippling in these compounds. Our band structure calculations for both unstrained and strained 2D h-BX are consistent with previous calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Exchange–correlation entropy from the generalized thermal adiabatic connection

Warm dense matter is a highly energetic phase characterized by strong correlations, thermal effects, and quantum mechanical electrons. Thermal density functional theory is commonly used in simulations of this challenging phase, driving the development of temperature-dependent approximations to the exchange–correlation free energy. Approaches using the adiabatic connection formula are well known at zero temperature and have been recently leveraged at non-zero temperatures as well. In this work, a generalized thermal adiabatic connection (GTAC) formula is proposed, introducing a fictitious temperature parameter. This allows extraction of the exchange–correlation entropy S XC using simulated interaction strength scaling. Here, this procedure uses a Hellmann–Feynman approach to express the exchange–correlation entropy in terms of a temperature- and interaction strength-dependent exchange–correlation potential energy. In addition, analysis of S XC as a function of interaction strength suggests new forms for approximations, and GTAC itself offers a new framework for exploring both the exact and approximate interplay of temperature, density, and interaction strength across a wide range of conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Microscopic origin of the charge transfer in single crystals based on thiophene derivatives: A combined NEXAFS and density functional theory approach

We have investigated the charge transfer mechanism in single crystals of DTBDT-TCNQ and DTBDT-F{sub 4}TCNQ (where DTBDT is dithieno[2,3-d;2′,3′-d′] benzo[1,2-b;4,5-b′]dithiophene) using a combination of near-edge X-ray absorption spectroscopy (NEXAFS) and density functional theory calculations (DFT) including final state effects beyond the sudden state approximation. In particular, we find that a description that considers the partial screening of the electron-hole Coulomb correlation on a static level as well as the rearrangement of electronic density shows excellent agreement with experiment and allows to uncover the details of the charge transfer mechanism in DTBDT-TCNQ and DTBDT-F{sub 4} TCNQ, as well as a reinterpretation of previous NEXAFS data on pure TCNQ. Finally, we further show that almost the same quality of agreement between theoretical results and experiment is obtained by the much faster Z+1/2 approximation, where the core hole effects are simulated by replacing N or F with atomic number Z with the neighboring atom with atomic number Z+1/2.

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

Investigation of the Thermodynamic and Electronic Properties of Double Perovskite Ca{sub 2}CoNbO{sub 6}

In the present work, a self-consistent ab initio calculation using the full-potential linearized augmented plane wave (FP-LAPW) method within the framework of the spin-polarized density functional theory (DFT) was used to study the structural, electronic, and thermodynamic properties of Ca{sub 2}CoNbO{sub 6} double perovskite compound. The generalized gradient approximation (GGA) described by Perdew–Burke–Ernzerhof (PBE) and GGA + U were used. The results obtained for the electronic properties show a ferrimagnetic and half-metallic behavior of the compound. The novelty of our work is the study of the thermodynamic properties of Ca{sub 2}CoNbO{sub 6} double perovskite such as heat capacity and Debye temperature which showed an important effect of pressure compared to the temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Regarding the use and misuse of retinal protonated Schiff base photochemistry as a test case for time-dependent density-functional theory

The excited-state relaxation of retinal protonated Schiff bases (PSBs) is an important test case for biological applications of time-dependent (TD) density-functional theory (DFT). While well-known shortcomings of approximate TD-DFT might seem discouraging for application to PSB relaxation, progress continues to be made in the development of new functionals and of criteria allowing problematic excitations to be identified within the framework of TD-DFT itself. Furthermore, experimental and theoretical ab initio advances have recently lead to a revised understanding of retinal PSB photochemistry, calling for a reappraisal of the performance of TD-DFT in describing this prototypical photoactive system. Here, we re-investigate the performance of functionals in (TD-)DFT calculations in light of these new benchmark results, which we extend to larger PSB models. We focus on the ability of the functionals to describe primarily the early skeletal relaxation of the chromophore and investigate how far along the out-of-plane pathways these functionals are able to describe the subsequent rotation around formal single and double bonds. Conventional global hybrid and range-separated hybrid functionals are investigated as the presence of Hartree-Fock exchange reduces problems with charge-transfer excitations as determined by the Peach-Benfield-Helgaker-Tozer Λ criterion and by comparison with multi-reference perturbation theory results. While we confirm that most functionals cannot render the complex photobehavior of the retinal PSB, do we also observe that LC-BLYP gives the best description of the initial part of the photoreaction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A systematic way for the cost reduction of density fitting methods

We present a simple approach for the reduction of the size of auxiliary basis sets used in methods exploiting the density fitting (resolution of identity) approximation for electron repulsion integrals. Starting out of the singular value decomposition of three-center two-electron integrals, new auxiliary functions are constructed as linear combinations of the original fitting functions. The new functions, which we term natural auxiliary functions (NAFs), are analogous to the natural orbitals widely used for the cost reduction of correlation methods. The use of the NAF basis enables the systematic truncation of the fitting basis, and thereby potentially the reduction of the computational expenses of the methods, though the scaling with the system size is not altered. The performance of the new approach has been tested for several quantum chemical methods. It is demonstrated that the most pronounced gain in computational efficiency can be expected for iterative models which scale quadratically with the size of the fitting basis set, such as the direct random phase approximation. The approach also has the promise of accelerating local correlation methods, for which the processing of three-center Coulomb integrals is a bottleneck.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗