Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dimer”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Assembly of Uranyl Peroxides from Ball Milled Solids

Mechanochemistry enables transformations of highly insoluble materials such as uranium dioxide or the mineral studtite [(UO 2 )(O 2 )(H 2 O) 2 ]·(H 2 O) 2 into uranyl triperoxide compounds that can subsequently assemble into hydroxide-bridged uranyl peroxide dimers in the presence of lithium hydroxide. Dissolution of these solids in water yields uranyl peroxide nanoclusters including U 24 , Li 24 [(UO 2 )(O 2 )(OH)] 24 . Insoluble uranium solids can transform into highly soluble uranyl peroxide phases in the solid state with miniscule quantities of water. Furthermore, such reactions are potentially applicable to uranium processing in the front and back end of the nuclear fuel cycle.

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

Tethering Effects in Oligomer-Based Metal–Organic Frameworks

Metal–organic frameworks (MOFs) can be con- structed using conventional molecular linkers or polymeric linkers (polyMOFs), but the relationship and relative properties of these related materials remain understudied. As an intermediate between these two extremes, a library of oligomeric ligand precursors (dimers, trimers) was used to prepare a series of oligomeric-linker MOFs (oligoMOFs) based on the prototypical IRMOF-1 system. IRMOF-1 was found to be remarkably tolerant to a wide variety of oligomeric linkers, the use of which greatly enhanced the MOF yield and prevented framework interpenetration. Tether length-dependent ordering of ligand and metal cluster orientations was also observed in these oligoMOFs. Improved low-humidity stability was found in oligoIRMOF-1 samples, with surface area preservation varying as a function of tether length and a complete suppression of crystalline hydrolysis products for all oligoIRMOF-1 materials. Furthermore, these findings pave the way toward a better understanding of the structure–function relationships between monomeric, oligomeric, and polymeric MOFs and highlight an underutilized strategy for tuning MOF properties.

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

Two Distinct Cu(II)-V(IV) Superexchange Interactions with Similar Bond Angles in a Triangular ?CuV2? Fragment

The strength and sign of superexchange interactions are often predicted on the basis of the bond angles between magnetic ions, but complications may arise in situations with a nontrivial arrangement of the magnetic orbitals. We report on a novel molecular tetramer compound [Cu(H2O)dmbpy]2[V2O2F8] (dmbpy = 4,4 '-dimethyl-2,2 ' bipyridyl) that is composed of triangular "CuV2" fragments and displays a spin gap behavior. By combining first-principles calculations and electronic models, we reveal that superexchange Cu-V interactions carry drastically diHerent coupling strengths along two Cu-F-V pathways with comparable bond angles in the triangular "CuV2" fragment. Counterintuitively, their strong disparity is found to originate from the restricted symmetry of the half-filled Cu dx2-y2 orbital stabilized by the crystal field, leading to one dominating antiferromagnetic Cu-V coupling in each fragment. We revisit the magnetic properties of the reported spin-gapped chain compound [enH2]Cu(H2O)2[V2O2F8] (enH2 = ethylene diammonium) containing similar triangular "CuV2" fragments, and the magnetic behavior of the molecular tetramer and the chain compounds is rationalized as that of weakly coupled spin dimers and spin trimers, respectively. This work demonstrates that fundamentally diHerent magnetic couplings can be observed between magnetic ions with similar bond angles in a single spin motif, thus providing a strategy to introduce various exchange interactions combined with low dimensionality in heterometallic Cu(II)-V(IV) compounds.

Complexes↗
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

Elucidating Actinide–Pertechnetate and Actinide–Perrhenate Bonding via a Family of Th–TcO 4 and Th–ReO 4 Frameworks and Solutions

