应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » [悬赏交流]大家讨论一下vasp中HF以及杂化泛函的应用吧
351/1返回列表
查看: 13761  |  回复: 34
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 1ST强帖 ，点击这里进行查看

uuv2010

荣誉版主 (职业作家)

[交流] [悬赏交流]大家讨论一下vasp中HF以及杂化泛函的应用吧

vasp手册上关于HF的内容链接如下：
http://cms.mpi.univie.ac.at/vasp/vasp/HF_type_calculations.html
大家一起讨论一下HF以及杂化泛函在VASP中的应用吧。
主要的参数包括
LHFCALC= .TRUE. | .FALSE.    （这个参数是必须的）

AEXX = [real] (fraction of exact exchange)
ALDAC= [real] (fraction of LDA correlation energy)
AGGAX= [real] (fraction of gradient correction to exchange)
AGGAC= [real] (fraction of gradient correction to correlation)

ENCUTFOCK= [real]
PRECFOCK= Low | Medium | Fast | Normal | Accurate
LMAXFOCK= [integer]                           
LMAXFOCKAE= [integer] (maximum L quantum number for accurate charge augmentation in HF routines)

HFSCREEN= [real]
LTHOMAS= .TRUE. | .FALSE.

NKRED= [integer]

NKREDX= [integer]  NKREDY= [integer]   NKREDZ= [integer]

EVENONLY= [logical]  ODDONLY= [logical]

大家有应用vasp中提供的HF和杂化泛函计算的好的例子，以及对上述参数的理解，欢迎讨论~

[ Last edited by uuv2010 on 2011-7-11 at 15:50 ]
回复此楼

» 收录本帖的淘帖专辑推荐

VASP first principle 表界面计算 vasp光学
HSE & 杂化函数 杂化泛函方面的研究 vasp VASP blabla
第一性原理 DFT+U，HSE等 个人收藏

» 本帖已获得的红花（最新10朵）

enola 骆驼忍着

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼2011-07-11 13:44:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

亚历山大王子

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵
引用回帖:
18楼: Originally posted by uuv2010 at 2011-07-13 15:34:23
是不是参数有问题？能把您的INCAR贴上来看看吗？

ISPIN=2
PREC=medium  
NSW=200
ISIF=3
ISMEAR=0
SIGMA=0.2
NEDOS=2000
这是我计算的INCAR文件，但是总无法继续进行计算，不知道问题出在哪里，请多多指教。
IBRION=2
NELM=200
ISTART=1
LHFCALC=.TURE.
HFSCREEN=0.2
ALGO=ALL
TIME=0.4
PRECFOCK=Fast
NKRED=2     
ISYM=0
EDIFFG=0.00001
ENCUT=200
赞一下(1人) 回复此楼
31楼2013-07-02 05:38:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

uuv2010

荣誉版主 (职业作家)
对于下面第一个问题悬赏50金币，给予大家认为最合理的答案提供者，欢迎根据自己的理解进行讨论


先提一个问题供大家讨论：
问题一：
对于对于d电子的体系，LMAXFOCK应该怎样设置？
手册上的说明是
LMAXFOCK  sets the maximum angular quantum number L for the augmentation of charge densities in Hartree-Fock type routines.
是不是对于d电子的体系需要LMAXFOCK =2 ？

[ Last edited by uuv2010 on 2011-7-22 at 21:47 ]
赞一下 回复此楼
2楼2011-07-11 13:50:03
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)
对于下面这个问题悬赏50金币，给予大家认为最合理的答案提供者，欢迎根据自己的理解进行讨论

在vasp手册中提供了一个例子

For metals and small gap semiconductors it is recommended to use.
ISTART = 1
LHFCALC = .TRUE. ; HFSCREEN = 0.2
ALGO = Damped ; TIME = 0.4
ENCUTFOCK = 0  ! omit flag for high quality calculations
NKRED     = 2  ! omit flag for high quality calculations

不知道这里为什么ENCUTFOCK = 0， NKRED     = 2 ？为什么不需要FOCK部分？

[ Last edited by uuv2010 on 2011-7-22 at 19:32 ]
赞一下 回复此楼
3楼2011-07-11 13:52:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)
★ ★ ★
uuv2010(金币+1):谢谢参与
uuv2010(金币+5): 是的，大体系的确很耗时，小体系还可以试一试。多谢提示！ 2011-07-11 20:51:25
zzy870720z(金币+2): 老前辈出马，一个顶俩，哈哈 2011-07-12 08:16:20
关键是k点的设置。
计算dos和能带结构时候的k点的设置。
如果k点小了dos图很难看，如果k点大了很难算。
HSE的关键还是太耗时
赞一下(2人) 回复此楼
4楼2011-07-11 20:46:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)

