Solvent effects on the S(N)2 reaction: Application of the density functional theory-based effective fragment potential method

The performance of the density functional theory (DFT)-based effective fragment potential (EFP) method is assessed using the S(N)2 reaction: Cl- + nH2O + CH3Br = CH3Cl + Br- + nH2O. The effect of the systematic addition of water molecules on the structures and relative energies of all species involved in the reaction has been studied. The EFP1 method is compared with second-order perturbation theory (MP2) and DFT results for n = 1, 2, and 3, and EFP1 results are also presented for four water molecules. The incremental hydration effects on the barrier height are the same for all methods. However, only full MP2 or MP2 with EFP1 solvent molecules are able to provide an accurate treatment of the transition state (TS) and hence the central barriers. Full DFT and DFT with EFP1 solvent molecules both predict central barriers that are too small. The results illustrate that the EFP1-based DFT method gives reliable results when combined with an accurate quantum mechanical (QM) method, so it may be used as an efficient alternative to fully QM methods in the treatment of larger microsolvated systems.     

A comparison of quantum chemical models for calculating NMR shielding parameters in peptides: mixed basis set and ONIOM methods combined with a complete basis set extrapolation

This article compares several quantum mechanical approaches to the computation of chemical shielding tensors in peptide fragments. First, we describe the effects of basis set quality up to the complete basis set (CBS) limit and level of theory (HF, MP2, and DFT) for four different atoms in trans N-methylacetamide. For both isotropic shielding and shielding anisotropy, the MP2 results in the CBS limit show the best agreement with experiment. The HF values show quite a different tendency to MP2, and even in the CBS limit they are far from experiment for not only the isotropic shielding of carbonyl carbon but also most shielding anisotropies. In most cases, the DFT values differ systematically from MP2, and small basis-set (double- or triple-zeta) results are often fortuitously in better agreement with the experiment than the CBS ones. Second, we compare the mixed basis set and ONIOM methods, combined with CBS extrapolation, for chemical shielding calculations at a DFT level using various model peptides. From the results, it is shown that the mixed basis set method provides better results than ONIOM, compared to CBS calculations using the nonpartitioned full systems. The information studied here will be useful in guiding the selection of proper quantum chemical models, which are in a tradeoff between accuracy and cost, for shielding studies of peptides and proteins.     

The dispersion interaction between quantum mechanics and effective fragment potential molecules

A method for calculating the dispersion energy between molecules modeled with the general effective fragment potential (EFP2) method and those modeled using a full quantum mechanics (QM) method, e.g., Hartree-Fock (HF) or second-order perturbation theory, is presented. C(6) dispersion coefficients are calculated for pairs of orbitals using dynamic polarizabilities from the EFP2 portion, and dipole integrals and orbital energies from the QM portion of the system. Dividing by the sixth power of the distance between localized molecular orbital centroids yields the first term in the commonly employed London series expansion. A C(8) term is estimated from the C(6) term to achieve closer agreement with symmetry adapted perturbation theory values. Two damping functions for the dispersion energy are evaluated. By using terms that are already computed during an ordinary HF or EFP2 calculation, the new method enables accurate and extremely rapid evaluation of the dispersion interaction between EFP2 and QM molecules.     

Infrared and Raman spectra and theoretical study of methyl trifluoromethyl sulfone, CF3SO2CH3

The infrared and Raman spectra were obtained for liquid CF3SO2CH3, as well as the infrared spectrum of the gaseous substance. The molecular geometry was optimized by means of the Hartree-Fock (HF), second order electron correlation (MP2) and density functional theory (DFT) procedures of quantum chemistry, resulting in a structure with Cs symmetry. The wavenumbers corresponding to the normal modes of vibration were calculated using the DFT (B3LYP/6-31G**) approximation and their agreement with the measured values improved after scaling of the associated force field. An assignment of bands is proposed on the basis of such calculations and the comparison with related molecules.     

De novo prediction of the ground state structure of transition metal complexes using semiempirical and ab initio quantum mechanics. Coordination isomerism