Technetium-99, a β-emitter produced from 235 U fission, poses a challenge for the nuclear industry due to co-extraction of pertechnetate (TcO 4 – ) with the actinides (An) during nuclear fuel reprocessing. Previous studies suggested that direct coordination of pertechnetate with An plays an important role in the coextraction process. However, few studies have provided direct evidence for An–TcO 4 – bonding in the solid state, and even fewer in solution. The present study describes synthesis and structural elucidation of a family of thorium(IV)-pertechnetate/perrhenate (ReO 4 – , nonradioactive surrogate) compounds, which is obtained by dissolution of thorium oxyhydroxide in perrhenic/pertechnic acid followed by crystallization, with or without heating. For reaction ratios of 3:1, 4:1, and 6:1 MO 4 – /Th(IV) (M = Tc, Re), the crystallized compounds reflect the same ratio, suggesting facile and flexible coordination. Furthermore, nine structures reveal 1-dimensional and 2-dimensional frameworks with varying topologies. While a multitude of compounds isolated from 4:1 (and 6:1) reaction solutions feature Th monomers linked by MO 4 – , the 3:1 reaction solution yielded the well-known dihydroxide-bridged thorium dimer, linked, and capped by MO 4 – . Density functional theory calculations on ReO 4 – /TcO 4 – isomorphs suggest similar bonding characteristics in the solid state, but experimental solution characterization noted differences. Specifically, small-angle X-ray scattering studies suggest the bonding of Th–TcO 4 – persists in solution, while Th–ReO 4 – bonding is less apparent.

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

Structural Investigation of Technetium Dibutylphosphate Species Using X-ray Absorption Fine Structure Spectroscopy

The speciation of Tc after the extraction of Tc(IV) from H 2 O and 1 M HNO 3 by dibutylphosphoric acid (HDBP) in dodecane has been studied by X-ray absorption fine structure (XAFS) spectroscopy. Results show the formation of dimeric species with Tc 2 O 2 and Tc 2 O units, and the formulas [Tc 2 O 2 (DBP·HDBP) 4 ] (1) and [Tc 2 O(NO 3 ) 2 (DBP) 2 (DBP·HDBP) 2 ] (2) were, respectively, proposed for the species extracted from H 2 O and 1 M HNO 3 . The interatomic Tc–Tc distances found in the Tc 2 O 2 and Tc 2 O units [2.55(3) and 3.57(4) Å, respectively] are similar to the ones found in Tc(IV) dinuclear species. It is likely that the speciation of Tc(IV) in dodecane is due to the extraction of a species with a Tc 2 O unit for (2) and to the redissolution of a Tc(IV)-DBP solid for (1). The XAFS results for (1) and (2) were compared to that obtained for the extraction of Tc(IV) with TBP/HDBP/dodecane from 0.5 M HNO 3 , (3) which highlight the formation of Tc mononuclear nitrate species {i.e. [Tc(NO 3 ) 3 (DBP)] or [Tc(NO 3 ) 2 (DBP·HDBP)]}. These results confirm the importance of the preparation and speciation of the Tc(IV) aqueous solutions prior to extraction and how much this influences and drives the final Tc speciation in organic extraction. Here, these studies outline the complexity of Tc separation chemistry and provide insights into the behavior of Tc during the reprocessing of used nuclear fuel.

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

Inducing Ferrimagnetic Exchange in 1D-FeSe 2 Chains Using Heteroleptic Amine Complexes: [Fe(en)(tren)][FeSe 2 ] 2

[Fe(en)(tren)][FeSe 2 ] 2 (en = ethylenediamine, C 2 H 8 N 2 , tren = tris(2-aminoethyl)amine, C 6 H 18 N 4 ) has been synthesized by a mixed-ligand solvothermal method. Its crystal structure contains heteroleptic [Fe(en)(tren)] 2+ complexes with distorted octahedral coordination, incorporated between 1D-FeSe 2 chains composed of edge–sharing FeSe 4 tetrahedra. The twisted octahedral coordination environment of the Fe–amine complex leads to partial dimerization of Fe–Fe distances in the FeSe 2 chains so that the FeSe 4 polyhedra deviate strongly from the regular tetrahedral geometry. 57 Fe Mössbauer spectroscopy reveals oxidation states of +3 for the Fe chain atoms and +2 for the Fe complex atoms. The close proximity of Fe atoms in the chains promotes ferromagnetic nearest neighbor interactions, as indicated by a positive Weiss constant, θ = +53.8(6) K, derived from the Curie–Weiss fitting. Magnetometry and heat capacity reveal two consecutive magnetic transitions below 10 K. DFT calculations suggest that the ordering observed at 4 K is due to antiferromagnetic intrachain interactions in the 1D-FeSe 2 chains. The combination of two different ligands creates an asymmetric coordination environment that induces changes in the structure of the Fe–Se fragments. Finally, this synthetic strategy opens new ways to explore the effects of ligand field strength on the structure of both Fe-amine complexes and surrounding Fe–Se chains.

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

The Role of the Organic Solvent Polarity in Isolating Uranyl Peroxide Capsule Fragments

