24小时热门版块排行榜    

查看: 2215  |  回复: 26

aylixmut

铜虫 (小有名气)

[求助] 很奇怪的问题,vasp计算ENCUT设大就不让计算,高手帮忙啊!

第一次碰到这么奇怪的问题,在集群上算,ENCUT超过400,就提示错误如下:
Tue Jul 31 16:36:12 CST 2012
running on    8 nodes
distr:  one band on    1 nodes,    8 groups
vasp.5.2.12 26May11 complex
POSCAR found type information on POSCAR  Fe F
POSCAR found :  2 types and      24 ions
LDA part: xc-table for Pade appr. of Perdew
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...(           1 )
WAVECAR not read
WARNING: random wavefunctions but no delay for mixing, default for NELMDL
entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
mpiexec: Warning: tasks 0-7 died with signal 11 (Segmentation fault).
Tue Jul 31 16:36:22 CST 2012

下面是我INCAR设置文件:
SYSTEM = bulk MoF3
               
# ISTART = 0     ! start from scratch (WAVECAR)
#  GGA = 91
  VOSKOWN = 1    ! used when GGA = 91
# ICHARG =  1    ! rho read from CHGCAR.

  ENCUT =500  
  PREC = Accurate
  LREAL= Auto  
# NPAR = 1

# NBANDS = 10
# RWIGS  = 2.20   ! unit in A


# ISMEAR = -5   ! Tetrahedron method with Bloechel corrections.
  ISMEAR = 0    ! Tetrahedron method with Bloechel corrections.
  SIGMA = 0.2   ! in eV

#MAG related values:
  ISPIN =2
  MAGMOM = 2 -2 2 -2 2 -2 18*0     
  LORBIT = 11

#LDA+U
   LDAU =.TRUE.
   LDAUTYPE = 2
   LDAUL = 2 -1
   LDAUU = 9 0
   LDAUJ = 1 0
   LDAUPRINT = 2
#Ionic Relaxation
     ISIF = 3        ! Relax ions only, no change of cell shape and volume
     IBRION = 2      ! ionic relax: 0-MD; 1-quasi-New; 2-CG.
#    POTIM = 0.2     ! POTIM = 0.2
     POTIM = 0.5
#    IALGO  = 38     ! algorithm
     NSW = 30        ! number of ionic step
     EDIFFG = -0.01  ! all forces are smaller than 0.005, eV/A unit?
   ENCUT改为300即让算,但是结果误差比较大。到底是什么原因

