The structure and binding energies of the van der Waals complexes of Ar and N2 with phenol and its cation, studied by high level ab initio and density functional theory calculations

We have investigated, using both ab initio and density functional theory methods, the minimum energy structures and corresponding binding energies of the van der Waals complexes between phenol and argon or the nitrogen molecule, and the corresponding complexes involving the phenol cation. Structures were obtained at the MP2 level using a large basis, and the corresponding energies were corrected for basis set superposition error (BSSE), higher order electron correlation effects, and for basis set size. The structures of the global minima were further refined for the effects of BSSE and the corresponding binding energies were evaluated. For each neutral species, we find only a single true minimum, pi bonded for argon and OH bonded for nitrogen. For both cationic species, we find that the OH-bonded complex is preferred over other minima which we have identified as having Ar or N(2) between exogeneous atoms. The ab initio calculations are generally in excellent agreement with experimental binding energies and rotational constants. We find that the B3LYP functional is particularly poor at describing these complexes, while a density functional theory (DFT) method with an empirical correction for dispersive interactions (DFT-D) is very successful, as are some of the new functionals proposed by Zhao and Truhlar [J. Phys. Chem. A 109, 5656 (2005); J. Chem. Theory Comput. 2, 1009 (2006); Phys. Chem. Chem. Phys. 7, 2701 (2005); J. Phys. Chem. A 108, 6908 (2004)]. Both the ab initio and DFT-D methods accurately predict the intermolecular vibrational modes.     

Nucleic acid reactivity: Challenges for next‐generation semiempirical quantum models

Semiempirical quantum models are routinely used to study mechanisms of RNA catalysis and phosphoryl transfer reactions using combined quantum mechanical (QM)/molecular mechanical methods. Herein, we provide a broad assessment of the performance of existing semiempirical quantum models to describe nucleic acid structure and reactivity to quantify their limitations and guide the development of next‐generation quantum models with improved accuracy. Neglect of diatomic differential overlap and self‐consistent density‐functional tight‐binding semiempirical models are evaluated against high‐level QM benchmark calculations for seven biologically important datasets. The datasets include: proton affinities, polarizabilities, nucleobase dimer interactions, dimethyl phosphate anion, nucleoside sugar and glycosidic torsion conformations, and RNA phosphoryl transfer model reactions. As an additional baseline, comparisons are made with several commonly used density‐functional models, including M062X and B3LYP (in some cases with dispersion corrections). The results show that, among the semiempirical models examined, the AM1/d‐PhoT model is the most robust at predicting proton affinities. AM1/d‐PhoT and DFTB3‐3ob/OPhyd reproduce the MP2 potential energy surfaces of 6 associative RNA phosphoryl transfer model reactions reasonably well. Further, a recently developed linear‐scaling “modified divide‐and‐conquer” model exhibits the most accurate results for binding energies of both hydrogen bonded and stacked nucleobase dimers. The semiempirical models considered here are shown to underestimate the isotropic polarizabilities of neutral molecules by approximately 30%. The semiempirical models also fail to adequately describe torsion profiles for the dimethyl phosphate anion, the nucleoside sugar ring puckers, and the rotations about the nucleoside glycosidic bond. The modeling of pentavalent phosphorus, particularly with thio substitutions often used experimentally as mechanistic probes, was problematic for all of the models considered. Analysis of the strengths and weakness of the models suggests that the creation of robust next‐generation models should emphasize the improvement of relative conformational energies and barriers, and nonbonded interactions. © 2015 Wiley Periodicals, Inc.     

Theoretical study of salicylaldehyde conformal isomers and their intramolecular oxygen and hydrogen relations

To verify the semiempirical‐type localized hydrogen bonding analysis methods introduced by us several years ago, the intramolecular oxygen and hydrogen relations within salicylaldehyde are selected as the major topic in this theoretical study. The B3LYP/6‐31G** density functional method is chosen for both the full‐optimization and frequency‐type calculations. Four ortho‐type planar conformal isomers are proven to be local minima, and four internal rotation transition states are found by QST3‐type calculation. The special interpretations of  CHO and  OH characteristic frequencies, energy barriers, and thermal chemical results are discussed. In the semiempirical scheme, both local hydrogen bonding population analysis and localized hydrogen bond energy breaking procedures are applied to five pairs of related oxygen and hydrogen atoms in each isomer. The explanations for the strong or weak hydrogen bonds and intra‐CHO repulsion relationships are discussed. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 74: 395–404, 1999     

