A domain optimization problem related to potential theory is studied by means of a relaxation approach where a concave functional on a given convex set is defined. The functional has a minimizing point that is the characteristic function of an optimal domain. As a consequence of the necessary conditions of optimality, the domain is the solution of a free boundary value problem.