[ Last edited by aylixmut on 2012-8-10 at 07:53 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylixmut

铜虫 (小有名气)

测试Sibulk(原胞2个原子时)就很正常,计算结果也很好,就是系统一大,比如32个原子,ENCUT只能很小算,一大就不让计算,提示错误如上,都是同样的位置停掉。 困扰好久了。
集群是学校管理员在管。
2楼2012-08-10 08:18:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

花无刺

木虫 (著名写手)

二师兄


WDD880227: 金币+1, 感谢帮助~~ 2012-08-13 19:02:32
先把自选轨道耦合关了试试

凭我的经验,结构弛豫可以不加自旋轨道耦合,影响不大的
加了,后面的计算反而会报错
这里user guide 完全是误导

PS:vasp5.1.12除了光棍节版,还有vasp.5.2.12 26May11 这个版本啊
3楼2012-08-10 08:21:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

花无刺

木虫 (著名写手)

二师兄

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
aylixmut: 金币+5, ★★★很有帮助 2012-08-10 10:01:10
liliangfang: 金币+1, 谢谢交流 2012-08-10 18:09:48
引用回帖:
2楼: Originally posted by aylixmut at 2012-08-10 08:18:49
测试Sibulk(原胞2个原子时)就很正常,计算结果也很好,就是系统一大,比如32个原子,ENCUT只能很小算,一大就不让计算,提示错误如上,都是同样的位置停掉。 困扰好久了。
集群是学校管理员在管。

回帖不能编辑的啊

请看看这个帖子:

http://muchong.com/bbs/viewthread.php?tid=4700360
4楼2012-08-10 08:31:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylixmut

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by 花无刺 at 2012-08-10 08:21:42
先把自选轨道耦合关了试试

凭我的经验,结构弛豫可以不加自旋轨道耦合,影响不大的
加了,后面的计算反而会报错
这里user guide 完全是误导

PS:vasp5.1.12除了光棍节版,还有vasp.5.2.12 26May11 这个版本 ...

我试过很多例子,有没算磁性的,都是一样报错哦。
5楼2012-08-10 10:01:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylixmut

铜虫 (小有名气)

引用回帖:
4楼: Originally posted by 花无刺 at 2012-08-10 08:31:20
回帖不能编辑的啊

请看看这个帖子:

http://muchong.com/bbs/viewthread.php?tid=4700360...

太好了,有人跟我一样的错误,只是他们提供的几种解决办法,我都试过了没用啊。
6楼2012-08-10 10:02:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cavediger

金虫 (正式写手)


liliangfang: 金币+1, 谢谢交流 2012-08-10 18:09:56
大概使用INTEL编译器搞的,几年前就碰到过,截断能大一点或K点多一点一进入主循环就退出来了,加-heap-arrays  64有可以解决的,也有不行的,而我就是后者,最后就改换编译器,用PGI重编一下就OK了
人比较懒,后来也没有再去研究什么原因了......
我们都老得太快,却聪明得太迟。当我一路跌跌撞撞登上山巅,蓦然发现,本来就有一条笔直的路,直达山顶....
7楼2012-08-10 17:01:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylixmut

铜虫 (小有名气)

引用回帖:
7楼: Originally posted by cavediger at 2012-08-10 17:01:12
大概使用INTEL编译器搞的,几年前就碰到过,截断能大一点或K点多一点一进入主循环就退出来了,加-heap-arrays  64有可以解决的,也有不行的,而我就是后者,最后就改换编译器,用PGI重编一下就OK了
人比较懒,后来 ...

谢谢回复,加 -heap-arrays 64 具体是怎么加的?
是加在makefile里,具体是哪个位置
,还是在命令行 make -heap-arrays 64?
谢谢!
8楼2012-08-10 22:27:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cavediger

金虫 (正式写手)


xueht987: 金币+1, 鼓励交流! 2012-08-12 10:26:46
引用回帖:
8楼: Originally posted by aylixmut at 2012-08-10 22:27:54
谢谢回复,加 -heap-arrays 64 具体是怎么加的?
是加在makefile里,具体是哪个位置
,还是在命令行 make -heap-arrays 64?
谢谢!...

FFLAGS =  -FR -lowercase -assume byterecl -heap-arrays 64
我们都老得太快,却聪明得太迟。当我一路跌跌撞撞登上山巅,蓦然发现,本来就有一条笔直的路,直达山顶....
9楼2012-08-11 16:59:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylixmut

铜虫 (小有名气)

引用回帖:
9楼: Originally posted by cavediger at 2012-08-11 16:59:37
FFLAGS =  -FR -lowercase -assume byterecl -heap-arrays 64...

谢谢,我试了也不行,看来我也是属于后者,有的折腾啊。谢谢
10楼2012-08-11 17:18:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 aylixmut 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +5 Ncdx123456 2026-03-13 6/300 2026-03-15 23:39 by lovewei0727
[考研] 0703化学调剂,求各位老师收留 +7 秋有木北 2026-03-14 7/350 2026-03-15 17:30 by 小物理化学
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 313分生物学求调剂 +6 Yyt杨1 2026-03-09 8/400 2026-03-14 03:00 by JourneyLucky
[考研] 306求调剂 +4 唐薏薏 2026-03-09 4/200 2026-03-14 01:19 by JourneyLucky
[考研] 材料工程专硕,一志愿中国矿业大学,总分314,求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 求调剂,一志愿江南大学环境工程085701 +3 Djdjj12 2026-03-10 4/200 2026-03-14 00:31 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 一志愿江南大学085701环境工程专硕总分287求调剂 +5 18266118446 2026-03-09 5/250 2026-03-11 16:51 by 2020015
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子!
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
信息提示
请填处理意见