Complete Basis Set Limit Values Research Articles

Explicitly-correlated double ionization potentials and double electron attachment equation-of-motion coupled cluster methods

Double ionization and double electron attachment equation-of-motion methods, based on linearly approximated explicitly correlated coupled-cluster singles and doubles [CCSD(F12)] are formulated and implemented. An extension of double electron attachment operator is introduced for proper account of short-range correlation effects in states with two additional electrons. Numerical tests for set of doubly ionized and doubly electron attached states of several molecules have shown a good agreement between obtained explicitly-correlated results and the corresponding complete basis set limit values already at double-ζ level.

Structure and binding energy of the H2S dimer at the CCSD(T) complete basis set limit.

This study presents results for the binding energy and geometry of the H2S dimer which have been computed using Møller-Plesset perturbation theory (MP2, MP4) and coupled cluster (CCSD, CCSD(T)) calculations with basis sets up to aug-cc-pV5Z. Estimates of De, EZPE, Do, and dimer geometry have been obtained at each level of theory by taking advantage of the systematic convergence behavior toward the complete basis set (CBS) limit. The CBS limit binding energy values of De are 1.91 (MP2), 1.75 (MP4), 1.41 (CCSD), and 1.69 kcal/mol (CCSD[T]). The most accurate values for the equilibrium S-S distance rSS (without counterpoise correction) are 4.080 (MP2/aug-cc-pV5Z), 4.131 (MP4/aug-cc-pVQZ), 4.225 (CCSD/aug-cc-pVQZ), and 4.146 Å (CCSD(T)/aug-cc-pVQZ). This study also evaluates the effect of counterpoise correction on the H2S dimer geometry and binding energy. As regards the structure of (H2S)2, MPn, CCSD, and CCSD(T) level values of rSS, obtained by performing geometry optimizations on the counterpoise-corrected potential energy surface, converge systematically to CBS limit values of 4.099 (MP2), 4.146 (MP4), 4.233 (CCSD), and 4.167 Å (CCSD(T)). The corresponding CBS limit values of the equilibrium binding energy De are 1.88 (MP2), 1.76 (MP4), 1.41 (CCSD), and 1.69 kcal/mol (CCSD(T)), the latter in excellent agreement with the measured binding energy value of 1.68 ± 0.02 kcal/mol reported by Ciaffoni et al. [Appl. Phys. B 92, 627 (2008)]. Combining CBS electronic binding energies De with EZPE predicted by CCSD(T) vibrational second-order perturbation theory calculations yields Do = 1.08 kcal/mol, which is around 0.6 kcal/mol smaller than the measured value of 1.7 ± 0.3 kcal/mol. Overall, the results presented here demonstrate that the application of high level calculations, in particular CCSD(T), in combination with augmented correlation consistent basis sets provides valuable insight into the structure and energetics of the hydrogen sulfide dimer.

Developments in stochastic coupled cluster theory: The initiator approximation and application to the uniform electron gas.

We describe further details of the stochastic coupled cluster method and a diagnostic of such calculations, the shoulder height, akin to the plateau found in full configuration interaction quantum Monte Carlo. We describe an initiator modification to stochastic coupled cluster theory and show that initiator calculations can at times be extrapolated to the unbiased limit. We apply this method to the 3D 14-electron uniform electron gas and present complete basis set limit values of the coupled cluster singles and doubles (CCSD) and previously unattainable coupled cluster singles and doubles with perturbative triples (CCSDT) correlation energies for up to r(s) = 2, showing a requirement to include triple excitations to accurately calculate energies at high densities.

Many-Body Basis Set Superposition Effect.

The basis set superposition effect (BSSE) arises in electronic structure calculations of molecular clusters when questions relating to interactions between monomers within the larger cluster are asked. The binding energy, or total energy, of the cluster may be broken down into many smaller subcluster calculations and the energies of these subsystems linearly combined to, hopefully, produce the desired quantity of interest. Unfortunately, BSSE can plague these smaller fragment calculations. In this work, we carefully examine the major sources of error associated with reproducing the binding energy and total energy of a molecular cluster. In order to do so, we decompose these energies in terms of a many-body expansion (MBE), where a "body" here refers to the monomers that make up the cluster. In our analysis, we found it necessary to introduce something we designate here as a many-ghost many-body expansion (MGMBE). The work presented here produces some surprising results, but perhaps the most significant of all is that BSSE effects up to the order of truncation in a MBE of the total energy cancel exactly. In the case of the binding energy, the only BSSE correction terms remaining arise from the removal of the one-body monomer total energies. Nevertheless, our earlier work indicated that BSSE effects continued to remain in the total energy of the cluster up to very high truncation order in the MBE. We show in this work that the vast majority of these high-order many-body effects arise from BSSE associated with the one-body monomer total energies. Also, we found that, remarkably, the complete basis set limit values for the three-body and four-body interactions differed very little from that at the MP2/aug-cc-pVDZ level for the respective subclusters embedded within a larger cluster.

