Comparison of the experimental, semi-experimental and ab initio equilibrium structures of acetylene: influence of relativisitic effects and of the diagonal Born-Oppenheimer corrections

The equilibrium structure of acetylene (also named ethyne) has been reinvestigated to resolve the small discrepancies noted between different determinations. The size of the system as well as the large amount of available experimental data provides the quite unique opportunity to check the magnitude and relevance of various contributions to equilibrium structure as well as to verify the accuracy of experimental results. With respect to pure theoretical investigation, quantum-chemical calculations at the coupled-cluster level have been employed together with extrapolation to the basis set limit, consideration of higher excitations in the cluster operator, inclusion of core correlation effects as well as relativistic and diagonal Born-Oppenheimer corrections. In particular, it is found that the extrapolation to the complete basis set limit, the inclusion of higher excitations in the electronic-correlation treatment and the relativistic corrections are of the same order of magnitude. It also appears that a basis set as large as a core-valence quintuple-zeta set is required for accurately accounting for the inner-shell correlation contribution. From a pure experimental point of view, the equilibrium structure has been determined using very accurate rotational constants recently obtained by a "global analysis" (that is to say that all non-negligible interactions are explicitly included in the Hamiltonian matrix) of rovibrational spectra. Finally, a semi-experimental equilibrium structure (where the equilibrium rotational constants are obtained from the experimental ground state rotational constants and computed rovibrational corrections) has been obtained from the available experimental ground-state rotational constants for ten isotopic species corrected for computed vibrational corrections. Such a determination led to the revision of the ground-state rotational constants of two isotopologues, thus showing that structural determination is a good method to identify errors in experimental rotational constants. The three structures are found in a very good agreement, and our recommended values are r(CC) = 120.2958(7) pm and r(CH) = 106.164(1) pm.     

Equilibrium structure and spectroscopic constants of maleic anhydride

The quadratic, cubic, and semi-diagonal quartic force fields of maleic anhydride have been calculated at the MP2 level of theory employing the cc-pVTZ basis set. The spectroscopic constants derived from the force field are in excellent agreement with the corresponding experimental values. The semi-experimental equilibrium structure has been derived from experimental ground state rotational constants and rovibrational corrections calculated from the cubic force field. This semi-experimental equilibrium structure is in excellent agreement with the ab initio structures computed at the CCSD(T) level of theory and it is closer to the ab initio structure than the purely experimental (or empirical) structures r ₀, r _m⁽¹⁾, and r _m⁽²⁾ obtained by microwave spectroscopy as well as the equilibrium structure derived from gas-phase electron diffraction data.     

The case of the weak N-X bond: Ab initio, semi-experimental, and experimental equilibrium structures of XNO (X = H, F, Cl, OH) and FNO2

The equilibrium structures of FNO, ClNO, HONO, and FNO2 have been determined using three different, somewhat complementary methods: a completely experimental, a semi-experimental (where the equilibrium rotational constants are derived from the experimental effective ground-state rotational constants and an ab initio cubic force field), and an ab initio, where geometry optimizations are usually performed at the coupled cluster level of nonrelativistic electronic structure theory using small to very large Gaussian basis sets. For the sake of comparison, the equilibrium structures of HNO and N2O have also been redetermined, confirming and extending earlier results. The semi-experimental method gives structural parameters in good agreement with the reliable experimental results for each compound investigated. Because of inadequate treatment of electron correlation, the single-reference CCSD(T) method gives N-X (X[double bond]F, Cl, OH) bonds that are too strong and associate bond lengths that are significantly too short. The discrepancy increases with increase in the size of the basis set. A much more elaborate treatment of electron correlation at the CCSDTQ level solves this problem and results in increased bond lengths, correctly representing the weakness of the N-X bond in these XNO and XNO2 species. The equilibrium structures determined are accurate to better than 0.001 A and 0.1 degrees .     

Molecular structure of thiourea

The molecular structure of thiourea has been investigated under C(s), C(2), and C(2v) symmetry constraints. At the coupled-cluster level in conjunction with a triple-ζ basis set, only the C(2) conformer has been found to be a real minimum on the potential energy surface. Its equilibrium structure has therefore been accurately evaluated using both theoretical and experimental data. With respect to the former, high-level quantum-chemical calculations at the coupled-cluster level in conjunction with correlation-consistent basis sets ranging in size from triple- to quintuple-zeta have been carried out. Extrapolation to the complete basis-set limit as well as core-correlation effects and inclusion of full treatment of triple excitations in the cluster operator have been considered. On the basis of the vibrational ground-state rotational constants available for five isotopic species and the corresponding computed vibrational corrections, the semiexperimental equilibrium geometry of thiourea has also been determined for the first time.     

