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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Phase Stability, Band Gap Tuning, and Rashba Splitting in Selenium-Alloyed Bournonite: CuPbSb(S 1– x Se x ) 3

Recently, bournonite (CuPbSbS 3 ) has been identified as a potential ferroelectric photovoltaic (PV) material with a 1.3 eV band gap, which falls in an appropriate range for single junction PV devices. Progress in applying bournonite has yielded several studies reporting successful thin-film processing, the most recent of which demonstrated a 2.65% power conversion efficiency for a bournonite-based PV device. In an effort to explore bournonite band gap engineering for PV and other prospective applications (e.g., thermoelectric, spintronic), we here consider the solid-state synthesis of selenium-alloyed bournonite, CuPbSb(S 1-x Se x ) 3 , across the full range of x (0.0 ≤ x ≤ 1.0), and report phase purity for 0.0 ≤ x ≤ 0.5. We characterize the crystal structure and band gap of the samples using X-ray diffraction and diffuse reflectance spectroscopy, determining a band gap decrease for single-phase samples from 1.25 eV at x = 0.0 to 1.06 eV at x = 0.5. Formation energies determined using dispersion corrected hybrid density functional theory (DFT) support experimental findings and are consistent with the stability (instability) of the x = 0.5 (x = 1.0) structure relative to starting materials. The computed band structures show a decreasing trend in the band gap with increasing x, again consistent with experiment, with the states at the conduction (valence) band minima (maxima) being principally derived from Pb (S/Se) states. Spin texture analysis and detailed comparison of spin splitting parameters show that Se-alloying modifies the band gap while maintaining Rashba spin splitting character within the electronic structure. Rashba splitting is most pronounced for the conduction band minimum along the Γ-X path in the Brillouin zone. Energy separation between spin states is maintained at ΔE~0.2 eV for the various Se contents, with σz as the spin direction. Furthermore, a smaller spin splitting occurs along Γ-Z (decreasing value with increasing x), with σ x as the spin direction.

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

Computational Study of an Iron(II) Polypyridine Electrocatalyst for CO 2 Reduction: Key Roles for Intramolecular Interactions in CO 2 Binding and Proton Transfer

A solar-driven conversion of CO 2 into fuels by artificial photosynthesis would not only mitigate the greenhouse effect but also provide an alternative to obtain fuels in a renewable fashion. To this end, the new iron polypyridine catalyst [Fe(bpy NHEt PY2Me)L 2 ] 2+ (L = H 2 O, CH 3 CN) was recently developed for the electrochemical reduction of CO 2 to CO. Here, we performed density functional theory (DFT) electronic structure calculations to shed light on a possible pathway for CO 2 reduction and the origin of the selectivity between CO 2 reduction versus the hydrogen evolution reaction. The metal center remains Lewis acidic throughout the reduction process due to ligand loss and mainly ligand-based reduction stabilized by antiferromagnetic coupling to a high-spin Fe(II) center. This results in a high barrier for hydride formation but a facile addition and activation of CO 2 via an η 2 coordination and stabilizing hydrogen bonding by the amine group. The second unoccupied equatorial coordination site opens up the possibility for an intramolecular protonation with a coordinated water ligand. This facilitates protonation because not only CO 2 but also the proton source H 2 O is activated and properly aligned for a proton transfer due to the Fe-OH 2 bond; consequently, both protonation steps are facile. The moderate ligand field allows a rapid ligand exchange for a second intramolecular protonation step and facilitates an exergonic CO release. The lower selectivity of the related [Fe(bpy OH PY2Me)L 2 ] 2+ complex can be related to its more acidic second coordination sphere, which enables an intramolecular proton transfer that is kinetically competitive with CO 2 addition.

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

Metal–insulator transition in composition-tuned nickel oxide films

