Calculation of bond dissociation energies for oxygen containing molecules by ab initio and density functional theory methods

Two ab initio (ROHF and MP2), one local (SVWN), four hybrid (BHandH, BHandHLYP, Becke3LYP, and Becke3P86), and two nonlocal (BLYP and BP86) density functional theory (DFT) methods are used for calculating the dissociation energies of molecules that contain H(SINGLE BOND)O, O(SINGLE BOND)O and O(SINGLE BOND)C bonds. The sensitivity to the basis set of the prediction of bond dissociation energies with DFT methods was tested with Becke3LYP on the H(SINGLE BOND)O dissociation energy of water. The 6–31 + G(d) methods are chosen as the smallest basis set which produces reasonable results. The calculated values for all other ab initio and DFT methods were performed with these basis sets and then compared with the experimental data. The suitability of DFT methods for computing reliable bond dissociation energies of oxygen containing molecules is discussed. © 1996 John Wiley & Sons, Inc.     

Computation of bond dissociation energy for sulfides and disulfides with ab initio and density functional theory methods

The geometries and S-H, S-S, and S-C bond dissociation energies for hydrogen sulfide, hydrogen disulfide, methanethiol, dimethyl disulfide, and dimethyl disulfide were calculated with both ab initio (ROHF and MP2), hybrid (BHandH, BHandHLYP, Becke3LYP and Becke3P86), and nonlocal (BLYP and BP86) density functional theory (DFT) methods. In all studies the 6–31 + G(d) basis set is used. The computed results are compared to the experimentally obtained values, targeting the selection of a suitable ab initio or DFT method for the study of these systems. © 1997 John Wiley & Sons, Inc.     

Electron localization function and electron localizability indicator applied to study the bonding in the peroxynitrous acid HOONO

The ground‐state electronic structure of peroxynitrous acid (HOONO) and its singlet biradicaloid form (HO ··· ONO) have been studied using topological analysis of the electron localization function (ELF), together with the electron localizability indicator (ELI‐D), at the DFT (B3LYP, M05, M052X, and M06), CCSD, and CASSCF levels. Three isomers of HOONO (cis‐cis, cis‐perp, and trans‐perp) have been considered. The results show that from all functionals applied, only B3LYP yields the correct geometrical structure. The ELF and ELI‐D‐topology of the O? O and central N? O bonds strongly depends on the wave function used for analysis. Calculations carried out at CAS (14,12)/aug‐cc‐pVTZ//CCSD(T)/aug‐cc‐pVTZ level reveal two bonds of the charge‐shift type: a protocovalent N? O bond with a basin population of 0.82–1.08e, and a more electron depleted O? O bond with a population of 0.66–0.71e. The most favorable dissociation channel (HOONO → HO + ONO) corresponds to breaking of the most electron‐deficient bond (O? O). In the case of cis‐ and trans‐HO ··· ONO, the ELF, ELI‐D, and electron density fields results demonstrate a closed‐shell O ··· O interaction. The α‐spin electrons are found mainly (0.64e) in the lone pairs of oxygen V_{i = 1,2} (O) from the OH group. The β‐spin electrons are delocalized over the ONO group, with the largest concentration (0.34e) on the lone pair of nitrogen V(N). © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Density functional theory study on molecular structure and vibrational IR spectra of isocytosine

