On the performance of various hierarchized bases in extrapolating the correlation energy to the complete basis set limit

High Accuracy Ab Initio Potential Energy Curves and Dipole Moment Functions for the X1Σ+ and a3Π Spin States of the CF+ Diatomic Molecule

Potential energy curves and dipole moment functions constructed using high-accuracy ab initio methods allow for an in-depth examination of the electronic structure of diatomic molecules. Ab initio computations serve as a valuable complement to experimental data, offering insights into the nature of short-lived molecules such as those encountered within the interstellar medium (ISM). While laboratory experiments provide critical groundwork, the ISM's conditions often permit longer lifetimes for lower stability molecules, enabling unique observations. The CF+ diatomic molecule is one such molecule that has been observed spectroscopically in the ISM. Previous experimental and theoretical work have examined different spectroscopic aspects of the CF+ molecule, but the development of newer, more complete potential energy curves and dipole moment functions allows for even greater insight. We constructed both potential energy curves and dipole moment functions for the ground X1Σ+ and first excited a3Π states of CF+ for both the 12C and 13C isotopologues. The potential energy curves were constructed using coupled cluster with single, double, and perturbative triple excitations (CCSD(T)) at the complete basis set limit with corrections from full triple, quadruple, quintuple, and hextuple excitations within a finite-basis coupled cluster wave function as well as corrections from full configuration interaction and relativistic effects. Rovibrational wave functions were calculated using a vibrational Hamiltonian matrix, which moves beyond the harmonic oscillator approximation. The equilibrium bond length, vibrational constant, and rotational constant were reproduced to within 0.00013 Å, 0.28 cm-1, and 0.00045 cm-1, respectively, of experimental values. Experimental transition energies from rovibrational spectra were reproduced with an error of no larger than 0.63 cm-1. The triplet excited state (a3Π) was found to have a longer equilibrium bond length at 1.21069 Å, a vibrational constant of 1611.29 cm-1, and a rotational constant of 1.56376 cm-1. Rovibrational line lists for the 12C and 13C isotopologues for both the X1Σ+ and the excited a3Π states were generated.

Reconciling spectroscopy with dynamics in global potential energy surfaces: the case of the astrophysically relevant SiC_{2}

SiC_{2} is a fascinating molecule due to its unusual bonding and astrophysical importance. In this work, we report the first global potential energy surface (PES) for ground-state SiC_{2} using the combined-hyperbolic-inverse-power-representation (CHIPR) method and accurate ab initio energies. The calibration grid data is obtained via a general dual-level protocol developed afresh herein that entails both coupled-cluster and multireference configuration interaction energies jointly extrapolated to the complete basis set limit. Such an approach is specially devised to recover much of the spectroscopy from the PES, while still permitting a proper fragmentation of the system to allow for reaction dynamics studies. Besides describing accurately the valence strongly-bound region that includes both the cyclic global minimum and isomerization barriers, the final analytic PES form is shown to properly reproduce dissociation energies,diatomic potentials, and long-range interactions at all asymptotic channels, in addition to naturally reflect the correct permutational symmetry of the potential. Bound vibrational state calculations have been carried out, unveiling an excellent match of the available experimental data on c-SiC_{2}(^{1}A_{1}). To further exploit the global nature of the PES, exploratory quasi-classical trajectory calculations for the endothermic C2 +Si → SiC+C reaction are also performed, yielding thermalized rate coefficients for temperatures up to 5000 K. The results hint for the prominence of this reaction in the innermost layers of the circumstellar envelopes around carbon-rich stars, thence conceivably playing therein a key contribution to the gas-phase formation of SiC, and eventually, solid SiC dust.

Optimized Structural Data at the Complete Basis Set Limit via Successive Quadratic Minimizations

Two variants of a successive quadratic minimization method (SQM and c-SQM) are suggested to calculate the structural properties of molecular systems at the complete basis set (CBS) limit. When applied to H₃⁺, H₂O, CH₂O, SH₂, and SO₂, they revealed CBS/(x₁, x₂) structural parameters that significantly surpass the raw ones calculated at the x₂ basis set level. Such a performance has also been verified for the intricate case of the water dimer. Because the c-SQM method is system specific, thus showing somewhat enhanced results relative to the general SQM protocol, it can be of higher cost depending on the level of calibration used. Yet, it hardly surpasses the general quality of the results obtained with the cost-effective SQM method. Since the number of cycles required to reach convergence is relatively small, both schemes are simple to use and easily adaptable to any of the existing extrapolation schemes for the Hartree-Fock and correlation energies.