Abstract Thin films of the solid solution Nd 1 − x La x NiO 3 are grown in order to study the expected 0 K phase transitions at a specific composition. We experimentally map out the structural, electronic and magnetic properties as a function of x and a discontinuous, possibly first order, insulator–metal transition is observed at low temperature when x = 0.2. Raman spectroscopy and scanning transmission electron microscopy show that this is not associated with a correspondingly discontinuous global structural change. On the other hand, results from density functional theory (DFT) and combined DFT and dynamical mean field theory calculations produce a 0 K first order transition at around this composition. We further estimate the temperature-dependence of the transition from thermodynamic considerations and find that a discontinuous insulator–metal transition can be reproduced theoretically and implies a narrow insulator–metal phase coexistence with x . Finally, muon spin rotation ( µ SR) measurements suggest that there are non-static magnetic moments in the system that may be understood in the context of the first order nature of the 0 K transition and its associated phase coexistence regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Influence of temperature, doping, and amorphization on the electronic structure and magnetic damping of iron

Hybrid magnonic quantum systems have drawn increased attention in recent years for coherent quantum information processing, but too large magnetic damping is a persistent concern when metallic magnets are used. Their intrinsic damping is largely determined by electron-magnon scattering induced by spin-orbit interactions. In the low scattering limit, damping is dominated by intra-band electronic transitions, which has been theoretically shown to be proportional to the electronic density of states at the Fermi level. In this work, we focus on body-centered-cubic iron as a paradigmatic ferromagnetic material. We comprehensively study its electronic structure using first-principles density functional theory simulations and account for finite lattice temperature, boron (B) doping, and structure amorphization. Our results indicate that temperature induced atomic disorder and amorphous atomic geometries only have a minor influence. Instead, boron doping noticeably decreases the density of states near the Fermi level with an optimal doping level of 6.25%. In addition, we show that this reduction varies significantly for different atomic geometries and report that the highest reduction correlates with a large magnetization of the material. Furthermore, this may suggest materials growth under external magnetic fields as a route to explore in experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Spin Alignment of Vector Mesons in Heavy-Ion Collisions

Polarized quarks and antiquarks in high-energy heavy-ion collisions can lead to the spin alignment of vector mesons formed by quark coalescence. Here, using the relativistic spin Boltzmann equation for vector mesons derived from Kadanoff-Baym equations with an effective quark-meson model for strong interaction and quark coalescence model for hadronizaton, we calculate the spin density matrix element $ρ_{00}$ for Φ mesons and show that anisotropies of local field correlations with respect to the spin quantization direction lead to Φ meson's spin alignment. We propose that the local correlation or fluctuation of Φ fields is the dominant mechanism for the observed Φ meson's spin alignment and its strength can be extracted from experimental data as functions of collision energies. The calculated transverse momentum dependence of $ρ_{00}$ agrees with STAR's data. We further predict the azimuthal angle dependence of $ρ_{00}$ which can be tested in future experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Local strain-dependent electronic structure and perpendicular magnetic anisotropy of a MnGaN 2D magnetic monolayer

Local strain-dependent spin-polarized electronic structure of a two-dimensional (2D) magnetic layer is an exciting property for practical applications. For example, it holds the promise for advanced ultrathin spintronic nanodevices with customized electronic and magnetic properties by local strain engineering. Here, we demonstrate that the spin-polarized electronic structure of a 2D manganese gallium nitride (MnGaN-2D) magnetic monolayer is sensitive to intrinsic local lattice strain, as proven by first-principles calculations and indicated by scanning tunneling spectroscopy measurements. Atomic resolution images reveal a highly non-Gaussian lattice spacing/strain distribution, while the spectroscopy reveals variations in the electronic density of states. Simulations of the MnGaN-2D monolayer based on first-principles calculations, including both isotropic and anisotropic strains, confirm a highly strain-dependent manganese partial density of states. Spin-orbit coupling is included which indicates either out-of-plane perpendicular magnetic anisotropy (PMA) or in-plane magnetic anisotropy, dependent on the type of strain whether compressive or tensile, suggesting that MnGaN-2D is magnetoelastic. The MnGaN-2D PMA is further supported by superconducting quantum interference device magnetometry measurements which reveal a high spin polarization of ~ 79 % at room temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Band structures and Z 2 invariants of two-dimensional transition metal dichalcogenide monolayers from fully relativistic Dirac-Kohn-Sham theory using Gaussian-type orbitals

