Finding Critical Points and Reconstruction of Electron Densities on Grids

Alberto Otero-de-la-Roza
DOI: https://doi.org/10.1063/5.0090232
IF: 3.5
2022-05-27
Physics Today
Abstract:The Journal of Chemical Physics, Volume 0, Issue ja, -Not available-. Bader's quantum theory of atoms in molecules (QTAIM) is one of the most popular ways of extracting chemical insight from the results of quantum mechanical calculations. One of the basic tasks in QTAIM is to locate the critical points of the electron density and calculate various quantities (density, Laplacian,...) on them, since these have been found to correlate with molecular properties of interest. If the electron density is given analytically, this process is relatively straightforward. However, locating the critical points is more challenging if the density is known only on a three-dimensional uniform grid. A density grid is common in periodic solids because it is the natural expression for the electron density in plane-wavecalculations. In this article, I explore the reconstruction of the electron density from a grid and its use in critical point location. The proposed reconstruction method employs polyharmonic spline interpolation combined with a smoothing function based on the promolecular density. The critical point search based on this reconstruction is fast, accurate, trivially parallelizable, works for periodic and non-periodic systems, does not present directional lattice bias when the grid is non-orthogonal, and mostly locates allcritical points of the underlying electron density. The proposed method also provides an accurate reconstruction of the electron density over the space spanned by the grid, which may be useful in other contexts besides critical point localization. Refining the grid in the few systems where the proposed method fails to find the complete critical point topology is a viable solution.
physics, multidisciplinary
What problem does this paper attempt to address?