Evaluation of the enthalpy of formation, proton affinity, and gas-phase basicity of gamma-butyrolactone and 2-pyrrolidinone by isodesmic reactions

The knowledge of thermochemical parameters such as the enthalpy of formation, gas-phase basicity, and proton affinity may be the key to understanding molecular reactivity. The obtention of these thermochemical parameters by theoretical chemical models may be advantageous when experimental measurements are difficult to accomplish. The development of ab initio composite models represents a major advance in the obtention of these thermochemical parameters, but these methods do not always lead to accurate values. Aiming at achieving a comparison between the ab initio models and the hybrid models based on the density functional theory (DFT), we have studied gamma-butyrolactone and 2-pyrrolidinone with a goal of obtaining high-quality thermochemical parameters using the composite chemical models G2, G2MP2, MP2, G3, CBS-Q, CBS-4, and CBS-QB3; the DFT methods B3LYP, B3P86, PW91PW91, mPW1PW, and B98; and the basis sets 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d), 6-311+G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ. Values obtained for the enthalpies of formation, proton affinity, and gas-phase basicity of the two target molecules were compared to the experimental data reported in the literature. The best results were achieved with the use of DFT models, and the B3LYP method led to the most accurate data.     

Theoretical study for the HNSi–HSiN isomerization: ab initio and DFT calculations

Ab initio and density functional theory (DFT) calculations using the GAUSSIAN 94 program have been performed to investigate the molecular structures of HNSi and HSiN in the ground state as well as the transition state for the HNSi–HSiN isomerization reaction at the 6-311G(d,p), 6-311+G(2d,p) and 6-311+G(2df,p) basis sets. The results show that DFT calculations at higher levels of theory reproduce experimental vibrational frequencies of both HNSi and HSiN better than ab initio methods including electron correlation effects. Those calculated geometries are accurate enough to predict the rotational constant of HNSi. The barrier height for the isomerization reaction is found to be about 10 kcal/mol.     

Structure and conformational stability of CH2CHCH2X (X=F, Cl and Br) molecules—post Hartree–Fock and density functional theory methods

The molecular structure and conformational stability of CH₂CHCH₂X (X=F, Cl and Br) molecules were studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 3-fluoropropene were optimized employing BLYP and B3LYP levels of theory of DFT method implementing 6-311+G(d,p) basis set. The MP2/6-31G^*, BLYP and B3LYP levels of theory of ab initio and DFT methods were used to optimize the 3-chloropropene and 3-bromopropene molecules. The structural and physical parameters of the molecules are discussed with the available experimental values. The rotational potential energy surface of the above molecules were obtained at MP2/6-31G^* and B3LYP/6-311+G(d,p) levels of theory. The Fourier decomposition of the rotational potentials were analyzed. The HF/6-31G^* and MP2/6-31G^* levels of theory have predicted the cis conformer as the minimum energy structure for 3-fluoropropene, which is in agreement with the experimental values, whereas the BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory reverses the order of conformation. The ΔE values calculated for 3-chloropropene at MP2/6-31G^*, BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory show that the gauche form is more stable than the cis form, which is in agreement with the experimental value. The same levels of theory have also predicted that the gauche form is stable than cis for 3-bromopropene molecule. The maximum hardness principle has been able to predict the stable conformer of 3-fluoropropene at HF/6-31G^* level of theory, but the same level of theory reverses the conformational stability of 3-chloropropene and 3-bromopropene molecules and MP2/6-31G^* level of theory predicted the stable conformer correctly.     

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     

Vibrational frequencies and infrared intensities of the hydrogen-bonded complexes of nitrous acid with ethers: ab initio and DFT studies

The vibrational spectra of the binary complexes formed by HONO-trans and HONO-cis with dimethyl and diethyl ethers have been investigated using ab initio calculations at the SCF and MP2 levels with 6-311++G(d,p) basis set and B3LYP calculations with 6-31G(d,p) and 6-31+G(d,p) basis sets. Full geometry optimisation was made for the complexes studied. The accuracy of the ab initio calculations have been estimated by comparison between the predicted values of the vibrational characteristics (vibrational frequencies and infrared intensities) and the available experimental data. It was established, that the methods, used in this study are well adapted to the problem under examination. The predicted values with the B3LYP calculations are very near to the results, obtained with 6-311++G(d,p)/MP2. The ab initio and DFT calculations show that the changes in the vibrational characteristics (vibrational frequencies and infrared intensities) upon hydrogen bonding for the hydrogen-bonded complex (CH3)2O...HONO-trans are larger than for the complex (CH3)2O...HONO-cis.     

