计算化学公社

标题: 求助：通过QE跑CPMD得到DOS [打印本页]

作者
Author: wth1219 时间: 2021-12-9 21:20
标题: 求助：通过QE跑CPMD得到DOS
各位老师好，在我个人比较熟悉的CP2K中，DOS和PDOS可以实现在跑FPMD的同时每隔N步直接输出；但对于我目前正在学习的QE，似乎只能通过projwfc.x后处理得到CPMD中一帧的DOS？

如果我开展长时间的CPMD，要计算几万帧的DOS，是不是要设置保存几万个波函数文件，然后挨个提交计算之后再求平均？感觉这即使编写了脚本降低了人力也有点过分啊……

不知道有没有其它我还没有发现抑或想到的方案呢？谢谢！

作者
Author: 卡开发发 时间: 2021-12-9 21:40
本帖最后由卡开发发于 2021-12-9 21:41 编辑

我能想到的办法是，输出的out每一步似乎能拿到能带占据的数目，然后自己想办法展宽，要不就改程序。qe在软件设计的角度来说，有些地方是不太好，模块过于零散。

作者
Author: wth1219 时间: 2021-12-9 21:48

卡开发发发表于 2021-12-9 21:40
我能想到的办法是，输出的out每一步似乎能拿到能带占据的数目，然后自己想办法展宽，要不就改程序。qe在软 ...

谢谢老师指教！
我对于QE完全是个初学者，改程序就不想了。（手动捂脸。）

我看到QE是可以输出wannier center的，不知道有没有办法从这个输出中计算得到DOS？

作者
Author: 卡开发发 时间: 2021-12-9 23:14

wth1219 发表于 2021-12-9 21:48
谢谢老师指教！
我对于QE完全是个初学者，改程序就不想了。（手动捂脸。）

似乎不大行，得有占据数的信息才行。另外CP那个模块我不大熟悉，有点吃不准。

作者
Author: wth1219 时间: 2021-12-10 00:25

卡开发发发表于 2021-12-9 23:14
似乎不大行，得有占据数的信息才行。另外CP那个模块我不大熟悉，有点吃不准。

好的，谢谢老师！

作者
Author: get-it 时间: 2021-12-11 06:12
最好先思考一下，真的有必要去计算几万帧的DOS吗？

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