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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】我对表面驰豫,算完了原子却一个都没有动
101/1返回列表
查看: 2043  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

liuhengjin

银虫 (小有名气)

[交流] 【求助】我对表面驰豫,算完了原子却一个都没有动

这是我的INCAR文件

SYSTEM = Ge(100)_(2x1)
ENCUT =400
ISMEAR = 0; SIGMA = 0.05
ISTART = 0; ICHARG = 2
GGA=PE
ALGO =Fast
EDIFF = 1E-5; EDIFFG = -1.0E-2
NELMIN=6
ISIF=3
IBRION=2
LREAL=AUTO

POSCAR为
Ge(110)_(2x1)                           
   1.02200000000000     
     4.0004566145628928    0.0000000000000000    0.0000000000000000
     0.0000000000000000    8.0009132291257874    0.0000000000000000
     0.0000000000000000    0.0000000000000000   29.1437500000000007
  22
Selective dynamics
Direct
  0.0000000000000000  0.0000000000000000  0.0970619772679981   F   F   F
  0.0000000000000000  0.0019822828384809  0.2967260685518983   T   T   T
  0.0000000000000000  0.9276824881117479  0.4932595612465049   T   T   T
  0.5000000000000000  0.2500000000000000  0.0000000000000000   F   F   F
  0.5000000000000000  0.2500000000000000  0.1971038563673539   T   T   T
  0.5000000000000000  0.2500000000000000  0.3964169710606655   T   T   T
  0.5000000000000000  0.0000000000000000  0.0485309886339991   F   F   F
  0.5000000000000000 -0.0029335760595449  0.2466868371805478   T   T   T
  0.5000000000000000  0.0058225771951400  0.4521603543234448   T   T   T
  0.0000000000000000  0.2500000000000000  0.1455929659000006   F   F   F
  0.0000000000000000  0.2500000000000000  0.3481028386587142   T   T   T
  0.0000000000000000  0.5000000000000000  0.0970619772679981   F   F   F
  0.0000000000000000  0.4980177171615190  0.2967260685518983   T   T   T
  0.0000000000000000  0.5723175118862555  0.4932595612465049   T   T   T
  0.5000000000000000  0.7500000000000000  0.0000000000000000   F   F   F
  0.5000000000000000  0.7500000000000000  0.1964311093813827   T   T   T
  0.5000000000000000  0.7500000000000000  0.3954745052361157   T   T   T
  0.5000000000000000  0.5000000000000000  0.0485309886339991   F   F   F
  0.5000000000000000  0.5029335760595451  0.2466868371805478   T   T   T
  0.5000000000000000  0.4941774228048601  0.4521603543234448   T   T   T
  0.0000000000000000  0.7500000000000000  0.1455929659000006   F   F   F
  0.0000000000000000  0.7500000000000000  0.3478604452255454   T   T   T



请大家指导，应该怎么设置?
我要驰豫的表面是Ge(100)_(2x1)
谢谢了

[ Last edited by liuhengjin on 2010-11-21 at 19:43 ]
回复此楼

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2010-11-21 19:42:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)
★ ★
mazuju028(金币+2):谢谢提示 2010-11-21 22:36:43
liuhengjin(金币+1):谢谢提醒 2010-11-22 13:58:28
你的NSW打开了吗？默认是NSW＝0，也就是离子步为0，所以肯定是计算一步之后立即停止计算，所以离子不会发生移动
赞一下(1人) 回复此楼
2楼2010-11-21 19:53:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)
★ ★
mazuju028(金币+2):谢谢交流 2010-11-21 22:37:46
表面优化你还用ISIF=3 难道你有优化晶胞的体积和形状？


本文来自: 小木虫论坛 http://muchong.com/bbs/viewthread.php?tid=2618313&fpage=1
赞一下(1人) 回复此楼
3楼2010-11-21 20:27:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhengjin

银虫 (小有名气)
2楼是对的，不过我还想就3楼的疑问再问下，表面优化不可以设置isif=3吗？这个应该比isif=2更好吧？

[ Last edited by liuhengjin on 2010-11-22 at 13:57 ]
赞一下 回复此楼
4楼2010-11-22 13:46:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢指教 2010-11-22 18:31:17
表面优化ISIF=2合理一些吧。
我是从全驰豫的体材料（ISIF=3)构建表面的。也可以这样理解：表面是个二维结构，只有一个方向不具有周期性，其它两个方向（最下面几层）都还是和体材料的环境一样的，它们还应该和体材料的晶格常数一致，所以驰豫的时候就ISIF=2了。

[ Last edited by stractor on 2010-11-22 at 15:12 ]
赞一下(2人) 回复此楼
5楼2010-11-22 15:07:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xianggui7895

银虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢指教 2010-11-22 18:31:29
一般先ISIF=3优化体材料单胞，而用ISIF=2优化表面，构建表面的单元是体材料的优化结果。
赞一下(2人) 回复此楼
6楼2010-11-22 16:17:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhengjin

银虫 (小有名气)
感谢各位的回答，谢谢。不过就6楼的观点，我个人认为，不应该先用ISIF=3来全驰豫，再用ISIF=2，这样会花掉很多时间，因为做全驰豫时底下的原子还是没有变，只是上层的原子改变了

[ Last edited by liuhengjin on 2010-11-23 at 13:39 ]
赞一下 回复此楼
7楼2010-11-23 13:26:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ustbwgs

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by liuhengjin at 2010-11-21 19:42:14:
这是我的INCAR文件

SYSTEM = Ge(100)_(2x1)
ENCUT =400
ISMEAR = 0; SIGMA = 0.05
ISTART = 0; ICHARG = 2
GGA=PE
ALGO =Fast
EDIFF = 1E-5; EDIFFG = -1.0E-2
NELMIN=6
ISIF=3
IBRION=2
LREAL=AUTO ...

为什么要每隔两个原子设置一个F  F  F呢
赞一下(1人) 回复此楼
8楼2010-11-23 14:38:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
sunyang1988(金币+1):谢谢回复 2010-11-23 16:27:23
本帖仅楼主可见
赞一下
9楼2010-11-23 16:15:29
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

xianggui7895

银虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
sunyang1988(金币+1):谢谢交流 2010-11-27 11:28:18
引用回帖:
Originally posted by liuhengjin at 2010-11-23 13:26:40:
感谢各位的回答，谢谢。不过就6楼的观点，我个人认为，不应该先用ISIF=3来全驰豫，再用ISIF=2，这样会花掉很多时间，因为做全驰豫时底下的原子还是没有变，只是上层的原子改变了

[ Last edited by liuhengjin ...

你可能误解我的意思了，全弛豫的是体材料单胞，不是表面模型。
赞一下(2人) 回复此楼
10楼2010-11-27 10:23:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liuhengjin 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 12:20 by n0P4fnsXKbBK
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 YXM1FDCSE8Bn 2026-10-31 5/250 2026-11-02 09:42 by dcc2yOIRGGTT
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 09:29 by dcc2yOIRGGTT
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 09:29 by dcc2yOIRGGTT
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-02 06:29 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 06:13 by xzcv0Tehd8lQ
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-02 06:07 by xzcv0Tehd8lQ
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:07 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 21/1050 2026-11-02 05:41 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-02 03:39 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 03:26 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:47 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-02 02:34 by xzcv0Tehd8lQ
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:27 by xzcv0Tehd8lQ
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
信息提示
关闭
请填处理意见
关闭