Accurate description of van der Waals complexes by density functional theory including empirical corrections

An empirical method to account for van der Waals interactions in practical calculations with the density functional theory (termed DFT-D) is tested for a wide variety of molecular complexes. As in previous schemes, the dispersive energy is described by damped interatomic potentials of the form C6R(-6). The use of pure, gradient-corrected density functionals (BLYP and PBE), together with the resolution-of-the-identity (RI) approximation for the Coulomb operator, allows very efficient computations for large systems. Opposed to previous work, extended AO basis sets of polarized TZV or QZV quality are employed, which reduces the basis set superposition error to a negligible extend. By using a global scaling factor for the atomic C6 coefficients, the functional dependence of the results could be strongly reduced. The "double counting" of correlation effects for strongly bound complexes is found to be insignificant if steep damping functions are employed. The method is applied to a total of 29 complexes of atoms and small molecules (Ne, CH4, NH3, H2O, CH3F, N2, F2, formic acid, ethene, and ethine) with each other and with benzene, to benzene, naphthalene, pyrene, and coronene dimers, the naphthalene trimer, coronene. H2O and four H-bonded and stacked DNA base pairs (AT and GC). In almost all cases, very good agreement with reliable theoretical or experimental results for binding energies and intermolecular distances is obtained. For stacked aromatic systems and the important base pairs, the DFT-D-BLYP model seems to be even superior to standard MP2 treatments that systematically overbind. The good results obtained suggest the approach as a practical tool to describe the properties of many important van der Waals systems in chemistry. Furthermore, the DFT-D data may either be used to calibrate much simpler (e.g., force-field) potentials or the optimized structures can be used as input for more accurate ab initio calculations of the interaction energies.     

Effect of the damping function in dispersion corrected density functional theory

It is shown by an extensive benchmark on molecular energy data that the mathematical form of the damping function in DFT-D methods has only a minor impact on the quality of the results. For 12 different functionals, a standard "zero-damping" formula and rational damping to finite values for small interatomic distances according to Becke and Johnson (BJ-damping) has been tested. The same (DFT-D3) scheme for the computation of the dispersion coefficients is used. The BJ-damping requires one fit parameter more for each functional (three instead of two) but has the advantage of avoiding repulsive interatomic forces at shorter distances. With BJ-damping better results for nonbonded distances and more clear effects of intramolecular dispersion in four representative molecular structures are found. For the noncovalently-bonded structures in the S22 set, both schemes lead to very similar intermolecular distances. For noncovalent interaction energies BJ-damping performs slightly better but both variants can be recommended in general. The exception to this is Hartree-Fock that can be recommended only in the BJ-variant and which is then close to the accuracy of corrected GGAs for non-covalent interactions. According to the thermodynamic benchmarks BJ-damping is more accurate especially for medium-range electron correlation problems and only small and practically insignificant double-counting effects are observed. It seems to provide a physically correct short-range behavior of correlation/dispersion even with unmodified standard functionals. In any case, the differences between the two methods are much smaller than the overall dispersion effect and often also smaller than the influence of the underlying density functional.     

Fully ab initio protein‐ligand interaction energies with dispersion corrected density functional theory

Dispersion corrected density functional theory (DFT‐D3) is used for fully ab initio protein‐ligand (PL) interaction energy calculation via molecular fractionation with conjugated caps (MFCC) and applied to PL complexes from the PDB comprising 3680, 1798, and 1060 atoms. Molecular fragments with n amino acids instead of one in the original MFCC approach are considered, thereby allowing for estimating the three‐body and higher many‐body terms. n > 1 is recommended both in terms of accuracy and efficiency of MFCC. For neutral protein side‐chains, the computed PL interaction energy is visibly independent of the fragment length n. The MFCC fractionation error is determined by comparison to a full‐system calculation for the 1060 atoms containing PL complex. For charged amino acid side‐chains, the variation of the MFCC result with n is increased. For these systems, using a continuum solvation model with a dielectricity constant typical for protein environments (? = 4) reduces both the variation with n and improves the stability of the DFT calculations considerably. The PL interaction energies for two typical complexes obtained ab initio for the first time are found to be rather large (?30 and ?54 kcal/mol). © 2012 Wiley Periodicals, Inc.     

Implementation of empirical dispersion corrections to density functional theory for periodic systems

