where the sum runs over all particles and periodic images and is evaluated using the same Ewald method as (1). The locations of local maxima in the field intensity, i.e. "hot-spots", are approximated as x i,hot = x i ± dc 2 pi , and the collection of 2N values for |E (x i,hot )| make up the field distribution for a single polarization and configuration.