A theoretical and experimental study of Sb4O6: vibrational analysis, infrared, and Raman spectra

The first ab initio theoretical study of tetraantimony hexoxide (Sb4O6) is reported. The normal mode frequencies, intensities, and the corresponding vibrational assignments of Sb4O6 in T(d) symmetry were calculated using the GAUSSIAN 98 set of quantum chemistry codes at the Hartree-Fock (HF)/CEP-121G, M?ller-Plesset (MP2)/CEP-121G, and density functional theory (DFT)/B3LYP/CEP-121G levels of theory. By comparison to experimental data deduced by our laboratory and others, correction factors for the calculated vibrational frequencies were determined and compared. Normal modes were decomposed into three non-redundant motions (Sb-O-Sb stretch, Sb-O-Sb bend, and Sb-O-Sb wag). Percent relative errors found for the HF, DFT, and MP2 corrected frequencies when compared to experiment are 5.8, 6.1, and 5.7 cm(-1), respectively. Electron distributions for selected molecular orbitals are also considered.     

Vibrational frequencies and structural determination of trimethylarsine oxide

The normal mode frequencies and corresponding vibrational assignments of trimethylarsine oxide are examined theoretically using the Gaussian 98 set of quantum chemistry codes. All normal modes were successfully assigned to one of eight types of motion (As-C stretch, As=O stretch, C-H stretch, C-As-C bend, As=O bend, H-C-H bend, CH3 wag, and CH3 twist) utilizing the C3v symmetry of the molecule. Calculations were performed at the Hartree-Fock, DFT(B3LYP), and MP2 levels of theory using the standard 6-311G** basis. Calculated infrared intensities and Raman activities are reported.     

Why does argon bind to deuterium? Isotope effects and structures of Ar x H5O2(+) complexes

Recently, we reported the spectrum of Ar x D4HO2(+) [McCunn; et, al. J. Phys. Chem. B 2008, 112, 321], and here, we extend that work to include the Ar x H4DO2(+) isotopologue in order to explore why the Ar atom has a much greater propensity for attachment to a dangling OD group than it does for OH, even when many more of the latter binding sites are available. Calculated (MP2/6-311+G(d,p) level of theory/basis) harmonic frequencies reproduce the observed multiplet patterns of OH and OD stretches and confirm the presence of various isomers arising from the different Ar binding sites. The preferential bonding of Ar to OD is traced to changes in the frequencies of the wag and rock modes of the H5O2(+) moiety rather than to shifts in the oscillator that directly binds the Ar atom.     

C24团簇结构与稳定性的理论研究

采用量子化学HF, B3LYP和MP2方法,选用6-31G^*, 6-311G^*, cc-pVDZ和cc-pVTZ基组,对C₂₄团簇的6种异构体进行了优化,并对它们的几何构型、振动频率、核独立化学位移(NICS)和稳定性进行了讨论,比较C₂₄团簇各种异构体的稳定性.研究表明:在6-311G^*和cc-pVDZ水平上,B3LYP方法给出的稳定性大小顺序分别为c＞f＞b＞e＞a＞d和c＞b＞f＞a＞e＞d, MP2方法给出的稳定性大小顺序为b＞c＞a＞e＞f＞d.     

Intermolecular potentials of the methane dimer calculated with Møller-Plesset perturbation theory and density functional theory

