24小时热门版块排行榜    

查看: 14345  |  回复: 34
本帖产生 1 个 1ST强帖 ,点击这里进行查看

[交流] [悬赏交流]大家讨论一下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朵)

» 猜你喜欢

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

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

查看全部散金贴

已阅   回复此楼   关注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
31楼2013-07-02 05:38:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖
对于下面第一个问题悬赏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的回帖
对于下面这个问题悬赏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的回帖
★ ★ ★
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的关键还是太耗时
4楼2011-07-11 20:46:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
5楼2011-07-11 20:49:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
送鲜花一朵
引用回帖:

多谢!您用HSE计算了很多了,能否做个专题,给大家科普一下?多谢多谢!
6楼2011-07-11 20:55:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
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的回帖
引用回帖:
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点能快点 试了下 发现依然很慢啊 求高手指点啊
12楼2011-07-12 09:27:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by wenyan0517 at 2011-07-12 09:27:22:
最近 也在用杂化算结构性质 但慢的要命啊 看一篇文献里说减少K点能快点 试了下 发现依然很慢啊 求高手指点啊

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

uuv2010(金币+1):谢谢参与
vasp论坛上的资料
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的回帖
是不是参数有问题?能把您的INCAR贴上来看看吗?

» 本帖已获得的红花(最新10朵)

18楼2011-07-13 15:34:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
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的回帖
★ ★ ★
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 ]

» 本帖已获得的红花(最新10朵)

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

youzhizhe(金币+1): 谢谢交流。 2011-07-14 13:06:27
vasp论坛上的资料比较权威。
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的回帖
★ ★
uuv2010(金币+1):谢谢参与
youzhizhe(金币+1): 谢谢交流。请引用回复,这样方便知道回复来源。 2011-07-14 13:07:02
我用MS算过,, 就算是gamma 点, 也很慢,,,
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即可。
27楼2011-07-22 19:38:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)



uuv2010(金币+1):谢谢参与
送鲜花一朵
uuv2010(金币+1): 如有好的见解或者测试的例子,欢迎贴出来和大家一起讨论,谢谢! 2011-07-26 12:58:45
恩,很多人想了解下。。。恳请高手赐教
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):谢谢参与
2011-07-12 13:57   回复  
uuv2010(金币+1):谢谢参与
祝福祝福
wyh62120楼
2011-07-13 15:53   回复  
uuv2010(金币+1):谢谢参与
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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 276求调剂 +3 YNRYG 2026-03-21 4/200 2026-03-23 08:31 by 醉在风里
[考研] 384求调剂 +3 子系博 2026-03-22 4/200 2026-03-22 11:04 by 搏击518
[考研] 化学调剂 +5 yzysaa 2026-03-21 5/250 2026-03-21 22:12 by peike
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿华南师大 070300(化学)304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
信息提示
请填处理意见