We propose an efficient, accurate method to integrate the basins of attraction of a smooth function defined on a general discrete grid and apply it to the Bader charge partitioning for the electron charge density. Starting with the evolution of trajectories in space following the gradient of charge density, we derive an expression for the fraction of space neighboring each grid point that flows to its neighbors. This serves as the basis to compute the fraction of each grid volume that belongs to a basin (Bader volume) and as a weight for the discrete integration of functions over the Bader volume. Compared with other grid-based algorithms, our approach is robust, more computationally efficient with linear computational effort, accurate, and has quadratic convergence. Moreover, it is straightforward to extend to nonuniform grids, such as from a mesh-refinement approach, and can be used to both identify basins of attraction of fixed points and integrate functions over the basins.

Original languageEnglish (US)
Article number064111
JournalJournal of Chemical Physics
Issue number6
StatePublished - Feb 14 2011

ASJC Scopus subject areas

  • General Physics and Astronomy
  • Physical and Theoretical Chemistry


Dive into the research topics of 'Accurate and efficient algorithm for Bader charge integration'. Together they form a unique fingerprint.

Cite this