We have calculated the intermolecular interaction potentials of the methane dimer at the minimum-energy D(3d) conformation using the Hartree-Fock (HF) self-consistent theory, the correlation-corrected second-order M?ller-Plesset (MP2) perturbation theory, and the density functional theory (DFT) with the Perdew-Wang (PW91) functional as the exchange or the correlation part. The HF calculations yield unbound potentials largely due to the exchange-repulsion interaction. In the MP2 calculations, the basis set effects on the repulsion exponent, the equilibrium bond length, the binding energy, and the asymptotic behavior of the calculated intermolecular potentials have been thoroughly studied. We have employed basis sets from the Slater-type orbitals fitted with Gaussian functions (STO-nG) (n=3-6) [Quantum Theory of Molecular and Solids: The Self-Consistent Field for Molecular and Solids (McGraw-Hill, New York, 1974), Vol. 4], Pople's medium size basis sets of Krishnan et al. [J. Chem. Phys. 72, 650 (1980)] [up to 6-311++G(3df,3pd)] to Dunning's correlation consistent basis sets [J. Chem. Phys. 90, 1007 (1989)] (cc-pVXZ and aug-cc-pVXZ) (X=D, T, and Q). With increasing basis size, the repulsion exponent and the equilibrium bond length converge at the 6-31G** basis set and the 6-311++G(2d,2p) basis set, respectively, while a large basis set (aug-cc-pVTZ) is required to converge the binding energy at a chemical accuracy (approximately 0.01 kcal/mol). Up to the largest basis set used, the asymptotic dispersion coefficient has not converged to the destined C6 value from molecular polarizability calculations. The slow convergence could indicate the inefficacy of using the MP2 calculations with Gaussian-type functions to model the asymptotic behavior. Both the basis set superposition error (BSSE) corrected and uncorrected results are presented to emphasize the importance of including such corrections. Only the BSSE corrected results systematically converge to the destined potential curve with increasing basis size. The DFT calculations generate a wide range of interaction patterns, from purely unbound to strongly bound, underestimating or overestimating the binding energy. The binding energy calculated using the PW91PW91 functional and the equilibrium bond length calculated using the PW91VP86 functional are close to the MP2 results at the basis set limit.     

Neutron spectroscopic study of hydrogen bonding dynamics in L-serine

Inelastic incoherent neutron scattering (IINS) spectra were obtained at 10 K for normal and deuterated L-serine. The geometry of L-serine molecule was optimized for the zwitterion form using ab initio HF, MP2 and DFT (B3LYP) levels with 6-31G* and 6-311 + +G4** basis sets. The theoretical frequencies of normal and d4-L-serine were compared with IINS spectra. Normal coordinate analysis and band assignments based on ab initio calculations and experimental data were presented. IINS frequencies due to the out-of-plane gamma(N-H...O) hydrogen bond motions were observed and identified.     

An evaluation of harmonic vibrational frequency scale factors

Scale factors for obtaining fundamental vibrational frequencies, low-frequency vibrational frequencies, zero-point vibrational energies (ZPVEs), and thermal contributions to enthalpy and entropy have been derived through a least-squares approach from harmonic frequencies determined at more than 100 levels of theory. Wave function procedures (HF, MP2, QCISD, QCISD(T), CCSD, and CCSD(T)) and a large and representative range of density functional theory (DFT) approaches (B3-LYP, BMK, EDF2, M05-2X, MPWB1K, O3-LYP, PBE, TPSS, etc.) have been examined in conjunction with basis sets such as 6-31G(d), 6-31+G(d,p), 6-31G(2df,p), 6-311+G(d,p), and 6-311+G(2df,p). The vibrational frequency scale factors were determined by a comparison of theoretical harmonic frequencies with the corresponding experimental fundamentals utilizing a standard set of 1066 individual vibrations. ZPVE scale factors were generally obtained from a comparison of the computed ZPVEs with experimental ZPVEs for a smaller standard set of 39 molecules, though the effect of expansion to a 48 molecule data set was also examined. In addition to evaluating the scale factors for a wide range of levels of theory, we have also probed the effect on scale factors of varying the percentage of incorporated exact exchange in hybrid DFT calculations using a modified B3-LYP functional. This has revealed a near-linear relationship between the magnitude of the scale factor and the proportion of exact exchange. Finally, we have investigated the effect of basis set size on HF, MP2, B3-LYP, and BMK scale factors by deriving values with basis sets ranging from 6-31G(d) up to 6-311++G(3df,3pd) as well as with basis sets in the cc-pVnZ and aug-cc-pVnZ series and with the TZV2P basis.     

Tricyclo[2.2.0.0]hexane: a new hypothetical molecule which should have only one inverted carbon atom

HF and MP2 calculations with the 6-31G** and 6-311G** basis sets and those at MP2/cc-pVTZ level were carried out for the hypothetical tricyclo[2.2.0.0^1,3]hexane. The results indicate that the molecule under study should have one carbon atom with highly unusual inverted configuration. The analysis of the vibrational frequencies of this molecule as well as the analysis of its plausible decomposition routes performed at the DFT level indicate that this unique molecule could be a plausible synthetic target.     

Comparative theoretical study of transition structures,barrier heights,and reaction energies for the intramolecular tautomerization in acetaldehyde/vinyl alcohol and acetaldimine/vinylamine systems

