高斯计算PBC结构的能带该怎么画出能带图呢

cqu_xiaoq · 发表于 Post on 2019-10-4 11:01:15

老师您好，首先我要道歉这个问题我在群里提问过，但是没有人回复，我下来有找了很久的资料却也没查到，所以我想发帖看看有没有没有注意到的，想在寻求下机会。问题如果：我用高斯计算了一个PBC结构的能带，K点选取了20个，计算完之后不知道该怎么来画出能带图，查阅了Multiwfn手册发现没有这一块的应用，又下查了，但是也没查到，所以抱着希望发帖在求助一下，虽然违规，但还望社长见谅。

kyuu · 发表于 Post on 2019-10-4 16:56:24

量子化学？综合交流？先告诉我你是怎样在PBC下用高斯选取20个K点的，不然根本就不知道该如何交流

biogon · 发表于 Post on 2019-10-4 17:42:49

kyuu 发表于 2019-10-4 16:56
量子化学？综合交流？先告诉我你是怎样在PBC下用高斯选取20个K点的，不然根本就不知道该如何交流

用NKPoint=N近似选择K点

cqu_xiaoq · 发表于 Post on 2019-10-4 17:48:09

kyuu 发表于 2019-10-4 16:56
量子化学？综合交流？先告诉我你是怎样在PBC下用高斯选取20个K点的，不然根本就不知道该如何交流

使用pbc下的选项pbc（NKpoint=20），此时就是选取了20个K点数，然后使用IOP(5/103=20)便会打印出相应的20个占据轨道和20个空轨道，在使用IOP(5/98=1)来保存所有的K点的特征值和轨道，问题是我算完之后log文件里的这些数据用啥程序来画图，我现在没法提取出来这些数据。还请指导一下

cqu_xiaoq · 发表于 Post on 2019-10-4 17:48:34

biogon 发表于 2019-10-4 17:42
用NKPoint=N近似选择K点

是这样，请问你做过类似的吗，你是怎么画的能带图呢？还请指导一下，谢谢。

卡开发发 · 发表于 Post on 2019-10-4 21:09:28

cqu_xiaoq 发表于 2019-10-4 17:48
使用pbc下的选项pbc（NKpoint=20），此时就是选取了20个K点数，然后使用IOP(5/103=20)便会打印出相应的20 ...

画能带用的是高对称点及其路径上的点，并不是给个均匀的k采样就行的。
如果有这些信息，即便gv不能做，要写个程序去画或者冒充其他程序输出应该不太难。

cqu_xiaoq · 发表于 Post on 2019-10-4 22:58:22

卡开发发发表于 2019-10-4 21:09
画能带用的是高对称点及其路径上的点，并不是给个均匀的k采样就行的。
如果有这些信息，即便gv不能做， ...

我功力不够，不懂程序，还有您讲的这个高对称的点及其路径上的点，这个我还不怎么懂，我下来在看看，还有别法子吗

卡开发发 · 发表于 Post on 2019-10-4 23:19:55

cqu_xiaoq 发表于 2019-10-4 22:58
我功力不够，不懂程序，还有您讲的这个高对称的点及其路径上的点，这个我还不怎么懂，我下来在看看，还有 ...

Gaussian没办法某个特定k点的计算方式（http://scuseria.rice.edu/gau/gdv_pbc/pbc_guide.html#band），虽然链接提到一种做法，但使用很麻烦，另外计算效率也比较低，PBC计算建议更换程序。

sobereva · 发表于 Post on 2019-10-5 04:48:15

Multiwfn没有对Gaussian的PBC计算画能带的功能
Gaussian官方PBC文档有绘制例子，即8L那个文档，但Gaussian的这个功能极不好用，而且那个文档里很多功能都是开发版或者开发者自己的code才有的，公开售卖版没法用
基本上Gaussian那个PBC功能只是作为泛函开发者的试验田，用来发理论方法的文章用的，实用性很低

biogon · 发表于 Post on 2019-10-5 15:27:00

sobereva 发表于 2019-10-5 04:48
Multiwfn没有对Gaussian的PBC计算画能带的功能
Gaussian官方PBC文档有绘制例子，即8L那个文档，但Gaussian ...

这个里面画出来的DOS图和multiwfn的有区别吗

卡开发发 · 发表于 Post on 2019-10-5 16:11:02

biogon 发表于 2019-10-5 15:27
这个里面画出来的DOS图和multiwfn的有区别吗

本质上没区别，但需要对不同的k态累加。

sobereva · 发表于 Post on 2019-10-6 07:07:53

biogon 发表于 2019-10-5 15:27
这个里面画出来的DOS图和multiwfn的有区别吗

Multiwfn是通过直接对轨道能级进行展宽得到的DOS，是近似模拟出来的；而PBC下计算的DOS是基于能带严格得到的，原理不同

kyuu · 发表于 Post on 2019-10-6 17:37:30

biogon 发表于 2019-10-4 17:42
用NKPoint=N近似选择K点

你能指定K点坐标？大神求带!

卡开发发 · 发表于 Post on 2019-10-6 17:49:36

kyuu 发表于 2019-10-6 17:37
你能指定K点坐标？大神求带!

选定的是均匀网格，而不是某个特定的k点。

kyuu · 发表于 Post on 2019-10-6 21:01:54

卡开发发发表于 2019-10-6 17:49
选定的是均匀网格，而不是某个特定的k点。

sarcasm

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 高斯计算PBC结构的能带该怎么画出能带图呢

浏览过的版块