24小时热门版块排行榜    

查看: 2632  |  回复: 21

fyuewen

金虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
sunyang1988: 金币+1, 谢谢交流 2012-07-19 23:09:34
引用回帖:
9楼: Originally posted by wujing110 at 2012-07-19 07:23:20
怎么会是一样的时候收敛啊??
那弛豫半天干嘛用的...

等弛豫好的时候你继续弛豫当然一样了。这就是你弛豫好的标准啊。。否则你怎么知道弛豫好了?
Physicsloveukulele!Aloha!
11楼2012-07-19 16:23:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wujing110

禁虫 (正式写手)

本帖内容被屏蔽

12楼2012-07-19 17:55:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

花无刺

木虫 (著名写手)

二师兄


sunyang1988: 金币+1, 谢谢交流 2012-07-19 23:08:38
1、ENCUT=520
520数字好听,其实400就足够了

2、ZnO的K点是不能取10 10 10的
你要是觉得不放心,取 24 24 13 试试

当然,要是你做了收敛性测试的话,没有必要取这么大的

13楼2012-07-19 19:03:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

argo

铁杆木虫 (著名写手)

上善若水

谢谢楼上各位的关注,我在三楼中所说的第一次和第二次是指INCAR中分别设成LREAL=Auto和LREAL=.FALSE.的计算。INCAR设置确实比较混乱,因为是刚开始学,这些参数都是别人告诉我怎么设我就怎么用,还没真正去理解为什么要这么设,所以,本帖的主要目的就是希望前辈们能够帮忙指点,应该如何准备这四个文件比较合适。
居善地,心善渊,与善仁,言善信,正善治,事善能,动善时。
14楼2012-07-20 08:41:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

argo

铁杆木虫 (著名写手)

上善若水

引用回帖:
13楼: Originally posted by 花无刺 at 2012-07-19 19:03:07
1、ENCUT=520
520数字好听,其实400就足够了

2、ZnO的K点是不能取10 10 10的
你要是觉得不放心,取 24 24 13 试试

当然,要是你做了收敛性测试的话,没有必要取这么大的


520也是别人这么设的,所以我跟着这么设而已
为什么k点不能取10 10 10,有原因吗?真心不知道,真心求教。谢谢
居善地,心善渊,与善仁,言善信,正善治,事善能,动善时。
15楼2012-07-20 08:43:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)


sunyang1988: 金币+1, 专家考核, 谢谢帮助 2012-07-20 18:24:48
引用回帖:
15楼: Originally posted by argo at 2012-07-20 08:43:07
520也是别人这么设的,所以我跟着这么设而已
为什么k点不能取10 10 10,有原因吗?真心不知道,真心求教。谢谢...

1, ZnO是c/a比较大的六方体系,Z轴方向的K点不需要取那么多~
2, 奇数的K点过原点,偶数的不过~估计你的LS是这个意思~不过这个可以通过设置gamma或者做k-shift来实现~
3, 10个左右的K点,对ZnO应该差不多也够了,不过对于驰豫嘛,还是建议多加些K点~不过ZnO是典型的半导体,K点不必非常密~
关于K点的取舍一是经验,二是要搞清楚K点是什么,为什么要取K点
思想重于技巧,内涵重于表象
16楼2012-07-20 11:27:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

argo

铁杆木虫 (著名写手)

上善若水

引用回帖:
16楼: Originally posted by souledge at 2012-07-20 11:27:15
1, ZnO是c/a比较大的六方体系,Z轴方向的K点不需要取那么多~
2, 奇数的K点过原点,偶数的不过~估计你的LS是这个意思~不过这个可以通过设置gamma或者做k-shift来实现~
3, 10个左右的K点,对ZnO应该差不多也够 ...

谢谢。已经设置gamma取点了
其实这些设置都是前面做过的同学帮我取的,我其实真的都还不清楚是在做什么,只是这样去做而已,但自己算出来的结果有点问题,所以应该是设置上什么地方不太合适了。其实个人认为不能这样什么都不懂就开始算,完全算是体力劳动,但……没办法
居善地,心善渊,与善仁,言善信,正善治,事善能,动善时。
17楼2012-07-20 11:54:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)