Liquid chloroform structure from computer simulation with a fullab initiointermolecular interaction potential

We have calculated the intermolecular interaction energies of the chloroform dimer in 12 orientations using the second-order Møller-Plesset perturbation theory. Single point energies of important geometries were calibrated by the coupled cluster with single and double and perturbative triple excitation method. Dunning's correlation consistent basis sets up to aug-cc-pVQZ have been employed in extrapolating the interaction energies to the complete basis set limit values. With the ab initio potential data we constructed a 5-site force field model for molecular dynamics simulations. We compared the simulation results with recent experiments and obtained quantitative agreements for the detailed atomwise radial distribution functions. Our results were also consistent with previous results using empirical force fields with polarization effects. Moreover, the calculated diffusion coefficients reproduced the experimental data over a wide range of thermodynamic conditions. To the best of our knowledge, this is the first ab initio force field which is capable of competing with existing empirical force fields for liquid chloroform.

An integrated experimental and quantum-chemical investigation on the vibrational spectra of chlorofluoromethane.

The vibrational analysis of the gas-phase infrared spectra of chlorofluoromethane (CH2ClF, HCFC-31) was carried out in the range 200-6200 cm(-1). The assignment of the absorption features in terms of fundamental, overtone, combination, and hot bands was performed on the medium-resolution (up to 0.2 cm(-1)) Fourier transform infrared spectra. From the absorption cross section spectra accurate values of the integrated band intensities were derived and the global warming potential of this compound was estimated, thus obtaining values of 323, 83, and 42 on a 20-, 100-, and 500-year horizon, respectively. The set of spectroscopic parameters here presented provides the basic data to model the atmospheric behavior of this greenhouse gas. In addition, the obtained vibrational properties were used to benchmark the predictions of state-of-the-art quantum-chemical computational strategies. Extrapolated complete basis set limit values for the equilibrium geometry and harmonic force field were obtained at the coupled-cluster singles and doubles level of theory augmented by a perturbative treatment of triple excitations, CCSD(T), in conjunction with a hierarchical series of correlation-consistent basis sets (cc-pVnZ, with n = T, Q, and 5), taking also into account the core-valence correlation effects and the corrections due to diffuse (aug) functions. To obtain the cubic and quartic semi-diagonal force constants, calculations employing second-order Møller-Plesset perturbation (MP2) theory, the double-hybrid density functional B2PLYP as well as CCSD(T) were performed. For all anharmonic force fields the performances of two different perturbative approaches in computing the vibrational energy levels (i.e., the generalized second order vibrational treatment, GVPT2, and the recently proposed hybrid degeneracy corrected model, HDCPT2) were evaluated and the obtained results allowed us to validate the spectroscopic predictions yielded by the HDCPT2 approach. The predictions of the deperturbed second-order perturbation approach, DVPT2, applied to the computation of infrared intensities beyond the double-harmonic approximation were compared to the accurate experimental values here determined. Anharmonic DFT and MP2 corrections to CCSD(T) intensities led to a very good agreement with the absorption cross section measurements over the whole spectral range here analysed.

A Comparison of One-Particle Basis Set Completeness, Higher-Order Electron Correlation, Relativistic Effects, and Adiabatic Corrections for Spectroscopic Constants of BH, CH+, and NH

To investigate the relative importance of various small sources of error in theoretical predictions of molecular properties, we report spectroscopic constants for the ground electronic states of BH, CH+, and NH, which are nearly converged to the adiabatic ab initio limit. Computations are performed using full configuration interaction and coupled-cluster singles, doubles, and perturbative triples methods with correlation-consistent basis sets of double- to sextuple-ζ quality. The equilibrium bond lengths, re, harmonic vibrational frequencies, ωe, anharmonicity constants, ωexe, centrifugal distortion constants, De, and other quantities are compared with experiment for each species. The systematic dependence of spectroscopic constants on the one-particle basis is used to estimate the complete basis set limit values by using a two-point linear extrapolation scheme. The importance of core correlation, scalar relativistic corrections, higher-order electron correlation, and basis set completeness are carefully...

TiCl, TiH, and TiH + bond energies: a test of a correlation-consistent Ti basis set

Correlation-consistent basis sets are developed for the Ti atom. The polarization functions are optimized for the average of the 3F and 5F states. One series of correlation-consistent basis sets is for 3d and 4s correlation, while the second series includes 3s and 3p correlation as well as 3d and 4s correlation. These basis sets are tested using the Ti 3F–5F separation and the dissociation energies of TiCl X4Φ, TiH X4Φ, and TiH+ X3Φ. The CCSD(T) complete basis set limit values are determined by extrapolation. The Douglas–Kroll approach is used to compute the scalar relativistic effect. Spin-orbit effects are taken from experiment and/or are computed at the CASSCF level. The Ti 3F–5F separation is in excellent agreement with experiment, while the TiCl, TiH, and TiH+ bond energies are in good agreement with experiment. Extrapolation with the valence basis set is consistent with other atoms, while including 3s and 3p correlation appears to make extrapolation more difficult.