计算化学公社

标题: vasp中采用杂化泛函计算能带结构求助 [打印本页]

作者
Author: didi_dudu 时间: 2016-9-18 14:31
标题: vasp中采用杂化泛函计算能带结构求助
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用于计算带隙是否会偏小？为了增加计算带隙宽度，可以采用什么方法呢？
（问题可能比较小白，求大神指导明路谢谢！（简单讲解，文章，或者其他链接资料，不胜感激））

作者
Author: 卡开发发 时间: 2016-9-18 19:35
本帖最后由卡开发发于 2016-9-19 08:14 编辑

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

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

作者
Author: didi_dudu 时间: 2016-9-19 10:55
本帖最后由 didi_dudu 于 2016-9-19 11:01 编辑

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

附上搜到的资料~

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

作者
Author: 卡开发发 时间: 2016-9-19 12:01

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

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

作者
Author: didi_dudu 时间: 2016-9-29 11:00

卡开发发发表于 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里面的某些值？在网上实在没搜到相应的具体介绍。如果知道的话，求解答！

作者
Author: 卡开发发 时间: 2016-9-29 18:08

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

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

作者
Author: luzihen 时间: 2016-11-11 17:50

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

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

作者
Author: 卡开发发 时间: 2016-11-11 20:24

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

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

作者
Author: ljh123 时间: 2025-5-5 12:55

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

老师您好，我有一个100个原子的体系想算带隙。优化是用PBE泛函。我能否结构优化用PBE泛函，算带隙时用PBE+U呢？

作者
Author: 卡开发发 时间: 2025-5-5 19:09

ljh123 发表于 2025-5-5 12:55
老师您好，我有一个100个原子的体系想算带隙。优化是用PBE泛函。我能否结构优化用PBE泛函，算带隙时用PBE ...

操作上肯定可以这么做，但没有什么道理。

作者
Author: ljh123 时间: 2025-5-5 19:47

卡开发发发表于 2025-5-5 19:09
操作上肯定可以这么做，但没有什么道理。

作为材料计算新手，确实拿不准DFT+U。于是准备跟着这篇教程，改用mBJ方法来算带隙了https://github.com/chouyoudou/mBJ-VASP-tutorial

作者
Author: 卡开发发 时间: 2025-5-6 00:03

ljh123 发表于 2025-5-5 19:47
作为材料计算新手，确实拿不准DFT+U。于是准备跟着这篇教程，改用mBJ方法来算带隙了https://github.com/c ...

mBJ处理带隙确实还可以，不过就是太难收敛了。

作者
Author: ljh123 时间: 2025-5-6 12:06

卡开发发发表于 2025-5-6 00:03
mBJ处理带隙确实还可以，不过就是太难收敛了。

老师您好，关于这个教程最后一步：删除EIGENVAL中有权重的k点。请问下老师EIGENVAL中该删除哪些内容呢（我看不太懂这个文件）。受限于文件大小，EIGENVAL文件的部分内容如下：
  114  114 1 2
  0.2278146E+02  0.1157963E-08  0.1006487E-08  0.2341959E-08  0.5000000E-15
  1.000000000000000E-004
  CAR
unknown system
685    3 468

  0.2500000E+00  0.0000000E+00  0.0000000E+00  0.3333333E+00
   1       -19.546239    -19.531942 1.000000 1.000000
   2       -18.823599    -18.806610 1.000000 1.000000
   3       -18.791804    -18.767915 1.000000 1.000000
   4       -18.560393    -18.291891 1.000000 1.000000
....
299       2.479529       6.607787 1.000000 1.000000
300       2.498071       6.352898 1.000000 0.992424
301       2.530856       6.442620 1.000000 0.951995
302       2.566631       6.523340 1.000000 0.415183
303       2.579190       6.693675 1.000000 0.035421
304       2.595249       6.735825 1.000000 0.000371

作者
Author: 卡开发发 时间: 2025-5-6 12:23

ljh123 发表于 2025-5-6 12:06
老师您好，关于这个教程最后一步：删除EIGENVAL中有权重的k点。请问下老师EIGENVAL中该删除哪些内容呢（ ...

有权重的一般应该是靠前的几个，比如你给的例子中第八行有四个数的最后一个是权重，你可以把权重不为0这个部分带着这部分和下面能级信息的都去掉（我不确定他的教程是咋样的）。

作者
Author: ljh123 时间: 2025-5-6 12:34

卡开发发发表于 2025-5-6 12:23
有权重的一般应该是靠前的几个，比如你给的例子中第八行有四个数的最后一个是权重，你可以把权重不为0这 ...

老师，由于我文件里每个k点权重都是0.3333E+00，我就把所有k点以及下面的内容都删了。最后只剩下开头这点内容：
  114  114 1 2
  0.2278146E+02  0.1157963E-08  0.1006487E-08  0.2341959E-08  0.5000000E-15
  1.000000000000000E-004
  CAR
unknown system
685    3 468

这时候用vaspkit处理能带数据，会出现报错：
forrtl: severe (24): end-of-file during read, unit 11, file /path_to_file/EIGENVAL
Image             PC             Routine          Line       Source
vaspkit          00000000008FEFDB  Unknown             Unknown  Unknown
vaspkit          00000000009378D8  Unknown             Unknown  Unknown
vaspkit          00000000006F312C  read_eigenval_          61  readeigenval.F90
vaspkit          0000000000631D4D  get_band_                18  writeband.F90
vaspkit          00000000006EFC36  task_                   237  task.F90
vaspkit          00000000005EB092  mainmenu_                41  mainmenu.F90
vaspkit          00000000006B4E26  MAIN__                   27  main.F90
vaspkit          00000000004085A2  Unknown             Unknown  Unknown
libc-2.28.so    000014B1E663ACF3  __libc_start_main    Unknown  Unknown
vaspkit          00000000004084A9  Unknown             Unknown  Unknown

作者
Author: 卡开发发 时间: 2025-5-8 05:04

ljh123 发表于 2025-5-6 12:34
老师，由于我文件里每个k点权重都是0.3333E+00，我就把所有k点以及下面的内容都删了。最后只剩下开头这点 ...

我不太清楚vaspkit的操作流程，我这边都是自己写python脚本提取数据，你可能得好好看看他们教程。

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