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 55 records · Page 3

Extending Orbital-Optimized Density Functional Theory to L-Edge XPS and Beyond: Spin–Orbit Coupling via Nonorthogonal Quasi-Degenerate Perturbation Theory

Quantum mechanical calculations of core electron binding energies (CEBEs) are relevant to interpreting X-ray photoelectron spectroscopy (XPS). Orbital-optimized density functional theory (OO-DFT) accurately predicts K-edge CEBEs but is challenged by the presence of significant spin–orbit coupling (SOC) at L- and higher edges involving inner-shell orbitals with nonzero angular momentum. Here, to extend OO-DFT to L-edges and higher, our method utilizes scalar-relativistic, spin-restricted open-shell OO-DFT to construct a minimal, quasi-degenerate basis of core-hole states corresponding to a chosen inner-shell (e.g., ionizing all six possible 2p spin orbitals). Nonorthogonal configuration interaction (NOCI) is then used to obtain the matrix elements of the full Hamiltonian including SOC in this quasi-degenerate model space of determinants. Using a screened 1-electron SOC operator parametrized with the Dirac-Coulomb-Breit (DCB) Hamiltonian results in doublet splitting (DS) values for third row elements that are nearly in quantitative agreement with experiment. The resulting NOCI eigenvalues are shifted by the average of the (scalar) OO-DFT CEBEs to yield CEBEs (split by SOC) corrected for dynamic correlation. Comparing calculations on gas phase molecules with experimental results establishes that NO-QDPT with the SCAN functional (NO-QDPT/SCAN), using the DCB screened 1-electron SOC operator is accurate to about 0.2 eV for L-edge CEBEs of molecules containing third row atoms. However, this NO-QDPT approach becomes less accurate for fourth-row elements starting in the middle of the 3d transition metal series, with errors increasing as atomic number increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)

This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting the oxidation states of Mn ions in the oxygen-evolving complex of photosystem II using supervised and unsupervised machine learning

Abstract Serial Femtosecond Crystallography at the X-ray Free Electron Laser (XFEL) sources enabled the imaging of the catalytic intermediates of the oxygen evolution reaction of Photosystem II (PSII). However, due to the incoherent transition of the S-states, the resolved structures are a convolution from different catalytic states. Here, we train Decision Tree Classifier and K-means clustering models on Mn compounds obtained from the Cambridge Crystallographic Database to predict the S-state of the X-ray, XFEL, and CryoEM structures by predicting the Mn’s oxidation states in the oxygen-evolving complex. The model agrees mostly with the XFEL structures in the dark S 1 state. However, significant discrepancies are observed for the excited XFEL states (S 2 , S 3, and S 0 ) and the dark states of the X-ray and CryoEM structures. Furthermore, there is a mismatch between the predicted S-states within the two monomers of the same dimer, mainly in the excited states. We validated our model against other metalloenzymes, the valence bond model and the Mn spin densities calculated using density functional theory for two of the mismatched predictions of PSII. The model suggests designing a more optimized sample delivery and illumiation systems are crucial to precisely resolve the geometry of the advanced S-states to overcome the noncoherent S-state transition. In addition, significant radiation damage is observed in X-ray and CryoEM structures, particularly at the dangler Mn center (Mn4). Our model represents a valuable tool for investigating the electronic structure of the catalytic metal cluster of PSII to understand the water splitting mechanism.

Plant Sciences↗

Dirac nodal arc in 1T-VSe 2

Transition metal dichalcogenides exhibit many fascinating properties including superconductivity, magnetic orders, and charge density wave. The combination of these features with a non-trivial band topology opens the possibility of additional exotic states such as Majorana fermions and quantum anomalous Hall effect. Here, we report on photon-energy and polarization dependent spin-resolved angle-resolved photoemission spectroscopy experiments on single crystal 1T-VSe 2 , revealing an unexpected band inversion and emergent Dirac nodal arc with spin-momentum locking. Density functional theory calculations suggest a surface lattice strain could be the driving mechanism for the topologically nontrivial electronic structure of 1T-VSe 2 .

36 MATERIALS SCIENCE↗

Identifying hydroxylated copper dimers in SSZ-13 via UV-vis-NIR spectroscopy

Cu-Exchanged zeolites are promising materials for the selective conversion of methane to methanol. Their activity is attributed to the presence of small Cu-oxo and Cu-hydroxy clusters, but the nature of active centers in various zeolite structures is still under debate. In this contribution, we combine time dependent density functional theory with spin–orbit coupling to predict the optical spectra of various Cu monomers and dimers in SSZ-13. As a result, we furthermore compare theoretical results to experimental measurements and find that the presence of Cu-hydroxy dimers and Cu monomers could potentially explain the experimentally observed UV-vis-NIR spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Native point defects in HgCdTe infrared detector material: Identifying deep centers from first principles

