Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “covalency”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Americium Oxalate: An Experimental and Computational Investigation of Metal–Ligand Bonding

A novel actinide-containing coordination polymer, [Am(C 2 O 4 )(H 2 O) 3 Cl] (Am-1), has been synthesized and structurally characterized. The crystallographic analysis reveals that the structure is two-dimensional and comprised of pseudo-dimeric Am 3+ nodes that are bridged by oxalate ligands to form sheets. Each metal center is nine-coordinate, forming a distorted capped square antiprism geometry with a C 1 symmetry, and features bound oxalate, aqua, and chloro ligands. The Am 3+ –ligand bonds were probed computationally using the quantum theory of atoms in molecules nd natural localized molecular orbital approaches to investigate the underlying mechanisms and hybrid atomic orbital contributions therein. The analyses indicate that the bonds within Am-1 are predominantly ionic and the 5f shell of the Am 3+ metal centers does not add a significant covalent contribution to the bonds. Further, our bonding assessment is supported by measurements on the optical properties of Am-1 using diffuse reflectance and photoluminescence spectroscopies. The position of the principal absorption band at 507 nm ( 5 L 6' ← 7 F 0' ) is notable because it is consistent with previously reported americium oxalate complexes in solution, indicating similarities in the electronic structure and ionic bonding. Compound Am-1 is an active phosphor, featuring strong bright-blue oxalate-based luminescence with no evidence of metal-centered emission.

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

Ligand Control of Oxidation and Crystallographic Disorder in the Isolation of Hexavalent Uranium Mono-Oxo Complexes

Here, the development of high-valent transuranic chemistry requires robust methodologies to access and fully characterize reactive species. We have recently demonstrated that the reducing nature of imidophosphorane ligands supports the two-electron oxidation of U 4+ to U 6+ and established the use of this ligand to evaluate the inverse-trans-influence (ITI) in actinide metal-ligand multiple bond (MLMB) complexes. To extend this methodology and analysis to transuranic complexes, new small-scale synthetic strategies and lower-symmetry ligand derivatives are necessary to improve crys-tallinity and reduce crystallographic disorder. To this end, the synthesis of two new imidophosphorane ligands, [N=P t Bu(pip) 2 ]¯ (NPC 1 ) and [N=P t Bu(pyrr) 2 ]¯ (NPC 2 ) (pip = piperidinyl; pyrr = pyrrolidinyl) is presented, which break pseudo-C 3 axes in the tetravalent complexes, U[NPC 1 ] 4 and U[NPC 2 ] 4 . The reaction of these complexes with two-electron oxygen atom transfer reagents (N 2 O, trimethylamine N-oxide (TMAO) and 2,3:5,6-dibenzo-7-azabicyclo[2.2.1]hepta-2,5-diene (dbabhNO)) yields the U 6+ mono-oxo complexes U(O)[NPC 1 ] 4 and U(O)[NPC 2 ] 4 . This methodology is optimized for direct translation to transuranic elements. Of the two ligands, the NPC 2 framework is most suitable for facilitating detailed bonding analysis and assessment of the ITI. Theoretical evaluation of the U–(NPC) bonding confirms a substantial difference between axially and equatorially-bonded N atoms, revealing mark-edly more covalent U–N ax interactions. The U 6d+5f combined contribution for U–N ax is nearly double that of U–N eq , accounting for ITI shortening and increased bond order of the axial bond. Two distinct N-atom hybridizations in the pyrrolidine/piperidine rings are noted across the complexes, with approximate sp 2 and sp 3 configurations describing the slightly shorter P–N “planar” and slightly longer P–N “pyramidal” bonds, respectively. In all complexes, the NPC 2 ligands feature more planar N atoms than NPC 1 , in accordance with a higher electron-donating capacity of the former.

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

Synthesis and Surface Attachment of Molecular Re(I) Complexes Supported by Functionalized Bipyridyl Ligands