Glycine Dimers: Structure,Stability, and Medium Effects

We report a study on different ionization states and conformations of the bimolecular (Gly)₂ system by means of quantum mechanical calculations. Optimized geometries for energy minima of the glycine dimer, as well as relative energies and free energies were computed as a function of the medium: gas phase, nonpolar polarizable solvent, and aqueous solution. The polarizable continuum model was employed to account for solvation effects. Energy calculations were done using the MP2/aug‐cc‐pVTZ and B3LYP/6‐311+G(2df,2p) methods on B3LYP/6‐31+G(d,p) optimized structures (some single‐point energy calculations were also done using the B3PW91 and PBE1KCIS methods). Ionized forms of the glycine dimer (either zwitterion–zwitterion or neutral–zwitterion) are predicted to exist in all media, in contrast to amino acid monomers. In aqueous solution, dimerization is an exergonic process (?4 kcal mol^?1). Thus, according to our results, zwitterion–zwitterion Gly dimers might be abundant in supersaturated glycine aqueous solutions, a fact that has been connected with the structure of α‐glycine crystals but that remains controversial in the literature. Another noticeable result is that zwitterion–zwitterion interactions are substantially underestimated when computed using methods based on density functional theory. For comparison, some calculations for the dimer of the simplest chiral amino acid alanine were done as well and differences to the glycine dimer are discussed.     

Performance of semiempirical methods in fullerene chemistry: relative energies and nucleus-independent chemical shifts

Semiempirical MNDO, AM1, and PM3 calculations are reported for 153 fullerene isomers in an attempt to assess the reliability of these methods through comparisons with ab initio and density functional results. B3LYP/6-31G^* relative energies are generally reproduced quite well by these calculations. Qualitative trends in ab initio nucleus-independent chemical shifts at the cage centers are captured by the semiempirical GIAO-MNDO approach while underestimating their absolute values. The agreement between the semiempirical results and the ab initio or density functional reference data is generally better for the larger fullerenes (C₆₀–C₁₀₂) than for the smaller ones (C₂₀–C₅₀). These systematic comparisons clarify the accuracy that may be expected from semiempirical computations in fullerene chemistry.     

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     

Prediction of activation energies for aromatic oxidation by cytochrome P450

We have estimated the activation energy for aromatic oxidation by compound I in cytochrome P450 for a diverse set of 17 substrates using state-of-the-art density functional theory (B3LYP) with large basis sets. The activation energies vary from 60 to 87 kJ/mol. We then test if these results can be reproduced by computationally less demanding methods. The best methods (a B3LYP calculation of the activation energy of a methoxy-radical model or a partial least-squares model of the semiempirical AM1 bond dissociation energies and spin densities of the tetrahedral intermediate for both a hydroxyl-cation and a hydroxyl-radical model) give correlations with r(2) of 0.8 and mean absolute deviations of 3 kJ/mol. Finally, we apply these simpler methods on several sets of reactions for which experimental data are available and show that we can predict the reactive sites by combining calculations of the activation energies with the solvent-accessible surface area of each site.     

An ab initio study of water clusters in gas phase and bulk aqueous media: (H2O)n,n=1–12

