On the bond order of C2 and N2

The commonly adopted bond order values of C₂ and N₂ are critically investigated with a new bond order concept. Ab initio calculations with extended basis sets suggest that C₂ can be described by a double to triple bond closer to acetylene than to ethylene and N₂ by a triple bond. The basis set dependence is discussed. Also a relation between the number of basis functions, MO's and non-vanishing eigenvalues of the bond order parts of the density matrix is presented.     

Bond functions for AB initio calculations. MCSCF results for CH,NH, OH and FH

Separate optimized s and p bond functions (BFs) were added to the corresponding Dunning basis sets (11s,6p15s) → (5s,4p13s) for the four hydrides. Properties calculated with these basis sets are quasi-identical to those obtained with conventional polarization functions (d_lp). The computer time ratios found are: t(BF)/t(dlp) = $\text{1}\text{2}$ for UHF calculations, and $\text{2}\text{3}$ for MC SCF calculations.     

Comparison between 6-31 + bond functions and 6-31 G* basis sets in UMP calculations of dissociation energies of the first row element compounds. I

A basis set with bond functions (6-31G + BF) has been tested for its applicability to calculation of dissociation energies of single and multiple bonds by Moeller-Plesset perturbation theory at the second and third orders. Results have been compared with those calculated in the 6-31G* basis set. The 6-31G + BF basis at the MP 2 and MP 3 levels yields better results than 6-31G* basis and the time consumption is less as well. Consideration of the bond functions on the bonds neighboring the bond being broken has no significant influence on the dissociation energies either at the SCF or at the MP 2 levels. If both reactants and products can be characterized by two-center bonds, the 6-31G + BF basis and UMP2 variant of perturbation theory can be recommended for practical calculation of D_e values, especially for the systems where the use of more exact bases is rather difficult.     

Basis set and correlation effects on computed positive ion hydrogen bond energies of the complexes AHn · AHn + 1+1: AHn = NH3, OH2, and FH

Basis set expansion and correlation effects on computed hydrogen bond energies of the positive ion complexes AH_n · AH_{n + 1}⁺¹, for AH_n = NH₃, OH₂ and FH, have been evaluated. The addition of diffuse functions on nonhydrogen atoms is the single most important enhancement of split-valence plus polarization basis sets for computing hydrogen bond energies. Basis set enhancement effects appear to be additive in these systems. The correlation energy contribution to the stabilization energies of these complexes is significant, with the second order term being the largest term and having a stabilizing effect. The third order term is smaller and of opposite sign, while the fourth order term is smaller yet and stabilizing. As a result, computed MP4 stabilization energies are bracketed by the MP2 and MP3 energies. The overall effect of basis set enhancement is to decrease hydrogen bond energies, whereas the addition of electron correlation increases stabilization energies.     

The use of biorthogonal sets in valence bond calculations

The “non-orthogonality problem” in valence bond theory is avoided by using two sets of basis functions, with a biorthogonality property, together with the method of moments. The basis set limit is reached whenall structures are included, but even truncated sets give good (though not monotonic) convergence. Calculations are reported for H₂, LiH, and H₂O.     

Analysis of bonding properties in molecular ground and excited states by a Cohen-type bond order

We propose a Cohen-type bond order analysis in terms of orthogonalized atomic basis functions which can be used to analyze NDO wave functions of large organic and metal–organic molecules. It is shown that for small molecules the results gained with this method are in excellent agreement with the same analysis based on ab initio STO -3G wavefunctions. For large planar aromatic systems these all-valence electron bond orders are found to be a consistent generalization of the π-bond order. A simple relation between these bond orders and the corresponding covalent bond energies is established. The method can be easily extended to study excited state multiconfiguration wave functions. We present calculations for C₂H₂, C₂H₄, C₂H₆, and Mn₂(CO)₁₀. The results indicate that the method can be used to discuss the photochemistry of organic and metal–organic compounds.     

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.     

Molecule‐adapted basis sets optimized with a quantum Monte Carlo method