Angular-Momentum Extrapolations to the Complete Basis Set Limit: Why and When They Work

The leading L^-3 dependence of the errors in the energies computed with nuclei-centered basis sets comprising functions with angular momenta not exceeding L is rigorously proven for the ¹Σ states of linear molecules and molecular ions with arbitrary even numbers of electrons. This major expansion of the domain of applicability over that offered by the routinely cited Hill asymptotic expression, which is valid only for the helium isoelectronic series, is accomplished with a formalism in which the off-diagonal cusp conditions for the one- and two-electron reduced density matrices play the central role. Despite being provided by these results with theoretical foundations more solid than ever before, the angular-momentum extrapolations to the complete basis set limit appear to work more by happenstance than mathematical rigor due to the poorly predictable variability in the prefactor multiplying the L^-3 term and the far from negligible contributions from the terms involving higher powers of L^-1.

Post-complete-basis-set extrapolation of conventional and explicitly correlated coupled-cluster energies: can the convergence to the CBS limit be diagnosed?

We assess benchmark correlation energies for 130 systems in test sets A24 and TS-106 both with the canonical CCSD(T) and explicitly correlated CCSD(T)-F12 methods. Aiming at enhanced accuracy, the calculated raw energies from both sets are CBS extrapolated to the complete basis set (CBS) limit and subsequently post-CBS extrapolated. Attention is focused at total energies, since their accuracy reflects on that of the interaction energies. Using up to triple-ζ basis sets for CBS and an additional quadruple-ζ for post-CBS, the mean and standard unsigned deviations with canonical CCSD(T) theory are 0.257 ± 0.25 kcal mol^-1, while the corresponding values for CCSD(T)-F12 in its F12a and F12b variants with specialized basis sets up to VQZ-F12 are 0.170 ± 0.13 kcal mol^-1 and 0.048 ± 0.04 kcal mol^-1. Although these show gains at post-CBS level that vary from 0.08 to 0.20 kcal mol^-1 relative to their CCSD(T)/VXZ analogues, the convergence is somewhat less clear when extending the basis up to V5Z-F12, the highest-rung available: 0.220 ± 0.17 kcal mol^-1 and 0.142 ± 0.08 kcal mol^-1, in the same order. An explanation for the up to one order of magnitude smaller deviations in energy differences is detailed. Based on energy differences involving basis set pairs employed for extrapolating to the CBS limit, a convergence diagnostic is also suggested. Arising from irregularities in the basis set that directly correlate with non-dynamical correlation, the new diagnostic may complement popular ones that feature other aspects of correlation.

Canonical and explicitly-correlated coupled cluster correlation energies of sub-kJ mol-1 accuracy via cost-effective hybrid-post-CBS extrapolation

Cost-effectiveness and accuracy are two basic pillars in electronic structure calculations. While cost-effectiveness enhances applicability, high accuracy is sustained when employing advanced computational tools. With the gold standard method of ab initio quantum chemistry at the focal point, canonical CCSD(T) and modern explicitly correlated CCSD(T)-F12 calculations are employed hand in hand to develop accurate hybrid post-CBS extrapolation schemes, which are validated using popular training sets involving a total of 130 molecules. By using raw valence-only calculations at CCSD(T)/VDZ and CCSD(T)/VQZ-F12 levels of theory, the novel scheme leads to the prediction of absolute energies that differ on average (-0.170 ± 0.224) kcal mol-1 from the highest affordable CCSD(T)-F12b/V(Q,5)Z-F12 extrapolations, but only (-0.048 ± 0.228) kcal mol-1 from the post-CBS extrapolated values based on CBS(D,T), CBS(D,Q) and CBS(T,Q) energies. From the cost-effectiveness standpoint, the approach is a kind of pseudo one-point extrapolation scheme since its cost is basically that of the highest-rung raw energy where it is based. Variants that imply no additional cost are also discussed, emerging h-pCBS(dt,dq)ab as the most effective. The approach can also be used with PNO-based local correlation methods that gained popularity due to allowing coupled-cluster calculations even for large molecules at reduced computational cost, namely local PNO-CCSD(T) and PNO-CCSD(T)-F12b. To gauge the approach performance, both the hydrogen molecule and the O-C2H5 torsion path of ethyl-methyl-ether, an extra molecule here considered with presupposed existence in astrophysical objects, are also studied. Additionally, the nonbonding interactions in the A24 test set are revisited per se. The results show that the title approach may be useful in high-accuracy quantum chemistry, with further improvements requiring the inclusion of contributions beyond the theory here employed such as the ones due to relativistic and nonadiabatic effects.

