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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Educating local radial basis functions using the highest gradient of interest in three dimensional geometries: EDUCATING LOCALIZED RBF-SCHEMES

We present a novel methodology to effectively localize radial basis function approximation methods in three dimensions. The local scheme requires shape parameter-dependent functions that can be used to approximate gradients of scattered data and to solve partial differential operators. The optimum shape parameter is obtained from the highest gradient of interest, where a known analytical function, when boundary conditions are not present, or a shape parameter-free global approximation are used to educate the localized scheme. The later option is applicable to problems where the operator needs to be solved multiple times, like in time evolution or stochastic integration. Past shape parameter's optimizations, for two-dimensional domains, based on the condition number of the interpolant matrix, were unable to provide satisfactory approximations. The applicability of our method is illustrated in the context of an analytical expression interpolation and during a Ginzburg-Landau relaxation of a free energy functional. In general, the optimum shape parameter depends on geometry, node distribution, and density, whereas the approximation errors decrease as the node density and the local stencil size increase. The effective localization of radial basis functions motivates its use in moving boundary problems and accelerates solutions through sparse matrix solvers.

free energy functional↗
DOE OSTI2014

Structural and electronic properties of Au{sub n−x}Pt{sub x} (n = 2–14; x ⩽ n) clusters: The density functional theory investigation

The structural evolutions and electronic properties of bimetallic Au{sub n–x}Pt{sub x} (n = 2–14; x ⩽ n) clusters are investigated by using the density functional theory (DFT) with the generalized gradient approximation (GGA). The monatomic doping Au{sub n–1}Pt clusters are emphasized and compared with the corresponding pristine Au{sub n} clusters. The results reveal that the planar configurations are favored for both Au{sub n–1}Pt and Au{sub n} clusters with size up to n = 13, and the former often employ the substitution patterns based on the structures of the latter. The most stable clusters are Au{sub 6} and Au{sub 6}Pt, which adopt regular planar triangle (D{sub 3h}) and hexagon-ring (D{sub 6h}) structures and can be regarded as the preferential building units in designing large clusters. For Pt-rich bimetallic clusters, their structures can be obtained from the substitution of Pt atoms by Au atoms from the Pt{sub n} structures, where Pt atoms assemble together and occupy the center yet Au atoms prefer the apex positions showing a segregation effect. With respect to pristine Au clusters, Au{sub n}Pt clusters exhibit somewhat weaker and less pronounced odd-even oscillations in the highest occupied and lowest unoccupied molecular-orbital gaps (HOMO-LUMO gap), electron affinity (EA), and ionization potential (IP) due to the partially released electron pairing effect. The analyses of electronic structure indicate that Pt atoms in AuPt clusters would delocalize their one 6s and one 5d electrons to contribute the electronic shell closure. The sp-d hybridizations as well as the d-d interactions between the host Au and dopant Pt atoms result in the enhanced stabilities of AuPt clusters.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Assessing the source of error in the Thomas–Fermi–von Weizsäcker density functional

We investigate the source of error in the Thomas–Fermi–von Weizsäcker (TFW) density functional relative to Kohn–Sham density functional theory (DFT). In particular, through numerical studies on a range of materials, for a variety of crystal structures subject to strain and atomic displacements, we find that while the ground state electron density in TFW orbital-free DFT is close to the Kohn–Sham density, the corresponding energy deviates significantly from the Kohn–Sham value. Here, we show that these differences are a consequence of the poor representation of the linear response within the TFW approximation for the electronic kinetic energy, confirming conjectures in the literature. In so doing, we find that the energy computed from a non-self-consistent Kohn–Sham calculation using the TFW electronic ground state density is in very good agreement with that obtained from the fully self-consistent Kohn–Sham solution.

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

Dirac Nodes and Magnetic Order in M 2 X 2 Transition‐Metal Chalcogenides

In this study, we perform a computational analysis of the M 2 X 2 transition‐metal chalcogenides (TMCs). Using density functional theory with a spin‐polarized generalized gradient approximation, we examine the magnetic and electronic properties for the antiferromagnetic and ferromagnetic states with M = Cr, Mn, and Fe and X = S and Se. After optimizing the geometric structure for stability, we examine the spin‐polarized electronic structure, density of states, and Mulliken population. It is discovered that these materials are quasi‐two‐dimensional honeycomb lattices with metallic antiferromagnetic ground states. The structures consist of a distorted tetrahedral crystal‐field symmetry that has a distinct magnetic moment. An analysis of the electronic structure shows the presence of nodal points that resemble Dirac nodes for all cases, which leads to the possibility of the realization of magnetic Dirac materials.

La Martina, Thomas↗
NASA NTRS1985

Local density maxima - Progenitors of structure

