计算化学公社

标题: 求助，Multiwfn如何根据cp2k输出信息计算偶极矩 [打印本页]

作者
Author: 云霭 时间: 2023-2-6 20:08
标题: 求助，Multiwfn如何根据cp2k输出信息计算偶极矩
大家好，由于cp2k中偶极矩和K点不能同时使用，我先使用K点跑MD，并输出HIRSHFELD电荷、盒子以及坐标文件。如果要得到偶极矩需要重跑一遍，请问Multiwfn能否直接根据已知文件信息计算相关偶极矩，如有提示，不胜感激，谢谢！

作者
Author: sobereva 时间: 2023-2-7 02:09
CP2K没法导出非gamma点以外的波函数信息到molden文件，Multiwfn也不支持k点。
我不知道你为什么在MD过程中考虑k点，有什么意义。必须考虑k点是晶胞较小的场合，而动力学会破坏原有晶体周期性，又不可能用小晶胞做动力学。极少有做动力学还考虑k点的。

作者
Author: 云霭 时间: 2023-2-7 07:57

sobereva 发表于 2023-2-7 02:09
CP2K没法导出非gamma点以外的波函数信息到molden文件，Multiwfn也不支持k点。
我不知道你为什么在MD过程中 ...

哦哦，谢谢老师指导！我只不过是看很多文献当中跑动力学都考虑了K点，所以觉得这个必不可少。

作者
Author: sobereva 时间: 2023-2-8 05:30

云霭发表于 2023-2-7 07:57
哦哦，谢谢老师指导！我只不过是看很多文献当中跑动力学都考虑了K点，所以觉得这个必不可少。

不要盲目效仿文献。文献里的体系特征和你可能有明显不同，导致适合的设置不同，而且文献里也未必用的就是对他们的体系而言真正合理的设置，此时更没借鉴意义

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