Equilibrium structure and torsional barrier of BH3NH3

Born-Oppenheimer equilibrium structures, r(e)(BO), of the electronic ground state of the borazane (BH3NH3) molecule of C3v point-group symmetry are computed ab initio using the CCSD(T) method with basis sets up to quintuple-zeta quality. Inclusion of the counterpoise correction and extrapolation of the structural parameters to the complete basis set limit yield a best estimate of r(e)(BO) of BH3NH3. The anharmonic force field of BH3NH3, computed at the CCSD(T) level of theory with a basis set of triple-zeta quality, allows the determination of semi-experimental equilibrium rotational constants, which in turn result in a semi-experimental equilibrium structure, r(e)(SE). The r(e)(BO) and r(e)(SE) structures are in excellent agreement, indicating the validity of the methods used for their determination. The empirical mass-dependent structure, r(m)(1), of BH3NH3 is also determined. Although it is inferior in quality to the previous two structures, it is much more accurate than the standard empirical r0 and r(s) structures reported earlier for BH3NH3. The semi-experimental r(e)(SE) as well as the empirical r(m)(1) structures determined are based on experimental ground-state rotational constants available from the literature for nine isotopologues of borazane. The effective barrier to the internal rotation of BH3NH3, a molecule isoelectronic with CH3CH3, has been computed ab initio, employing the focal-point analysis (FPA) approach, to be 699 +/- 11 cm(-1). This compares favorably with an empirical redetermination of the effective barrier based on the above r(e)(SE) structure, V3 = 718(17) cm(-1).     

A systematic ab initio study of the equilibrium geometry and vibrational wave numbers of bismuthine

The equilibrium structure and the harmonic and anharmonic force fields of BiH(3) are determined by high-level ab initio calculations using a variety of correlation treatments, basis sets, and pseudopotentials, partly in combination with core polarization potentials. Spin-orbit effects are included by a configuration interaction treatment. This systematic study serves to establish a reliable computational protocol for such calculations and, in particular, to minimize basis set superposition errors through an improved new basis set and/or counterpoise corrections. Using the recommended procedures, the best ab initio results for the equilibrium geometry and the fundamental vibrational wave numbers are in good agreement with the available experimental data, which further supports the recent spectroscopic identification of BiH(3). The ground-state total atomization energy of BiH(3) is predicted to be 153.1 kcal/mol.     

硝胺及其甲基衍生物的从头计算研究 I: 平衡几何构型的全优化

我们用移植的TEXAS梯度法从头计算程序, 选取STO-4-21G基组, 对硝胺、甲硝胺和二甲硝胺等分子的平衡几何构型进行了全优化计算。其C-N, N-H和C-H等键长的计算结果以4-21G基组的标准校正值校正之后, 所得理论预测值与实验数据良好相符, 对N-N,N-O键长以及键角、两面角等构型参数, 对照计算和实验结果进行了讨论。     

Determination of the “Privileged Structure” of 8-Hydroxyquinoline

An accurate semi-experimental equilibrium structure of 8-hydroxyquinoline (8-HQ) has been determined combining experiment and theory. The cm-wave rotational spectrum of 8-HQ was recorded in a pulsed supersonic jet using broadband dual-path reflection and narrowband Fabry-Perot-type resonator Fourier-transform microwave spectrometers. Accurate rotational and quartic centrifugal distortion constants and ¹⁴N quadrupole coupling constants are determined. Rotational constants of all ¹³C, ¹⁸O and ¹⁵N singly substituted isotopologues in natural abundance and those of a chemically synthesized OD isotopologue were used to obtain geometric parameters for all the heavy atoms and the hydroxyl hydrogen from a number of structure determination models. Theoretical approaches allowed for the determination of a semi-experimental equilibrium structure, in which computed rovibrational and electronic corrections were utilized to convert vibrational ground state constants into equilibrium constants. Despite the molecule having only a horizontal plane of symmetry and possessing 11 individual heavy atoms, microwave spectroscopy has allowed for a reliable and accurate structure determination. A mass dependent, structure was determined and proved to be equally reliable by comparison with the B3LYP-D3(BJ)/aVTZ equilibrium structure.     

