Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density gradient function”

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

Comparing GGA, GGA+ U , and meta-GGA functionals for redox-dependent binding at open metal sites in metal–organic frameworks

Metal–organic frameworks (MOFs) with open metal sites have been widely investigated for the selective adsorption of small molecules via redox mechanisms where charge transfer can take place between the binding site and the adsorbate of interest. Quantum-chemical screening methods based on density functional theory have emerged as a promising route to accelerate the discovery of MOFs with enhanced binding affinities toward various adsorbates. However, the success of this approach is linked to the accuracy of the underlying density functional approximations (DFAs). In this work, we compare commonly used generalized gradient approximation (GGA), GGA+U, and meta-GGA exchange-correlation functionals in modeling redox-dependent binding at open metal sites in MOFs using O 2 and N 2 as representative small molecules. We find that the self-interaction error inherent to the widely used Perdew, Burke, and Ernzerhof (PBE) GGA predicts metal sites that are artificially redox-active, as evidenced by their strong binding affinities, short metal–adsorbate bond distances, and large degree of charge transfer. The incorporation of metal-specific, empirical Hubbard U corrections based on the transition metal oxide literature systematically reduces the redox activity of the open metal sites, often improving agreement with experiment. Additionally, the binding behavior shifts from strong chemisorption to weaker physisorption as a function of U. The M06-L meta-GGA typically predicts binding energies between those of PBE-D3(BJ) and PBE-D3(BJ)+U when using empirically derived U values from the transition metal oxide literature. Despite the strong sensitivity of the binding affinities toward a given DFA, the GGA, GGA+U, and meta-GGA approaches often yield the same qualitative trends and structure–property relationships.

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

A Drift‐Kinetic Method for Obtaining Gradients in Plasma Properties From Single‐Point Distribution Function Data

Abstract In this paper, we derive a new drift‐kinetic method for estimating gradients in the plasma properties through a velocity space distribution at a single point. The gradients are intrinsically related to agyrotropic features of the distribution function. This method predicts the gradients in the magnetized distribution function and can predict gradients of arbitrary moments of the gyrotropic background distribution function. The method allows for estimates on density and pressure gradients on the scale of a Larmor radius, proving to resolve smaller scales than any method currently available to spacecraft. The model is verified with a set of fully kinetic VPIC particle‐in‐cell simulations.

Wetherton, B. A.↗
DOE OSTI2021

Kohn-Sham Solver (KSSOLV) v2.0

KSSOLV is a MATLAB toolbox for solving Kohn-Sham density functional theory based electronic structure eigenvalue problems. It uses an object oriented features of MATLAB to represent atom, molecules, wavefunctions and Hamiltonians and their operations. It is designed to make it easier for users to prototype and test new algorithms for solving the Kohn-Sham problem. KSSOLV2.0 contains significant improvement over the original KSSOLV described in a paper published in ACM Transaction on Mathematical Software (attached). In addition to performing ground state calculation for small molecules, it can also perform geometry optimization for both molecules and solids. It uses standard pseudopotentials and implements local density approximation, generalized gradient approximation and hybrid functionals. Future releases will also include time-dependent DFT and post DFT calculations such as the GW quasi-particle energy calculation and Bethe-Salpeter equation solver for optical absorption.

Yang, Chao↗
DOE OSTI2024

Tailoring the Electronic Structures and Spectral Properties of ZnO with Irradiation Defects Generated Under Intense Electronic Excitation: A Combined Experimental and DFT Approach

Here, the relationship between the composition of the internal defect states, spectral properties, and correlated electronic structures of wurtzite zinc oxide (ZnO) crystals under 645 MeV Xe 35+ irradiation is systematically investigated, employing experimental characterizations combined with first-principle calculations. Based on the ion irradiation-induced thermal expansion and relaxation processes, the high concentration of vacancy/interstitial defects produced from the transient disordered phase in molten track states trigger photoelectric changes, as follows: i) the generation of internal defect states effectively reduces the intrinsic bandgap (3.25 eV → 2.66 eV); ii) a large number of defective active sites inhibits the recombination between electron–hole pairs, causing dark conductance and photoconductance to increase with increasing damage levels until optimal fluence is achieved. Based on the density functional theory (DFT) with the GGA + U (GGA = generalized gradient approximation) method, the defective models associated with the different electronic structures, density of states, formation energy, and the nature of the chemical bonding are established. The narrowing of the bandgap observed experimentally and the enhancement of carrier concentration originating from the internal electron defect states are qualitatively verified, therefore laying the foundation for designing future nanoscale photoelectronic devices and microelectronics applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fully thermal meta-GGA exchange correlation free-energy density functional