Peptide models. XXXIII. Extrapolation of low-level Hartree-Fock data of peptide conformation to large basis set SCF,MP2, DFT,and CCSD(T) results. The Ramachandran surface of alanine dipeptide computed at various levels of theory

At the dawn of the new millenium, new concepts are required for a more profound understanding of protein structures. Together with NMR and X-ray-based 3D-structure determinations in silico methods are now widely accepted. Homology-based modeling studies, molecular dynamics methods, and quantum mechanical approaches are more commonly used. Despite the steady and exponential increase in computational power, high level ab initio methods will not be in common use for studying the structure and dynamics of large peptides and proteins in the near future. We are presenting here a novel approach, in which low- and medium-level ab initio energy results are scaled, thus extrapolating to a higher level of information. This scaling is of special significance, because we observed previously on molecular properties such as energy, chemical shielding data, etc., determined at a higher theoretical level, do correlate better with experimental data, than those originating from lower theoretical treatments. The Ramachandran surface of an alanine dipeptide now determined at six different levels of theory [RHF and B3LYP 3-21G, 6-31+G(d) and 6-311++G(d,p)] serves as a suitable test. Minima, first-order critical points and partially optimized structures, determined at different levels of theory (SCF, DFT), were completed with high level energy calculations such as MP2, MP4D, and CCSD(T). For the first time three different CCSD(T) sets of energies were determined for all stable B3LYP/6-311++G(d,p) minima of an alanine dipeptide. From the simplest ab initio data (e.g., RHF/3-21G) to more complex results [CCSD(T)/6-311+G(d,p)//B3LYP/6-311++G(d,p)] all data sets were compared, analyzed in a comprehensive manner, and evaluated by means of statistics.     

丁烯和丁烷自由基阳离子的分子结构和超精细结构的 理论研究

用密度函数理论B3LYP方法和6-31G(d,p),6-311G(d,p)及6-311+G(d,p)基组，分别对1-C4H^+~8,2-C4H^+~8和C4H^+~10进行了构型优化和频率分析计算，预言1-C4H^+~8具有非平面构型，与以往报道的从头算和密度函数理论计算结果不同。在各自由基阳离子的B3LYP构型上，进行了B3LYP、MP2及MRSDCI方法的超精细偶合常数计算，得到了比以往更好的结果，特别是MP2/B3LYP计算值是至今与实验值符合得最好的理论计算结果。     

Density functional calculations of difluorodiazete structures with Gaussian-orbital-type approach

The comparison of local nonlocal and hybrid DFT methods with RHF , MP 2, CCSD , and CPF ab initio methods in generating geometries and relative energies of cis- and trans-difluorodiazete, SVWN , BLYP , BP 86, BECKE 3LYP , and BECKE 3P 86 DFT methods with 6-311 + g(2d) and 6-311 + + g(3df) basis sets. The geometries generated with RHF ab initio models are quite different from experimental values and energy evaluation prefers the wrong isomer. The hybrid methods give superior geometries while energies evaluated with nonlocal DFT methods are better than the one obtained with MP 2 or CCDS ab initio methods. The results suggest DFT as the method of choice of studying similar systems. © 1996 John Wiley & Sons, Inc.     

The structure and energetics of pyrrolidinones, tetrahydrofuranones, piperidinones, and tetrahydropyranones: a computational study

We have performed high level ab initio quantum chemical calculations for 2- and 3-pyrrolidinone; tetrahydro-2- and -3-furanone; 2-, 3-, and 4-piperidinone; and tetrahydro-2-, -3-, and -4-pyranone. The most stable molecular structures were obtained from DFT calculations using the B3LYP density functional and the 6-31G(d) and 6-311+G(3df,2p) basis sets. The respective enthalpies of formation have been computed by the G3(MP2)//B3LYP composite method and appropriately chosen reactions. The calculated results are in excellent agreement with experimental data reported in the literature.     

The evaluation of bond dissociation energies for simple sulphur containing molecules using ab initio and density functional methods

The S–H and C–S bond dissociation energies for simple alkylthiols and dialkylsulphides, along with the S–S bond dissociation energy for dimethyl disulphide, compounds which have been used in the metal–organic chemical vapour deposition (MOCVD) growth of wide band gap II–VI (12–16) Zn- and Cd-based compound semiconductors, have been computed using the ab initio (ROHF and MP2) and density functional theory (DFT) methods (BHandH, BHandHLYP, B3LYP, B3P86, B3PW91, BLYP and BP86) with the 6-311+G(2d,p) basis set along with high accuracy complete basis set, CBS-4 and CBS-Q energy computations. The computed energies are compared with experimental results and the suitability of the DFT methods, for the computational study of these systems, is discussed.     

FTIR, FT-Raman spectra and ab initio DFT vibrational analysis of 2-bromo-4-methyl-phenylamine

The FTIR and FT-Raman spectra of 2-bromo-4-methyl-phenylamine (BMP) have been recorded. From the standard geometrical parameters the geometry of BMP was optimized at ab initio and DFT levels of theory with complete relaxation in the potential energy surface using 6-311+g(d,p) and 6-311+g(2df,2p) basis sets. Several thermodynamic parameters were also calculated for the minimum energy conformer at ab initio and DFT level of theories. The harmonic vibrational frequencies were calculated and the scaled values have been compared with experimental FTIR and FT-Raman spectra. Majority of the computed wavenumbers were found to be in good agreement with the experimental observations. The experimental spectra also coincide satisfactorily with those of theoretically constructed bar type spectrograms.     

FT-IR, FT-Raman vibrational spectra and molecular structure investigation of 2-chloro-4-methylaniline: A combined experimental and theoretical study

In this work, the experimental and theoretical vibrational spectra of 2-chloro-4-methylaniline (2Cl4MA, C₇H₈NCl) were studied. FT-IR and FT-Raman spectra of 2Cl4MA in the liquid phase have been recorded in the region 4000–400 cm⁻¹ and 3500–50 cm⁻¹, respectively. The structural and spectroscopic data of the molecule in the ground state have been calculated by using Hartree-Fock (HF) and density functional method (B3LYP) with the 6-31G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p) and 6-311G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p) basis sets. The vibrational frequencies have been calculated and scaled values have been compared with experimental FT-IR and FT-Raman spectra. The observed and calculated frequencies are found to be in good agreement. The complete assignments were performed on the basis of the total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. The DFT-B3LYP/6-311++G(d,p) calculations have been found more reliable than the ab initio HF/6-311++G(d,p) calculations for the vibrational study of 2Cl4MA. The optimized geometric parameters (bond lengths and bond angles) were compared with experimental values of aniline and p-methylaniline molecules.     

