Accurate MRCI and CC study of the most relevant stationary points and other topographical attributes for the ground-state C(2)H(2) potential energy surface

Extensive ab initio calculations have been performed to determine the energy, geometry, vibrational frequencies, and relative energetics of all stationary points of the C(2)H(2) ground-state potential-energy surface. The geometries of acetylene and vinylidene minima as well as all transition states are reported at the CASSCF, MRCI, and CCSD(T) levels with aug-cc-pVXZ basis sets. Other more advanced levels of CC theory have also been utilized where judged adequate, mostly for check purposes. Also reported are theoretical limiting values of the energetics of the reaction, deduced from series of computations using the USTE extrapolation method. The data here reported should be valuable for modeling a single-sheeted global potential energy surface for the title system.     

Polarization consistent basis sets. V. The elements Si-Cl

Polarization consistent basis sets, optimized for density functional calculations, are proposed for the elements Si-Cl. Their performance for atomization energies, equilibrium geometries, harmonic vibrational frequencies, and associated infrared intensities is compared with other commonly used basis sets. Atomization energies can be predicted to within 0.01 kJ/mol per atom of the basis set limit by extrapolation of the pc-2, -3, and -4 results. Equilibrium bond distances and harmonic vibrational frequencies can be calculated to within 10(-5) A and 0.5 cm(-1), respectively, of the basis set limit. The pc-n basis sets are shown to give comparable or better accuracy than other alternatives, while containing fewer or equal number of primitive basis functions.     

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.     

Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs

MP2 and CCSD(T) complete basis set (CBS) limit interaction energies and geometries for more than 100 DNA base pairs, amino acid pairs and model complexes are for the first time presented together. Extrapolation to the CBS limit is done by using two-point extrapolation methods and different basis sets (aug-cc-pVDZ - aug-cc-pVTZ, aug-cc-pVTZ - aug-cc-pVQZ, cc-pVTZ - cc-pVQZ) are utilized. The CCSD(T) correction term, determined as a difference between CCSD(T) and MP2 interaction energies, is evaluated with smaller basis sets (6-31G** and cc-pVDZ). Two sets of complex geometries were used, optimized or experimental ones. The JSCH-2005 benchmark set, which is now available to the chemical community, can be used for testing lower-level computational methods. For the first screening the smaller training set (S22) containing 22 model complexes can be recommended. In this case larger basis sets were used for extrapolation to the CBS limit and also CCSD(T) and counterpoise-corrected MP2 optimized geometries were sometimes adopted.     

The s-homodesmotic method for the computation of conventional strain energies of bicyclic systems and individual rings within these systems

The s-homodesmotic method for computing conventional strain energies (CSE) has been extended for the first time to bicyclic systems and to individual rings within these systems. Unique isodesmic, homodesmotic, and hyperhomodesmotic reactions originate from the s-homodesmotic method. These are used to investigate 12 bicyclic systems comprising cyclopropane and cyclobutane and how the CSE of each system compares to the sum of the individual rings within each. Equilibrium geometries, harmonic vibrational frequencies, and the corresponding electronic energies and zero point vibrational energy corrections are computed for all relevant molecules using second-order perturbation theory and density functional theory (B3LYP) with the correlation consistent basis sets cc-pVDZ and cc-pVTZ. Single-point CCSD(T) energies are computed at the MP2/cc-pVTZ optimized geometries to ascertain the importance of higher order correlation effects. Results indicate that CSEs are additive when the two rings are separated by one or two bonds and somewhat additive in other cases.     

Density functional theory studies on the structures and vibrational spectra of 3,6-dichlorocarbazole and 3,6-dibromocarbazole

Becke 3-Lee-Yang-Parr density functional theory (DFT) calculations using 6-311G** and 6-311G(2df,p) basis sets were carried out to study molecular structures and vibrational spectra of 3,6-dichlorocarbazole and 3,6-dibromocarbazole. The optimized geometries, vibrational frequencies, IR intensities, and Raman activities have been obtained. On the basis of B3LYP calculations, a normal mode analysis was performed to assign the vibrational fundamental frequencies according to the potential energy distributions. The computational frequencies are in good agreement with the observed results.     

Approaching the basis set limit for transition metal compounds with highly polar bonds: a benchmark coupled-cluster study of the ScF3 and FeF3 molecular structures and spectra

The molecular equilibrium geometries, quadratic and cubic force constants, vibrational frequencies, and infrared intensities of scandium and iron trifluorides are determined ab initio with a sequence of atomic natural orbital basis sets using the CCSD(T) treatment of electron correlation. The largest basis set of spdf ghi quality contains 462 contracted Gaussian functions. Relativistic corrections are applied to compute the equilibrium geometries and vibrational frequencies. The cubic force constants are used to estimate vibrational corrections to the effective r(g) internuclear distances determined in the gas electron diffraction experiments. The computed molecular properties are extrapolated to the complete basis-set limit. The predicted values are compared to the available experimental data; uncertainties and inconsistencies in these data are then discussed.     

