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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1975

Guidance of an aeromaneuvering orbit-to-orbit shuttle through a statistically varying atmosphere

One candidate for a reusable upper stage to be carried by the space shuttle is an aeromaneuvering orbit-to-orbit shuttle (AMOOS). This concept uses the drag of the vehicle during a pass through the atmosphere rather than the propulsion system to slow the vehicle on a return from a high energy orbit. The nature and magnitude of the sensitivity of AMOOS to uncertainties in the properties of the atmosphere are shown. Various guidance schemes for correcting for the effects that the unpredictable variations in the atmosphere have on the trajectory are discussed. For the mission studied here, a payload retrieval from geosynchronous orbit with aerodynamic plane change, a linear feedback guidance scheme was developed. A relatively simple heuristic law was used to demonstrate the concept. Using optimal control theory a feedback law was developed analytically. Testing with a large number of different atmospheres showed this law to be a feasible means of controlling the AMOOS trajectory. Refinements to the technique offer promise of significant improvement, and these are discussed.

Rehder, J. J.↗
NASA NTRS2007

Semianalytical Propagation of Satellite Orbits about an Arbitrary Central Body

Precision mean element (PME) satellite theories play a key role in orbit dynamics analyses. These theories employ: nonsingular orbital elements comprehensive force models Generalized Method of Averaging Numerical interpolation concepts The Draper Semianalytical Satellite Theory (DSST) (Refs. 1 - 6), whose development was led by the author, and the independently-developed Universal Semianalytical Method (USM) (Ref. 7) are examples of such theories. These theories provide the capability to tailor the force modeling to meet the desired computational speed vs. accuracy trade-off. The flexibility of such theories is demonstrated by their ability to include complicated atmosphere density models and spacecraft models in the perturbation theory context. The value of high speed satellite theories, in this era of computational plenty, is that they allow new ways of looking at astrodynamical problems such as orbit design (Refs. 8, 9) and atmosphere density updating (Refs. 10, 11). In the mid to late-1980 s, the geodynamics community led the development of very precise geopotential models such as GEM T2 and GEM T3 (Ref. 12), and with the subsequent analysis of the TOPEX flight data, JGM-2 and JGM-3 (Ref. 13). These were high degree and order geopotentials, at least 50 x 50. In 1993, the DSST implementation in the GTDS program was extended to include the 50 x 50 geopotential models (Ref. 14). The 50 x 50 geopotential, J2000 integration coordinate system, and solid Earth tide capabilities were integrated in GTDS by Scott Carter (Ref. 15). This capability demonstrated 1 m accuracy versus the TOPEX Precise Orbit Ephemerides. Subsequently the DSST Standalone program was also extended to include high degree and order geopotential models (Ref. 5). More recently GTDS has been hosted in the Linux PC environment. However, all of these efforts have been limited to modeling the motion of an artificial Earth satellite. They did not consider the additional complexities associated with lunar, planetary, or other natural satellite orbiters. Such complexities include: additional coordinate systems (associated with the direction of the north pole of rotation and the prime meridian of the new central bodies) (Ref. 16) normalized gravity model coefficients (desirable for high degree and order fields) (Ref. 17) indirect oblateness

Cefola, Paul J.↗
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↗
DOE OSTI2024

Revisiting Artifacts of Kohn–Sham Density Functionals for Biosimulation

We revisit the problem of unphysical charge density delocalization/fractionalization induced by the self-interaction error of common approximate Kohn–Sham (KS) density functional theory functionals on simulation of small to medium-sized proteins in a vacuum. Aside from producing unphysical electron densities and total energies, the vanishing of the HOMO–LUMO gap associated with the unphysical charge delocalization leads to an unphysical low-energy spectrum and catastrophic failure of most popular solvers for the KS self-consistent field (SCF) problem. We apply a robust quasi-Newton SCF solver to obtain solutions for some of these difficult cases. The anatomy of the charge delocalization is revealed by the natural deformation orbitals obtained from the density matrix difference between the Hartree–Fock and KS solutions; the charge delocalization not only can occur between charged fragments (such as in zwitterionic polypeptides) but also involves neutral fragments. The vanishing-gap phenomenon and troublesome SCF convergence are both attributed to the unphysical KS Fock operator eigenspectra of molecular fragments (e.g., amino acids or their side chains). Analysis of amino acid pairs suggests that the unphysical charge delocalization can be partially ameliorated by the use of some range-separated hybrid functionals but not by semilocal or standard hybrid functionals. Last, we demonstrate that solutions without the unphysical charge delocalization can be located even for semilocal KS functionals highly prone to such defects, but such solutions have non-Aufbau character and are unstable with respect to mixing of the non-overlapping “frontier” orbitals. Caution should be exercised when unexpectedly small (or vanishing) HOMO–LUMO gaps and atypical SCF convergence patterns (e.g., oscillatory) are observed in KS DFT simulations in any context (bio or otherwise).

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

Toward Quantum Confinement in Graphitic Carbon Nitride-Based Polymeric Monolayers

raphitic carbon nitride (g-C 3 N 4 ) has garnered much attention due to its potential as an efficient metal-free photocatalyst. This study examines the evolution of properties in zero-dimensional quantum dots up to sizable clusters that mimic extended g-C 3 N 4 monolayers. We employ density functional theory to investigate systematically the structural, electronic, and optical properties of the g-C 3 N 4 -based melamine and heptazine building blocks using a “bottom-up” construction of polymeric monolayers. The results from our computations indicate that the melamine- and heptazine-based polymeric g-C 3 N 4 systems must be reduced to at least 2.74 and 4.00 nm, respectively, to observe an increase of its optical gap with a size reduction. The present study also examines the nature of the electronic transitions exhibited by g-C 3 N 4 -based monolayers through full natural transition orbital and density of state analyses. Furthermore, the most promising sites for water splitting and subsequent chemical doping studies are identified, which generally correspond to the nitrogen and carbon atoms, respectively.

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

Theoretical and Experimental Study of the Spectroscopy and Thermochemistry of UC +/0/–

A combination of high-level ab initio calculations and anion photoelectron detachment (PD) measurements is reported for the UC, UC – , and UC + molecules. To better compare the theoretical values with the experimental photoelectron spectrum (PES), a value of 1.493 eV for the adiabatic electron affinity (AEA) of UC was calculated at the Feller–Peterson–Dixon (FPD) level. The lowest vertical detachment energy (VDE) is predicted to be 1.500 eV compared to the experimental value of 1.487 ± 0.035 eV. A shoulder to lower energy in the experimental PD spectrum with the 355 nm laser can be assigned to a combination of low-lying excited states of UC – and excited vibrational states. The VDEs calculated for the low-lying excited electronic states of UC at the SO-CASPT2 level are consistent with the observed additional electron binding energies at 1.990, 2.112, 2.316, and 3.760 eV. Potential energy curves for the Ω states and the associated spectroscopic properties are also reported. Compared to UN and UN + , the bond dissociation energy (BDE) of UC (411.3 kJ/mol) is predicted to be considerably lower. The natural bond orbitals (NBO) calculations show that the UC 0/+/– molecules have a bond order of 2.5 with their ground-state configuration arising from changes in the oxidation state of the U atom in terms of the 7s orbital occupation: UC (5f 2 7s 1 ), UC – (5f 2 7s 2 ), and UC + (5f 2 7s 0 ). Furthermore, the behavior of the UN and UC sequence of molecules and anions differs from the corresponding sequences for UO and UF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