Uranyl peroxide capsules are a fascinating class of polyoxometalates (POMs), discovered only in the 21st century. Understanding the reactivity between peroxide, alkali cations, and uranyl in alkaline solutions is important in nuclear science disciplines including mineralogy, nuclear energy, and legacy nuclear wastes. In this work, we have developed a general procedure to isolate different fragments of the uranyl-peroxide POM capsules, using organic solvents to partially remove K + salts from crude solids of the monomer building block UO 2 (O 2 ) 3 4– (K–U 1 ), leading to stabilization of these reactive fragments. Higher polarity organic solvents remove more K + salts from the crude solid, owed to higher solubility, resulting in more extensive linking of uranyl peroxide building units. By this strategy we have isolated and structurally characterized a dimer K 6 [(UO 2 ) 2 (O 2 ) 4 (OH) 2 ]·7H 2 O (K–U 2 ) and a hexamer face frequently observed in the capsules, K 12 [(UO 2 ) 6 (O 2 ) 9 (OH) 6 ]$·x$H 2 O(K–U 6 ). Comparing experimental and computed Raman spectra shows that these intermediates crystallize by a solid-to-solid transformation, via polymerization of the monomer building block. By small-angle X-ray scattering (SAXS), we track the conversion of the fragments to POM capsules; the reaction rate increases from K–U 1 (days) < K–U 2 (hours) < K–U 6 (instantaneous). This study provides a general synthetic procedure to isolate metastable uranyl peroxide oligomers and control the oligomerization, which will be later applied to systems with the heavier alkalis that are even less stable.

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

High-Pressure, High-Temperature Synthesis and Characterization of Polar and Magnetic LuCrWO 6

We report a new polar and magnetic oxide, LuCrWO 6 , was synthesized under high pressure (6 GPa) and high temperature (1673 K). LuCrWO 6 is isostructural with the previously reported polar YCrWO 6 (SG: Pna21, no. 33). The ordering of CrO 6 and WO 6 octahedra in the edge-shared dimers induce the polar structure. The effective size of rare earth, Ln cation does not seem to affect the symmetry of LnCrWO 6 . Second harmonic generation measurements of LuCrWO 6 confirmed the noncentrosymmetric character and strong piezoelectric domains are observed from piezoresponse force microscopy at room temperature. LuCrWO 6 exhibits antiferromagnetic behavior, T N , of ~18 K with a Weiss temperature of -30.7 K.

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

New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions

Insensitive high explosives based on TATB (1,3,5-triamino-2,4,6-trinitrobenzene) are needed in applications when safety is of paramount importance, but the basic material properties that give rise to its insensitivity are not fully understood. Molecular dynamics modeling using empirical force fields (FFs) has been the main route to characterize many complicated dynamical properties of TATB single crystal, but these FFs have not been comprehensively tested at extreme conditions typical of detonation. We collect a benchmark data set of (quasi)static TATB physical properties as determined by experiments and electronic structure calculations and apply this data set to validate four existing TATB FFs along with a new TATB FF that we develop here and denote as the CEA-LLNL-Missouri (CLM) FF. Benchmark data include vibrational spectra, the TATB crystal temperature–pressure–volume equation of state and lattice parameters, properties of TATB crystal polymorphs and transitions to the gaseous and liquid states, dimer energy landscapes, the pressure-dependent elastic tensor, and the energy landscape for inelastic deformation via sliding of TATB crystal layers. As a general assessment, we find that the two existing nonreactive FFs are more accurate in describing TATB’s physical properties compared to the two variants of the ReaxFF reactive FF considered. The new CLM FF is found to consistently yield similar or better agreement with experiments and electronic structure theory than any of the existing FF models, and it presents a distinct improvement in accurately modeling TATB elasticity and equation of state. So this work is expected to help improve the accuracy of FF-based modeling of complicated dynamic responses that ultimately govern the safety and performance characteristics of this material.

36 MATERIALS SCIENCE↗
DOE OSTI2022

All-Atom Simulations Uncover Structural and Dynamical Properties of STING Proteins in the Membrane System