The ground state coordination isomers for 30 different trigonal bipyramidal transition metal complexes have been predicted using different levels of quantum mechanics: semiempirical (PM3(tm)), ab initio (MP2//HF), pure (BPW91) and hybrid (B3PW91) density functional theory (DFT) methods. For species where these methods failed to reproduce crystallographic data, hybrid quantum mechanics/molecular mechanics (QM/MM) methods were used to study more exact experimental models. Literature deficiencies regarding ground state multiplicity of these species were supplemented by spin predictions using previously tested PM3(tm) methods. Geometry optimization calculations were performed for each possible coordination isomer. The predicted ground state minima provided by the different methods are compared to each other and with crystallographic data. Pure DFT functionals outperformed hybrid functionals and MP2//HF. The very rapid PM3(tm) parameterization method provided accurate predictions in comparison to other levels of theory. An integrated MM/PM3(tm)/DFT de novo scheme accurately reproduced crystallographic data for species where the individual methods failed.     

A combined effective fragment potential-fragment molecular orbital method. II. Analytic gradient and application to the geometry optimization of solvated tetraglycine and chignolin

The gradient for the fragment molecular orbital (FMO) method interfaced with effective fragment potentials (EFP), denoted by FMO∕EFP, was developed and applied to polypeptides solvated in water. The structures of neutral and zwitterionic tetraglycine immersed in water layers of 2.0, 2.5, 3.0, 3.5, 4.0, and 4.5 A? are investigated by performing FMO∕EFP geometry optimizations at the RHF∕cc-pVDZ level of theory for the solutes. The geometries optimized with FMO-RHF∕EFP are compared to those from the conventional RHF∕EFP and are found to be in very close agreement. Using the optimized geometries, the stability of the hydrated zwitterionic and neutral structures is discussed structurally and in terms of energetics at the second-order M?ller-Plesset theory (MP2)∕cc-pVDZ level. To demonstrate the potential of the method for proteins, the geometry of hydrated chignolin (protein data bank ID: 1UAO) was optimized, and the importance of the inclusion of water was examined by comparing the solvated and gas phase structures of chignolin with the experimental NMR structure.     

Vibrational frequencies of hydrazoic acid and methyl azide: density functional theory study

Harmonic vibrational frequencies of HN3 and CH3N3 molecules and their several isotopomers are calculated using HF, MP2 and five popular density functional theory (DFT) methods. On the basis of the comparison between calculated and experimental results, assignments of fundamental vibrational modes arc examined. HF and MP2 results are in bad agreement with experimental values. Of the five DFT methods, BLYP reproduces the observed fundamental frequencies the most satisfactorily. Two hybrid DFT methods are found to yield frequencies generally higher than the observed fundamental frequencies. The results indicate that BLYP calculation is a very promising approach for understanding the observed spectral features.     

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     

An assessment of density functional methods for studying molecular adsorption in cluster models of zeolites

The performance of six density functional theory (DFT) methods has been tested for a zeolite cluster containing three tetrahedral atoms (3T) and the complexes it forms with water and methanol molecules. The DFT methods (BLYP, BP86, BPW91, B3LYP, B3P86, B3PW91) give results in good agreement with second-order perturbation theory (MP2). The results in this paper provide evidence of the suitability of DFT methods for studying hydrogen-bonded adsorption complexes in zeolites. Generally, the hybrid DFT methods are in closer agreement with experiment and MP2 than the pure DFT methods for geometrical parameters. The only exception is the Z⁻ geometry, perhaps due to its anionic character. All DFT methods give results in good overall agreement with MP2 for intramolecular geometrical parameters of the adsorption complexes, intramolecular vibrational frequencies, and adsorption energies. The B3LYP method gives intermolecular geometries and intermolecular vibrational frequencies which are closest to those obtained from the MP2 method. Thus, the B3LYP method seems to be the best choice for a density functional treatment of molecular adsorption in zeolite systems.     

Geometry, dipole moment, polarizability and first hyperpolarizability of polymethineimine: an assessment of electron correlation contributions

We have investigated the geometries as well as the longitudinal dipole moment (micro), polarizability (alpha), and first hyperpolarizability (beta) of polymethineimine oligomers using different approaches [Hartree-Fock (HF), second-order M?ller-Plesset (MP2), and hybrid density functional theory (DFT) methods (B3LYP and PBE0)] for evaluating the geometries and the nonlinear optical properties. It turns out that (i) HF and the selected DFT methods provide the incorrect sign for beta of short and medium size oligomers. (ii) The B3LYP and PBE0 electron correlation correction are too small for micro, too large for alpha, and for some oligomer lengths, they are in the wrong direction for beta. (iii) On the contrary to polyacetylene, the hybrid-DFT geometries are in poor agreement with MP2 geometries; the former showing much smaller bond length alternations.     

A gauge-independent zeroth-order regular approximation to the exact relativistic Hamiltonian--formulation and applications

A simple modification of the zeroth-order regular approximation (ZORA) in relativistic theory is suggested to suppress its erroneous gauge dependence to a high level of approximation. The method, coined gauge-independent ZORA (ZORA-GI), can be easily installed in any existing nonrelativistic quantum chemical package by programming simple one-electron matrix elements for the quasirelativistic Hamiltonian. Results of benchmark calculations obtained with ZORA-GI at the Hartree-Fock (HF) and second-order Moller-Plesset perturbation theory (MP2) level for dihalogens X(2) (X=F,Cl,Br,I,At) are in good agreement with the results of four-component relativistic calculations (HF level) and experimental data (MP2 level). ZORA-GI calculations based on MP2 or coupled-cluster theory with single and double perturbations and a perturbative inclusion of triple excitations [CCSD(T)] lead to accurate atomization energies and molecular geometries for the tetroxides of group VIII elements. With ZORA-GI/CCSD(T), an improved estimate for the atomization energy of hassium (Z=108) tetroxide is obtained.     

Reassessment of methyl rotation barriers and conformations by correlated quantum chemistry methods

Internal rotations of the methyl group in ortho‐substituted and 2,6‐disubstituted toluenes in their ground state have been investigated by means of various ab initio quantum chemistry methods. Computed barriers at the Hartree‐Fock (HF) level using medium sized basis sets agreed reasonably with experimental results in the case of the studied ortho‐substituted toluenes. However, this agreement worsens when using very large basis sets. Furthermore, the determination of the conformation and barriers of more weakly hindered methyl groups, that is, for 2,6‐dihalogenotoluenes or toluene itself, necessitates high level correlated computations, because of a possible failure of HF calculations in this case. Density functional theory (DFT) techniques required, in several cases, much more extended basis sets than the post‐HF Møller‐Plesset perturbation (MP2, MP4) ones, to insure the convergence of the computed barriers. Non‐negligible variations of the computed barriers when using different DFT functionals are observed for some systems. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 2093–2100, 2003     

Application of a semi-empirical dispersion correction for modeling water clusters

The Grimme-D3 semi-empirical dispersion energy correction has been implemented for the original effective fragment potential for water (EFP1), and for systems that contain water molecules described by both correlated ab initio quantum mechanical (QM) molecules and EFP1. Binding energies obtained with these EFP1-D and QM/EFP1-D methods were tested using 27 benchmark species, including neutral, protonated, deprotonated, and auto-ionized water clusters and nine solute–water binary complexes. The EFP1-D and QM/EFP1-D binding energies are compared with those obtained using fully QM methods: second-order perturbation theory, and coupled cluster theory, CCSD(T), at the complete basis set (CBS) limit. The results show that the EFP1-D and QM/EFP1-D binding energies are in good agreement with CCSD(T)/CBS binding energies with a mean absolute error of 5.9 kcal/mol for water clusters and 0.8 kcal/mol for solute–water binary complexes. © 2018 Wiley Periodicals, Inc.     

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

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

Interfacing Q-Chem and CHARMM to perform QM/MM reaction path calculations

A hybrid quantum mechanical/molecular mechanical (QM/MM) potential energy function with Hartree-Fock, density functional theory (DFT), and post-HF (RIMP2, MP2, CCSD) capability has been implemented in the CHARMM and Q-Chem software packages. In addition, we have modified CHARMM and Q-Chem to take advantage of the newly introduced replica path and the nudged elastic band methods, which are powerful techniques for studying reaction pathways in a highly parallel (i.e., parallel/parallel) fashion, with each pathway point being distributed to a different node of a large cluster. To test our implementation, a series of systems were studied and comparisons were made to both full QM calculations and previous QM/MM studies and experiments. For instance, the differences between HF, DFT, MP2, and CCSD QM/MM calculations of H2O...H2O, H2O...Na+, and H2O...Cl- complexes have been explored. Furthermore, the recently implemented polarizable Drude water model was used to make comparisons to the popular TIP3P and TIP4P water models for doing QM/MM calculations. We have also computed the energetic profile of the chorismate mutase catalyzed Claisen rearrangement at various QM/MM levels of theory and have compared the results with previous studies. Our best estimate for the activation energy is 8.20 kcal/mol and for the reaction energy is -23.1 kcal/mol, both calculated at the MP2/6-31+G(d)//MP2/6-31+G(d)/C22 level of theory.     

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.     

Protonation of glycine and alanine: proton affinities, intrinsic basicities and proton transfer path

Proton affinities and intrinsic basicities for nitrogen and oxygen protonation in the gas phase of the amino acids glycine and alanine were calculated using density functional theory (DFT) and ab initio methods at different levels of theory from Hartree-Fock (HF) to G2 approximations. All methods gave good agreement for proton affinities for nitrogen protonation for both amino acids. However, dramatic differences were found between DFT, MP4//MP2, and G2 results on one hand, and MP4//HF results on the other to the calculation of structural and energetic characteristics of oxygen protonation in glycine and alanine. An investigation into the source of these differences revealed that electron correlation effects are chiefly responsible for the differences in calculated oxygen proton affinities between the various methods. It has been found that proton transfer between nitrogen and oxygen protonation sites in both amino acids occurs without a transfer path barrier when correlated methods were used to calculate the path energetics.     

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

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

Study on the conformational equilibrium of the alanine dipeptide in water solution by using the averaged solvent electrostatic potential from molecular dynamics methodology

The ASEP/MD method has been employed for studying the solvent effect on the conformational equilibrium of the alanine dipeptide in water solution. MP2 and density functional theory (DFT) levels of theory were used and results were compared. While in gas phase cyclic structures showing intramolecular hydrogen bonds were found to be the most stable, the stability order is reversed in water solution. Intermolecular interaction with the solvent causes the predominance of extended structures as the stabilizing contacts dipeptide-water are favoured. Free-energy differences in solution were calculated and PPII, α(R), and C5 conformers were identified as the most stable at MP2 level. Experimental data from Raman and IR techniques show discrepancies about the relative abundance of α(R) y C5, our results support the Raman data. The DFT level of theory agrees with MP2 in the location and stability of PPII and α(R) forms but fails in the location of C5. MP2 results suggest the possibility of finding traces of C7eq conformer in water solution, in agreement with recent experiments.     

Comparison of some representative density functional theory and wave function theory methods for the studies of amino acids

Energies of different conformers of 22 amino acid molecules and their protonated and deprotonated species were calculated by some density functional theory (DFT; SVWN, B3LYP, B3PW91, MPWB1K, BHandHLYP) and wave function theory (WFT; HF, MP2) methods with the 6-311++G(d,p) basis set to obtain the relative conformer energies, vertical electron detachment energies, deprotonation energies, and proton affinities. Taking the CCSD/6-311++G(d,p) results as the references, the performances of the tested DFT and WFT methods for amino acids with various intramolecular hydrogen bonds were determined. The BHandHLYP method was the best overall performer among the tested DFT methods, and its accuracy was even better than that of the more expensive MP2 method. The computational dependencies of the five DFT methods and the HF and MP2 methods on the basis sets were further examined with the 6-31G(d,p), 6-311++G(d,p), aug-cc-pVDZ, 6-311++G(2df,p), and aug-cc-pVTZ basis sets. The differences between the small and large basis set results have decreased quickly for the hybrid generalized gradient approximation (GGA) methods. The basis set convergence of the MP2 results has been, however, very slow. Considering both the cost and the accuracy, the BHandHLYP functional with the 6-311++G(d,p) basis set is the best choice for the amino acid systems that are rich in hydrogen bonds.