We prove that the hexagonal lattice is a local minimizer, among all point configurations, of the interaction energy per unit volume for pair potentials that are completely monotonic functions of the square distance. This includes Gaussian interactions and power laws.
