Medium-size polarized basis sets for high-level-correlated calculations of molecular electric properties

Summary The basis set polarization method is used to derive the first-order polarized basis sets for Ge through Br for calculations of atomic and molecular electric properties. The performance of the [15.12.9/9.7.4] GTO/CGTO basis sets generated in this study is verified in calculations of atomic dipole polarizabilities and dipole moments and polarizabilities of the third-row atom hydrides. Whenever accurate reference data are available for comparison, the excellent performance of the derived first-order polarized basis sets is demonstrated. The role of the core polarization and relativistic contributions to atomic and molecular is also investigated. The detailed basis set data for Ge through Br are given in Appendix.     

Medium-size polarized basis sets for high-level-correlated calculations of molecular electric properties

Summary The first-order polarized basis sets for the use in high-level-correlated investigations of molecular electric properties have been generated for Pb, Bi, Po, and At. The performance of the standard [10.17.14.5/13.11.8.2] and extended [20.17.14.9/13.11.8.4] basis sets has been examined in nonrelativistic and quasirelativistic calculations for atoms and simple closed-shell hydrides. The relativistic contributions to electric dipole properties of those systems have been evaluated by using the recently developed quasirelativistic scheme. The predicted dipole polarizability of At is in good agreement with the results of other relativistic calculations. The calculated quasirelativistic dipole moments of BiH₃ (–0.499 a.u.), PoH₂ (–0.207 a.u.), and AtH (+0.036 a.u.) involve a significant relativistic contribution which amounts to —0.230 a.u., –0.177 a.u., and –0.097 a.u., respectively. The basis set details append this paper. They are also available as a part of the basis set library of the MOLCAS system.     

Reduced-size polarized basis sets for calculations of molecular electric properties. I. The basis set generation

Following the recent studies of basis sets explicitly dependent on oscillatory external electric field we have investigated the possibility of some further truncation of the so-called polarized basis sets without any major deterioration of the computed data for molecular dipole moments, dipole polarizabilities, and related electric properties of molecules. It has been found that basis sets of contracted Gaussian functions of the form [3s1p] for H and [4s3p1d] for the first-row atoms can satisfy this requirement with particular choice of contractions in their polarization part. With m denoting the number of primitive GTOs in the contracted polarization function, the basis sets devised in this article will be referred to as the ZmPol sets. In comparison with earlier, medium-size polarized basis sets (PolX), these new ZmPol basis sets are reduced by 2/3 in their size and lead to the order of magnitude computing time savings for large molecules. Simultaneously, the dipole moment and polarizability data remain at almost the same level of accuracy as in the case of the PolX sets. Among a variety of possible applications in computational chemistry, the ZmPolX are also to be used for calculations of frequencies and intensities in the Raman spectra of large organic molecules (see Part II, this issue).     

Medium-size polarized basis sets for high-level-correlated calculations of molecular electric properties

The basis set polarization method is used to generate the first-order polarized basis sets for Sn, Sb, Te, and I. The standard (spd) and extended (spdf) versions of those basis sets are derived for the purpose of calculations of dipole moments and dipole polarizabilities for molecules involving the fourth-row atoms. The verification of the performance of the generated polarized basis sets is achieved mainly by a cross-examination of different atomic and molecular results calculated in this paper. The role of the core-polarization and relativistic effects is investigated. It is shown that the relativistic contribution to dipole moments of the fourth-row hydrides is commensurate with the contribution due to electron correlation and must be explicitly considered in accurate calculations. The detailed basis set data for Sn through I are presented in the Appendix.     

Reduced-size polarized basis sets for calculations of molecular electric properties. II. Simulation of the Raman spectra