Determination of Molecular Structure in Terms of Potential Energy Functions from Gas-Phase Electron Diffraction Supplemented by Other Experimental and Computational Data

Different approaches of equilibrium structure determinations by the gas-phase electron diffraction (ED) method or by its combination with other relevant techniques have been reviewed. Some problems and limitations of these approaches are discussed. Special attention is paid to various potential energy function models. Different types of equilibrium bond lengths obtained by the optimization of ED data or their combination with experimental and computational spectroscopic data are compared in tables. Relations between different types of vibrational corrections are discussed. Structure data determined by other methods or approaches are given for comparison.     

On basis set superposition error corrected stabilization energies for large n-body clusters

In this contribution, we propose an approximate basis set superposition error (BSSE) correction scheme for the site-site function counterpoise and for the Valiron-Mayer function counterpoise correction of second order to account for the basis set superposition error in clusters with a large number of subunits. The accuracy of the proposed scheme has been investigated for a water cluster series at the CCSD(T), CCSD, MP2, and self-consistent field levels of theory using Dunning's correlation consistent basis sets. The BSSE corrected stabilization energies for a series of water clusters are presented. A study regarding the possible savings with respect to computational resources has been carried out as well as a monitoring of the basis set dependence of the approximate BSSE corrections.     

Basis set limit CCSD(T) harmonic vibrational frequencies

Benchmark, frozen-core CCSD(T) equilibrium harmonic vibrational frequencies of 12 closed-shell and five open-shell molecules are computed to within 1 cm-1 of the basis set limit using the explicitly correlated CCSD(T)-R12 method. The convergence of the standard CCSD(T) method with the one-particle basis sets of Dunning and co-workers is examined and found to be slow, with mean and maximum absolute errors of 1.3 and 3.5 cm-1 remaining at the cc-pV6Z level. Finite basis set effects do not appear to introduce systematic errors in equilibrium harmonic frequencies, and mean absolute errors reduce by a factor of 2 for each basis set cardinal number increment. The convergence of individual equilibrium harmonic frequencies is not guaranteed to be monotonic due to the associated shift in the equilibrium structure. The inclusion of computed scalar relativistic effects and previously available corrections for core-valence correlation and higher-order excitations in the cluster operator results in an agreement with experimentally derived harmonic frequencies of 0.1, 0.3, and -0.4 cm-1 for HF, N2, and CO, respectively. F2 continues to present a challenge to computational chemistry with an error of 3.2 cm-1, primarily resulting from the high basis set dependence of the quadruples contribution.     

Revised vibrational band assignments for the experimental IR and Raman spectra of 2,3,4-trifluorobenzonitrile based on ab initio, DFT and normal coordinate calculations

In the present study, a systematic vibrational spectroscopic investigation for the experimental IR and Raman spectra of 2,3,4-trifluorobenzonitrile (TFB), aided by electronic structure calculations has been carried out. The electronic structure calculations – ab initio (RHF) and hybrid density functional methods (B3LYP) – have been performed with 6-31G* basis set. Molecular equilibrium geometries, electronic energies, IR intensities, harmonic vibrational frequencies, depolarization ratios and Raman activities have been computed. The results of the calculations have been used to simulate IR and Raman spectra for TFB that showed excellent agreement with the observed spectra. Potential energy distribution (PED) and normal mode analysis have also been performed. The assignments proposed based on the experimental IR and Raman spectra have been reviewed. A complete assignment of the observed spectra has been proposed.     

Analysis of the vibrational spectra of 1-methyluracil and its isotopic derivatives by ab initio calculations

The vibrational spectrum of 1-methyluracil trapped in an argon matrix has been analysed based on ab initio Hartree—Fock SCF calculations with a split-valence 4–21 basis set. The directly computed theoretical harmonic force field was scaled with empirical scale factors which were transferred from uracil (except for the methyl vibrational modes) to provide an a priori prediction of fundamental frequencies and intensities. The average deviations between experiment and prediction were 9.8 cm⁻¹ for the in-plane vibrations and 18.3 cm⁻¹ for the ring out-of-plane modes. After a few corrections of assignment of the observed spectrum, a new set of scale factors was optimized to give the best force field available from combined consideration of the experimental and theoretical information. These scale factors reduced the average deviations to 6.7 cm⁻¹ for the in-plane modes and to 11.7 cm⁻¹ for the out-of-plane ones. The vibrational spectra of seven isotopic derivatives (-2¹⁸0, -4¹⁸0, -3d, -5d, -6d, -5, 6d₂ and -CD₃) of 1-methyluracil were predicted with the force field resulting from the optimized set of scale factors, and compared with the crystal-phase experimental data. A few misassignments in the experimental isotopic spectra have been corrected. Vibrational absorption intensities have been computed and compared with experiment and with an earlier calculation.     

