Implementation and Performance of DFT-D with Respect to Basis Set and Functional for Study of Dispersion Interactions in Nanoscale Aromatic Hydrocarbons

The Accuracy of Density Functional Theory in the Description of Cation-π and π-Hydrogen Bond Interactions

Cation-π and π-hydrogen bond interactions are ubiquitous in protein folding, molecular recognition, and ligand-receptor associations. As such systems are routinely studied at the DFT level, it becomes essential to understand the underlying accuracy of the plethora of density functionals currently available for the description of these interactions. For that purpose, we carried out theoretical calculations on two small model systems (benzene-Na(+) and benzene-H2O) that represent a paradigm for those intermolecular interactions and systematically tested 46 density functionals against the results of high-level post-HF methods, ranging from MP2 to extrapolated CCSD(T)/CBS. A total of 13 basis sets were also tested to examine the convergence of the interaction energy with basis set size. The convergence was surprisingly fast, with deviations below 0.2 kcal/mol for double-ζ polarized basis sets with diffuse functions. Concerning functional benchmarking, the Truhlar group functionals were particularly well suited for the description of the π-hydrogen bond interactions. In the case of cation-π interactions, there was not a clear correlation between accuracy and functional sophistication. Despite the large number of functionals predicting interaction energies within chemical accuracy (five for π-hydrogen bond and 20 for cation-π interactions), not a single functional has shown chemical accuracy in both cases. Moreover, if we calculate the average error for these two interactions, only two density functionals resulted in an average error below 1.0 kcal/mol (M06 and HCTH, with average errors of 0.6 and 0.8 kcal/mol). The obtained results serve as a guide for future computer simulations on this kind of system.

Aromatic rings in chemical and biological recognition: energetics and structures

This review describes a multidimensional treatment of molecular recognition phenomena involving aromatic rings in chemical and biological systems. It summarizes new results reported since the appearance of an earlier review in 2003 in host-guest chemistry, biological affinity assays and biostructural analysis, data base mining in the Cambridge Structural Database (CSD) and the Protein Data Bank (PDB), and advanced computational studies. Topics addressed are arene-arene, perfluoroarene-arene, S⋅⋅⋅aromatic, cation-π, and anion-π interactions, as well as hydrogen bonding to π systems. The generated knowledge benefits, in particular, structure-based hit-to-lead development and lead optimization both in the pharmaceutical and in the crop protection industry. It equally facilitates the development of new advanced materials and supramolecular systems, and should inspire further utilization of interactions with aromatic rings to control the stereochemical outcome of synthetic transformations.

Theoretical study of polaron formation in poly(G)-poly(C) cations

Polaron formation in poly(G)-poly(C) cations is investigated with density functional theory (DFT) and molecular mechanics (MM) employing a two-layer ONIOM method. In these calculations, the high layer, composed of all complementary base pairs, is treated by a DFT method, while the low layer, which includes the sugar-phosphate backbone, counterions and water molecules, is described by the AMBER force field. The high layer is the model system in which the charge transfer takes place. According to our calculations, three or four guanines move in a paddle-like fashion when an electron is removed from the neutral model system. In the cation model system, about 80% of the charge is delocalized onto the guanine residues, and the remaining charge is delocalized onto the cytosine residues. This happens because guanine has a lower ionization potential (IP) than cytosine. The counterions and water molecules in the low layer are important in the geometry optimization. The optimized geometry of the model system is closer to the standard B-form structure when counterions and water molecules are included than when they are omitted. Comparison of the optimized neutral and cationic model systems reveals a polaron in poly(G)-poly(C) cations extending from the first to the third guanine. It is demonstrated that the position of counterions and the number of surrounding water molecules can affect polaron formation.

Assessment of DFT and DFT-D for Potential Energy Surfaces of Rare Gas Trimers-Implementation and Analysis of Functionals and Extrapolation Procedures

Given the recent developments in methodology associated with the accurate computation of molecular systems with weak interactions, it is of particular interest to revisit systems that are notoriously challenging for determining reliable potential energy surface (PES) descriptions. Additionally, challenges associated with carrying out complete basis set extrapolation procedures and treatment of basis set superposition error (BSSE) are of importance in these descriptions. In this work, investigation into the ability to accurately predict the potential energy surfaces of the main Rg3 molecules (Rg = He, Ne, Ar) is made across a range of wave function types and large basis sets, including the use of several established extrapolation procedures and counterpoise corrections. Wave function types span most classes of density functional types, including the newest DFT-D schemes, and are benchmarked against high accuracy CCSD(T)/CBS methodology. Study of such systems is valuable, as they serve as simple models for many complex properties, most importantly n-body weak interaction energies.

Benchmark calculations of water-acene interaction energies: Extrapolation to the water-graphene limit and assessment of dispersion-corrected DFT methods

In a previous study (J. Phys. Chem. C, 2009, 113, 10242-10248) we used density functional theory based symmetry-adapted perturbation theory (DFT-SAPT) calculations of water interacting with benzene (C(6)H(6)), coronene (C(24)H(12)), and circumcoronene (C(54)H(18)) to estimate the interaction energy between a water molecule and a graphene sheet. The present study extends this earlier work by use of a more realistic geometry with the water molecule oriented perpendicular to the acene with both hydrogen atoms pointing down. We also include results for an intermediate C(48)H(18) acene. Extrapolation of the water-acene results gives a value of -3.0 +/- 0.15 kcal mol(-1) for the binding of a water molecule to graphene. Several popular dispersion-corrected DFT methods are applied to the water-acene systems and the resulting interacting energies are compared to results of the DFT-SAPT calculations in order to assess their performance.

On the Structure and Geometry of Biomolecular Binding Motifs (Hydrogen-Bonding, Stacking, X-H···π): WFT and DFT Calculations

The strengths of noncovalent interactions are generally very sensitive to a number of geometric parameters. Among the most important of these parameters is the separation between the interacting moieties (in the case of an intermolecular interaction, this would be the intermolecular separation). Most works seeking to characterize the properties of intermolecular interactions are mainly concerned with binding energies obtained at the potential energy minimum (as determined at some particular level of theory). In this work, in order to extend our understanding of these types of noncovalent interactions, we investigate the distance dependence of several types of intermolecular interactions, these are hydrogen bonds, stacking interactions, dispersion interactions, and X-H···π interactions. There are several methods that have traditionally been used to treat noncovalent interactions as well as many new methods that have emerged within the past three or four years. Here we obtain reference data using estimated CCSD(T) values at the complete basis set limit (using the CBS(T) method); potential energy curves are also produced using several other methods thought to be accurate for intermolecular interactions, these are MP2/cc-pVTZ, MP2/aug-cc-pVDZ, MP2/6-31G*(0.25), SCS(MI)-MP2/cc-pVTZ, estimated MP2.5/CBS, DFT-SAPT/aug-cc-pVTZ, DFT/M06-2X/6-311+G(2df,2p), and DFT-D/TPSS/6-311++G(3df,3pd). The basis set superposition error is systematically considered throughout the study. It is found that the MP2.5 and DFT-SAPT methods, which are both quite computationally intensive, produce potential energy curves that are in very good agreement to those of the reference method. Among the MP2 techniques, which can be said to be of medium computational expense, the best results are obtained with MP2/cc-pVTZ and SCS(MI)-MP2/cc-pVTZ. DFT-D/TPSS/6-311++G(3df,3pd) is the DFT-based method that can be said to give the most well-balanced description of intermolecular interactions.