Two-dimensional (2D) materials exhibit a wide range of remarkable phenomena, many of which owe their existence to the relativistic spin-orbit coupling (SOC) effects. To understand and predict properties of materials containing heavy elements, such as the transition-metal dichalcogenides (TMDs), relativistic effects must be taken into account in first-principles calculations. We present an all-electron method based on the four-component Dirac Hamiltonian and Gaussian-type orbitals (GTOs) that overcomes complications associated with linear dependencies and ill-conditioned matrices that arise when diffuse functions are included in the basis. Until now, there has been no systematic study of the convergence of GTO basis sets for periodic solids either at the nonrelativistic or the relativistic level. Here, in this study, we provide such a study of relativistic band structures of the 2D TMDs in the hexagonal (2H), tetragonal (1T), and distorted tetragonal (1T') structures, along with a discussion of their SOC-driven properties (Rashba splitting and Z 2 topological invariants). We demonstrate the viability of our approach even when large basis sets with multiple basis functions involving various valence orbitals (denoted triple- and quadruple-ζ) are used in the relativistic regime. Our method does not require the use of pseudopotentials and provides access to all electronic states within the same framework. Our study paves the way for direct studies of material properties, such as the parameters in spin Hamiltonians, that depend heavily on the electron density near atomic nuclei where relativistic and SOC effects are the strongest.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2012

Particle Acceleration in Dissipative Pulsar Magnetospheres

Pulsar magnetospheres represent unipolar inductor-type electrical circuits at which an EM potential across the polar cap (due to the rotation of their magnetic field) drives currents that run in and out of the polar cap and close at infinity. An estimate ofthe magnitude of this current can be obtained by dividing the potential induced across the polar cap V approx = B(sub O) R(sub O)(Omega R(sub O)/c)(exp 2) by the impedance of free space Z approx eq 4 pi/c; the resulting polar cap current density is close to $n {GJ} c$ where $n_{GJ}$ is the Goldreich-Julian (GJ) charge density. This argument suggests that even at current densities close to the GJ one, pulsar magnetospheres have a significant component of electric field $E_{parallel}$, parallel to the magnetic field, a condition necessary for particle acceleration and the production of radiation. We present the magnetic and electric field structures as well as the currents, charge densities, spin down rates and potential drops along the magnetic field lines of pulsar magnetospheres which do not obey the ideal MHD condition $E cdot B = 0$. By relating the current density along the poloidal field lines to the parallel electric field via a kind of Ohm's law $J = sigma E_{parallel}$ we study the structure of these magnetospheres as a function of the conductivity $sigma$. We find that for $sigma gg OmegaS the solution tends to the (ideal) Force-Free one and to the Vacuum one for $sigma 11 OmegaS. Finally, we present dissipative magnetospheric solutions with spatially variable $sigma$ that supports various microphysical properties and are compatible with the observations.

Kazanas, Z.↗
DOE OSTI2026

Ferromagnetic Interlayer Exchange Coupling in a Few Layers of CrSBr on a Gold Thin Film

Abstract The two-dimensional character of van der Waals magnets allows for efficient control of their properties via proximity effects and electrical stimuli, making them promising candidates for application in spin electronics. We use spin-polarized low-energy electron microscopy to directly image the magnetic texture of thin CrSBr on top of a Au film, discovering a ferromagnetic ground state. We argue that the stabilization of the ferromagnetic ordering─as compared to the conventional antiferromagnetic one─is obtained via electron transfer from the Au film to the CrSBr flakes, in agreement with ab initio density functional theory calculations. Reflected-electron spectroscopy shows clear differences in the unoccupied density of states between a few layers of CrSBr on Au and bulk CrSBr, pointing toward electronic band structure modification in thin CrSBr. This work sheds light on the possibility of tuning the magnetic properties of two-dimensional magnets via substrate engineering.

Bosma, Rixt [University of Groningen , , ,]↗
DOE OSTI2021

Nuclei with Up to $$\varvec{A=6}$$ Nucleons with Artificial Neural Network Wave Functions

The ground-breaking works of Weinberg have opened the way to calculations of atomic nuclei that are based on systematically improvable Hamiltonians. Solving the associated many-body Schrödinger equation involves non-trivial difficulties, due to the non-perturbative nature and strong spin-isospin dependence of nuclear interactions. Artificial neural networks have proven to be able to compactly represent the wave functions of nuclei with up to $A=4$ nucleons. In this work, we extend this approach to $^6$Li and $^6$He nuclei, using as input a leading-order pionless effective field theory Hamiltonian. We successfully benchmark their binding energies, point-nucleon densities, and radii with the highly-accurate hyperspherical harmonics method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Homoleptic Perchlorophenyl “Ate” Complexes of Thorium(IV) and Uranium(IV)

The reaction of AnCl 4 (DME) n (An = Th, n = 2; U, n = 0) with 5 equiv of LiC 6 Cl 5 in Et 2 O resulted in the formation of homoleptic actinide-aryl “ate” complexes [Li(DME) 2 (Et 2 O)] 2 [Li(DME) 2 ][Th(C 6 Cl 5 ) 5 ] 3 ([Li][1]) and [Li(Et 2 O) 4 ][U(C 6 Cl 5 ) 5 ] ([Li][2]). Similarly, the reaction of AnCl 4 (DME)n (An = Th, n = 2; U, n = 0) with 3 equiv of LiC 6 Cl 5 in Et 2 O resulted in the formation of heteroleptic actinide-aryl “ate” complexes [Li(DME) 2 (Et 2 O)][Li(Et 2 O) 2 ][ThCl 3 (C 6 Cl 5 ) 3 ] ([Li][3]) and [Li(Et 2 O) 3 ][UCl 2 (C 6 Cl 5 ) 3 ] ([Li][4]). Density functional calculations show that the An–C ipso σ-bonds are considerably more covalent for the uranium complexes vs the thorium analogues, in line with past results. Additionally, good agreement between experiment and calculations is obtained for the 13 C ipso NMR chemical shifts in [Li][1] and [Li][3]. Here, the calculations demonstrate a deshielding by ca. 29 ppm from spin–orbit coupling effects originating at Th, which is a direct consequence of 5f orbital participation in the Th–C bonds.

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

Geometry, spin coupling, and dielectric control of redox potentials in [4Fe–4S] Clusters

Iron–sulfur (Fe–S) clusters are common biological cofactors that facilitate vital redox reactions. Despite extensive research, the molecular basis of redox potential tuning in ferredoxin-like proteins remains an active area of debate. In this study, we combine statistical analysis of over one thousand [4Fe–4S]-containing protein structures from the Protein Data Bank (PDB) with broken-symmetry and extended broken-symmetry density functional theory to examine how cysteine ligand orientations and environmental screening affect redox properties of the clusters. We identified five main ligand configurations, three of which are predominant in natural structures. Among these, the adiabatic electron affinity differs by less than 0.1 V, indicating that, while geometry plays a secondary role, it allows localized fine-tuning of redox properties. In contrast, electrostatic and solvation effects primarily determine the overall potential range.

Computational Chemistry↗
DOE OSTI2024

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

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

Inelastic neutron scattering study of the anisotropic S = 1 spin chain [ Ni ( HF 2 ) ( 3 – Clpyridine ) 4 ] BF 4

[Ni(HF 2 )(3–Clpyridine) 4 ]BF 4 (NBCT) is a one-dimensional, S = 1 spin chain material that shows no magnetic neutron Bragg peaks down to temperatures of 0.1 K. Previous work identified NBCT as being in the Haldane phase and near a quantum phase transition as a function of D/J to the large-D quantum paramagnet phase (QPM), where D is the axial single-ion anisotropy and J is the intrachain superexchange. Herein, inelastic neutron scattering results are presented on partially deuterated, 11B-enriched NBCT polycrystalline samples in zero magnetic field and down to temperatures of 0.3 K. Comparison to density matrix renormalization group calculations yields D/J = 1.51 and a significant rhombic single-ion anisotropy E (E/D ≈ 0.03, E/J ≈ 0.05). Here, these D, J, and E values place NBCT in the large-D QPM phase but precipitously near a quantum phase transition to a long-range ordered phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A benchmark dataset for Hydrogen Combustion

The generation of reference data for deep learning models is challenging for reactive systems, and more so for combustion reactions due to the extreme conditions that create radical species and alternative spin states during the combustion process. Here, we extend intrinsic reaction coordinate (IRC) calculations with ab initio MD simulations and normal mode displacement calculations to more extensively cover the potential energy surface for 19 reaction channels for hydrogen combustion. A total of ~290,000 potential energies and ~1,270,000 nuclear force vectors are evaluated with a high quality range-separated hybrid density functional, ωB97X-V, to construct the reference data set, including transition state ensembles, for the deep learning models to study hydrogen combustion reaction.

08 HYDROGEN↗
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)↗
NASA NTRS2020

