Search NASA⌕ Search

SEARCH · Search NASA

Results for “Classical density functional theory”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Phase Equilibrium of Water with Hexagonal and Cubic Ice Using the SCAN Functional

Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye toward studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ices Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Furthermore, our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic.

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

Density-related properties from self-interaction corrected density functional theory calculations

The Perdew–Zunger self-interaction correction (PZ-SIC) removes unphysical electron self-interaction from calculations employing standard density functional approximations. Doing so improves many computed properties, bringing them into better agreement with experimental observations or with results from high-level quantum chemistry calculations. However, while PZ-SIC generally corrects in the right direction relative to the corresponding reference values, in many cases, it over-corrects. For this reason, scaled-down versions of PZ-SIC have been proposed and investigated. These approaches have mostly employed exterior scaling in which SIC correction terms are scaled in the same way at every point in space. Recently, a new local, or interior, scaling SIC method was proposed on non-empirical grounds to restore a property of the exact, but unknown, density functional that is broken in PZ-SIC. Here, the scaling at each point depends on the character of the charge density at that point. However, the local scaling can be done in various ways while still restoring the behavior of the exact functional. In this work, we compare and contrast the performance of various interior scaling approaches for addressing over-corrections of calculated molecular dipole moments and atomic polarizabilities and properties that reflect the nature of the electronic charge density.

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

Photoelectron spectra of early 3 d -transition metal dioxide molecular anions from GW calculations

Photoelectron spectra of early 3 d -transition metal dioxide anions, Sc O 2 - , Ti O 2 - , V O 2 - , Cr O 2 - , and Mn O 2 - , are calculated using semilocal and hybrid density functional theory (DFT) and many-body perturbation theory within the GW approximation using one-shot perturbative and eigenvalue self-consistent formalisms. Different levels of theory are compared with each other and with available photoelectron spectra. We show that one-shot GW with a PBE0 starting point ( G 0 W 0 @PBE0) consistently provides very good agreement for all experimentally measured binding energies (within 0.1 eV–0.2 eV or less). We attribute this to the success of PBE0 in mitigating self-interaction error and providing good quasiparticle wave functions, which renders a first-order perturbative GW correction effective. One-shot GW calculations with a Perdew–Burke–Ernzerhof (PBE) starting point do poorly in predicting electron removal energies by underbinding orbitals with typical errors near 1.5 eV. A higher exact exchange amount of 50% in the DFT starting point of one-shot GW does not provide very good agreement with experiment by overbinding orbitals with typical errors near 0.5 eV. While not as accurate as G 0 W 0 @PBE0, the G -only eigenvalue self-consistent GW scheme with W fixed to the PBE level provides a reasonably predictive level of theory (typical errors near 0.3 eV) to describe photoelectron spectra of these 3 d -transition metal dioxide anions. Adding eigenvalue self-consistency also in W , on the other hand, worsens the agreement with experiment overall. Overall, our findings on the performance of various GW methods are discussed in the context of our previous studies on other transition metal oxide molecular systems.

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

6⁢𝐻-perovskite dimer lattice with antiferromagnetic interactions: Ba 3⁢ Zn 1−𝑥⁢ Ca 𝑥 ⁢Ru 2 ⁢O 9

Here, we investigate the magnetic behavior of the 6⁢𝐻-perovskite dimer lattice Ba 3 ⁢Zn 1−𝑥 ⁢Ca 𝑥 ⁢Ru 2 ⁢O 9 using analytical theory, density functional theory, inelastic neutron scattering, and modeling of historical magnetization and neutron-scattering data. A dimer mean-field theory built upon classical Luttinger–Tisza analysis generates a phase diagram revealing a transition from a nonmagnetic singlet to a finite-moment ground state as interdimer couplings increase. A (generalized) linear spin-wave theory captures multiplet mixing, excitation gap closing, and fluctuation-induced moment suppression. Density-functional-theory calculations on selected compounds, together with neutron spectroscopy of dilute Ba 3 ⁢Zn⁢(Ru 1−𝑥 ⁢Sb 𝑥 ) 2 ⁢O 9 , confirm the exchange hierarchy, enabling quantification of previously published experiments within this framework. Our results identify three mechanisms for magnetic moment suppression—quantum fluctuations, ligand hybridization, and nonmagnetic-singlet/magnetic-multiplet mixing.

