Minimal basis sets in calculations of intermolecular interaction energies

Several minimal (7, 3/3) Gaussian basis sets have been used to calculate the energies and some other properties of CH₄ and H₂O. Improved basis sets developed for these molecules have been extended to NH₃ and HF and employed to H₂CO and CH₃OH. Interaction energies between XH_n molecules have been calculated using the old and the new minimal basis sets. The results obtained with the new basis sets are comparable in accuracy to those calculated with significantly more extended basis sets involving polarization functions. Binding energies calculated using the counterpoise method are not much different for the new and the old minimal basis sets, and are likely to be more accurate than the results of much more extended calculations.     

Development and testing of a compact basis set for use in effective core potential calculations on rhodium complexes

We present a set of effective core potential (ECP) basis sets for rhodium atoms which are of reasonable size for use in electronic structure calculations. In these ECP basis sets, the Los Alamos ECP is used to simulate the effect of the core electrons while an optimized set of Gaussian functions, which includes polarization and diffuse functions, is used to describe the valence electrons. These basis sets were optimized to reproduce the ionization energy and electron affinity of atomic rhodium. They were also tested by computing the electronic ground state geometry and harmonic frequencies of [Rh(CO)₂μ‐Cl]₂, Rh(CO)₂ClPy, and RhCO (neutral and its positive, and negative ions) as well as the enthalpy of the reaction of [Rh(CO)₂μ‐Cl]₂ with pyridine (Py) to give Rh(CO)₂ClPy, at different levels of theory. Good agreement with experimental values was obtained. Although the number of basis functions used in our ECP basis sets is smaller than those of other ECP basis sets of comparable quality, we show that the newly developed ECP basis sets provide the flexibility and precision required to reproduce a wide range of chemical and physical properties of rhodium compounds. Therefore, we recommend the use of these compact yet accurate ECP basis sets for electronic structure calculations on molecules involving rhodium atoms. © 2012 Wiley Periodicals, Inc.     

Uniform quality gaussian basis sets for molecular calculations,I. C1 hydrocarbons

The optimality of MO basis sets of Gaussian functions, when constructed from AO basis sets optimized for the neutral atom or for atom ions, is investigated. A formal charge parameter Q is defined and used to adjust the AO basis sets to the molecular environment, by virtue of a simple quadratic expression. Calculations on a series of C₁ hydrocarbons (CH₂, CH₃, CH₃⁺, CH₃^?, CH₄) using 3G basis sets indicate considerable variations in the optimum Q value with the molecular species. The proposed method offers a simple alternative technique to a full molecular basis set optimization.     

Correlation balance in basis sets for electronic structure calculations

The balanced addition of polarization functions to the 6–31G and 6–311G basis sets for correlated wave functions is evaluated using bond energy predictions at the MP 2 and full MP 4 levels as a measure of correlation-balanced basis sets. The homolytic dissociations of the XH bonds in H₂, CH₄, NH₃, H₂O, and HF and the XY bonds in C₂H₆, NH₂NH₂, HOOH, and CH₃OH are used as the basis for the evaluation. It is found that correlation balance is achieved for HH, XH, and XY bonds, particularly at the MP 2 level, only if at least as many polarization sets, and sometimes more, are added to the hydrogens as are added to the heavy atoms.     

K-shell ionization energies in molecules by SCF GF calculations

Energies of CH₄, NH₃, H₂O and C₂H₄ K-ionized molecules are calculated by means of a Group Function method using minimal or near minimal basis sets of STO's. Further results from very large basis sets are reported for CH₄, NH₃, and H₂O. Results seemingly do not suffer the shortcomings of a previous SCF MO treatment.     

A Systemic DFT Study on Several 5<Emphasis Type="Italic">d</Emphasis>-Electron Element Dimers: Hf<Subscript>2</Subscript>, Ta<Subscript>2</Subscript>, Re<Subscript>2</Subscript>, W<Subscript>2</Subscript>, and Hg<Subscript>2</Subscript>

Eleven kinds of density functionals in conjunction with three different basis sets are employed to investigate the homonuclear 5d-electron dimers: Hf₂, Ta₂, Re₂, W₂ and Hg₂. The computed bond lengths, vibrational frequencies and dissociation energies of these molecules are used to compare with available experimental data to find the appropriate combination of functional and basis set. The different functionals and basis sets favor different ground electronic state for Hf₂ and Re₂ molecules, indicating that these two dimers are sensitive to the functionals used. The molecular properties of Hg₂ dimer depend strongly on both functionals and basis sets used. It is found that the BP86 and PBEPBE functionals are generally successful in describing the 5d-electron dimers. For the ground states of these dimers, the bonding patterns are determined by natural bond orbital (NBO) analysis. Natural electron configurations show that the 6s and 5d orbitals in the bonding atoms hybrid with each other for the studied dimers except for Hg₂.     