Inverse Compton Cooling in the Coronae of Simulated Black Hole Accretion Flows

We present a formulation for a local cooling function to be employed in the diffuse, hot corona region of 3D GRMHD simulations of accreting black holes. This new cooling function calculates the cooling rate due to inverse Compton scattering by considering the relevant microphysics in each cell in the corona and approximating the radiation energy density and Compton temperature there by integrating over the thermal seed photon flux from the disk surface. The method either assumes ion and electron tem- peratures are equal (1T), or calculates them separately (2T) using an instantaneous equilibrium approach predicated on the actual relevant rate equations (Coulomb and Compton). The method is shown to be consistent with a more detailed ray-tracing calculation where the bulk of the cooling occurs, but is substantially less costly to per- form. As an example, we apply these methods to a harm3d simulation of a 10M⊙, non-spinning black hole, accreting at nominally 1% the Eddington value. Both 1T and 2T approaches lead to increased radiative efficiency and a larger fraction of total cooling in the corona as compared to the original target-temperature cooling function used by harm3d, especially in the 1T case. Time-averaged post-processing reveals that the continuum spectral observations predicted from these simulations are qualitatively similar to actual X-ray binary data, especially so for the 1T approach which yields a harder power-law component (Γ = 2.25) compared to the 2T version (Γ = 2.53).

