Performance of theoretical methods and basis sets on the molecular structure, atomisation and ionisation energies, electron affinity, and vibrational spectrum of silylene

This work compares the performance of theoretical methods and basis sets on the molecular structure, atomisation and ionisation energies, electron affinity, and vibrational spectrum of silylene. Silylene, its cation and anion have been studied in ¹ A ₁, ² A ₁ and ² B ₁ states, respectively, in the gas phase and C_2v symmetry. The methods considered are second-order Møller-Plesset perturbation theory (MP2), the density functional theory (DFT), Gaussian-2 (G2) and complete basis set methods (CBS-4M and CBS-Q). The basis sets used are 6-31G(d,p), 6-311G(d,p), 6-31++G(d,p) and 6-311++G(d,p). The functional used for the DFT method is B3LYP. Silylene and its cation and anion have been optimised using the MP2 and DFT methods and the named basis sets. Single-point energy calculations (G2, CBS-4M and CBS-Q) were performed using MP2/6-311++G(d,p) structures and these energies have been used to calculate atomisation energy, ionisation energy and adiabatic electron affinity. Frequency calculations were also done and the raw vibrational frequencies were assigned. It is interesting to note the close similarity between the predicted parameters and some of the available literature values. The results obtained are consistent and converge with different basis sets with improved size and quality. However, the parameters obtained are very much method dependent.     

Structural and spectroscopic studies on 2-pyranones

Experimental methods of infrared, Raman and electronic absorption spectroscopy and DFT calculations using B3LYP functionals and 6-31G** and 6-311++G** basis sets have been used to understand the structural and spectral characteristics of 2-pyranones, 6-phenyl-4-methylsulfanyl-2-oxo-2H-pyran and 6-phenyl-4-methylsulfanyl-2-oxo-2H-pyran-3-carbonitrile in the electronic ground (S₀) and first excited (S₁) states. Information about the size, shape, charge density distribution and site of chemical reactivity of the molecules has been obtained by mapping electron density isosurface with electrostatic potential surfaces (ESP). Based on TD-DFT calculations using 6-31+G**5D basis set, an assignment of absorption peaks in the UV–VIS region has been suggested. The S₁ state is found to be a ¹(π,π*) state. A complete vibrational analysis has been attempted on the basis of experimental infrared and Raman spectra and calculated frequency and intensity of the vibrational bands and potential energy distribution over the internal coordinates. Characteristic vibrational bands of the 2-pyranone ring and methylsulfanyl and carbonyl groups have been identified.     

Comparative study of DFT methods applied to small titanium/oxygen compounds

The performance of a number of different local and nonlocal density functional theory (DFT) methods has been investigated for some small titanium—oxygen systems. Equilibrium geometries, ionization potentials, dipole moments, atomization energies, and harmonic vibrational frequencies have been calculated for the TiO, TiO₂, and Ti₂ molecules, and the results are compared with experimental data and ab initio calculations. It is shown that most DFT methods perform much better than the ab initio Hartree—Fock (HF), second-order perturbation theory (MP2), and configuration interaction including single and double excitations (CISD) treatments. For good agreement with experimental data, gradient corrections to the exchange part of the DFT functional are needed, as well as some type of correction for the errors in the calculated energy splittings between different atomic states of titanium. Hybrid methods including a mixture of HF exchange with DFT exchange correlation do not perform as well as “pure” DFT methods for the studied systems. © 1996 John Wiley & Sons, Inc.     

Theoretical study of structure,vibrational and electronic spectra of isomers of methyl-3-methoxy-2-propenoate

A systematic quantum mechanical study of the possible conformations, their relative stabilities, vibrational and electronic spectra and thermodynamic parameters of methyl-3-methoxy-2-propenoate has been reported for the electronic ground (S₀) and first excited (S₁) states using time-dependent and time-independent Density Functional Theory (DFT) and RHF methods in extended basis sets. Detailed studies have been restricted to the E-isomer, which is found to be substantially more stable than the Z-isomer. Four possible conformers c′Cc, c′Tc, t′Cc, t′Tc, of which the first two are most stable, have been identified in the S₀ and S₁ states. Electronic excitation to S₁ state is accompanied with a reversal in the relative stability of the c′Cc and c′Tc conformers and a substantial reduction in the rotational barrier between them, as compared with the S₀ state. Optimized geometries of these conformers in the S₀ and S₁ states are being reported. Based on suitably scaled RHF/6-31G^** and DFT/6-311G^** calculations, assignments have been provided to the fundamental vibrational bands of both these conformers in terms of frequency, form and intensity of vibrations and potential energy distribution across the symmetry coordinates in the S₀ state. A complete interpretation of the electronic spectra of the conformers has been provided.     

