介电常数和介电损耗计算求助

ZZU_SCU

我用要lammps计算介电常数和介电损耗。第一步，通过https://mp.weixin.qq.com/s/0cRizbPrGSU8HO3A71Tgzw 帖子的

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

复制代码

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

1. # Time-correlated data for fix 4
   
2. # Timestep Number-of-time-windows
   
3. # 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
   
4. 0 10000
   
5. 1 0 1 2.79478e-10 4.35704e-10 1.58788e-05 6.79257e-10 2.47549e-05 0.902166
   
6. 2 10 0 0.0 0.0 0.0 0.0 0.0 0.0
   
7. 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 两个代码有很多共同点。
具体公式为

我想请教偶极矩的矩阵和偶极矩的关系。



问题一：
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再开根号么？

ZZU_SCU

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

呵呵小怪

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

ZZU_SCU

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

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

呵呵小怪

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

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

ZZU_SCU

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

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

呵呵小怪

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

非常感谢

xkui

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

sharklmh

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

wildon

楼主现在进展怎么样

wyw2000

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

suzhe

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

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

suzhe

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

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

huangyoyo

想请教一下为什么我计算出来的偶极矩自相关联函数不是振荡的。
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。