We investigate the native point defects in the long-wavelength infrared (LWIR) detector material Hg0.75Cd0.25Te using a dielectric-dependent hybrid density functional combined with spin–orbit coupling. Characterizing these point defects is essential as they are responsible for intrinsic doping and nonradiative recombination centers in the detector material. The dielectric-dependent hybrid functional allows for an accurate description of the bandgap (Eg) for Hg1−xCdxTe (MCT) over the entire compositional range, a level of accuracy challenging with standard hybrid functionals. Our comprehensive examination of the native point defects confirms that cation vacancies VHg(Cd) are the primary sources of p-type conductivity in the LWIR material given their low defect formation energies and the presence of a shallow acceptor level (−/0) near the valence-band maximum. In addition to the shallow acceptor level, the cation vacancies exhibit a deep charge transition level (2−/−) situated near the midgap, characteristic of nonradiative recombination centers. Our results indicate that Hg interstitial could also be a deep center in the LWIR MCT through a metastable configuration under the Hg-rich growth conditions. While an isolated Te antisite does not show deep levels, the formation of VHg–TeHg defect complex introduces a deep acceptor level within the bandgap.

Physics↗

Ab Initio Evaluation of Uranium Carbide S(α,β) and Thermal Neutron Cross Sections

Uranium Carbide (UC) is a nuclear fuel material which offers better neutron economy and lower fuel-cycle costs compared to conventional mixed-oxide fuels. UC’s lattice binding and dynamical properties impact thermal neutron scattering and low temperature epithermal resonance absorption. The Thermal Scattering Law (TSL) describes the scattering system available energy and momentum transfer states. There is no TSL evaluation for UC in the ENDF/B-VIII.0 database; herein, ab-initio lattice dynamics (AILD) techniques are invoked to calculate the phonon spectrum for UC using spin-orbit-coupling density functional theory (DFT). The TSLs, inelastic and elastic thermal scattering cross sections for Uranium and Carbon in UC, respectively, are calculated in FLASSH for use in higher fidelity reactor design calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ab Initio Evaluation of Plutonium Dioxide $S(α,β)$ and Thermal Neutron Cross Sections

Plutonium Dioxide (PuO 2 ) is an advanced fuel-cycle nuclear fuel material for thermal and fast reactors, with a melting point at 3261.8 K and thermal conductivity above 3 Wm -1 K -1 for temperatures below 1500 K. Plutonium from spent UO 2 is reprocessed to form chemically stable mixed oxide (MOX) fuel. In contrast to metallic fuels, the low PuO 2 conductivity results in high thermal gradients between fuel center and surface, resulting in structural annealing and efficient fission product migration over burnup. Crystal binding affects the doppler broadening of epithermal resonances for nuclear fuels and scattering behavior of thermal neutrons. Currently, there is no ENDF/B Thermal Scattering Law, i.e., $S(α,β)$, evaluation for PuO 2 , only free atom cross sections for Plutonium and Oxygen, respectively, which neglect lattice contributions to total cross section. Herein, ab initio lattice dynamics (AILD) techniques are employed to calculate the phonon density of states (DOS) using spin-orbit-coupling density functional theory (DFT) to predict the PuO2 paramagnetic ground-state structure. These results highlight the evaluation of $S(α,β)$ for Plutonium and Oxygen, respectively in PuO 2 and consequential generation of thermal neutron scattering cross sections for high fidelity criticality safety analysis and reactor calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

One-dimensional Rashba states with unconventional spin texture in Bi chains

Spin-polarized electrons confined in low-dimensional structures are of high interest for spintronics applications. Here, in this study, we investigate the electronic structure of an ordered array of Bi monomer and dimer chains on the Ag(110) surface. By means of spin-resolved photoemission spectroscopy, we find Rashba-Bychkov split bands crossing the Fermi level with one-dimensional constant energy contours. These bands are up-spin polarized for positive wave vectors and down-spin polarized for negative wave vectors, at variance with the Rashba-Bychkov model that predicts a pair of states with opposite spin in each half of the surface Brillouin zone. Density functional theory shows that spin-selective hybridization with the Ag bulk bands originates this unconventional spin texture.

1-dimensional spin chains↗

Magnetic interactions and spin excitations in van der Waals ferromagnet VI 3

