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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Monte Carlo Explicitly Correlated Second-Order Many-Body Green’s Function Calculations of Semiconductor Band Gaps

A systematically converging series of ab initio, post-density-functional, size-consistent, electron-correlated approximations is desired for predictive computing of felectronic band structures of insulating, semiconducting, and metallic solids. A series that meets all of these desiderata (except the applicability to metals) is ab initio many-body Green's function theory based on Gaussian-type-orbital (GTO) basis sets. Here, its leading-order approximation, the second-order Green's function (GF2) method in the diagonal and frequency-independent approximations with the aug-cc-pVDZ basis set, is applied to the fundamental band gaps of three semiconductors (diamond, silicon, and silicon carbide in the zincblende structure) using cluster models. Corrections are made to the basis-set-incompleteness errors by the explicit-correlation (F12) ansatz (GF2-F12) for the valence band edges. The crystals are modeled as surface-passivated clusters of increasing sizes, whose wave functions are expanded by up to 2709 GTO basis functions. Immense computational costs of these calculations are overcome by the highly scalable stochastic algorithm of the Monte Carlo GF2-F12 method, whose operation cost per state increases only as a cubic power of system size, which has a tiny memory footprint and easily achieves near-perfect parallel efficiency on thousands of CPUs or on hundreds of GPUs. The correlated, F12-corrected highest-occupied and lowest-unoccupied molecular-orbital energy (HOMO-LUMO) gap is 5.78 ± 0.07 eV for C 87 H 76 as compared with the experimental value of the fundamental (indirect) band gap of bulk diamond at 5.48 eV. The correlated, F12-corrected HOMO-LUMO gaps for Si 75 H 76 and Si 32 C 43 H 76 are 2.56 ± 0.15 eV and 3.50 ± 0.12 eV, respectively, which are expected to decrease further with increasing cluster sizes. As a result, the experimental fundamental (indirect) band gaps of bulk silicon and silicon carbide are 1.17 eV and 2.42 eV, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Impact of nonadiabatic charge transfer on the rate of redox chemistry of carbon oxides on rutile TiO2 (110) surface

We present the results of a density functional theory (DFT) within the LDA+U approximation on large models of partially reduced TiO2 (110) rutile surface to investigate the nature of charge transfer and the role of non-adiabatic effects on three prototypical redox reactions: (i) O2 adsorption (ii) CO oxidation and (iii) CO2 reduction. Charge-constrained DFT (cDFT) is used to estimate kinetic parameters for a Marcus theory rate law that accounts for adiabatic coupling effects on reaction rates. We find that for O2 adsorption, the coupling between adiabatic states is strong, leading to fast charge transfer rates. The lowest energy structures at high coverage consist of two chemisorbed O2-, one adsorbed at a VO site and the other adsorbed at an adjacent Ti5C site. For CO oxidation, however, all reactions are kinetically hindered on the ground state due to the weak adiabatic coupling at the state crossing, such that one has to overcome two kinetically unfavorable charge transfer events to drive the process (non-adiabatically) on the thermal ground state. The process can be driven by photochemical means but would result in an adsorbed radical [OCOO-] intermediate species. Similarly, CO2 reduction also proceeds via a non-adiabatic charge transfer to form an adsorbed CO2- species followed by a second non-adiabatic charge transfer to produce CO. Our analysis provides important computational guidelines for modeling these types of processes. We thank Z. Dohnalek, M. Hendersen, G. Kimmel, H. Metieu, and N. Petrik for invaluable discussions. This work was supported by the US Department of Energy, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences and performed at the Pacific Northwest National Laboratory (PNNL). PNNL is a multi-program national laboratory operated for DOE by Battelle. Computational resources were provided at W. R. Wiley Environmental Molecular Science Laboratory (EMSL), a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research located at PNNL and the National Energy Research Scientific Computing Center (NERSC) at Lawrence Berkeley National Laboratory.

charge transfer↗
DOE OSTI2026

Machine-learning guided search for phonon-mediated superconductivity in boron and carbon compounds

We present a workflow that iteratively combines ab-initio calculations with a machine-learning (ML) guided search for superconducting compounds with both dynamical stability and instability from imaginary phonon modes, the latter of which have been largely overlooked in previous studies. Electron-phonon coupling (EPC) properties and critical temperature (T c ) of 417 boron, carbon, and borocarbide compounds have been calculated with density functional perturbation theory (DFPT) and isotropic Eliashberg approximation. Our study addresses T c convergence of Brillouin zone sampling with an ansatz test, stabilizing imaginary phonon modes for significant EPC contributions, and comparing the performance of two ML models, especially when including compounds of dynamical instability. We predict a few promising superconducting compounds with formation energy just above the ground state convex hull, such as Ca 5 B 3 N 6 (35 K), TaNbC 2 (28.4 K), Nb 3 B 3 C (16.4 K), Y 2 B 3 C 2 (4.0 K), Pd 3 CaB (7.0 K), MoRuB 2 (15.6 K), RuVB 2 (15.0 K), RuSc 3 C 4 (6.6 K) among others.

