Following Becke, 102 we define the Voronoi polyhedra using a two-center coordinate system, via confocal elliptical coordinates r = (λ AB , µ AB , φ AB ) for the pair of nuclei situated at R A and R B . In this coordinate system, φ AB is the angle between r and the internuclear axis (A10)