First-principles calculation of 17O and 25Mg NMR shieldings in MgO at finite temperature: rovibrational effect in solids

The temperature dependence of (17)O and (25)Mg NMR chemical shifts in solid MgO have been calculated using a first-principles approach. Density functional theory, pseudopotentials, a plane-wave basis set, and periodic boundary conditions were used both to describe the motion of the nuclei and to compute the NMR chemical shifts. The chemical shifts were obtained using the gauge-including projector augmented wave method. In a crystalline solid, the temperature dependence is due to both (i) the variation of the averaged equilibrium structure and (ii) the fluctuation of the atoms around this structure. In MgO, the equilibrium structure at each temperature is uniquely defined by the cubic lattice parameters, which we take from experiment. We evaluate the effect of the fluctuations within a quasiharmonic approximation. In particular, the dynamical matrix, defining the harmonic Hamiltonian, has been computed for each equilibrium volume. This harmonic Hamiltonian was used to generate nuclear configurations that obey quantum statistical mechanics. The chemical shifts were averaged over these nuclear configurations. The results reproduce the previously published experimental NMR data measured on MgO between room temperature and 1000 degrees C. It is shown that the chemical shift behavior with temperature cannot be explained by thermal expansion alone. Vibrational corrections due to the fluctuations of atoms around their equilibrium position are crucial to reproduce the experimental results.     

Vibrational spectra and structure of 5,6-diamino uracil and 5,6-dihydro-5-methyl uracil by density functional theory calculations

The molecular vibrations of 5,6-diamino uracil and 5,6-dihydro-5-methyl uracil were investigated in polycrystalline sample, at room temperature, by FT-IR and FT-Raman spectroscopies. The spectra were interpreted with the aid of normal coordinate analysis following a full structure optimization and force field calculations based on the density functional theory (DFT) using standard B3LYP/6-31G* and B3LYP/6-311+G** methods and basis set combinations. The DFT force field transformed to natural internal coordinates was corrected by a well-established set of scale factors that were found to be transferable to the title compounds. The infrared and Raman spectra were also predicted from the calculated intensities. Comparison of the simulated spectra with the experimental spectra provides important information about the ability of the computational method to describe the vibrational modes.     

First principles based group additive values for the gas phase standard entropy and heat capacity of hydrocarbons and hydrocarbon radicals

In this work a complete and consistent set of 95 Benson group additive values (GAVs) for standard entropies S(o) and heat capacities C(p)(o) of hydrocarbons and hydrocarbon radicals is presented. These GAVs include 46 groups, among which 25 radical groups, which, to the best of our knowledge, have not been reported before. The GAVs have been determined from a set of B3LYP/6-311G(d,p) ideal gas statistical thermodynamics values for 265 species, consistently with previously reported GAVs for standard enthalpies of formation. One-dimensional hindered rotor corrections for all internal rotations are included. The computational methodology has been compared to experimental entropies (298 K) for 39 species, with a mean absolute deviation (MAD) between experiment and calculation of 1.2 J mol(-1) K(-1), and to 46 experimental heat capacities (298 K) with a resulting MAD = 1.8 J mol(-1) K(-1). The constructed database allowed evaluation of corrections on S(o) and C(p)(o) for non-nearest-neighbor effects, which have not been determined previously. The group additive model predicts the S(o) and C(p)(o) within approximately 5 J mol(-1) K(-1) of the ab initio values for 11 of the 14 molecules of the test set, corresponding to an acceptable maximal deviation of a factor of 1.6 on the equilibrium coefficient. The obtained GAVs can be applied for the prediction of S(o) and C(p)(o) for a wide range of hydrocarbons and hydrocarbon radicals. The constructed database also allowed determination of a large set of hydrogen bond increments, which can be useful for the prediction of radical thermochemistry.     

