Search NASA⌕ Search

SEARCH · Search NASA

Results for “Natural orbitals”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Actinide 2-metallabiphenylenes that satisfy Hückel’s rule

Aromaticity and antiaromaticity, as defined by Hückel’s rule, are key ideas in organic chemistry, and are both exemplified in biphenylene—a molecule that consists of two benzene rings joined by a four-membered ring at its core. Biphenylene analogues in which one of the benzene rings has been replaced by a different (4n + 2) π-electron system have so far been associated only with organic compounds. In addition, efforts to prepare a zirconabiphenylene compound resulted in the isolation of a bis(alkyne) zirconocene complex instead. Here we report the synthesis and characterization of, to our knowledge, the first 2-metallabiphenylene compounds. Single-crystal X-ray diffraction studies reveal that these complexes have nearly planar, 11-membered metallatricycles with metrical parameters that compare well with those reported for biphenylene. Nuclear magnetic resonance spectroscopy, in addition to nucleus-independent chemical shift calculations, provides evidence that these complexes contain an antiaromatic cyclobutadiene ring and an aromatic benzene ring. Finally, spectroscopic evidence, Kohn–Sham molecular orbital compositions and natural bond orbital calculations suggest covalency and delocalization of the uranium f 2 electrons with the carbon-containing ligand.

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

Experimental and quantum mechanical characterization of an oxygen-bridged plutonium(IV) dimer

We report the synthesis and characterization of K 4 {[PuCl 2 (NO 3 ) 3 ] 2 (μ 2 -O)}·H 2 O, which contains the first known μ 2 -oxo bridge between two Pu IV metal centers. Adding to its uniqueness is the Pu-(μ 2 -O) bond length of 2.04 Å, which is the shortest of other analogous compounds. The Pu-(μ 2 -O)-Pu bridge is characterized by the mixing of s -, d -, and p -orbitals from Pu with the p -orbitals of O; the 5 f -orbitals do not participate in bonding. Natural bond orbital analysis indicates that Pu and O interact through one 3c-2e σ Pu-O-Pu and two 3c-2e π Pu-O-Pu bonding orbitals and that the electron density is highly polarized on the μ 2 -O. Bond topology properties analysis indicates that the Pu-(μ 2 -O) bond shares both ionic and covalent character. Quantum mechanical calculations also show that the dimer has multiconfigurational ground states, where the nonet, septet, quintet, triplet, and singlet are close in energy. This work demonstrates the interplay between experimental and computational efforts that is required to understand the chemical bonding of Pu compounds.

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

Magnetic Compton Scattering Study of Li-Rich Battery Materials

The redox process in a lithium-ion battery occurs when a conduction electron from the lithium anode is transferred to the redox orbital of the cathode. Understanding the nature of orbitals involved in anionic as well as cationic redox reactions is important for improving the capacity and energy density of Li-ion batteries. In this connection, we have obtained magnetic Compton profiles (MCPs) from the Li-rich cation-disordered rock-salt compound LixTi0.4Mn0.4O2 (LTMO). The MCPs, which involved the scattering of circularly polarized hard X-rays, are given by the momentum density of all the unpaired spins in the material. The net magnetic moment in the ground state can be extracted from the area under the MCP, along with a SQUID measurement. Our analysis gives insight into the role of Mn 3d magnetic electrons and O 2p holes in the magnetic redox properties of LTMO.

Suzuki, Kosuke (ORCID:0000000317110604)↗
DOE OSTI2022

Interaction of Th with H 0/–/+ : Combined Experimental and Theoretical Thermodynamic Properties