Density functional theory with the combined Becke3-LYP exchange-correlation energy functional [DFT(B3-LYP) method] using the 6-31G(d, p) basis set is applied to predict molecular parameters (geometries, rotational constants, dipole moments) and vibrational IR spectra (harmonic wavenumbers, absolute intensities) of six tautomers of the isocytosine molecule. The results are compared with the corresponding data calculated at the conventional ab initio Hartree-Fock (HF) level using the same basis set and with available experimental data. Calculations show that (a) three amino tautomers are slightly nonplanar species with, evidently, a distorted amino group, (b) the DFT (B3-LYP)/6-31G(d, p) method predicts better molecular parameters, than do the HF calculations, and (c) the DFT(B3-LYP)-calculated vibrational IR spectra of isocytosine agree well with the available recorded IR spectra, and they show marked improvement over the IR spectra predicted at the HF/6-31G(d, p) level. Tautomeric stabilities of isocytosine are discussed on the basis of computed electronic energies by the DFT(B3-LYP) and ab initio approaches [including the MP2 and MP4(SDQ) calculations of electronic energies] and predicted zero-point vibrational energies by DFT(B3-LYP) and HF methods. This relative energies at 0 K of the tautomeric forms of isocytosine predicted by both conventional ab initio and DFT(B3-LYP) methods correlate well with the experimental data, showing the predominance of the aminohydroxy tautomer of isocytosine for an isolated molecule. © 1997 John Wiley & Sons, Inc.     

Density functional Gaussian-type orbital approach in theoretical study of S2F2 isomerization

The structures of two isomers, difluorodisulfane (FSSF) and thiothionylfluoride (SSF₂), and the corresponding transition structure were generated with density functional theory (DFT) methods. Three groups of DFT methods were used: local (Local Spin Density Approximation, LSDA), nonlocal (local with gradient corrections; BLYP and BP86), and hybrid methods that include a mixture of Hartree-Fock (HF) exchange with nonlocal correlation (Becke3BLYP, Becke3P86). An extended basis set [6-311 + + G(3df)] was used for all calculations, although satisfactory results can be obtained with the 6-311G(d) basis set. The geometries obtained were compared with both restricted Hartree-Fock (RHF) calculated and experimentally obtained values. The energy outcome and the activation barrier for the isomerization were evaluated. It was determined that excellent geometries can be obtained with the Becke3B86 hybrid method, whereas for reasonable energies MP2 single-point calculations on these geometries are necessary. © 1996 by John Wiley & Sons, Inc.     

Theoretical investigation of one- and two-photon absorption properties of platinum acetylide chromophores

In this paper, we have theoretically investigated bis((4-phenylethynyl)phenyl) ethynyl)bis(trimethylphosphine)platinum(II) (PE2) and its analogs three platinum acetylide complexes (1-3) that feature highly pi-conjugated ligands (alkynyl-dimethylfluorene substituted with electron-donating or -withdrawing moieties). The geometrical and electronic structures are calculated at the ECP60MWB//6-31G*(H, C, P, N, S) basis set level by the density functional theory (DFT) method; one-photon absorption properties have been calculated by using time-dependent DFT (TDDFT) and Zerner's intermediate neglect of differential overlap (ZINDO) methods, and two-photon absorption (TPA) properties are obtained with the ZINDO/sum-over-states method. The values of beta(sp) and beta(d) for Pt are adjusted to -1 eV and -28.5 eV, respectively, to make one-photon absorption spectra calculated by ZINDO closest to the experimental data and TDDFT results. The calculated results indicate that all molecules in this work (involving cis isomers of molecules 1-3) take on two TPA peaks in the 600-800 nm region. The peak at 700-750 nm should not be simply attributed to the appearance of noncentrosymmetric cis isomers in solution, although trans and cis isomers adhere to a different selection rule. Every TPA peak results from its transition character. Molecules 1-3 show greater two-photon absorption strength compared with PE2 and retain good transparency.     

ClC(O)NCS分子振动光谱的理论研究

采用密度泛函理论DFT(B3LYP)方法,以6-31G^*为基组对ClC(O)NCS的反式和顺式两种构型的几何结构、振动谐性力场和红外光谱进行了研究.B3LYP/6-31G^*的理论力场由适用于B3LYP/6-31G^*计算水平和大多数有机分子的一套固定标度因子进行标度.根据标度后的理论力场进行简正坐标分析得到的势能分布(PED)和红外光谱强度值对ClC(O)NCS分子的顺式和反式两种构型的振动基频进行了理论归属.     

Theoretical investigation of the structure and vibrational spectra of carbamoyl azide

