Engineering topics

Schaefer, III, Henry F. [University of Georgia, Athens, GA (United States)] (ORCID:0000000302522083)

Publications and source records attributed to Schaefer, III, Henry F. [University of Georgia, Athens, GA (United States)] (ORCID:0000000302522083).

9 records

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Linear-Scaling Asymmetric Triples Correction through the Solution of the DLPNO–CCSD Lambda Equations: DLPNO–CCSD(T) Λ

In this research, we derive equations for solving for the stationary points of the DLPNO–CCSD Lagrangian, in the t 1 -transformed formalism introduced earlier and as currently implemented in the P SI 4 quantum chemistry software package. These lambda equations in the local pair natural orbital basis allow for the evaluation of CCSD(T) Λ energetics with linear-scaling computational effort, also known as the asymmetric triples correction. This DLPNO–CCSD(T) Λ method allows for accurate triples contributions to be computed for larger molecules, especially in cases that CCSD(T) is known to be insufficient, such as with multireference systems and bond-breaking systems. We showcase the accuracy of our code on reaction energies, barrier heights, and noncovalent interaction energies. Also showcased are the capabilities of our code by evaluating DLPNO–CCSD(T) Λ energetics on large noncovalent dimers up to 112 atoms, as well as a rhodium catalyst complex containing 66 atoms.

Cluster chemistry↗

DOE OSTI2026

Concordant Mode Approach (CMA): Vibrational Analysis of New and Upgraded Intermolecular Benchmarks for Noncovalent Bonding

The Concordant Mode Approach (CMA) is a novel method that offers tremendous potential for increasing the system size and the level of theory attainable in quantum chemical computations of molecular vibrational frequencies. To investigate the extension of CMA to intermolecular vibrations, computations with coupled cluster singles and doubles with perturbative triples theory [CCSD(T)] using two augmented correlation-consistent polarized-valence triple-ζ basis sets (aug-cc-pVTZ or h-aug-ccpVTZ) were performed on 17 prototypical loosely bound complexes of hydrogen-bonded, dispersion, and mixed character. These Level A results provide new and upgraded benchmarks for noncovalent bonding and a severe test for CMA vibrational analyses. The Level A target frequencies were recovered remarkably well using second-order Møller−Plesset perturbation theory (MP2) with h-aug-cc-pVTZ for generating the underlying (Level B) normal modes of the CMA scheme. Employing this Level B within the lowest-rung CMA-0A method reproduces the 435 benchmark frequencies with a mean absolute error (MAE) of 0.23 cm −1 and a corresponding standard deviation (σ) of 0.84 cm −1 ; strikingly, the corresponding subset of 106 interfragment frequencies exhibits MAE = 0.34 cm −1 and σ = 0.90 cm −1 . Subsequent application of the higher-rung CMA-2A scheme eliminates all outliers and reduces the overall MAE to a minuscule 0.08 cm −1 with the inclusion of only 3.0% of the off-diagonal couplings not accounted for by CMA-0A. Accordingly, the highly efficient CMA methodology proves to be robust even for vibrations on flat potential energy surfaces.

Aromatic compounds↗

DOE OSTI2025

Vibrational Signatures of Unrealized Phosphorus Suboxide Intermediates in White Phosphorus Oxidation Reactions

White phosphorus ignition notoriously produces the phosphoric acid anhydride P 4 O 10 , yet the intermediate oxidation steps remain undetermined. Here, we report the first geometric and vibrational characterization of two P 4 O 2 isomers, P 3 OPO and P 3 PO 2 , and substantiate a previously proposed cyclic P 4 O 2 isomer. We formally assign the infrared bands observed by Andrews and Mielke at 898 and 891 cm −1 to the antisymmetric P− O−P vibrations of P 3 OPO species. Additional bands corresponding to terminal P=O and −PO 2 stretches of P 3 OPO and P 3 PO 2 discussed herein also went unrecognized due to the peculiar bonding of P 4 O 2 species compared to oxo-bridged P 4 O x (x = 3−6) species. Sequential addition of oxygen atoms to the P 4 tetrahedron appears to form P 3 OPO and P 3 PO 2 , while cyclic P 4 O 2 is formed from P 2 O dimerization. CCSD(T) geometries, CCSD(T) + MP2[δVPT2] fundamental frequencies, and enthalpies of formation extrapolated using focal-point analysis are reported. The predicted enthalpies of formation relative to tetrahedral P 4 plus 3 Σ g − O 2 for bent-P 3 OPO, extended-P 3 OPO, P 3 PO 2 , and cyclic-P 4 O 2 are −93, −88, −85, and −63 kcal mol −1 , respectively.

