Gromacs求介电常数的两个功能:gmx current和gmx diploes的区别以及原理？

liu840397562

请问各位老师同学们，gmx current和gmx diploes计算介电常数的原理区别是什么呢？查了手册也没找到，网上找到了gmx current是计算的静态介电常数，使用公式为Epsilon(0) = 1 + 4 Pi/3 * (<M^2> - <M>^2)/(V k_B T)，而gmx dipoles使用的公式以及原理是什么呢？求各位赐教。

sobereva

gmx current输出两个epsilon，一个是纯粹基于偶极矩波动情况计算的，和gmx dipoles的一致；还一个似乎是依赖于流密度算的，和gmx dipoles给的有轻微差别。用前者就完了

liu840397562

sobereva 发表于 2023-11-10 07:26
gmx current输出两个epsilon，一个是纯粹基于偶极矩波动情况计算的，和gmx dipoles的一致；还一个似乎是依 ...

谢谢sob老师解答，我在找资料过程中发现，似乎还有一种公式形式与上文中所述很像，仅系数不一样为：Epsilon(0) = 1 + 1/3 * (<M^2> - <M>^2)/(V k_B T) 这两者似乎完全不一样，但均是通过偶极矩的涨落计算的，请问有什么联系吗老师?

nauxi

liu840397562 发表于 2023-11-10 10:45
谢谢sob老师解答，我在找资料过程中发现，似乎还有一种公式形式与上文中所述很像，仅系数不一样为：Epsil ...

请问你找到区别了吗，1/3这个分母貌似还有一个真空介电常数ε0