In this dissertation, we mainly discuss the image approximation methods to the reaction field and their applications to electrostatic interactions in molecular dynamics simulations. The Poisson-Boltzman (PB) equation, a three-dimensional second or...