The application of density functional theory to materials in the warm dense matter regime has motivated the development of exchange-correlation functionals which incorporate proper, explicit temperature dependence. Previous work has yielded fully-thermal exchange-correlation free energy functionals at the local density approximation (LDA) and generalized gradient approximation (GGA) levels of refinement. Recently an additive thermal correction scheme was utilized to construct a meta-GGA exchange-correlation (XC) functional in which thermal effects are treated at the GGA level. Here, the f TSCAN free-energy XC functional presented here includes thermal effects through the meta-GGA level in the context of the SCAN (strongly constrained and appropriately normed) ground-state functional. The f TSCAN functional provides generality while achieving similar performance to a thermal GGA functional at high temperatures, e.g. pressures within 1% of path integral Monte Carlo simulations of warm dense hydrogen, and a significant improvement over ground-state functionals. At low temperatures, f TSCAN demonstrates improvements in accuracy relative to lower-level and deorbitalized functionals, indicating that calculations using f TSCAN may be expected to perform well across experimentally relevant densities and pressures.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Electronic transport, thermal transport, thermal expansion, and magnetization in the strongly correlated metal LaNi⁢O 3

Perovskite structured LaNiO 3 is a strongly correlated metal with intriguing thermal and magnetic properties. The volume dependence of calculated and measured physical properties can add additional critical information to develop a more in-depth understanding of this strongly correlated phenomenon. Taking advantage of recent single crystal LaNiO 3 growth using the floating-zone method, we have measured the thermal expansion, the magnetostriction, and the pressure dependence of the magnetic susceptibility, which then allows derivation of the Grüneisen parameters γ e = $\frac{dlnN(E_F)}{d lnV}$, γ χ = $\frac{dlnχ}{d lnV}$, as well as of electric and thermal transport properties. We simulate the volume dependence of structural and magnetic properties using Density Functional Theory calculations at the Generalized Gradient Approximation level. A large discrepancy between experimental values and calculated ones suggests that strong correlations are likely to be dynamic in nature. This study also provides a side-by-side comparison of measurements in single crystal and polycrystalline samples of LaNiO 3 to elucidate intrinsic materials properties. A broad hump at high temperatures in the temperature dependence of magnetization found in the single crystal sample of LaNiO 3 has been rationalized by a model that includes the influence of electron correlations on the Landau diamagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Magnetic contributions to corundum-eskolaite and corundum-hematite phase equilibria: A DFT cluster expansion study

Magnetic contributions have the potential to significantly influence predicted phase stability within alloy and mineral mixing phase diagrams, yet have been historically challenging to incorporate due to a significant increase to phase space sampling. Here, in this work, we employ a computational protocol that includes spin orientation as an additional configurational component within multi-component cluster expansions between magnetic and non-magnetic metal oxide alloys [calculated using density functional theory (DFT) and the generalized gradient approximation]. This approach was used to determine the effect of magnetic contributions to corundum-eskolaite and corundum-hematite phase equilibria from first principles. Two-component cluster expansions of the magnetic components of eskolaite and hematite were first performed showing the ability of this method to properly calculate their respective magnetic properties. Two-component cluster expansions were then performed for non-magnetic Al(III) and ferromagnetic Cr(III) and Fe(III), and phase diagrams were calculated for later comparison. Finally, a non-magnetic Al(III) and “up” and “down” magnetic configurations for anti-ferromagnetic Cr(III) and Fe(III) were performed. Magnetic contributions to the calculated phase diagram for the corundumeskolaite system were shown to be inconsequential but are vital for accurate determination of the corundum-hematite solvus.

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

The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm–Dancoff approximation

