Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spin density functionals”

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 37 records · Page 2

Spin-state gaps and self-interaction-corrected density functional approximations: Octahedral Fe(II) complexes as case study

Accurate prediction of a spin-state energy difference is crucial for understanding the spin crossover phenomena and is very challenging for density functional approximations, especially for local and semi-local approximations due to delocalization errors. Here, we investigate the effect of the self-interaction error removal from the local spin density approximation (LSDA) and Perdew–Burke–Ernzerhof generalized gradient approximation on the spin-state gaps of Fe(II) complexes with various ligands using recently developed locally scaled self-interaction correction (LSIC) by Zope et al. [J. Chem. Phys. 151, 214108 (2019)]. The LSIC method is exact for one-electron density, recovers the uniform electron gas limit of the underlying functional, and approaches the well-known Perdew–Zunger self-interaction correction (PZSIC) as a particular case when the scaling factor is set to unity. Our results, when compared with reference diffusion Monte Carlo results, show that the PZSIC method significantly overestimates spin-state gaps favoring low spin states for all ligands and does not improve upon density functional approximations. The perturbative LSIC-LSDA using PZSIC densities significantly improves the gaps with a mean absolute error of 0.51 eV but slightly overcorrects for the stronger CO ligands. Finally, the quasi-self-consistent LSIC-LSDA, such as coupled-cluster single double and perturbative triple [CCSD(T)], gives a correct sign of spin-state gaps for all ligands with a mean absolute error of 0.56 eV, comparable to that of CCSD(T) (0.49 eV).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Small bimetallic clusters Agn-1M (M = Au, Co, Cu, Ni, Pd, Pt; n = 3, 9, 15): Density functional theory and genetic algorithm

We investigated the effect of size and composition on the properties of bimetallic nanoclusters. The geometric structures, stabilities, and electronic properties of size-selected Ag n-1 M (M = Au, Co, Cu, Ni, Pd, Pt; n = 3, 9, 15) bimetallic nanoclusters are systematically analyzed using spin-polarized density functional theory (DFT) within the generalized gradient approximation (GGA). We determine the most stable geometries for these clusters using a genetic algorithm (GA) in combination with DFT. Our results show that doping pure silver clusters with an M atom (transition metal), referred to as a “guest atom”, increases the stability as compared to pure Ag n (n = 3, 9, 15) clusters. The results for various properties including formation energy per atom, electronic structure, magnetic moments, and vibrational density of states (VDOS) are evaluated as a function of both size and composition of the system. The adsorption of selected bimetallic clusters on hydroxylated alumina substrate shows weak binding and minor changes in geometric properties except for Ag 8 Pt.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Adsorption-Induced Surface Magnetism

We report the emergence of adsorption-induced magnetism from heterohelicene molecules on a nonmagnetic Cu(100) surface. Spin-polarized low-energy electron microscopy measurements reveal spin-dependent electron reflectivity for enantiopure 7,12,17-trioxa[11]helicene (TO[11]H) monolayers, indicating the formation of a spin-polarized state localized in the topmost copper layer. Control experiments on clean Cu(100) and TO[11]H on highly oriented pyrolytic graphite show no such effect, excluding artifacts and chirality-induced spin selectivity as origins. Spin-polarized density functional theory calculations with hybrid functionals attribute the magnetism to strong chemisorption, which induces hybridization between the molecular HOMO and copper s- and d-states, driving asymmetric spin-polarized charge redistribution at the interface. An extended Newns-Anderson-Grimley model incorporating on-site Coulomb repulsion in Cu d-orbitals reproduces the emergence of interfacial spin polarization above a threshold interaction strength, highlighting the key roles of hybridization parameters and Coulomb correlation. These findings reveal a mechanism for inducing magnetism at molecule-metal interfaces without inherently magnetic components, offering avenues for engineering spin-polarized states in organic-inorganic hybrid systems.

Hubbard model↗

Tunneling magnetoresistance in MgO tunnel junctions with Fe-based leads in empirically corrected density functional theory

