极客学术

 找回密码
 立即注册
查看: 50|回复: 10

[交流] 如何用VASP优化二维材料结构

[复制链接]

25

主题

40

帖子

125

积分

注册会员

Rank: 2

积分
125
发表于 2018-11-29 15:49:59 | 显示全部楼层 |阅读模式
一般来说,二维材料会有一个非常大的真空层作为z轴,而VASP的ISIF参数常用的有3个,ISIF=2表示不动晶格只优化原子位置,ISIF=3表示晶格,体积与原子位置都优化,而ISIF=4表示晶格与原子位置优化,体积不发生变化。
二维材料比如说从实验的范德华式的层状材料剥离出来,因此需要优化原子位置与晶格参数,请问各位虫友都是怎么样解决这个问题,我所了解的至少有以下方案:
1. 使用ISIF=3优化,在PBE+D3泛函下,我发现真空层会减少很多,甚至直接减少到零,然后报错。精度不是很高的情况下(10e-5 eV与 0.005 eV/A)能够优化出结果。
2. 使用ISIF=4优化,基本上由于体积限制真空层不会变化太多,然后对于优化后的平衡晶格常数进行0.5% 1%的拉伸,画出E-V曲线,找到能量最低点。
3. 使用ISIF=3优化,但是重新编译VASP,使得其只改变xy平面内的参数,z轴真空层能够保持不变,z轴也始终与xy平面垂直。或许先粗糙再精确的收敛条件能够得到比较好的结果。
不知各位是怎么看的?二维材料也是比较热的材料体系,还是想拿到一个尽可能精确的结构!


回复

使用道具 举报

26

主题

50

帖子

152

积分

注册会员

Rank: 2

积分
152
发表于 2018-11-29 15:50:48 | 显示全部楼层
修改vasp源文件
实现拉伸x方向,优化y的晶格矢量,x和z方向晶格矢量不优化
原子位置优化
回复

使用道具 举报

22

主题

36

帖子

110

积分

注册会员

Rank: 2

积分
110
发表于 2018-11-29 15:51:04 | 显示全部楼层
3最方便
然后还可以手动优化,稍嫌麻烦
回复

使用道具 举报

26

主题

50

帖子

152

积分

注册会员

Rank: 2

积分
152
发表于 2018-11-29 15:51:21 | 显示全部楼层
3不是只适用于正交晶系吗?还有手动优化是不是手动改晶格常数,看能量最低。
回复

使用道具 举报

22

主题

36

帖子

110

积分

注册会员

Rank: 2

积分
110
发表于 2018-11-29 15:51:26 | 显示全部楼层
我觉得3不只是正交吧?
手动优化你说的对
回复

使用道具 举报

26

主题

50

帖子

152

积分

注册会员

Rank: 2

积分
152
发表于 2018-11-29 15:51:45 | 显示全部楼层
麻烦再请教你2个问题,
1.我看你在这个贴子(http://muchong.com/t-9268257-1)里说你分别用脚本和改原文件优化了蓝磷烯,蓝磷烯是六方格子,那改源文件应该对于所有晶系都可行的。
2..对于二维材料的优化,如果用脚本优化应该和优化块状晶体然后拟合状态方程的脚本不太一样吧,二维材料应该不能用拟合状态方程画E-V曲线来找最优的晶格参数吧。那这样脚本的意思就只是按比例改变a,b轴(一般a=b),看能量最低,其实和手动改晶格常数一个道理吧。但为什么楼主说的方法2,用ISLF=4优化过以后得到的结构进一步加小的伸缩画的是E-V曲线呢。
回复

使用道具 举报

22

主题

36

帖子

110

积分

注册会员

Rank: 2

积分
110
发表于 2018-11-29 15:51:59 | 显示全部楼层
跟3D没有本质区别,主要根据对称性看有几个独立晶格常数。比如块体MoS2,只有两个,平面内是相同的啊
回复

使用道具 举报

26

主题

50

帖子

152

积分

注册会员

Rank: 2

积分
152
发表于 2018-11-29 15:52:05 | 显示全部楼层
好的,意思就是说我用侯老师教程里地那个脚本改变晶格缩放系数优化单层MoS2也是适用的,但优化二维的材料c方向不应该固定的吗?
回复

使用道具 举报

22

主题

36

帖子

110

积分

注册会员

Rank: 2

积分
110
发表于 2018-11-29 15:52:17 | 显示全部楼层
那就更方便了啊,单层mos2的c轴固定,只拉伸一个平面内方向完事。
回复

使用道具 举报

26

主题

50

帖子

152

积分

注册会员

Rank: 2

积分
152
发表于 2018-11-29 15:52:22 | 显示全部楼层
好的,多谢。突然想到那楼主说的方法2中第一部先用 ISLF=4优化得到一个较准确的结构后一步对得到的结构改用ISLF=2,然后手动改晶格常数进行微凉的拉伸过挤压是否也可行。
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|Archiver|手机版|小黑屋|极客学术

GMT+8, 2018-12-13 11:34 , Processed in 0.242415 second(s), 6 queries , File On.

Powered by Discuz! X3.3

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表