The primary computational challenge when simulating nonadiabatic ab initio molecular dynamics is the unfavourable compute costs of electronic structure calculations with molecular size. Simple electronic structure theories, like time-dependent density functional theory within the Tamm–Dancoff approximation (TDDFT/TDA), alleviate this cost for moderately sized molecular systems simulated on realistic time scales. Although TDDFT/TDA does have some limitations in accuracy, an appealing feature is that, in addition to including electron correlation through the use of a density functional, the cost of calculating analytic nuclear gradients and nonadiabatic coupling vectors is often computationally feasible even for moderately sized basis sets. Here in this work, some of the benefits and limitations of TDDFT/TDA are discussed and analysed with regard to its applicability as a ‘back-end’ electronic structure method for the symmetric quasi-classical Meyer–Miller model (SQC/MM). In order to investigate the benefits and limitations of TDDFT/TDA, SQC/MM is employed to predict and analyse a prototypical example of excited-state hydrogen transfer in gas-phase malonaldehyde. Then, the ring-opening dynamics of selenophene are simulated, which highlight some of the deficiencies of TDDFT/TDA. Additionally, some new algorithms are proposed that speed up the calculation of analytic nuclear gradients and nonadiabatic coupling vectors for a set of excited electronic states.

molecular dynamics↗
DOE OSTI2020

Localization in the SCAN meta-generalized gradient approximation functional leading to broken symmetry ground states for graphene and benzene

Density functional theory calculations play a central role in understanding chemical and solid-state systems. Progress depends on density functionals that accurately reproduce both energies, for thermochemistry, and properly describe ground states and other properties that are of interest. The Cr dimer, benzene and graphene are especially important benchmark systems for quantum chemistry and condensed matter physics. The Strongly Constrained and Appropriately Normed (SCAN) functional, which is an advanced meta-generalized gradient approximation functional that significantly improves molecular energies is shown to perform poorly for the Cr dimer. This is connected with its poor performance for itinerant solid-state magnets and is a consequence of over localization of electrons, thus illustrating an analogy between the Cr dimer and itinerant magnets. The Cr dimer is a notoriously difficult system for density functionals. However, we additionally find that SCAN predicts an incorrect symmetry broken ground state for 2D graphene and for the benzene molecule, which is surprising considering that ground states of these are known to be well described even by the simplest local density approximation. We show that SCAN overly favors localized spin polarized states, which is a serious deficiency of this approach. Thus, the challenge of finding density functionals that accurately treat both localized and delocalized electronic systems remains.

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

Effect of XC functionals and dispersion corrections on the DFT‐computed structural and vibrational properties of SrCl 2 –NaCl and ZrF 4 –LiF

Density functional theory (DFT) calculations were performed to examine the impact of exchange–correlation (XC) functionals and van der Waals corrections (specifically the D3 method) on the structural and vibrational properties of the SrCl 2 –NaCl and ZrF 4 –LiF salt systems. Multiple XC functionals, including the local density approximation (LDA), the generalized gradient approximation using the Perdew–Burke–Ernzerhof (PBE) model, and its modified form suitable for solids (PBEsol), the dispersion-corrected PBE-D3 and PBEsol-D3, were considered. Of these functionals, LDA was found to exhibit the highest degree of error, while PBEsol and PBE-D3 displayed the least error. Underestimated lattice parameters compared with experimental values were observed to result in higher force constants, leading to an overprediction of vibrational frequencies. Conversely, an overestimation of lattice parameters was associated with lower vibrational frequencies. The methodology presented in this study yielded results that are in good agreement with experiment, irrespective of the method (finite differences vs. density functional perturbation theory) employed for calculating infrared and Raman spectra. It was further demonstrated that for alkali halides with weak Raman scattering, utilizing a supercell constructed from primitive cells better predicts Raman features than does the use of conventional cells.

Raman↗
DOE OSTI2020

Thermal hybrid exchange-correlation density functional for improving the description of warm dense matter

Finite-temperature density functional theory (FT-DFT) has become a standard tool for first principles calculations of the properties of warm dense matter (WDM) relevant to high-energy density physics (HEDP) applications. In this work we introduce theoretical grounds of thermal hybrid exchange-correlation (XC) functionals within the generalized Mermin-Kohn-Sham scheme for an improved description of WDM. Building on the previously developed KSDT (Karasiev-Sjostrom-Dufty-Trickey) [Karasiev et al., Phys. Rev. Lett. 112, 076403 (2014)] local density approximation (LDA) and the KDT16 (Karasiev-Dufty-Trickey 2016) [Karasiev et al., Phys. Rev. Lett. 120, 076402 (2018)] generalized-gradient approximation (GGA) XC free-energy density functionals, we construct a new thermal hybrid XC functional, referred to here as KDT0. The KDT0 model at low temperature reduces to the popular ground-state PBE0 hybrid due to properties of the used KDT16 density functional approximation. Application to static calculations of electronic band gap and band structure at a wide range of temperatures for various systems of interest to HEDP show that KDT0 provides a significant improvement to the lower LDA and GGA rung XC functionals and to the ground-state PBE0 hybrid.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Spiral Galaxies Flat Rotational Velocity Curve Explained by the Constant Group Velocity of a Nonlinear Density Wave