BSSE‐corrected geometry and harmonic and anharmonic vibrational frequencies of formamide–water and formamide–formamide dimers

The basis set superposition error (BSSE) influence in the geometry structure, interaction energies, and intermolecular harmonic and anharmonic vibrational frequencies of cyclic formamide–formamide and formamide–water dimers have been studied using different basis sets (6‐31G, 6‐31G**, 6‐31++G**, D95V, D95V**, and D95V++**). The a posteriori “counterpoise” (CP) correction scheme has been compared with the a priori “chemical Hamiltonian approach” (CHA) both at the Hartree–Fock (HF) and second‐order Møller–Plesset many‐body perturbation (MP2) levels of theory. The effect of BSSE on geometrical parameters, interaction energies, and intermolecular harmonic vibrational frequencies are discussed and compared with the existing experimental data. As expected, the BSSE‐free CP and CHA interaction energies usually show less deep minima than those obtained from the uncorrected methods at both the HF and MP2 levels. Focusing on the correlated level, the amount of BSSE in the intermolecular interaction energies is much larger than that at the HF level, and this effect is also conserved in the values of the force constants and harmonic vibrational frequencies. All these results clearly indicate the importance of the proper BSSE‐free correlation treatment with the well‐defined basis functions. At the same time, the results show a good agreement between the a priori CHA and a posteriori CP correction scheme; this agreement is remarkable in the case of large and well‐balanced basis sets. The anharmonic frequency correction values also show an important BSSE dependence, especially for hydrogen bond stretching and for low frequencies belonging to the intermolecular normal modes. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Singlet Nitric Oxide Dimer Revisited: An Unexpected High Stability Isomer

A new form of the asymmetric cis isomer of the nitric oxide dimer, a cyclic parallelogram, has been found and characterized as the second most stable conformation by ab initio computational methods. Optimum equilibrium geometries, energetics, and harmonic frequencies are computed for the cyclic structure and the more common symmetric cis and trans isomers at the levels of second-order perturbation theory (MP2) and coupled-cluster theory limited to double excitations (CCD) using the 6-311G(2d) and 6-311+G(2d) basis sets. Single-point, more highly correlated calculations are then performed at these CCD optimum geometries. In addition, second-order geometries and frequencies are computed using the larger 6-311G(3df) basis, and single-point MP2 calculations are performed at these geometries using the Aug-cc-pVTZ and Aug-cc-pVQZ basis sets. At the highest levels of theory considered, the cyclic isomer lies within 17 kJ/mol of the symmetric cis global minimum and is more stable than the often-studied symmetric trans conformation.     

Coupled-cluster study of isomers of H2SO2

A theoretical study has been made on six isomers of H2SO2 using coupled-cluster singles and doubles with noniterative triple excitations (CCSD(T)). The isomers studied are sulfoxylic acid (S(OH)2; C2 and Cs conformers), sulfinic acid (HS(=O)OH; 2 C1 conformers), dihydrogen sulfone (H2SO2; C2v), sulfhydryl hydroperoxide (HSOOH; C1), thiadioxirane (Cs), and dihydrogen persulfoxide (H2SOO; Cs). Molecular geometries, harmonic vibrational frequencies, and infrared intensities of all species were obtained using the CCSD(T) method and the 6-311++G(2d,2p) basis set. All aforementioned species were found to be local minima, with the exception of thiadioxirane, which has one imaginary frequency. A prior possible infrared observation of sulfinic acid was reassessed on the basis of the present data. In agreement with previous MP2 results, the present CCSD(T) data provide support for at most 4 of the 8 observed frequencies. The CCSD(T) frequencies and intensities should be of assistance in future identification of H2SO2 isomers by vibrational spectroscopy. Relative energies were calculated using the CCSD(T) method and several larger basis sets. As found previously, the lowest energy species is C2 S(OH)2, followed by Cs S(OH)2, HS(=O)OH, H2SO2, HSOOH, thiadioxirane, and H2SOO. Expanding the basis set significantly reduces the relative energies of HS(=O)OH and H2SO2. The CCSD(T) method was used with extended basis sets (up to aug-cc-pV(Q+d)Z) and basis set extrapolation in two reaction schemes to calculate the DeltaH degrees t (25 degrees C) of C2 S(OH)2. The two reaction schemes gave -285.8 and -282.7 kJ mol-1, which are quite close to a prior theoretical estimate (-290 kJ mol-1).     