A recently developed empirical dispersion correction (Grimme et al., J. Chem. Phys. 2010, 132, 154104) to standard density functional theory (DFT‐D3) is implemented in the plane‐wave program package VASP. The DFT‐D3 implementation is compared with an implementation of the earlier DFT‐D2 version (Grimme, J. Comput. Chem. 2004, 25, 1463; Grimme, J. Comput. Chem. 2006, 27, 1787). Summation of empirical pair potential terms is performed over all atom pairs in the reference cell and over atoms in shells of neighboring cells until convergence of the dispersion energy is obtained. For DFT‐D3, the definition of coordination numbers has to be modified with respect to the molecular version to ensure convergence. The effect of three‐center terms as implemented in the original molecular DFT‐D3 version is investigated. The empirical parameters are taken from the original DFT‐D3 version where they had been optimized for a reference set of small molecules. As the coordination numbers of atoms in bulk and surfaces are much larger than in the reference compounds, this effect has to be discussed. The results of test calculations for bulk properties of metals, metal oxides, benzene, and graphite indicate that the original parameters are also suitable for solid‐state systems. In particular, the interlayer distance in bulk graphite and lattice constants of molecular crystals is considerably improved over standard functionals. With the molecular standard parameters (Grimme et al., J. Chem. Phys. 2010, 132, 154104; Grimme, J. Comput. Chem. 2006, 27, 1787) a slight overbinding is observed for ionic oxides where dispersion should not contribute to the bond. For simple adsorbate systems, such as Xe atoms and benzene on Ag(111), the DFT‐D implementations reproduce experimental results with a similar accuracy as more sophisticated approaches based on perturbation theory (Rohlfing and Bredow, Phys. Rev. Lett. 2008, 101, 266106). © 2012 Wiley Periodicals, Inc.     

Correlated ab initio quantum chemical calculations of di- and trisaccharide conformations

High level correlated quantum chemical calculations, using MP2 and local MP2 theory, have been performed for conformations of the disaccharide, beta-maltose, and the trisaccharide, 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranose. For beta-maltose, MP2 and local MP2 calculations using the 6-311++G** basis set are in good agreement, predicting a global minimum gas-phase conformation with a counterclockwise hydrogen bond network and the experimentally-observed intersaccharide hydrogen bonding arrangement. For conformations of 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranose, MP2/6-311++G**, and local MP2/6-311++G** calculations do not provide a consensus prediction of relative energetics, with the MP2 method finding large differences in stability between extended and folded trisaccharide conformations. Local MP2 calculations, less susceptible to intramolecular basis set superposition errors, predict a narrower range of trisaccharide energetics, in line with estimates from Hartree-Fock theory and B3LYP and BP86 density functionals. All levels of theory predict compact, highly hydrogen-bonded conformations as lowest in energy on the in vacuo potential energy surface of the trisaccharide. These high level, correlated local MP2/6-311++G** calculations of di- and trisaccharide energetics constitute potential reference data in the development and testing of improved empirical and semiempirical potentials for modeling of carbohydrates in the condensed phase.     

Interaction energies for the purine inhibitor roscovitine with cyclin-dependent kinase 2: correlated ab initio quantum-chemical, DFT and empirical calculations

The interaction between roscovitine and cyclin-dependent kinase 2 (cdk2) was investigated by performing correlated ab initio quantum-chemical calculations. The whole protein was fragmented into smaller systems consisting of one or a few amino acids, and the interaction energies of these fragments with roscovitine were determined by using the MP2 method with the extended aug-cc-pVDZ basis set. For selected complexes, the complete basis set limit MP2 interaction energies, as well as the coupled-cluster corrections with inclusion of single, double and noninteractive triples contributions [CCSD(T)], were also evaluated. The energies of interaction between roscovitine and small fragments and between roscovitine and substantial sections of protein (722 atoms) were also computed by using density-functional tight-binding methods covering dispersion energy (DFTB-D) and the Cornell empirical potential. Total stabilisation energy originates predominantly from dispersion energy and methods that do not account for the dispersion energy cannot, therefore, be recommended for the study of protein-inhibitor interactions. The Cornell empirical potential describes reasonably well the interaction between roscovitine and protein; therefore, this method can be applied in future thermodynamic calculations. A limited number of amino acid residues contribute significantly to the binding of roscovitine and cdk2, whereas a rather large number of amino acids make a negligible contribution.     

Molecular structure and vibrational spectra of phenobaraitone by density functional theory and ab initio hartree-Fock calculations