The Monte Carlo simulated annealing method is adapted to optimize correlated Gaussian‐type functions in nonrelativistic molecular environments. Starting from an atom‐centered atomic Gaussian basis set, the uncontracted functions are reoptimized in the molecular environments corresponding to the H₂O, CN^?, N₂, CO, BF, NO⁺, CO₂, and CS systems. These new molecular adapted basis sets are used to calculate total energies, harmonic vibrational frequencies, and equilibrium geometries at a correlated level of theory. The present methodology is a simple and effective way to improve molecular correlated wave functions, without the need to enlarge the molecular basis set. Additionally, this methodology can be used to generate hierarchical sequences of molecular basis sets with increasing size, which are relevant to establish complete basis set limits. © 2014 Wiley Periodicals, Inc.     

Combined bond polarization function basis sets for accurate ab initio calculation of the dissociation energies of AHn molecules (A=LI to F)

An alternative route toward developing basis sets for post-Hartree-Fock calculations, the hybrid bond polarization function method, is investigated. Two new basis sets, denoted 6-31G(d, p)+ B and 6-31 + G(d,p)+B, are defined for the first-row hydrides. The dissociation energies of the first-row hydride species in their respective ground states are computed using full fourth-order Møller-Plesset theory, and compared with results obtained with large polarized basis sets containing no bond functions. It is shown that results are competitive even with basis sets as large as 6-311++G(3df,3pd), while computation times are reduced by a factor of 4 to 20. On empirical grounds, the basis set superposition error should be neglected entirely.     

Bond functions,covalent potential curves,and the basis set superposition error

In the current practice of quantum chemistry, it is not clear whether corrections for basis set superposition errors should be applied to the calculation of potential energy curves, in order to improve agreement with experimental data. To examine this question, spectroscopic parameters derived from theoretical potential curves are reported for the homonuclear diatomics C₂, N₂, O₂, and F₂, using a configuration interaction method. Three different basis sets were used, including double zeta plus polarization, triple zeta plus double polarization, and double zeta polarization augmented by bond functions. The bond function basis sets, which were optimized in the preceding paper to obtain accurate dissociation energies, also gave the most accurate parameters. The potential curves were then corrected for basis set superposition error using the counterpoise correction, and the spectroscopic parameters were computed again. The BSSE-corrected curves showed worse agreement with experiment for all properties than the original (uncorrected) curves. The reasons for this finding are discussed. In addition to the numerical results, some problems in the application of the BSSE correction to basis sets containing bond functions are shown. In particular, there is an overcounting of the lowering due to the bond functions, regardless of which type of correction is applied. Also, genuine BSSE affects cannot be separated from energy-lowering effects due to basis set incompleteness, and we postulate that it is the latter which is strongly dominant in the calculation of covalent potential curves. Based on these arguments, two conclusions follow: (1) application of BSSE corrections to potential curves should not be routinely applied in situations where the bonding is strong, and (2) appropriate use of bond functions can lead to systematic improvement in the quality of potential curves.     

An improved generator coordinate Hartree–Fock method applied to the choice of contracted Gaussian basis sets for first‐row diatomic molecules

Accurate Gaussian basis sets (18s for Li and Be and 20s11p for the atoms from B to Ne) for the first‐row atoms, generated with an improved generator coordinate Hartree–Fock method, were contracted and enriched with polarization functions. These basis sets were tested for B₂, C₂, BeO, CN⁻, LiF, N₂, CO, BF, NO⁺, O₂, and F₂. At the Hartree–Fock (HP), second‐order Møller–Plesset (MP2), fourth‐order Møller–Plesset (MP4), and density functional theory (DFT) levels, the dipole moments, bond lengths, and harmonic vibrational frequencies were studied, and at the MP2, MP4, and DFT levels, the dissociation energies were evaluated and compared with the corresponding experimental values and with values obtained using other contracted Gaussian basis sets and numerical HF calculations. For all diatomic molecules studied, the differences between our total energies, obtained with the largest contracted basis set [6s5p3d1f], and those calculated with the numerical HF methods were always less than 3.2 mhartree. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 15–23, 2000     

Combined bond-polarization basis sets for accurate determination of dissociation energies. II. Basis set superposition error as a function of the parent basis set