Beginning with a short review of the observational data and of some previous theoretical ideas and numerical simulations, the question is addressed of how the large-scale structure that emerges around local density maxima depends on the nature of the primordial density perturbations field. The density contrast profile around local maxima is given, to a good approximation, by the primordial two-point correlation function. The mean number density of objects of a given core mass is calculated as a function of the primordial power spectrum, p(k). In an open universe, rich clusters should have halos steeper than galactic haloes. The observed structure is found to be consistent with omega-sub-zero less than 1.0 and n = -1.

Hoffman, Y.↗
DOE OSTI2015

The interaction between hexagonal boron nitride and water from first principles

The use of hexagonal boron nitride (h-BN) in microfluidic and nanofluidic applications requires a fundamental understanding of the interaction between water and the h-BN surface. A crucial component of the interaction is the binding energy, which is sensitive to the treatment of electron correlation. In this work, we use state of the art quantum Monte Carlo and quantum chemistry techniques to compute the binding energy. Compared to high-level many-body theory, we found that the second-order Møller-Plesset perturbation theory captures the interaction accurately and can thus be used to develop force field parameters between h-BN and water for use in atomic scale simulations. On the contrary, density functional theory with standard dispersion corrections tends to overestimate the binding energy by approximately 75%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Many-body perturbation theory with hybrid density functional theory starting points accelerated by adaptively compressed exchange

We report on the use of the adaptively compressed exchange (ACE) operator to accelerate many-body perturbation theory (MBPT) calculations, including G 0 W 0 and the Bethe–Salpeter equation (BSE), for hybrid density functional theory starting points. We show that by approximating the exact exchange operator with the low-rank ACE operator, substantial computational savings can be achieved with systematically controllable errors in the quasiparticle energies computed with full-frequency G 0 W 0 and the optical absorption spectra and vertical excitation energies computed by solving the BSE within density matrix perturbation theory. Our implementation makes use of the ACE-accelerated electronic Hamiltonian to carry out both G 0 W 0 and BSE without explicitly computing empty states. We show the robustness of the approach and present the computational gains obtained on both the central processing unit and graphics processing unit nodes. In conclusion, our work will facilitate the exploration and evaluation of fine-tuned hybrid starting points aimed at enhancing the accuracy of MBPT calculations without involving computationally demanding self-consistency in Hedin’s equations.

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

Nuclear kinetic density from ab initio theory

The nuclear kinetic density is one of many fundamental, nonobservable quantities in density functional theory (DFT) dependent on the nonlocal nuclear density. Often, approximations may be made when computing the density that may result in spurious contributions in other DFT quantities. With the ability to compute the nonlocal nuclear density from ab initio wave functions, it is now possible to estimate effects of such spurious contributions. Herein, we derive the kinetic density using ab initio nonlocal scalar one-body nuclear densities computed within the no-core shell model (NCSM) approach, utilizing two- and three-nucleon chiral interactions as the sole input. The ability to compute translationally invariant nonlocal densities allows us to gauge the impact of the spurious center-of-mass (c.m.) contributions in DFT quantities, such as the kinetic density, and provide ab initio insight into refining energy density functionals. The nonlocal nuclear densities are derived from the NCSM one-body densities calculated in second quantization. We present a review of c.m. contaminated and translationally invariant nuclear densities. We then derive an analytic expression for the kinetic density using these nonlocal densities, producing an ab initio kinetic density. The ground-state nonlocal densities of He 4 , 6 , 8 , C 12 , and O 16 are used to compute the kinetic densities of the aforementioned nuclei. The impact of c.m. removal techniques in the density are discussed and compared to a procedure applied in DFT. The results of this work can be extended to other fundamental quantities in DFT. The use of a general nonlocal density allows for the calculation of fundamental quantities taken as input in theories such as DFT. This allows benchmarking c.m. removal procedures and provides a bridge for comparison between ab initio and DFT many-body techniques.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Perspectives on van der Waals Density Functionals: The Case of TiS 2

The van der Waals interaction is of foundational importance for a wide variety of physical systems. In particular, van der Waals forces lie at the heart of potential device technologies that may be realized from the functional organization of layered two-dimensional (2D) nanomaterials. Furthermore, for intermediate to large-scale applications modeling, van der Waals density functionals have become the de facto choice for first-principles calculations. In particular, the vdW-DF family of functionals have provided a systematic approach to this theoretically challenging problem. While much progress has been made, there remains room for improvement in the microscopic description of vdW forces from these density functionals. In this work, we compute benchmark results for the binding energy and the electronic density response to binding in TiS2 via accurate diffusion quantum Monte Carlo calculations. We compare these benchmark data to results obtained from local, semilocal, and van der Waals functionals. In particular, we gauge the quality of the original vdW-DF/vdW-DF2 functionals, as well as updated variants such as vdW-DF-C09, vdW-DF-optB88, vdW-DF-optB86b, and vdW-DF2-B86R. We find a close relationship between the accuracy of predicted interlayer separation distances and binding energies for TiS 2 , with the vdW-DF-optB88 functional performing very well in terms of both quantities. In general, the more recently developed functionals are systematic improvements over older ones. However, when considering the response of the electron density to binding, we find that local-density approximation (LDA) and PBEsol generally outperform the vdW-DF functionals in describing the interlayer charge accumulation with vdW-DF-C09 variants performing the best overall.

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