Eleven 2,2'-bipyridine (bpy) ligands functionalized with attachment groups for covalent immobilization on silicon surfaces were prepared. Five of the ligands feature silatrane functional groups for attachment to metal oxide coatings on the silicon surfaces, while six contain either alkene or alkyne functional groups for attachment to hydrogen-terminated silicon surfaces. The bpy ligands were coordinated to Re(CO) 5 Cl to form complexes of the type Re(bpy)(CO) 3 Cl, which are related to known catalysts for CO 2 reduction. Six of the new complexes were characterized using X-ray crystallography. As proof of principle, four molecular Re complexes were immobilized on either a thin layer of TiO 2 on silicon or hydrogen-terminated silicon. The surface-immobilized complexes were characterized using X-ray photoelectron spectroscopy, IR spectroscopy, and cyclic voltammetry (CV) in the dark and for one representative example in the light. The CO stretching frequencies of the attached complexes were similar to those of the pure molecular complexes, but the CVs were less analogous. For two of the complexes, comparison of the electrocatalytic CO 2 reduction performance showed lower CO Faradaic efficiencies for the immobilized complexes than the same complex in solution under similar conditions. In particular, a complex containing a silatrane linked to bpy with an amide linker showed poor catalytic performance and control experiments suggest that amide linkers in conjugation with a redox-active ligand are not stable under highly reducing conditions and alkyl linkers are more stable. So a conclusion of this work is that understanding the behavior of molecular Re catalysts attached to semiconducting silicon is more complicated than related complexes, which have previously been immobilized on metallic electrodes.

14 SOLAR ENERGY↗
DOE OSTI2023

Magnetic Properties of Tetravalent Pu in the Perovskites BaPuO 3 and SrPuO 3

BaPuO 3 and SrPuO 3 were synthesized, and their structures were refined in the orthorhombic space group Pbnm , a common distortion from the classic $Pm\bar{3}$m over barm cubic perovskite. Magnetic-susceptibility measurements, obtained as a function of temperature over the range of 1.8-320 K, exhibit temperature-dependent behavior, with evidence of long-range magnetic order at temperatures higher than their lanthanide and actinide analogues: BaPuO 3 below 164(1) K and SrPuO 3 below 76(1) K. Effective moments of 1.66(10)μ B for BaPuO 3 and 1.84(8)μ B for SrPuO 3 were obtained by fitting their paramagnetic susceptibilities using the Curie-Weiss law. Further, both are below the free-ion value of 2.68 μ B expected for a Pu 4+ 5 $\mid$ 4 ground level. Ab initio wave function calculations, performed at the relativistic complete active space level including spin-orbit coupling and with an embedded cluster approach that neglects interactions between Pu centers, were used to generate embedded-cluster Pu 4+ magnetic susceptibilities. The calculations agree well with experimental data at higher temperatures, providing evidence that a single-ion representation is sufficient to account for the observed paramagnetic behavior without the need to invoke charge transfer, disproportionation, strong covalent bonding, or other more complex electronic behavior.

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

A Spectroscopic and Computational Evaluation of Uranyl Oxo Engagement with Transition Metal Cations

Here, we report the synthesis and characterization of five novel Cd 2+ /UO 2 2+ heterometallic complexes that feature Cd-oxo distances ranging from 78 to 171% of the sum of the van der Waals radii for these atoms. This work marks an extension of our previously reported Pb 2+ /UO 2 2+ and Ag + /UO 2 2+ complexes, yet with much more pronounced structural and spectroscopic effects resulting from Cd-oxo interactions. We observe a major shift in the U═O symmetric stretch and significant uranyl bond length asymmetry. The ρbcp values calculated using Quantum Theory of Atoms in Molecules (QTAIM) support the asymmetry displayed in the structural data and indicate a decrease in covalent character in U═O bonds with close Cd-oxo contacts, more so than in related compounds containing Pb 2+ and Ag + . Second-order perturbation theory (SOPT) analysis reveals that O sp x → Cd s is the most significant orbital overlap and U═O bonding and antibonding orbitals also contribute to the interaction (U═O σ/π → Cd d and Cd s → U═O σ/π*). The overall stabilization energies for these interactions were lower than those in previously reported Pb 2+ cations, yet larger than related Ag + compounds. Analysis of the equatorial coordination sphere of the Cd 2+ /UO 2 2+ compounds (along with Pb 2+ /UO 2 2+ complexes) reveals that 7-coordinate uranium favors closer, stronger M n+ -oxo contacts. These results indicate that U═O bond strength tuning is possible with judicious choice of metal cations for oxo interactions and equatorial ligand coordination.

cations↗
DOE OSTI2024

Effects of Structural Constraints on Excited-State Properties in Dimeric Cu(I) Diimine Complexes

