Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ORBIT CALCULATION”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Structural and Electronic Complexities of a Sulfur‐Bridged Di‐Iron Complex Composed of Mono‐ and Di‐Nitrosyl Units

The delocalized, thermodynamically stable cation, [(N 2 S 2 )Fe(NO)•Fe(NO) 2 ] + , an adduct of mono-nitrosyl and dinitrosyl iron units, is analyzed to address the unusual stability of the sulfur-bridged diiron complex in its three overall redox levels, +, 0, and −. X-ray diffraction and myriad spectroscopic techniques probe products of sequential electron uptake in the corresponding neutral and anionic species. Conundrums include unified blueshifts of the overall 3-band, ν(NO), pattern with added electrons. One-electron reduction changes the anti-ferromagnetically coupled, S = 0, cationic diiron species to the neutral analog, S = ½, with unpaired spin mainly localized on the MNIU, which decreases its ∠Fe–N–O angle by 10 degrees in response to the extra electron density. Subsequent reduction to the anionic species, S = 1, involves a major geometric change at the MNIU, which moves the Fe in {Fe(NO)} 8 out of the N 2 S 2 plane. Site-specific 15 N labeling of nitrosyl in the MNIU confirms the IR analysis and shows rapid NO exchange between the MNIU/DNIU (mono-nitrosyl iron unit/dinitrosyl iron unit) pairs during its synthesis at RT. Mössbauer spectroscopy, S K-edge XAS, and molecular orbital calculations confirm the ability of NO and the versatility of sulfur bridges to buffer and distribute electrons, a key to their major importance in metalloenzymes.

DNIC↗
DOE OSTI2024

Secondary electron emission measurements from imidazolium-based ionic liquids

The electron-induced secondary electron emission (SEE) yields of imidazolium-based ionic liquids are presented for primary electron beam energies between 30 and 1000 eV. These results are important for understanding plasma synthesis of nanoparticles in plasma discharges with an ionic liquid electrode. Due to their low vapor pressure and high conductivity, ionic liquids can produce metal nanoparticles in low-pressure plasmas through reduction of dissolved metal salts. In this work, the low vapor pressure of ionic liquids is exploited to directly measure SEE yields by bombarding the liquid with electrons and measuring the resulting currents. The ionic liquids studied are [BMIM][Ac], [EMIM][Ac], and [BMIM][BF 4 ]. The SEE yields vary significantly over the energy range, with maximum yields of around 2 at 200 eV for [BMIM][Ac] and [EMIM][Ac], and 1.8 at 250 eV for [BMIM][BF 4 ]. Molecular orbital calculations indicate that the acetate anion is the likely electron donor for [BMIM][Ac] and [EMIM][Ac], while in [BMIM][BF 4 ], the electrons likely originate from the [BMIM] + cation. The differences in SEE yields are attributed to varying ionization potentials and molecular structures of the ionic liquids. These findings are essential for accurate modeling of plasma discharges and understanding SEE mechanisms in ionic liquids.

ionic liquids↗
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 OSTI2026

Disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y alloys

This dataset contains input and output files from density functional theory (DFT) simulations used to study the disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y (0 ≤ x + y ≤ 1) alloys and their binary end members MnTe, MnSb, and MnBi. The alloys adopt the hexagonal NiAs-type (nickeline) structure and span ternary (MnTexSb1-x, MnTexBi1-x, MnBixSb1-x), and quaternary compositions across the full MnTe–MnSb–MnBi composition triangle. For each alloy composition, the dataset provides DFT calculations in three magnetic configurations: A-type antiferromagnetic (AFM), C-type AFM, and ferromagnetic (FM). Every magnetic configuration folder contains the fully relaxed crystal structure (CONTCAR), VASP input parameters (INCAR), and the main VASP output file (OUTCAR), from which total electronic energies, Mn magnetic moments, lattice parameters, and percent volume changes between magnetic states are extracted. These data are used to construct compositional phase diagrams, evaluate thermodynamic stability (formability), and map magnetoelastic responses across the alloy space. For A-type AFM and FM configurations, additional data are provided as follows: (i) FORCE_CONSTANTS and thermal_properties.yaml files at the top level of A-type_AFM/ and FM/ folders — present only for compositions marked with an asterisk (*) in Table I of the main text. These are derived from Phonopy finite-displacement calculations on full disordered 128-atom supercells and provide vibrational free energy, entropy (Svib)contribution from explicit disorder calculations. (Table I of the associated main manuscript) (ii) A VCA/ subfolder within A-type_AFM/ and FM/, containing FORCE_CONSTANTS and thermal_properties.yaml from Virtual Crystal Approximation phonon calculations (without spin-orbit coupling). VCA data are available for all compositions and are used to estimate vibrational contributions to the Gibbs free energy across the full composition space. (iii) A SOC/ subfolder containing CONTCAR, INCAR, and OUTCAR from spin-orbit coupling calculations, providing relativistic corrections to electronic energies and lattice parameters (Tables S2–S3 of the SM, and Table I of the main manuscript). (iv) A SOC/VCA/ subfolder containing FORCE_CONSTANTS and thermal_properties.yaml from VCA phonon calculations performed within the SOC framework, combining relativistic and vibrational thermodynamic corrections. The computed properties are used to map the AFM–FM magnetic crossover near MnTe0.75Sb0.25, demonstrate disorder- and spin-induced phonon broadening, identify a semiconductor-to-metal crossover, and quantify the pronounced magnetoelastic volume response near the magnetic phase boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Photoluminescence of a Uranium(IV) Alkoxide Complex