The transition structures associated with the possible intramolecular tautomerization for acetaldehyde/vinyl alcohol and acetaldimine/vinylamine systems as models of keto/enol and imine/enamine interconversion processes, respectively, were characterized. The relative stabilities of the tautomers and the associated barrier heights were calculated. Ab initio analytical gradients and second derivatives at the HF level of theory and 3-21G, 6-31G, 6-31G**, 6-31++G**, and 6-311++G** basis-set, DFT (BP86/6-311++G** and BLYP/6-311++G**), and semiempirical (AM1 and PM3) procedures were used to identify the stationary points. Correlation effects were estimated using the perturbational approach at MP2/6-31G**, MP2/6-311++G**, and MP2/6-311++G (3df,2p) levels. The geometry, electronic structure, harmonic vibrational frequencies, and transition vector associated with the transition structures as well as the relative stabilities of different isomers and barrier heights were analyzed. The dependence of these properties upon theoretical methods is analyzed and discussed. The transition structures are four-membered rings and the corresponding transition vectors are associated to collective fluctuations. The 1,3 intramolecular hydrogen migration is much more advanced than are the hybridization changes on donor and acceptor centers at the transition structure. The corresponding barrier heights can be related to the change of bond orders and acid/base properties of these centers. A comparison of the results obtained with different methods renders that the nature of the transition structure seems to be a rather robust entity. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 9–24, 1998     

Anharmonic vibrational frequency shifts upon interaction of phenol+ with the open shell ligand O2. The performance of DFT methods versus MP2

Anharmonic vibrational frequency shifts of the phenol(+) O-H stretching mode upon complex formation with the open-shell ligand O(2) were computed at several DFT and MP2 levels of theory, with various basis sets, up to 6-311++G(2df,2pd). It was found that all DFT levels of theory significantly outperform the MP2 method with this respect. The best agreement with the experimental frequency shift for the hydrogen-bonded minimum on the potential energy surfaces was obtained with the HCTH/407 functional (-93.7 cm(-1) theoretical vs -86 cm(-1) experimental), which is a significant improvement over other, more standard DFT functionals (such as, e.g., B3LYP, PBE1PBE), which predict too large downshifts (-139.9 and -147.7 cm(-1), respectively). Good agreement with the experiment was also obtained with the mPW1B95 functional proposed by Truhlar et al. (-109.2 cm(-1)). We have attributed this trend due to the corrected long-range behavior of the HCTH/407 and mPW1B95 functionals, despite the fact that they have been designed primarily for other purposes. MP2 method, even with the largest basis set used, manages to reproduce only less than 50% of the experimentally detected frequency downshift for the hydrogen-bonded dimer. This was attributed to the much more significant spin contamination of the reference HF wave function (compared to DFT Kohn-Sham wave functions), which was found to be strongly dependent on the O-H stretching vibrational coordinate. All DFT levels of theory outperform MP2 in the case of computed anharmonic OH stretching frequency shifts upon ionization of the neutral phenol molecule as well. Besides the hydrogen-bonded minimum, DFT levels of theory also predict existence of two other minima, corresponding to stacked arrangement of the phenol(+) and O(2) subunits. mPW1B95 and PBE1PBE functionals predict a very slight blue shift of the phenol(+) O-H stretching mode in the case of stacked dimer with the nearly perpendicular orientation of oxygen molecule with respect to the phenolic ring, which is entirely of electrostatic origin, in agreement with the experimental observations of an additional band in the IR photodissociation spectra of phenol(+)-O(2) dimer [Patzer, A.; Knorke, H.; Langer, J.; Dopfer, O. Chem. Phys. Lett. 2008, 457, 298]. The structural features of the minima on the studied PESs were discussed in details as well, on the basis of NBO and AIM analyses.     

Density functional theory and post-Hartree-Fock studies on molecular structure and harmonic vibrational spectrum of formaldehyde