Copper(I) bis-diimine complexes have played important roles in light-activated processes that can lead to their potential applications in photocatalysis and chemical sensing. Their metal-to-ligand charge-transfer (MLCT) excited-state properties are tunable by various structural factors. Dimeric Cu(I) complexes with connecting diimine derivative ligands offer another structural tuning platform for the excited-state properties. Here, we investigate excited-state properties in two covalently connected dimeric Cu(I)'s with varying structural constraints exerted by the number of carbons in the polyethylene bridge (C0 and C4) connecting the two copper(I) diimine moieties. An interesting feature of Cu(I) diimine complexes is their ability to flatten following a photoinduced structural change. Herein, we observe larger structural constraints and more structural rearrangement required upon excitation of the longer bridged complex C4 to achieve a conformation toward a more flattened tetrahedral coordination geometry compared to the shorter bridged C0. Vibrational wavepacket analysis of these complexes further supports the effect of these structural constraints where we observe a more rapid dephasing of the C0 complex, as opposed to the C4 complex, despite similar normal mode vibrations. The experimental results were supplemented by TDDFT calculations. In conclusion, the studies provide insight into using metal-metal interactions through constraints to tune excited-state dynamics and pathways.

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

Distinguishing Desirable and Undesirable Reactions in Multicomponent Systems for Redox Activation of the Uranyl Ion

Although it has been established that covalent functionalization of the U–O bonds in the uranyl dication (UO 2 2+ ) generally requires use of strong reductants and electrophiles, little work has examined how interactions between the individual reaction components could affect final outcomes in solution. Here, the patterns of such reactivity have been studied in a UO 2 2+ -containing model system supported by a workhorse pentadentate ligand, 2,2′-[(methylimino)bis(2,1-ethanediylnitrilomethylidyne)]bis-phenol. Oxo activation and functionalization have been tested with (i) electrochemical and chemical reduction, and (ii) coordinating and noncoordinating solvents. In acetonitrile, uranyl reduction was achieved cleanly, but treatment of the reduced species with tris(pentafluorophenyl)borane (BCF) resulted in a mixture of products arising from direct electron transfer to BCF. In dichloromethane (CH 2 Cl 2 ), electrochemical reduction of uranyl was achieved cleanly, but clean chemical reactivity was inaccessible. Despite these challenges, one trinuclear and oxo-deficient uranium-containing product was crystallized from CH 2 Cl 2 solution and characterized; thus, desirable electrophilic reactivity can proceed to some degree in CH 2 Cl 2 with BCF. Computational studies were used to investigate the properties of the trinuclear uranium product and the changes that could be inducible by further reduction. Here, taken together, the reactivity patterns identified here could inform design of improved systems for actinyl oxo functionalization.

cations↗
DOE OSTI2025

Affinity for OH – Produces Four-Coordinated Zn 2+ Impurities in Hydrated Amorphous Calcium Carbonate

Using ab initio based molecular dynamics and electronic structure calculations, we show that Zn impurities in hydrated amorphous calcium carbonate (ACC) have a much lower coordination number than other divalent impurities due to covalent interactions between the 3d Zn shell and the oxygen atoms of the carbonate and water groups. Further, the local structure around Zn in ACC, including the predicted low coordination number, is confirmed by X-ray absorption spectroscopy of synthetic Zn-bearing ACC. The strong Zn–O chemical interaction leads to substantial water dissociation and slightly disrupts the hydrogen bonding network. Implications of Zn 2+ incorporation for ACC stability are discussed.

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

Actinide 5f Occupations: The Case of PuO 2

In actinide chemistry, the formal number of open-shell 5f electrons, n open , is a well-defined quantity with an integer value. The effective 5f occupation, n f , additionally takes donation and back-donation into account, generally has a non-integer value, and has varying numerical definitions. The present study explores the important distinction between n f and n open in actinide chemistry with the example of PuO2, by using electronic structure methods with a relativistic Hamiltonian in combination with experimental Pu M 5 -edge high-energy-resolution X-ray absorption and emission spectroscopic data. The total donation to the metal in PuO 2 is between 3.1 and 2.4 electrons, depending on the type of calculation, most of which is to the Pu 6d and 5f shells. The donation into 5f is sensitive to the approximations in the electronic structure model but likely amounts to 1.6/0.8 electrons when the diffuse regions of the 5f shell are included/excluded. Valence band resonant inelastic X-ray scattering experiments demonstrate that Pu 5f electron density is present in the valence band; thus, there is a clear experimental signature of covalent bonding in PuO 2 . Pu M 5 -edge and M 3 -edge high-energy-resolution X-ray absorption near-edge structures for Pu 3+ and Pu 4+ in an aqueous solution are compared to PuO 2 , showing that Pu in PuO 2 has an nf closer to Pu 4+ (aq).

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

