A new density functional (DF) of the generalized gradient approximation (GGA) type for general chemistry applications termed B97-D is proposed. It is based on Becke's power-series ansatz from 1997 and is explicitly parameterized by including damped atom-pairwise dispersion corrections of the form C(6) x R(-6). A general computational scheme for the parameters used in this correction has been established and parameters for elements up to xenon and a scaling factor for the dispersion part for several common density functionals (BLYP, PBE, TPSS, B3LYP) are reported. The new functional is tested in comparison with other GGAs and the B3LYP hybrid functional on standard thermochemical benchmark sets, for 40 noncovalently bound complexes, including large stacked aromatic molecules and group II element clusters, and for the computation of molecular geometries. Further cross-validation tests were performed for organometallic reactions and other difficult problems for standard functionals. In summary, it is found that B97-D belongs to one of the most accurate general purpose GGAs, reaching, for example for the G97/2 set of heat of formations, a mean absolute deviation of only 3.8 kcal mol(-1). The performance for noncovalently bound systems including many pure van der Waals complexes is exceptionally good, reaching on the average CCSD(T) accuracy. The basic strategy in the development to restrict the density functional description to shorter electron correlation lengths scales and to describe situations with medium to large interatomic distances by damped C(6) x R(-6) terms seems to be very successful, as demonstrated for some notoriously difficult reactions. As an example, for the isomerization of larger branched to linear alkanes, B97-D is the only DF available that yields the right sign for the energy difference. From a practical point of view, the new functional seems to be quite robust and it is thus suggested as an efficient and accurate quantum chemical method for large systems where dispersion forces are of general importance.

Semiempirical GGA-type density functional constructed with a long-range dispersion correction.

Molecular structure does not easily identify the intricate noncovalent interactions that govern many areas of biology and chemistry, including design of new materials and drugs. We develop an approach to detect noncovalent interactions in real space, based on the electron density and its derivatives. Our approach reveals the underlying chemistry that compliments the covalent structure. It provides a rich representation of van der Waals interactions, hydrogen bonds, and steric repulsion in small molecules, molecular complexes, and solids. Most importantly, the method, requiring only knowledge of the atomic coordinates, is efficient and applicable to large systems, such as proteins or DNA. Across these applications, a view of nonbonded interactions emerges as continuous surfaces rather than close contacts between atom pairs, offering rich insight into the design of new and improved ligands.

/pdf/revealing-noncovalent-interactions-1oqknfvmfi.pdf

Revealing noncovalent interactions.

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

QM/MM Methods for Biomolecular Systems

Dispersion corrections to standard Kohn–Sham density functional theory (DFT) are reviewed. The focus is on computationally efficient methods for large systems that do not depend on virtual orbitals or rely on separated fragments. The recommended approaches (van der Waals density functional and DFT-D) are asymptotically correct and can be used in combination with standard or slightly modified (short-range) exchange–correlation functionals. The importance of the dispersion energy in intramolecular cases (conformational problems and thermochemistry) is highlighted. © 2011 John Wiley & Sons, Ltd. WIREs Comput Mol Sci 2011 1 211-228 DOI: 10.1002/wcms.30

Density functional theory with London dispersion corrections

Introduction
PaDEL-Descriptor is a software for calculating molecular descriptors and fingerprints. The software currently calculates 797 descriptors (663 1D, 2D descriptors, and 134 3D descriptors) and 10 types of fingerprints. These descriptors and fingerprints are calculated mainly using The Chemistry Development Kit. Some additional descriptors and fingerprints were added, which include atom type electrotopological state descriptors, McGowan volume, molecular linear free energy relation descriptors, ring counts, count of chemical substructures identified by Laggner, and binary fingerprints and count of chemical substructures identified by Klekota and Roth.

Methods
PaDEL-Descriptor was developed using the Java language and consists of a library component and an interface component. The library component allows it to be easily integrated into quantitative structure activity relationship software to provide the descriptor calculation feature while the interface component allows it to be used as a standalone software. The software uses a Master/Worker pattern to take advantage of the multiple CPU cores that are present in most modern computers to speed up calculations of molecular descriptors.

Results
The software has several advantages over existing standalone molecular descriptor calculation software. It is free and open source, has both graphical user interface and command line interfaces, can work on all major platforms (Windows, Linux, MacOS), supports more than 90 different molecular file formats, and is multithreaded.