Theoretical Gas Phase Study of the Gauche and Trans Conformers of 1-Bromo-2-Iodoethane and Solvent Effects

This is a gas-phase study of the gauche and trans conformers of 1-bromo-2-iodoethane. The methods used are the second-order Møller-Plesset theory (MP2) and density functional theory (DFT). The functional used for the DFT method is B3LYP and the basis sets used are 6-311++G(d,p) for all atoms except that different basis sets, namely 3-21G, LANECP, CRENBL ECP, Stuttgart RLC ECP and 6-311G(d,p), have been explored for the iodine atom. The results indicate that the trans conformer is preferred. The energy difference between the gauche and trans conformers (ΔE _g?t) and related thermodynamic parameters are reported. The ΔE _g?t values are 12.50 kJ?mol^?1 (B3LYP) and 10.00 kJ?mol^?1 (MP2) with the basis sets being 6-311++G(d,p)[C,H,Br]/6-311G(d,p)[I]. The conformers of 1-bromo-2-iodoethane have also been subjected to vibrational analysis. The results from the two theoretical levels are in good agreement but they are not much affected by the basis set of the iodine atom. The study has been extended to explore solvent effects using Self-Consistent Reaction Field methods. The structural parameters of the conformers are little affected by the polarity of the solvent but ΔE _g?t decreases and the solvation Gibbs energy increases with increasing polarity of the solvent.     

Synthesis, spectroscopic and DFT investigation of dimethyl-2-(5-acetyl-2,2-dimethyl-4,6-dioxo-1,3-dioxan-5-yl)-3-(triphenylphosphinylidene)succinate

Dimethyl-2-(5-acetyl-2,2-dimethyl-4,6-dioxo-1,3-dioxan-5-yl)-3-(triphenylphosphinylidene)succinate has been synthesized and characterized by elemental analysis, FT-IR and ¹H, ¹³C and ³¹P NMR. The vibrational wavenumbers, gauge including atomic orbital (GIAO) ¹H and ¹³C chemical shift values of title compound in the ground state have been computed with density functional theory method (DFT) and the B3LYP functional. The basis sets used are 6-311G(d,p) and 6-31G(d). The harmonic vibrational wavenumbers have been computed and the scaled values have been compared with the experimental FT-IR spectra. The complete assignments have been performed on basis of the total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. Most of the computed wavenumbers are found to be in good agreement with the observed spectrum.     

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     

Quantum DFT and DF–DFT study of vibrational spectra of sulfuric acid,sulfuric acid monohydrate,formic acid and its cyclic dimer

Rigorous theoretical treatment of vibrational frequencies is critically important for the interpretation of unassigned experimental vibrational spectra and accurate determination of thermodynamic properties of molecular clusters. IR spectra of trans monomers of sulfuric and acetic acids, sulfuric acid monohydrate and cyclic dimer of the formic acid have been studied using DFT and DF–DFT methods using BLYP, B3LYP and PW91 with 12 different Pople and Dunning basis sets. New data for above-mentioned structures have been reported, scaling factors have been calculated and a comprehensive analysis of the performance of BLYP, B3LYP and PW91 methods has been performed. Comparison of the obtained results with experiments shows that results of pure PW91 and BLYP are better than predictions of well-established hybrid B3LYP method. Our analysis shows on the existence of the considerable difference in scaling factors weighted to high and low frequencies. In the case of formic acid dimer, the deviation the predicted low frequencies from the experimental data is considerable that leads to quite large (∼6–7 kcal mol⁻¹) uncertainties in the absolute values of the cluster Gibbs free energy. In order to determine an efficient computational strategy that comprises accuracy and reasonable computational costs, the effect of density fitting (DF) and basis set superposition error (BSSE) on the vibration frequencies has been investigated. We found that application of the DF that substantially (2.5–3.5 times) increases the performance of pure PW91 and BLYP methods gives excellent results, which are very close to those of conventional DFT. This suggests that DF–DFT is a viable low-cost alternative to conventional DFT in predicting vibrational spectra. It has been found that while vibrational spectra obtained using the counterpoise correction (CP) for the BSSE do not deviate much from uncorrected ones, the difference in absorption intensities between corrected and uncorrected spectra obtained using small and medium-sized basis sets is considerable. This suggests that application of DF–DFT uncorrected for the BSSE with large basis sets is a more efficient strategy of predicting vibrational spectra than the application of conventional DFT with small basis sets.     

