Ab initio Study of Anharmonic Force Field and Spectroscopic Constants for Germanium Dichloride

Ab initio study of the equilibrium structure, spectroscopy constants, and anharmonic force field for several isotopomers of germanium dichloride (⁷⁰GeCl₂, ⁷²GeCl₂, and ⁷⁶GeCl₂) have been carried out at the MP2 and CCSD(T) levels of theory using cc-pVTZ basis set. The calculated geometries, rotational constants, vibration-rotation interaction constants, harmonic frequencies, anharmonic constants, quartic and sextic centrifugal distortion constants, cubic and quartic force constants are compared with experimental data. For small mass differences of the Ge isotopes, the isotopic effects for germanium dichloride are much weaker. The agreements are satisfactory for these two methods, but the deviations of CCSD(T) results are slightly larger than that of MP2, because of CCSD(T)'s inadequate treatment of electron correlation in hypervalent Cl atom.     

Calculations of vibrational energy levels by using a hybrid ab initio and DFT quartic force field: application to acetonitrile

A hybrid quartic force field with quadratic force constants calculated at the CCSD(T)/cc-pVTZ level and cubic and quartic terms determined by a B3LYP/cc-pVTZ treatment is proposed to compute the vibrational energy levels of acetonitrile from a variational method. Fundamentals and overtones calculated in the range of 300-3200 cm(-1) are in fair agreement with the 31 observed data, with an absolute mean deviation of less than 0.8%. These results allow us to explain several Fermi resonances.     

A coupled-cluster study of the structure and vibrational spectra of pyrazole and imidazole

Harmonic force fields were calculated at the corresponding optimized geometries for pyrazole and imidazole at the HF, B3LYP, MP2, CCSD and CCSD(T) levels using the 6-31G* basis set and at the HF and B3LYP levels using the cc-pVTZ basis set. The agreement between the calculated and experimental geometries by the CCSD and CCSD(T) methods was generally similar to that obtained with the B3LYP and MP2 methods. The force fields were scaled using one-scale-factor (1SF), 3SF and 7SF scaling schemes. The scale factors were varied with respect to the experimental frequencies. Using 7SF scaling, the root-mean-square (RMS) deviation of the calculated frequencies from the experimental frequencies by the HF, B3LYP, MP2, CCSD and CCSD(T) methods and the 6-31G* basis set was 16, 7, 13, 11 and 11 cm(-1), respectively. This shows that the B3LYP method is preferred for force field calculations over the perturbative MP2, CCSD and CCSD(T) methods. Using 1SF scaling, the CCSD(T) scale factor was 0.931, the highest among the five methods used but close to that obtained with the B3LYP method and the cc-pVTZ basis set with lower RMS deviation.     

Scaled quantum chemical force fields for 1,1-difluorocyclopropane and the influence of vibrational anharmonicity

Potential functions and harmonic (omega(i)) and anharmonic (nu(i)) fundamental frequencies have been calculated for 1,1-difluorocyclopropane (DFCP) and its d4 and d2 isotopomers using the program Gaussian 03. B3LYP and MP2 models were employed, each with the bases 6-311++G** and cc-pVTZ. Anharmonicity corrections Delta(i) = omega(i) - nu(i) are listed and shown to be different for symmetric and antisymmetric CH stretching modes in situations where Fermi resonance appears to be absent. The same effect is missing in C2H4, for which similar calculations were made. The quadratic force fields for DFCP have been scaled in symmetry coordinate space with 15 scale factors both to observed frequencies nu(obsd)and also to omega (obsd), where omega(obsd) = nu(obsd) + Delta. With nu(obsd) especially, different scale factors are needed for the symmetric and antisymmetric CH stretching force constants due to their differing anharmonicities. The source of the latter in the quartic and cubic force field is explored. MP2 calculations of valence interaction force constants involving the stretching of bonds on a common carbon atom are preferred to those from a B3LYP model. In either model, scaling to omega(obsd) rather than to nu(obsd) does not remove the necessity of varying scale factors for differing types of motion in the same group. Theoretical values of the five quartic centrifugal distortion constants are listed for the normal species and compared with new experimental data. The predictions are sufficiently good to be useful in fitting pure rotational transitions. A weakness is identified in the current Gaussian 03 code for the calculation of vibration-rotation quantities, and limitations are noted in the manner in which Fermi resonance is handled.     