Geometries of several clusters of water molecules including single minimum energy structures of n‐mers (n=1–5), several hexamers and two structures of each of heptamer to decamer derived from hexamer cage and hexamer prism were optimized. One structural form of each of 11‐mer and 12‐mer were also studied. The geometry optimization calculations were performed at the RHF/6‐311G* level for all the cases and at the MP2/6‐311++G** level for some selected cases. The optimized cluster geometries were used to calculate total energies of the clusters in gas phase employing the B3LYP density functional method and the 6‐311G* basis set. Frequency analysis was carried out in all the cases to ensure that the optimized geometries corresponded to total energy minima. Zero‐point and thermal free energy corrections were applied for comparison of energies of certain hexamers. The optimized cluster geometries were used to solvate the clusters in bulk water using the polarized continuum model (PCM) of the self‐consistent reaction field (SCRF) theory, the 6‐311G* basis set, and the B3LYP density functional method. For the cases for which MP2/6‐311++G** geometry optimization was performed, solvation calculations in water were also carried out using the B3LYP density functional method, the 6‐311++G** basis set, and the PCM model of SCRF theory, besides the corresponding gas‐phase calculations. It is found that the cage form of water hexamer cluster is most stable in gas phase among the different hexamers, which is in agreement with the earlier theoretical and experimental results. Further, use of a newly defined relative population index (RPI) in terms of successive total energy differences per water molecule for different cluster sizes suggests that stabilities of trimers, hexamers, and nonamers in gas phase and those of hexamers and nonamers in bulk water would be favored while those of pentamer and decamer in both the phases would be relatively disfavored. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 90–104, 2001     

Theoretical studies of dissociative phosphoryl transfer in interconversion of phosphoenolpyruvate to phosphonopyruvate: solvent effects, thio effects, and implications for enzymatic reactions

The conversion of phosphoenolpyruvate (PEP) to phosphonopyruvate (P-pyr) is catalyzed by PEP mutase via a dissociative mechanism. In this work, we investigate the uncatalyzed reaction using ab initio methods, density functional theory, and the semiempirical MNDO/d method. Comparisons of geometries and relative energies of stationary points (minima and transition states) with density functional results indicate that the semiempirical method is reasonably accurate. Solvent effects are examined using implicit solvent models, including the recently extended smooth conductor-like screening model. Due to the large negative charge carried by the system, solvation is found to drastically alter the location and energy of stationary points along the dissociative reaction pathways. The influence of substituting a nonbridging phosphoryl oxygen by sulfur (thio effects) was also investigated. Implications of these results for the enzymatic reaction are discussed.     

Comparison of aromatic NH···π, OH···π, and CH···π interactions of alanine using MP2, CCSD, and DFT methods

A comparison of the performance of various density functional methods including long‐range corrected and dispersion corrected methods [MPW1PW91, B3LYP, B3PW91, B97‐D, B1B95, MPWB1K, M06‐2X, SVWN5, ωB97XD, long‐range correction (LC)‐ωPBE, and CAM‐B3LYP using 6‐31+G(d,p) basis set] in the study of CH···π, OH···π, and NH···π interactions were done using weak complexes of neutral (A) and cationic (A⁺) forms of alanine with benzene by taking the Møller–Plesset (MP2)/6‐31+G(d,p) results as the reference. Further, the binding energies of the neutral alanine–benzene complexes were assessed at coupled cluster (CCSD)/6‐31G(d,p) method. Analysis of the molecular geometries and interaction energies at density functional theory (DFT), MP2, CCSD methods and CCSD(T) single point level reveal that MP2 is the best overall performer for noncovalent interactions giving accuracy close to CCSD method. MPWB1K fared better in interaction energy calculations than other DFT methods. In the case of M06‐2X, SVWN5, and the dispersion corrected B97‐D, the interaction energies are significantly overrated for neutral systems compared to other methods. However, for cationic systems, B97‐D yields structures and interaction energies similar to MP2 and MPWB1K methods. Among the long‐range corrected methods, LC‐ωPBE and CAM‐B3LYP methods show close agreement with MP2 values while ωB97XD energies are notably higher than MP2 values. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010     

PDDG/PM3 and PDDG/MNDO: improved semiempirical methods