Molecular structure↗

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 OSTI2025

Tests of the DFT Ladder for the Fulminic Acid Challenge

Properties of the historically pivotal fulminic acid (HCNO) molecule have been computed with a panoply of 473 density functionals of all varieties, providing a snapshot of the performance of contemporary density functional theory (DFT) for a challenging chemical system. Exhaustive tabulations and statistical analyses have been carried out for geometric parameters, vibrational frequencies, barriers to linearity, and the HCN–O dissociation energy. As the DFT ladder is climbed, confusion rather than consensus ensues regarding the details of the distinctive, extremely flat H–C–N bending potential of fulminic acid and whether the equilibrium structure is linear or bent. While high-ranking DFT functionals produce the smallest errors for the HCN + O( 3 P) → HCNO reaction energy, lower rungs emerge as the best performers for many of the bond distances and harmonic vibrational frequencies. This research shows that the current DFT zoo of approximations does not constitute a transparent ladder of increasingly accurate methods that consistently converges on definitive predictions for various properties of HCNO. Additional analyses are performed on the side effects of popular dispersion corrections on the covalently bonded properties and thermochemistry of HCNO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Linear-Scaling Local Natural Orbital-Based Full Triples Treatment in Coupled-Cluster Theory

We present an efficient, asymptotically linear-scaling implementation of the canonically O(N 8 ) coupled-cluster method with singles, doubles, and full triples excitations (CCSDT) method. We apply the domain-based local pair natural orbital (DLPNO) approach for computing CCSDT amplitudes. Our method, called DLPNO–CCSDT, uses the converged coupled-cluster amplitudes from a preceding DLPNO–CCSD(T) computation as a starting point for the solution of the CCSDT equations in the local natural orbital basis. To simplify the working equations, we t1-dress our two-electron integrals and Fock matrices, allowing our equations to take on the form of CCDT. With appropriate parameters, our method can recover more than 99.99% of the total canonical CCSDT correlation energy. In addition, we demonstrate that our method consistently yields sub-kJ mol –1 errors in relative energies when compared to canonical CCSDT, and, likewise, when computing the difference between CCSDT and CCSD(T). Finally, to highlight the low scaling of our algorithm, we present timings on linear alkanes (up to 30 carbons and 730 basis functions) and water clusters (up to 131 water molecules and 3144 basis functions).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Stabilization of [(N 5 ) 2 BX] 2– and [(N 5 ) 2 B 2 X 2 ] 2– (X = H, F, Cl, Br) by Conjugation and Hyperconjugation Effects

The isolation of nucleophilic boron bases has led to a paradigm shift in boron chemistry. Previous studies of the bis(carbene) borylene complexes revealed that these compounds possess strong donor abilities, and their reaction inertness is due to the large steric hindrance between boron reagents and reactant. In the present study, we have theoretically studied the [(N 5 ) 2 BX] 2– and [(N 5 ) 2 B 2 X 2 ] 2– compounds (X = H, F, Cl, Br). Their electronic structures and properties are discussed by using the NBO, LOL, and ELF methods. We found that both π-conjugation and hyperconjugation effects can effectively stabilize the substituted nucleophilic anionic boron compounds [(N 5 ) 2 BX] 2– and [(N 5 ) 2 B 2 X 2 ] 2– . Substituents, especially X = H, stabilize the boron center through highly delocalized π-bonding, involving the formally “empty” in-plane p orbitals of the boron atom. While the halogen substituents have high electron withdrawal ability, leading to systems being less stable, we suggest the borinium anions [(N 5 ) 2 BH] 2– and [(N 5 ) 2 B 2 H 2 ] 2– as possible synthetic targets of novel environmentally friendly catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Why Are 5-Thioglycopyranosyl Donors More Axially Selective than their Glycopyranosyl Counterparts? A Low and Variable Temperature NMR Spectroscopy and Computational Study