The influence of basis sets on wave function tails

Wave function tails are analyzed quantitatively by investigating the dependence of exterior electron density (EED ) on basis sets; the EED is defined as the integrated electron density outside the repulsive molecular surface. Ab initio MO calculations with large scale basis sets were performed to establish the benchmark order of EED values for valence orbitals of some simple molecules. It is found that very popular basis sets, such as 4-31G, which are determined by energy optimization, are inferior in describing the wave function tails to some similar size basis sets, such as MIDI -4, which are obtained by least-squares fit to near Hartree-Fock atomic functions. Further the EED values for atomic 2s functions are shown to be unfavorably smaller than those for atomic 2p functions when the same value is used for the exponent α in the GTO basis sets. This indicates that the frequently used constraint α_s = α_p is not appropriate for describing wave function tails with medium-size basis sets. Deficiencies in the energy-optimized basis sets are found to become more serious for molecules including heavier atoms.     

Uniform quality gaussian basis sets for molecular calculations. II. C2 hydrocarbons

This investigation is a continuation of a study on the optimality of MO basis sets of Gaussian functions, when constructed from AO basis sets optimized for the neutral atom or for ions. A formal charge parameter Q is used to adjust AO basis sets to the molecular environment, by virtue of a simple quadratic equation. Calculations are performed on a series of seven C₂ hydrocarbons (C₂H₂, C₂H₄, C₂H₆, C₂H₃⁺ (open), C₂H₃⁺ (bridged), C₂H₅⁺ (bridged), and C₂H₄^? radical anion). A simple rule is formulated to give approximate values of the charge parameter Q.     

NLO‐X (X = I–III): New Gaussian basis sets for prediction of linear and nonlinear electric properties

New adjusted Gaussian basis sets are proposed for first and second rows elements (H, B, C, N, O, F, Si, P, S, and Cl) with the purpose of calculating linear and mainly nonlinear optical (L–NLO) properties for molecules. These basis sets are new generation of Thakkar‐DZ basis sets, which were recontracted and augmented with diffuse and polarization extrabasis functions. Atomic energy and polarizability were used as reference data for fitting the basis sets, which were further applied for prediction of L–NLO properties of diatomic, H₂, N₂, F₂, Cl₂, BH, BF, BCl, HF, HCl, CO, CS, SiO, PN, and polyatomic, CH₄, SiH₄, H₂O, H₂S, NH₃, PH₃, OCS, NNO, and HCN molecules. The results are satisfactory for all electric properties tested; dipole moment (µ), polarizability (α), and first hyperpolarizability (β), with an affordable computational cost. Three new basis sets are presented and called as NLO‐I (ADZP), NLO‐II (DZP), and NLO‐III (VDZP). The NLO‐III is the best choice to predict L–NLO properties of large molecular systems, because it presents a balance between computational cost and accuracy. The average errors for β at B3LYP/NLO‐III level were of 8% for diatomic molecules and 14% for polyatomic molecules that are within the experimental uncertainty. © 2014 Wiley Periodicals, Inc.     

Valence bond studies of the D2h isomer of O4: An interim report

Single-zeta and π-electron double-zeta basis sets are used to examine some theories of the origin of the stability of the D_2h isomer of O₄, using ab initio valence-bond procedures. With these basis sets, resonance between covalent-type (i.e., O₂ ·· O₂) valence-bond structures does not lead to a stabilization of the dimer relative to the separated monomers. When basis sets of the same size are used to construct wave functions for covalent and ionic structures, covalent-ionic resonance (i.e., O₂ ·· O₂ ↔ O₂⁺ ·· O₂⁻ ↔ O₂⁻ ·· O₂⁺) is also unable to stabilize the dimer. Without consideration of the basis-set superposition error, stability is obtained when the size of the AO basis is increased for the dimer relative to the monomer, either via the basis for the ionic structures or by the inclusion of midbond functions. Brief consideration is given to an increased-valence structure for the dimer. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 547–555, 1997     

H2O,H2, HF,F2 and F2O nuclear magnetic shielding constants and indirect nuclear spin–spin coupling constants (SSCCs) in the BHandH/pcJ‐n and BHandH/XZP Kohn–Sham limits

Good performance of segmented contracted basis sets XZP, where X = D, T, Q and 5, for obtaining H₂O, H₂, HF, F₂ and F₂O nuclear isotropic shielding constants in the BHandH Kohn–Sham basis set limit was shown. The results of two‐ and three‐parameter complete basis set limit extrapolation schemes were compared with experimental results, earlier literature data and benchmark ab initio results. Similar convergence patterns of shieldings obtained from calculations using general purpose XZP basis sets and from polarization‐consistent basis sets pcS‐n and pcJ‐n, where n = 0, 1, 2, 3 and 4, designed to accurately predict magnetic properties were observed. On the contrary, the SSCCs were more sensitive to the XZP basis set size and generally less accurate than those estimated using pcJ‐n basis set family. The BHandH density functional markedly outperforms B3LYP method in predicting heavy atom shieldings and SSCCs values in the studied systems. Copyright © 2009 John Wiley & Sons, Ltd.     