Brooks E. Kinch↗
DOE OSTI2020

Inverse Compton Cooling in the Coronae of Simulated Black Hole Accretion Flows

Here we present a formulation for a local cooling function to be employed in the diffuse, hot corona region of 3D GRMHD simulations of accreting black holes. This new cooling function calculates the cooling rate due to inverse Compton scattering by considering the relevant microphysics in each cell in the corona and approximating the radiation energy density and Compton temperature thereby integrating over the thermal seed photon flux from the disk surface. The method either assumes the ion and electron temperatures are equal (1T) or calculates them separately (2T) using an instantaneous equilibrium approach predicated on the actual relevant rate equations (Coulomb and Compton). The method is shown to be consistent with a more detailed ray-tracing calculation where the bulk of the cooling occurs, but is substantially less costly to perform. As an example, we apply these methods to a harm3d simulation of a 10M⊙, non-spinning black hole, accreting at nominally 1% the Eddington value. This new approach leads to radiative efficiency values sime65% above Novikov–Thorne, with a larger fraction of total cooling in the corona as compared to simulations performed using the original target-temperature cooling function. Time-averaged post-processing reveals that the continuum spectral observations predicted from these simulations are qualitatively similar to actual X-ray binary data, especially so for the 1T approach, which yields a harder power-law component (Γ = 2.25) compared to the 2T version (Γ = 2.53).

79 ASTRONOMY AND ASTROPHYSICS↗