24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3933  |  回复: 9

兜兜有颗糖

铜虫 (初入文坛)

[交流] VASP计算表面吸附 但是能量总是不收敛/。。。已有5人参与

我要计算的体系是GaN的半极性面(1-10)
但是算的时候能量是正值,一共弛豫了99步,能量是正的170多电子伏

poscar如下:
test
1
3.2074 0.0000 0.0000
0.0000 5.2137 0.0000
0.0000 0.0000 14.4812
Ga N H
5 5 2
Direct
0.00000   0.38260   0.25575                          
0.50000   0.38260   0.06394                          
0.50000   0.38260   0.44756                          
0.00000   0.88260   0.38362                          
0.50000   0.88260   0.19181                          
0.00000   0.00000   0.25575                          
0.50000   0.00000   0.06394                          
0.50000   0.00000   0.44756                          
0.00000   0.50000   0.38362                          
0.50000   0.50000   0.19181                          
1.00000  -0.11740   0.00000                          
1.00000   0.50000   0.00000   


INCAR:
   SYSTEM = GaN
   Start parameter for this Run:
   #  LWAVE   = F
   #  LCHARG  = F      
   #  MAXMIX =  40
   #  ALGO = "Fast"
   #  ISYM = 2
      PREC   = Normal
      IALGO   = 48
   #  LPANE   = .TRUE.
   #  NPAR    = 4
   #  NSIM    = 4
   Electronic  Relaxation:

   #  NELMIN = 2
   #  NBLOCK = 1 ; KBLOCK =   1
   Ionic Relaxation:
      NSW = 99
   #  ISIF = 0 ! relax ions except changing volume
      ISIF = 2
   #  POTIM = 0.1
      IBRION = 2  
   DOS related values:
      ISMEAR  = -5
      ISTART = 0
      ICHARG = 2


另外 kpoints我用的是331  真空层的厚度8唉

这个是一个很小的体系,之前试过一个超大的体系,一共100多个原子,然后进行选择性弛豫,不过似乎服务器跑不动,所以就想要先做个小一点的,但是就是收敛不了
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

兜兜有颗糖

铜虫 (初入文坛)

我是新手,,,。,希望有大牛来交流下。。。
2楼2013-04-17 17:24:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengangelo

禁虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

3楼2013-04-17 17:43:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

兜兜有颗糖

铜虫 (初入文坛)

引用回帖:
3楼: Originally posted by fengangelo at 2013-04-17 17:43:44
你要不找个范本INCAR算一下

我很好奇的是。。。加了真空层之后INCAR要注意些什么,,,,
4楼2013-04-17 17:55:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengangelo

禁虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

5楼2013-04-17 18:50:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

兜兜有颗糖

铜虫 (初入文坛)

引用回帖:
5楼: Originally posted by fengangelo at 2013-04-17 18:50:58
我的日常INCAR, ISIF不考虑, IALGO  =    38
你要不要关注下 LDIPOL

NWRITE =      2
ISTART =      0
GGA= RP

Electronic Relaxation
ENCUT  = 400.0
EDIFF  = 0.1E-03
LREAL=.FALSE.
NELM=100
LDIP ...

恩恩 谢谢
6楼2013-04-17 23:32:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wgj501

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by fengangelo at 2013-04-17 18:50:58
我的日常INCAR, ISIF不考虑, IALGO  =    38
你要不要关注下 LDIPOL

NWRITE =      2
ISTART =      0
GGA= RP

Electronic Relaxation
ENCUT  = 400.0
EDIFF  = 0.1E-03
LREAL=.FALSE.
NELM=100
LDIP ...

请问下,这里DIPOL = 不需要指定吗?是不是默认的0.5,0.5,0.
5
沉下去!!!!
7楼2013-04-18 14:06:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yawhohaha

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
7楼: Originally posted by wgj501 at 2013-04-18 14:06:05
请问下,这里DIPOL = 不需要指定吗?是不是默认的0.5,0.5,0.
5...

要指定...
eg.

LDIPOL=.TRUE.
IDIPOL=3
DIPOL= 0.5 0.5 0.23 (center of thin film)
8楼2015-05-22 04:30:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

490105179

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
没设定POTIM如何弛豫?
9楼2015-05-22 09:16:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyan6930

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
固定下面几层原子

[ 发自小木虫客户端 ]
10楼2015-05-23 20:09:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 兜兜有颗糖 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见