wbn 发表于 2019-11-15 02:53
前面的那些系数是根据Born model 来的，应该等于e^2/4\pi\epsilon_0 ， 你可以自己验证一下数字对不对。分 ...

非常感谢老师的指导

前面的那些系数是根据Born model 来的，应该等于e^2/4\pi\epsilon_0 ， 你可以自己验证一下数字对不对。分子的solvation radius目前并没有什么好的方法可以计算，你可以计算单个分子的solvation energy (http://sobereva.com/333)然后反推出solvation radius， 或者算分子的范德华半径 （http://sobereva.com/190）， 只要方法能说的通就行。因为这个模型本来就是很粗糙的，分子也不可能都是球形。内重组能一般情况下确实比外重组能小一个数量级，但是我还是建议算一下加上去，反正也不麻烦