Unveiling Unusual Reactivity of SO 2 and Unusual Type of S–X Long Bonds

A new reactivity of SO 2 to form unusual stable M−X−SO 2 (X = F, Cl, H) complexes is unveiled in this study. Moreover, a new type of S−X long bonds, which are significantly longer than traditional S−X covalent bonds, has been discovered. The P,N-ligated Ni−F complex model 1A can bind a SO 2 molecule through the new F−S long bond (2.207 Å), and a stable Ni−F−SO 2 complex 1B is generated, being exergonic by 2.2 kcal/mol. According to natural localized molecular orbital analysis, the new S−F long bond has a unique p(F) → π*(O=S=O) bonding interaction, which is shown to arise from the long S−F length. In comparison, the strength of the new F−S long bond (−2.2 kcal/mol) is found to be significantly stronger than common noncovalent interactions such as the hydrogen and halogen bond. The substituent modulations suggest that the electron-donating groups can increase the strength of new F−S bonds and enhance binding free energies ΔG bind . The scope of possible M−X complexes was explored, and various metals and X (F, Cl, and H) ligands were found to form stable M−X−SO 2 complexes. Specifically, the anionic M−X complexes display much higher ΔG bind values, ranging from −8 to −10 kcal/mol. The study paves the way for a green, recyclable, and adjustable SO 2 absorption method.

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

Reactions of U(DMSO) 8 (ClO 4 ) 4 with Terpyridine Yield Dimeric Hydrolysis Products and Induce C–C Coupling

Reactions have been carried out using UIV(DMSO) 8 (ClO 4 ) 4 with 2,2′:6′,2″-terpyridine (terpy) under nonaqueous conditions. At room temperature in acetonitrile, the combination of the U(IV) starting material with terpy resulted in a mixture containing [UO 2 (DMSO) 2 terpy][ClO 4 ] 2 ·MeCN, while increasing the water content led to the hydrolysis products [(UO 2 (DMSO)terpy) 2 (μ 2 –O)][ClO 4 ] 2 and [(UO 2 terpy) 2 (μ 2 –OH) 2 ][ClO 4 ] 2 ·MeCN·H 2 O. Performing the reaction at slightly elevated temperature with no added water led to the formation of [UO 2 sexipyridine][ClO 4 ] 2 ·MeCN. This new uranyl complex contains the hexadentate ligand 2,2′:6′,2″:6″,2″:6‴,2⁗:6⁗,2⁗′-sexipyridine, which formed in situ from the tetravalent uranium starting material, where photoexcited uranyl or in situ generated peroxide could have induced C–C coupling. Analysis of bonding in the dimeric uranyl species via quantum chemical methods revealed a small increase in covalency of the bridging oxo unit and a slightly greater stability compared to the bridging hydroxo compound, which causes a significant shift in the uranyl symmetric stretch in the Raman spectrum. Structural, spectroscopic, and computational comparisons are made across the series of compounds, providing insight into the bonding and reactivity of uranium in nonaqueous media.

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

From Serendipity to Exploration: Synthesis and Properties of PtSiSb

PtSiSb is a compound with covalent Si–Sb bonds, which are relatively rare instances in solid-state materials. The original report on PtSiSb indicates synthetic challenges for the isolation of single-phase samples, as well as potential discrepancies in the predicted electronic structure and experimental transport properties. In this work, we considered PtSiSb as a case study with the objective of demonstrating the importance of synthesis as a bridge between theoretical predictions and experimental properties. Using guidance from in situ powder X-ray diffraction, the synthesis of single-phase polycrystalline samples was developed, allowing further experimental characterization of transport and catalytic properties. Thermal conductivity and heat capacity demonstrated anomalous behavior, indicating an unconventional phonon scattering mechanism in PtSiSb. Charge transport measurements indicated semimetallic behavior. Hydrogen evolution reaction (HER) tests in acidic electrolyte revealed that the PtSiSb cathode needs an overpotential η 10 = 120 mV to achieve a current density j = −10 mA/cm 2 . This is lower activity than elemental Pt but comparable to known and active transition-metal phosphide HER catalysts, such as CoP, with η 10 = 93 mV when measured at identical conditions.

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

