Towards Improving The Accuracy of Implicit Solvent Models and Understanding Electrostatic Catalysis in Complex Solvent Environment