The permanent dipole moment of gas-phase para-amino benzoic acid revisited

The permanent dipole moment of para-amino benzoic acid has been calculated at various theoretical levels, including Hartree-Fock, second-order M?ller-Plesset perturbation (MP2), coupled cluster singles and doubles (CCSD), and triples corrections CCSD(T), and hybrid density functional theory at B3LYP level. It is found that the B3LYP method fails to provide correct results for the geometry and the permanent dipole moment. These results are significantly improved by MP2 calculations. Our best estimated dipole moments obtained at CCSD and CCSD(T) levels are in good agreement with experiment.     

Semispectroscopic and quantitative structure-property relationship estimates of the equilibrium and vibrationally averaged structure and dipole moment of 1-buten-3-yne

Systematic quantum chemical calculations have been performed to obtain precise estimates of the equilibrium and vibrationally averaged molecular structure and electric dipole moment of vinylacetylene (VA, 1-buten-3-yne). Anharmonic (cubic and semi-diagonal quartic) MP2/cc-pVTZ force fields in normal coordinates were computed to account for anharmonic vibrational effects, including zero-point contributions to the rotational constants and the electric dipole moment. A simultaneous weighted least-squares structural refinement was performed, resulting in the best semispectroscopic estimate of the re structure of VA. The refinement was based on experimentally measured ground-state rotational constants of two isotopologs of VA corrected to equilibrium values using MP2/cc-pVTZ vibration-rotation interaction constants and all-electron CCSD(T)/aug-cc-pVTZ structural constraints. The semispectroscopic re structure of VA agrees excellently with the high-level CCSD(T)/aug-cc-pVTZ ab initio structure. The most dependable, CCSD(T)/cc-pVQZ//CCSD(T)/aug-cc-pVTZ equilibrium electric dipole moment of VA, in D, is mua= 0.4088, mub= 0.0004, and muc= 0. The vibrationally corrected a-component of 0.4214 D is in excellent agreement with one of the available experimental values. The present analysis shows that mub is negligible even after vibrational correction. A simple quantitative structure-property relationship (QSPR) model resulted in a highly similar estimate, 0.45 D, for the electric dipole moment of VA.     

Anharmonic force field of cis- and trans-formic acid from high-level ab initio calculations, and analysis of resonance polyads

The quadratic, cubic, and semidiagonal quartic force fields of cis- and trans-formic acid have been calculated using three different levels of theory. They all give satisfactory results, including the one at the lowest level of theory which is the MP2 method employing a basis set of triple-zeta quality. The results are used to theoretically analyze resonance polyads, including the one involving the 4(1), 5(1), 6(1), 8(1), 7(1)9(1), and 9(2) vibrational states. A semiexperimental equilibrium structure is derived from experimental ground state rotational constants and rovibrational interaction parameters calculated from the ab initio force field. The ab initio structure calculated at the CCSD(T) level of theory using a basis set of quintuple-zeta quality is in excellent agreement with the semiexperimental structure.     

Vibrational anharmonicity and harmonic force fields for dichloromethane from quantum-chemical calculations