An ab initio theoretical study of the optical spectrum of CF2

The ab initio calculation methods have been used to calculate the spectral and electronic characteristics of difluorocarbene in the ground electronic state (¹A₁), the lowest-lying singlet (¹B₁) and triplet (³B₁) states. The optimized equilibrium geometries, rotational constants, harmonic vibrational frequencies and energy gaps, electronic charges, dipole moments of these states have been computed with different basis sets. The calculated vibrational frequency of ³B₁ state (v₂=522 cm^?1) and the energy separation (2.26 eV) between ³B₁ and ¹A₁ states are in good agreement with the experimental results (519 cm^?1, 2.46 eV respectively). According to the calculations the previous assignment of vibrational symmetries of ¹B₁ state was corrected, and some experimentally undetermined vibrational frequencies were predicted.     

Vibrational Spectra and Quantum Calculations of Ethylbenzene

Normal vibrations of ethylbenzene in the first excited state have been studied using resonant two-photon ionization spectroscopy. The band origin of ethylbenzene of S₁←S₀ transition appeared at 37586 cm^-1. A vibrational spectrum of 2000 cm^-1 above the band origin in the first excited state has been obtained. Several chain torsions and normal vibrations are obtained in the spectrum. The energies of the first excited state are calculated by the time-dependent density function theory and configuration interaction singles (CIS) methods with various basis sets. The optimized structures and vibrational frequencies of the S₀ and S₁ states are calculated using Hartree-Fock and CIS methods with 6-311++G(2d,2p) basis set. The calculated geometric structures in the S₀ and S₁ states are gauche conformations that the symmetric plane of ethyl group is perpendicular to the ring plane. All the observed spectral bands have been successfully assigned with the help of our calculations.     

Density functional and multiconfigurational ab initio study of the ground and excited states of Os2

Multiconfigurational ab initio methods predict that the ⁵Π_u state as the ground state instead of the ⁷Δ_u state. Although multiconfigurational perturbation theory correctly predicts the ground state, they overestimate the bond dissociation energy (BDE). Only multireference configuration interaction method can reasonably calculate the BDE. The spin‐orbit effect on the spectroscopic constants is not significant. The results calculated by density functional theory (DFT) vary significantly depending on the selection of a DFT functional. No DFT functional gives the same energy ordering as calculated by the second‐order multiconfigurational perturbation theory (CASPT2). The old generalized gradient approximations functionals are well suited for predicting the ground state and calculating the bond length and the vibrational frequency of Os₂. According to the CASPT2 calculation, the ground state of Os₂ has a quadruple bond. © 2014 Wiley Periodicals, Inc.     

FT-IR and FT-Raman, UV spectroscopic investigation of 1-bromo-3-fluorobenzene using DFT (B3LYP, B3PW91 and MPW91PW91) calculations

The FT-IR and FT-Raman spectra of 1-bromo-3-fluorobenzene (C₆H₄FBr) molecule have been recorded using Bruker IFS 66 V spectrometer in the range of 4000–100 cm⁻¹. The molecular geometry and vibrational frequencies in the ground state are calculated using the DFT (B3LYP, B3PW91 and MPW91PW91) methods with 6-31++G(d,p) and 6-311++G(d,p) basis sets. The computed values of frequencies are scaled using a suitable scale factor to yield good coherence with the observed values. The isotropic DFT (B3LYP, B3PW91 and MPW1PW91) analysis showed good agreement with the experimental observations. Comparison of the fundamental vibrational frequencies with calculated results by B3LYP methods. The complete data of this molecule provide the information for future development of substituted benzene. The influence of bromine and fluorine atom on the geometry of benzene and its normal modes of vibrations has also been discussed. A study on the electronic properties, such as absorption wavelengths, excitation energy, dipole moment and frontier molecular orbital energies, was performed by time dependent DFT (TD-DFT) approach. The electronic structure and the assignment of the absorption bands in the electronic spectra of steady compounds were discussed. The calculated HOMO and LUMO energies show that charge transfer occurs within the molecule. On the basis of the thermodynamic properties of the title compound at different temperatures have been calculated in gas phase, revealing the correlations between standard heat capacities (C) standard entropies (S), standard enthalpy changes (H) and temperatures.     

