Accurate p$K_a$ Computation Using Matched Interface and Boundary (MIB) Method Based Poisson-Boltzmann Solver. Commun. Comput. Phys. 2018, 23 (2), 520-539. https://doi.org/10.4208/cicp.OA-2017-0078.