Molecular structure and vibrational frequencies of carbamoyl azide NH2CO-NNN have been investigated with ab initio and density functional theory (DFT) methods. The molecular geometries for all the possible conformers of the molecule were optimized using DFT-B3LYP, DFT-BLYP and MP2 applying the standard 6-311++G** basis set. From the calculations, the molecule was predicted to exist predominantly in cis conformation with the cis-trans rotational barrier of about 7.91-9.10 kcal/mol depending on the level of theory applied. The vibrational frequencies and the corresponding vibrational assignments of carbamoyl azide in Cs symmetry were examined theoretically and the calculated Infrared and Raman spectra of the molecule in the cis conformation were plotted. Observed frequencies for normal modes were compare 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.     

Accurate interaction energies of hydrogen-bonded nucleic acid base pairs

Hydrogen-bonded nucleic acids base pairs substantially contribute to the structure and stability of nucleic acids. The study presents reference ab initio structures and interaction energies of selected base pairs with binding energies ranging from -5 to -47 kcal/mol. The molecular structures are obtained using the RI-MP2 (resolution of identity MP2) method with extended cc-pVTZ basis set of atomic orbitals. The RI-MP2 method provides results essentially identical with the standard MP2 method. The interaction energies are calculated using the Complete Basis Set (CBS) extrapolation at the RI-MP2 level. For some base pairs, Coupled-Cluster corrections with inclusion of noniterative triple contributions (CCSD(T)) are given. The calculations are compared with selected medium quality methods. The PW91 DFT functional with the 6-31G basis set matches well the RI-MP2/CBS absolute interaction energies and reproduces the relative values of base pairing energies with a maximum relative error of 2.6 kcal/mol when applied with Becke3LYP-optimized geometries. The Becke3LYP DFT functional underestimates the interaction energies by few kcal/mol with relative error of 2.2 kcal/mol. Very good performance of nonpolarizable Cornell et al. force field is confirmed and this indirectly supports the view that H-bonded base pairs are primarily stabilized by electrostatic interactions.     

An ab initio MO study on the structures and electronic states of hydrogen-bonded O3- HF and SO2-HF complexes