Quantum‐chemical,IR, NMR,and X‐ray diffraction studies on 2‐(4‐chlorophenyl)‐1‐methyl‐ 1H‐benzo[d]imidazole

The title molecule, 2‐(4‐chlorophenyl)‐1‐methyl‐1H‐benzo[d]imidazole (C₁₄H₁₁ClN₂), was prepared and characterized by ¹H NMR, ¹³C NMR, IR, and single‐crystal X‐ray diffraction. The molecular geometry, vibrational frequencies, and gauge including atomic orbital (GIAO) ¹H and ¹³C NMR chemical shift values of the title compound in the ground state have been calculated by using the Hartree‐Fock (HF) and density functional theory (DFT/B3LYP) method with 6‐31G(d) basis sets, and compared with the experimental data. The calculated results show that the optimized geometries can well reproduce the crystal structural parameters, and the theoretical vibrational frequencies and GIAO ¹H and ¹³C NMR chemical shifts show good agreement with experimental values. The energetic behavior of the title compound in solvent media has been examined using B3LYP method with the 6‐31G(d) basis set by applying the Onsager and the polarizable continuum model (PCM). Besides, molecular electrostatic potential (MEP), frontier molecular orbitals (FMO) analysis, and nonlinear optical (NLO) properties of the title compound were investigated by theoretical calculations. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

An improved generator coordinate Hartree–Fock method applied to the choice of contracted Gaussian basis sets for first‐row diatomic molecules

Accurate Gaussian basis sets (18s for Li and Be and 20s11p for the atoms from B to Ne) for the first‐row atoms, generated with an improved generator coordinate Hartree–Fock method, were contracted and enriched with polarization functions. These basis sets were tested for B₂, C₂, BeO, CN⁻, LiF, N₂, CO, BF, NO⁺, O₂, and F₂. At the Hartree–Fock (HP), second‐order Møller–Plesset (MP2), fourth‐order Møller–Plesset (MP4), and density functional theory (DFT) levels, the dipole moments, bond lengths, and harmonic vibrational frequencies were studied, and at the MP2, MP4, and DFT levels, the dissociation energies were evaluated and compared with the corresponding experimental values and with values obtained using other contracted Gaussian basis sets and numerical HF calculations. For all diatomic molecules studied, the differences between our total energies, obtained with the largest contracted basis set [6s5p3d1f], and those calculated with the numerical HF methods were always less than 3.2 mhartree. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 15–23, 2000     

Basis-set effects in quantum-mechanical calculations of translational-to-vibrational energy transfer in O(3P) + CO2 collisions

Cross sections for the vibrational excitation of CO₂ in collisions with O(³P) from the ground state to the (00⁰1) state are computed using vibrational-close-coupling with the infinite-order-sudden approximation for rotation (VCC IOS). For a relative translational energy of 3 eV, several vibrational basis sets are tested for convergence of this cross section. Using the best converged vibrational basis set for 3 eV, the (00⁰1) → (00⁰1) cross section is computed over the relative translational energy range of 2 to 6 eV. Previous discrepancies of factors of 2 to 7 between the VCC IOS method and quasiclassical trajectory studies are reduced to factors of 1.1 to 2.5.     

Reaction mechanism of HSH and CH3SH with NH2CH2COCH2X (X = F and Cl) molecules

The reaction mechanism of model compounds H₂S and CH₃SH for cysteine proteases with NH₂CH₂COCH₂X (X = F and Cl) molecules has been investigated using DFT methods with B3LYP and B3PW91 hybrid density functionals at 6‐31+G* basis sets. The single point energy has been calculated for the above reactions with B3LYP and B3PW91 functionals using aug‐cc‐PVDZ infinite basis set in both gas and solution phases. The intrinsic reaction coordinates calculations have been performed to confirm that each transition state is linked by the desired reactants and products. The geometries and relative energies for various stationary points have been determined and discussed. The zero point vibrational energy corrections have been made to predict the reliable energy. The negative value of reaction energy indicates that the overall reaction profile is found to be exothermic. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Theoretical investigation of the isomers photolysis reaction for chlorine nitrate

