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