The rotation velocity curves of stars in galaxies, the motions of pairs of galaxies, and the behavior of galaxies in clusters and super-clusters all indicate that there is a lack of mass on different scales in the universe. In this paper, we derive the expression for rotational velocity using the nonlinear density wave theory considering only stellar components and we show that such theory can support the observed flat rotational velocity curve due to the main property of the soliton wave, which is a constant group velocity of the wave. The surface mass density (SMD) function, used to derive gravitational potential gradient and rotational velocity, is not assumed but rather derived as a solution of the nonlinear Srödinger equation, on the contrary to the widely used, in the literature, exponential disk approximation. Three parameters relevant to the curve shape are the intensities of equilibrium SMD, the amplitude of the wave, and total angular velocity or differential rotation, equivalently. Since the shape of the rotational velocity is highly sensitive to the mentioned parameters, this theory eventually provides a method for a very accurate estimation of galaxy mass and angular velocity as well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Tunable noninteracting free-energy density functionals for high-energy-density physics applications

In this work, we introduce the concept of a tunable noninteracting free-energy density functional and present two examples realized: (i) via a simple one-parameter convex combination of two existing functionals and (ii) via the construction of a generalized gradient approximation (GGA) enhancement factor that contains one free parameter and is designed to satisfy a set of incorporated constraints. Functional (i), constructed as a combination of the local Thomas–Fermi and a pseudopotential-adapted GGA for the noninteracting free-energy, has already demonstrated its practical usability for establishing the high temperature end of the equation of state of deuterium [Phys. Rev. B 104, 144104 (2021)] and CHON resin [Phys. Rev. E 106, 045207 (2022)] for inertial confinement fusion applications. Hugoniot calculations for liquid deuterium are given as another example of how the application of computationally efficient orbital-free density functional theory (OF-DFT) can be utilized with the employment of the developed functionals. Once the functionals have been tuned such that the OF-DFT Hugoniot calculation matches the Kohn–Sham solution at some low-temperature point, agreement with the reference Kohn–Sham results for the rest of the high temperature Hugoniot path is very good with relative errors for compression and pressure on the order of 2% or less.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Polynomial Scaling Localized Active Space Unitary Selective Coupled Cluster Singles and Doubles

