On the stability domain of systems of four unit charges