24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2795  |  回复: 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娜 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +7 张.1 2026-04-05 7/350 2026-04-05 20:40 by 啵啵啵0119
[考研] 338求调剂 +3 我想上岸ii 2026-04-05 3/150 2026-04-05 19:59 by nepu_uu
[考研] 080200学硕,机械工程专业277分,求带走! +7 瓶子PZ 2026-03-31 7/350 2026-04-05 17:49 by liucky
[考研] 288求调剂 一志愿哈工大 材料与化工 +13 洛神哥哥 2026-04-03 13/650 2026-04-05 17:27 by zzx2138
[考研] 283求调剂 +10 A child 2026-04-04 10/500 2026-04-05 08:22 by qlm5820
[论文投稿] 求文献 5+3 ys879651$ 2026-04-02 3/150 2026-04-04 17:22 by bobvan
[考研] 0710生物学336分求调剂 +6 kiyy 2026-04-01 8/400 2026-04-04 10:10 by kiyy
[考研] 本科985,专业0812分336求调剂 +4 莫莫很行 2026-04-03 4/200 2026-04-03 21:31 by zhq0425
[考研] 286求调剂 +8 lim0922 2026-04-02 8/400 2026-04-03 20:19 by rzh123456
[考研] 建环,能源,土木老师路过看一看!!! +5 嘿嘿uu 2026-04-01 5/250 2026-04-03 11:47 by znian
[基金申请] 请问共同通讯和共同一作的认可度问题 10+4 psa1234 2026-04-01 10/500 2026-04-03 11:08 by Kittylucky
[考研] 材料考研调剂 +10 Gs大王 2026-04-02 10/500 2026-04-03 09:47 by 遗忘消失的灆
[考研] 调剂 +7 祉岷. 2026-04-02 7/350 2026-04-03 09:11 by 花呗还欠600
[考研] 318求调剂,计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 271求调剂 +15 勒布朗@ 2026-03-31 20/1000 2026-04-02 11:24 by Sammy2
[考研] 335求调剂 +3 321* 2026-03-31 4/200 2026-04-01 00:00 by 321*
[考研] 274求调剂 +6 xiao爱同学 2026-03-30 6/300 2026-03-31 10:04 by cal0306
[考研] 323分 食品与营养调剂 +3 嘿ooo 2026-03-31 3/150 2026-03-31 09:38 by longlotian
[考研] 370求调剂 +3 080700调剂 2026-03-30 3/150 2026-03-31 01:09 by A_Zhe
[考研] 材料专硕 085600求调剂 +7 BBQ233 2026-03-30 7/350 2026-03-30 17:44 by oooqiao
信息提示
请填处理意见