Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ORBITAL ELEMENT”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Comparing the fluorination of Pb+ and Po+: insights into p-block relativistic effects for superheavy elements

It is expected that the chemical properties of the heaviest of the superheavy elements (SHEs, = 113 – 118) do not align with what is suggested by their current positions on the periodic table. Specifically, the onset of significant relativistic effects, including increased spin–orbit splitting of the p-orbitals may lead to enhanced stability of low-oxidation states. Notably, it is predicted that flerovium ( = 114) may exhibit pseudo-noble gas behavior from a electron configuration that acts as an ‘inert pair’. Even though it is expected that the – splitting becomes pronounced toward the end of the p-block’s sixth row, there is currently limited experimental evidence to confirm its extent or impact. Here, the production of gas-phase lead ( = 82) and polonium ( = 84) fluoride cations (PbF and PoF ) were compared to elucidate differences in accessible oxidation states. The PoF and PbF species were produced and identified with the FIONA spectrometer at the Lawrence Berkeley National Laboratory 88-Inch Cylctron Facility. Polonium showed notably different fluorination chemistry than lead, producing PoF and PoF as primary products compared to lead’s PbF and PbF . The distribution of PoF products observed offer insights as to the role of spin–orbit splitting for polonium. Similar studies of superheavy elements would elucidate the accessibility of their low-oxidation states as well as to inform chemical predictions for eighth-row elements not yet discovered.