Quantum chemistry calculations have been performed using Gaussian03 program to compute optimized geometry, harmonic vibrational frequency along with intensities in IR and Raman spectra at RHF/6-31++G** and B3LYP/6-31++G** levels for phenobarbitone (C₁₂H₁₂N₂O₃) in the ground state. The scaled harmonic vibrational frequencies have been compared with experimental FT-IR and FT-Raman spectra. Theoretical vibrational spectra of the title compound were interpreted by means of potential energy distributions (PEDs) using MOLVIB program. A detailed interpretation of the infrared spectra of the title compound is reported. On the basis of the agreement between the calculated and observed results, the assignments of fundamental vibrational modes of phenobarbitone were examined and some assignments were proposed. The theoretical spectrograms for FT-IR and FT-Raman spectra of the title compound have been constructed.     

Reparameterization of a meta-generalized gradient approximation functional by combining TPSS exchange with τ1 correlation

We present a new local meta-GGA exchange-correlation density functional by combining the TPSS meta-GGA exchange and the τ1 meta-GGA correlation functionals. The TPSS meta-GGA exchange-correlation and the τ1 meta-GGA correlation functionals have been implemented in the deMon code. The parameters in the τ1 meta-GGA correlation model are reoptimized in a synchronized way to match the original TPSS meta-GGA exchange counterpart. This reparametrized meta-GGA functional is referred to as “TPSSτ3”. The TPSSτ3 and TPSSτ1 meta-GGAs are validated using a test set that consists of covalent molecules, hydrogen-bonded complexes, and van der Waals interactions. The calculated results from TPSSτ1 and TPSSτ3 are analyzed and compared with reliable experimental data and theoretical data, as well as with those from Bmτ1 and TPSS calculations. The τ1 correlation model describes the aromatic compounds better than TPSS. TPSSτ3 yields satisfactory results for the covalent molecules, the hydrogen-bonded complexes, and the van der Waals complexes in the test set compared with TPSS, Bmτ1 and TPSSτ1. Contribution to the Serafin Fraga Memorial Issue.     

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.     

低阶煤中氢键作用的色散校正密度泛函理论研究

本研究以苯酚…苯酚、苯酚…苯、苯酚…二苯醚、苯酚…喹啉和苯甲酸…苯甲酸为对象，采用色散校正的密度泛函理论分别研究褐煤中自缔合OH、OH-π、OH-醚O、OH-N和COOH-COOH之间形成的氢键。此外，还研究了氢键供体中取代基(CH₃-、CH₃O-、OH-、NH₂-、COOH-和NO₂-)对氢键的影响。对上述复合物进行了几何优化，并计算了能量、Mulliken电荷分布及振动频率。从优化的结构中可以看出上述复合物之间都存在氢键，所有复合物中O-H键键长都比苯酚中自由羟基的长，这表明这些复合物之间存在相互作用。其中，羧酸…羧酸复合物中O-H键的键长最长。此外，通过Mulliken电荷分布可看出上述复合物之间存在电荷转移。基于振动频率分析，所有的O-H键伸缩振动都发生了红移，尤其是羧酸…羧酸和苯酚…喹啉复合物，这可为煤中羟基振动的红外光谱分析提供依据。根据键能不同氢键强度按以下顺序依次递减：COOH-COOH>OH-N > 自缔合OH≈OH-醚O > OH-π，这与振动频率的分析结果一致。此外，不同取代基对氢键作用的影响不同。     

Improved parameterization of the quantum harmonic oscillator model based on localized wannier functions to describe Van der Waals interactions in density functional theory

Recently, the quantum harmonic oscillator model has been combined with maximally localized Wannier functions to account for long‐range dispersion interactions in density functional theory calculations (Silvestrelli, J. Chem. Phys. 2013, 139, 054106). Here, we present a new, improved set of values for the three parameters involved in this scheme. To test the new parameter set we have computed the potential energy curves for various systems, including an isolated Ar2 dimer, two N2 dimers interacting within different configurations, and a water molecule physisorbed on pristine graphene. While the original set of parameters generally overestimates the interaction energies and underestimates the equilibrium distances, the new parameterization substantially improves the agreement with experimental and theoretical reference values. © 2016 Wiley Periodicals, Inc.     

Equilibrium structure of the carbon dioxide-water complex in the gas phase: an ab initio and density functional study

