The RESP-dPol model assigns partial charges with a RESP-like fitting method in the context of the typed polarizabilities discussed above. It differs from standard RESP because the induced dipoles also contribute to the computed ESPs. Pioneering work by Cieplak et al. used a similar methodology and achieved reasonable agreement with experimental solvation free energies. In RESP-dPol, we adopt the two-stage hyperbolic restraints used in RESP to reduce nonphysical variations of partial charges. The first stage is used to make sure all polar regions are well-fitted, and a weak restraint (a = 0.005 a.u., b = 0.1 a.u.)