Jingzhen Hu, Shan Zhao and Weihua Geng. (2018). Accurate p$K_a$ Computation Using Matched Interface and Boundary (MIB) Method Based Poisson-Boltzmann Solver. Communications in Computational Physics. 23 (2). 520-539. doi:10.4208/cicp.OA-2017-0078