Performance of numerical atom‐centered basis sets in the ground‐state correlated calculations of noncovalent interactions: Water and methane dimer cases

Numerical atom‐centered basis sets (orbitals) (NAO) are known for their compactness and rapid convergence in the Hartree–Fock and density‐functional theory (DFT) molecular electronic‐structure calculations. To date, not much is known about the performance of the numerical sets against the well‐studied Gaussian‐type bases in correlated calculations. In this study, one instance of NAO [Blum et al., The Fritz Haber Institute ab initio Molecular Simulations Package (FHI‐aims), 2009] was thoroughly examined in comparison to the correlation‐consistent basis sets in the ground‐state correlated calculations on the hydrogen‐bonded water and dispersion‐dominated methane dimers. It was shown that these NAO demonstrate improved, comparing to the unaugmented correlation‐consistent based, convergence of interaction energies in correlated calculations. However, the present version of NAO constructed in the DFT calculations on covalently‐bound diatomics exhibits enormous basis‐set superposition error (BSSE)—even with the largest bases. Moreover, these basis sets are essentially unable to capture diffuse character of the wave function, necessary for example, for the complete convergence of correlated interaction energies of the weakly‐bound complexes. The problem is usually treated by addition of the external Gaussian diffuse functions to the NAO part, what indeed allows to obtain accurate results. However, the operation increases BSSE with the resulting hybrid basis sets even further and breaks down the initial concept of NAO (i.e., improved compactness) due to the significant increase in their size. These findings clearly point at the need in the alternative strategies for the construction of sufficiently‐delocalized and BSSE‐balanced purely‐numerical bases adapted for correlated calculations, possible ones were outlined here. For comparison with the considered NAOs, a complementary study on the convergence properties of the correlation‐consistent basis sets, with a special emphasis on BSSE, was also performed. Some of its conclusions may represent independent interest. © 2013 Wiley Periodicals, Inc.     

Model studies of the optical rotation,and theoretical determination of its sign for β‐pinene and trans‐pinane

We have carried out extensive studies on the basis set dependence of the calculated specific optical rotation (OR) in molecules at the level of the time–dependent Hartree–Fock and density functional approximations. To reach the limits of the basis set saturation, we have devised an artificial model, the asymmetrically deformed (chiral) methane (CM) molecule. This small system permits to use basis sets which are prohibitively large for real chiral molecules and yet shows all the important features of the basis set dependence of the OR values. The convergence of the OR has been studied with n‐aug‐cc‐pVXZ basis sets of Dunning up to the 6–ζ. In a parallel series of calculations, we have used the recently developed large polarized (LPolX) basis sets. The relatively small LPolX sets have been shown to be competitive to very large n‐aug‐cc‐pVXZ basis sets. The conclusions reached in calculations of OR in CM concerning the usefulness of LPolX basis sets have been further tested on (S)‐methyloxirane and (S)‐fluoro‐oxirane. The smallest set of the LPolX family (LPol–ds) has been found to yield OR values of similar quality as those obtained with much larger Dunning's aug‐cc‐pVQZ basis set. These results have encouraged us to carry out the OR calculations with LPol–ds basis sets for systems as large as β‐pinene and trans‐pinane. In both cases, our calculations have lead to the correct sign of the OR value in these molecules. This makes the relatively small LPol–ds basis sets likely to be useful in OR calculations for large molecules. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Improved efficiency of focal point conformational analysis with truncated correlation consistent basis sets

It has been suggested that the computational cost of correlated ab initio calculations could be reduced efficiently by using truncated basis sets on hydrogen atoms (Mintz et al., J Chem Phys 2004, 121, 5629). We now explore this proposal in the context of conformational analysis of small molecules, such as hydrogen peroxide, dimethyl ether, ethyl methyl ether, formic acid, methyl formate, and several small alcohols. It is found that truncated correlation consistent basis sets that lack certain higher angular momentum functions on hydrogen atoms offer accuracy similar to traditional Dunning's basis sets for conformational analysis. Combination of such basis sets with the basis set extrapolation technique to estimate Hartree-Fock and M?ller-Plesset second order energies provides composite extrapolation model chemistries that are significantly more accurate and faster than analogous single point calculations with traditional correlation consistent basis sets. Root mean square errors of best composite extrapolation model chemistries on the used set of molecules are within 0.03 kcal/mol of traditional focal point conformational energies. The applicability of composite extrapolation methods is illustrated by performing conformational analysis of tert-butanol and cyclohexanol. For comparison, conformational energies calculated with popular molecular mechanics force fields are also given.     