sunyang1988: 金币+1, 专家考核, 这个有点贵 2012-07-20 18:25:43
引用回帖:
17楼: Originally posted by argo at 2012-07-20 11:54:44
谢谢。已经设置gamma取点了
其实这些设置都是前面做过的同学帮我取的,我其实真的都还不清楚是在做什么,只是这样去做而已,但自己算出来的结果有点问题,所以应该是设置上什么地方不太合适了。其实个人认为不能这 ...

看看周围有没有人有MedeA,这个上手远快过直接用VASP,虽然计算也同样是VASP来做。不过至少可以学习一下基本参数的设置,比较直观一点~

计算后的结构如图:

计算的输入如下:
CODE:
    PREC = Accurate
   ENCUT = 520
  IBRION = 1
     NSW = 100
    ISIF = 3
    ALGO = Very Fast (RMM-DIIS)
    NELM = 100
  NELMIN = 2
   EDIFF = 1.0e-06
  EDIFFG = -0.01
VOSKOWN = 1
  NBLOCK = 1
   ISPIN = 1
  INIWAV = 1
  ISTART = 0
  ICHARG = 2
   LWAVE = .FALSE.
  LCHARG = .FALSE.
ADDGRID = .FALSE.
  ISMEAR = 1
   SIGMA = 0.1
   LREAL = .FALSE.
RWIGS = 0.73 1.25

CODE:
    1.0
            3.24270000 0.00000000 0.00000000
            -1.62134993 2.80826062 0.00000000
            0.00000010 0.00000016 5.19480000
  2 2
Direct
        0.33333333 0.66666667 0.38250000
        0.66666667 0.33333333 0.88250000
        0.33333333 0.66666667 0.00000000
        0.66666667 0.33333333 0.50000000

CODE:
Automatic mesh
0
Gamma
  15  15  11
  0.  0.  0.

思想重于技巧,内涵重于表象
18楼2012-07-20 12:30:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

argo

铁杆木虫 (著名写手)

上善若水

引用回帖:
18楼: Originally posted by souledge at 2012-07-20 12:30:49
看看周围有没有人有MedeA,这个上手远快过直接用VASP,虽然计算也同样是VASP来做。不过至少可以学习一下基本参数的设置,比较直观一点~
c1/41/941486_1342757825_373.png
计算后的结构如图:
f2/77/941486_134 ...

谢谢souledge,但回复 不是应助回复无法评分。

好象周边没人使用MedeA,或者他们有使用但我不知道。
居善地,心善渊,与善仁,言善信,正善治,事善能,动善时。
19楼2012-07-20 19:29:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

X.Wang

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
franch: 金币+2, 谢谢回帖交流。 2012-07-29 21:45:45
argo: 金币+5, 有帮助, 谢谢 2012-07-30 10:24:06
你可以试一试:
INCAR:
SYSTEM= ZnO
ENCUT=520
NSW = 60
IBRION = 2
ISIF  =  3
EDIFFG=-0.02
ISMEAR=0
SIGMA=0.1  (这个参数一般来说是默认值比较好,但是如果你测试过是0.1比较好就按照测试结果设置)
LWAVE=.FALSE.


POSCAR:
ZnO
1.0
         2.8148423098        -1.6251499653         0.0000000000
         0.0000000000        3.2502999306          0.0000000000
         0.0000000000         0.0000000000         5.2072000504
     2    2
Direct
      0.333333343         0.666666687         0.000000000
      0.666666627         0.333333313         0.500000000
      0.333333343         0.666666687         0.382900000
      0.666666627         0.333333313         0.882900000

KPOINTS:
Automatic mesh
0
Mon
6 6 4
0 0 0
(PS:K空间的建立是需要测试的但是比例应该基本如此,即与体系的基矢倒数比值保持一致)


POTCAR:
potcar的设置只要保证和原子顺序一致即可
在你的POSCAR设置的前提下为先Zn后O那么也就是说POTCAR中也是先Zn后O
20楼2012-07-29 20:45:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 argo 的主题更新
信息提示
请填处理意见