THEORETICAL PREDICTIONS OF STRUCTURES AND VIBRATIONAL INFRARED FREQUENCIES. I. MERCAPTANS AND SULFIDES

Abstract

We report the computed equilibrium geometries and vibrational infrared frequencies of a group of thirteen mercaptans and sulfides. The computations were based on the Gaussian 86 Program Package utilizing 3-21G basis sets. The theoretical bond distances and bond angles are in agreement with the available experimental data. The agreement between computed frequencies and available experimental values seems reasonable. We also used the computations as guidelines for the assignment of available experimental infrared frequencies. We believe that it is possible to make reliable frequency predictions by combining computations with available experimental data for groups of similar molecules. However, it is necessary to introduce different correction factors for different types of vibrational modes if we use 3-21G basis sets.     

Ab initio spectroscopic characterization of the HNNO and ONHN radicals

A composite coupled cluster methodology is used with systematic sequences of correlation consistent basis sets to accurately determine the structure, vibrational frequencies, and isotopic shifts for trans-HNNO ((2)A'), cis-HNNO ((2)A'), and ONHN ((2)A'). Anharmonic corrections to the vibrational frequencies and rotational constants are obtained using density functional theory. With basis sets larger than double-zeta, large differences between restricted open-shell Hartree-Fock (ROHF)-based and unrestricted Hartree-Fock (UHF)-based coupled cluster harmonic frequencies are calculated, with the UHF-based ones judged to be more reliable based on an analysis of the orbital hessian eigenvalues. The final calculated anharmonic vibrational band origins are generally in good agreement with the experimental values measured in rare gas matrices. The calculation of the vibrational band origins of the isovalent NO(2) molecule at similar levels of theory exhibits an agreement with experiment to within a few wavenumbers. In the latter case, however, a ROHF treatment was required since the UHF approach failed to provide realistic frequencies for the antisymmetric stretching mode. The heat of formation at 0 K of trans-HNNO is calculated to be 50.5 ± 0.5 kcal∕mol using a very similar composite coupled cluster methodology as in the structure and harmonic frequency determinations.     

Theoretical studies on vibrational spectra of some halides of group IVB elements

The vibrational spectra of group IVB elements halides MX4 (M=Ti(IV), Zr(IV), Hf(II); X=F, Cl, Br and I), have been investigated by ab initio RHF, MP2 and density functional theory B3LYP method with LanL2DZ basis sets. The optimized geometries, calculated vibrational frequencies and Far-IR intensities of MX4 are evaluated via comparison with experimental data. The vibrational frequencies, calculated by these methods, are compared to each other. The results indicate that B3LYP method is more reliable than RHF and MP2 methods for the frequencies calculations for these compounds. With this method, some vibrational frequencies of M2X6(2+)(M=Ti(IV), Zr(IV) and Hf(II); X=F, Cl, Br and I) are also predicted.     

羰基镍簇Ni(CO)n(n=1-4)的结构和Ni-CO键解离性质的密度泛函理论研究

在密度泛函理论框架下, 应用不同泛函计算了配合物Ni(CO)n(n=1~4)的平衡几何构型和振动频率. 考察了泛函和基组重叠误差对预测Ni—CO键解离能的影响. 计算结果表明, 用杂化泛函能得到与实验一致的优化几何构型和较合理的振动频率. 对Ni(CO)n(n=2~4)体系, 用“纯”泛函, 如BP86和BPW91, 可得到与CCSD(T)更符合、 并与实验值接近的解离能. 当解离产物出现单个金属原子或离子(如金属羰基配合物的完全解离)时, BSSE校正项的计算中应保持金属部分的电子结构一致. 只有考虑配体基组和不考虑配体基组两种情况下金属的电子构型与配合物中金属的构型一致时, 才能得到合理的BSSE校正, 从而预测合理的解离能.     

Conventional strain energies of 1,2-dihydroazete, 2,3-dihydroazete, 1,2-dihydrophosphete, and 2,3-dihydrophosphete

The conventional strain energies of 1,2-dihydroazete, 2,3-dihydroazete, 1,2-dihydrophosphete, and 2,3-dihydrophosphete are determined within the isodesmic, homodesmotic, and hyperhomodesmotic models. Optimum equilibrium geometries, harmonic vibrational frequencies, and corresponding electronic energies and zero-point vibrational energies are computed for all pertinent molecular systems using SCF theory, second-order perturbation theory, and density functional theory and employing the correlation consistent basis sets cc-pVDZ, cc-pVTZ, and cc-pVQZ. Single-point fourth-order perturbation theory, CCSD, and CCSD(T) calculations employing the cc-pVTZ and the cc-pVQZ basis sets are computed using the MP2/cc-pVTZ and MP2/cc-pVQZ optimized geometries, respectfully, to ascertain the contribution of higher order correlation. Three DFT functionals, B3LYP, wB97XD, and M06-2X, are employed to determine whether they can yield results similar to those obtained at the CCSD(T) level.     

Ab initio and DFT studies on vibrational spectra of some halides of group IIIB elements

The vibrational spectra of some group IIIB elements halides MX(3) and their dimmers, M(2)X(6) (M=Sc(III), Y(III), La(III); X=F, Cl, Br, I), have been systematically investigated by ab initio restricted Hartree-Fock (RHF) and density functional B3LYP methods with LanL2DZ and SDD basis sets. The optimized geometries and calculated vibrational frequencies are evaluated via comparison with experimental values. The vibrational frequencies, calculated by two methods with different basis sets, are compared to each other. The effect of the methods and the basis sets used on the calculated vibrational frequencies are discussed. Some vibrational frequencies of these complexes are also predicted.     

[M(Im)2X2]型配合物的Far-IR和Raman光谱的理论研究

本文用从头计算RHF和密度泛函B3LYP方法以及LanL2DZ，SDD和6-31G(d)基组计算了配合物M(Im)₂X₂ (Im=imidazole；M=Zn(Ⅱ)，Pd(Ⅱ)，Pt(Ⅱ)；X=F，Cl，Br，I)的几何构型以及Far-IR和Raman振动频率。计算结果表明，对Zn(Ⅱ)配合物而言,B3LYP/6-31G(d)方法得到的几何参数与实验值吻合得最好，B3LYP/SDD次之。在计算Far-IR和Raman振动频率时，发现采用6-31G(d)基组，两种方法计算的结果差别不大。对LanL2DZ和SDD基组而言，对计算结果影响较大的是理论方法，基组影响甚微，个别的振动频率基组影响较大，相比较而言，SDD基组得到的结果更好一些。本文所使用的两种计算方法都能得到与实验值比较吻合的结果，而用从头计算RHF方法计算的结果与实验值更接近一些。在此基础上，预测了Pd(Ⅱ)和Pt(Ⅱ)配合物的Far-IR和Raman振动频率。     

Accuracy and efficiency of atomic basis set methods versus plane wave calculations with ultrasoft pseudopotentials for DNA base molecules

Recent results from Preuss et al. (J Comput Chem 2004, 25, 112) on DNA base molecules, obtained by plane wave density functional calculations using ultrasoft pseudopotentials, are compared with calculations using Gaussian basis sets. Bond lengths and angles agree closely, but dihedral angles and vibrational frequencies show significant differences. The Gaussian basis calculations are at least an order of magnitude more efficient than the plane wave/ultrasoft pseudopotential calculations at a similar level of accuracy; the advantage is even larger if the Fourier Transform Coulomb method is used. To obtain definite benchmark values, the geometries of the four DNA bases were optimized at the MP2 level with large basis sets, up to cc-pVQZ and aug-cc-pVTZ.     

Benchmarking computational chemistry approaches on iminodiacetic acid

In the present study, theoretical harmonic vibrational frequencies (IR and Raman), carbon and proton NMR chemical shifts, geometric parameters, atomic charges (only for heteroatoms), reactivity indices (eLUMO, eHOMO, electronegativity, and hardness), and thermodynamical data (inner energy, enthalpy, Gibbs free energy, and entropy) of iminodiacetic acid (IDA) molecule have been investigated. We utilized ORCA software for B3LYP and HF (combined with Pople and Karlsruhe basis sets) calculations and MOPAC2016 software for semi-empirical calculations (AM1, PM3, and PM6). Theoretical vibrational frequencies and carbon and proton NMR chemical shifts have been compared with the corresponding experimental data. Although there was a strong correlation between the experimental and computational vibrational frequencies at low frequencies (<2200 cm^?1), the computational predictions of vibrational frequencies were unsuccessful at high frequencies (>2200 cm^?1). Distinctly, the studied computational approaches appeared to perform better in the prediction of carbon and proton NMR chemical shifts. Theoretical vibrational frequencies were also compared to each other to understand the impact of method choice (HF vs B3LYP D3 vs semi-empirical methods), dispersion correction (B3LYP D3 vs B3LYP), water solvation (SMD supplemented vs non-supplemented calculations), the family of basis set (Pople vs Karlsruhe basis sets), numbers of zeta (double vs triple zeta), polarization function (polarized vs nonpolarized basis sets), and diffusion function (diffusion supplemented vs non-supplemented basis sets). Moreover, geometric parameters, heteroatom charges, reactivity indices, and thermodynamical data produced by distinct computational approaches, as well, were compared to each other. Based on these comparisons, we detected critical factors (such as water solvation) acting on the computation of geometries, energies, and charges.