Efficient Evaluation of Molecular Electrostatic Potential in Large Systems

An algorithm for the efficient computation of molecular electrostatic potential is reported. It is based on the partition/expansion of density into (pseudo) atomic fragments with the method of Deformed Atoms in Molecules, which allows to compute the potential as a sum of atomic contributions. These contributions are expressed as a series of irregular spherical harmonics times effective multipole moments and inverse multipole moments, including short-range terms. The problem is split into two steps. The first one consists of the partition/expansion of density accompanied by the computation of multipole moments, and its cost depends on the size of the basis set used in the computation of electron density within the Linear Combination of Atomic Orbitals framework. The second one is the actual computation of the electrostatic potential from the quantities calculated in the first step, and its cost depends on the number of computation points. For a precision in the electrostatic potential of six decimal figures, the algorithm leads to a dramatic reduction of the computation time with respect to the calculation from electron density matrix and integrals involving basis set functions.

[1]  James J. P. Stewart,et al.  Optimization of parameters for semiempirical methods VI: more modifications to the NDDO approximations and re-optimization of parameters , 2012, Journal of Molecular Modeling.

[2]  J. Murray,et al.  An Overview of Strengths and Directionalities of Noncovalent Interactions: σ-Holes and π-Holes , 2019, Crystals.

[3]  J. Rico,et al.  Translation of real solid spherical harmonics , 2013 .

[4]  Timothy Clark,et al.  Halogen bonding: the σ-hole , 2007 .

[5]  Shridhar R. Gadre,et al.  An appraisal of Poincaré-Hopf relation and application to topography of molecular electrostatic potentials. , 2008, The Journal of chemical physics.

[6]  Jacopo Tomasi,et al.  The electrostatic molecular potential as a tool for the interpretation of molecular properties , 1973 .

[7]  Shridhar R. Gadre,et al.  Topology of molecular electron density and electrostatic potential with DAMQT , 2017, Comput. Phys. Commun..

[8]  I. S. Gradshteyn,et al.  Table of Integrals, Series, and Products , 1976 .

[9]  W. H. Henneker,et al.  Molecular Charge Distributions and Chemical Binding , 1967 .

[10]  J. Murray,et al.  Computational characterization of nucleotide bases: Molecular surface electrostatic potentials and local ionization energies, and local polarization energies , 2001 .

[11]  Rafael López,et al.  Electrostatic potentials and fields from density expansions of deformed atoms in molecules , 2004, J. Comput. Chem..

[12]  S. R. Gadre,et al.  Exploring the Gradient Paths and Zero Flux Surfaces of Molecular Electrostatic Potential. , 2016, Journal of chemical theory and computation.

[13]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[14]  Rafael López,et al.  Analytical method for the representation of atoms‐in‐molecules densities , 2004, J. Comput. Chem..

[15]  M. Alderton,et al.  Distributed multipole analysis , 2006 .

[16]  Jacopo Tomasi,et al.  Electronic Molecular Structure, Reactivity and Intermolecular Forces: An Euristic Interpretation by Means of Electrostatic Molecular Potentials , 1978 .

[17]  R. N. Shirsat,et al.  TLC Detection and Theoretical Structure Elucidation of Nitrogen Containing Compounds with Cobalt Thiocyanate , 2019, Analytical Chemistry Letters.

[18]  I. Shrivastava,et al.  Shapes and sizes of molecular anions via topographical analysis of electrostatic potential , 1991 .

[19]  R. López,et al.  Analysis of the molecular density , 1999 .

[20]  J. Elguero,et al.  Complexes of O=C=S with Nitrogen Bases: Chalcogen Bonds, Tetrel Bonds, and Other Secondary Interactions. , 2018, Chemphyschem : a European journal of chemical physics and physical chemistry.

[21]  R. S. Mulliken.  Electronic Population Analysis on LCAO–MO Molecular Wave Functions. I , 1955 .

[22]  Gilbert N. Lewis,et al.  The Atom and the Molecule , 1916, Resonance.

[23]  Peter Politzer,et al.  Relationships between computed molecular properties and solute-solvent interactions in supercritical solutions , 1993 .

[24]  W. Heisenberg,et al.  Zur Quantentheorie der Molekeln , 1924 .

[25]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[26]  Shridhar R. Gadre,et al.  Basis set dependence of the molecular electrostatic potential topography. A case study of substituted benzenes , 1995 .

[27]  S. R. Gadre,et al.  Molecular Electrostatic Potential-Based Atoms in Molecules: Shielding Effects and Reactivity Patterns , 2016 .

[28]  Single-centre expansion of Gaussian basis functions and the angular decomposition of their overlap integrals , 1989 .

[29]  W. J. Orville-Thomas.  Atoms in Molecules — a Quantum Theory , 1996 .

[30]  Peter Politzer,et al.  Use of the electrostatic potential at the molecular surface to interpret and predict nucleophilic processes , 1990 .

[31]  Shridhar R. Gadre,et al.  Maximal and minimal characteristics of molecular electrostatic potentials , 1990 .

[32]  J. Murray,et al.  Molecular surface electrostatic potentials in the analysis of non-hydrogen-bonding noncovalent interactions , 1993, Proceedings / Indian Academy of Sciences.

[33]  Julio Caballero,et al.  Rationalizing the stability and interactions of 2,4-diamino-5-(4-chlorophenyl)-6-ethylpyrimidin-1-ium 2-hydroxy-3,5-dinitrobenzoate salt , 2019, Journal of Molecular Structure.

[34]  J. Murray,et al.  Relationships of critical constants and boiling points to computed molecular surface properties , 1993 .