In this work, high-level electronic structure calculations of the lowlying energy electronic states for ThH, ThH – , and ThH + are reported and compared to experimental measurements. The inclusion of spin–orbit coupling is critical to predict the ground-state ordering as inclusion of spin–orbit switches the coupled-cluster CCSD(T) ordering of the two lowest energy states for ThH and ThH + . At the multireference spin–orbit SO-CASPT2 level, the ground states of ThH, ThH – , and ThH + are predicted to be the 2 Δ 3/2 , 3 Φ 2 , and 3 Δ 1 states, respectively. The adiabatic electron affinity is calculated to be 0.820 eV, and the vertical detachment energy is calculated to be 0.832 eV in comparison to an experimental value of 0.87 ± 0.02 eV. The observed ThH – photoelectron spectrum has many transitions, which approximately correlate with excitations of Th + and/or Th. The adiabatic ionization energy of ThH including spin–orbit corrections is calculated to be 6.181 eV. The natural bond orbital results are consistent with a significant contribution of the Th + H – ionic configuration to the bonding in ThH. The bond dissociation energies for ThH, ThH – , and ThH + using the Feller–Peterson–Dixon approach were calculated to be similar for all three molecules and lie between 259 and 280 kJ/mol.

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

Overdestabilization vs Overstabilization in the Theoretical Analysis of f-Orbital Covalency

The complex nature of the f-orbital electronic structures and their interaction with the chemical environment pose significant computational challenges. Advanced computational techniques that variationally include scalar relativities and spin–orbit coupling directly at the molecular orbital level have been developed to address this complexity. Among these, variational relativistic multiconfigurational multireference methods stand out for their high accuracy and systematic improvement in studies of f-block complexes. Additionally, these advanced methods offer the potential for calibrating low-scaling electronic structure methods such as density functional theory. However, studies on the Cl K-edge X-ray absorption spectra of the [Ce(III)Cl 6 ] 3– and [Ce(IV)Cl 6 ] 2– complexes show that time-dependent density functional theory with approximate exchange–correlation kernels can lead to inaccuracies, resulting in an overstabilization of 4f orbitals and incorrect assessments of covalency. In contrast, approaches utilizing small active space wave function methods may understate the stability of these orbitals. The results herein demonstrate the need for large active space, multireference, and variational relativistic methods in studying f-block complexes.

chemical calculations↗
DOE OSTI2021

Imaging orbital ferromagnetism in a moiré Chern insulator

Capturing exotic magnetism Ferromagnetism is associated with the breaking of time-reversal symmetry, most frequently by the spin degree of freedom. Although the orbital motion of electrons can also contribute to ferromagnetism, in most materials, it is small relative to the spin contribution. Tschirhart et al. showed that the reverse is true in an unusual magnetic state hosted by twisted bilayer graphene. Their scanning magnetometry measurements were consistent with ferromagnetism of predominantly orbital origin. Science , abd3190, this issue p. 1323

Science & Technology - Other Topics↗
DOE OSTI2026

