Thanks to the favorable properties of the Delaunay triangulation (see Sec. II C 1), topological properties, such as the number of maxima and saddle points etc., converge almost instantly. Topographical properties (such as Bader charges) also converge quickly, as can be seen in Figs. DFT grids well before 1 million grid points. Results using the uniform-grid Multiwfn package, v3.8 with charge densities calculated using QChem, v5.0 are also shown (the same numbers are obtained if Psi4 v1.6.1 is used in place of QChem). Such uniform-grid methods routinely use tens of millions of grid points -we quote results obtained using Multiwfn's "Lunatic quality grid" (of order 10 million points) and using an even larger custom grid (of order 100 million points, obtained by specifying a custom spacing), which we denote as extralunatic, which was necessary to achieve convergence in all cases.