24小时热门版块排行榜    

查看: 2764  |  回复: 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娜 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +4 素年祭语 2026-03-15 4/200 2026-03-16 09:42 by 闲人终南山
[考研] 326求调剂 +4 上岸的小葡 2026-03-15 5/250 2026-03-16 08:39 by Linda Hu
[考研] 调剂 +8 调剂的考研学生 2026-03-09 8/400 2026-03-15 22:14 by Winj1e
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 304求调剂 +5 小熊joy 2026-03-14 5/250 2026-03-14 21:07 by peike
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 求调剂 +5 鹤遨予卿 2026-03-09 5/250 2026-03-14 02:44 by JourneyLucky
[考研] 化学工程321分求调剂(南京工业,浙江工业) +3 大米饭! 2026-03-09 4/200 2026-03-14 02:34 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考研] 296求调剂 +3 大口吃饭 身体健 2026-03-13 3/150 2026-03-13 10:31 by 学员8dgXkO
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
请填处理意见