Two new semiempirical methods employing a Pairwise Distance Directed Gaussian modification have been developed: PDDG/PM3 and PDDG/MNDO; they are easily implemented in existing software, and yield heats of formation for compounds containing C, H, N, and O atoms with significantly improved accuracy over the standard NDDO schemes, PM5, PM3, AM1, and MNDO. The PDDG/PM3 results for heats of formation also show substantial improvement over density functional theory with large basis sets. The PDDG modifications consist of a single function, which is added to the existing pairwise core repulsion functions within PM3 and MNDO, a reparameterized semiempirical parameter set, and modified computation of the energy of formation of a gaseous atom. The PDDG addition introduces functional group information via pairwise atomic interactions using only atom-based parameters. For 622 diverse molecules containing C, H, N, and O atoms, mean absolute errors in calculated heats of formation are reduced from 4.4 to 3.2 kcal/mol and from 8.4 to 5.2 kcal/mol using the PDDG modified versions of PM3 and MNDO over the standard versions, respectively. Several specific problems are overcome, including the relative stability of hydrocarbon isomers, and energetics of small rings and molecules containing multiple heteroatoms. The internal consistency of PDDG energies is also significantly improved, enabling more reliable analysis of isomerization energies and trends across series of molecules; PDDG isomerization energies show significant improvement over B3LYP/6-31G* results. Comparison of heats of formation, ionization potentials, dipole moments, isomer, and conformer energetics, intermolecular interaction energies, activation energies, and molecular geometries from the PDDG techniques is made to experimental data and values from other semiempirical and ab initio methods.     

Quantum study on the electrocyclic reactions of CH2CHCHCHX (XH,F, Cl)

The microscopic mechanisms of the electrocyclic reactions for cis‐1,3‐butadiene and its monofluoro‐, monochloroderivatives have been studied by density functional theory (DFT), using the B3LYP method and 6‐311++G** basis sets. We optimized the geometric configurations of reactants, transition states, and products; verified all the probable transition states through vibrational analysis; and calculated the relative single‐point energies at the QCISD(T)/6‐311++G**//B3LYP/6‐311++G**. The results show that the monofluoro‐, monochloroderivatives of cis‐1,3‐butadiene both have two conformers; the reactant favors the electrocyclic reaction when one outboard hydrogen atom of the CH₂ groups is substituted by the fluorine or chlorine atom. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

A systematic quantum chemical investigation of the C-H bond activation in methane by gas phase vanadium oxide cation VO+

The interaction between a methane molecule and the VO(+) cation in the gas phase has been investigated by means of single reference density functional (B3LYP) and wave function-based multireference (MR) correlation calculations. For the latter, an extrapolation technique is used to evaluate correlation energies at the basis set limit. A comprehensive picture for the C-H activation features a variety of molecular structures corresponding to both minima and transition states. Possible reaction paths are discussed, also taking into account change of the spin multiplicity. Activation of the methane molecule by VO(+) is always an endothermic process. Competing reaction paths might be expected. An evaluation of miscellaneous computational methods is performed using calculated energy differences for various molecular structures. Results obtained from the MR calculations exhibit no systematic convergence with increasing size of the active space used, and for two largest active spaces relative energies still differ by up to 25 kJ/mol. Simple mean difference between the B3LYP results and the best MR values is -50 +/- 19 kJ/mol.     

Probing furanose ring conformation by gas-phase computational methods: energy profile and structural parameters in methyl beta-D-arabinofuranoside as a function of ring conformation

The potential energy surface of methyl beta-D-arabinofuranoside (3) has been studied by ab initio molecular orbital (HF/6-31G) and density functional theory (B3LYP/6-31G) calculations via minimization of the 10 possible envelope conformers. The partial potential energy surface identified that the global minimum and lowest energy northern conformer was E(2). In the HF calculations, (2)E was the most stable southern conformer, while the density functional theory methods identified (4)E as the local minimum in this hemisphere. Additional calculations at higher levels of theory showed that the B3LYP-derived energies of many of the envelope conformers of 3 are dependent upon the basis set used. It has also been demonstrated that B3LYP/6-31+G//B3LYP/6-31G single point energies are essentially the same as those obtained from full geometry optimizations at the B3LYP/6-31+G level. The northern and southern minima of the B3LYP/6-31+G surface are, respectively, the E(2) and (2)E conformers. The B3LYP/6-31G geometries were used to study the relationship between ring conformation and various structural parameters including bond angles, dihedral angles, bond lengths, and interatomic distances.     

A semiempirical long-range corrected exchange correlation functional including a short-range Gaussian attenuation (LCgau-B97)