Density functional theory (DFT) with the Becke's three-parameter exchange correlation functional and the functional of Lee, Yang and Parr, gradient-corrected functionals of Perdew, and Perdew and Wang [the DFT(B3LYP), DFT(B3P86) and DFT(B3PW91) methods, respectively], and several levels of conventional ab initio post-Hartree-Fock theory (second- and fourth-order perturbation theory M?ller-Plesset MP2 and MP4(SDTQ), coupled cluster with the single and double excitations (CCSD), and CCSD with perturbative triple excitation [CCSD(T)], configuration interaction with the single and double excitations [CISD], and quadratic configuration interaction method [QCISD(T)], using several basis sets [ranging from a simple 6-31G(d,p) basis set to a 6-311+ +G(3df, 2pd) one], were applied to study of the molecular structure (geometrical parameters, rotational constants, dipole moment) and harmonized infrared (IR) spectrum of formaldehyde (CH₂O). High-level ab initio methods CCSD(T) and QCISD(T) with the 6-311+ +G(3df, 2pd) predict correctly molecular parameters, vibrational harmonic wavenumbers and the shifts of the harmonic IR spectrum of ¹²CH₂ ¹⁶O upon isotopic substitution. Received: 30 January 1997 / Accepted: 7 May 1997     

Ab initio and DFT studies of the molecular structures and vibrational spectra of succinonitrile

The Molecular structure, conformational stability and vibrational frequencies of succinonitrile NCCH2CH2CN have been investigated with ab initio and density functional theory (DFT) methods implementing the standard 6-311++G* basis set. The potential energy surfaces (PES) have been explored at DFT-B3LYP, HF and MP2 levels of theory. In agreements with previous experimental results, the molecule was predicted to exist in equilibrium mixture of trans and gauche conforms with the trans form being slightly lower in energy. The vibrational frequencies and the corresponding vibrational assignments of succinonitrile in both C2h and C2 symmetry were examined theoretically and the calculated Infrared and Raman spectra of the molecule were plotted. Observed frequencies for normal modes were compared with those calculated from normal mode coordinate analysis carried out on the basis of ab initio and DFT force fields using the standard 6-311++G* basis set of the theoretical optimized geometry. Theoretical IR intensities and Raman activities are reported.     

DFT and experimental studies of the structure and vibrational spectra of curcumin

The potential energy surface of curcumin [1,7‐bis(4‐hydroxy‐3‐methoxyphenyl)‐1,6‐heptadiene‐3,5‐dione] was explored with the DFT correlation functional B3LYP method using 6‐311G* basis. The single‐point calculations were performed at levels up to B3LYP/6‐311++G**//B3LYP/6‐311G*. All isomers were located and relative energies determined. According to the calculation the planar enol form is more stable than the nonplanar diketo form. The results of the optimized molecular structure are presented and compared with the experimental X‐ray diffraction. In addition, harmonic vibrational frequencies of the molecule were evaluated theoretically using B3LYP density functional methods. The computed vibrational frequencies were used to determine the types of molecular motions associated with each of the experimental bands observed. Our vibrational data show that in both the solid state and in all studied solutions curcumin exists in the enol form. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Comparison of ab initio and density functional methods for vibrational analysis of TeCl4

Vibrational analysis of tellurium tetrachloride, TeCl₄, was performed with Hartree–Fock (HF), MP2, and generalized gradient approximation density functional theory (DFT) methods supplemented with polarized double-zeta split valence (DZVP) basis sets and relativistic effective core potentials (RECP) of Hay and Wadt. The molecular geometry is best reproduced at the HF and MP2/RECP+DZVP [polarized Hay and Wadt RECP for Te and 6–31G(d) basis set for Cl] levels of theory. The DFT methods gave rise to poorer results, especially those using Becke's 1988 exchange functional. Generally, the vibrational frequencies calculated by the MP2 and B3-type DFT methods with the all electron and RECP+DZVP basis sets as well as at the HF/RECP level were in satisfactory accord with the experimental data. The agreement was good enough to assist the assignment of the measured vibrational spectra. The best agreement with the experimental vibrational frequencies was achieved with the scaled HF/RECP force field. Consistent results were obtained for the unobserved A₂ (ν₄) fundamental, where the results of the best methods were within 4 cm⁻¹. The best force fields were obtained with the following methods: Becke3–Lee–Yang–Parr and Becke3–Perdew/all electron basis, MP2 and Becke3-Perdew/RECP+DZVP, and HF/RECP. The methods using RECPs are advantageous for large-scale computations. The RECP basis set effectively compensates the errors of the HF method for TeCl₄; however, it provides poor results with correlated methods. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 308–318, 1998     

G3 and density functional theory investigations of the structures and energies of SF(n)Cl (n=0-5) and their anions

