24小时热门版块排行榜    

查看: 1482  |  回复: 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的回帖

happyjwx

木虫 (小有名气)


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

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的回帖

limaojlu

金虫 (正式写手)

我后来查了 这100个K点只是原来 14个在不可约布里渊区里面的 K点的在整个个(第一)布里渊区展开。不知道是不是K点选太多的原因,我们继续试一试。
5楼2010-01-08 09:23:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳变了,能看出什么问题? +7 基诺咪客 2026-08-17 10/500 2026-08-18 16:24 by Lullabygxb
[基金申请] 重要消息,中午系统在维护 +7 yuleib84 2026-08-18 8/400 2026-08-18 16:05 by jnhyjjm
[基金申请] 明天放榜? +4 Shxjjxjkx 2026-08-18 4/200 2026-08-18 15:59 by abcabc1133
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[教师之家] 两姐妹,一个辛辛苦苦30岁工作当大学教师,另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 5/250 2026-08-18 14:55 by LNP@mRNA
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见