A systematic preparation of new contracted Gaussian-type orbital sets. VIII. MINI-1 and MIDI-1 sets for Ga through Cd

Minimal contracted Gaussian basis sets are presented for Ga through Cd. Characteristically these Gaussian-based minimal sets give far better d orbital energies than those by minimal STO basis sets. These new basis sets were tested on Br₂ for which a new benchmark calculation was also performed. The test result is satisfactory in that these basis sets produce good general agreement with the near Hartree–Fock calculation with respect to the molecular spectroscopic constants.     

Ab initio calculations on sulfur-containing compounds. I. Uniform quality basis sets for sulfur: Total energies and geometries of H2S

Uniform quality basis sets (UQ-NG ; N=3, 4, 5), with s = p and s ≠ p, and a 6-31 G^* basis set have been optimized for the sulfur atom. These uniform quality basis sets in their uncontracted and contracted forms were used, together with other basis sets reported in the literature (a total of 40 basis sets), to study their accuracy in predicting the bond length and bond angle of H₂S.     

Basis set dependence of correlation corrections to protonation energies

The protonation energies of NH₃ and H₂O have been computed using a variety of basis sets. It is found that the effect of election correlation on these energies cannot reliably be determined without the use of large (triple-split and polarized) basis sets.     

Substituent effects in second-row molecules: Basis set performance in calculations of normal valency phosphorus and sulfur compounds

Ab inito molecular orbital calculations of the phosphorus- and sulfur-containing series PH₂X, PH₃X⁺, SHX, and SH₂X⁺ (X = H, CH₃, NH₂, OH, F) have been carried out over a range of Gaussian basis sets and the results (optimized geometrical structures, relative energies, and electron distributions) critically compared. As in first-row molecules there are large discrepancies between substituent interaction energies at different basis set levels, particularly in electron-rich molecules; use of basis sets lower than the supplemented 6-31G basis incurs the risk of obtaining substituent stabilizations with large errors, including the wrong sign. Only a small part of the discrepancies is accounted for by structural differences between the optimized geometries. Supplementation of low level basis sets by d functions frequently leads to exaggerated stabilization energies for π-donor substituents. Poor performance also results from the use of split valence basis sets in which the valence shell electron density is too heavily concentrated in diffuse component of the valence shell functions, again likely to occur in electron-rich molecules. Isodesmic reaction energies are much less sensitive to basis set variation, but d function supplementation is necessary to achieve reliable results, suggesting a marginal valence role for d functions, not merely polarization of the bonding density. Optimized molecular geometries are relatively insensitive to basis set and electron population analysis data, for better-than-minimal bases, are uniform to an unexpected degree.     

Comparative study of BSSE correction methods at DFT and MP2 levels of theory

A comparative study of intermolecular potential energy curves is performed on the complexes H₂O(SINGLE BOND)HF, H₂O(SINGLE BOND)H₂O, H₂O(SINGLE BOND)H₂S, and H₂S(SINGLE BOND)H₂S using nine different basis sets at the MP2 and DFT (BLYP and B3LYP) levels of theory. The basis set superposition error is corrected by means of the counterpoise scheme and based on the “chemical Hamiltonian approach.” The counterpoise and CHA-corrected DFT curves are generally close to each other. Using small basis sets, the B3LYP functional cannot be favored against the BLYP one because the BLYP results sometimes get closer to the MP2 values than those of B3LYP. From the results—including the available literature data—we suggest that one has to use at least polarized-valence triple-zeta-quality basis sets (TZV, 6-311G) for the investigation of hydrogen-bonded complexes. Special attention must be paid to the physical nature of the binding. If the dispersion forces become significant DFT methods are not able to describe the interaction. Proper correction for the basis set superposition error is found to be mandatory in all cases. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 575–584, 1998     

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 effect of bond functions on molecular dissociation energies

Multi-reference Cl calculations are reported for the ground states of HCl and N₂ at their equilibrium distances, and for their separated atoms. Basis sets of double-zeta and double-zeta plus polarization quality are systematically augmented by additional sets of functions located at the bond centers. It is shown that use of bond functions can lead to either an underestimate or an overestimate of the the bond energy. Optimum basis sets for each molecule were obtained, giving D_e values of 4.59 eV for HCl (expt. 4.62 eV) and 9.96 eV for N₂ (expt. 9.905 eV) at the estimated full Cl level. The quality of the potential curves obtained with these basis sets is discussed.