to the applied field. The geometries of the inclusions were generated by finding the bestfit ellipsoids to the Bader zero flux surfaces of the inclusions from the quantum dynamicsoptimized composites, where the ellipsoid volumes were constrained to be equal to the Bader volumes. Using the particle permittivity values derived in Sec. 3.1.2, we calculate the permittivities of composites with rod-and disk-shaped inclusions using the FEM [see Sec.