uuv2010(金币+5): 多谢提示 2011-07-11 20:56:12
http://muchong.com/bbs/viewthread.php?tid=2945479

http://muchong.com/bbs/viewthread.php?tid=2654952

http://muchong.com/bbs/viewthread.php?tid=2624954

http://muchong.com/bbs/viewthread.php?tid=2605650
这些全是没有解决的问题。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

uuv2010
5楼2011-07-11 20:49:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)
送鲜花一朵
引用回帖:

多谢！您用HSE计算了很多了，能否做个专题，给大家科普一下？多谢多谢！
赞一下 回复此楼
6楼2011-07-11 20:55:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)
引用回帖:
Originally posted by uuv2010 at 2011-07-11 20:55:21:
多谢！您用HSE计算了很多了，能否做个专题，给大家科普一下？多谢多谢！

我也是一知半解。
害怕误导大家。
http://cms.mpi.univie.ac.at/vasp-forum/forum.php
vasp论坛上的资料比较权威。
大家可以看看。
赞一下 回复此楼
7楼2011-07-11 20:59:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)
引用回帖:
Originally posted by xuzhl at 2011-07-11 21:05:42:

希望您能就具体的问题进行讨论，免得帖子水掉。非常感谢您的参与！
赞一下 回复此楼
9楼2011-07-11 21:08:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuliuxuan

至尊木虫 (职业作家)
★
uuv2010(金币+1):谢谢参与
多谢
回复此楼
10楼2011-07-11 21:16:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuzhl

木虫 (著名写手)
引用回帖:
Originally posted by uuv2010 at 2011-07-11 21:08:36:
希望您能就具体的问题进行讨论，免得帖子水掉。非常感谢您的参与！

这个真不会，隔行如隔山啊。。。只能祝福了
赞一下 回复此楼
11楼2011-07-11 21:21:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenyan0517

新虫 (初入文坛)
★
uuv2010(金币+1):谢谢参与
最近 也在用杂化算结构性质 但慢的要命啊 看一篇文献里说减少K点能快点 试了下 发现依然很慢啊 求高手指点啊
赞一下(1人) 回复此楼
12楼2011-07-12 09:27:22
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)
引用回帖:
Originally posted by wenyan0517 at 2011-07-12 09:27:22:
最近 也在用杂化算结构性质 但慢的要命啊 看一篇文献里说减少K点能快点 试了下 发现依然很慢啊 求高手指点啊

的确很慢，看来只能计算小体系了。
赞一下 回复此楼
13楼2011-07-12 13:10:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leimiao_hit

木虫之王 (文学泰斗)
★
uuv2010(金币+1):谢谢参与
vasp论坛上的资料
赞一下(1人) 回复此楼
15楼2011-07-12 13:30:07
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenyan0517

新虫 (初入文坛)
引用回帖:
Originally posted by uuv2010 at 2011-07-12 13:10:40:
的确很慢，看来只能计算小体系了。

我的体系不算大啊 单胞里就14个原子 唉 用杂化计算的伤不起啊
赞一下 回复此楼
17楼2011-07-13 09:15:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)
是不是参数有问题？能把您的INCAR贴上来看看吗？
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

亚历山大王子
18楼2011-07-13 15:34:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leimiao_hit

木虫之王 (文学泰斗)
ALDAC= [real] (fraction of LDA correlation energy)
AGGAX= [real] (fraction of gradient correction to exchange)
AGGAC= [real] (fraction of gradient correction to correlation)
赞一下 回复此楼
19楼2011-07-13 15:42:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)
★ ★ ★
youzhizhe(金币+3, 1ST强帖+1): 谢谢交流。 2011-07-14 13:06:19
引用回帖:

对于杂化泛函计算band的的撒点问题，即如下帖子所述的其中的一个：
http://muchong.com/bbs/viewthread.php?tid=2654952
就我自己粗浅的理解谈谈：
- Add the points at the end of the KPOINTS file, but set the weights of these added k-points to 0.
-the number of the k-points used in the standard mesh PLUS the number of the k-points along the lines
-NKRED can not be used.
上述三点是最重要的，当然这样需要自己产生band所需要的一系列k点，并且会给出如下的警告信息：
| Your generating k-point grid is not commensurate to the symmetry |
| of the lattice. This can cause slow convergence with respect | | to k-points for HF type calculations |
| suggested SOLUTIONS: |
| ) if not already the case, use automatic k-point generation |
| ) shift your grid to Gamma (G) (e.g. required for hex or fcc lattice) |
这个信息中红色标记的是重要的提示，就是说这样的设置对于HF类型的计算是正确的，而对于其他的比如纯DFT的计算这样的设置是有问题的。
所以这个警告信息并不用担心。
这种方法做出的band需要将前面那些有权重的K点去掉，然后才能做出正确的band图。