Pore, Jennifer L [Nuclear Science Division, Lawren↗
DOE OSTI2025

4f-orbital covalency enables a single-crystal-to-single-crystal ring-opening isomerization in a CeIV–cyclopropenyl complex

Metal–ligand bonding interactions for f-element compounds are typically highly polarized with only minor covalent character. Whereas the 5d/6d orbitals are known to be chemically accessible for dative bonding, recent quantum chemical and spectroscopic analyses have indicated appreciable 4f/5f-orbital involvement in certain metal–ligand bonds. However, 4f-orbital covalency has not been compellingly linked to distinctive modes of chemical reactivity via rigorous comparative study and mechanistic investigation. Here a series of MIV–cyclopropenyl complexes (M = Ti, Zr, Ce, Hf, Th) are described, wherein the cerium congener exhibits a 4f-covalent Ce=Cα interaction, causing a ring-opening isomerization reaction through a single-crystal-to-single-crystal transformation. The results provide evidence for 4f-orbital covalency by demonstrating its expression in the reactivity of an f-element complex within an isostructural series of tetravalent d- and f-block metal complexes. They also provide new directions for the study of orbital covalency effects of molecular compounds in solid-state chemical transformations.

Vincenzini, Brett D↗
DOE OSTI2024

Efficient Parameterization of Density Functional Tight-Binding for 5 f -Elements: A Th–O Case Study

Density functional tight binding (DFTB) models for f-element species are challenging to parametrize owing to the large number of adjustable parameters. The explicit optimization of the terms entering the semiempirical DFTB Hamiltonian related to f orbitals is crucial to generating a reliable parametrization for f-block elements, because they play import roles in bonding interactions. However, since the number of parameters grows quadratically with the number of orbitals, the computational cost for parameter optimization is much more expensive for the f-elements than for the main group elements. In this work we present a set of efficient approaches for mitigating the hurdle imposed by the large size of the parameter space. A novel group-by-orbital correction functions for two-center bond integrals was developed. With this approach the number of parameters is reduced, and it grows linearly with the number of elements, maintaining the accuracy and the number of parameters, in the case of f elements, by more than 40%. The parameter optimization step was accelerated by means of the mini-batch BFGS method. This method allows parameter optimizations with much larger training sets than other single batch methods. A stochastic optimizer was employed that helped overcome shallow local minima in the objective function. The proposed algorithm was used to parametrize the DFTB Hamiltonian for the Th–O system, which was subsequently applied to the study of ThO 2 nanoparticles. The training set consisted of 6322 unique structures, which is barely feasible with conventional optimization methods. The optimized parameter set, LANL-ThO, displays good agreement with DFT-calculated properties such as energies, forces, and structures for both clusters and bulk ThO 2 . Benefiting from the fewer number of parameters and lower computational costs for objective function evaluations, this new approach shows its potential applications in DFTB parametrization for elements with high angular momentum, which present a challenge to conventional methods.

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

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗
DOE OSTI2026

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗
DOE OSTI2026

One-Step Relativistic Driven Similarity Renormalization Group Multireference Perturbation Theory

We present an efficient implementation of a one-step relativistic second-order multireference perturbation theory based on the multireference driven similarity renormalization group (MR-DSRG) using the exact two-component (X2C) Hamiltonian, which we denote X2C-DSRG-MRPT2. We show that the X2C-DSRG-MRPT2 method can accurately capture spin–orbit coupling (SOC) effects in the electronic structure of strongly correlated systems containing elements across the periodic table. We further demonstrate that the X2C-DSRG-MRPT2 method, through its variational treatment of SOC effects, can yield spin–orbit splittings with mean absolute percentage errors consistently below 7% with respect to experimental values for systems containing up to sixth row elements. With its modest computational scaling (fourth power in system size for the perturbative step) and high accuracy, X2C-DSRG-MRPT2 provides a promising avenue for the routine treatment of relativistic effects in strongly correlated molecular systems.

Hamiltonians↗
DOE OSTI2025

The balance of orbital overlap and orbital energy in the activation of methane by actinide cations: insights from inductively coupled plasma tandem mass spectrometry

The actinides present a unique challenge to chemical theory. The classical view of covalent bonding is driven by the extent of spatial overlap of valence orbitals. Modern theory has expanded assessments of covalency to include considerations of orbital energy degeneracy to assess orbital energy mixing between metal and ligand valence orbitals. Actinide–ligand (An–L) bonding has more recently been described as a balance between orbital overlap and orbital energy mixing, where 5f and L valence orbital overlap decreases while energy mixing between An 5f and L valence orbitals increases across the series. To test these existing views, we employed inductively coupled plasma tandem mass spectrometry to examine the kinetic energy dependences of reactions of actinide cations, Th + –Am + , with methane. Further, this is the first experimental report of the energy dependences of methane activation reactions involving the cations of Pa, Np, Pu, and Am and the first experimental determination of transuranic An + –D, An + –CD 2 , An + –CD 3 , and An + –CD bond dissociation energies. The correlation of the measured An + –CD 2 bond energies with E p (6d 2 ) indicates that An + 6d orbitals are the dominant contributors in the An + –CD 2 bonds. Close examination of the relative reactivities of An + offers additional support that the balance of classical and modern views of molecular bonding may lie between Np + and Pu + and that the increased reactivity of Th + –Np + may be attributed to the increased spatial extension of the 5f orbitals whereas covalent An + bond formation may be more driven by the decreasing energies of the 5f orbitals across the actinide series.

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

Electronic structure simulations in the cloud computing environment

The transformative impact of modern computational paradigms and technologies, such as high-performance computing, quantum computing, and cloud computing, has opened up profound new opportunities for scientific simulations. Scalable computational chemistry is one beneficiary of this technological progress. The main focus of this paper is on the performance of various quantum chemical formulations, ranging from low-order methods to high-accuracy approaches, implemented in different computational chemistry packages, such as NWChem, NWChemEx, SPEC, ExaChem, and FLOSIC codes on the Azure Quantum Element (AQE) Microsoft cloud services. We pay particular attention to the intricate workflows for performing composite chemistry simulations, associated data curation, and mechanisms for accuracy assessment, as defined by the enabling cloud Computational Chemistry as a Service (CCaaS). Our focus also extends to Arrows' automated workflow for high throughput simulations. Finally, we provide a perspective on the role of cloud computing in supporting the mission of leadership computational facilities (LCFs).

computational chemistry, electronic structure, Clo↗
DOE OSTI2025

Spinor-based coupled cluster thermochemistry: RgF and CnF 0/+ as compared to AuF and HgF 0/+

A relativistic coupled cluster approach that includes spin-orbit variationally in the eXact 2-component (X2C) approximation with the inclusion of the Gaunt interaction (X2Cg) was used to probe the thermochemistry and ground state spectroscopic constants of RgF, CnF, and CnF + . Utilizing large sequences of correlation consistent basis sets at the X2Cg-CCSD(T) level of theory, this work reports 0 K bond dissociation energies (BDEs) of AuF, RgF, HgF 0/+ , and CnF 0/+ , as well as ionization energies of Au, Rg, Hg, Cn, HgF, and CnF.

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

One Million Open-source Cislunar Orbits

Cislunar space, encompassing the region from geosynchronous orbit to beyond the Moon, is poised to become a cornerstone for future exploration, scientific discovery, and national security. Missions in this region, spanning durations from weeks to decades, require robust infrastructure and reliable transit capabilities. The complex gravitational influences of the Moon, Sun, and planets, along with thermal radiation from Earth and the Sun, lead to significant trajectory deviations, resulting in kilometer-scale errors within days. Leveraging the high-performance computing resources at Lawrence Livermore National Laboratory (LLNL), we have simulated one million high-fidelity cislunar trajectories, now publicly available via LLNL’s Green Data Oasis and the Unified Data Library. Generated using the open-source Space Situational Awareness Python package, these trajectories match the precision of commercial tools such as AGI’s Systems Tool Kit and NASA’s General Mission Analysis Tool. This data set is a valuable resource for reference, statistical analysis of cislunar orbit populations, and training machine learning models for rapid orbit classification with minimal observational input. Preliminary analysis reveals stable bands in Keplerian element space, particularly around five geosynchronous radii across a range of inclinations and eccentricities. Beyond this threshold, the Moon’s influence disrupts most unassisted orbits, though co-orbiting L4/L5 Lunar Trojans persist throughout the six-year simulation.

Astronomy and AstroPhysics↗
DOE OSTI2024

A new computational framework for spinor-based relativistic exact two-component calculations using contracted basis functions

Here, a new computational framework for spinor-based relativistic exact two-component (X2C) calculations is developed using contracted basis sets with a spin–orbit contraction scheme. Generally contracted, j-adapted basis sets of p-block elements using primitive functions in the correlation-consistent basis sets are constructed for the X2C Hamiltonian with atomic mean-field spin–orbit integrals (the X2CAMF scheme). The contraction coefficients are taken from atomic X2CAMF Hartree–Fock spinors, thereby following the simple concept of a linear combination of atomic orbitals. Benchmark calculations of spin–orbit splittings, equilibrium bond lengths, and harmonic vibrational frequencies demonstrate the accuracy and efficacy of the j-adapted spin–orbit contraction scheme.

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

Understanding the optoelectronic properties of doped 2D organic-inorganic halide perovskite quantum wells: towards efficient quantum well IR photodetectors

Metal halide perovskite (MHP) multiple quantum wells which consist of multilayers of alternate organic and inorganic layers exhibit large exciton binding energies due to the dielectric confinement between the inorganic and organic layers. These naturally formed multiple quantum wells have strong spin-orbit coupling (SOC) due to the presence of heavy elements in their crystal structures. Although the fundamental properties of 2D MHPs are far from being entirely understood, it is widely accepted that their band edge absorption coefficient results from strong exciton interactions. However, studies demonstrating how different exciton interactions and doping effects influence electronic traps and disorder on the band edge absorption coefficient of 2D MHPs have not been demonstrated. Understanding these interactions in MHPs will allow us to access low energy optical transitions for the fabrication of solution processable short-to-mid-wavelength IR photodetectors (1 – 8 μm). Moreover, upon doping, it is possible to move the Fermi energy into the conduction band (CB) to favorably promote the transport of charges in a working device. Herein, we study the development of 2D MHPs having strong SOC, high carrier mobility, and tunable quantum well structures. Our studies shed light on the design and modulation of fundamental physical phenomena by carefully elucidating the role of dopants (n-type and p-type), exciton heterogeneity, orientation, structure, and bias stress effects on the performance of MHPs as potential IR photodetectors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

SSAPy - Space Situational Awareness for Python

SSAPy is a fast and flexible orbit modeling and analysis tool for orbits spanning from low-Earth into the cislunar regime. Orbits can be flexibly specified from common input formats such as Keplerian elements or two-line element (TLE) data files. SSAPy allows users to model satellites and specify parameters such as satellite area, mass, and drag coefficients. SSAPy includes a customizable force-propagation with a range of Earth, Lunar, radiation, atmospheric, and maneuvering models. SSAPy makes use of various community integration methods and can calculate time-evolved orbital quantities, including satellite magnitudes and state vectors. Users can specify various space- and ground-based observation models with support for multiple coordinate and reference frames. SSAPy also supports orbit analysis and propagation methods such as multiple hypothesis tracking and has built-in uncertainty quantification. The majority of SSAPy’s methods are vectorized and parallelizable, allowing for effective use of high-performance computer (HPC) systems. Finally, SSAPy has plotting functionality, allowing users to visualize orbits and trajectories. Examples are shown in Figure 1 and Figure 2.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Lattice Symmetry-Guided Charge Transport in 2D Supramolecular Polymers Promotes Triplet Formation

Singlet-to-triplet intersystem crossing (ISC) in organic molecules is intimately connected with their geometries: by modifying the molecular shape, symmetry selection rules pertaining to spin-orbit coupling can be partially relieved, leading to extra matrix elements for increased ISC. As an analog to this molecular design concept, the study finds that the lattice symmetry of supramolecular polymers also defines their triplet formation efficiencies. A supramolecular polymer self-assembled from weakly interacting molecules is considered. Its 2D oblique unit cell effectively renders it as a coplanar array of 1D molecular columns weakly bound to each other. Using momentum-resolved photoluminescence imaging in combination with Monte Carlo simulations, the study found that photogenerated charge carriers in the supramolecular polymer predominantly recombine as spin-uncorrelated carrier pairs through inter-column charge transfer states. This lattice-defined recombination pathway leads to a substantial triplet formation efficiency (≈60%) in the supramolecular polymer. These findings suggest that lattice symmetry of micro-/macroscopic structures relying on intermolecular interactions can be strategized for controlled triplet formation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

New precious metal containing normal spinels: LiRhRu 1-x Ir x O 4 , LiFeIr 1-x Ru x O 4 , and LiCoIr 1-x Ru x O 4

Large spin-orbit-coupled cations in geometrically frustrated crystal structures have the most suitable setting for exploring novel exotic states of matter. Spinel oxides (AM 2 O 4 ) are well-known examples of geometrically frustrated systems. In this study, we report for the first time the synthesis of compositions LiRhRu 1-x Ir x O 4 (x = 0–0.5), LiFeIr 1-x Ru x O 4 (x = 0–0.5), and LiCoIr 1-x Ru x O 4 (x = 0–0.3) containing precious metal cations on edge-sharing octahedral M-sites, and systematically investigate their magnetic and electrical properties. 57 Fe Mössbauer spectroscopy revealed that iron is trivalent in all LiFeIr 1-x Ru x O 4 solid solutions. Magnetic measurements indicate deviations from theoretical spin-only magnetic moment values, indicating the influence of spin-orbit coupling owing to the presence of 4d and 5d block elements. The LiFeIr 1-x Ru x O 4 series shows spin-glass-like freezing behavior with T g ≈ 20 K, and a small frustration index (f ≈ 1-2), indicating that the frustration originates from site disorder. LiRhRu 1-x Ir x O 4 and LiCoIr 1-x Ru x O 4 exhibit strongly geometrically frustrated magnetism. Electrical resistivity measurements as a function of temperature indicate that all phases are semiconducting. Seebeck coefficient measurements show that LiRhRu 1-x Ir x O 4 and LiFeIr 1-x Ru x O 4 are p-type semiconductors with holes as the major charge carriers. A sign reversal of the Seebeck coefficient indicates both holes and electrons as carriers for LiCoIr 1-x RuxO 4 (x = 0–0.2), but only holes as major carriers for x = 0.3. Here, the Seebeck coefficient and power factor increase drastically in the LiRhRu 1-x Ir x O 4 solid solution with Ir substitution, reaching a maximum of ≈ +125 μV/K and ≈2.3×10 -6 W/mK 2 at ∼650 K for x = 0.5.

Charge carriers↗
DOE OSTI2025

High-pressure reduction of carbon dioxide in reactive liquid mixtures

Density-functional theory based molecular-dynamics simulations were used to investigate high-pressure chemical reactions in liquid mixtures of CO 2 with several elements (Si, Mn, and Fe) at high temperatures of 2000-3000 K. Our ab initio simulations indicate that these reactant elements can reduce CO 2 to C at high pressures (20 GPa) leading to the formation of C-C chains, with Si by far the most effective carbon-reducing agent. A combined chemical analysis using Bader charge analysis and Crystal Orbital Hamilton Population (COHP) on simulation snapshots shows that significant charge transfer from the reducing element to the C atoms creates instability in the C-O covalent bonds. COHP analysis further shows that Mn/Fe-O and Mn/Fe-C bonding interactions are weaker compared to the Si counterparts. These results further our understanding of the redox chemistry of CO 2 at conditions relevant to planetary mantle interiors and demonstrate the effectiveness of high pressure in the reduction of CO 2 directly to solid carbon.

Ab initio molecular dynamics↗
DOE OSTI2025

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↗