Anharmonic and related constants have been calculated for CH2Cl2, CD2Cl2, and CHDCl2 by using the program Gaussian03 and B3LYP and MP2 models. Bases used were 6-311++G** and cc-pVTZ. The size of grid used in the B3LYP/6-311++G** model had a noticeable effect on resulting data. Features of the MP2/6-311++G** calculations suggested a deleterious effect of the absence of f functions in this basis set. The need for the replacement of second-order terms in the perturbation theory formulas for the vibrational anharmonic constants x ij in the presence of Fermi resonance was explored, and minor resonances were found associated with the cubic constants varphi 122 and varphi 299 (d 0 isotopomer), phi122 and phi849 (d2), and phi278 (d1). Computed xij values for nuCH and nuCD motions agree quite well with earlier experimental data. Observed anharmonic frequencies, nu obsd, were corrected to "observed" harmonic frequencies, omega obsd, by using computed differences Delta = omegaQC-nuQC. These differences Delta are larger for the antisymmetric nuasCH2 mode than for symmetric nusCH2 motion. This fact made it necessary to use differing scale factors for the two kinds of CH stretching force constants in a subsequent scaling of the harmonic force field to nuobsd. Force field scaling was also carried out by refining to omega obsd. In both approaches, the B3LYP models required differing scale factors for symmetric and antisymmetric CCl stretching force constants, indicating a failure to compute an accurate C-Cl stretch-stretch interaction force constant. The MP2/cc-pVTZ force field was preferred. Both scaled and unscaled harmonic force fields were used to calculate centrifugal distortion constants (CDCs) and contributions to the vibrational dependence of the rotational constants (alphas). Variations in the CDCs can, in part, be explained by the magnitudes of the frequencies used in the scaling process.     

Equilibrium structure of sulfuric acid

The equilibrium structure of the more stable conformer of H2SO4, of C2 symmetry, has been calculated ab initio using the CCSD(T) method and taking into account the core correlation correction. The accuracy of this structure has been checked by comparing it to that of similar molecules and by estimating the effects of basis set enlargement and of diffuse functions. Furthermore, the quadratic, cubic, and quartic force fields have been calculated at the MP2 level of theory using a basis set of triple-zeta quality. The spectroscopic constants derived from the force field are in satisfactory agreement with the experimental ones. The resulting band origins are compared to literature infrared values, including those for overtone and combination bands. Normal modes of vibration are pictured. Using this force field, semiexperimental equilibrium rotational constants are determined which allows us to check the accuracy of the ab initio structure and to refine it using a mixed regression method.     

Anharmonic Effect of the Unimolecular Dissociation of Ethyl Radical

The anharmonic and harmonic rate constants have been calculated for the unimolecular dissociation of ethyl radical using the method proposed by Yao and Lin (YL method) at both B3LYP/6‐311++G** and MP2/6‐311++G** levels. The different rate constants indicate that the results obtained from B3LYP and MP2 method are very close. The anharmonic and tunneling effect of the title reaction has also been examined. The comparison shows that, both in microcanonical and canonical systems, the anharmonic rate constants are higher than those for harmonic cases, especially in the case of high total energies and temperatures, which indicates that anharmonic effect of the unimolecular dissociation of ethyl into C₂H₄ and H is so significant that cannot be neglected. The tunneling effect is very small for the decomposition of C₂H₅ radical.     

Anharmonic Effect in CH<Subscript>3</Subscript>CH<Subscript>2</Subscript>C(=O)OCH<Subscript>2</Subscript>CH<Subscript>3</Subscript> Decomposition

In this paper, using the B3LYP functional and CCSD(T) method with 6-311++G** basis set, the harmonic and anharmonic rate constants in the unimolecular dissociation of ethyl propanoate have been calculated using Rice–Ramsperger–Kassel–Marcus theory. The anharmonic rate constants of the title reaction have also been examined, the comparison shows that, the anharmonic effect especially in the case of high total energies and temperature for channels 3 to 6 is significant, so that the anharmonic effect cannot be neglected for unimolecular dissociation reaction of CH₃CH₂C(=O)OCH₂CH₃ both in microcanonical and canonical systems.     

Theoretic studies on the kinetics and mechanism of multi-channel gas-phase unimolecular reactions of 1-chloropropane and 2-chloropropane

