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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP计算表面吸附 但是能量总是不收敛/。。。
101/1返回列表
查看: 4399  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

兜兜有颗糖

铜虫 (初入文坛)

[交流] 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多个原子，然后进行选择性弛豫，不过似乎服务器跑不动，所以就想要先做个小一点的，但是就是收敛不了
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-04-17 17:23:57
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 兜兜有颗糖 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 68gx6bh3ei 2026-07-09 3/150 2026-07-10 06:33 by pbuc7it0i8
[硕博家园] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 ph7bxbho8f 2026-07-09 3/150 2026-07-10 05:53 by j4frgf6zcw
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 al9t2f50x3 2026-07-09 4/200 2026-07-10 05:33 by j4frgf6zcw
[考研] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 al9t2f50x3 2026-07-09 4/200 2026-07-10 05:33 by j4frgf6zcw
[基金申请] b口会评 +6 当空照 2026-07-08 6/300 2026-07-09 20:49 by tonyliu0401
[基金申请] 信息学部会评时间 +3 xuel2011 2026-07-06 3/150 2026-07-09 19:05 by mollyzhang_2003
[基金申请] 祈祷青基必中 50+3 hadengry 2026-07-09 9/450 2026-07-09 17:37 by 文杰猪猪
[考研] 在职考研 +6 wu5d 2026-07-08 7/350 2026-07-09 16:40 by 化工小霸王呀
[论文投稿] MSER送审了还被拒稿 +3 S778950906 2026-07-08 5/250 2026-07-09 15:46 by S778950906
[基金申请] 关于如何从代码看上不上会 +11 且听虎啸 2026-07-09 14/700 2026-07-09 15:27 by libeny
[有机交流] 硝基还原 20+3 暮年飞雪 2026-07-07 5/250 2026-07-09 13:30 by tfang
[基金申请] 化学口会评？应用基础研究和基础研究，希望能分开 +3 Tide man 2026-07-08 3/150 2026-07-09 11:36 by dede5556
[基金申请] 今年E04面上 +15 布布和一二 2026-07-05 21/1050 2026-07-09 10:25 by 学员fOzRO9
[硕博家园] 大龄残疾硕士的一点执念 +14 3776737374 2026-07-06 22/1100 2026-07-09 06:07 by lovechizi
[基金申请] 生命口会评 +4 wuke100666 2026-07-07 4/200 2026-07-08 16:38 by 逍遥不够自在
[教师之家] 网络举报学术不端成为了男女之间情杀的手段了，真荒唐。。。 +3 zju2000 2026-07-04 3/150 2026-07-07 21:00 by hyx552200
[有机交流] Suzuki偶联中硼酸酯原料掉下的硼酸酯是否可以自身偶联，生成B2Pin2杂质 500+3 天下大同 2026-07-06 3/150 2026-07-07 10:31 by martyr_oyx
[论文投稿] 刑法学论文投稿求助 +3 高源678 2026-07-04 5/250 2026-07-06 12:20 by 平天尺
[论文投稿] Journal of Environmental Chemical Engineering +3 雪花天涯 2026-07-05 3/150 2026-07-06 09:00 by bobvan
[基金申请] 最后一年，祈求好运 +3 pyw1992 2026-07-06 3/150 2026-07-06 08:53 by 苏_frp
信息提示
关闭
请填处理意见
关闭