Recent studies have shown that the stimulator of interferon gene (STING) protein plays a central role in the immune system by facilitating the production of type I interferons in cells. The STING signaling pathway is also a prominent activator of cancer-killing T cells that initiate a powerful adaptive immune response. Since biomolecular signaling pathways are complicated and not easily identified through traditional experiments, molecular dynamics (MD) has often been used to study structural and dynamical responses of biological pathways. In this work, we carried out MD simulations for full-length chicken and human STING (chSTING and hSTING) proteins. Specifically, we investigated ligand-bound closed (holo) and ligand-unbound open (apo) forms of STING in the membrane system by comparing their conformational and dynamical differences. Our research provides clues for understanding the mechanism of the STING signaling pathway by uncovering detailed insights for the examined systems: the residues from each chain in the binding pocket are strongly correlated to one another in the open STING structure compared with those in the closed STING structure. Ligand-bound closed STING displays ~174° rotation of the ligand-binding domain (LBD) relative to the open STING structure. The dynamical analysis of residue Cys148 located in the linker region of hSTING does not support the earlier hypothesis that Cys148 can form disulfide bonds between adjacent STING dimers. We also demonstrate that using the full-length proteins is critical, since the MD simulations of the LBD portion alone cannot properly describe the global conformational properties of STING.

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

Free Energy and Flexibility Analysis of Autoinhibited Human BRAF

The RAF serine/threonine protein kinases function as direct effectors of RAS in the intracellular transmission of extracellular growth signals, and they are key targets for drug discovery, given the high incidence of oncogenic mutations in RAF and other components of this signaling pathway. In its inactive state, RAF is held in an autoinhibited conformation in the cytosol through a combination of intramolecular interactions and binding to a regulatory 14−3−3 protein dimer. Activation of RAF is initiated by its interaction with membrane-localized GTP-bound RAS, which induces conformational changes that release RAF from its autoinhibited state. However, the molecular mechanisms governing RAF activation remain incomplete, largely due to the challenges in experimentally capturing the intermediate conformational states in this process. To address this gap, we developed a comprehensive all-atom model of BRAF based on existing cryo-EM structures. Using this model, we performed extensive molecular dynamics simulations to evaluate the stability and free energy landscape of autoinhibited BRAF in solution. Our analysis reveals conformational flexibility within the autoinhibited complex, suggesting that this dynamic behavior may play a role in facilitating BRAF activation upon engagement with the membrane-bound RAS.

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

Chemistry on Quantum Computers with Virtual Quantum Subspace Expansion

Several novel methods for performing calculations relevant to quantum chemistry on quantum computers have been proposed but not yet explored experimentally. Virtual quantum subspace expansion is one such algorithm developed for modeling complex molecules using their full orbital space and without the need for additional quantum resources. Here, we implement this method on the IBM Q platform and calculate the potential energy curves of the hydrogen and lithium dimers using only two qubits and simple classical post-processing. A comparable level of accuracy would require twenty qubits with previous approaches. We also develop an approach to minimize the impact of experimental noise on the stability of a generalized eigenvalue problem that is a crucial component of the algorithm. Our results demonstrate that virtual quantum subspace expansion works well in practice.

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

Confirmation of the PPLB Derivative Discontinuity: Exact Chemical Potential at Finite Temperatures of a Model System

The landmark 1982 paper of Perdew, Parr, Levy, and Balduz (often called PPLB) laid the foundation for our modern understanding of the role of the derivative discontinuity in density functional theory, which drives much development to account for its effects. A simple model for the chemical potential at vanishing temperature played a crucial role in their argument. We investigate the validity of this model in the simplest non-trivial system to which it can be applied and which can be easily solved exactly, the Hubbard dimer. We find exact agreement in the crucial zero-temperature limit, and show the model remains accurate for a significant range of temperatures. We identify how this range depends on the strength of correlations. We extend the model to approximate free energies accounting for the derivative discontinuity, a feature missing in standard semilocal approximations. Here, we provide a correction to this approximation to yield even more accurate free energies. We discuss the relevance of these results for warm dense matter.

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

Determining Partial Atomic Charges for Liquid Water: Assessing Electronic Structure and Charge Models

Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, natural population analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree–Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 Å of the central water) in the quantum chemical calculation converges the charges of the central water molecule. Furthermore, our final pragmatic quantum chemical charge-assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 Å.

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

Efficient Four-Component Dirac–Coulomb–Gaunt Hartree–Fock in the Pauli Spinor Representation

