24小时热门版块排行榜    

查看: 1327  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 281求调剂(0805) +5 烟汐忆海 2026-03-16 13/650 2026-03-18 14:30 by stone_128
[考研] 0703化学调剂 ,六级已过,有科研经历 +10 曦熙兮 2026-03-15 10/500 2026-03-18 14:19 by 007_lilei
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +6 Liwangman 2026-03-15 6/300 2026-03-18 13:21 by 尽舜尧1
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 267一志愿南京工业大学0817化工求调剂 +6 SUICHILD 2026-03-12 6/300 2026-03-17 09:24 by 雾散后相遇lc
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
信息提示
请填处理意见