Vibrational frequency scale factors for density functional theory and the polarization consistent basis sets

Calculated harmonic vibrational frequencies systematically deviate from experimental vibrational frequencies. The observed deviation can be corrected by applying a scale factor. Scale factors for: (i) harmonic vibrational frequencies [categorized into low (<1000 cm^?1) and high (>1000 cm^?1)], (ii) vibrational contributions to enthalpy and entropy, and (iii) zero‐point vibrational energies (ZPVEs) have been determined for widely used density functionals in combination with polarization consistent basis sets (pc‐n, n = 0,1,2,3,4). The density functionals include pure functionals (BP86, BPW91, BLYP, HCTH93, PBEPBE), hybrid functionals with Hartree‐Fock exchange (B3LYP, B3P86, B3PW91, PBE1PBE, mPW1K, BH&HLYP), hybrid meta functionals with the kinetic energy density gradient (M05, M06, M05‐2X, M06‐2X), a double hybrid functional with Møller‐Plesset correlation (B2GP‐PLYP), and a dispersion corrected functional (B97‐D). The experimental frequencies for calibration were from 41 organic molecules and the ZPVEs for comparison were from 24 small molecules (diatomics, triatomics). For this family of basis sets, the scale factors for each property are more dependent on the functional selection than on basis set level, and thus allow for a suggested scale factor for each density functional when employing polarization consistent basis sets (pc‐n, n = 1,2,3,4). A separate scale factor is recommended when the un‐polarized basis set, pc‐0, is used in combination with the density functionals. © 2012 Wiley Periodicals, Inc.     

Non‐empirical calculations of NMR indirect carbon–carbon coupling constants: 2. Strained polycarbocycles

High‐level non‐empirical calculations of carbon–carbon spin–spin coupling constants in a series of strained polycarbocycles have been carried out, in excellent agreement with available experimental data. The utmost importance of electronic correlation effects in this case has been demonstrated and it has been shown that the Second‐Order Polarization Propagator Approach (SOPPA) is an adequate method to account for those effects. It has been demonstrated that the most reliable basis sets to calculate J(C,C) at the SOPPA level are the correlation‐consistent basis sets of Dunning and co‐workers augmented with inner core s‐functions or decontracted in their s‐parts. The nature of the unusual bridgehead–bridgehead bonds in bicyclobutane and propellane in terms of s‐characters of bonding hybrids and also the hybridization effects in spiropentane are discussed based on the arguments derived from the current calculations of J(C,C) in the title compounds. The values of the unknown J(C,C) in propellane and spiropentane are predicted with high reliability. Copyright © 2003 John Wiley & Sons, Ltd.     

Ab initio calculations of the Ar–ethane intermolecular potential energy surface using bond function basis sets

The intermolecular potential energy surface (PES) of argon with ethane has been studied by ab initio calculations at the levels of second‐order Møller–Plesset perturbation (MP2) theory and coupled‐cluster theory with single, double, and noniterative triple configurations (CCSD(T)) using a series of augmented correlation‐consistent basis sets. Two sets of bond functions, bf1 (3s3p2d) and bf2 (6s6p4d2f), have been added to the basis sets to show a dramatic and systematic improvement in the convergence of the entire PES. The PES of Ar–ethane is characterized by a global minimum at a near T‐shaped configuration with a well depth of 0.611 kcal mol^?1, a second minimum at a collinear configuration with a well depth of 0.456 kcal mol^?1, and a saddle point connecting the two minima. It is shown that an augmented correlation‐consistent basis set with a set of bond functions, either bf1 or bf2, can effectively produce results equivalent to the next larger augmented correlation‐consistent basis set, that is, aug‐cc‐pVDZ‐bf1 ≈ aug‐cc‐pVTZ, aug‐cc‐pVTZ‐bf1 ≈ aug‐cc‐pVQZ. Very importantly, the use of bond functions improves the PES globally, resulting accurate potential anisotropy. Finally, MP2 method is inadequate for accurate calculations, because it gives a potentially overestimated well depth and, more seriously, a poor potential anisotropy. © 2012 Wiley Periodicals, Inc.     