An Algorithm for Atom-Centered Lossy Compression of the Atomic Orbital Basis in Density Functional Theory Calculations

Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-ζ (QZ) size are required to adequately converge Kohn–Sham density functional theory (DFT) calculations toward the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. Here, in this work, we develop and test a threshold-based natural atomic orbital (NAO) scheme in which ϵ-NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. This enables compression of the AO basis that is optimal for a given threshold, 10 –ϵ , by discarding NAOs with occupation numbers below that threshold. Extensive pilot test calculations using the Hartree–Fock functional and taking the converged density matrix as input suggest that a threshold of 10 –5 can yield a compression factor (ratio of AO to compressed ϵ-NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold 10 –7 , while the compression factor is typically between 2 and 2.5.

basis sets↗
DOE OSTI2024

Properties of NoO and YbO

Bond dissociation energies (BDEs) and spectroscopic parameters for the ground states of YbO and NoO were calculated at the Brueckner doubles (BD(T)) level and for NoO at the coupled cluster CCSD(T) level based on the closed-shell configurations Yb 2+ (4f 14 )O 2 – (2p 6 ) and No 2+ (5f 14 )O 2 – (2p 6 ). For YbO, the BD(T) BDE (3.93 eV) and vibrational frequency (686.3 cm –1 ) are consistent with the experiment; the BD(T) bond distance was 1.8492 Å, about 0.04 Å longer than the experiment. For NoO, the calculated BDEs are 3.07 eV (BD(T)) and 3.08 eV (CCSD(T)) and are expected to be good to ±3 kcal/mol (±0.13 eV). Potential energy curves (PECs) for the spin–orbit (Ω) states associated with the f 13 s 1 configurations of YbO and NoO were calculated at the SA-CASSCF/SO-CASPT2/aQ-DK level. PECs considering both the f 13 s 1 and f 14 configurations were calculated at this level, as the ground states are expected to contain both configurations. Due to a bias in the SA-CASSCF to the f 13 s 1 configurations, the f 14 configurations could not be properly described simultaneously with the f 13 s 1 ones. Finally, the natural population analysis based on the natural bond orbitals shows that YbO and NoO are quite ionic, with 5d orbitals being more important in YbO than the 6d orbitals in NoO.

basis sets↗
DOE OSTI2022

Molecular Properties of Thorium Hydrides: Electron Affinities and Thermochemistry

In this work, high-level electronic structure calculations of the ground and low-lying energy electronic states for ThH x and ThH x – for x = 2–5 are reported and compared to available anion photoelectron detachment experiments. The adiabatic electron affinities (EAs) are predicted to be 0.82, 0.88, 0.51, and 2.36 eV for x = 2 to 5, respectively, at the Feller– Peterson–Dixon (FPD) level. The vertical detachment energies (VDEs) are predicted to be 0.84, 0.88, 0.81, and 4.38 eV for x = 2–5, respectively. The corresponding experimental VDEs are 0.871 eV for x = 2, 0.88 eV for x = 3, and 4.09 eV for x = 5. As for ThH, there is a significant spin–orbit (SO) correction for the EA of ThH 2 , and this correction decreases substantially for x > 2. The observed ThH 2 – photoelectron spectrum has many transitions as predicted at the CASPT2-SO level. The FPD bond dissociation energies (BDEs) increase from 67 to 75 kcal/mol for x = 2 to x = 4 at the FPD level. The BDE for ThH 5 is much lower as it is a complex of H 2 with ThH 3 . The hydride affinities for x = 2 to 4 are all comparable and near 70 kcal/mol. A natural bond orbital analysis is consistent with a significant Th + –H– ionic contribution to the Th–H bonds. There is very little participation of the 5f orbitals in the bonding and the valence electrons on the Th are dominated by 7s and 6d for the neutrals and anions except for ThH 2 – where there is a significant contribution from the 7p.

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

Energetic and Electronic Properties of UX +/0/– for X = Li and Be and Comparison of the Properties of the Uranium Atom Binding to 2nd Row Elements Li–F

The bonding and spectroscopic properties of ULi +/0/– and UBe +/0/– to complete the series for UX +/0/– for X = Li to F were investigated by high-level ab initio SO-CASPT2 and CCSD(T) electronic structure calculations. The low-lying spin–orbit states were obtained at the SA-CASPT2/aQ-PP level; bond dissociation energies (BDEs), ionization energies (IEs), adiabatic electronic affinities (AEAs), and vertical detachment energies (VDEs) were calculated at the Feller-Peterson-Dixon (FPD) level. A dense manifold of low-lying states was predicted for ULi +/0/– and UBe +/0/– . Here, the calculated BDEs for ULi (37.7 kJ/mol) and UBe (8.0 kJ/mol) show that UBe is weakly bound. For redox processes, the BDEs increased for ULi + (109.3 kJ/mol), ULi – (47.4 kJ/mol), UBe + (35.6 kJ/mol), and UBe – (72.3 kJ/mol). The IE(ULi) = 4.650 eV is lower than IE(Li); the IE(UBe) = 5.901 eV is close to the IE(U) and to the IEs of UB, UC, UN, UO, and UF. The AEAs of ULi (0.708 eV) and UBe (0.989 eV) are lower than those for UB, UC, UN, and UO but higher than that for EA(UF). Natural bond orbital (NBO) calculations show that ULi has the 5f 3 6d 1 7s 2 configuration for U and 2s 1 for Li, with a small partial negative charge slightly delocalized on U. UBe arises from the U(5f 3 6d 1 7s 2 ) and Be(2s 2 ) electron configurations with no charge separation. The same calculations were made for WX (X = Li, Be, C–F) to enable detailed comparisons of the properties for UX with WX (X = Li–F). For WX, BDE(WX) is higher than that for UX for X = Li to N and lower than BDE(UX) for X = O and F, mostly due to the higher IE of W than U as ionic character becomes more important going from Li to F.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Experimental and Computational Description of the Interaction of H and H – with U

In this work, the results of ab initio correlated molecular orbital theory electronic structure calculations for low-lying electronic states are presented for UH and UH – and compared to photoelectron spectroscopy measurements. The calculations were performed at the CCSD(T)/CBS and multireference CASPT2 including spin–orbit effects by the state interacting approach levels. The ground states of UH and UH – are predicted to be 4 Ι 9/2 and 5 Λ 6 , respectively. The spectroscopic parameters T e , r e , ω e , ω e x e , and Be were obtained, and potential energy curves were calculated for the low energy Ω states of UH. The calculated adiabatic electron affinity is 0.468 eV in excellent agreement with an experimental value of 0.462 ± 0.013 eV. The lowest vertical detachment energy was predicted to be 0.506 eV for the ground state, and the adiabatic ionization energy (IE) is predicted to be 6.116 eV. The bond dissociation energy (BDE) and heat of formation values of UH were obtained using the IE calculated at the Feller–Peterson–Dixon level. For UH, UH – , and UH + , the BDEs were predicted to be 225.5, 197.9, and 235.5 kJ/mol, respectively. The BDE for UH is predicted to be ~20% lower in energy than that for ThH. The analysis of the natural bond orbitals shows a significant U + H – ionic component in the bond of UH.

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

Structure and Bonding of 1,2,4-Triazole Thiones Derived from Nitron

The molecular structures of two 1,2,4-triazole thiones derived from Nitron, namely Nitron(S) and Nitron(S) Me , have been determined by single crystal X-ray diffraction. Comparison with the structure of Nitron indicates that the formation of the thiones is accompanied by a shortening of the C–N bond within the 5-membered ring that is opposite to the thione group, and a lengthening of the exocyclic C–N bond. Analysis of the structures of Nitron, Nitron(S) and Nitron(S) Me , together with hypothetical versions in which the phenyl substituents are replaced by hydrogen, indicates that these changes are in accord with bond order differences predicted by natural bond orbital methods. Here, a similar evaluation of the thione groups of Nitron(S) and Nitron(S)Me indicates that the C–S bonds are composed of both covalent and ionic components and are best described by a combination of resonance structures that feature C=S double and C–S single bonds.

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

Structure−Property Relationships in a Series of Hybrid Halopyridinium Ruthenium(IV) Halides

Reported is the synthesis and characterization of a family of hybrid ruthenium halides consisting of haloruthenium octahedra (X = Cl, Br) charge balanced with halopyridinium (XPy; X= H, Cl, Br, I) organic cations which assemble via noncovalent interactions between ion pairs. Diffuse reflectance spectroscopy showed that compounds containing RuBr 6 2- octahedra displayed a lower bandgap (1.05 eV< x < 1.08 eV) compared to compounds with RuCl 6 2- octahedra (1.22 eV < x < 1.43 eV). Additionally, computational density functional theory (DFT) based natural bonding orbital (NBO) analysis and density of state (DOS) methods were used to characterize second-sphere non-covalent interaction strengths and elucidate molecular orbital perturbations to elucidate their influence on Ru-X orbital constructs and in turn rationalize band gap trends. Analyses showed ligand to metal charge transfer is responsible for the small bandgap energies and are unperturbed by second-sphere interactions. Here, this report serves as a platform for probing the relationship between the structural and photophysical properties within the haloruthenium family of low dimensional transition metal halide perovskites.

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

Theoretical Evaluation of Metal–Ligand Bonding in Neptunium Compounds in Relation to 237 Np Mössbauer Spectroscopy

The 237 Np Mössbauer isomer shift and quadrupole splitting (QS) are powerful probes for the metal–ligand bonding of neptunium, a 5f-element of vital importance in the nuclear fuel cycle. A large set of Np compounds with different oxidation states (III) to (VII) is studied to investigate, by first-principles calculations, isomer shifts and the QS trends in relation to the Np oxidation state. Natural Bond Orbital analysis reveals that in addition to donation bonding to the 5f shell, participation of the 6d and 7s neptunium shells in covalent (donation) bonding substantially impacts the isomer shifts. The isomer shift cannot be interpreted solely by the 5f shell electron count. The isomer shift for Np(II) compounds is estimated to be in the range of 31–34 mm/s, less positive than for Np(III) compounds. For the QS, density functional calculations fail to reproduce the quadrupole splitting for some Np(VI) ionic solids. A multiconfigurational wave function approach reproduces the observed QS trends. The calculations give a semiquantitative interpretation of the trends for Np oxidation states (V) to (VII). Furthermore, the contrasting QS for standard and “reverse” neptunyl(VI), at the opposite extremes of the observed QS scale, arises predominantly from the different crystal environments.

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

Dinitrogen Binding and Functionalization from a Low-Coordinate Alkynyliron Complex

Alkynyl complexes of low-coordinate transition metals offer a sterically open environment and interesting bonding opportunities. Here, in this work, we explore the capacity of iron(I) alkynyl complexes to bind N 2 and isolate a N 2 complex including its X-ray crystal structure. Silylation of the N 2 complex gives an isolable, formally iron(IV) complex with a disilylhydrazido(2−) ligand, but natural bond orbital analysis indicates that an iron(II) formulation is preferable. The structure of this compound is similar to an earlier reported phenyl complex in which phenyl migration forms a new N–C bond, but the alkynyl group does not migrate. DFT calculations are used to test the possible reasons why the alkynyl is resistant to migration, and these show that the large Fe–C bond energy in the alkynyl complex is a factor that could contribute to the lack of migration.

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

Linear versus Bent Uranium(II) Metallocenes—A Local Vibrational Mode Study

Uranium metallocenes have recently attracted attention driven by their use as catalysts in organometallic synthesis. In addition to bent U(IV) and U(III), an U(II) metallocene [(η 5 -C 5 i Pr 5 ) 2 U] was synthesized with an unusual linear Cp–U–Cp angle. Here, in this work, we investigated 22 U(II) metallocenes, (i) assessing the intrinsic strength of the U-ring interactions in these complexes with a novel bond strength measure based on our local vibrational mode analysis and (ii) systematically exploring what makes these U(II) metallocenes bent. We included relativistic effects through the NESCau Hamiltonian and complemented the local mode analysis with natural bonding orbital (NBO) and quantum theory of atoms in molecules (QTAIM) data. Our study led to the following results: (i) reduction of bulky U-ring ligand substituents does not lead to bent complexes for alkyl substituents ( i Pr and i Bu) in contrast to SiMe 3 ring substituents, which are all bent. (ii) The most bent complexes are [(η 5 -C 5 H 4 SiMe 3 ) 2 U] (130°) and [η 5 -P 5 H 5 ) 2 U] (143°). (iii) Linear complexes showed one hybridized NBO with s/d character, while bent structures were characterized by s/d/f mixing. (iv) We did not observe a correlation between the strength of the U-ring interaction and the amount of the ring-U-ring bend; the strongest interaction was found for [η 5 -Cp) 2 U] and the weakest for [η 5 -P 5 H 5 ) 2 U]. In conclusion, our results provide a foundation for the design of U(II) metallocenes with specific physicochemical properties and increased reactivity.

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

N 5 -Scaling Excited-State-Specific Perturbation Theory

In this paper, we show that by working in a basis similar to that of the natural transition orbitals and using a modified zeroth-order Hamiltonian, the cost of a recently introduced perturbative correction to excited-state mean field theory can be reduced from seventh to fifth order in the system size. The (occupied) 2 (virtual) 3 asymptotic scaling matches that of ground-state second-order Møller-Plesset theory but with a significantly higher prefactor because the bottleneck is iterative: it appears in the Krylov-subspace-based solution of the linear equation that yields the first-order wave function. Here, we discuss the details of the modified zeroth-order Hamiltonian we use to reduce the cost and the automatic code generation process we used to derive and verify the cost scaling of the different terms. Overall, we find that our modifications have little impact on the method's accuracy, which remains competitive with singles and doubles equation-of-motion coupled cluster.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