The potential energy surfaces of the gas-phase unimolecular decomposition reactions of 1-chloropropane and 2-chloropropane are investigated by various quantum chemical methods including CCSD(T), CCSD, GBS-QB3, C3B3, MP4, MP2, and B3LYP. Modified strong collision/RRKM theory was used to calculate the unimolecular rate constants as a function of pressure and temperature. It is found that the major reaction pathway is the HCl elimination. It is predicted that bond dissociation reactions have negligible contribution to the overall rate constant. The computed rate constants are compared with the available experimental data.     

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.     

甲基肼自由基在高温燃烧分解反应下的非谐振效应

In this work, the harmonic and anharmonic rate constants of the decomposition reaction of monomethylhydrazine (MMH) radicals have been calculated by using transition state (TS) and Rice-Ramsperger-Kassel-Marcus (RRKM) theories with either MP2 or B3LYP method at 6-311++G (3df, 2p) basis set, respectively. The reaction mechanism and anharmonic effect of the MMH radicals are studied in detail and both of the harmonic and anharmonic rate constants increase sharply with increasing temperature in the canonical system. In the microcanonical system, these constants also show sharp increase with the energies. Overall, the anharmonic effect becomes more pronounced with the increasing temperature or energy in the canonical and microcanonical systems, respectively. These results indicate that the anharmonic effect of the decomposition reaction of MMH radicals is quite significant and cannot be ignored.     

Predicting Rate Constants for Nucleophilic Reactions of Amines with Diarylcarbenium Ions Using an ONIOM Method

The rate constants of the nucleophilic reactions between amines and benzhydrylium ions were calculated using first-principles theoretical methods. Solvation models including PCM, CPCM, and COSMORS, as well as different types of atomic radii including UA0, UAKS, UAHF, Bondi, and UFF, and several single-point energy calculation methods (B3LYP, B3P86, B3PW91, BHANDH, PBEPBE, BMK, M06, MP2, and ONIOM method) were exam-ined. By comparing the correlation between experimental rate constants and the calculated values, the ONIOM(CCSD(T)/6-311++G(2df,2p):B3LYP/6-311++G(2df,2p))// B3LYP/6-31G(d)/PCM/UFF) method was found to perform the best. This method was then employed to calculate the rate constants of the reactions between diverse amines and diarylcarbenium ions. The calculated rate constants for 65 reactions of amines with diarylcarbenium ions are in agreement with the experimental values, indicating that it is feasible to predict the rate constant of a reaction between an amine and a diarylcarbenium ion through ab initio calculation.     

Density functional theory and ab initio studies of the structure and energetics of digallium tetraoxide, Ga2O4, in the gas phase