Pajerowski, Daniel M. [Oak Ridge National Laborato↗
DOE OSTI2023

Room temperature colossal superparamagnetic order in aminoferrocene–graphene molecular magnets

Intensive studies are published for graphene-based molecular magnets due to their remarkable electric, thermal, and mechanical properties. However, to date, most of all produced molecular magnets are ligand based and subject to challenges regarding the stability of the ligand(s). The lack of long-range coupling limits high operating temperature and leads to a short-range magnetic order. Herein, we introduce an aminoferrocene-based graphene system with room temperature superparamagnetic behavior in the long-range magnetic order that exhibits colossal magnetocrystalline anisotropy of 8 × 10 5 and 3 × 10 7 J/m 3 in aminoferrocene and graphene-based aminoferrocene, respectively. These values are comparable to and even two orders of magnitude larger than pure iron metal. Aminoferrocene [C 10 H 11 FeN] + is synthesized by an electrophilic substitution reaction. It was then reacted with graphene oxide that was prepared by the modified Hammers method. The phase structure and functionalization of surface groups were characterized and confirmed by XRD, FT-IR, and Raman spectroscopy. To model the behavior of the aminoferrocene between two sheets of hydroxylated graphene, we have used density functional theory by placing the aminoferrocene molecule between two highly ordered hydroxylated sheets and allowing the structure to relax. Here, the strong bowing of the isolated graphene sheets suggests that the charge transfer and resulting magnetization could be strongly influenced by pressure effects. In contrast to strategies based on ligands surface attachment, our present work that uses interlayer intercalated aminoferrocene opens routes for future molecular magnets as well as the design of qubit arrays and quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Nonrad: Computing nonradiative capture coefficients from first principles

Point defects in semiconductor crystals provide a means for carriers to recombine nonradiatively. This recombination process impacts the performance of devices. We present the Nonrad code that implements the first-principles approach of Alkauskas et al. (2014) [8] for the evaluation of nonradiative capture coefficients based on a quantum-mechanical description of the capture process. An approach for evaluating electron-phonon coupling within the projector augmented wave formalism is presented. We also show that the common procedure of replacing Dirac delta functions with Gaussians can introduce errors into the resulting capture rate, and implement an alternative scheme to properly account for vibrational broadening. Finally, we assess the accuracy of using an analytic approximation to the Sommerfeld parameter by comparing with direct numerical evaluation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Improvements in the orbitalwise scaling down of Perdew–Zunger self-interaction correction in many-electron regions

The Perdew-Zunger (PZ) method offers a way to remove the self-interaction (SI) error from density functional approximations on an orbital by orbital basis. The PZ method provides significant improvements for the properties such as barrier heights or dissociation energies but results in over-correcting the properties well described by SI-uncorrected semi-local functional. One cure to rectify the over-correcting tendency is to scale down the magnitude of SI-correction of each orbital in the many electron region. We have implemented the orbitalwise scaled down SI-correction (OSIC) scheme of Vydrov et al. [J. Chem. Phys. 124, 094108 (2006)] using the Fermi-L owdin SI-correction method. After validating the OSIC implementation with previously reported OSICLSDA results, we examine its performance with the most successful non-empirical SCAN meta-GGA functional. Using different forms of scaling factors to identify one-electron regions, we assess the performance of OSIC-SCAN for a wide range of properties: total energies, ionization potentials and electron affinities for atoms, atomization energies, dissociation and reaction energies, and reaction barrier heights of molecules. Our results show that OSIC-SCAN provides superior results than the previously reported OSIC-LSDA, -PBE, and -TPSS results. Furthermore, we propose selective scaling of OSIC (SOSIC) to remove its major shortcoming that destroys the -1/r asymptotic behavior of the potentials. The SOSIC method gives the highest occupied orbital eigenvalues practically identical to those in PZSIC and unlike OSIC provides bound atomic anions even with larger powers of scaling factors. Furthermore, SOSIC compared to PZSIC or OSIC provides more balanced description of total energies and barrier heights.

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

Phonons of metallic hydrogen with quantum Monte Carlo

Here, we describe a simple scheme to perform phonon calculations with quantum Monte Carlo (QMC) methods and demonstrate it on metallic hydrogen. Because of the energy and length scales of metallic hydrogen and the statistical noise inherent to QMC methods, the conventional manner of calculating force constants is prohibitively expensive. We show that our alternate approach is nearly 100 times more efficient in resolving the force constants needed to calculate the phonon spectrum in the harmonic approximation. This requires only the calculation of atomic forces, as in the conventional approach, and otherwise little or no programmatic modification.

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

AlScO 3 perovskite—An ~8 eV bandgap oxide predicted to exhibit low small hole polaron ionization energies and p -type conductivity at elevated temperatures

Here, we investigate electronic structure and dopability of an ultrawide bandgap (UWBG) AlScO 3 perovskite, a known high-pressure and long-lived metastable oxide. From first-principles electronic structure calculations, HSE06(+G 0 W 0 ), we find this material to exhibit an indirect bandgap of around 8.0 eV. Defect calculations point to cation and oxygen vacancies as the dominant intrinsic point defects limiting extrinsic doping. While acceptor behaving Al and Sc vacancies prevent n-type doping, oxygen vacancies permit the Fermi energy to reach ~0.3 eV above the valence band maximum, rendering AlScO 3 p-type dopable. Furthermore, we find that both Mg and Zn could serve as extrinsic p-type dopants. Specifically, Mg is predicted to have achievable net acceptor concentrations of ~10 17 cm -3 with ionization energy of bound small hole polarons of ~0.49 eV and free ones below 0.1 eV. These values place AlScO 3 among the UWBG oxides with lowest bound small hole polaron ionization energies, which, as we find, is likely due to large ionic dielectric constant that correlates well with low hole polaron ionization energies across various UWBG oxides.

14 SOLAR ENERGY↗
DOE OSTI2023

Mechanistic origins of methyl-driven Overhauser DNP

The Overhauser effect in the dynamic nuclear polarization (DNP) of non-conducting solids has drawn much attention due to the potential for efficient high-field DNP as well as a general interest in the underlying principles that enable the Overhauser effect in small molecules. We recently reported the observation of 1 H and 2 H Overhauser effects in H 3 C- or D 3 C-functionalized Blatter radical analogs, which we presumed to be caused by methyl rotation. In this work, we look at the mechanism for methyl-driven Overhauser DNP in greater detail, considering methyl librations and tunneling in addition to classical rotation. Here, we predict the temperature dependence of these mechanisms using density functional theory and spin dynamics simulations. Comparisons with results from ultralow-temperature magic angle spinning-DNP experiments revealed that cross-relaxation at temperatures above 60 K originates from both libration and rotation, while librations dominate at lower temperatures. Due to the zero-point vibrational nature of these motions, they are not quenched by very low temperatures, and methyl-driven Overhauser DNP is expected to increase in efficiency down to 0 K, predominantly due to increases in nuclear relaxation times.

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

Dispersion interactions in proposed covalent superhydride superconductors

Recent developments in high temperature superconducting materials under high pressure have made numerical evaluation of the superconducting transition temperature (T c ) of predicted materials critically important as a means to identification. Existing methods of calculating T c often do not agree with each other or experiment, often due to the large number of complex factors that contribute to this property; among them is the neglect of dispersion interactions in commonly used density functionals. We evaluate the effect of including dispersion interactions on the predicted superconducting properties of two examples of the covalent superhydride class of very high-T c superconducting materials. In both cases, dispersion is found to have sizable effects, increasing the electron-phonon coupling as compared to the reference case of elemental niobium. A detailed investigation traces the origin of this effect in a 270 GPa $\mathrm {R3m CSH_7}$ structure to structural distortions driven by long-range electron-phonon interactions rather than novel bonding networks.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Substrate-Dependent Band Structures in Trilayer Graphene / h - BN Heterostructures

The tight-binding model has been spectacularly successful in elucidating the electronic and optical properties of a vast number of materials. Within the tight-binding model, the hopping parameters that determine much of the band structure are often taken as constants. Here, using ABA-stacked trilayer graphene as the model system, we show that, contrary to conventional wisdom, the hopping parameters and therefore band structures are not constants, but are systematically variable depending on their relative alignment angle between h-BN. Moreover, the addition or removal of the h-BN substrate results in an inversion of the K and K' valley in trilayer graphene’s lowest Landau level. Our work illustrates the oft-ignored and rather surprising impact of the substrates on band structures of 2D materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Spectral quadrature for the first principles study of crystal defects: Application to magnesium

In this work, we present an accurate and efficient finite-difference formulation and parallel implementation of Kohn-Sham Density (Operator) Functional Theory (DFT) for non periodic systems embedded in a bulk environment. Specifically, employing non-local pseudopotentials, local reformulation of electrostatics, and truncation of the spatial Kohn-Sham Hamiltonian, and the Linear Scaling Spectral Quadrature method to solve for the pointwise electronic fields in real-space and the non-local component of the atomic force, we develop a parallel finite difference framework suitable for distributed memory computing architectures to simulate non-periodic systems embedded in a bulk environment. Choosing examples from magnesium-aluminum alloys, we first demonstrate the convergence of energies and forces with respect to spectral quadrature polynomial order, and the width of the spatially truncated Hamiltonian. Next, we demonstrate the parallel scaling of our framework, and show that the computation time and memory scale linearly with respect to the number of atoms. Next, we use the developed framework to simulate isolated point defects and their interactions in magnesium-aluminum alloys. Our findings conclude that the binding energies of divacancies, Al solute-vacancy and two Al solute atoms are anisotropic and are dependent on cell size. Furthermore, the binding is favorable in all three cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Atomic Fingerprinting of Heteroatoms Using Noncontact Atomic Force Microscopy

We report immense strides have been made in increasing the resolution of scanning probe microscopy. Noncontact atomic force microscopy (n-AFM) now offers one the ability to characterize and visualize single molecules with subatomic resolution. Specifically, nc-AFM with a carbon monoxide (CO) functionalized tip has the ability to discriminate functional groups (–C≡C–, –CH 2 , –C$=$O, …), although discriminating atomic species often remains as an ongoing challenge. Here, real-space pseudopotentials constructed within density functional theory are employed to accurately simulate nc-AFM images of molecules containing heteroatoms (S, I, and N) within dibenzothiophene (DBT), 2-iodotriphenylene (ITP), acridine (ACR) and ferrous phthalocyanine (FePc). It is found that S and I atoms can be easily identified from C based on their unique features. For N atoms, a use of tip functionalization is proposed to effectively discriminate them from C atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Comparing the performance of density functionals in describing the adsorption of atoms and small molecules on Ni(111)

Reaction networks calculated using density functional theory (DFT) play a crucial role in understanding and predicting heterogeneous catalyst performance. However, results obtained using different density functionals can vary significantly. Accordingly, it is important to assess the accuracy of available functionals in capturing adsorption energetics derived from single-crystal adsorption microcalorimetry experiments. Here, we evaluate how different density functionals describe the molecular and dissociative adsorption of CH 3 I, the adsorption of CH 3 , I, and H, and the dissociative adsorption of CH 4 on Ni(111). First, we identify energetically-preferred adsorbed states for each adsorbate at three surface coverages and calculate the corresponding enthalpies of adsorption at 160 K. Then, we compare the results to experimental measurements published in the literature. For each adsorption process we find at least one functional that is quantitatively accurate within experimental error, but no functional is accurate (within experimental error) for all adsorption processes studied. Even when assuming an additional DFT error of ±20 kJ/mol on top of the experimental error, only the PBE-D3 and RPBE-D3 functionals are accurate for all considered adsorption systems. Overall, these results indicate that quantitative agreement between density functional theory calculations and experimental measurements is both system- and functional-dependent. Importantly, single-crystal adsorption microcalorimetry experiments will continue to play an important role for benchmarking density functionals.

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

Comparative study of Minnesota functionals performance on ferroelectric BaTiO 3 and PbTiO 3

Density functional theory based simulations are the leading tool for the computational investigation of ferroelectrics and the parametrization of their classical potentials. However, the predictions often depend strongly on the exchange correlation functional. The most popular choices, LDA and GGA, tend to underestimate or overestimate some structural, electric, and energy properties. These impede development of highly accurate classical potentials that extend the reach of first-principles simulations to finite temperatures and realistic sizes. In this work, we investigate the performance of recently developed Minnesota exchange correlation functionals on prototypical ferroelectrics BaTiO 3 , and PbTiO 3 , in comparison with some popular ones. We find that there exists strong correlation between predictions for some properties (tetragonality, phase energy difference, and polarization) by different functionals. Along the correlation line, we find a range of functionals (including some from the Minnesota suite) whose predictions fall in between those of LDA and GGA and, therefore, offer an improvement. A way to relatively rank functional performance with respect to chosen benchmarks is proposed and applied to identify the top performers for BaTiO 3 and PbTiO 3 . Here, the performance is found to be material dependent. Therefore, we propose that the performance assessment carried out in this work is employed to "screen" functionals prior to their use on ferroelectrics. Quick and computationally inexpensive, it is likely to lead to improved descriptions, especially for classical potential parametrization. The Minnesota functional suite is found to be suitable for this task.

36 MATERIALS SCIENCE↗