计算化学公社

标题: 介电常数和介电损耗计算求助 [打印本页]

作者
Author: ZZU_SCU 时间: 2022-4-19 22:07
标题: 介电常数和介电损耗计算求助
我用要lammps计算介电常数和介电损耗。第一步，通过https://mp.weixin.qq.com/s/0cRizbPrGSU8HO3A71Tgzw 帖子的

variable chunkID atom 1
compute mychunk all chunk/atom v_chunkID nchunk once ids once
compute dipole all dipole/chunk mychunk
variable dipole_x equal c_dipole[1][1]
variable dipole_y equal c_dipole[1][2]
variable dipole_z equal c_dipole[1][3]
fix 1 all ave/correlate 5 1000 10000 v_dipole_x v_dipole_y v_dipole_z type auto/upper ave running lx.txt

复制代码

计算时间偶极矩时间关联函数，输出文件的开头

# Time-correlated data for fix 4
# Timestep Number-of-time-windows
# Index TimeDelta Ncount v_dipole_x*v_dipole_x v_dipole_x*v_dipole_y v_dipole_x*v_dipole_z v_dipole_y*v_dipole_y v_dipole_y*v_dipole_z v_dipole_z*v_dipole_z
0 10000
1 0 1 2.79478e-10 4.35704e-10 1.58788e-05 6.79257e-10 2.47549e-05 0.902166
2 10 0 0.0 0.0 0.0 0.0 0.0 0.0
3 20 0 0.0 0.0 0.0 0.0 0.0 0.0

复制代码

我不会python处理时间偶极矩时间关联函数，参考https://zhuanlan.zhihu.com/p/391385899 热导率的python代码。准备修改用公式使用，公式来源于文献Nanoscale, 2021, 13, 10933 两个代码有很多共同点。
具体公式为
(, 下载次数 Times of downloads: 22)
我想请教偶极矩的矩阵和偶极矩的关系。

(, 下载次数 Times of downloads: 24)

问题一：
lammps的输出文件的 v_dipole_x*v_dipole_x v_dipole_x*v_dipole_y v_dipole_x*v_dipole_z v_dipole_y*v_dipole_y v_dipole_y*v_dipole_z v_dipole_z*v_dipole_z对应于
<MxMx><MxMy><MxMz>
<MyMy><MyMz>
<MzMz>么？
为什么没有<MyMx><MzMx>和<MzMy>
问题二：
某一帧结构文件的偶极矩等于（<MxMx>）**2+（<MyMy>）**2+（<MzMz>）**2再开根号么？

作者
Author: ZZU_SCU 时间: 2022-4-21 19:28
有没有人可以帮我解答啊，不要沉啊。

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

作者
Author: ZZU_SCU 时间: 2022-6-29 10:36

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

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

作者
Author: 呵呵小怪 时间: 2022-6-29 22:33

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

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

作者
Author: ZZU_SCU 时间: 2022-6-30 09:20

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

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

作者
Author: 呵呵小怪 时间: 2022-6-30 09:28

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

非常感谢

作者
Author: xkui 时间: 2022-11-7 17:53
楼主现在处理出来介电常数了吗？

作者
Author: sharklmh 时间: 2023-3-10 10:09
楼主现在处理出来介电常数了吗，可以交流一下吗

作者
Author: wildon 时间: 2023-5-11 16:27
楼主现在进展怎么样

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

作者
Author: suzhe 时间: 2023-11-15 15:34

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

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

作者
Author: suzhe 时间: 2023-11-17 14:08

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

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

作者
Author: huangyoyo 时间: 2024-4-7 13:40
想请教一下为什么我计算出来的偶极矩自相关联函数不是振荡的。
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。

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