应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 很奇怪的问题，vasp计算ENCUT设大就不让计算，高手帮忙啊！
271/3返回列表123下一页
查看: 2591  |  回复: 26
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-08-10 07:32:30
已阅   回复此楼   关注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 的主题更新
271/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 03:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 03:31 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[教师之家] 跳槽后在研项目怎么办？ +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] 什么时候开奖？ +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
信息提示
关闭
请填处理意见
关闭