“Accurate p$K_a$ Computation Using Matched Interface and Boundary (MIB) Method Based Poisson-Boltzmann Solver” (2018) Communications in Computational Physics, 23(2), pp. 520–539. doi:10.4208/cicp.OA-2017-0078.