Assessment of DFT Functionals for QTAIM Topological Analysis of Halogen Bonds with Benzene

Halogen bonding, a noncovalent interaction between a halogen atom and a nucleophilic site, is receiving a growing attention in the chemical community stimulating a large number of theoretical investigations. The density functional theory (DFT) approach revealed to be one of the most suitable methods owing to its accuracy and low computational cost. We report here a detailed analysis of the performance of an extensive set of DFT functionals in reproducing accurate binding energies and topological properties for the halogen-bonding interaction of either NCX or PhX molecules (X = F, Cl, Br, I) with the aromatic system of benzene in the T-shaped configuration. It was found that the better performance for both sets of properties is provided by a small subset of functionals able to take into account, implicitly or explicitly (by inclusion of an additive pairwise potential), the dispersion contribution, that is, ωB97X, M06-2X, M11, mPW2PLYP-D, and B2PLYP-D3.

Association of symmetrical alkane diols with pyridine: DFT/GIAO calculation of 1 H NMR chemical shifts

Proton nuclear magnetic resonance (NMR) shifts of the free diol and of its 1 : 1 and 1 : 2 hydrogen-bonded complexes with pyridine have been computed for five symmetrical alkane diols on the basis of density functional theory, by applying the gauge-including atomic orbital method to geometry-optimized conformers. For certain conformers, intramolecular OH···OH interactions, evidenced by high NMR OH proton shifts, are further enhanced on going from the free diol to the corresponding 1 : 1 diol/pyridine complex. This is confirmed by atoms-in-molecules and non-covalent interaction plots. The computed OH and CH proton shifts for the diol and the two complexes correlate well with values obtained by analysing data from the NMR titration of the diols in benzene against pyridine. Shift values for the diols in neat pyridine are calculated by weighting the shifts of the various protons in the three forms (free diol, 1 : 1 and 1 : 2 diol/pyridine complexes) according to the experimentally determined association constants. The results are in good agreement with those observed, and after empirical scaling, the root mean square difference is 0.18 ppm. Copyright © 2016 John Wiley & Sons, Ltd.

GenLocDip: A Generalized Program to Calculate and Visualize Local Electric Dipole Moments

Local dipole moments (i.e., dipole moments of atomic or molecular subsystems) are essential for understanding various phenomena in nanoscience, such as solvent effects on the conductance of single molecules in break junctions or the interaction between the tip and the adsorbate in atomic force microscopy. We introduce GenLocDip, a program for calculating and visualizing local dipole moments of molecular subsystems. GenLocDip currently uses the Atoms-In-Molecules (AIM) partitioning scheme and is interfaced to various AIM programs. This enables postprocessing of a variety of electronic structure output formats including cube and wavefunction files, and, in general, output from any other code capable of writing the electron density on a three-dimensional grid. It uses a modified version of Bader's and Laidig's approach for achieving origin-independence of local dipoles by referring to internal reference points which can (but do not need to be) bond critical points (BCPs). Furthermore, the code allows the export of critical points and local dipole moments into a POVray readable input format. It is particularly designed for fragments of large systems, for which no BCPs have been calculated for computational efficiency reasons, because large interfragment distances prevent their identification, or because a local partitioning scheme different from AIM was used. The program requires only minimal user input and is written in the Fortran90 programming language. To demonstrate the capabilities of the program, examples are given for covalently and non-covalently bound systems, in particular molecular adsorbates. © 2016 Wiley Periodicals, Inc.

Natural occupation numbers in two-electron quantum rings

Natural orbitals (NOs) are central constituents for evaluating correlation energies through efficient approximations. Here, we report the closed-form expression of the NOs of two-electron quantum rings, which are prototypical finite-extension systems and new starting points for the development of exchange-correlation functionals in density functional theory. We also show that the natural occupation numbers for these two-electron paradigms are in general non-vanishing and follow the same power law decay as atomic and molecular two-electron systems.

Experimental and theoretical charge density distribution in Pigment Yellow 101

The charge density distribution in 2,2'-dihydroxy-1,1'-naphthalazine (Pigment Yellow 101; P.Y.101) has been determined using high-resolution X-ray diffraction and multipole refinement, along with density functional theory calculations. Topological analysis of the resulting densities highlights the localisation of single/double bonds in the central C=N-N=C moiety of the molecule in its ground state. The density in the N-N is examined in detail, where we show that very small differences between experiment and theory are amplified by use of the Laplacian of the density. Quantification of hydrogen bonds highlights the importance of the intramolecular N-H···O interaction, known to be vital for retention of fluorescence in the solid state, relative to the many but weak intermolecular contacts located. However, a popular method for deriving H-bond strengths from density data appears to struggle with the intramolecular N-H···O interaction. We also show that theoretical estimation of anisotropic displacements for hydrogen atoms brings little benefit overall, and degrades agreement with experiment for one intra-molecular contact.

Prediction of Intramolecular Polarization of Aromatic Amino Acids Using Kriging Machine Learning

Present computing power enables novel ways of modeling polarization. Here we show that the machine learning method kriging accurately captures the way the electron density of a topological atom responds to a change in the positions of the surrounding atoms. The success of this method is demonstrated on the four aromatic amino acids histidine, phenylalanine, tryptophan, and tyrosine. A new technique of varying training set sizes to vastly reduce training times while maintaining accuracy is described and applied to each amino acid. Each amino acid has its geometry distorted via normal modes of vibration over all local energy minima in the Ramachandran map. These geometries are then used to train the kriging models. Total electrostatic energies predicted by the kriging models for previously unseen geometries are compared to the true energies, yielding mean absolute errors of 2.9, 5.1, 4.2, and 2.8 kJ mol(-1) for histidine, phenylalanine, tryptophan, and tyrosine, respectively.

Unusual bonding modes of perfluorobenzene in its polymeric (dimeric, trimeric and tetrameric) forms: entirely negative fluorine interacting cooperatively with entirely negative fluorine

The Laplacian of the charge density plot for the tetramer of perfluorobenzene, displaying the attraction between entirely negative fluorine atoms, as well as the windmill type intermolecular topology between those atoms.

Atomic decomposition of conceptual DFT descriptors: application to proton transfer reactions

We present an atomic decomposition of the molecular energy, reaction force and reaction flux, which is based on Bader's atoms-in-molecules theory.