In this report, we describe the photoluminescence of a homoleptic uranium(IV) alkoxide complex. Excitation of [Li(THF)] 2 [U IV (O t Bu) 6 ] leads to the first example of photoluminescence from a well-defined actinide complex originating from an f–f excitation, supported by second order multiconfigurational electronic structure calculations including spin–orbit coupling. These calculations show strong spin–orbit coupling between the excited triplet and singlet states for the 5f-orbital manifold, which leads to a long-lived excited state lifetime of 0.85 s at low temperature. The photophysical properties of homoleptic uranium(V) and uranium(VI) tertbutoxide complexes are also presented; we find that oxidation of the uranium(IV) alkoxide results in quenching of luminescence in [Li(THF)][U V (O t Bu) 6 ] and [U VI (O t Bu) 6 ]. This is attributed to competing ligand to metal charge transfer absorption processes shifted to lower energy upon oxidation of the actinide center, which mask the relevant f–f transitions in the visible region of the electronic absorption spectrum.

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

Investigation into Scalable and Detection-Enhanced Satellite Conjunction Assessment

Imaging opportunities (viewable conjunctions) of Resident Space Objects (RSOs) by satellites are not continuously discovered. We propose to continuously produce and report viewable conjunctions among objects in orbit. Viewable conjunctions are events in space and time when a satellite may favorably view a Resident Space Object (RSO). Favorability is defined by a set of constraints, e.g., solar illumination, distance between observer and target, orbital location for viewable event. Computing viewable conjunctions requires calculation of orbital propagation while considering constraints based on the state vectors of position, velocity, with covariance for both satellite and RSO. We propose two parallel lanes of effort: acceleration and research. The objective of acceleration is to avoid missed opportunities and reduce latency for satellite maneuver requests through continuous prediction and reporting of viewable conjunctions. The effort will begin by deploying currently available software on dedicated systems and continue with optimizing the code for high performance computing hardware. The research lane aims to expand RSO inspection and modeling capabilities. Among our current research ideas are spectral characterization of RSO materials and planning multiple observations to recover RSO 3D form. Computing resources at Oak Ridge National Laboratory (ORNL) are available for the acceleration work. Laika, Maxar conjunction prediction dashboard software, and Bluesim, Maxar orbital propagation software, are expected to be the first software in the acceleration lane. Laike and Bluesim are to be provided by the sponsor, and output will be made accessible through its dashboard. Deliverables will follow a gated schedule to the sponsor. ORNL will provide progressively more robust viewable conjunction assessments from both modelled and actual ephemerides.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Overset-Grid Method with Smooth Orbital Partitioning for Molecular Scattering Calculations

To solve molecular photoionization and electron scattering problems, we use an overset-grid representation of electronic continuum functions, which has an extended central spherical grid that overlaps small spherical grids (subgrids) centered on each atom of a polyatomic molecule. Here, in this work, we present an improved algorithm that smoothly partitions the total wave function between the central grid and the atomic subgrids. The smooth partitioning allows one to use approximately one-fourth the number of partial waves on the central grid compared to our previous implementation with switching functions. The resulting numerical method for treating electron scattering and photoionization of polyatomic molecules combines the accuracy and flexibility of pure numerical grid representations with the rapid convergence of hybrid combinations of atom-centered basis-set expansions and grid methods. The overset-grid representation is implemented using the complex Kohn variational principle for scattering and photoionization amplitudes. The faster convergence with respect to the number of central grid partial waves is demonstrated and accuracy is verified by comparisons with the previous implementation and with far more computationally demanding single-center numerical expansions in electron-molecule scattering and photoionization calculations on the neon dimer (Ne 2 ) system, carbon tetrafluoride (CF 4 ) molecule, and the pyridine (C 5 H 5 N) molecule in the static-exchange approximation.