Adsorption‐induced changes of intramolecular optical transitions: PTCDA/NaCl and PTCDA/KCl

Structural and optical properties of isolated perylene‐3,4,9,10‐tetracarboxylic acid dianhydride molecules adsorbed on (100) oriented NaCl and KCl surfaces were studied theoretically to analyze the recently observed red‐shift of the optical excitation spectrum after adsorption (Müller et al., Phys. Rev. B, 2011, 83, 241203; Paulheim et al. Phys. Chem. Chem. Phys., 2013, 15, 4906). The ground‐state structures were obtained by periodic dispersion‐corrected density functional theory (DFT) calculations. For the excited‐state calculations, nonperiodic time‐dependent DFT methods were applied for a cluster model embedded in point charges. The range‐separated hybrid functional CAM‐B3LYP was used. Correlation‐consistent basis sets were used and the calculated excitation energies were extrapolated to the complete basis set limit. The shift of the first optical excitation energy was analyzed in terms of electronic and geometric contributions. It was found that both the distortion of the molecule due to the interaction with the surface and the electrostatic potential of the surface play an important role. © 2015 Wiley Periodicals, Inc.     

Performance of the elongation method with larger basis sets

The elongation method proposed by Imamura serves as a theoretical model for polymerization processes. It can now be used together with larger basis sets, Hartree–Fock and density functional methods from the Gaussian 94 package with direct self‐consistent field (SCF). This allows electronic structure calculation of elongating clusters with an efficiency superior to full cluster calculations and a precision superior to previous versions of our elongation method. Performance and accuracy compared with full cluster calculations on a regular polymer using the BLYP/6‐31G(d, p) method. Interaction energies of water and hydrogen fluoride polymers of increasing length are compared between HF, BLYP methods and 4‐31G, 6‐31G(d, p) basis sets: Diffuse and polarization functions have a large influence on the interaction energy on both polymers. Local density of states are calculated for different cluster lengths. They are in good agreement with full cluster calculations. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 74: 35–47, 1999     

Revisiting the extrapolation of correlation energies to complete basis set limit

The extrapolation scheme of correlation energy is revisited to evaluate the complete basis set limit from double‐zeta (DZ) and triple‐zeta levels of calculations. The DZ level results are adjusted to the standard asymptotic behavior with respect to the cardinal number, observed at the higher levels of basis sets. Two types of adjusting schemes with effective scaling factors, which recover errors in extrapolations with the DZ level basis set, are examined. The first scheme scales the cardinal number for the DZ level energy, while the second scheme scales the prefactor of the extrapolation function. Systematic assessments on the Gaussian‐3X and Gaussian‐2 test sets reveal that these calibration schemes successfully and drastically reduce errors without additional computational efforts. © 2015 Wiley Periodicals, Inc.     

Compatibility of correlation‐consistent basis sets with a hybrid Hartree–Fock/density functional method

The present study examines the feasibility of combining the correlation‐consistent basis sets developed by Dunning and coworkers with the hybrid Hartree–Fock/density functional method B3LYP. Furthermore, extrapolation to the complete basis set (CBS) limit minimizes errors due to the presence of an incomplete basis set and can act as a rigorous test of the limitations of the B3LYP method. Equilibrium geometries, energies, and harmonic vibrational frequencies were determined for a series of well‐studied, yet computationally challenging, small inorganics and their respective ions. The results were then extrapolated to the CBS limit, where applicable, and compared to experiment. It was found that a union between the hybrid Hartree–Fock/density functional B3LYP method and Dunning's augmented correlation‐consistent basis sets gave results that were comparable to molecular orbital methods that explicitly account for electron correlation. Furthermore, the minimum basis set necessary to attain reasonable results for the systems studied was aug‐cc‐pVTZ. Upgrading to the aug‐cc‐pVQZ level and subsequent extrapolation to the CBS limit further improved the overall agreement with the experiment. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 207–216, 1999     

