24小时热门版块排行榜    

查看: 3218  |  回复: 15

xiewei娜

银虫 (小有名气)

[求助] 请教MS建表面slab模型在转为POSCAR文件过程的若干问题 已有3人参与

学生想用MS建表面slab模型,并用于vasp的计算,遇到几个疑难问题,希望有经验的老师前辈帮忙解答,感激不尽!
①优化slab模型时固定中间几层的原子,而弛豫最接近表面的几层原子,是否上下表面都需要弛豫?
②在是否固定原子坐标这个设置上我用的方法很笨,对着MS的可视化模型手动找坐标然后对应到POSCAR里面设置,这样做很容易出错,不知道有没有更简便的方法?
③对于slab的厚度在ms建模的时候就可以设置的,但对于弛豫的原子层厚度,可以从哪里读取数据还是如何测得呢?如果需要同时弛豫上下表面的原子的话,是不是还要保证上下弛豫的原子层厚度一样呢?
④建表面模型加真空层后还需要建超胞吗,怎么考虑要不要建表面超胞呢?
⑤以下展示我建表面slab模型的具体操作,错误或遗漏的地方希望前辈们能加以指正:

请教MS建表面slab模型在转为POSCAR文件过程的若干问题


请教MS建表面slab模型在转为POSCAR文件过程的若干问题-1
I[MAK0LMZDQ3V%XLQ%1N0{4.png


请教MS建表面slab模型在转为POSCAR文件过程的若干问题-2
FBT8{S4CKCJREC9G7G${3YA.png


请教MS建表面slab模型在转为POSCAR文件过程的若干问题-3
LNKU~JHY8U5SAN6SZ`C}$8G.png


请教MS建表面slab模型在转为POSCAR文件过程的若干问题-4
~NPHI%L~95~G39KBWC5M8EA.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

大海12345

木虫 (正式写手)

小白

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
xiewei娜: 金币+4, 有帮助 2015-09-15 08:32:31
liliangfang: 金币+1, 谢谢交流 2015-09-18 07:55:57
关于上下表面的问题可以参考文献,大部分算的是上表面,固定底层几层原子。固定可以通过工具选中在modify里面有个constrait就是固定位置。导出vasp只要看z方向坐标查层数即可。个人认为周期性结构加真空层不需要扩胞吧。个人愚见,期待大神指点。

发自小木虫Android客户端
前途只能靠自己,记着你的责任与承诺
2楼2015-09-15 07:12:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

magicmonk

至尊木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
xiewei娜: 金币+10, ★★★很有帮助 2015-09-15 08:29:26
liliangfang: 金币+2, 谢谢交流 2015-09-18 07:56:04
①优化slab模型时固定中间几层的原子,而弛豫最接近表面的几层原子,是否上下表面都需要弛豫?
建立表面的话,最好上下对称,即上下都弛豫,这样可以避免体系产生偶极。
②在是否固定原子坐标这个设置上我用的方法很笨,对着MS的可视化模型手动找坐标然后对应到POSCAR里面设置,这样做很容易出错,不知道有没有更简便的方法?
自己编一个shell脚本,大致算法是,确定坐标开始的行数---每行坐标的第三列是z坐标----制定z坐标范围条件,并在其后面添加相应固定符号
③对于slab的厚度在ms建模的时候就可以设置的,但对于弛豫的原子层厚度,可以从哪里读取数据还是如何测得呢?
弛豫的厚度测试其实就是固定层数的测试,如果你固定的那几层在弛豫后键长层间距变化不大的话,说明是体相性质,固定得比较合适。因为表面是体相的延伸。
如果需要同时弛豫上下表面的原子的话,是不是还要保证上下弛豫的原子层厚度一样呢?
是的

④建表面模型加真空层后还需要建超胞吗,怎么考虑要不要建表面超胞呢?
一般不需要,除非你要做界面模型。具体请看文献。
⑤以下展示我建表面slab模型的具体操作,错误或遗漏的地方希望前辈们能加以指正:
1、真空层最好取15-20埃,具体需要测试
2、表面层数太少
我们最深的恐惧不是我们能力不够,我们最深的恐惧是能力超越了极限。
3楼2015-09-15 07:28:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiewei娜

银虫 (小有名气)

引用回帖:
3楼: Originally posted by magicmonk at 2015-09-15 07:28:00
①优化slab模型时固定中间几层的原子,而弛豫最接近表面的几层原子,是否上下表面都需要弛豫?
建立表面的话,最好上下对称,即上下都弛豫,这样可以避免体系产生偶极。
②在是否固定原子坐标这个设置上我用的方法 ...

感谢你的解答!
用脚本文件的话感觉还是很困难,有没有相应模板可以学习呢?
4楼2015-09-15 08:32:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiewei娜

银虫 (小有名气)

引用回帖:
2楼: Originally posted by 大海12345 at 2015-09-15 07:12:56
关于上下表面的问题可以参考文献,大部分算的是上表面,固定底层几层原子。固定可以通过工具选中在modify里面有个constrait就是固定位置。导出vasp只要看z方向坐标查层数即可。个人认为周期性结构加真空层不需要扩胞 ...

我看constrait里面固定的都是总体结构的参数,选择性固定原子的,可以再详细解说一下怎么操作吗
请教MS建表面slab模型在转为POSCAR文件过程的若干问题-5
K{P~7]RCGD]E[7_}ZRSCS]Q.png


请教MS建表面slab模型在转为POSCAR文件过程的若干问题-6
L$M7HI}Q9S~_Q9B2EE7O89R.png


请教MS建表面slab模型在转为POSCAR文件过程的若干问题-7
`FYY7[04D5SXQ$SUO]%]1T0.png

5楼2015-09-15 08:36:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiewei娜

银虫 (小有名气)

引用回帖:
3楼: Originally posted by magicmonk at 2015-09-15 07:28:00
①优化slab模型时固定中间几层的原子,而弛豫最接近表面的几层原子,是否上下表面都需要弛豫?
建立表面的话,最好上下对称,即上下都弛豫,这样可以避免体系产生偶极。
②在是否固定原子坐标这个设置上我用的方法 ...

我看文献也会算到极性表面的,那种稳定的极性表面也是可以存在的,但好像极性面由于存在偶极矩算表面能会比较麻烦,请问这个具体需要注意什么问题吗?
6楼2015-09-15 08:41:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magicmonk

至尊木虫 (著名写手)

【答案】应助回帖


liliangfang: 金币+1, 谢谢交流 2015-09-18 07:56:14
引用回帖:
6楼: Originally posted by xiewei娜 at 2015-09-15 08:41:32
我看文献也会算到极性表面的,那种稳定的极性表面也是可以存在的,但好像极性面由于存在偶极矩算表面能会比较麻烦,请问这个具体需要注意什么问题吗?...

需要进行偶极校正,vasp手册里有讲
我们最深的恐惧不是我们能力不够,我们最深的恐惧是能力超越了极限。
7楼2015-09-15 09:05:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magicmonk

至尊木虫 (著名写手)

【答案】应助回帖


liliangfang: 金币+1, 谢谢交流 2015-09-18 07:56:20
引用回帖:
5楼: Originally posted by xiewei娜 at 2015-09-15 08:36:17
我看constrait里面固定的都是总体结构的参数,选择性固定原子的,可以再详细解说一下怎么操作吗

K{P~7]RCGD]EQ.png

L$M7HI}Q9S~_Q9B2EE7O89R.png

`FYY7%]1T0.png
...

我主要用vasp,对MS不太熟,抱歉!
我们最深的恐惧不是我们能力不够,我们最深的恐惧是能力超越了极限。
8楼2015-09-15 09:06:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大海12345

木虫 (正式写手)

小白

引用回帖:
5楼: Originally posted by xiewei娜 at 2015-09-15 08:36:17
我看constrait里面固定的都是总体结构的参数,选择性固定原子的,可以再详细解说一下怎么操作吗

K{P~7]RCGD]EQ.png

L$M7HI}Q9S~_Q9B2EE7O89R.png

`FYY7%]1T0.png
...

固定分数位置即可,要是想让某个方向固定就选上面的
前途只能靠自己,记着你的责任与承诺
9楼2015-09-15 09:27:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiewei娜

银虫 (小有名气)

引用回帖:
9楼: Originally posted by 大海12345 at 2015-09-15 09:27:28
固定分数位置即可,要是想让某个方向固定就选上面的...

按你这样做的话会不会所有原子的设置都一样,我想达到的效果是上下表面原子不固定,中间的固定。
10楼2015-09-15 15:29:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiewei娜 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 范进中举一文的中心思想 +5 炎黄贵胄 2026-08-22 6/300 2026-08-24 06:34 by sincosx
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见