A first-order system least-squares finite element method for the Poisson-Boltzmann equation