Ab initio molecular orbital (MO) calculations have been carried out for base-hydrogen fluoride (HF) complexes (base = O3 and SO2) in order to elucidate the structures and energetics of the complexes. The ab initio calculations were performed up to the QCISD(T)/6-311++G(d,p) level of theory. In both complexes, hydrogen-bonded structures where the hydrogen of HF orients toward one of the oxygen atoms of bases were obtained as stable forms. The calculations showed that cis and trans isomers exist in both complexes. All calculations for the SO2-HF complex indicated that the cis form is more stable in energy than the trans form. On the other hand, in O3-HF complexes, the stable structures are changed by the ab initio levels of theory used, and the energies of the cis and trans forms are close to each other. From the most sophisticated calculations (QCISD(T)/6-311++G(d,p)//QCISD/6-311+G(d) level), it was predicted that the complex formation energies for cis SO2-HF, trans SO2-HF, cis O3-HF, and trans O3-HF are 6.1, 5.7, 3.4, and 3.6 kcal/mol, respectively, indicating that the binding energy of HF to SO2 is larger than that of O3. The harmonic vibrational frequencies calculated for cis O3-HF and cis SO2-HF complexes were in good agreement with the experimental values measured by Andrews et al. Also, the calculated rotation constants for cis SO2-HF agreed with the experiment.     

The determination of the equilibrium structures of oxygen,ozone, and hydrogen peroxide using the ab initio and density functional theory methods

The geometries and energies of small oxygen containing molecules are studied by both the ab initio and density functional theory (DFT) methods. The RHF, MP2, and QCISD(T) ab initio methods, BHandH, BHandHLYP, BeckeSLYP, Becke3P86 DFT hybrid methods, BLYP, and the BP86 non-local DFT methods with the 3-21G¹, 6-31G(d,p), 6-311 + G(2d,2p) and 6-311 + + G(3df,3pd) basis sets were used for the computational study. The obtained results from the different methods were compared to the experimental values. The suitability of the DFT methods for reproducing experimental data were discussed.     

Study of the formamide–methanol dimer with ab initio and density functional theory methods

For the first time, the structures and energies for the hydrogen bonding of a 1:1 complex formed between formamide and methanol molecules have been computed with various pure and hybrid density functional theory (DFT) and ab initio methods at varied basis set levels from 6‐31g to 6‐31+g(d,p). Five reasonable geometries on the potential energy surface of methanol and formamide system are considered and their relative stability is discussed. The infrared (IR) spectrum frequencies, IR intensities, and vibrational frequency shifts are reported. From the systematic studies, it is found that all the DFT methods selected here correctly compute the dimerization energies and geometries, with the B3P86 method predicting the hydrogen bond lengths relatively shorter and BPW91 yielding the interaction energies relatively lower. Finally, the solvent effects on the geometries of the formamide–methanol complexes have also been investigated using self‐consistent reaction field (SCRF) calculations with five different DFT methods at the 6‐31+g(d,p) basis set level. The results indicate that the polarity of the solvent has played an important role on the structures and relative stabilities of different isomers. Moreover, the basis set superposition error correction is critical to the interaction energies in the polar solvents. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Theoretical Study on the Spectroscopic Properties of CO3.−.nH2O Clusters: Extrapolation to Bulk

Vertical detachment energies (VDE) and UV/Vis absorption spectra of hydrated carbonate radical anion clusters, CO₃^.?.n H₂O (n=1–8), are determined by means of ab initio electronic structure theory. The VDE values of the hydrated clusters are calculated with second‐order Moller–Plesset perturbation (MP2) and coupled cluster theory using the 6‐311++G(d,p) set of basis functions. The bulk VDE value of an aqueous carbonate radical anion solution is predicted to be 10.6 eV from the calculated weighted average VDE values of the CO₃^.?.n H₂O clusters. UV/Vis absorption spectra of the hydrated clusters are calculated by means of time‐dependent density functional theory using the Becke three‐parameter nonlocal exchange and the Lee–Yang–Parr nonlocal correlation functional (B3LYP). The simulated UV/Vis spectrum of the CO₃^.?.8 H₂O cluster is in excellent agreement with the reported experimental spectrum for CO₃^.? (aq), obtained based on pulse radiolysis experiments.     

Vibrational spectroscopic study of pyrrole and its deuterated derivatives: comparison of the quality of the applicability of the DFT/Becke3P86 and the DFT/Becke3LYP functionals

The aim of our work is a reassignment of the vibrational spectra of pyrrole based on high-level quantum-chemical calculations. This is a continuation of our earlier works on one side of pyrrolydine and N-methylpyrrolidine and on the other side of five-membered ring parent molecules with two-ring nitrogens.

Infrared and Raman spectra are reported for the isotopic species of pyrrole, 1-deuteropyrrole and pentadeuteropyrrole.

The molecular structure and the harmonic force field were calculated applying the ab initio density functional theory (DFT) level with both the Becke3P86 and the Becke3LYP functionals with the 6-311G(d,p) basis set. The force fields were fitted to the experimental fundamentals using in both cases eight scale factors. Though some scaled frequencies show larger deviations from the experimental ones, the percentage deviations of the calculated frequencies from the experimental ones are less than 1.0% for pyrrole and less than 1.2% for the deuterated derivatives in the case of both applied functionals.     

Density functional theory and molecular clusters

Density functional theory (DFT) methods, including nonlocal density gradient terms in the exchange and correlation energy functionals, were applied to various types of molecular clusters: H-bonded, ionic, electrostatic, and London. Reliable results on the structure and stabilization energy were obtained for the first two types of cluster as long as Becke3LYP and Becke3P86 functionals and basis sets of at least DZ + P quality were used. DFT methods with currently available functionals failed completely, however, for London-type clusters, for which no minimum was found on the potential energy surfaces. DFT interaction energy exhibits the same basis set extension dependence as the Hartree-Fock (HF) interaction energy. Therefore, the Boys-Bernardi function counterpoise procedure should be employed for elimination of the DFT basis set superposition error. © 1995 John Wiley & Sons, Inc.     

A comparative study of DFT and traditional ab initio methodologies on the OsO4 molecule

The performance of different conventional ab initio methodologies and density functional procedures is compared through its application to the theoretical calculation of the bond distance and harmonic vibrational frequencies of the OsO₄ molecule. The problem of the basis set is first considered, with up to nine different basis sets being tested in calculations using the hybrid Becke3LYP density functional, and the most appropriate basis set is used in the comparison of Hartree–Fock, post‐Hartree–Fock, and density functional methods. The post‐Hartree–Fock methods analyzed are MP2, CISD, and CCSD(T), and the density functionals tested are SVWN, BLYP, BPW91, and Becke3LYP. The results show that for this particular system density functional methods perform better than do HF‐based methods with the exception of CCSD(T), which gives the best overall results. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 544–551, 2000     

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-chemical Study on the Structures and Properties of Uracil-BX3 （X = F, Cl） Complexes

1 INTRODUCTION The intermolecular interaction of bases in DNA or RNA is of immense interest and significance to che- mists and biologists alike. The interactions of these bases with metal cations, solvent molecules and other small molecules or ions would affect the struc- ture and biological properties or recognition process,which has been investigated widely[1～8]. Boron contained compounds are electron deficient com- pounds and have been extensively used as catalysts in chemical react…     

S2O2分子稳定结构的密度泛函理论研究

采用密度泛函理论中的B3LYP方法,在6-311+G(3df)和Aug-cc-pVTZ水平上,研究了单态S2O2分子的各种可能的异构体及其相对稳定性。结果表明,具有C2v对称性的三角平面分叉异构体的热力学稳定性要高于目前实验上唯一被发现的具有C2v对称性的cis-OSSO异构体,同时trans-OSSO的稳定性与顺式异构体十分接近,这2种异构体应该可以在实验上被观察到。同时本文还讨论了3个最稳定构型的前线分子轨道和链型OSSO的内扭转势能。     

Theoretical Study of Haloacetonitrile Anions: CH2XCN- (X=F,Cl)

Haloacetonitrile anions CH2XCN- (X=F, Cl) were studied by HF-SCF, Becke3-LYP, and MP2 methods together with the Dunning's basis set aug-cc-PVTZ. The vertical electron attachments to the neutral are endothermic. The geometrically optimized CH2FCN- is mainly a valence-bounded anion and CH2FCN-→CH2CN+F- is a nonadiabatic dissociation. This theoretical study in good agreement with the experimental results shows that the Becke3-LYP method is reasonable in describing the electronic structures of anions and dissociative attachment dynamics, while significant differences between MP2 and Becke3-LYP results are shown for the dissociation potential curves of CH2ClCN-→CH2CN+Cl-.2,70-(Ethylene)-bis-8-hydroxyquinoline was optimized with DFT/B3LYP and ab initio HF methods, so ionization potential and electron affinity could be determined. Absorption spectrum was calculated by ZINDO and TD-DFT. CIS method was used to calculate the S1 excited states of the compound and afterwards the emission spectrum was computed. When the solvent effect was taken into account, the computed results show encouraging agreement with known experimental data. The results of analyzing the relationship between the energies and absorption spectra indicate that the ability to transporting electrons is strengthened compared with 8-hydroxyquinoline and that absorption and emission spectra are red-shifted. The intramolecular reorganization energy of tris(2,70-(ethylene)-bis-8-hydroxyquinoline)-aluminum implies its electron transporting property is worse than tris(8-hydroxyquinoline)-aluminum. The predicted maximum emission wavelength is red-shifted compared with tris(8-hydroxyquinoline)-aluminum.