Using a combination of density functional theory (DFT) and spin-wave theory methods, we investigate the magnetic interactions and spin excitations in semiconducting VI 3 . Exchange parameters of monolayer, bilayer, and bulk forms are evaluated by mapping the magnetic energies of various spin configurations, calculated using DFT + U, onto the Heisenberg model. The intralayer couplings remain largely unchanged in three forms of VI 3 , while the interlayer couplings show stronger dependence on the dimensionality of the materials. We calculate the spin-wave spectra within a linear spin-wave theory and discuss how various exchange parameters affect the magnon bands. Here, the magnon–magnon interaction is further incorporated, and the Curie temperature is estimated using a self-consistently renormalized spin-wave theory. To understand the roles of constituent atoms on magnetocrystalline anisotropy energy (MAE), we resolve MAE into sublattices and find that a strong negative V-I inter-sublattice contribution is responsible for the relatively small easy-axis MAE in VI 3 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

One-bond 13 C– 13 C spin-coupling constants in saccharides: a comparison of experimental and calculated values by density functional theory using solid-state 13 C NMR and X-ray crystallography

Methyl aldohexopyranosides were 13 C-labeled at contiguous carbons, crystallized, and studied by single-crystal X-ray crystallography and solid-state 13 C nuclear magnetic resonance (NMR) spectroscopy to examine the degree to which density functional theory (DFT) can calculate one-bond 13C–13C spin-coupling constants ( 1 J CC ) in saccharides with sufficient accuracy to permit their use in MA'AT analysis, a newly-reported hybrid DFT/NMR method that provides probability distributions of molecular torsion angles in solution (Zhang et al., J. Phys. Chem. B, 2017, 121, 3042–3058; Meredith et al., J. Chem. Inf. Model., 2022, 62, 3135–3141). Experimental 1 J CC values in crystalline samples of the doubly 13 C-labeled compounds were measured by solid-state 13 C NMR and compared to those calculated from five different DFT models: (1) 1 J CC values calculated from single structures identical to those observed in crystalline samples by X-ray crystallography (all atom refinement); (2) 1 J CC values calculated from the single structures in (1) but after Hirshfeld atom refinement (HAR); (3) 1 J CC values calculated from the single structures in (1) after DFT-optimization of hydrogen atoms only; and (4 and 5) 1 J CC values calculated in rotamers of torsion angle θ 2 (C1–C2–O2–O2H) or ω(C4–C5–C6–O6) from which either specific or generalized parameterized equations were obtained and used to calculate 1 J CC values in the specific θ 2 or ω rotamers observed in crystalline samples. Good qualitative agreement was observed between calculated 1 J CC values and those measured by solid-state 13 C NMR regardless of the DFT model, but in no cases were calculated 1 J CC values quantitative, differing (over-estimated) on average by 4–5% from experimental values. These findings, and those reported recently from solution NMR studies (Tetrault et al., J. Phys. Chem. B 2022, 126, 9506–9515), indicate that improvements in DFT calculations are needed before calculated 1 J CC values can be used directly as reliable constraints in MA'AT analyses of saccharides in solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of random antiferromagnetic exchange on the spin waves in a three-dimensional Heisenberg ferromagnet

Neutron scattering is used to study spin waves in the three-dimensional Heisenberg ferromagnet YTiO 3 , with spin-spin exchange disorder introduced via La substitution at the Y site. No significant changes are observed in the spin-wave dispersion up to a La concentration of 20%. However, a strong broadening of the spectrum is found, indicative of shortened spin-wave lifetimes. Density-functional theory calculations predict minimal changes in exchange constants as a result of average structural changes due to La substitution, in agreement with the data. The absence of significant changes in the spin-wave dispersion, the considerable lifetime effect, and the reduced ordered magnetic moment previously observed in the La-substituted system are qualitatively captured by an isotropic, nearest-neighbor, three-dimensional Furthermore, Heisenberg ferromagnet model with random antiferromagnetic exchange. We therefore establish Y 1-$x$ La $x$ TiO 3 as a model system for the study of antiferromagnetic spin-exchange disorder in a three-dimensional Heisenberg ferromagnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Assessing spin-density wave formation in La 3 Ni 2 O 7 from electronic structure calculations