FTIR and FTRaman spectra,assignments, ab initio HF and DFT analysis of 4-nitrotoluene

In this work, the experimental and theoretical study on molecular structure and vibrational spectra of 4-nitrotoluene are studied. The FTIR and FTRaman experimental spectra of the molecule have been recorded in the range of 4000–100 cm^?1. Making use of the recorded data, the complete vibrational assignments are made and analysis of the observed fundamental bands of molecule is carried out. The experimental determinations of vibrational frequencies are compared with those obtained theoretically from ab initio HF and DFT quantum mechanical calculations using HF/6-31G (d, p), B3LYP/6-31++G* (d, p) and B3LYP/6-311++G* (d, p) methods. The differences between the observed and scaled wave number values of most of the fundamentals are very small in B3LYP than HF. The geometries and normal modes of vibrations obtained from ab initio HF and B3LYP calculations are in good agreement with the experimentally observed data. Comparison of the simulated spectra provides important information about the ability of the computational method (B3LYP) to describe the vibrational modes. The vibrations of NO₂ and CH₃ groups coupled with skeletal vibrations are also investigated.     

戊烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p)、6-31＋G(d,p)、6-311G(d,p)及6-311＋G(d,p)基组，分别对2-C5H10＋和1-C5H10＋的各种构象进行了几何构型优化，并用B3LYP/6-311G(d,p)进行了频率分析计算.计算预言1-C5H10＋具有非平面构型，与以往报导的从头算计算结论相反.在两个自由基阳离子的各种构象的B3LYP几何构型上，进行了B3LYP和UMP2(full)方法的超精细偶合常数计算，得到了比以往更好的结果.     