Investigation of New d{sup 0} Half-metallic Heusler Alloy LiCsN{sub 2} Using First-Principle Calculations

Using the self-consistent full-potential linearized augmented plane wave (FPLAPW) method, we predicated a new d{sup 0} half-metallic Heusler compound LiCsN{sub 2} with integer total magnetic moment 4.00 μ B per formula unit based on density functional theory (DFT). Our calculation shows that this alloy has a huge half-metallic gap (1.06 eV) and a wide band gap (4.52 eV). We found the LiCsN{sub 2} alloy is a robust half-metallic material with respect to the lattice compression. It maintains the half-metallic character until the lattice constant is contracted to 6.10 Å. The Curie temperature of LiCsN{sub 2} was estimated to be 6145.4 K within the mean field approximations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Bayesian Monte Carlo Evaluation Framework for Cross Sections Nuclear Data and Integral Benchmark Experiments

The new Bayesian Monte Carlo (MC) evaluation framework described in this abstract has been conceived as an attempt to improve nuclear data evaluations of differential crosssection data by removing the following two approximations conventionally employed for nuclear data evaluations: all probability density functions (PDFs) of all data and model parameters, both prior and posterior, are assumed to be normal (i.e., Gaussian) PDFs, and all uncertainties and covariances are propagated using a linear approximation. With these approximations removed, the Bayesian MC (BMC) framework could be used to account for nonlinear effects and would enable improved evaluations of differential cross sections and IBE data that are presently performed based on the assumptions itemized above. The BMC would also improve upon the uniform sampling of IBE parameters from within ranges defined by their evaluated uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

On the theory of electric double layer with explicit account of a polarizable co-solvent

We present a continuation of our theoretical research into the influence of co-solvent polarizability on a differential capacitance of the electric double layer. We formulate a modified Poisson-Boltzmann theory, using the formalism of density functional approach on the level of local density approximation taking into account the electrostatic interactions of ions and co-solvent molecules as well as their excluded volume. We derive the modified Poisson-Boltzmann equation, considering the three-component symmetric lattice gas model as a reference system and minimizing the grand thermodynamic potential with respect to the electrostatic potential. We apply present modified Poisson-Boltzmann equation to the electric double layer theory, showing that accounting for the excluded volume of co-solvent molecules and ions slightly changes the main result of our previous simplified theory. Namely, in the case of small co-solvent polarizability with its increase under the enough small surface potentials of electrode, the differential capacitance undergoes the significant growth. Oppositely, when the surface potential exceeds some threshold value (which is slightly smaller than the saturation potential), the increase in the co-solvent polarizability results in a differential capacitance decrease. However, when the co-solvent polarizability exceeds some threshold value, its increase generates a considerable enhancement of the differential capacitance in a wide range of surface potentials. We demonstrate that two qualitatively different behaviors of the differential capacitance are related to the depletion and adsorption of co-solvent molecules at the charged electrode. We show that an additive of the strongly polarizable co-solvent to an electrolyte solution can shift significantly the saturation potential in two qualitatively different manners. Namely, a small additive of strongly polarizable co-solvent results in a shift of saturation potential to higher surface potentials. On the contrary, a sufficiently large additive of co-solvent shifts the saturation potential to lower surface potentials. We obtain that an increase in the co-solvent polarizability makes the electrostatic potential profile longer-ranged. However, increase in the co-solvent concentration in the bulk leads to non-monotonic behavior of the electrostatic potential profile. An increase in the co-solvent concentration in the bulk at its sufficiently small values makes the electrostatic potential profile longer-ranged. Oppositely, when the co-solvent concentration in the bulk exceeds some threshold value, its further increase leads to decrease in electrostatic potential at all distances from the electrode.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1976

A kinetic-theory approach to turbulent chemically reacting flows

The paper examines the mathematical and physical foundations for the kinetic theory of reactive turbulent flows, discussing the differences and relation between the kinetic and averaged equations, and comparing some solutions of the kinetic equations obtained by the Green's function method with those obtained by the approximate bimodal method. The kinetic method described consists essentially in constructing the probability density functions of the chemical species on the basis of solutions of the Langevin stochastic equation for the influence of eddies on the behavior of fluid elements. When the kinetic equations are solved for the structure of the diffusion flame established in a shear layer by the bimodal method, discontinuities in gradients of the mean concentrations at the two flame edges appear. This is a consequence of the bimodal approximation of all distribution functions by two dissimilar half-Maxwellian functions, which is a very crude approximation. These discontinuities do not appear when the solutions are constructed by the Green's function method described here.