High-level ab initio (MP2/6-311++G(2d,2p) geometry, Gaussian-2, MP4(SDTQ) and QCISD(T) binding energies) and density-functional (Becke3LYP/6-311++G(2df,2pd)) calculations have been performed on the charge-transfer complex between water and carbon dioxide. The complex appears to have two equivalent non-planar minima of C_s symmetry. Minima are separated by transition states with C₁ symmetry, whereas the totally planar structure with C_2v symmetry is a second-order transition state. All the critical points lie at approximately the same energy (less than 0.05 Kj mol⁻¹ difference). Therefore, the experimentally observable structure should be planar. The best equilibrium intermolecular distance for this complex calculated at the MP2/6-311++G(2d,2p) level is 2.800 Å. Our best estimate of the observable intermolecular distance (corrected for anharmonicity) is 2.84 Å, in agreement with the experimentally derived value of 2.836 Å. Our best estimate of the binding energy at the QCISD(T) level, taking into account the variation of the distance owing to anharmonicity and the use of more sophisticated theoretical treatments, is −12.0 ± 0.2 kJ mol⁻¹. Our best estimate of the barrier to internal rotation, also at the MP2/6-311++G(2d,2p) level, is 4.0 kJ mol⁻¹, outside the error limits of the experimental determination (3.64 ± 0.04 kJ mol⁻¹). Density functional theory at the level employed here gives an equilibrium intermolecular distance that is too large (2.857 Å), a binding energy that is too small (8.1 kJ mol⁻¹), attributable neither to geometry nor to the basis set, and also a barrier to internal rotation that is slightly too small (3.39 kJ mol⁻¹). The overall picture is, however, reasonably good.     

Effect of differently terminal groups of poly(amido‐amine) dendrimers on dispersion stability of nano‐silica and ab initio calculations

To investigate the effect of differently terminal groups of the lowest‐order generation poly(amido‐amine) dendrimers on dispersion stability of nano‐silica, the four types of G0‐CH₂CH₃ (G0E), G0‐CH₂CH₂CH₂CH₃ (G0B), G0‐NH₂ (G0N), and G0‐COOH (G0C) dendrimer molecules are used to modify the silica based on the dry modification. The zeta potential, the surface charge density, and the storage stability of kinds of modified SiO₂ dispersion systems have been studied. The results show that the effect of carboxyl groups on dispersion stability is stronger than that of the other groups such as the amine and alkyl groups. Mulliken charge distributions of the main active sites are analyzed through the conductor‐like polarizable calculation model (CPCM) on basis of the density functional theory (DFT) method, indicating the formation of chemical bonding between the modifiers and SiO₂ particles. The most stable SiO₂ dispersion system modified by G0‐COOH dendrimer molecule is obtained due to the combined effect including the hydrogen bonding, electrostatic repulsive force, and the steric hindrance of the terminal groups. Copyright © 2015 John Wiley & Sons, Ltd.     

Theoretical molecular structures for partially bonded complexes of trimethylamine with SO2 and SO3: ab initio and density functional theory calculations

The self-consistent reaction field (SCRF) method based on Onsager's reaction field theory is applied to investigate the effect of polar media on molecular structures of complexes of trimethylamime (TMA) with SO_x (x=2,3). The calculated SCRF N–S bond lengths at the MPW1PW91/6-311+G(3df) level are in satisfactory agreement with the experimental N–S bond lengths for the TMA–SO_x upon crystallization. The results are enough to demonstrate the usefulness of the reaction field theory in providing qualitative understanding of the medium effect on the partially bonded system such as TMA–SO_x.     

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     

Noncovalent interactions in inorganic supramolecular chemistry based in heavy metals. Quantum chemistry point of view

Complexity is a concept that is being considered in chemistry as it has shown potential to reveal interesting phenomena. Thus, it is possible to study chemical phenomena in a new approach called systems chemistry. The systems chemistry has an organization and function, which are regulated by the interactions among its components. At the simplest level, noncovalent interactions between molecules can lead to the emergence of large structures. Consequently, it is possible to go from the molecular to the supramolecular systems chemistry, which aims to develop chemical systems highly complex through intra- and intermolecular forces. Proper use of the interactions previously mentioned allow a glimpse of supramolecular system chemistry in many tasks such as structural properties reflecting certain behaviors in the chemistry of materials, for example, electrical and optical, processes of molecular recognition and among others. In the last time, within this area, inorganic supramolecular systems chemistry has been developed. Those systems have a structural orientation which is defined by certain forces that predominate in the associations among molecules. It is possible to recognize these forces as hydrogen bonding, π-π stacking, halogen bonding, electrostatic, hydrophobic, charge transfer, metal coordination, and metallophilic interactions. The presence of these forces in supramolecular system yields certain properties such as light absorption and luminescence. The quantum theoretical modeling plays an important role in the designing of the supramolecular system. The goal is to apply supramolecular principles in order to understand the associated forces in many inorganic molecules that include heavy metals for instance gold, platinum, and mercury. Relevant systems will be studied in detail, considering functional aspects such as enhanced coordination of functionalized molecular self-assembly, electronic and optoelectronic properties.