求助：受限水的垂直介电常数计算问题

aj1685889 · 发表于 Post on 2025-1-11 23:27:09

本帖最后由 aj1685889 于 2025-2-11 11:27 编辑

我的导师最进让我做关于受限体系的垂直介电常数计算，体系非常简单，就是水在双层石墨烯之间，参考的文献我放在附件里面了，编写的in文件我放在尾端的附件了也有我的相关疑问，关于密度问题，我也是使用了NPT,1bar条件处理得到的适合的狭缝水密度填充的，文献中的计算方法就说是根据原子的电荷密度进行积分得到局部极化密度，之后再次积分便得到了总极化密度，但是我在尝试了很多次运行，总是会出现一个问题就是计算的垂直介电常数体相区的值在负值，这一点我没有明白为什么，尝试过修改lj参数，尝试过续跑，基本上都没有成功，包括使用文献中的SF以及非周期Ewald3dc的BH方法，基本上都是处于0刻度线下，而不是体相介电常数的倒数，如水=1/71=0.013，现在非常头大，我的编程计算是使用的Fortran，关于电荷密度使用的是固定电荷乘以z轴位置，积分方法使用的梯形积分法，如果有大佬懂，我可以发给您查看程序，非常希望有人能够解答，苦恼了很久了！！

wudiazhu · 发表于 Post on 2025-1-13 14:41:35

我觉得你可以拿你的脚本算一个没有confinement的水盒子试试，看下算出来是不是1/71。如果是的话，检查一下介电常数里的偶极波动项，把每一步的偶极平方打印出来和偶极系综平均的平方对比看看。

aj1685889 · 发表于 Post on 2025-1-13 16:57:37

wudiazhu 发表于 2025-1-13 14:41
我觉得你可以拿你的脚本算一个没有confinement的水盒子试试，看下算出来是不是1/71。如果是的话，检查一下 ...

好的，谢谢大佬指导，我马上提交任务试一试

aj1685889 · 发表于 Post on 2025-2-9 13:24:02

本帖最后由 aj1685889 于 2025-2-11 11:27 编辑

1111

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Lammps] 求助：受限水的垂直介电常数计算问题