[ Last edited by uuv2010 on 2011-7-14 at 09:04 ]
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

gaozhibin
22楼2011-07-13 18:32:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leimiao_hit

木虫之王 (文学泰斗)
★
youzhizhe(金币+1): 谢谢交流。 2011-07-14 13:06:27
vasp论坛上的资料比较权威。
赞一下(1人) 回复此楼
23楼2011-07-13 18:50:22
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyzn

铁杆木虫 (知名作家)
★
uuv2010(金币+1):谢谢参与
看看
回复此楼
24楼2011-07-13 18:57:07
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljq1988114

木虫 (正式写手)
★
uuv2010(金币+1):谢谢参与
困惑
回复此楼
25楼2011-07-14 10:49:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

franch

版主 (著名写手)
★ ★
uuv2010(金币+1):谢谢参与
youzhizhe(金币+1): 谢谢交流。请引用回复，这样方便知道回复来源。 2011-07-14 13:07:02
我用MS算过，， 就算是gamma 点， 也很慢，，，
赞一下(2人) 回复此楼
26楼2011-07-14 11:05:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

caihe605163

新虫 (初入文坛)
★
uuv2010(金币+1):谢谢参与
uuv2010(金币+50): 欢迎参加交流 2011-07-22 19:38:36
引用回帖:
Originally posted by uuv2010 at 2011-07-11 13:52:02:
对于下面这个问题悬赏50金币，给予大家认为最合理的答案提供者，欢迎根据自己的理解进行讨论

在vasp手册中提供了一个例子

For metals and small gap semiconductors it is recommended  ...

这是简单的HSE06的计算，ENCUTFOCK = 0没有考虑FOCK部分。为了考虑FOCK部分，在vasp5.2中高精度的计算设这ENCUTFOCK =F即可。
赞一下(1人) 回复此楼
27楼2011-07-22 19:38:03
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)
★
uuv2010(金币+1):谢谢参与
送鲜花一朵
uuv2010(金币+1): 如有好的见解或者测试的例子，欢迎贴出来和大家一起讨论，谢谢！ 2011-07-26 12:58:45
恩，很多人想了解下。。。恳请高手赐教
赞一下(1人) 回复此楼
28楼2011-07-25 21:46:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
对这个NKRED的值不太清楚, 手册的例子上给出的是NKRED=2。莫非大家都用这个值？
赞一下 回复此楼
33楼2014-09-25 11:35:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

骆驼忍着

铜虫 (小有名气)
送红花一朵
收益多多
回复此楼
34楼2015-01-06 16:11:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaozhibin

铜虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵
引用回帖:
22楼: Originally posted by uuv2010 at 2011-07-13 18:32:26
对于杂化泛函计算band的的撒点问题，即如下帖子所述的其中的一个：
http://muchong.com/bbs/viewthread.php?tid=2654952
就我自己粗浅的理解谈谈：
- Add the points at the end of the KPOINTS file, but s ...

你说的对，只有红色的字体： This can cause slow convergence with respect | | to k-points for HF type calculations |  才说明我们打开了 HF 计算。

学习了！！
赞一下 回复此楼
35楼2015-11-22 16:24:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
xuzhl8楼
2011-07-11 21:05   回复  
uuv2010(金币+1):谢谢参与
iamacr14楼
2011-07-12 13:24   回复  
uuv2010(金币+1):谢谢参与
wangwenqiang16楼
2011-07-12 13:57   回复  
uuv2010(金币+1):谢谢参与
祝福祝福
wyh62120楼
2011-07-13 15:53   回复  
uuv2010(金币+1):谢谢参与
太白药谷21楼
2011-07-13 15:54   回复  
uuv2010(金币+1):谢谢参与
ZHUHAPPYRR29楼
2013-04-23 10:14   回复  
ZHUHAPPYRR30楼
2013-05-07 08:35   回复  
cskx_hxj32楼
2013-09-03 10:55   回复  
顶一下
相关版块跳转 我要订阅楼主 uuv2010 的主题更新
351/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