The minority-spin Fe/MgO interface states are at the Fermi level in density functional theory (DFT), but experimental evidence and GW calculations place them slightly higher in energy. This small shift can strongly influence tunneling magnetoresistance (TMR) in junctions with a thin MgO barrier and its dependence on the concentration of Co in the electrodes. Here, in this study, an empirical potential correction to DFT is introduced to shift the interface states up to match the tunnel spectroscopy data. With this shift, TMR in Fe/MgO/Fe junctions exceeds 800% and 3000% at 3 and 4 monolayers (ML) of MgO, respectively. We further consider the effect of alloying of the Fe electrodes with up to 30% Co or 10% V, treating them in the coherent potential approximation (CPA). Alloying with Co broadens the interface states and brings a large incoherent minority-spin spectral weight to the Fermi level. Alloying with V brings the minority-spin resonant states close to the Fermi level. However, in both cases the minority-spin spectral weight at the Fermi level resides primarily at the periphery of the Brillouin zone, which is favorable for spin filtering. Using convolutions of k II -resolved barrier densities of states calculated in CPA, it is found that TMR is strongly reduced by alloying with Co or V but still remains above 500% at 4 ML of MgO up to 30% of Co or 5% V. At 5 ML, the TMR increases above 1000% in all systems considered. However, while TMR declines sharply with increasing bias up to 0.2 eV in the tunnel junctions with pure Fe leads, it remains almost constant up to 0.5 eV if leads are alloyed with Co.

74 ATOMIC AND MOLECULAR PHYSICS↗

Nonconventional NMR Spin-Coupling Constants in Oligosaccharide Conformational Modeling: Structural Dependencies Determined from Density Functional Theory Calculations

Nonconventional NMR spin-coupling constants were investigated to determine their potential as conformational constraints in MA’AT modeling of the O-glycosidic linkages of oligosaccharides. Four ( 1 J C1',H1' , 1 J C1',C2' , 2 J C1',H2' , and 2 J C2',H1' ) and eight ( 1 J C4,H4 , 1 J C3,C4 , 1 J C4,C5 , 2 J C3,H4 , 2 J C4,H3 , 2 J C5,H4 , 2 J C4,H5 , and 2 J C3,C5 ) spin-couplings in methyl β-d-galactopyranosyl-(1→4)-β-d-glucopyranoside (methyl β-lactoside) were calculated using density functional theory (DFT) to determine their dependencies on O-glycosidic linkage C–O torsion angles, $\phi$ and ψ, respectively. Long-range 4 J H1',H4 was also examined as a potential conformational constraint of either $\phi$ or ψ. Secondary effects of exocyclic (hydroxyl) C–O bond rotation within or proximal to these coupling pathways were investigated. Based on the findings of methyl β-lactoside, analogous J-couplings were studied in five additional two-bond O-glycosidic linkages [βGlcNAc-(1→4)-βMan, 2-deoxy-βGlc-(1→4)-βGlc, αMan-(1→3)-βMan, αMan-(1→2)-αMan, and βGlcNAc(1→2)-αMan] to determine whether the coupling behaviors observed in methyl β-lactoside were more broadly observed. Of the 13 nonconventional J-couplings studied, 7 exhibit properties that may be useful in future MA’AT modeling of O-glycosidic linkages, none of which involve coupling pathways that include the linkage C–O bonds. The findings also provide new insights into the general effects of exocyclic C–O bond conformation on the magnitude of experimental spin-couplings in saccharides and other hydroxyl-containing molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Layer-dependent spin-resolved electronic structure of ferromagnetic triple-layered ruthenate Sr4Ru3O10

High-resolution angle- and spin-resolved photoemission spectroscopy (ARPES) of the triple-layered ruthenate Sr4Ru3O10 reveals features of the electronic structure that extend our understanding of the layered strontium ruthenates. The spectra near the Fermi energy are very different from the nonmagnetic analogues Sr2RuO4 and Sr3Ru2O7 with distinct Fermi surfaces for wide electronlike minority spin bands around the zone center and narrow holelike majority spin Fermi surface contours around the zone corners. The most dramatic results are two narrow spectral peaks ∼30 meV below the Fermi level, a spin-minority holelike band at the Brillouin zone center, and a spin-majority saddle-band van Hove singularity at the zone edge, which exhibits almost 100% spin polarization at low temperature, and a strong temperature dependent coherence-incoherence crossover attributed to Hund metal correlations. Quantitative comparison of the ARPES to spin-polarized density functional theory (DFT) calculations identify the specific antibonding and nonbonding orbital origins of the narrow bands, with a prediction of different spatial localization in the central and outer layers. This is shown to be consistent with experimental ARPES multizone matrix element intensity variations, and implicates outer-layer-specific control of the in-plane metamagnetism. The renormalization of the bands relative to the mean-field DFT, the demonstration of spin-polarized oxygen bands, and of spin-minority and spin-majority band-crossing hybridization provide a more complete picture of the magnetism which displays aspects of both delocalized and local moment behavior.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Noncollinear ground states of solids with a source-free exchange correlation functional