Four-component Dirac Hartree--Fock is an accurate mean-field method for treating molecular systems where relativistic effects are important. However, the computational cost and complexity of the two-electron interaction makes this method less common, even though we can consider the Dirac Hartree--Fock Hamiltonian the ``ground truth" of electronic structure, barring explicit quantum-electrodynamical effects. Being able to calculate these effects is then vital to the design of lower scaling methods for accurate predictions in computational spectroscopy and properties of heavy element complexes that must include relativistic effects for even qualitative accuracy. In this work, we present a Pauli quaternion formalism of maximal component- and spin-separation for computing the Dirac-Coulomb-Gaunt Hartree-Fock ground state, with a minimal floating-point-operation count algorithm. Furthermore, this approach also allows one to explicitly separate different spin physics from the two-body interactions, such as spin-free, spin-orbit, and the spin-spin contributions. Additionally, we use this formalism to examine relativistic trends in the periodic table, and analyze the basis set dependence of atomic gold and gold dimer systems.

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

A Massively Parallel Implementation of the CCSD(T) Method Using the Resolution-of-the-Identity Approximation and a Hybrid Distributed/Shared Memory Parallelization Model

In this work, a parallel algorithm is described for the coupled-cluster singles and doubles method augmented with a perturbative correction for triple excitations [CCSD(T)] using the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs). The algorithm bypasses the storage of four-center ERIs by adopting an integral-direct strategy. The CCSD amplitude equations are given in a compact quasi-linear form by factorizing them in terms of amplitude-dressed three-center intermediates. A hybrid MPI/OpenMP parallelization scheme is employed, which uses the OpenMP-based shared memory model for intranode parallelization and the MPI-based distributed memory model for internode parallelization. Parallel efficiency has been optimized for all terms in the CCSD amplitude equations. Two different algorithms have been implemented for the rate-limiting terms in the CCSD amplitude equations that entail and -scaling computational costs, where N O and N V denote the number of correlated occupied and virtual orbitals, respectively. One of the algorithms assembles the four-center ERIs requiring N V 4 and N O 2 N V 2 -scaling memory costs in a distributed manner on a number of MPI ranks, while the other algorithm completely bypasses the assembling of quartic memory-scaling ERIs and thus largely reduces the memory demand. It is demonstrated that the former memory-expensive algorithm is faster on a few hundred cores, while the latter memory-economic algorithm shows a better strong scaling in the limit of a few thousand cores. The program is shown to exhibit a near-linear scaling, in particular for the compute-intensive triples correction step, on up to 8000 cores. The performance of the program is demonstrated via calculations involving molecules with 24–51 atoms and up to 1624 atomic basis functions. As the first application, the complete basis set (CBS) limit for the interaction energy of the π-stacked uracil dimer from the S66 data set has been investigated. This work reports the first calculation of the interaction energy at the CCSD(T)/aug-cc-pVQZ level without local orbital approximation. The CBS limit for the CCSD correlation contribution to the interaction energy was found to be -8.01 kcal/mol, which agrees very well with the value -7.99 kcal/mol reported by Schmitz, Hättig, and Tew [ Phys. Chem. Chem. Phys. 2014 , 16 , 22167-22178]. The CBS limit for the total interaction energy was estimated to be -9.64 kcal/mol.

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

Comprehensive Basis-Set Testing of Extended Symmetry-Adapted Perturbation Theory and Assessment of Mixed-Basis Combinations to Reduce Cost

Hybrid or “extended” symmetry-adapted perturbation theory (XSAPT) replaces traditional SAPT’s treatment of dispersion with better performing alternatives while at the same time extending two-body (dimer) SAPT to a many-body treatment of polarization using a self-consistent charge embedding procedure. The present work presents a systematic study of how XSAPT interaction energies and energy components converge with respect to the choice of Gaussian basis set. Errors can be reduced in a systematic way using correlation-consistent basis sets, with aug-cc-pVTZ results converged within <0.1 kcal/mol. Similar (if slightly less systematic) behavior is obtained using Karlsruhe basis sets at much lower cost, and we introduce new versions with limited augmentation that are even more efficient. Pople-style basis sets, which are more efficient still, often afford good results if a large number of polarization functions are included. The dispersion models used in XSAPT afford much faster basis-set convergence as compared to the perturbative description of dispersion in conventional SAPT, meaning that “compromise” basis sets (such as jun-cc-pVDZ) are no longer required and benchmark-quality results can be obtained using triple-ζ basis sets. The use of diffuse functions proves to be essential, especially for the description of hydrogen bonds. As a result, the “δ(Hartree–Fock)” correction for high-order induction can be performed in double-ζ basis sets without significant loss of accuracy, leading to a mixed-basis approach that offers 4× speedup over the existing (cubic scaling) XSAPT approach.

74 ATOMIC AND MOLECULAR PHYSICS↗