The accuracies of the calculated vibrational frequencies and Raman intensities given by two new, highly compact Pol-type basis sets, Z2PolX and Z3PolX, have been determined and compared to the 6-31G(d), PolX, and aug-cc-pVTZ basis sets. Calculation of accurate Raman intensities has previously required large basis sets, but the ZmPolX basis sets are smaller even than PolX, which are the most compact basis sets able to calculate accurate Raman intensities. For the largest compound studied, C5H10O2, Z3PolX required more than an order of magnitude less CPU time than PolX, which has been shown to be 10 times faster than aug-cc-pVTZ. Two sets of test molecules were studied: one was a series of small molecules for which experimental values for absolute Raman activities were available; the second was a series of medium-sized molecules (mainly common organic solvents) where only relative Raman band intensities were available. The accuracies of the Raman intensities given by both of the ZmPolX basis sets were good compared to those of the PolX and aug-cc-pVTZ sets, and much better than the 6-31G(d) values. The errors in even unscaled frequency values <2000 cm(-1) were also acceptable and were slightly lower for Z3PolX than Z2PolX (30 cm(-1) vs. 48 cm(-1)). The combination of good intensity and frequency data meant that for the medium-sized organic molecules there was a close correspondence between the simulated Raman spectra and experimental data, and that the observed bands could easily be assigned on the basis of these calculations. Achieving this level of accuracy in the simulations at modest computational cost should now allow computational methods to be combined with experimental Raman studies much more widely than is currently the case.     

Reduced-size polarized basis sets for calculations of molecular electric properties. IV. First-row transition metals

Recent studies of the perturbation-dependent basis sets have indicated the possibility of a significant reduction of the size of the usual CGTO sets without considerable loss of accuracy in calculations of molecular electric properties. The resulting (ZPolX) basis sets have been developed for several atoms of the first and second row of the Periodic Table. The same method of the ZPolX basis set generation is extended for the first-row transition metals and the corresponding contracted ZPolX basis sets of the size [6s5p3d1f] are determined for both nonrelativistic and scalar relativistic calculations. The performance of the ZPolX basis sets is verified in calculations on the first-row transition metal oxides at the level of the ROHF, ROHF/CASPT2, and ROHF/CCSD(T) approximations. Also the study of the dipole polarizability of TiCl₄ confirms the excellent features of these very compact basis sets. The ZPolX basis sets for nonrelativistic and relativistic calculations of molecular electric properties are available on the web page http://www.chem.uni.torun.pl/zchk/basis-sets.html.     

Efficient polarized basis sets for evaluation of static and dynamic molecular electric properties

New medium size Gaussian‐type basis set R‐ORP for evaluation of static and dynamic electric properties in molecular systems is presented. It is obtained in a close resemblance to the original ORP basis set, from the source basis set through addition of two first‐order polarization functions whose exponent values are optimized with respect to the finite field restricted open‐shell Hartree–Fock (ROHF) atomic polarizabilities. As the source set the VTZ basis set of Ahlrichs and coworkers, augmented with additional diffuse functions and contracted to the form [6s/3s] for hydrogen and [11s7p/4s3p] for carbon through fluorine, is chosen. The resulting basis set is of the form [6s2p/3s2p] for hydrogen and [11s7p2d/4s3p2d] for other atoms. Presented basis set is next tested in the CCSD static and dynamic molecular polarizability and hyperpolarizability calculations for a set of ten and four test molecules, respectively, for which very accurate reference data exist. Additionally, the recently developed ORP basis set is employed in the calculations to examine the limits of its applicability. Results are compared to the literature data obtained in both, large and diffuse, as well as reduced‐size basis sets. In the case of polarizability calculations, the aug‐pc‐1 and R‐ORP are the optimal choices among the investigated smaller basis sets, with the overall performance of the aug‐pc‐1 set being better. Among the larger sets, the ORP performs better in the case of average polarizability, while the RMSE values for polarizability anisotropy are practically identical for d‐aug‐cc‐pVDZ and ORP sets. Finally, the R‐ORP and ORP basis sets compete other small bases in the evaluation of the first hyperpolarizability in investigated systems. © 2016 Wiley Periodicals, Inc.     

Optimization of Gaussian basis sets for Dirac-Hartree-Fock calculations

We investigate the optimization of Gaussian basis sets for relativistic calculations within the framework of the restricted Dirac-Hartree-Fock (DHF) method for atoms. We compare results for Rn of nonrelativistic and relativistic basis set optimizations with a finite nuclear-size. Optimization of separate sets for each spin-orbit component shows that the basis set demands for the lower j component are greater than for the higher j component. In particular, the p _1/2 set requires almost as many functions as the s _1/2 set. This implies that for the development of basis sets for heavy atoms, the symmetry type for which a given number of functions is selected should be based on j, not on l, as has been the case in most molecular calculations performed to date.     

