How can I minimize a potential energy function in 2D with periodic boundary conditions using python?

Question

Right now, I am trying to minimize an electrostatic potential energy function in 2D with periodic boundary conditions and two charges q1 and q2. Because of the conditions, the t