Estimation of isotropic nuclear magnetic shieldings in the CCSD(T) and MP2 complete basis set limit using affordable correlation calculations

A linear correlation between isotropic nuclear magnetic shielding constants for seven model molecules (CH₂O, H₂O, HF, F₂, HCN, SiH₄ and H₂S) calculated with 37 methods (34 density functionals, RHF, MP2 and CCSD(T)), with affordable pcS‐2 basis set and corresponding complete basis set results, estimated from calculations with the family of polarization‐consistent pcS‐n basis sets is reported. This dependence was also supported by inspection of profiles of deviation between CBS estimated nuclear shieldings and shieldings obtained with the significantly smaller basis sets pcS‐2 and aug‐cc‐pVTZ‐J for the selected set of 37 calculation methods. It was possible to formulate a practical approach of estimating the values of isotropic nuclear magnetic shielding constants at the CCSD(T)/CBS and MP2/CBS levels from affordable CCSD(T)/pcS‐2, MP2/pcS‐2 and DFT/CBS calculations with pcS‐n basis sets. The proposed method leads to a fairly accurate estimation of nuclear magnetic shieldings and considerable saving of computational efforts. Copyright © 2013 John Wiley & Sons, Ltd.     

The application of cholesky decomposition in valence bond calculation

The Cholesky decomposition (CD) technique, used to approximate the two‐electron repulsion integrals (ERIs), is applied to the valence bond self‐consistent field (VBSCF) method. Test calculations on ethylene, C_2nH_2n+2, and C_2nH_4n?2 molecules (n = 1–7) show that the performance of the VBSCF method is much improved using the CD technique, and thus, the integral transformation from basis functions to VB orbitals is no longer the bottleneck in VBSCF calculations. The errors of the CD‐based ERIs and of the total energy are controlled by the CD threshold, for which a value of 10^?6 ensures to control the total energy error within 10^?6 Hartree. © 2016 Wiley Periodicals, Inc.     

Complete basis set prediction of methanol isotropic nuclear magnetic shieldings and indirect nuclear spin–spin coupling constants (SSCC) using polarization‐consistent and XZP basis sets and B3LYP and BHandH density functionals

Efficient B3LYP and BHandH density functionals were used to estimate methanol's nuclear magnetic isotropic shieldings and spin–spin coupling constants in the basis set limit. Polarization‐consistent pcS‐n and pcJ‐n (n = 0, 1, 2, 3 and 4), and segmented contracted XZP, where X = D, T, Q and 5, basis sets were used and the results fitted with simple mathematical formulas. The performance of the methods was assessed from comparison with experiment and higher level calculations. ¹J(CH) and ³J(HH) values were determined from very diluted solutions in deuterochloroform and compared with theoretical predictions. The agreement between complete basis set (CBS) density functional theory (DFT) predicted isotropic shieldings and spin–spin values and experiment was good. The BHandH/pcS‐n methanol shieldings obtained using structures optimized at the same level of theory are approaching the accuracy of the advanced coupled‐cluster‐singles‐doubles‐approximate triples (CCSD(T)) calculations. Copyright © 2009 John Wiley & Sons, Ltd.     

Computer simulation of trifluoromethane properties with ab initio force field