Conclusion
PaDEL-Descriptor is a useful addition to the currently available molecular descriptor calculation software. The software can be downloaded at http://padel.nus.edu.sg/software/padeldescriptor. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011

https://onlinelibrary.wiley.com/doi/pdf/10.1002/jcc.21707

PaDEL‐descriptor: An open source software to calculate molecular descriptors and fingerprints

Evaluation of human intestinal absorption data and subsequent derivation of a quantitative structure-activity relationship (QSAR) with the Abraham descriptors.

Additive models for the estimation of Abraham's molecular descriptors R2, π2H, Σα2H, Σβ2H, Σβ2O, and log L16 have been developed. For five of the six descriptors, one set of 81 atom and functional group fragments is capable of reproducing experimentally derived results with correlation coefficients ranging from 0.95 to 0.99. However, one descriptor, Σα2H, required an entirely separate set of 51 fragments to be developed, resulting in a correlation coefficient of 0.97. Of particular importance is the speed of calculation (approximately 700 molecules/min), allowing so-called “high-throughput screening”. Several applications of this model for molecules containing intramolecular interactions are discussed.

Estimation of Molecular Linear Free Energy Relation Descriptors Using a Group Contribution Approach

The suitability of a hybrid density functional to qualitatively reproduce geometric and energetic details of parallel pi-stacked aromatic complexes is presented. The hybrid functional includes an ad hoc mixture of half the exact (HF) exchange with half of the uniform electron gas exchange, plus Lee, Yang, and Parr's expression for correlation energy. This functional, in combination with polarized, diffuse basis sets, gives a binding energy for the parallel-displaced benzene dimer in good agreement with the best available high-level calculations reported in the literature, and qualitatively reproduces the local MP2 potential energy surface of the parallel-displaced benzene dimer. This method was further critically compared to high-level calculations recently reported in the literature for a range of pi-stacked complexes, including monosubstituted benzene-benzene dimers, along with DNA and RNA bases, and generally agrees with MP2 and/or CCSD(T) results to within +/-2 kJ mol(-1). We also show that the resulting BH&H binding energy is closely related to the electron density in the intermolecular region. The net result is that the BH&H functional, presumably due to fortuitous cancellation of errors, provides a pragmatic, computationally efficient quantum mechanical tool for the study of large pi-stacked systems such as DNA.

Hybrid density functional theory for π-stacking interactions : application to benzenes, pyridines, and DNA bases

The ability of functional groups to act as hydrogen bond acids and bases can be obtained from either equilibrium constants for 1:1 hydrogen bonding or overall hydrogen bond constants. Either method leads to structural constants for hydrogen bonding that in some way are analogous to substituent constants. Extensive lists of these functional group constants are reported. It is shown that those derived from overall hydrogen bond constants are the more useful in analyses of physicochemical and biochemical properties.

Hydrogen bond structural group constants.

Potential energy curves for the parallel-displaced, T-shaped and sandwich structures of the benzene dimer are computed with density fitted local second-order Moller–Plesset perturbation theory (DF-LMP2) as well as with the spin-component scaled (SCS) variant of DF-LMP2. While DF-LMP2 strongly overestimates the dispersion interaction, in common with canonical MP2, the DF-SCS-LMP2 interaction energies are in excellent agreement with the best available literature values along the entire potential energy curves. The DF-SCS-LMP2 dissociation energies for the three structures are also compared with new complete basis set estimates of the interaction energies obtained from accurate coupled cluster (CCSD(T)) and DF-SCS-MP2 calculations. Since LMP2 is essentially free of basis set superposition errors, counterpoise corrections are not required. As a result, DF-SCS-LMP2 is computationally inexpensive and represents an attractive method for the study of larger π-stacked systems such as truncated sections of DNA.

James Alexis Platts

Papers

Evaluation of human intestinal absorption data and subsequent derivation of a quantitative structure-activity relationship (QSAR) with the Abraham descriptors.

Estimation of Molecular Linear Free Energy Relation Descriptors Using a Group Contribution Approach

Hybrid density functional theory for π-stacking interactions : application to benzenes, pyridines, and DNA bases

Hydrogen bond structural group constants.

Calculation of intermolecular interactions in the benzene dimer using coupled-cluster and local electron correlation methods