Extrapolation in quantum chemistry: Insights on energetics and reaction dynamics

Since there is no exact solution for problems in physics and chemistry, extrapolation methods may assume a key role in quantitative quantum chemistry. Two topics where it bears considerable impact are addressed, both at the heart of computational quantum chemistry: electronic structure and reaction dynamics. In the first, the problem of extrapolating the energy obtained by solving the electronic Schrödinger equation to the limit of the complete one-electron basis set is addressed. With the uniform-singlet-and-triplet-extrapolation (USTE) scheme at the focal point, the emphasis is on recent updates covering from the energy itself to other molecular properties. The second topic refers to extrapolation of quantum mechanical reactive scattering probabilities from zero total angular momentum to any of the values that it may assume when running quasiclassical trajectories, QCT/QM-[Formula: see text]J. With the extrapolation guided in both cases by physically motivated asymptotic theories, realism is seeked by avoiding unsecure jumps into the unknown. Although, mostly review oriented, a few issues are addressed for the first time here and there. Prospects for future work conclude the overview.

Effective basis set extrapolations for CCSDT, CCSDT(Q), and CCSDTQ correlation energies

It is well established that extrapolating the coupled-cluster single double triple [CCSD and (T)] correlation energies using empirically motivated extrapolation exponents can accelerate the basis set convergence. Here, we consider the extrapolation of coupled-cluster expansion terms beyond the CCSD(T) level to the complete basis set (CBS) limit. We obtain reference CCSDT-CCSD(T) [T₃-(T)], CCSDT(Q)-CCSDT [(Q)], and CCSDTQ-CCSDT(Q) [T₄-(Q)] contributions from cc-pV{5,6}Z extrapolations for a diverse set of 16 first- and second-row systems. We use these basis-set limit results to fit extrapolation exponents in conjunction with the cc-pV{D,T}Z, cc-pV{T,Q}Z, and cc-pV{Q,5}Z basis set pairs. The optimal extrapolation exponents result in noticeable improvements in performance (relative to α = 3.0) in conjunction with the cc-pV{T,Q}Z basis set pair; however, smaller improvements are obtained for the other basis sets. These results confirm that the basis sets and basis set extrapolations used for obtaining post-CCSD(T) components in composite thermochemical theories such as Weizmann-4 and HEAT are sufficiently close to the CBS limit for attaining sub-kJ/mole accuracy. The fitted extrapolation exponents demonstrate that the T₃-(T) correlation component converges more slowly to the CBS limit than the (Q) and T₄ terms. A systematic investigation of the effect of diffuse functions shows that it diminishes (i) in the order T₃-(T) > (Q) > T₄-(Q) and (ii) with the size of the basis set. Importantly, we find that diffuse functions tend to systematically reduce the T₃-(T) contribution but systematically increases the (Q) contribution. Thus, the use of the cc-pVnZ basis sets benefits from a certain degree of error cancellation between these two components.

On the Interplay between Charge-Shift Bonding and Halogen Bonding

The nature of halogen-bond interactions has been analysed from the perspective of the astatine element, which is potentially the strongest halogen-bond donor. Relativistic quantum calculations on complexes formed between halide anions and a series of Y₃ C-X (Y=F to X, X=I, At) halogen-bond donors disclosed unexpected trends, e. g., At₃ C-At revealing a weaker donating ability than I₃ C-I despite a stronger polarizability. All the observed peculiarities have their origin in a specific component of C-Y bonds: the charge-shift bonding. Descriptors of the Quantum Chemical Topology show that the halogen-bond strength can be quantitatively anticipated from the magnitude of charge-shift bonding operating in Y₃ C-X. The charge-shift mechanism weakens the ability of the halogen atom X to engage in halogen bonds. This outcome provides rationales for outlier halogen-bond complexes, which are at variance with the consensus that the halogen-bond strength scales with the polarizability of the halogen atom.

