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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助完毕】vasp-slab-边界附近原子顶位吸附小分子的结构优化很难收敛吗？
61/1返回列表
查看: 1557  |  回复: 5
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

cenwanglai

荣誉版主 (知名作家)

[交流] 【求助完毕】vasp-slab-边界附近原子顶位吸附小分子的结构优化很难收敛吗？

我的表面是p(3*3)大小的，在表面中间吸附、掺杂都能100步内收敛。但是在靠近边界的最表层原子顶位上的吸附（如吸附SO2）结构,力就很难收敛，都300多步了，而且能量在小数点后第三位上上下波动。

同样在边界附近去掉一个最表层原子（O），做一个空位结构，也能较快收敛。但同一系列的自洽计算（ISMEAR=0,SIGMA=0.05）,scf的EENTRO~=.03。一般中间缺陷或者吸附结构在scf计算中第六电子步就是0.0000了。

问题：
（1）边界上的缺陷和吸附结构优化有什么要注意的呢？
（2）要增大表面大小，把边界外推，让吸附结构靠近表面中心？
（3）超胞结构优化（cell fixed)边界与表面中心的处理不一样吗？

[ Last edited by cenwanglai on 2011-4-25 at 20:13 ]
回复此楼

» 收录本帖的淘帖专辑推荐

VASP 杂七杂八

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

查看全部散金贴
1楼2010-12-13 16:17:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)
★ ★
zzy870720z(金币+2):鼓励交流 2010-12-13 22:59:04
如果有图会比较容易理解，你的表面是周期边界条件的吗？边界是指在cell的边界处吗？
按我的理解，在边界处的计算应该和中心计算没有太多的区别，唯一可能的原因就是原子一离开cell就会做一个加法或减法回到cell，这种加法减法可能会导致误差。
顶位的吸附也是一个问题，有可能会有多个暂稳态使得原子在多个暂稳态之间跳跃
赞一下(1人) 回复此楼
2楼2010-12-13 21:09:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)
引用回帖:
Originally posted by bingmou at 2010-12-13 21:09:26:
如果有图会比较容易理解，你的表面是周期边界条件的吗？边界是指在cell的边界处吗？
按我的理解，在边界处的计算应该和中心计算没有太多的区别，唯一可能的原因就是原子一离开cell就会做一个加法或减法回到cell， ...

是周期性边界条件。你完全明白了我的问题所在。

我的POTIM调到0.1了。是不是可以再调小一点？
赞一下 回复此楼
3楼2010-12-14 09:10:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)
★
sunyang1988(金币+1):谢谢交流 2010-12-14 09:17:37
cenwanglai(金币+2):是氧化物表面。中间是有杂原子的。除此之外，没有差别。杂原子对表面有一定扭曲。 2010-12-14 11:28:57
氧化物表面么？

能确定边界和中心的环境是一样的么？
赞一下(1人) 回复此楼
4楼2010-12-14 09:16:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)
★ ★
zzy870720z(金币+2):谢谢指教 2010-12-14 09:33:26
cenwanglai(金币+2):表面被吸附的小分子受力无法收敛。其它结构能很好的收敛到0.02.边界附近吸附势只能到0.05. 2010-12-14 11:30:14
引用回帖:
Originally posted by cenwanglai at 2010-12-14 09:10:04:


是周期性边界条件。你完全明白了我的问题所在。

我的POTIM调到0.1了。是不是可以再调小一点？

我觉得你可以仔细研究是一下是哪些原子的力无法收敛，以及他们为什么在那里无法收敛。POTIM已经很小了，我觉得没有必要修改，而且VASP在弛豫的过程中是修改POTIM的，如果发现不对会缩小POTIM，如果发现正确会加大POTIM
赞一下(1人) 回复此楼
5楼2010-12-14 09:24:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

葵cc_

木虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
你好，请问你是在根据什么来判断是否收敛的？查看什么量，怎样才算收敛？
赞一下 回复此楼
6楼2014-04-03 19:06:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cenwanglai 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 我感觉当老师好累啊。。 +19 fairy1122 2024-11-15 20/1000 2024-11-16 14:16 by lx7991
[基金申请] 求助！请问海外博后依托内地单位申请青基的优劣？ 80+3 2280999712 2024-11-15 10/500 2024-11-16 13:41 by dxcharlary
[精细化工] 同一个反应相同的反应条件，是不是反应结果相差不大？ +6 青霉素 2024-11-11 6/300 2024-11-16 13:35 by linjunc
[第一性原理] 有没有要一起买VASP版权的，可以6个人 5+3 变成小神 2024-11-14 9/450 2024-11-16 13:33 by dxcharlary
[论文投稿] 论文返修状态变成了awaiting AE recommendation 10+4 猪小耍 2024-11-13 15/750 2024-11-16 09:24 by guxue
[教师之家] 这种人痛苦吗 +12 2671 2024-11-12 14/700 2024-11-15 22:15 by 仇恨毁灭财富
[论文投稿] chemical science和advanced science哪个好一点 +6 yly150 2024-11-12 6/300 2024-11-15 21:06 by hcj9909053
[论文投稿] 找不到审稿人 +6 吉祥如意141913 2024-11-13 6/300 2024-11-15 16:48 by LEyblyt
[论文投稿] 核心期刊那个阶段最难 +4 公共茅斯 2024-11-11 4/200 2024-11-15 16:44 by sugarmei
[基金申请] 广东省自然科学基金-面上项目消息有了么？ +4 pachang 2024-11-14 6/300 2024-11-15 09:39 by zhoulei7170
[基金申请] 求助 +4 Enenenene 2024-11-15 4/200 2024-11-15 09:36 by 榨菜拌青椒
[论文投稿] 投稿系统中的通讯作者和文章中的通讯作者不一样，文章目前被录用了？ +4 babybabygo 2024-11-12 5/250 2024-11-14 19:13 by 走了002
[硕博家园] 研究生的生活该是什么样 +4 lqy0719 2024-11-14 4/200 2024-11-14 16:45 by 阿荣喝酒
[有机交流] 同一个反应回流情况不同 20+3 1853846 2024-11-12 3/150 2024-11-14 15:42 by 太阳谷
[硕博家园] 大龄已婚想读博如何 +15 米娅阳 2024-11-11 18/900 2024-11-14 14:18 by xiaomi0401
[基金申请] 76批博后基金 +3 feiyi3986 2024-11-14 3/150 2024-11-14 11:50 by puly
[论文投稿] 要不要撤稿另投 10+4 wangzhesd 2024-11-09 8/400 2024-11-14 09:09 by LIU_V
[硕博家园] 如何看清华北大硕博的毕业论文 +6 2025lucky 2024-11-10 6/300 2024-11-12 22:46 by 楚扬长流
[论文投稿] 论文接受后，编辑部发邮件提醒缺少附件 15+4 上善若水明泪 2024-11-10 4/200 2024-11-12 14:24 by 北京莱茵润色
[论文投稿] 爱思唯尔投稿系统里的通讯作者可以和文章里的通讯作者标注不同吗 +7 Omnissiah 2024-11-10 7/350 2024-11-12 14:07 by holypower
信息提示
关闭
请填处理意见
关闭