Here, we employ correlated density-functional theory methods (DFT + Hubbard U) to investigate the spin-density wave state of the bilayer Ruddlesden-Popper (RP) nickelate La 3 Ni 2 O 7 which becomes superconducting under pressure. We predict that the ground state of this bilayer RP material has traits of both the double spin-stripe and the single spin-charge stripe phases proposed in the literature as it corresponds to in-plane up/up′/down/down′ diagonal stripes with up/down being high spin (formally Ni 2+ : d 8 ), and up′ /down′ being low spin (formally Ni 3+ : d 7 ). The main feature of this solution (that is insulating even at U = 0) is the dominant role of the d x 2 −y 2 bands around the Fermi level, which would become doped with the introduction of electrons via oxygen vacancies. In spite of the similarity with cuprates in terms of the dominant role of d x 2 −y 2 bands, some differences are apparent in the magnetic ground state of La 3 Ni 2 O 7 : the antiferromagnetic out-of-plane coupling within the bilayer (linked to the d z 2 orbitals forming a spin-singlet-like configuration) is found to be the dominant one while in-plane interactions are reduced due to the stripe order of the ground state. With pressure, this striped magnetic ground state remains similar in nature but the increase in bandwidth quickly transitions La 3 Ni 2 O 7 into a metallic state with all the activity close to the Fermi level involving, to a large extent, d x 2 −y 2 orbitals. This is reminiscent of the cuprates and may provide key insights into how superconductivity arises in this material under pressure.

LaBollita, Harrison↗

CHARGED POINT DEFECTS AND DEFECT COMPLEXES IN PHOSPHORENE AND STRUCTURAL PREDICTION OF TAU-11, A DENSITY FUNCTIONAL THEORY APPROACH

With the advent of powerful computers, first principle calculations have become a robust and reliable tool to complement experiments in the study of materials. First principle calculations such as density functional theory can be used to guide experiments to save time and effort when it comes to new and advanced materials discovery. Study of charge defects in materials and crystal structure prediction are such areas in which density functional theory calculations facilitate the prediction of defect formation energy, charge transition levels and elemental site occupations. Experimentally, studying charge defects requires inferring results from various techniques. Likewise, determining site occupancies for similar elements using experimental techniques is challenging and sometimes impossible. In this work, we study the charged intrinsic and extrinsic defects and defect complexes in phosphorene and determine the site occupancies of Al, Fe and Si in the τ11 phase using density functional theory. In the defect work, we calculate the defect formation energy, charge transition level, binding energies and Stokes shift for vacancy, dopant substitution and dopant-vacancy defect complexes in phosphorene. We found that vacancy defect in phosphorene becomes negatively charged in n-type doping and may passivate the dopants and reduce carrier concentration and mobility. For non-metal dopants in phosphorene, we predict that O, S and Si prefer to form dopant-vacancy complexes removing the vacancy defect states from the band gap. Mn dopant-vacancy defect complex exhibits possibility of switching between two magnetic spin states. Lastly, using density functional theory, we complement the results from neutron powder diffraction to determine the site occupancies of Al and Si in τ11 phase which was used to determine the alloying element to develop a high temperature and low density alloy.

36 MATERIALS SCIENCE↗

Coherent Spins in van der Waals Semiconductor GeS 2 at Ambient Conditions

Optically active spin defects in van der Waals (vdW) materials have emerged as versatile quantum sensors, enabling applications for a wide range of quantum phenomena in condensed matter systems. Their ease of exfoliation and compatibility with device integration make them promising candidates for future quantum technologies. Here we report the observation and room-temperature coherent control of ensemble spin defects in the high-temperature crystalline phase of germanium disulfide (β-GeS 2 ), a two-dimensional (2D) semiconductor with low nuclear spin density. The defects exhibit spin-1/2 behavior, and their dynamics can be explained by a weakly coupled spin-pair model. We implement dynamical decoupling techniques to extend the coherence time (T 2 ) by a factor of 20. Finally, we use density functional theory (DFT) calculations to estimate the structures and spin densities of two possible spin defect candidates. This work will help to expand the field of quantum sensing with spin defects in 2D materials.

2D materials↗

Modeling surface spin polarization on ceria-supported Pt nanoparticles

In this work, we employ density functional theory simulations to investigate possible spin polarization of CeO 2 -(111) surface and its impact on the interactions between a ceria support and Pt nanoparticles. With a Gaussian type orbital basis, our simulations suggest that the CeO 2 -(111) surface exhibits a robust surface spin polarization due to the internal charge transfer between atomic Ce and O layers. In turn, it can lower the surface oxygen vacancy formation energy and enhance the oxide reducibility. We show that the inclusion of spin polarization can significantly reduce the major activation barrier in the proposed reaction pathway of CO oxidation on ceria-supported Pt nanoparticles. For metal-support interactions, surface spin polarization enhances the bonding between Pt nanoparticles and ceria surface oxygen, while CO adsorption on Pt nanoparticles weakens the interfacial interaction regardless of spin polarization. However, the stable surface spin polarization can only be found in the simulations based on the Gaussian type orbital basis. Furthermore, given the potential importance in the design of future high-performance catalysts, our present study suggests a pressing need to examine the surface ferromagnetism of transition metal oxides in both experiment and theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