A systematic investigation on the neutral and anionic digallium tetraoxide, Ga(2)O(4) has been carried out by using density functional theory (DFT), second-order M?ller-Plesset perturbation theory (MP2), and the coupled cluster approach with single and double substitutions and a perturbative treatment of the triple excitations [CCSD(T)]. The geometry of neutral Ga(2)O(4) has been proposed earlier, from an experimental study, to adopt a D(2d) symmetry (J. Phys. Chem. 1979, 83, 656). However, the current research reveals that, out of the several isomers considered for neutral and anionic digallium tetraoxide, the (3)B(1u) and (2)B(3g) of the planar D(2h) geometry (7a-D(2h)) are the lowest-energy states for Ga(2)O(4) and Ga(2)O(4)(-). Our computations rule out the D(2d) geometry (3-D(2d)) as a viable contender for neutral Ga(2)O(4). The (3)B(2) (3-D(2d)) state is located above the (3)B(1u) (7-D(2h)) state by at least 4.26 eV. The energies of the low lying states, geometrical parameters, and energetic features (VEDE, AEDE, and AEA) are reported. The AEA of Ga(2)O(4) is calculated to be 3.94 eV (B3LYP), 3.24 eV (MP2), 3.42 eV [CCSD(T)//B3LYP], and 3.38 eV [CCSD(T)//MP2], respectively. In addition, the dissociation energy, D(e), for the process Ga(2)O(4) ((3)B(1u)) → 2GaO(2) ((2)A(2)) is 3.59 eV (B3LYP), 5.08 eV (MP2), 4.82 eV [CCSD(T)//B3LYP], and 4.80 eV [CCSD(T)//MP2]. The results obtained in this work are critically analyzed, discussed, and compared with those of the analogous metal oxides.     

Quantum chemical study of the hydrogen-bonded HXeOH–H2O complex

A combined MP2 and DFT/B3LYP study of the HXeOH–H₂O complex is presented. These computational methods have been used to extract information on the structural, energetical and vibrational properties of the complex. Additionally, we have applied anharmonic vibrational calculations based on the MP2-computed intermolecular potential energy surface. Large perturbations both on the subunit structures and their fundamental vibrational modes are found upon complexation. Large changes of anharmonicity of the HXeOH subunit reflects the perturbation of the molecule's electronic structure. The computed BSSE-corrected interaction energies are −40.23 and −38.94 kJ mol⁻¹ at the CCSD(T)//MP2 and CCSD(T)//B3LYP levels of theory, respectively. The estimated deformation energy contribution to the interaction energy is about 5%, which is very large compared with classical hydrogen-bonded complexes. The topological analysis of the Electron Localization Function (ELF) was applied to study further the hydrogen-bonded interaction between the two complex partners. The obtained interaction pattern suggests that the interaction between HXeOH and H₂O is a typical hydrogen bond interaction driven mainly by electrostatic interactions.     

Anharmonic analysis of the vibrational spectra of some cyanides and related molecules of astrophysical importance

A detailed analysis of the vibrational spectra of carbonyl cyanide, diethynyl ketone and acetyl cyanide has been conducted in harmonic and anharmonic approximations. RHF, MP2 and density functional theory (DFT) methods with 6-311++G(2df,2p) basis sets and B3LYP functionals have been employed. Spectroscopic constants such as anharmonicity constants, rotational and centrifugal distortion constants, rotation-vibration coupling constants and Coriolis coupling coefficients have been calculated for each molecule and compared with the experimental data, where available. A close agreement between the calculated and experimental values of the spectroscopic constants has been obtained. Complete assignments have been provided to the fundamental bands, overtones and combination tones of the molecules. Density functional theory based anharmonic frequencies compare well with the experimental frequencies within +/-18 cm(-1) on an average. RHF and MP2 methods, however, give much higher values for the frequencies that need scaling even in the anharmonic approximation.     

Microwave spectra and barriers to internal rotation of z- and e-1-propenyl isocyanide

A synthetic procedure yielding a mixture of Z- and E-1-propenyl isocyanide (CH(3)CH═CHNC) is described. The microwave spectrum of this mixture has been recorded in the 12-100 GHz spectral range, and the spectra of the Z and E isomers have been assigned for the first time. Most transitions of the Z form were split into two components of equal intensity due to tunneling of the methyl group, which allowed the barrier to internal rotation of this group to be determined as 4.0124(12) kJ/mol by fitting 568 transitions with a maximum value of J = 46 using the computer program Xiam. This fit had a root-mean-square deviation as large as 4.325. The same transitions were therefore fitted anew using the more sophisticated program Erham. This fit has a rms deviation marginally better (4.136) than the Xiam fit. No split MW lines were found for E-1-propenyl isocyanide. The absence of splittings is ascribed to a barrier to internal rotation of the methyl group that is significantly higher than the barrier of the Z isomer. It is concluded that the barrier must be larger than 6 kJ/mol for the E form. The experimental work was augmented by quantum chemical calculations at CCSD/cc-pVTZ, B3LYP/cc-pVTZ, and MP2/cc-pVTZ levels of theory. The CCSD method predicts rotational constants of the Z and E forms well. The B3LYP barriers to internal rotation of a series of substituted propenes were calculated and found to be in good agreement with experiments. Calculations of the quartic centrifugal distortion constants of the two 1-propenyl isocyanides by the B3LYP and MP2 methods were less successful.