We applied an improved long‐range correction scheme including a short‐range Gaussian attenuation (LCgau) to the Becke97 (B97) exchange correlation functional. In the optimization of LCgau‐B97 functional, the linear parameters are determined by least squares fitting. Optimizing μ parameter (0.2) that controls long‐range portion of Hartree‐Fock (HF) exchange to excitation energies of large molecules (Chai and Head‐Gordon, J Chem Phys 2008, 128, 084106) and additional short‐range Gaussian parameters (a = 0.15 and k = 0.9) that controls HF exchange inclusion ranging from short‐range to mid‐range (0.5–3 Å) to ground state properties achieved high performances of LCgau‐B97 simultaneously on both ground state and excited state properties, which is better than other tested semiempirical density functional theory (DFT) functionals, such as ωB97, ωB97X, BMK, and M0x‐family. We also found that while a small μ value (～0.2) in LC‐DFT is appropriate to the local excitation and intramolecular charge‐transfer excitation energies, a larger μ value (0.42) is desirable in the Rydberg excitation‐energy calculations. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

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     

Energetics of reactions involving radical species in solution: Calculation of relative electrode potentials for nitroimidazoles using density functional and continuum methods

Nonlocal density functional calculations and a semiempirical modified Born method for computing free energies of hydration were used to calculate the electrode potentials for a series of nitroimidazoles to a mean accuracy of about 80 mV. The density functional calculations used the nonlocal Becke '88 functional for exchange and either the nonlocal Lee-Yang-Parr or the local Vosko-Wilk-Nusair functionals for correlation and were performed at the HF/3-21G geometry. The most suitable geometry for these calculations was determined from a survey of various semiempirical, Hartree-Fock (HF) and density functional methods, with a variety of basis sets. The HF/3-21G method was found to yield a very favorable compromise between speed and accuracy in the determination of the geometry of 2-nitroimidazole, but the small basis set density functional calculations performed very badly. Density functional atom-optimized basis sets were found to give better overall results than traditional Pople-type basis sets. The free energy of hydration calculations employed the AM1 SM2 method. Both the gas-phase energies and the free energies of hydration made a significant contribution to the computed electrode potential. Indeed, an inverse relationship was found between the gas-phase electron affinity and the difference in free energy of hydration between the neutral nitroimidazole and its radical anion. The protocol established here may be useful for investigating novel bioreductive agents. © 1996 John Wiley & Sons, Inc.     

Activation Energy and Transition State Determination of the Olefin Insertion Process of Metallocene Catalysts Using a Semiempirical Molecular Orbital Calculation

The transition state of the olefin insertion process of metallocene catalysts can be determined by adopting the semiempirical PM3 model. In computational chemistry, the computational methods most employed are the ab initio method and density functional theory, which are very time consuming. The semiempirical molecular orbital method requires much less computational resources than the above methods. However, the accuracy and reliability of the semiempirical molecular orbital method remains to be determined. The PM3 model is the most recently developed the semiempirical molecular orbital method and can also be applied to transition metal calculations. This study is intended to investigate the reliability of computational results determined using semiempirical PM3 model on metallocene catalysts through comparison with published results on the density functional theory (DFT). The saddle point finding procedure is adopted to find the transition state of the ethylene insertion process of metallocene catalysts. Results on the geometry and energy trends of the ethylene insertion process of metallocene catalysts determined using the PM3 model are in good agreement with the DFT results. In addition, the saddle point of the potential energy surface of ethylene insertion is verified in accordance with the eigenvalue of the vibrational frequency spectrum. Correct eigenvalues indicate that the correct saddle point of the potential energy surface of ethylene insertion has been successfully located. Hence, the eigenvalue of the vibrational frequency spectrum is a valuable reference in terms of saddle point justification. Computational results and vibrational frequency spectrum analysis demonstrate that the PM3 model can be used to locate the correct saddle point of the potential energy surface. The results obtained using the PM3 model confirm that the eigenvalue of the transition state lies nearly on the vibrational frequency spectrum. The eigenvalues are also analyzed, providing a valuable reference for further studies of the transition state of olefin insertion of metallocene catalysts. The activation energies for the olefin insertion reaction are also studied for evaluation of the catalyst.