Computations of structures and total energies have been carried out for neutral and anionic SF(n)Cl (n=0-5), using the composite G3 method and density functional theory (DFT) at the B3LYP6-311+G(3df) level. The total energies and zero-point energies have been used here to derive electron affinities, bond dissociation energies, and heats of formation. In addition, vibrational frequencies, polarizabilities, and dipole moments are reported. Results are compared with earlier work for SF(m) (m=1-6) and demonstrate how the relatively weak S-Cl bond and reduced symmetry influence the properties of these molecules and anions. Comparisons are also made between G3 and DFT results for SF(n)Cl. Of particular interest is the alternating pattern of agreement between calculated electron affinity values with n. These calculations also provide critical energetic data needed to understand experimental measurements of electron attachment to SF(5)Cl [Van Doren et al., J. Chem. Phys. 128, 094309 (2008)] for which numerous ion products have been reported in the literature at low electron energy.     

Molecular structure and vibrational analysis of 3-Ethylpyridine using ab initio HF and density functional theory (B3LYP) calculations

The FT-Raman and FT-IR spectra for 3-Ethylpyridine (3-EP) have been recorded in the region 4000-100 cm(-1) and compared with the harmonic vibrational frequencies calculated using HF/DFT (B3LYP) method by employing 6-31G(d,p) and 6-311++G(d,p) basis set with appropriate scale factors. IR intensities and Raman activities are also calculated by HF and DFT (B3LYP) methods. Optimized geometries of the molecule have been interpreted and compared with the reported experimental values of some substituted benzene. The experimental geometrical parameters show satisfactory agreement with the theoretical prediction from HF and DFT. The scaled vibrational frequencies at B3LYP/6-311++G(d,p) seem to coincide with the experimentally observed values with acceptable deviations. The theoretical spectrograms (IR and Raman) have been constructed and compared with the experimental FT-IR and FT-Raman spectra. Some of the vibrational frequencies of the pyridine are effected upon profusely with the C2H5 substitutions in comparison to pyridine and these differences are interpreted.     

Structure and vibrational frequencies of 2-butanimine

The conformational behavior and structural stability of 2-butanimine were investigated by utilizing ab initio calculations with 6-311++G** basis set at HF, MP2, B3LYP and BLYP levels. The vibrational frequencies of 2-butanimine were computed. Complete vibrational assignments were made on the basis of normal coordinate calculations for stable conformer of the molecule. HF results without scaled quantum mechanical (SQM) force field procedure considered are in bad agreement with experimental values. Of the two DFT methods, BLYP reproduces the observed fundamental frequencies most satisfactorily with the mean absolute deviation of the non-CH stretching modes less than 21.3 cm(-1). The results indicate that BLYP calculation is a very promising approach for understanding the observed spectral features.     

去氢抗坏血酸分子振动光谱的理论研究

采用RHF, MP2, DFT(B3LYP)方法, 以6-311++G**为基组研究了去氢抗坏血酸分子(DHA)的平衡几何构型和振动光谱. 计算结果表明, 采用RHF, B3LYP以及MP2 方法优化得到的几何结构以及频率值是一致的. 采用B3LYP/6-311++G**计算了DHA分子平衡构型下的谐振动力场﹑振动频率和振动强度. 使用Wilson的GF矩阵方法对DHA分子进行了简正坐标分析, 依据所得的势能分布对DHA分子的振动基频进行了合理的理论归属.     

Vibrational spectroscopic (FTIR and FTRaman) investigation using ab initio (HF) and DFT (LSDA and B3LYP) analysis on the structure of Toluic acid

The FTRaman and FTIR spectra for Toluic acid (TA) have been recorded in the region 4000-100 cm(-1) and compared with the harmonic vibrational frequencies calculated using HF/DFT (LSDA and B3LYP) method BY employing 6-311G (d, p) basis set with appropriate scale factors. IR intensities and Raman activities are also calculated by HF and DFT (LSDA/B3LYP) methods. Optimized geometries of the molecule have been interpreted and compared with the reported experimental values for benzoic acid and some substituted benzoic acids. The experimental geometrical parameters show satisfactory agreement with the theoretical prediction from HF and DFT. The scaled vibrational frequencies at B3LYP/6-311G (d, p) seem to coincide with the experimentally observed values with acceptable deviations. The theoretical spectrograms (IR and Raman) have been constructed and compared with the experimental FT-IR and FT-Raman spectra. Some of the vibrational frequencies of the TA are effected upon profusely with the methyl substitutions in comparison to benzoic acid and these differences are interpreted.