想请教一下为什么我计算出来的偶极矩自相关联函数不是振荡的。
fix           3 all ave/correlate 5 10000 100000 v_dipole_x v_dipole_y v_dipole_z type auto file sic.dat ave run
分子动力学时间步长设置为0.005ps。

wyw2000 发表于 2023-11-3 11:08
楼主现在关于节电常数的计算过程，进展怎么样

你好啊，请问你最近也在用lammps计算介电常数吗，可以一起交流一下嘛

呵呵小怪 发表于 2022-6-29 22:33
**** 作者被禁止或删除 内容自动屏蔽 ****

我也困惑了半天，最后发现偶极矩张量是厄米的，所以只有上半部分。希望对你有帮助！

楼主现在关于节电常数的计算过程，进展怎么样

楼主现在进展怎么样

楼主现在处理出来介电常数了吗，可以交流一下吗

楼主现在处理出来介电常数了吗？

ZZU_SCU 发表于 2022-6-30 09:20
目前是这样的，偶极矩=（++）/3

非常感谢

呵呵小怪 发表于 2022-6-29 22:33
看来短时间是弄不完了，另外请教一下大佬对于原贴中提出的问题一有答案了吗，谢谢

目前是这样的，偶极矩=（<MxMx>+<MyMy>+<MzMz>）/3

ZZU_SCU 发表于 2022-6-29 10:36
还没有，停留在matlab或者python输入公式那一步。需要找个matlab工程师帮处理一下

看来短时间是弄不完了，另外请教一下大佬对于原贴中提出的问题一有答案了吗，谢谢

呵呵小怪 发表于 2022-6-29 00:40
请问解决了吗，我也找到了那篇文章，但是没法成功计算偶极矩

还没有，停留在matlab或者python输入公式那一步。需要找个matlab工程师帮处理一下

请问解决了吗，我也找到了那篇文章，但是没法成功计算偶极矩

有没有人可以帮我解答啊，不要沉啊。