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...

ver descrição completa

Detalhes bibliográficos
Autores: López Fernández, Rafael, Ema López, Ignacio de, García de la Vega, José Manuel, Ramírez Moreno, Guillermo, Martínez, Frank
Formato: artículo
Fecha de publicación:2019
País:España
Recursos:Universidad Autónoma de Madrid
Repositorio:Biblos-e Archivo. Repositorio Institucional de la UAM
Idioma:inglés
OAI Identifier:oai:repositorio.uam.es:10486/709307
Acesso em linha:http://hdl.handle.net/10486/709307
Access Level:acceso abierto
Palavra-chave:Atoms in molecules
Chemical reactivity
Electron density
Molecular electrostatic potential
Noncovalent interactions
Química
Descrição
Resumo: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