We present a polynomial-scaling algorithm for the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. In this approach, cluster excitations are selected based on a threshold ϵ determined by the absolute gradients of the LAS-UCCSD energy with respect to cluster amplitudes. Using the generalized Wick’s theorem for multireference wave functions, we derive the gradient expression as a polynomial function of one-, two-, and three-body reduced density matrices and 1- and 2-electron integrals, valid for any multireference wave function. The resulting gradient implementation exhibits a memory scaling of 𝒪(N 6 ), with N spin orbitals in the combined active space of all fragments. The variational quantum eigensolver is used to optimize the selected cluster excitations on a quantum simulator. Furthermore, by plotting the energy error, defined as the difference between the LAS-USCCSD and corresponding CASCI energies, against the inverse cluster amplitude selection threshold (ϵ –1 ) for polyene chains containing 2 to 5 π-bond units, we establish a relationship between the energy error and the threshold. To further validate the accuracy of LAS-USCCSD, we computed the cis–trans isomerization energy of stilbene (a 20-qubit system) and the magnetic coupling constant of the tris-hydroxo-bridged chromium dimer [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ (evaluated as both 12- and 20-qubit systems) using the Qiskit-Qulacs simulator. Assessing such examples is important to determine the practical feasibility of quantum simulations for chemically realistic systems. Toward this goal, with the LAS-USCCSD algorithm we estimated the quantum resources required for simulating an active space of (30e,22o) in [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ , a size that remains beyond the reach of current quantum simulators for accurate treatment.

Algorithms↗
DOE OSTI2022

Construction of meta-GGA functionals through restoration of exact constraint adherence to regularized SCAN functionals

The strongly constrained and appropriately normed (SCAN) meta-GGA exchange–correlation functional is constructed as a chemical environment-determined interpolation between two separate energy densities: one describes single-orbital electron densities accurately and another describes slowly varying densities accurately. To conserve constraints known for the exact exchange–correlation functional, the derivatives of this interpolation vanish in the slowly varying limit. While theoretically convenient, this choice introduces numerical challenges that degrade the functional’s efficiency. We have recently reported a modification to the SCAN meta-GGA, termed restored-regularized-SCAN (r 2 SCAN) , that introduces two regularizations into SCAN, which improve its numerical performance at the expense of not recovering the fourth order term of the slowly varying density gradient expansion for exchange. Here, we show the derivation of a progression of density functional approximations [regularized SCAN (rSCAN), r++SCAN, r 2 SCAN, and r 4 SCAN] with increasing adherence to exact conditions while maintaining a smooth interpolation. The greater smoothness of r 2 SCAN seems to lead to better general accuracy than the additional exact constraint of SCAN or r 4 SCAN does.

, correlation energy↗
DOE OSTI2020

PyFLOSIC: Python-based Fermi–Löwdin orbital self-interaction correction

We present pyflosic, an open-source, general-purpose python implementation of the Fermi–Löwdin orbital self-interaction correction (FLO-SIC), which is based on the python simulation of chemistry framework (pyscf) electronic structure and quantum chemistry code. Thanks to pyscf, pyflosic can be used with any kind of Gaussian-type basis set, various kinds of radial and angular quadrature grids, and all exchange-correlation functionals within the local density approximation, generalized-gradient approximation (GGA), and meta-GGA provided in the libxc and xcfun libraries. A central aspect of FLO-SIC is the Fermi-orbital descriptors, which are used to estimate the self-interaction correction. Importantly, they can be initialized automatically within pyflosic; they can also be optimized within pyflosic with an interface to the atomic simulation environment, a python library that provides a variety of powerful gradient-based algorithms for geometry optimization. Although pyflosic has already facilitated applications of FLO-SIC to chemical studies, it offers an excellent starting point for further developments in FLO-SIC approaches, thanks to its use of a high-level programming language and pronounced modularity.

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

Engineering local strain for single-atom nuclear acoustic resonance in silicon

Mechanical strain plays a key role in the physics and operation of nanoscale semiconductor systems, including quantum dots and single-dopant devices. In this report we describe the design of a nanoelectronic device, where a single nuclear spin is coherently controlled via nuclear acoustic resonance (NAR) through the local application of dynamical strain. The strain drives spin transitions by modulating the nuclear quadrupole interaction. We adopt an AlN piezoelectric actuator compatible with standard silicon metal–oxide–semiconductor processing and optimize the device layout to maximize the NAR drive. We predict NAR Rabi frequencies of order 200 Hz for a single 123Sb nucleus in a wide region of the device. Spin transitions driven directly by electric fields are suppressed in the center of the device, allowing the observation of pure NAR. Using electric field gradient-elastic tensors calculated by the density-functional theory, we extend our predictions to other high-spin group-V donors in silicon and to the isoelectronic 73Ge atom.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Leptodermic corrections to the TOV equations and nuclear astrophysics within the effective surface approximation

The macroscopic model for a neutron star (NS) as a liquid drop at the equilibrium is used to extend the Tolman-Oppenheimer-Volkoff (TOV) equations taking into account the gradient terms responsible for the system surface. The parameters of the Schwarzschild metric in the spherical case are found with these surface corrections to the known leading (zero) order of the leptodermic approximation a/R << 1, where a is the NS effective-surface (ES) thickness, and R is the effective NS radius. The energy density $\mathscr{E}$ is considered in a general form including the functions of the particle number density and of its gradient terms. The macroscopic gravitational component $Φ$(ρ) of the energy density is taken into account in the simplest form as expansion in powers of $ρ$ – $\overline{ρ}$, where $\overline{ρ}$ is the saturation density, up to second order, in terms of its contributions to the separation particle energy and incompressibility. Density distributions ρ across the NS ES in the normal direction to the ES, which are derived in the simple analytical form at the same leading approximation, was used for the derivation of the modified TOV (MTOV) equations by accounting for their NS surface corrections. As a result, the MTOV equations are analytically solved at first order and the results are compared with the standard TOV approach of the zero order.

Magner, A. G. [Institute for Nuclear Research, Kyi↗