The effect of the parent basis set on the basis set superposition error caused by bond functions is investigated systematically. An important difference between BSSE at the SCF and correlated levels is pointed out. Three new basis sets are defined, denoted 6-311 + G(d,p)B, 6-311 + G(2d,p)B, and 6-311 + G(2df,p)B. BSSE for the first-row hydrides seems to increase uniformly with increasing atomic number of the central atom. Expansion of the valence part of the basis set from 6-31G to 6-311G, as well as adding f functions, has a significant effect on the BSSE. Additional BSSEs incurred by bond functions are less than or equal to 1 kcal/mol for the 6-311 + G(2df,p)B basis set. For the dissociation energies of the first-row hydride species, agreement with experiment within only a few kcal/mol can be obtained even without resorting to isogyric reaction cycles. For high-quality calculations, adding bond functions seems to have definite advantages over expanding the polarization space beyond the [2d1f] level.     

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.     

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.     

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.     

Relation of the force constant of a bond to the electric field at a nucleus

The change in the electric field at a nucleus in a molecule due to bond stretch is related to the force constant of the stretched bond. The validity of this relationship using approximate wave functions at the SCF and MP2 levels of theory is tested for the diatomic molecules H₂, HF, CO, and N₂. The effect of basis set variation on H₂ is also investigated. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1664–1667, 1997     

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.     

Thermochemistry and NBO analysis of peptide bond: Investigation of basis sets and binding energy

Ab initio methods were used to analyze the structure, energetic and binding energy of the five began dipeptides with methionine, Met-Gly, Met-Ala, Met-Ser, Met-Cys, and Met-Thr dipeptides, in gas phase. The structures of the dipetides and involved amino acids in them were optimized by using Hartree-Fock and DFT methods and 3-21G(d), 6-31G(d), 6-311G, 6-311G(d), and 6-311+G(d) basis sets. The effect of basis sets and electron correlations were analyzed with special emphasis on the calculated binding energies and thermodynamic functions. All used methods revealed that Met-Thr has the highest binding energy among all of the five dipeptide molecules. These numerical results suggest that Thr donates the proton easier than other four amino acids and it has the most tendency to join with methionine and it forms the most strong bond with methionine. This fact may be the reason behind the obtained high binding energies for Met-Thr at all levels. From comparison of the values of binding energy for dipeptides in different levels of theory, we could identify that the order of tendency for joint with methionine is Thr > Gly > Ala > Cys > Ser. Also, these data represented that the highest binding energy provide in HF/6-311G level for all of the dipeptides (14.4202, 11.2387, 8.3267, 9.8853, 17.3362 kcal mol⁻¹ for dipeptides 1–5, respectively). Moreover, natural bond orbital (NBO) analysis demonstrated that the effect of basis sets and electron correlations on σ_N1-C2 bonding orbital occupancy is the same as the basis set and electron correlation effects on binding energy of dipeptides in all cases. The obtained results from studying the effect of basis sets and electron correlations on binding energy, NMR and NBO properties showed that the effect of basis sets is almost independent of molecular structure and computational method, while electron correlation effects are relatively dependent to molecular structure and basis set type. In investigating the effect of basis sets and electron correlations on binding properties, the NBO results are in good agreement with the energetic and thermochemistry data at all levels of calculations. The article is published in the original.     

Bond critical points in the electronic structures of binary hydrides

The bond critical points of the binary hydrides formed by the elements of the first two rows of the periodic table have been calculated. Particular attention has been paid to the basis-set dependence of the bond critical points at the experimental equilibrium geometries, or where necessary at model geometries. With the exception of H₂S, stepwise extension of the basis set leads to a smooth convergence of the bond critical points to a set of values which appear to converge to the Hartree–Fock limit. For H₂S it is shownb that the position of the bodn critical point is not only more sensitive to the presence of polarization functions in the basis set, but depends strongly on the orbital exponents of the polarization functions. Extensive optimizations of the exponents of the polarization functions have been carried out with the (12s9p/5s) basis set for second-row hydrides. The effects of contracting the Huzinaga basis sets have been examined.