In this paper, we expand upon the source-free (SF) exchange correlation (XC) functional developed by Sangeeta Sharma and coworkers to plane-wave density functional theory (DFT) based on the projector augmented wave (PAW) method. This constraint is implemented by the current authors within the VASP source code, using a fast Poisson solver that capitalizes on the parallel three-dimensional fast Fourier transforms (FFTs) implemented in VASP. Using this modified XC functional, we explore the improved convergence behavior that results from applying this constraint to the GGA-PBE+U+J functional. In the process, we compare the noncollinear magnetic ground state computed by each functional and their SF counterpart for a select number of magnetic materials in order to provide a metric for comparing with experimentally determined magnetic orderings. We observe significantly improved agreement with experimentally measured magnetic ground-state structures after applying the source-free constraint. Furthermore, we explore the importance of considering probability current densities in spin-polarized systems, even under no applied field. We analyze the XC torque as well, in order to provide theoretical and computational analyses of the net XC magnetic torque induced by the source-free constraint. Along these lines, we highlight the importance of properly considering the real-space integral of the source-free local magnetic XC field. Our analyses on probability currents, net torque, and constant terms draw additional links to the rich body of previous research on spin-current density functional theory (SCDFT), and pave the way for future extensions and corrections to the SF corrected XC functional.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Enhancing Nitrogen Activation in Electrochemical Reduction: The Role of Rare Earth Oxide Surface Configurations

The pursuit of sustainable ammonia synthesis has prompted the exploration of ambient electrochemical nitrogen reduction reaction (e-NRR) as an alternative to the energy-intensive Haber-Bosch process. Here, this study conducted a theoretical investigation into the use of rare earth oxide materials, specifically dysprosium oxide (Dy 2 O 3 ), as potential electrocatalysts for NRR. Utilizing spin-polarized density functional theory calculations, we explored the interaction between Dy 2 O 3 surfaces and nitrogen (N 2 ) molecules, examining the capability of Dy 2 O 3 to adsorb and activate N 2 under ambient conditions. The results indicate that Dy 2 O 3 surfaces exhibit diverse configurations and bonding environments, providing a variety of reactive sites that display different behaviors in N 2 adsorption and activation. The distinctive electronic structure and surface chemistry of a particular Dy 2 O 3 surface configuration were found to significantly enhance the activation of N 2 by promoting charge transfer, which facilitates the NRR process. This research provides deep insights into the mechanistic pathways of N 2 reduction over Dy 2 O 3 , highlighting the surface properties as pivotal in catalysis. These theoretical insights serve as a foundation for the development of novel rare earth-based electrocatalytic materials for efficient ambient e-NRR, potentially transforming ammonia production into a greener and more energy-efficient process.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

First-Principles Studies of Tritium Species Diffusivity Across the Interface of Nickel-Plated Zircaloy-4

In this study, we will calculate the 3H2 and 3H2O diffusion mechanism in the Ni coated layer and the dominant 3H species dissolution into the Zircaloy-4 getter. We will optimize the Ni/Zr interface and introduce Sn and O impurities. To achieve the objectives, we will perform the investigation by using spin-polarized density functional theory (DFT) methods to capture the effects of Ni magnetic properties on 3H diffusion. We will introduce NiO/Ni(OH) on the Ni side closer to the interface. Then, the 3H diffusion pathways and barriers from Ni side to the zircalory-4 will be explored. We will investigate how NiO/Ni(OH) affects the 3H diffusion across the interface.

Paudel, Hari P.↗

Nonprecious Single Atom Catalyst for Methane Pyrolysis