CoSn-type NiIn1–xSbx (0 ≤ x ≤ 0.17): Site-Selective Substitution, Electronic Structure, Chemical Bonding, and Structural Transformation

CoSn-type intermetallic compounds have emerged as a model platform for Kagome-derived flat-band physics, where subtle chemical perturbations can strongly influence electronic structure and phase stability. Here, we present a combined experimental and theoretical study of Sb-substitution in CoSn-type NiIn1–xSbx (0 ≤ x ≤ 0.17) to elucidate the interplay between site selectivity, solubility limit, chemical bonding, and electronic structure. Rietveld refinements on Neutron powder diffraction data confirmed the selective Sb-substitution at the electron-rich In2 (2d) site forming the honeycomb substructure, while the In1 (1a) site within the Kagome layer remains exclusively occupied by In. Density functional theory (DFT) calculations revealed that pristine CoSn-type NiIn hosts Ni 3d-dominated flat bands near the Fermi level (EF), originating from the Kagome-like Ni substructure. Partial replacement of In by Sb within the honeycomb layer alters these flat-band features below EF, reducing the density of states and suppressing the flat-band topology near the Fermi level. Orbital-resolved electronic structure and chemical-bonding analyses show that Sb-substitution enhances Ni-p-block (In/Sb) covalency and optimizes charge compensation, stabilizing the CoSn-type structure up to the solubility limit x ≈ 0.17. Beyond the limit, the higher-Sb compositions show satellite reflections consistent with an incommensurately modulated phase. These results establish a link between site-selective chemical substitution, bonding optimization, and flat-band electronic structure evolution, providing fundamental insights into how chemical substitution influences the electronic properties of Kagome-based intermetallic compounds.

