Exploring Higher Accuracy Poisson-Boltzmann Methods For Biomolecular Simulations