The development of a suitable catalytic system for methane pyrolysis reactions requires a detailed investigation of the activation energy of C-H bonds on catalysts, as well as their stability against sintering and coke formation. In this work, both single-metal Ni atoms and small clusters of Ni atoms deposited on titanium nitride (TiN) plasmonic nanoparticles were characterized for the C-H bond activation of a methane pyrolysis reaction using ab initio spin-polarized density functional theory (DFT) calculations. The present work shows the complete reaction pathway, including energy barriers for C-H bond activation and dehydrogenated fragments, during the methane pyrolysis reaction on catalytic systems. Interestingly, the C-H bond activation barriers were low for both Ni single-atom and Ni-clusters, showing the energy barriers of ~1.10 eV and ~0.88 eV, respectively. Additionally, single-atom Ni-TiN showed weaker binding to adsorbates, and a net endothermic reaction pathway indicated that the single-atom Ni-TiN was expected to resist coke formation on its surface. However, these Ni single-atom catalysts can sinter, aggregate into a small cluster, and form a coke layer from the highly exothermic reaction pathway that the cluster takes despite the facile reaction pathway.

08 HYDROGEN↗

Time-dependent density functional theory of charge, energy, and spin transport, and quantum continuum mechanics (Final Report)

The overarching goal of this project, which unfolded over a period of seventeen years, including two no-cost extensions in the last two years, was to extend and deepen our ability to use time-dependent density functional theory and related theoretical methods for the description of the transport of charge, spin, or energy in quantum electronic systems at the nanoscale. The output of our efforts is contained in 52 papers, listed below, which we have published since 2006. There have also been many presentations at international conferences, workshops, colloquia and seminars. This body of work has been carried out with the essential input of 8 postdoctoral students (who were supported by the grant) and several external senior collaborators and their postdocs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Generating Bright Emissive States by Modulating the Bandgap of Monolayer Tungsten Diselenide

Transition metal dichalcogenides (TMDs) are essential due to their fascinating electronic and optical properties, strong exciton binding energy, and layer-dependent bandgap. They can be tuned to function as a single-photon emitter, but the quantum yield of photoluminescence of single-layer WSe 2 is low. There is some evidence that a bright emissive state can be produced by introducing local defects through functionalization. In this paper we use spin-polarized periodic density functional theory (DFT) to study the effect of functionalization for the specific case of cyclic carbenes. We find that the simplest of these molecules, cyclopentadiene carbene (Cyc), binds to the surface by a covalent bond leading to a substantial change in the bandgap (1.24 eV compared to 1.64 eV for the pristine surface). Moreover, there are semiflat bands below the Fermi level that originate from σ bonding and the interaction between the lone-pair p orbitals on the carbon of Cyc that bonds with Se on the TMD. Cyc is found to form a Type IIa heterojunction before and after contact with the surface, where the “a” refers to ordering of the energy levels where the TMD levels have the larger splitting and “b” to the reverse. The effect of various electron-withdrawing and electron-donating groups on Cyc is investigated, and it is found that the direct bandgap and heterointerfaces can be chemically tuned with covalently bound functional groups. With an electron-withdrawing group (EWG) such as –CHO, –COCl, and –CN attached to Cyc, a Type III//Type IIa interface is formed, whereas all other EWGs used in this study form a Type IIa junction, before and after contact to the surface. The electron-donating groups (EDG) form Type Ib//Type Ia junctions, and in particular, we find that the system WSe 2 + Cyc-Me band structure consists of semiflat bands at valence band maxima localized on Cyc-Me and the conduction band minimum is coupled between the surface and functional group with a direct bandgap of 0.88 eV. Hence, we predict that the Cyc-Me-functionalized monolayer WSe 2 will produce a red-shifted bright emissive state. Furthermore, we also find that carbenes with the Cyc all-carbon ring have (1) a triplet ground state and (2) form covalent bonds to WSe 2 , while heterocyclic carbenes (1) have a singlet ground state and (2) do not form covalent bonds to WSe 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Excited State Properties of Point Defects in Semiconductors and Insulators Investigated with Time-Dependent Density Functional Theory

Here, we present a formulation of spin-conserving and spin-flip hybrid time-dependent density functional theory (TDDFT), including the calculation of analytical forces, which allows for efficient calculations of excited state properties of solid-state systems with hundreds to thousands of atoms. We discuss an implementation on both GPU- and CPU-based architectures along with several acceleration techniques. We then apply our formulation to the study of several point defects in semiconductors and insulators, specifically the negatively charged nitrogen-vacancy and neutral silicon-vacancy centers in diamond, the neutral divacancy center in 4H silicon carbide, and the neutral oxygen-vacancy center in magnesium oxide. Our results highlight the importance of taking into account structural relaxations in excited states in order to interpret and predict optical absorption and emission mechanisms in spin defects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dzyaloshinskii-Moriya interaction induced magnetoelectric coupling in a tetrahedral molecular spin-frustrated system