A global CHIPR potential energy surface for ground-state C3H and exploratory dynamics studies of reaction C2 + CH → C3 + H

A full-dimensional global potential-energy surface (PES) is first reported for ground-state doublet C₃H using the combined-hyperbolic-inverse-power-representation (CHIPR) method and accurate ab initio energies extrapolated to the complete basis set limit. The PES is based on a many-body expansion-type development where the two-body and three-body energy terms are from our previously reported analytic potentials for C₂H(²A') and C₃(¹A',³A'), while the effective four-body one is calibrated using an extension of the CHIPR formalism for tetratomics. The final form is shown to accurately reproduce all known stationary structures of the PES, some of which are unreported thus far, and their interconversion pathways. Moreover, it warrants by built-in construction the appropriate permutational symmetry and describes in a physically reasonable manner all long-range features and the correct asymptotic behavior at dissociation. Exploratory quasi-classical trajectory calculations for the reaction C₂ + CH → C₃ + H are also performed, yielding thermalized rate coefficients for temperatures up to 4000 K.

Complete Basis Set Extrapolation of Electronic Correlation Energies Using the Riemann Zeta Function

In this article, we demonstrate the effectiveness of the method of complete basis set (CBS) extrapolation of correlation energies based on the application of the Riemann zeta function. Instead of fitting the results obtained with a systematic sequence of one-electron bases with a certain functional form, an analytic resummation of the missing contributions coming from higher angular momenta, l, is performed. The assumption that these contributions vanish asymptotically as an inverse power of l leads to an expression for the CBS limit given in terms of the zeta function. This result is turned into an extrapolation method that is very easy to use and requires no "empirical" parameters to be optimized. The performance of the method is assessed by comparing the results with very accurate reference data obtained with explicitly correlated theories and with results obtained with standard extrapolation schemes. On average, the errors of the zeta-function extrapolation are several times smaller compared with the conventional schemes employing the same sequence of bases. A recipe for the estimation of the residual extrapolation error is also proposed.

Accurate CHIPR Potential Energy Surface for the Lowest Triplet State of C3

We report the first global ab initio-based potential energy surface (PES) for ground-state triplet C₃(³A') based on accurate energies extrapolated to the complete basis set (CBS) limit, and using the combined-hyperbolic-inverse-power-representation method for the analytical modeling. By relying on a cost-effective CBS(D,T) protocol, we ensure that the final form reproduces all topographical features of the PES, including its cyclic-linear isomerization barrier, with CBS(5,6)-quality. To partially account for the incompleteness of the N-electron basis and other minor effects, the available accurate experimental data on the relevant diatomics were used to obtain direct-fit curves that replace the theoretical ones in the many-body expansion. Besides describing properly long-range interactions at all asymptotic channels and permutational symmetry by built-in construction, the PES reported here reproduces the proper exothermicities at dissociation regions as well as the spectroscopy of the diatomic fragments. Bound vibrational state calculations in both linear and cyclic isomers have also been carried out, unveiling a good match of the available data on C₃(ã ³Π_u), while assisting with IR band positions for C₃(³A₂^') that may serve as a guide for its laboratory and astronomical detection.

CBS extrapolation of Hartree-Fock energy: Pople and Dunning basis sets hand-to-hand on the endeavour

The Hartree-Fock (HF) energy is shown to be extrapolatable from subminimal, minimal, and extended basis sets. Unprecedentedly, it can be reliably extrapolated to the complete basis set limit (CBS) from as small as Pople's STO-2G up to the largest aug-cc-pVXZ basis sets of Dunning's correlation-consistent type and even more complete variants of such ansatzes. The approach is tested on a panoply of small, medium and large molecular systems. As a case study, extrapolation is shown to be particularly reliable and cost-effective (typically an order of magnitude more economical for similar accuracy) when employing a hybrid (STO-2G,extended) basis set pair. The additional cost relative to a single-point calculation with the extended basis set is of no significance. After training the novel scheme with a meager set of 18 systems (TS-18), four test sets covering from neutrals (TS-105) to neutrals plus cations (TS-26) and anions (TS-17), these involving up to third-row elements, as well as tetrapeptide conformers (TS-10) were utilized for its validation. Although the answer to the question of whether all basis sets converge to the same HF/CBS limit is a clear yes, that of whether they show the same convergence rate is open, an issue also examined in the present work. For systems involving atoms up to the second-row of the periodic table the answer seems to be a clear yes, at least when using correlation-consistent basis sets, but it is shown to be a no when the study is extended to anionic systems containing third-row elements. For the latter, augmented correlation consistent basis sets display a faster, perhaps optimal, convergence rate. The implications of this in practice can hardly be anticipated. For example, it cannot be predicted a priori which of such basis sets will yield the best relative energies for those anions, although a logical expectation suggests it to be the augmented ones which turn out to display the fastest convergence. As a blind-test type illustration, the novel scheme is applied to the conformational analysis of 2 tetrapeptides with 5 conformers each, using from subminimal up to extended basis sets as large as cc-pV5Z and aug-cc-pVQZ. The (STO-2G,extended) results so obtained are predicted to be in very good agreement with the best data reported thus far.

Do CCSD and approximate CCSD-F12 variants converge to the same basis set limits? The case of atomization energies

While the title question is a clear "yes" from purely theoretical arguments, the case is less clear for practical calculations with finite (one-particle) basis sets. To shed further light on this issue, the convergence to the basis set limit of CCSD (coupled cluster theory with all single and double excitations) and of different approximate implementations of CCSD-F12 (explicitly correlated CCSD) has been investigated in detail for the W4-17 thermochemical benchmark. Near the CBS ([1-particle] complete basis set) limit, CCSD and CCSD(F12^*) agree to within their respective uncertainties (about ±0.04 kcal/mol) due to residual basis set incompleteness error, but a nontrivial difference remains between CCSD-F12b and CCSD(F12^*), which is roughly proportional to the degree of static correlation. The observed basis set convergence behavior results from the superposition of a rapidly converging, attractive, CCSD[F12]-CCSD-F12b difference (consisting mostly of third-order terms) and a more slowly converging, repulsive, fourth-order difference between CCSD(F12^*) and CCSD[F12]. For accurate thermochemistry, we recommend CCSD(F12^*) over CCSD-F12b if at all possible. There are some indications that the nZaPa family of basis sets exhibits somewhat smoother convergence than the correlation consistent family.

Cn (n=2-4): current status

The major aspects of the C₂, C₃ and C₄ elemental carbon clusters are surveyed. For C₂, a brief analysis of its current status is presented. Regarding C₃, the most recent results obtained in our group are reviewed with emphasis on modelling its potential energy surface which is particularly complicated due to the presence of multiple conical intersections. As for C₄, the most stable isomeric forms of both triplet and singlet spin states and their possible interconversion pathways are examined afresh by means of accurate ab initio calculations. The main strategies for modelling the ground triplet C₄ potential are also discussed. Starting from a truncated cluster expansion and a previously reported DMBE form for C₃, an approximate four-body term is calibrated from the ab initio energies. The final six-dimensional global DMBE form so obtained reproduces all known topographical aspects while providing an accurate description of the C₄ linear-rhombic isomerization pathway. It is therefore commended for both spectroscopic and reaction dynamics studies.This article is part of the theme issue 'Modern theoretical chemistry'.

Straightening the Hierarchical Staircase for Basis Set Extrapolations: A Low-Cost Approach to High-Accuracy Computational Chemistry

Because the one-electron basis set limit is difficult to reach in correlated post-Hartree-Fock ab initio calculations, the low-cost route of using methods that extrapolate to the estimated basis set limit attracts immediate interest. The situation is somewhat more satisfactory at the Hartree-Fock level because numerical calculation of the energy is often affordable at nearly converged basis set levels. Still, extrapolation schemes for the Hartree-Fock energy are addressed here, although the focus is on the more slowly convergent and computationally demanding correlation energy. Because they are frequently based on the gold-standard coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)], correlated calculations are often affordable only with the smallest basis sets, and hence single-level extrapolations from one raw energy could attain maximum usefulness. This possibility is examined. Whenever possible, this review uses raw data from second-order Møller-Plesset perturbation theory, as well as CCSD, CCSD(T), and multireference configuration interaction methods. Inescapably, the emphasis is on work done by the author's research group. Certain issues in need of further research or review are pinpointed.

CBS extrapolation in electronic structure pushed to the end: a revival of minimal and sub-minimal basis sets

The complete basis set (CBS) limit is secluded in calculations of electronic structure, and hence CBS extrapolation draws immediate attention.