Intermolecular interaction potentials of the trifluoromethane dimer in 15 orientations have been calculated using the Hartree‐Fock (HF) self‐consistent theory and the second‐order Møller‐Plesset (MP2) perturbation theory. Single point energies at important geometries were also calibrated by the coupled cluster with single and double and perturbative triple excitation [CCSD(T)] calculations. We have employed Pople's medium size basis sets [up to 6‐311++G(3df,3pd)] and Dunning's correlation consistent basis sets (up to aug‐cc‐pVQZ). Basis set limit potential values were obtained through well‐studied extrapolation methods. The calculated MP2 potential data were employed to parameterize a 5‐site force field for molecular simulations. We performed molecular dynamics simulations using the constructed ab initio force field and compared the simulation results with experiments. Quantitative agreements for the atom‐wise radial distribution functions and the self‐diffusion coefficients over a wide range of experimental conditions can be obtained, thus validating the ab initio force field without using experimental data a priori. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Consistent gaussian basis sets of double‐ and triple‐zeta valence with polarization quality of the fifth period for solid‐state calculations

Consistent basis sets of double‐ and triple‐zeta valence with polarization quality for the fifth period have been derived for periodic quantum‐chemical solid‐state calculations with the crystalline‐orbital program CRYSTAL. They are an extension of the pob‐TZVP basis sets, and are based on the full‐relativistic effective core potentials (ECPs) of the Stuttgart/Cologne group and on the def2‐SVP and def2‐TZVP valence basis of the Ahlrichs group. We optimized orbital exponents and contraction coefficients to supply robust and stable self‐consistent field (SCF) convergence for a wide range of different compounds. The computed crystal structures are compared to those obtained with standard basis sets available from the CRYSTAL basis set database. For the applied hybrid density functional PW1PW, the average deviations of calculated lattice constants from experimental references are smaller with pob‐DZVP and pob‐TZVP than with standard basis sets. © 2018 Wiley Periodicals, Inc.     

The ORP basis set designed for optical rotation calculations

Details of generation of the optical rotation prediction (ORP) basis set developed for accurate optical rotation (OR) calculations are presented. Specific rotation calculations carried out at the density functional theory (DFT) level for model chiral methane molecule, fluorooxirane, methyloxirane, and dimethylmethylenecyclopropane reveal that the ORP set outperforms larger basis sets, among them the aug‐cc‐pVTZ basis set of Dunning (J. Chem. Phys. 1989, 90, 1007) and the aug‐pc‐2 basis set of Jensen (J. Chem. Phys. 2002, 117, 9234; J. Chem. Theory Comput. 2008, 4, 719). It is shown to be an attractive choice also in the case of larger systems, namely norbornanone, β‐pinene, trans‐pinane, and nopinone. The ORP basis set is further used in OR calculations for 24 other systems, and the results are compared to the aug‐cc‐pVDZ values. Whenever large discrepancies of results are observed, the ORP values are in an excellent agreement with the aug‐cc‐pVTZ results. The ORP basis set enables accurate specific rotation calculations at a reduced cost and thus can be recommended for routine DFT OR calculations, also for large and conformationally flexible molecules. © 2013 Wiley Periodicals, Inc.     

Gas‐Phase Electron‐Diffraction Investigation and Quantum‐Chemical Calculations of the Structure of 1,5‐Dimethylsemibullvalene‐2,4,6,8‐tetracarboxylic Dianhydride

The bridged homotropilidines have been of interest for decades because their molecules offer the potential for homoaromaticity. Although many of these have been shown not to be homoaromatic, the energy differences of the delocalized (homoaromatic) forms and the localized (nonhomoaromatic) ones, and the barriers to the interconversion of the localized forms via a Cope rearrangement, have been found to vary greatly. The title compound is a strong candidate for homoaromaticity, and, since the structures of the possible localized and delocalized forms could differ significantly, we have carried out an electron‐diffraction investigation of it augmented by quantum‐mechanical calculations with different basis sets at several levels of theory. Three models were explored: one representing a localized form of C_s symmetry, one a delocalized form of C_2v symmetry, and one a 2 : 1 mixture of the localized/delocalized forms. Although none of the models could be ruled out, the experimental evidence slightly favors the C_s form. These results are consistent with those from the DFT B3PW91 calculations with basis sets ranging from 6‐31G(d) to cc‐pVTZ, which, surprisingly, predict essentially equal thermally corrected free energies for each. The results are discussed.