24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1343  |  回复: 6
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 limaojlu 的 6 个金币

limaojlu

金虫 (正式写手)

[交流] 【求助】请教大家:利用杂化泛函计算固体电子结构的问题。

在VASP5.2中利用杂化泛函算氧化硅的静态计算, 按照手册的说明一般先是做正常的DFT计算然后再做,得到WAVECAR后再做杂化泛函计算。当我做杂化计算时总出错:提示是:
running on    8 nodes
distr:  one band on    1 nodes,    8 groups
vasp.5.2.2 15Apr09 complex
POSCAR found :  2 types and    9 ions
LDA part: xc-table for Pade appr. of Perdew
found WAVECAR, reading the header
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: small aliasing (wrap around) errors must be expected
FFT: planning ...(           1 )
reading WAVECAR
the WAVECAR file was read sucessfully
initial charge from wavefunction
entering main loop
       N       E                     dE             d eps       ncg     rms          ort
internal error in FOCK_ACC: number of k-points incorrect          14
         100           2           2           2
internal error in FOCK_ACC: number of k-points incorrect          14
         100           2           2           2
internal error in FOCK_ACC: number of k-points incorrect          14
         100           2           2           2
internal error in FOCK_ACC: number of k-points incorrect          14
         100           2           2           2
internal error in FOCK_ACC: number of k-points incorrect          14
         100           2           2           2
internal error in FOCK_ACC: number of k-points incorrect          14
         100           2           2           2
internal error in FOCK_ACC: number of k-points incorrect          14
         100           2           2           2
internal error in FOCK_ACC: number of k-points incorrect          14
         100           2           2           2

从提示看是不是我的K网格选择有问题,我选用的是5x5x4的MP取点,这个K点的选取在做正常DFT计算中,没有问题的。这里14应该是指不可约的K点数目,我查看了OUTCAR, 但是用于杂化泛函处理的K点有100个。

大家 有没有做过类似计算的经验啊,谢谢指教啊。

下面是我的INCAR和KPOINTS

....................................
SYSTEM = SiO2-Primitive Cell

Start parameter for this run:
ISTART = 1
ICHARG = 2
NWRITE = 2
ISPIN = 2

Electronic Relaxation
ENCUT = 540
NELM = 100
EDIFF = 1E-05
LREAL = .FALSE.

Ionic Relaxation
EDIFFG = -0.03   Stopping-Criterion for IOM (HF-FORCE)
NSW = 0
IBRION = -1
ISIF = 2
POTIM = 0.5

LHFCALC = .TRUE.
HFSCREEN = 0.2
NBANDS = 40
ALGO = All
TIME = 0.4
ENCUTFOCK = 0
NKRED = 2


DOS related values:
ISMEAR = 0
SIGMA = 0.1

Write flags
LWAVE = .TRUE.
LCHARG = .TRUE.
...............................................
Auto
0
M-K
5 5 4
0 0 0
回复此楼

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

杂化泛函方面的研究

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

etree

铜虫 (初入文坛)

★ ★
limaojlu(金币+1,VIP+0): 1-6 09:53
ice_rain(金币+1,VIP+0):谢谢交流 1-6 12:00
可能跟编译的参数有关
2楼2010-01-06 01:42:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylayl08

荣誉版主 (文坛精英)

至尊天蝎

优秀版主

★ ★
limaojlu(金币+2,VIP+0): 1-8 09:20
引用回帖:
Originally posted by etree at 2010-1-6 01:42:
可能跟编译的参数有关

能否给个具体意见呢?哪里的参数出了问题?
珍惜拥有,快乐人生!欢迎光临计算模拟区(*^__^*)。。。。。。。。。。。。。。
3楼2010-01-07 15:50:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

happyjwx

木虫 (小有名气)


limaojlu(金币+1,VIP+0): 1-8 09:21
没用过5.2的新功能,一直当4.6用呢
跟贴学习一下
夜半追忆往昔
4楼2010-01-07 18:50:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

limaojlu

金虫 (正式写手)

我后来查了 这100个K点只是原来 14个在不可约布里渊区里面的 K点的在整个个(第一)布里渊区展开。不知道是不是K点选太多的原因,我们继续试一试。
5楼2010-01-08 09:23:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

etree

铜虫 (初入文坛)

贴一下你的makefile文件,已经用的库的信息。
引用回帖:
Originally posted by aylayl08 at 2010-1-7 15:50:

能否给个具体意见呢?哪里的参数出了问题?

6楼2010-01-09 00:31:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

limaojlu

金虫 (正式写手)

★ ★
wuchenwf(金币+2,VIP+0):谢谢 1-9 22:20
解决了,去掉ENCUTFOCK = 0 NKRED = 2,这两个参数就好了,可能是因为NKRED这个对屏蔽HF交换势减少展开项的比例可能在我的六方氧化硅胞中,不合适。
7楼2010-01-09 16:02:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 limaojlu 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085501机械英二77总分294求调剂,接受跨专业学习 +5 守法公民亓纪 2026-04-08 5/250 2026-04-08 21:19 by bljnqdcc
[考研] 264求调剂 +11 麦小叮当 2026-04-07 11/550 2026-04-08 16:05 by 一只好果子?
[考研] 274求调剂求调剂 +10 Jachenbingoo 2026-04-06 13/650 2026-04-08 14:25 by zhq0425
[考研] 求材料调剂,一志愿郑州大学289分 +21 硕星赴 2026-04-03 21/1050 2026-04-08 11:55 by 猪会飞
[考研] 285求调剂 +7 AZMK 2026-04-07 7/350 2026-04-08 11:10 by 逆水乘风
[考研] 0703总分331求调剂 +17 ZY-05 2026-04-04 21/1050 2026-04-08 10:16 by screening
[考研] 326分,一志愿沪9,求生物学调剂 +4 刘墨墨 2026-04-05 4/200 2026-04-08 06:22 by lijunpoly
[考研] 305求调剂 +3 77Qi 2026-04-07 3/150 2026-04-07 23:18 by jp9609
[考研] 生物调剂 +5 橙子橙子橙子啊 2026-04-05 9/450 2026-04-07 15:31 by 上岸快快
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
[考研] 材料专硕283求调剂 +17 试试看呗 2026-04-04 18/900 2026-04-06 09:24 by 286640313
[考研] 085600,320分求调剂 +16 大馋小子 2026-04-04 17/850 2026-04-06 07:58 by MOF_Catal
[考研] 0703求调剂383分 +9 W55j 2026-04-03 9/450 2026-04-06 06:50 by houyaoxu
[考研] 男生,一志愿沪9生物学071000,初试308求调剂 +3 刘墨墨 2026-04-04 3/150 2026-04-05 08:26 by barlinike
[考研] [调剂信息]085408光电信息 求调剂 总分291分数一英一 +3 iz11az 2026-04-02 3/150 2026-04-04 19:09 by 蓝云思雨
[考研] 求调剂 +3 ffyyu 2026-04-02 3/150 2026-04-04 19:03 by 蓝云思雨
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 14:46 by 湘农储能材料
[考研] 考研调剂 +3 李木子0120 2026-04-02 5/250 2026-04-02 21:45 by dongzh2009
[考研] 一志愿大工学硕,求调剂 +4 yub0811 2026-04-02 4/200 2026-04-02 21:36 by 百灵童888
信息提示
请填处理意见