Does Basis Set Superposition Error Significantly Affect Post-CCSD(T) Corrections?

We have investigated the title question for both a subset of the W4-11 total atomization energies benchmark, and for the A24x8 noncovalent interactions benchmark. Overall, counterpoise corrections to post-CCSD(T) contributions are about two orders of magnitude less important than those to the CCSD(T) interaction energy. Counterpoise corrections for connected quadruple substitutions (Q) are negligible, and or especially so.

Basis Set Extrapolation from the Vanishing Counterpoise Correction Condition.

Basis set extrapolations are typically rationalized either from analytical arguments involving the partial-wave or principal expansions of the correlation energy in helium-like systems or from fitting extrapolation parameters to reference energetics for a small(ish) training set. Seeking to avoid both, we explore a third alternative: extracting extrapolation parameters from the requirement that the BSSE (basis set superposition error) should vanish at the complete basis set limit. We find this to be a viable approach provided that the underlying basis sets are not too small and reasonably well balanced.

Exploring the influence of ( - 1)d subvalence correlation and of spin-orbit coupling on chalcogen bonding.

This article presents a comprehensive computational investigation into chalcogen bonding interactions, focusing specifically on elucidating the role of subvalence ( - 1)d and ( - 1)sp correlation. The incorporation of inner-shell ( - 1)d correlation leads to a decrease in interaction energies for chalcogen-bonded systems (at least those studied herein), contradicting the observations regarding halogen bonding documented by Kesharwani in , , 122 (8), 2184-2197. The significance of ( - 1)sp subvalence correlation appears to be lower by an order of magnitude.

W4Λ: Leveraging Λ Coupled-Cluster for Accurate Computational Thermochemistry Approaches.

High-accuracy composite wave function methods like Weizmann-4 (W4) theory, high-accuracy extrapolated ab initio thermochemistry (HEAT), and the Feller-Peterson-Dixon (FPD) approach enable sub-kJ/mol accuracy in gas-phase thermochemical properties. Their biggest computational bottleneck is the evaluation of the valence post-CCSD(T) correction term. We demonstrate here, for the W4-17 thermochemistry benchmark and subsets thereof, that the Λ coupled-cluster expansion converges more rapidly and smoothly than the regular coupled-cluster series.

On the sensitivity of computed partial charges toward basis set and (exchange-)correlation treatment.

Partial charges are a central concept in general chemistry and chemical biology, yet dozens of different computational definitions exist. In prior work [Cho et al., ChemPhysChem 21, 688-696 (2020)], we showed that these can be reduced to at most three 'principal components of ionicity'.

Correlation Consistent Basis Sets for Explicitly Correlated Theory: The Transition Metals.

We present correlation consistent basis sets for explicitly correlated (F12) calculations, denoted VZ(-PP)-F12-wis ( = D,T), for the -block elements. The cc-pVDZ-F12-wis basis set is contracted to [8s7p5d2f] for the 3-block, while its ECP counterpart for the 4 and 5-blocks, cc-pVDZ-PP-F12-wis, is contracted to [6s6p5d2f]. The corresponding contracted sizes for cc-pVTZ(-PP)-F12-wis are [98632] for the 3-block elements and [77632] for the 4 and 5-block elements.

A double-hybrid density functional based on good local physics with outstanding performance on the GMTKN55 database.

