计算化学公社

标题: 求助LAMMPS计算受限水的介电常数问题 [打印本页]

作者
Author: CSUDzh 时间: 2025-1-6 10:34
标题: 求助LAMMPS计算受限水的介电常数问题
本帖最后由 CSUDzh 于 2025-1-10 16:06 编辑

我正在尝试使用 lammps 计算石墨烯通道内不同距离处的介电常数。
根据文献中的公式，需要计算体系的总偶极矩和局域极化密度。我使用 compute chunk/atom bin/1d 命令将石墨烯通道沿z轴分成不同距离的块（其中还有一个整块是为了计算体系总的电偶极矩M，并划分了边界以去除石墨烯边界的效应）。由于关于chunk命令输出的是一个矢量，所以我编写了一个循环来提取矢量中的每个量来通过定义（单位体积内的电偶极矩）计算局部极化密度 P（z）。
文献中计算的介电常数在不同距离处与体介电常数具有相同的数量级，但我计算的介电常数却是10的15次方这么大的数量级。我尝试使用lammps自带的dipole/chunk命令计算每个chunk的电偶极矩、计算电偶极矩时减去块内静电荷和质心的乘积、重置 image 的值、并在计算电偶极矩时使用unwrapped coordinates，但结果仍然变化很小。
请问我对公式中物理量的计算方式有什么问题吗，或者是 lammps 命令的使用有问题吗？
(, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 9) (, 下载次数 Times of downloads: 4)

作者
Author: aj1685889 时间: 2025-2-9 13:29
我也在做这个计算，不过我的计算垂直介电常数的过程就是简单的积分，求出电荷密度就可以，但是脚本一直卡着无法进行下一步

欢迎光临计算化学公社 (http://bbs.keinsci.com/)