The Voronoi functional is maximized by the Delaunay triangulation in the plane
H. Edelsbrunner, A. Glazyrin, O. Musin, A. Nikitenko, Combinatorica 37 (2017) 887–910.