In two recent papers [A. D. Becke, J.

The Importance of Tight f Basis Functions for Heavy p-Block Oxides and Halides: A Parallel With Tight d functions in the Second Row.

It is well-known that both wave function ab initio and DFT calculations on second-row compounds exhibit anomalously slow basis set convergence unless the basis sets are augmented with additional "tight" (high-exponent) d functions, as in the cc-pV(n+d)Z and aug-cc-pV(n+d)Z basis sets. This has been rationalized as being necessary for a better description of the low-lying 3d orbital, which as the oxidation state increases sinks low enough to act as a back-donation acceptor from chalcogen and halogen lone pairs. This prompts the question whether a similar phenomenon exists for the isovalent compounds of the heavy p-block.

Performance of Localized-Orbital Coupled-Cluster Approaches for the Conformational Energies of Longer -Alkane Chains.

We report an update and enhancement of the ACONFL (conformer energies of large alkanes [, 3521-3535]) dataset. For the ACONF12 (-dodecane) subset, we report basis set limit canonical coupled-cluster with singles, doubles, and perturbative triples [i.e.

S66x8 noncovalent interactions revisited: new benchmark and performance of composite localized coupled-cluster methods.

The S66x8 noncovalent interactions benchmark has been re-evaluated at the "sterling silver" level, using explicitly correlated MP2-F12 near the complete basis set limit, CCSD(F12*)/aug-cc-pVTZ-F12, and a (T) correction from conventional CCSD(T)/sano-V{D,T}Z+ calculations. The revised reference values differ by 0.1 kcal mol RMS from the original Hobza benchmark and its revision by Brauer , but by only 0.

Reduced-Scaling Double Hybrid Density Functional Theory with Rapid Basis Set Convergence through Localized Pair Natural Orbital F12.

Following earlier work [Mehta, N.; Martin, J. M.

Explicitly Correlated Double-Hybrid DFT: A Comprehensive Analysis of the Basis Set Convergence on the GMTKN55 Database.

Double-hybrid density functional theory (DHDFT) offers a pathway to accuracy approaching composite wavefunction approaches such as G4 theory. However, the Görling-Levy second-order perturbation theory (GLPT2) term causes them to partially inherit the slow ∝ (with the maximum angular momentum) basis set convergence of correlated wavefunction methods. This could potentially be remedied by introducing F12 explicit correlation: we investigate the basis set convergence of both DHDFT and DHDFT-F12 (where GLPT2 is replaced by GLPT2-F12) for the large and chemically diverse general main-group thermochemistry, kinetics, and noncovalent interactions (GMTKN55) benchmark suite.

Benefits of Range-Separated Hybrid and Double-Hybrid Functionals for a Large and Diverse Data Set of Reaction Energies and Barrier Heights.

To better understand the thermochemical kinetics and mechanism of a specific chemical reaction, an accurate estimation of barrier heights (forward and reverse) and reaction energies is vital. Because of the large size of reactants and transition state structures involved in real-life mechanistic studies (e.g.

Thermally-Activated Tunneling in the Two-Water Bridge Catalyzed Tautomerization of Phosphinylidene Compounds.

Phosphinylidenes are an important class of organophosphorus compounds that can exhibit tautomerization between tricoordinated P(III) hydroxide (R R POH) and a pentacoordinated P(V) oxide (R R P(O)H) form. Herein we show, using the canonical variational transition state theory combined with multidimensional small-curvature tunneling approximation, the dominance of proton tunneling in the two-water-bridged tautomerizations of phosphinous acid and model phosphinylidenes comprising phosphosphinates, H-phosphonates, H-phosphinates and secondary phosphine oxides. Based on the studied system, the contribution of thermally-activated tunneling is predicted to speed up the semiclassical reaction rate by ca.

Automatic generation of complementary auxiliary basis sets for explicitly correlated methods.

Explicitly correlated calculations, aside from the orbital basis set, typically require three auxiliary basis sets: Coulomb-exchange fitting (JK), resolution of the identity MP2 (RI-MP2), and complementary auxiliary basis set (CABS). If unavailable for the orbital basis set and chemical elements of interest, the first two can be auto-generated on the fly using existing algorithms, but not the third. In this paper, we present a quite simple algorithm named autoCABS; a Python implementation under a free software license is offered at Github.

MP2-F12 Basis Set Convergence near the Complete Basis Set Limit: Are Functions Sufficient?

We have investigated the title question for the W4-08 thermochemical benchmark using -saturated truncations of a large reference (REF) basis set, as well as for standard F12-optimized basis sets. With the REF basis set, the root-mean-square (RMS) contribution of functions to the MP2-F12 total atomization energies (TAEs) is about 0.01 kcal/mol, the largest individual contributions being 0.

Do Double-Hybrid Functionals Benefit from Regularization in the PT2 Term? Observations from an Extensive Benchmark.

We put to the test a recent suggestion [Shee, J., et al. , (50), 12084-12097] that MP2 regularization might improve the performance of double-hybrid density functionals.

Heavy-Atom Tunneling in the Covalent/Dative Bond Complexation of Cyclo[18]carbon-Piperidine.

Recent quantum chemical computations demonstrated the electron-acceptance behavior of this highly reactive cyclo[18]carbon (C) ring with piperidine (pip). The C-pip complexation exhibited a double-well potential along the N-C reaction coordinate, forming a van der Waals (vdW) adduct and a more stable, strong covalent/dative bond (DB) complex by overcoming a low activation barrier. By means of direct dynamical computations using canonical variational transition state theory (CVT), including the small-curvature tunneling (SCT), we show the conspicuous role of heavy atom quantum mechanical tunneling (QMT) in the transformation of vdW to DB complex in the solvent phase near absolute zero.

The MOBH35 Metal-Organic Barrier Heights Reconsidered: Performance of Local-Orbital Coupled Cluster Approaches in Different Static Correlation Regimes.

We have revisited the MOBH35 (Metal-Organic Barrier Heights, 35 reactions) benchmark [Iron; , Janes, , 2019, 123 (17), 3761-3781; ibid. 2019, 123, 6379-6380] for realistic organometallic catalytic reactions, using both canonical CCSD(T) and localized orbital approximations to it. For low levels of static correlation, all of DLPNO-CCSD(T), PNO-LCCSD(T), and LNO-CCSD(T) perform well; for moderately strong levels of static correlation, DLPNO-CCSD(T) and (T) may break down catastrophically, and PNO-LCCSD(T) is vulnerable as well.

Pure and Hybrid SCAN, rSCAN, and rSCAN: Which One Is Preferred in KS- and HF-DFT Calculations, and How Does D4 Dispersion Correction Affect This Ranking?

Using the large and chemically diverse GMTKN55 dataset, we have tested the performance of pure and hybrid KS-DFT and HF-DFT functionals constructed from three variants of the SCAN meta-GGA exchange-correlation functional: original SCAN, rSCAN, and rSCAN. Without any dispersion correction involved, HF-SCANn outperforms the two other HF-DFT functionals. In contrast, among the self-consistent variants, SCANn and rSCANn offer essentially the same performance at lower percentages of HF-exchange, while at higher percentages, SCANn marginally outperforms rSCANn and rSCANn.

Coupled Cluster Benchmark of New DFT and Local Correlation Methods: Mechanisms of Hydroarylation and Oxidative Coupling Catalyzed by Ru(II, III) Chloride Carbonyls.

We have evaluated a set of accurate canonical CCSD(T) energies for stationary points on the potential energy surface for Ru(II, III) chloride carbonyl catalysis of two competing reactions between benzene and methyl acrylate (MA), namely, hydroarylation and oxidative coupling. We have then applied this set to evaluate the performance of localized orbital coupled-cluster methods and several new and common density functionals. We find that (a) DLPNO-CCSD(T) with TightPNO cutoffs is an acceptable substitute for full canonical CCSD(T) calculations on this system; (b) for the closed-shell systems where it could be applied, LNO-CCSD(T) with tight convergence criteria is very close to the canonical results; (c) the recent ωB97X-V and ωB97M-V functionals exhibit superior performance to commonly used DFT functionals in both closed- and open-shell calculations; (d) the revDSD-PBEP86 revision of the DSD-PBEP86 double hybrid represents an improvement over the original, even though transition metals were not involved in its parametrization; and (e) DSD-SCAN and DOD-SCAN show comparable efficiency.

Surprisingly Good Performance of XYG3 Family Functionals Using a Scaled KS-MP3 Correlation.

By adding a GLPT3 (third-order Görling-Levy perturbation theory, or KS-MP3) term to the XYG7 form for a double hybrid, we are able to bring down WTMAD2 (weighted total mean absolute deviation) for the very large and chemically diverse GMTKN55 benchmark to an unprecedented 1.17 kcal/mol, competitive with much costlier composite wave function ab initio approaches. Intriguingly, (a) the introduction of makes an empirical dispersion correction redundant; (b) generalized gradient approximation (GGA) or meta-GGA semilocal correlation functionals offer no advantage over the local density approximation (LDA) in this framework; (c) if a dispersion correction is retained, then simple Slater exchange leads to no significant loss in accuracy.

Exploring Avenues beyond Revised DSD Functionals: II. Random-Phase Approximation and Scaled MP3 Corrections.

For revDSD double hybrids, the Görling-Levy second-order perturbation theory component is an Achilles' heel when applied to systems with significant near-degeneracy ("static") correlation. We have explored its replacement by the direct random phase approximation (dRPA), inspired by the SCS-dRPA75 functional of Kállay and co-workers. The addition to the final energy of both a D4 empirical dispersion correction and of a semilocal correlation component lead to significant improvements, with DSD-PBEdRPA-D4 approaching the performance of revDSD-PBEP86-D4 and the Berkeley ωB97M(2).

Exploring Avenues beyond Revised DSD Functionals: I. Range Separation, with DSD as a Special Case.

We have explored the use of range separation as a possible avenue for further improvement on our revDSD minimally empirical double hybrid functionals. Such ωDSD functionals encompass the XYG3 type of double hybrid (, DSD) as a special case for ω → 0. As in our previous studies, the large and chemically diverse GMTKN55 benchmark suite was used for evaluation.