The atoms in molecules theory provides a rigorous definition of chemical bonds for all types of molecules and solids [1][2][3][4][5][6][7][8][9][10] and has proven useful in analysis of the physical properties of insulators, pure metals, and alloys [4][5][6][7][8]. In particular, it has been observed that the strength of the bonding between a given pair of atoms in a molecule correlates with the values of the electron density at the bond critical point, r b [1]. A simple relationship between the binding energy and r b in periodic solids has also been reported [4,6,8]. The AIM theory also leads to unique partition of three-dimensional space into a collection of chemically identifiable regions called atomic basins (the atoms in a molecule or in a crystal). These are the most transferable pieces one can define in exhaustive partitioning of the real space [1]. In this chapter we describe the implementation of an algorithm that uses the r(r) topological information to determine the main elements of the AIM theory for periodic systems, and discuss an application to nanocatalysis.
Implementing the Determination of the Topological Properties of r(r) from a Three-dimensional GridThe CPs are usually calculated using the Newton-Raphson (NR) technique [11]. The NR algorithm starts from a truncated Taylor expansion at a point r ¼ r 0 þ h, about r 0 of a multidimensional scalar function ð'rðrÞÞ:where H is the Hessian (the Jacobian of 'rðrÞÞ at point r 0 . The best step, h, to move from the initial r 0 to the CP is h ¼ ÀH À1 'rðr 0 Þ. This correction is then used to obtain a vector r new ¼ r old þ th (t is a small value) and the process is iter-
231The Quantum Theory of Atoms in Molecules. Edited