For the reactant ClONO₂ the equilibrium geometry and energy have been calculated by using various density functional (DFT) theory methods. The harmonic vibrational frequencies were computed using the B3LYP method for the ground state of the reactant; the results were in good agreement with those of the experiment. It was shown that the DFT methods were superior to other ab initio methods in optimizing geometry and predicting vibrational frequencies. In the process of forming the products from the dissociation reactions of ClONO₂, the results indicated it undergoes the TS₂, but not the TS₁ transition state. According to the relative energies of various species, the potential energy surface reflected intuitively the mechanisms of these dissociation reactions. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 44–50, 2003     

Theoretical study of the vibrational frequencies of carbon disulfide

A benchmark comparison for different computational methods and basis sets has been presented. In this study, five computational methods (Hartree–Fock (HF), MP2, B3LYP, MPW1MP91, and PBE1PBE) along with 18 basis sets have been applied to optimize the geometry of carbon disulfide (CS₂), and further calculate the vibrational frequencies of the optimized geometries. The differences between the calculated frequencies and corresponding experimental data are used to evaluate the efficiency of each combination of computational method and basis set. The comparison of frequency difference indicates that B3LYP generally gives the best prediction of frequencies for CS₂, whereas the other two density functional theory (DFT) methods, i.e., MPW1PW91 and PBE1PBE, often give parallel results. Although MP2 predicts the frequencies with accuracy almost as good as those from DFT methods, in a particular case, HF calculation outperforms MP2 as well as MPW1PW91 and PBE1PBE for prediction of the frequency of asymmetrical stretching for CS₂. © 2013 Wiley Periodicals, Inc.     

Synthesis,molecular conformation,vibrational, electronic transition,and chemical shift assignments of 4-(thiophene-3-ylmethoxy)phthalonitrile: a combined experimental and theoretical analysis

This work presents the synthesis and characterization of a novel compound, 4-(thiophene-3-ylmethoxy)phthalonitrile (TMP). The spectroscopic properties of the compound were examined by FT-IR, FT-Raman, NMR, and UV techniques. FT-IR and FT-Raman spectra in solid state were observed in the region 4000–400 cm⁻¹ and 3500–50 cm⁻¹, respectively. The ¹H and ¹³C NMR spectra were recorded in CDCl₃ solution. The UV absorption spectrum of the compound that dissolved in THF was recorded in the range of 200–800 nm. The structural and spectroscopic data of the molecule in the ground state were calculated using density functional theory (DFT) employing B3LYP exchange correlation and the 6-311++G(d,p) basis set. The vibrational wavenumbers were calculated and scaled values were compared with experimental FT-IR and FT-Raman spectra. The complete assignments were performed on the basis of the experimental results and total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. Isotropic chemical shifts (¹³C NMR and ¹H NMR) were calculated using the gauge-invariant atomic orbital (GIAO) method. A study on the electronic properties, such as HOMO and LUMO energies, were performed by time-dependent DFT (TD-DFT) approach. The HOMO and LUMO analyses have been used to elucidate information regarding charge transfer within the molecule. Comparison of the calculated frequencies, NMR chemical shifts, absorption wavelengths with the experimental values revealed that DFT method produces good results.     

Accurate theoretical IR and Raman spectrum of Al2O2 and Al2O3 molecules

The accurate harmonic vibration frequencies together with the infrared (IR) and Raman intensities of the most stable conformers of Al₂O₂ and Al₂O₃ molecules have been calculated by the density functional theory (DFT) method with B3LYP exchange–correlation potential and using a set of the augmented correlated consistent basis sets up to quintuple order. The anharmonic vibration frequencies of the non-linear Al₂O₂ molecule have also been calculated. The obtained equilibrium geometrical parameters, harmonic and anharmonic vibration frequencies along with the IR and Raman intensities good converge to their limits with increasing the size of the used basis set. A comparison of the calculated harmonic and anharmonic vibrational frequencies with the available experimental ones points out that the small differences between the calculated harmonic and experimental frequencies can be further substantially reduced when calculations of the anharmonic vibrational frequencies will be available for all types of molecular geometries.