We have investigated magnetoelectric (ME) coupling in the single-molecule magnet Mn 4 Te 4 (PEt 3 ) 4 with tetrahedral spin frustration. Our density functional studies found that an electric dipole moment can emerge with various noncollinear spin orderings. The forms of spin-dependent dipole are determined and consistent with that in noncentrosymmetric magnets driven by the Dzyaloshinskii-Moriya interaction. Writing a parameterized spin Hamiltonian, after solving for eigenvalues and eigenstates, we quantified the ME coupling by calculating the thermal average of the electric and magnetic susceptibilities, which can be influenced by external magnetic and electric fields, respectively. The quadratic relations are expected to be observable in experiments.

36 MATERIALS SCIENCE↗

Study of Self-Interaction Errors in Density Functional Calculations of Magnetic Exchange Coupling Constants Using Three Self-Interaction Correction Methods

We examine the role of self-interaction errors (SIE) removal on the evaluation of magnetic exchange coupling constants. In particular we analyze the effect of scaling down the self-interaction-correction (SIC) for three non-empirical density functional approximations (DFAs) namely, the local spin density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation, and recent SCAN family of meta-GGA functionals. To this end, we employ three one-electron SIC methods: Perdew-Zunger [Perdew, J. P.; Zunger, A. Phys. Rev. B, 1981, 23, 5048] SIC, the orbitalwise scaled SIC method [Vydrov, O. A. et al., J. Chem. Phys. 2006, 124, 094108], and the recent local scaling method [Zope, R. R. et al., J. Chem. Phys. 2019, 151, 214108]. We compute the magnetic exchange coupling constants using the spin projection and non projection approaches for sets of molecules composed of dinuclear and polynuclear H-He models, organic radical molecules, and chlorocuprate, and compare these results against accurate theories and experiment. Our results show that for the systems that mainly consist of single electron regions, PZSIC performs well but for more complex organic systems and the chlorcuprates, an overcorrecting tendency of PZSIC combined with the DFAs utilized in this work is more pronounced, and in such cases LSIC with kinetic energy density ratio performs better than PZSIC. Furthermore, analysis of the results in terms of SIC corrections to the density and to the total energy shows that both density and energy correction are required to obtain an improved prediction of magnetic exchange couplings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Real-space visualization of short-range antiferromagnetic correlations in a magnetically enhanced thermoelectric

Short-range magnetic correlations can significantly increase the thermopower of magnetic semiconductors, representing a noteworthy development in the decades-long effort to develop high-performance thermoelectric materials. Here, we reveal the nature of the thermopower-enhancing magnetic correlations in the antiferromagnetic semiconductor MnTe. Using magnetic pair distribution function analysis of neutron scattering data, we obtain a detailed, real-space view of robust, nanometer-scale, antiferromagnetic correlations that persist into the paramagnetic phase above the Neel temperature $T_N$ = 307 K. In this work, the magnetic correlation length in the paramagnetic state is significantly longer along the crystallographic c axis than within the ab plane, pointing to anisotropic magnetic interactions. Ab initio calculations of the spin-spin correlations using density functional theory in the disordered local moment approach reproduce this result with quantitative accuracy. These findings constitute the first real-space picture of short-range spin correlations in a magnetically enhanced thermoelectric and inform future efforts to optimize thermoelectric performance by magnetic means.

36 MATERIALS SCIENCE↗

Ab initio prediction of rapid kinetics of Fe impurities in δ -Pu

Here, we study the formation energies of iron impurities in δ-Pu within spin–orbital-polarized density functional theory (SOP-DFT). The thermodynamic solubility limit of iron in δ-Pu is calculated, indicating low miscibility. We show that surprisingly, Fe impurities at equilibrium are almost equally likely to occupy octahedral interstitial sites or substitutional sites, with slight preference for the former. In contrast, we find the energy of the tetrahedral interstitial Fe to be nearly 1 eV higher than the octahedral one. We explore the energy landscape for Fe impurity hopping diffusion and conclude that Fe impurities in δ-Pu are divided into two populations: (i) Immobile substitutional Fe impurities and (ii) highly mobile interstitial Fe impurities. The latter, (ii), migrate between octahedral interstitial sites with an energy barrier of around 0.2 eV. The energy barrier for exchange between the two populations is calculated to exceed 0.7 eV. Finally, we discuss the role of magnetic order on the impurity energetics.

Atomic structure↗