Stacked but not Stuck: Unveiling the Role of π→π* Interactions with the Help of the Benzofuran–Formaldehyde Complex

The 1:1 benzofuran–formaldehyde complex has been chosen as model system for analyzing π→π* interactions in supramolecular organizations involving heteroaromatic rings and carbonyl groups. A joint “rotational spectroscopy–quantum chemistry” strategy unveiled the dominant role of π→π* interactions in tuning the intermolecular interactions of such adduct. The exploration of the intermolecular potential energy surface led to the identification of 14 low-energy minima, with 4 stacked isomers being more stable than those linked by hydrogen bond or lone-pair→π interactions. All energy minima are separated by loose transition states, thus suggesting an effective relaxation to the global minimum under the experimental conditions. This expectation has been confirmed by the experimental detection of only one species, which was unambiguously assigned owing to the computation of accurate spectroscopic parameters and the characterization of 11 isotopologues. The large number of isotopic species opened the way to the determination of the first semi-experimental equilibrium structure for a molecular complex of such a dimension.     

Bond dissociation energies for removal of the hydroxyl group in some alcohols from quantum chemical calculations

In this theoretical work, 22 alcohols and their geometric structure properties have been investigated employing quantum chemical methods to calculate the C? OH equilibrium bond distances and bond dissociation energies (BDEs). Since DFT methods have been researched to have low basis sets sensitivity for small and medium molecules in our previous work (Zhao et al., J Mol Struct, 2006, 766, 87), 22 title compounds have been studied by employing the hybrid density functional theory (B3LYP, B3PW91, B3P86, PBE1PBE) in conjunction with the 6‐311G** basis set and the complete basis set (CBS–Q) method. Comparison with the available experimental data shows that CBS–Q and B3P86 methods calculated results agree very well with the experimental values, with the average absolute errors of 1.3 kcal/mol and 3.5 kcal/mol, respectively. So considering the expensive computational time, CBS–Q method can be chosen as a satisfactory method of predicting the accurate BDEs for removal of the OH group in small and medium size alcohols. And B3P86 method may give accurate BDEs for larger alcohols we haven't studied. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Counterpoise corrected interaction energies are not systematically better than uncorrected ones: comparison with CCSD(T) CBS extrapolated values

The effect of the inclusion of counterpoise corrections (CP) on the accuracy of interaction energies has been studied for different systems accounting for (1) intermolecular interactions, (2) intramolecular interactions and (3) chemical reactions. To minimize the error associated with the method of choice, the energy calculations were performed using CCSD(T) in all the cases. The values obtained using aug-cc-pVXZ basis sets are compared to CBS-extrapolated values. It has been concluded that at least for the tested systems CP corrections systematically leads to results that differ from the CBS-extrapolated ones to a larger extension than the uncorrected ones. Accordingly, from a practical point of view, we do not recommend the inclusion of such corrections in the calculation of interaction energies, except for CBS extrapolations. The best way of dealing with basis set superposition error (BSSE) is not to use CP corrections, but to make a computational effort for increasing the basis set. This approach does not eliminate BSSE but significantly decreases it, and more importantly it proportionally decreases all the errors arising from the basis set truncation.     

Interpretation of the vibrational spectra of matrix-isolated uracil from scaled ab initio quantum mechanical force fields

The vibrational spectrum of uracil trapped in an argon matrix has been interpreted based on ab initio Hartree–Fock SCF calculations with a split-valence 4?21 basis set. The directly computed theoretical general valence force field was scaled with empirical scale factors in order to correct for the systematic errors originating in the limitation of the theoretical model. Scale factors transferred from related molecules provided a priori prediction of fundamental frequencies and intensities, permitting several corrections to be proposed for earlier assignments. Using the observed spectrum with the few altered assignments, a new set of scale factors was optimized to give the best force field available from combined consideration of the experimental and the theoretical data. For unknown reasons, the out-of-plane force field predicted a spectrum agreeing slightly less well with experiment than did the in-plane force field. However, the overall agreement between theory and experiment provided additional support for the assumptions involved in the method. The computed force fields were compared with others available from previous work. The comparison demonstrated the importance of expanding the energy surface around the true energy minimum and of using a proper scaling procedure. Previous scaled CNDO /2 calculations were found to be surprisingly good despite the large corrections required and the fact that they were made at an incorrect geometry.