Molecules↗
DOE OSTI2025

Solving the “Coloring Problem” in InPd 3– x Ag x ( x = 0–0.7) by Phase Diagrams Modeling and Diffraction Experiments

Here, a series of InPd 3–x Ag x (x = 0–1) compositions were synthesized by conventional high-temperature synthesis, and as-synthesized samples were characterized by powder X-ray diffraction experiments. Up to x = 0.7, InPd 3–x Ag x adopts the ternary substitutional variant of the InPd 3 structure (TiAl 3 -type), when x > 0.7, elemental Ag starts to segregate along with the main phase. Accurate structural characterization in InPd 3–x Ag x faces a critical challenge due to the narrow X-ray scattering contrast among constituents In, Pd, and Ag and nearly identical neutron scattering lengths of Pd and Ag. To overcome this “coloring problem”, a combination of calculation of phase diagrams modeling (CALPHAD) and diffraction techniques (X-ray and neutron) was employed. In the compositional range 0 ≤ x ≤ 0.7, InPd 3–x Ag x presents a ternary variant of the TiAl 3 -type structure, where Ag atoms selectively substitute one (the 2b Wyckoff site) of the two Pd sites in InPd 3 . Notably, in contrast to the isologous InPd 3–x Cu x (x = 0–1) system, Ag substitution does not form an ordered VRh 2 Sn-type structure at the limiting composition. The distinct site preference in InPd 3–x Ag x is elucidated by charge population analysis, electronic structure calculations, and orbital-resolved chemical bonding investigations, and the extent of substitution is supported by formation free energy calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Correlation-driven topological band inversion in VSe 2

Recent spectroscopic studies have uncovered topological surface states and band inversion in 1⁢T−VSe 2 , positioning this material at the intersection of correlated electron physics and nontrivial band topology. While previous interpretations attribute these features to surface strain, the microscopic origin of the topological band structure remains unresolved. Here, we present an alternative explanation based on electronic correlations, showing that a negative effective Hubbard interaction (𝑈 eff < 0) applied to the Se 4p orbitals can reproduce the experimentally observed band inversion at both the Γ and 𝑀 points. Using density functional theory (DFT) calculations with orbital-selective interactions, we demonstrate that this approach naturally gives rise to topological surface states without invoking structural distortions. In conclusion, our results highlight the crucial role of ligand orbital correlations in shaping band topology and provide a novel framework for understanding and engineering topological phases in chalcogenide-based quantum materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Reduced fast-ion transport calculations of m = n = 1 fishbone-like instabilities in MAST-U

Fast-ion transport associated with an m = n = 1 fishbone-like burst in MAST-U discharge 47128 is investigated using a reduced guiding-center-based transport model (ORBIT-Kick) constrained by multi-diagnostic measurements. The two-dimensional beam-emission spectroscopy system provides measurements of the core poloidal mode structure and fluctuation amplitude, while EFIT++ reconstructions constrained by the motional Stark effect diagnostic indicate a flat q-profile with q 0 > 1⁠, indicating the absence of a resonant q = 1 surface and supporting a pressure-driven infernal-mode interpretation. Analytic m = n = 1 displacement profiles consistent with the measured core mode structure and equilibrium constraints are used as the mode structure inputs to ORBIT-Kick. The calculations show that the dominant resonances occur between the mode and co-passing fast ions, producing redistribution localized near the magnetic axis. Synthetic neutron camera signals from TRANSP-Kick recover up to 90% of the experimentally observed neutron deficit at the time of peak mode amplitude, indicating that the measured m = n = 1 mode is a dominant contributor to core fast-ion transport. However, the synthetic neutron signals recover rapidly, whereas the measured neutron emission continues to decrease after the peak amplitude. In conclusion, the remaining discrepancy may arise from contributions not included in the present single-harmonic model, including higher-m and higher-n harmonics, multi-harmonic interactions, and additional transport mechanisms, motivating future diagnostic development and modeling efforts to resolve and incorporate these additional contributions.