Electron affinities of metals computed by density functional theory and ab initio methods

An extensive computational study of the meal electron affinity was performed using the ab initio and density functional theory (DFT) methods. HF, MP2, MP3, MP4, QCISD, and QCISD(T) was used as computational methods, while the hybrid, local, and nonlocal DFT methods with the LYP, P86, PW91, and VWN correlation functionals were used. Two basis sets, one small and applicable to almost all metals (LanL2DZ) and one large [6-311 + + G(3df, 3 pd)] used only for small metals, were employed. The computed results were compared with the experimental data and the capabilities of the DFT methods to perform this study were discussed. © 1997 John Wiley & Sons, Inc.     

CH and Chalogen bond dissociation energies for fluorinated and chlorinated methane evaluated with hybrid B3LYP density functional theory methods and their comparison with experimental data and the CBS-Q ab initio computational approach

The CH and Chalogen bond dissociation energies (BDEs) were computed with the hybrid B3LYP/6-311 + G(2d,2p) theory model for chlorinated and fluorinated methane. All computed values were substantially lower (5–10 kcal mol⁻¹) than the experimental values. To obtain better agreement, a correlation factor was introduced. When this factor was applied, excellent agreement between the B3LYP/6-311 + G(2d,2p) computed energies and the experimental BDEs was observed. On the other hand, the CBS-Q ab initio computational approach generated BDEs which are in good agreement with experimental values without a correction factor.     

High level of ab initio and density functional theory computational study of structural and energetic properties of tetrafluorhydrazine rotamers

The HF, MP2, MP3, MP4, and QCISD ab initio methods were compared with local, hybrid, and gradient-corrected density functional theory (DFT) methods for computing structures and energies of N2F4 rotamers. In all DFT calculations 6-311 + G(2d) basis set was used. The generated structures energies of trans- and gauche-N₂F₄ rotamers, and their dissociation energies to nitrogen difluoride were compared with experimental data. Suitable hybrid and gradient-corrected DFT methods for determining structures and energies for these and similar molecular systems were discussed.     

Kinetics of the multichannel reaction of methanethiyl radical (CH3S*) with 3O2

The CH3S* + O2 reaction system is considered an important process in atmospheric chemistry and in combustion as a pathway for the exothermic conversion of methane-thiyl radical, CH3S*. Several density functional and ab initio computational methods are used in this study to determine thermochemical parameters, reaction paths, and kinetic barriers in the CH3S* + O2 reaction system. The data are also used to evaluate feasibility of the DFT methods for higher molecular weight oxy-sulfur hydrocarbons, where sulfur presents added complexity from its many valence states. The methods include: B3LYP/6-311++G(d,p), B3LYP/6-311++G(3df,2p), CCSD(T)/6-311G(d,p)//MP2/6-31G(d,p), B3P86/6-311G(2d,2p)//B3P86/6-31G(d), B3PW91/6-311++G(3df,2p), G3MP2, and CBS-QB3. The well depth for the CH3S* + 3O2 reaction to the syn-CH3SOO* adduct is found to be 9.7 kcal/mol. Low barrier exit channels from the syn-CH3SOO* adduct include: CH2S + HO2, (TS6, E(a) is 12.5 kcal/mol), CH3 + SO2 via CH3SO2 (TS2', E(a) is 17.8) and CH3SO + O (TS17, E(a) is 24.7) where the activation energy is relative to the syn-CH3SOO* stabilized adduct. The transition state (TS5) for formation of the CH3SOO adduct from CH3S* + O2 and the reverse dissociation of CH3SOO to CH3S* + O2 is relatively tight compared to typical association and simple bond dissociation reactions; this is a result of the very weak interaction. Reverse reaction is the dominant dissociation path due to enthalpy and entropy considerations. The rate constants from the chemical activation reaction and from the stabilized adduct to these products are estimated as functions of temperature and pressure. Our forward rate constant and CH3S loss profile are in agreement with the experiments under similar conditions. Of the methods above, the G3MP2 and CBS-QB3 composite methods are recommended for thermochemical determinations on these carbon-sulfur-oxygen systems, when they are feasible.