Nepal, Niraj K. [Ames Laboratory (AMES), Ames, IA ↗
DOE OSTI2020

Implications of the fractional charge of hydroxide at the electrochemical interface

Rational design of materials that efficiently convert electrical energy into chemical bonds will ultimately depend on a thorough understanding of the electrochemical interface at the atomic level. Towards this goal, the use of density functional theory (DFT) at the generalized gradient approximation (GGA) level has been applied widely in the past 15 years. In the calculation of electrochemical reaction energetics using GGA-DFT, it is frequently implicitly assumed that ions in the Helmholtz plane have unit charge. However, the ion charge is observed to be fractional near the interface through both a capacitor model and through Bader charge partitioning. In this work, we show that this spurious charge transfer can be effectively mitigated by continuum charging of the electrolyte. Finally, we show that, similar to hydronium, the observed fractional charge of hydroxide is not due to a GGA level self-interaction error, as the partial charge is observed even when using hybrid level exchange–correlation functionals.

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

On the choice of reference orbitals for linear-response calculations of solution-phase K-edge X-ray absorption spectra

The simplest response theory methods for computing vertical excitation spectra in condensed-phase are configuration interaction with single excitations (CIS) and linear-response time-dependent density functional theory (TDDFT) within the Tamm–Dancoff approximation. In applications to X-ray absorption spectroscopy (XAS), methods like CIS and TDDFT that codify only single excitations into the wave function are prone to catastrophic errors in main-edge and post-edge features whose shapes act as a crucial fingerprint in structural analyses of liquids. Here, we show that these errors manifest primarily due to a lack of orbital relaxation in conventional linear-response theories and that core-ionized (n - 1-electron) references, like those of electron-affinity TDDFT, can eliminate the errors in the spectral profile, even in the highest-energy parts of the post-edge. Crucially, we find that single excitations atop core-ionized references are sufficient to elucidate liquid-phase XAS spectra with semi-quantitative accuracy, opening the door for methods like electron-affinity CIS/TDDFT to be used as efficient alternatives to higher-order wave function approaches.

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

Magnetic properties of MnBi based alloys: First-principles calculations for MnBi-Co and MnBi-Co-Fe cases

First-principles calculations of fundamental magnetic properties were performed for ordered MnBi, MnBi-Co, and MnBi-Co-Fe alloys to evaluate maximum energy product (BH)(max). Full potential linear-augmented plane wave (FLAPW) and linear-muffin-tin-orbital (LMTO) calculations using density functional theory (DFT) within the local spin density approximation (LSDA) were used and found to give a reasonable description of saturation magnetization (M-s), effective anisotropy constant (K-eff), and Curie temperature (T-c) for NiAs-structured MnBi crystal. We found that upon addition of Co, the M-s and K-eff increased, while T-c reduced. The magnetic anisotropy changed from weak anisotropy easy plane for MnBi to the strong easy axis anisotropy for MnBi-Co and MnBi-Co-Fe. (C) 2013 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported License.

Hong, YK↗
DOE OSTI2014

Assessing electronic structure approaches for gas-ligand interactions in metal-organic frameworks: The CO{sub 2}-benzene complex

Adsorption of gas molecules in metal-organic frameworks is governed by many factors, the most dominant of which are the interaction of the gas with open metal sites, and the interaction of the gas with the ligands. Herein, we examine the latter class of interaction in the context of CO{sub 2} binding to benzene. We begin by clarifying the geometry of the CO{sub 2}–benzene complex. We then generate a benchmark binding curve using a coupled-cluster approach with single, double, and perturbative triple excitations [CCSD(T)] at the complete basis set (CBS) limit. Against this ΔCCSD(T)/CBS standard, we evaluate a plethora of electronic structure approximations: Hartree-Fock, second-order Møller-Plesset perturbation theory (MP2) with the resolution-of-the-identity approximation, attenuated MP2, and a number of density functionals with and without different empirical and nonempirical van der Waals corrections. We find that finite-basis MP2 significantly overbinds the complex. On the other hand, even the simplest empirical correction to standard density functionals is sufficient to bring the binding energies to well within 1 kJ/mol of the benchmark, corresponding to an error of less than 10%; PBE-D in particular performs well. Methods that explicitly include nonlocal correlation kernels, such as VV10, vdW-DF2, and ωB97X-V, perform with similar accuracy for this system, as do ωB97X and M06-L.

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

Photoionization of the iodine 3d, 4s, and 4p orbitals in methyl iodide

Ionization of the I 3d, 4s, and 4p orbitals in methyl iodide (CH 3 I) has been studied by using synchrotron radiation to measure the total ion yield and by recording photoelectron spectra with linearly polarized radiation in two polarization orientations. The complete photoelectron spectrum of CH 3 I has been recorded at several photon energies, and bands due to the C 1s, I 3d, 4s, 4p, and 4d atomic-like orbitals, as well as the molecular orbitals, have been observed and assigned. In the vicinity of the I 3d 5/2 and 3d 3/2 ionization thresholds at 626.8 and 638.3 eV, respectively, the ion yield displays weak structure in the pre-edge region due to transitions into valence or Rydberg states, and, at higher energies, a shoulder and a broad maximum attributed to the I 3d 5/2 → εf and the I 3d 3/2 → εf shape resonances, respectively. The absorption spectrum calculated using time-dependent density functional theory, within the Tamm-Dancoff approximation, has allowed assignments to be proposed for the valence and Rydberg states. The Stieltjes imaging technique has been used to simulate the absorption spectrum above the ionization threshold and indicates that transitions into the f(l = 3) continuum channel dominate. This conclusion has been corroborated by a Continuum Multiple Scattering–Xα (CMS–Xα) calculation. The asymmetric broadening of the photoelectron bands associated with the I 3d orbital, due to post collision interaction, is taken into account in our experimental analysis. Experimentally derived photoelectron anisotropy parameters for the I 3d orbital are in good agreement with the theoretical predictions obtained with the CMS–Xα approach. The I 3d shake-up/shake-off photoelectron spectrum has been recorded, and assignments have been proposed for several of the satellites. The M 4 N 45 N 45 and M 5 N 45 N 45 Auger electron yields have been measured, and that for the M 5 N 45 N 45 decay exhibits a maximum due to interchannel coupling between the 3d5/2 and 3d3/2 continua. The photoelectron band associated with the I 4p orbital has an unusual appearance. Based upon previous theoretical work for the analogous Xe 4p orbital, it appears that the initial I 4p -1 hole state decays rapidly through Coster-Kronig and super-Coster-Kronig transitions. This leads to a redistribution of the spectral intensity associated with the I 4p orbital and results in a photoelectron spectrum containing a single structured band together with an extended continuum. Another continuum is observed on the high binding energy side of the peak due to the 4s orbital, and we assign this to super-Coster-Kronig transitions into the 4p -1 4d -1 continuum.

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

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing the accuracy. Advances in real space methods include techniques for accurate computation of bandgaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods, including GW and density functional based techniques. To provide an improved foundation for these calculations, we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2020

Photoabsorption, photoionization, and Auger processes at the carbon K edge in C H 3 I

The dynamics of photoabsorption, photoionization, and the associated Auger decay have been investigated at the carbon K edge in methyl iodide (CH 3 I) using linearly polarized synchrotron radiation. Ion yield measurements were used to investigate transitions in the pre-edge region due to excitations into either unoccupied valence or Rydberg states. The assignment of these transitions was achieved through comparison with theoretical x-ray absorption spectra calculated using time-dependent density functional theory, within the Tamm-Dancoff approximation. Several of the Rydberg states belonging to series converging onto the C1s ionization limit exhibit significant vibrational structure that is also interpreted using theoretical calculations. The C1s in CH 3 I photoelectron spectrum was measured, and the observed vibrational structure was assigned with the aid of theoretical predictions. Polarization dependent, resonantly excited, valence shell photoelectron spectra were recorded at photon energies coinciding with the C1s → σ*, C1s → 6sa1 and C1s → 6pe transitions in CH 3 I , thereby allowing photoelectron angular distributions to be determined. The nonresonantly excited C(KVV) Auger electron spectrum was measured and some of the features observed at high kinetic energies were attributed to transitions into valence orbitals possessing significant iodine character. The contributions of participator and spectator Auger decay to the resonantly excited photoelectron spectra have been assessed. The influence of participator decay appears minor whereas spectator decay results in the enhanced population of satellite states.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Vibronic States and Their Effect on the Temperature and Strain Dependence of Silicon-Vacancy Qubits in 4$\textit{H}$-SiC

Silicon-vacancy qubits in silicon carbide ( Si C ) are emerging tools in quantum-technology applications due to their excellent optical and spin properties. In this paper, we explore the effect of temperature and strain on these properties by focusing on the two silicon-vacancy qubits, V1 and V2, in 4$\textit{H}$-SiC. We apply density-functional theory beyond the Born-Oppenheimer approximation to describe the temperature-dependent mixing of electronic excited states assisted by phonons. We obtain a polaronic gap of around 5 and 22 meV for the V1 and V2 centers, respectively, which results in a significant difference in the temperature-dependent dephasing and zero-field splitting of the excited states, which explains recent experimental findings. Finally, we also compute how crystal deformations affect the zero-phonon line of these emitters. Our predictions are important ingredients in any quantum applications of these qubits sensitive to these effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Angle-Resolved Spectroscopy of Positronium Emission from a Cu(110) Surface

We report the affinity A Ps of positronium (Ps) atoms for a metal is the negative of the maximum kinetic energy with which Ps is emitted into vacuum when thermalized positrons in a metal encounter the surface. When this quantity is measured by ground state Ps time of flight (TOF), the precision is severely limited by the short triplet state lifetime of 142 ns. By quickly converting the emitted Ps atoms into long-lived Rydberg states, we are able to dramatically increase the TOF to allow precision measurements of A Ps . From our measurements made on a Cu(110) sample at T = 128 K, we find A Ps (128 K) = (-2.476±0.010 stat ±0.013 syst ) eV, compared with the result A Ps (128 K) = (-2.545±0.010 stat ±0.010 syst ) eV found using highly accurate generalized gradient approximations for both electrons and positrons within density functional theory. Such precision opens up opportunities in the quest for an improved density functional.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

First-principles study of magnetism and electric field effects in 2D systems

This review article provides a bird's-eye view of what first-principles based methods can contribute to next-generation device design and simulation. After a brief overview of methods and capabilities in the area, the authors focus on published work by their group since 2015 and current work on CrI3. The authors introduce both single- and dual-gate models in the framework of density functional theory and the constrained random phase approximation in estimating the Hubbard U for 2D systems vs their 3D counterparts. A wide range of systems, including graphene-based heterogeneous systems, transition metal dichalcogenides, and topological insulators, and a rich array of physical phenomena, including the macroscopic origin of polarization, field effects on magnetic order, interface state resonance induced peak in transmission coefficients, spin filtration, etc., are covered. For CrI3, the authors present their new results on bilayer systems such as the interplay between stacking and magnetic order, pressure dependence, and electric field induced magnetic phase transitions. The authors find that a bare bilayer CrI3, graphene|bilayer CrI3|graphene, hexagonal boron nitride (h-BN)|bilayer CrI3|h-BN, and h-BN|bilayer CrI3|graphene all have a different response at high field, while at small field, the difference is small except for graphene|bilayer CrI3|graphene. The authors conclude with discussion of some ongoing work and work planned in the near future, with the inclusion of further method development and applications.

Cheng, Hai-Ping (ORCID:0000000159901725)↗
DOE OSTI2020

Isoscalar and Isovector Giant Resonances in Closed Shells Nuclei and Bulk Properties of Nuclear Matter

Centroid energies, E CEN , of the isoscalar (T = 0) and isovector (T = 1) giant resonances of multipolarities L = 0 – 3 in 40,48 Ca, 68 Ni, 90 Zr, 116 Sn, 144 Sm and 208 Pb, were calculated within the fully self-consistent spherical Hartree-Fock (HF)-based random phase approximation (RPA) theory, using 33 different energy density functionals associated with Skyrme-type effective nucleon-nucleon interactions of the standard form commonly employed in the literature. Here, we also calculate the Pearson linear correlation coefficients between each E CEN and each bulk property of nuclear matter (NM), associated with the Skyrme interactions used in the calculations, and determine the sensitivity of E CEN to bulk properties of NM. By comparing the calculated values of E CEN to the experimental data, we constrain the values of the bulk NM properties. We find that interactions associated with the values of the NM effective mass, m* / m = 0.70 to 0.90, incompressibility coefficient, K NM = 210 to 240 MeV, and the enhancement coefficient of the energy weighted sum rule of the isovector giant dipole resonance, κ = 0.25 to 0.70, best reproduce the experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Multigrid based First-Principles Molecular Dynamics

MGmol ls a First-Principles Molecular Dynamics code. It relies on the Born-Oppenheimer approximation and models the electronic structure using Density Functional Theory, either LDA or PBE. Norm-conserving pseudopotentials are used to model atomic cores.

Fattebert, Jean-Luc↗
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↗
NASA NTRS1973

The muscle spindle as a feedback element in muscle control

The muscle spindle, the feedback element in the myotatic (stretch) reflex, is a major contributor to muscular control. Therefore, an accurate description of behavior of the muscle spindle during active contraction of the muscle, as well as during passive stretch, is essential to the understanding of muscle control. Animal experiments were performed in order to obtain the data necessary to model the muscle spindle. Spectral density functions were used to identify a linear approximation of the two types of nerve endings from the spindle. A model reference adaptive control system was used on a hybrid computer to optimize the anatomically defined lumped parameter estimate of the spindle. The derived nonlinear model accurately predicts the behavior of the muscle spindle both during active discharge and during its silent period. This model is used to determine the mechanism employed to control muscle movement.

Andrews, L. T.↗