5-Thioglycopyranosyl donors differ in reactivity and selectivity from simple glycopyranosyl donors. An extensive study has been conducted on the nature and stability of the reactive intermediates generated on the activation of per-O-acetyl and per-O-methyl 5-thioglucopyranosyl donors and the corresponding glucopyranosyl donors. Variable temperature nuclear magnetic resonance (NMR) studies with per-O-methylated or per-O-acetyl glycosyl sulfoxides and trichloroacetimidates on activation with trifluoromethanesulfonic anhydride or trimethylsilyl triflate are reported. These show that following initial adduct formation with the promoter conversion of the 5-thioglucopyranosyl donors to the 5-thioglucopyranosyl triflates requires higher temperatures than that of the glucopyranosyl donors to the glucopyranosyl triflates. It is demonstrated that neighboring group participation is a less important phenomenon for the peracetylated thioglucosyl donors than for the peracetylated glucosyl donors. A simple thiocarbenium ion was generated by protonation of 2,3-dihydro-4H-thiopyran at low temperature and characterized by NMR spectroscopy. However, the corresponding 5-thioglucopyranosyl thenium ions were not observed in any of the NMR studies of the 5-thiopyranosyl donors: the electron-withdrawing C–O bonds around the thiopyranoside core discourage thiocarbenium ion formation, just as they discourage oxocarbenium ion formation. Density functional theory (DFT) calculations reveal the tetrahydrothiopyranyl thiocarbenium ion to be approximately 2.5 kcal/mol lower in energy than the corresponding tetrahydropyranyl oxocarbenium ion relative to the corresponding covalent triflates. However, the computations reveal a 5.8 kcal/mol activation barrier for conversion of the tetrahydrothiopyranyl triflate to the thiocarbenium ion, while formation of the oxocarbenium ion–triflate ion pair from tetrahydropyranyl triflate requires only 2.6 kcal·mol –1 . Overall, the greater axial selectivity of 5-thioglycopyranosyl donors compared to analogous glycopyranosyl donors derives from (i) the lower kinetic reactivity necessitating higher reaction temperatures, (ii) the greater stability of the thiocarbenium ion over the oxocarbenium ion facilitating equilibration under thermodynamic conditions, (iii) the greater magnitude of the anomeric effect in the 5-thiosugars, and (iv) decreased neighboring group participation in the per-esterified 5-thiosugars.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Convergent Concordant Mode Approach for Molecular Vibrations: CMA-2

The concordant mode approach (CMA) is a promising new scheme for dramatically increasing the system size and level of theory achievable in quantum chemical computations of molecular vibrational frequencies. Here, we achieve advances in the CMA hierarchy by computations targeting CCSD(T)/cc-pVTZ (coupled cluster singles and doubles with perturbative triples using a correlation-consistent polarized-valence triple-ζ basis set) benchmarks within the G2 molecular test set, executing a statistical analysis for 1501 frequencies from 111 compounds and then separately solving the refractory case of pyridine. First, MP2/cc-pVTZ (second-order Møller–Plesset perturbation theory with the same basis set) proves to be an excellent and preferred choice for generating the underlying (Level B) normal modes of the CMA scheme. Utilizing this Level B within the CMA-0A method reproduces the 1501 benchmark frequencies with a mean absolute error (MAE) of only 0.11 cm –1 and an attendant standard deviation of 0.49 cm –1 . Second, a convergent CMA-2 method is constituted that allows efficient computation of higher level (Level A) frequencies to any reasonable accuracy threshold by using only Hartree–Fock (HF) and MP2 or density functional theory (DFT) data to generate ξ parameters, which select the sparse off-diagonal force field elements for explicit evaluation at Level A. When Level B = MP2/cc-pVTZ, a cutoff of ξ = 0.02 provides an average maximum absolute error per molecule of only 0.17 cm –1 by incurring merely a 33% increase in average cost over CMA-0A. This CMA-2 method also eradicates the 4 problematic CMA-0A outliers of pyridine with even less effort (ξ = 0.04, 22% increase). Finally, the newly developed CMA procedures are shown to be highly successful when applied to 1-(1H-pyrrol-3-yl)ethanol, a new test molecule with diverse types of vibration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