24小时热门版块排行榜    

查看: 1329  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工 +5 安全上岸! 2026-03-16 5/250 2026-03-18 15:33 by cmz0325
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 一志愿西南交大,求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 材料专硕274一志愿陕西师范大学求调剂 +6 薛云鹏 2026-03-13 6/300 2026-03-18 14:14 by 脱颖而出
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 332求调剂 +6 Zz版 2026-03-13 6/300 2026-03-17 17:03 by ruiyingmiao
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
信息提示
请填处理意见