Chung, P. M.↗
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 OSTI2017

Assessment of electronic structure methods for the determination of the ground spin states of Fe( ii ), Fe( iii ) and Fe( iv ) complexes

Our ability to understand and simulate the reactions catalyzed by iron depends strongly on our ability to predict the relative energetics of spin states. In this work, we studied the electronic structures of Fe 2+ ion, gaseous FeO and 14 iron complexes using Kohn–Sham density functional theory with particular focus on determining the ground spin state of these species as well as the magnitudes of relevant spin-state energy splittings. The 14 iron complexes investigated in this work have hexacoordinate geometries of which seven are Fe(II), five are Fe(III) and two are Fe(IV) complexes. These are calculated using 20 exchange–correlation functionals. In particular, we use a local spin density approximation (LSDA) – GVWN5, four generalized gradient approximations (GGAs) – BLYP, PBE, OPBE and OLYP, two non-separable gradient approximations (NGAs) – GAM and N12, two meta-GGAs – M06-L and M11-L, a meta-NGA – MN15-L, five hybrid GGAs – B3LYP, B3LYP*, PBE0, B97-3 and SOGGA11-X, four hybrid meta-GGAs – M06, PW6B95, MPW1B95 and M08-SO and a hybrid meta-NGA – MN15. The density functional results are compared to reference data, which include experimental results as well as the results of diffusion Monte Carlo (DMC) calculations and ligand field theory estimates from the literature. For the Fe 2+ ion, all functionals except M11-L correctly predict the ground spin state to be quintet. However, quantitatively, most of the functionals are not close to the experimentally determined spin-state splitting energies. For FeO all functionals predict quintet to be the ground spin state. For the 14 iron complexes, the hybrid functionals B3LYP, MPW1B95 and MN15 correctly predict the ground spin state of 13 out of 14 complexes and PW6B95 gets all the 14 complexes right. Here, the local functionals, OPBE, OLYP and M06-L, predict the correct ground spin state for 12 out of 14 complexes. Two of the tested functionals are not recommended to be used for this type of study, in particular M08-SO and M11-L, because M08-SO systematically overstabilizes the high spin state, and M11-L systematically overstabilizes the low spin state.

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

Predicting band gaps and band-edge positions of oxide perovskites using density functional theory and machine learning

Density functional theory (DFT) within the local or semilocal density approximations, i.e., the local density approximation (LDA) or generalized gradient approximation (GGA), has become a workhorse in the electronic structure theory of solids, being extremely fast and reliable for energetics and structural properties, yet remaining highly inaccurate for predicting band gaps of semiconductors and insulators. The accurate prediction of band gaps using first-principles methods is time consuming, requiring hybrid functionals, quasiparticle GW, or quantum Monte Carlo methods. Efficiently correcting DFT-LDA/GGA band gaps and unveiling the main chemical and structural factors involved in this correction is desirable for discovering novel materials in high-throughput calculations. In this direction, we, in this study, use DFT and machine learning techniques to correct band gaps and band-edge positions of a representative subset of ABO 3 perovskite oxides. Relying on the results of HSE06 hybrid functional calculations as target values of band gaps, we find a systematic band-gap correction of ~1.5 eV for this class of materials, where ~1 eV comes from downward shifting the valence band and ~0.5 eV from uplifting the conduction band. The main chemical and structural factors determining the band-gap correction are determined through a feature selection procedure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Long-lived K isomer and enhanced γ vibration in the neutron-rich nucleus 172 Dy: Collectivity beyond double midshell

The level structure of 172Dy has been investigated for the first time by means of decay spectroscopy following in-flight fission of a 238U beam. A long-lived isomeric state with T1/2 = 0.71(5) s and Kπ = 8- has been identified at 1278 keV, which decays to the ground-state and γ -vibrational bands through hindered electromagnetic transitions, as well as to the daughter nucleus 172Ho via allowed β decays. The robust nature of the Kπ = 8- isomer and the ground-state rotational band reveals an axially-symmetric structure for this nucleus. Meanwhile, the γ -vibrational levels have been identified at unusually low excitation energy compared to the neighboring well-deformed nuclei, indicating the significance of the microscopic effect on the non-axial collectivity in this doubly mid-shell region. The underlying mechanism of enhanced γ vibration is discussed in comparison with the deformed Quasiparticle Random-Phase Approximation based on a Skyrme energy-density functional.

172Dy↗
DOE OSTI2023

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute–solvent interactions in modulating the excited-state dynamics. Here, to tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV–vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC) 5 Fe III (μ-CN)Ru II (pyridine) 4 (μ-NC)Fe III (CN) 5 ] 4– ). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute–solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm–Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