Wong, Henry H. [University of California, Los Ange↗
DOE OSTI2025

The Electronic Structure of Zirconium and Hafnium Monochalcogenides

High-level ab initio CCSD(T) and spin–orbit icMRCI+Q calculations were used to predict potential energy curves (PECs) for the lowest-lying states of ZrO, ZrS, HfO, and HfS. The prediction of the ground state is basis set dependent at the icMRCI+Q level for ZrO and ZrS due to the small singlet–triplet splitting between the lowest 1 Σ + and 3 Δ states. CCSD(T) with a spin orbit correction predicted the 1 Σ + ground state in agreement with experiment. New all-electron basis sets were developed for Hf to improve the results over those predicted by use of effective core potentials (ECPs) that subsume the 4f electrons into the definition of the core. The use of the new DK-4f basis sets rather than ECPs became more important for HfO and HfS where there is a lack of a good core–valence separation. icMRCI+Q, CCSD(T), and DFT calculations for the spectroscopic parameters of ZrO, ZrS, HfO, and HfS were benchmarked with available experimental data. Bond dissociation energies (BDEs) of these four systems were calculated at the Feller–Peterson–Dixon (FPD) level to be 762.1 (ZrO), 543.5 (ZrS), 803.8 (HfO), and 575.1 kJ/mol (HfS), in excellent agreement with experiment. The HfS BDE was remeasured using the R3PI method, providing an updated experimental measurement of D 0 (HfS) = 5.978 ± 0.002 eV = 576.8 ± 0.2 kJ/mol. This experimental value, combined with experimental measurements of the ionization energies of Hf and HfS, gives the cationic BDE of D 0 (Hf + -S) = 5.124 ± 0.002 eV = 494.4 ± 0.2 kJ/mol.

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

First-principles insights into Si substitution effects in Sm 2 (Fe,Si) 17 C x magnet

The partial substitution of Fe by Si enhances the phase stability of Sm 2 Fe 17 C x magnets with x > 1.0. Here, we elucidate the Si-substitution scheme and its impact on phase stability and magnetic properties in Sm 2 (Fe,Si) 17 C 3 from first-principles calculations and chemical bonding analysis. The calculated substitution energies for Si at various Fe sites are negative, indicating improved phase stability. Si preferentially substitutes Fe atoms at the 9d site in Sm 2 (Fe,Si) 17 C 3 while it tends to enter the Fe 18h site in Sm 2 (Fe,Si) 17 . This difference in site preference is attributed to the distinct chemical environments surrounding the Fe (Si) sites in the two compounds. Si substitution favors the formation of Sm–Si bonds while minimizing the Si–C and Si–Si interactions. Crystal orbital Hamilton populations and crystal orbital bond index calculations indicate that the partial replacement of Fe with Si strengthens the chemical bonding of Sm–Fe 3 (18f) and Sm–Fe 4 (18h) and improves overall phase stability in Sm 2 (Fe,Si) 17 C 3 . Beyond the dilution effect, Si substitution also reduces the magnetic moments of neighboring Fe atoms, a phenomenon linked to the strong Fe–Si bonding. These findings highlight the dual role of Si in modifying both the structural and magnetic characteristics of Sm 2 Fe 17 -based magnetic compounds.

Chemical bonding↗
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 OSTI2024

Self-consistent quasiparticle G W and hybrid functional calculations for Al/InAs/Al heterojunctions: Band offset and spin-orbit coupling effects

The electronic structure of surfaces and interfaces plays a key role in the properties of quantum devices. Here, we study the electronic structure of realistic Al/InAs/Al heterojunctions using a combination of density functional theory with hybrid functionals and state-of-the-art quasiparticle G W (QSGW) calculations. We find a good agreement between QSGW calculations and hybrid functional calculations, which themselves compare favorably well with angle-resolved photoemission spectroscopy experiments. Our paper confirms the need for well-controlled quality of the interfaces to obtain the needed properties of InAs/Al heterojunctions. A detailed analysis of the effects of spin-orbit coupling on the spin splitting of the electronic states shows a linear scaling in k space, related to the two-dimensional nature of some interface states. The good agreement by QSGW and hybrid functional calculations opens the door towards trustable use of an effective approximation to QSGW for studying very large heterojunctions. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2025

Elastic dijet production in electron scattering on a longitudinally polarized proton at small x : A portal to orbital angular momentum distributions

We calculate the elastic production of dijets from electron collisions with a longitudinally polarized proton target at small values of the Bjorken x variable. Building on the pioneering proposals of Hatta [Gluon orbital angular momentum at small- x , .] and Bhattacharya [Signature of the gluon orbital angular momentum, ., Probing the quark orbital angular momentum at electron-ion colliders using exclusive π 0 production, ., and Exploring orbital angular momentum and spin-orbit correlation for gluons at the Electron-Ion Collider, .] for measuring the quark and gluon orbital angular momentum (OAM) distributions, our focus is on both the longitudinal double spin asymmetry (DSA) and longitudinal single spin asymmetry (SSA). We compute the numerators of these asymmetries in the small- x formalism of the light-cone operator treatment. Utilizing the small- x expressions for the OAM distributions derived earlier in Kovchegov and Manley [Orbital angular momentum at small x revisited, .], we demonstrate that the DSA provides a robust probe for both the quark and gluon OAM distributions within the proton. In contrast, we find that while the SSA is also sensitive to the OAM distributions, extraction of the latter from the SSA would require new developments in small- x theory and phenomenology, and is probably not feasible at this point in time. These findings highlight the potential of DSA measurements in elastic dijet production at the future Electron-Ion Collider to provide the first-ever direct access to the quark and gluon OAM distributions at small x , paving the way for new insights into the proton spin puzzle. Published by the American Physical Society 2025

Kovchegov, Yuri V. (ORCID:0000000169909173)↗
DOE OSTI2026

A Cautionary Tale: Failure of the Valence CASSCF to Describe the Hallmark of Hydrogen Bonding

Valence CASSCF (vCAS) calculations for the hydrogen-bonded (H2O)2, (HF)2, and HF–H2O dimers fail to predict the lengthening of the hydrogen donor bond and the corresponding red shift in the donor bond vibrational stretching frequency. Analysis of the active orbitals in the vCAS calculation reveals that one of the weakly occupied orbitals in each dimer is used to correlate the electrons in the acceptor lone pairs rather than those in the donor bonding orbitals. This has significant implications for subsequent treatments of the dynamical correlation that are based on the vCAS wavefunction.

Nguyen, Long H.↗
DOE OSTI2026

Bond Dissociation Energies of the Actinide Halides AnX, An = Ac–Lr and X = F–I, Utilizing Relativistic Composite Coupled Cluster Approaches

Bond dissociation energies (BDEs) have been calculated for the set of actinide halides AnX with An=Ac, Pa, and Np-Lr and X=F-I. Two composite thermochemistry methods based on the Feller-Peterson-Dixon (FPD) approach have been utilized, one involving spinor-based relativistic CCSD(T) calculations where spin-orbit (SO) was included at the orbital level and another using scalar relativistic CCSD(T) with a posteriori SO contributions based on 2-component multireference configuration interaction calculations. The method that was chosen for a given actinide halide was based on which representation yielded the best single determinant reference determinant for the coupled cluster calculation. The spinor-based method was chosen for all cases except for AmX, CmX, and BkX. Both composite approaches included contributions accounting for basis set truncation, outer-core correlation, the Gaunt interaction, and QED. The scalar FPD results, as well as the spinor-based calculations for AcF, also included higher order electron correlation up through CCSDT(Q). In addition to BDEs, CCSD(T) equilibrium bond lengths, harmonic frequencies, and vibrational anharmonicity constants are reported for all species. Last, the FPD BDEs for the fluorides were used to confirm the trend across the actinide series previously predicted by Gibson using bonding models based atomic promotion energies that provide a single 6d electron for bonding. In particular the local minimum in the BDEs at AmF is confirmed in the present calculations. Furthermore, the BDEs for LrX are predicted to be slightly larger than those of AcX, making them the largest in the actinide halide series.

Actinides↗
DOE OSTI2025

Toward an Understanding of Linear Scaling Relations through Energy Decomposition Analysis

The discovery of linear scaling relations has fundamentally changed the field of heterogeneous catalysis. The scaling relations have been rationalized based on the d-band theory, specifically a separation of sp and d electron contributions to adsorption energies. Within the framework of energy decomposition analysis, a full understanding of such a separation would require one to further break down the adsorption energy into distinct energy components such as electrostatics, polarization, charge transfer, and van der Waals interactions, and to examine the sp and d contributions to each of them. As a step in this direction, we analyzed the interaction energy between CH x (x = 1–4) adsorbates and fcc(100) transition metal surfaces (M = Cu, Ag, Au, Rh, and Pt), with the surfaces represented both as slabs in plane-wave density functional theory (pw-DFT) calculations and as atomic clusters in atomic-orbital basis density functional theory (ao-DFT) calculations. Through an absolutely localized molecular orbital (ALMO) based energy decomposition analysis of the ao-DFT adsorption energy, each of the interaction energy components (electrostatics, polarization, van der Waals, and charge transfer) was found to follow its own scaling relations, with an intricate interplay among these energy components yielding the overall scaling relations for the total adsorption energies. Using the recently introduced ALMO-based polarization and charge-transfer analysis schemes, we further dissected polarization into metal surface and adsorbate contributions, and charge transfer into metal → adsorbate and adsorbate → metal contributions. The contributions from the sp and d electrons of the metal to these terms were further quantified, and the dominant role of the metal d electrons was reaffirmed. These results shed light on how CHx adsorbates interact with metal surfaces and further reveal the physical origin of the scaling relations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