Roy, Nilanjan [National Institute of Technology Si↗
DOE OSTI2020

A Computational Assessment of Actinide Dioxide Cations AnO 2 2+ for An = U to Lr: The Limited Stability Range of the Hexavalent Actinyl Moiety, [O=An=O] 2+

The isolated gas-phase actinide dioxide dications, AnO 2 2+ , were evaluated by DFT and coupled cluster CCSD(T) calculations for 12 actinides, An = U–Lr. CASSCF calculations were used to define the orbitals for the CCSD(T) calculations. The characteristic linear [O$=$An$=$O] 2+ hexavalent actinyl(VI) was found to be the lowest energy structure for An = U, Np, and Pu, which also form stable actinyl(VI) species in solution and possibly for Am when spin–orbit effects are included. For Am, there is a divalent [An II (O 2 )] 2+ structure where the dioxygen is an end-on physisorbed η 1 - 3 O 2 2 kcal/mol above the actinyl when spin–orbit effects are included which lower the energy of the actinyl structure. For An = Cm, Bk, and Lr, the lowest energy structure is trivalent [An III (O 2 – )] 2+ where the dioxygen is a side-on superoxide, η 2 -O 2 – . For Cm, the actinyl is close in energy to the ground state when spin–orbit effects are included. For An = Cf, Es, Fm, Md, and No, the lowest energy structure is divalent [An II (O 2 )] 2+ where the dioxygen is an end-on physisorbed η 1 - 3 O 2 . The relative energies suggest that curyl(VI) and berkelyl(VI), like well-known americyl(VI), might be stabilized by coordinating ligands in condensed phases. The results further indicate that for californyl and beyond, the actinyl(VI) moieties will probably be elusive even using strong donor ligands. The prevalence of low oxidation states (OSs) An(II) and An(III) for transplutonium actinides reflects stabilization of the 5f orbitals and validates established trends, including the remarkably high stability of divalent No. Finally, bond distances and other parameters suggest maximum bond covalency around plutonyl(VI), with a particularly substantial decrease in bond strength between americyl(VI) and curyl(VI).

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

Theoretical Prediction and Interpretation of 237 Np Mössbauer Isomer Shifts

Different theoretical approaches for the calculation of 237 Np Mössbauer isomer shifts are investigated. The traditional contact density route is compared to a previously proposed alternative approach that uses energy derivatives with respect to the nuclear radius. Both approaches yield similar results as long as suitable basis sets augmented with large exponents and relativistic methods are used. Density functional theory (DFT) calculations do not show a strong dependency of the 237 Np isomer shift on the chosen functional. Wavefunction calculations show that dynamic electron correlation can be important when covalent bonding influences the isomer shift. Effects from spin–orbit coupling are small. The isomer shifts of ionic solids and Np(III) organometallic complexes are largely governed by the oxidation state of Np. Isomer shifts of organometallic Np(IV) complexes are strongly affected by donation bonding. Furthermore, detailed analysis of the wavefunction results with different active spaces demonstrates that correlation among the outer core Np and occupied ligand frontier orbitals contributes significantly to isomer shifts of Np(IV) compounds.

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

237 Np Mössbauer Isomer Shifts: A Lesson About the Balance of Static and Dynamic Electron Correlation in Heavy Element Complexes

A large set of neptunium compounds with different oxidation states (III to VII) was assembled to study the Mössbauer isomer shift by wave function calculations and better understand covalency in f-elements complexes. The contact density approach was used to calculate the isomer shift using complete active space self-consistent field (CASSCF) multiconfiguration wave functions, as well as matrix product states [from Density Matrix Renormalization Group (DMRG) algorithms] for large active spaces. Dynamic correlation effects for the isomer shifts were treated via CASPT2 energy derivatives with respect to the nuclear radius. The CASSCF calculations appear to produce different orbital overlocalization errors for low and high Np oxidation states. For compounds with low Np oxidation numbers, the errors can be attributed to the overlocalization of the 5f orbitals. For the compounds with high Np oxidation numbers, the main errors arise from an overlocalization of ligand orbitals and concomitant to weak donation bonding. Attempts to mitigate the overlocalization errors with large active spaces using DMRG were only partially successful, showing that explicit treatment of dynamic correlation is necessary for accurate predictions of Mössbauer isomer shifts. The CASPT2 calculations perform very satisfactorily. For a subset of Np compounds, both static and dynamic correlation effects were substantial. Furthermore, a rational active space selection based on orbital entanglement diagrams proved beneficial for determining the optimal reference wave function.

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

The Good, the Bad, and the Ugly: Pseudopotential Inconsistency Errors in Molecular Applications of Density Functional Theory

The pseudopotential (PP) approximation is one of the most common techniques in computational chemistry. Despite its long history, the development of custom PPs has not tracked with the explosion of different density functional approximations (DFAs). As a result, the use of PPs with exchange/correlation models for which they were not developed is widespread, although this practice is known to be theoretically unsound. The extent of PP inconsistency errors (PPIEs) associated with this practice has not been systematically explored across the types of energy differences commonly evaluated in chemical applications. Here, we evaluate PPIEs for a number of PPs and DFAs across 196 chemically relevant systems of both transition-metal and main-group elements, as represented by the W4-11, TMC34, and S22 data sets. Near the complete basis set limit, these PPs are found to cleanly approach all-electron (AE) results for noncovalent interactions but introduce root-mean-squared errors (RMSEs) upwards of 15 kcal mol –1 into predictions of covalent bond energies for a number of popular DFAs. We achieve significant improvements through the use of empirical atom- and DFA-specific PP corrections, indicating considerable systematicity of the PPIEs. The results of this work have implications for chemical modeling in both molecular contexts and for DFA design, which we discuss.

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

The General Atomic and Molecular Electronic Structure System (GAMESS): Novel Methods on Novel Architectures

The primary focus of GAMESS over the last 5 years has been the development of new high-performance codes that are able to take effective and efficient advantage of the most advanced computer architectures, both CPU and accelerators. These efforts include employing density fitting and fragmentation methods to reduce the high scaling of well-correlated (e.g., coupled-cluster) methods as well as developing novel codes that can take optimal advantage of graphical processing units and other modern accelerators. Because accurate wave functions can be very complex, an important new functionality in GAMESS is the quasi-atomic orbital analysis, an unbiased approach to the understanding of covalent bonds embedded in the wave function. Finally, best practices for the maintenance and distribution of GAMESS are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