vasp中采用杂化泛函计算能带结构求助。

didi_dudu

http://cms.mpi.univie.ac.at/wiki ... dified_KPOINTS_file
跟着教程在学习vasp能带结构的Hybrid calculation.有些东西没看懂，请教一下
教程中说在进行hybrid calculate时候，把之前第一步的IBZKPT拷过来，然后增加k-points的总数
Then add the desired additional k-points with zero weight and change the total number of k-points

Explicit k-points list
      18 <--- CHANGE TOTAL NUMBER OF K-POINTS !!
Reciprocal lattice
    0.00000000000000    0.00000000000000    0.00000000000000             1
    0.25000000000000    0.00000000000000    0.00000000000000             8
    0.50000000000000    0.00000000000000    0.00000000000000             4
    0.25000000000000    0.25000000000000    0.00000000000000             6
    0.50000000000000    0.25000000000000    0.00000000000000            24
   -0.25000000000000    0.25000000000000    0.00000000000000            12
    0.50000000000000    0.50000000000000    0.00000000000000             3
   -0.25000000000000    0.50000000000000    0.25000000000000             6
0.00000000 0.00000000 0.00000000 0.000 <--- ZERO WEIGHT !!
0.00000000 0.05555556 0.05555556 0.000
0.00000000 0.11111111 0.11111111 0.000
0.00000000 0.16666667 0.16666667 0.000
0.00000000 0.22222222 0.22222222 0.000
0.00000000 0.27777778 0.27777778 0.000
0.00000000 0.33333333 0.33333333 0.000
0.00000000 0.38888889 0.38888889 0.000
0.00000000 0.44444444 0.44444444 0.000
0.00000000 0.50000000 0.50000000 0.000

==============
zero weight之后这十个k点是新增加的。这十个点有什么意义呢？
另，GGA用于计算带隙是否会偏小？为了增加计算带隙宽度，可以采用什么方法呢？
（问题可能比较小白，求大神指导明路谢谢！（简单讲解，文章，或者其他链接资料，不胜感激））

卡开发发

0权重的k点在计算能量的时候这些k点的本征值会被计算不会被计入到总能当中啊（偷懒下，公式就不写了，你可以看vasp官网的lecture），否则，比如你计算n次Gamma的话，和你计算1次Gamma的总能不就不同了？以往纯泛函计算的时候也是同一个原理，只是非自洽过程中，这些0权重的点能够通过自洽时的电荷密度以非自洽的方式进行计算罢了，这个过程电荷密度是不更新的。杂化泛函没办法这样做罢了。

GGA算带隙得看体系，一般非半满、全满的d、f电子的系统可以DFT+U，U可以通过线性响应或者DMFT来确定，其他情形只好乖乖使用杂化泛函，或者GWA之类。

didi_dudu

谢谢！！本来是想问下细节的还是想想不做伸手党自己先查查才对就改了~

附上搜到的资料~

http://muchong.com/html/201211/5237955.html  

卡开发发

didi_dudu 发表于 2016-9-19 10:55
谢谢！！本来是想问下细节的还是想想不做伸手党自己先查查才对就改了~附上搜到的资料~http:/ ...

上次听别人说vasp的LDAUTYPE=3可以做线性响应，我不确定具体怎么做，PWscf应该很容易实现，自己写个简单脚本就行了。

didi_dudu

卡开发发 发表于 2016-9-19 12:01
上次听别人说vasp的LDAUTYPE=3可以做线性响应，我不确定具体怎么做，PWscf应该很容易实现，自己写个简单 ...

呃~再问个问题，看vasp教程里，做完band计算之后要画band 有一个提醒
Mind: Remove from the bandstructure plot the eigenvalues corresponding to the the regular k-points mesh.   
这是什么意思呢？是不是要删除EIGENVAL里面的某些值？在网上实在没搜到相应的具体介绍。如果知道的话，求解答！

卡开发发

didi_dudu 发表于 2016-9-29 11:00
呃~再问个问题，看vasp教程里，做完band计算之后要画band 有一个提醒
Mind: Remove from the bandstruct ...

应该是保留那些零权重的点，这些才是沿着高对称点的path，其余的点事Monkhorst-Pack生成的点，画能带用不到。

luzihen

卡开发发 发表于 2016-9-18 19:35
0权重的k点在计算能量的时候这些k点的本征值会被计算不会被计入到总能当中啊（偷懒下，公式就不写了，你可 ...

"这个过程电荷密度是不更新的。杂化泛函没办法这样做罢了。"为啥?不能直接用杂化范函自洽么?

卡开发发

luzihen 发表于 2016-11-11 17:50
"这个过程电荷密度是不更新的。杂化泛函没办法这样做罢了。"为啥?不能直接用杂化范函自洽么?

可以自洽，但是不可以非自洽，具体的算法我没研究过。