Standardized basis sets for high-level-correlated relativistic calculations of atomic and molecular electric properties in the spin-averaged Douglas-Kroll (no-pair) approximation I. Groups Ib and IIb

Summary The technique developed earlier for the generation of the so-called first-order polarized basis sets for accurate non-relativistic calculations of molecular electric properties is used to obtain similar basis sets suitable for calculations in the Douglas-Kroll no-pair approximation. The corresponding (relativistic) basis sets are devised for atoms of the Groups Ib and IIb of the periodic table and tested in calculations of atomic polarizabilities and dipole moments of the coinage metal hydrides. Excellent performance of these basis sets has been found in the case of molecular calculations.     

Polarized basis sets and the calculation of infrared intensities from nuclear electric shielding tensors

The idea of the basis set polarization which follows from the known dependence of basis set functions on the perturbation strength is applied to the calculation of the dipole moment derivatives with respect to nuclear displacements. The differentiation of the dipole moment function is replaced by the straightforward evaluation of derivatives of the intramolecular electric field with respect to the external electric field strength. The method and its efficiency are illustrated by a series of calculations of the dipole moment derivatives for the water molecule. Already a polarized basis set of 26 CGTO's derived from the minimal CGTO basis set provides fairly reasonable results.     

Polarized basis sets for accurate calculations of static and dynamic electric properties of molecules

We report on the development and testing of large polarized basis sets (LPolX, where X is the element symbol) for accurate calculations of linear and nonlinear electric properties of molecules. The method used to generate LPolX sets is based on our studies of the analytic dependence of Gaussian functions on external time‐independent and time‐dependent electric fields. At variance with the earlier investigations of small, highly compact (ZPolX) basis sets for moderately accurate calculations of electric properties of large molecules, the present goal is to obtain basis sets that are nearly saturated with respect to the selected class of electric properties and can be used for accurate studies of interaction‐induced properties. This saturation makes the LPolX sets also useful in calculations of optical properties for chiral molecules. In this article, the LPolX sets are generated for X = H, C, N, O, and F, and examined in calculations of linear and nonlinear electric properties of four standard test systems: HF, N₂, CO, and HCN. The study of the performance of LPolX basis sets has been carried out at different levels of approximation ranging from the SCF HF method to highly correlated CCSD(T) approach. The results obtained in this study compare favorably with accurate reference data and show a high level of saturation of LPolX basis sets with respect to the polarization effect due to external electric fields. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

An augmented effective core potential basis set for the calculation of molecular polarizabilities

Calculations of molecular polarizabilities require basis sets capable of accurately describing the responses of the electrons to an external perturbation. Unfortunately, basis sets that yield suitable quantitative results have traditionally been all-electron sets with large numbers of primitives, making their use computationally intractable even for moderately sized systems. We present a systematic augmentation of the effective core potential basis set of Stevens et al. [J Chem Phys 81, 12 (1984), Can J Chem 70, 612 (1992)] for 39 main group elements based on the procedure used to construct diffuse and polarization functions in the well-known Sadlej basis sets [Collec Czech Chem Comm 53, 1995 (1988)]. Representative calculations have been performed and we have shown that results to within 1% of all-electron calculations using the Sadlej basis set can be obtained for <1-35% of the computational cost using this new basis set.     

New basis sets for the evaluation of interaction‐induced electric properties in hydrogen‐bonded complexes

Interaction‐induced static electric properties, that is, dipole moment, polarizability, and first hyperpolarizability, of the CO? (HF)_n and N₂? (HF)_n, n = 1–9 hydrogen‐bonded complexes are evaluated within the finite field approach using the Hartree–Fock, density functional theory, Møller–Plesset second‐order perturbation theory, and coupled cluster methods, and the LPol‐n (n = ds, dl, fs, fl) basis sets. To compare the performance of the different methods with respect to the increase of the complex size, we consider as model systems linear chains of the complexes. We analyze the results in terms of the many‐body and cooperative effects. © 2012 Wiley Periodicals, Inc.