求助使用GROMACS模拟计算电解液介电常数时有关g-factor的问题

18763285297

各位老师好，我在做MD模拟计算电解液的介电常数，发现与实验值偏差较大。经过查阅文献，发现在一篇文章（C. Zhang, M. Sprik, Computing the dielectric constant of liquid water at constant dielectric displacement, Phys. Rev. B 93 (2016), 144201. https://doi.org/10.1103/ PhysRevB.93.144201.）中提到了一种方法：分别对分子进行电场E=0和电位移场D=0的MD模拟，计算g-factor, 然后复合两种方式计算出来的Gk(r):
再代入公式： 计算出介电常数。
目前有文章采用此方法进行了模拟 (T. Yamaguchi, H. Yamada, T. Fujiwara, T. Fujiwara, K. Hori, Simulations of dielectric constants and viscosities of organic electrolytes by quantum mechanics and molecular dynamics, J. Mol. Liq. 312 (2020) 113288, https://doi.org/ 10.1016/j.molliq.2020.113288.)，且模拟结果与实验值非常吻合：
我想请教一下如果想要进行D=0的模拟，是否要根据公式D =ε0E + P 计算出电场强度E的值呢？其中的P为偶极矩的密度，需要从x,y,z三个方向分别计算出矢量和，并代入计算出三个方向的Ex, Ey, Ez吗？
再次感谢各位大佬