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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问怎么用VASP对二维结构进行加压呢 新人不知道如何入手 求助求助
71/1返回列表
查看: 2570  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

骆驼忍着

铜虫 (小有名气)

[求助] 请问怎么用VASP对二维结构进行加压呢 新人不知道如何入手 求助求助 已有2人参与

最近看了一篇PRB文献，对单层黑磷加单轴应变，不过用的是ABINT 我呢想用VASP模拟一下  但在网上找到的都是三维结构的计算方法  不知对二维结构要怎么做呢  求具体操作步骤
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 文献二.pdf
    • 2014-11-30 15:19:18, 2.67 M

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

技术贴

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-11-30 15:19:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★
感谢参与，应助指数 +1
franch: 金币+4, 谢谢回帖交流。。。 2014-11-30 16:40:01
简单的方法：
如果平面内晶格常数是a,b，法向方向上的晶格常数为c（带真空层）
施加应变可以通过改变晶格常数a, b的值来实现。
比如沿着a方向上的单轴应变，拉伸应变的话，将a增加，对每一个增加后的值a'，固定该值a'，优化晶格常数b的值以及原子坐标得出该单轴应变a'下的最优的b'值。对另一个a''值，重复优化晶格常数b的值以及原子坐标。
赞一下(3人) 回复此楼
2楼2014-11-30 15:41:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

骆驼忍着

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by valenhou001 at 2014-11-30 15:41:00
简单的方法：
如果平面内晶格常数是a,b，法向方向上的晶格常数为c（带真空层）
施加应变可以通过改变晶格常数a, b的值来实现。
比如沿着a方向上的单轴应变，拉伸应变的话，将a增加，对每一个增加后的值a'，固定该 ...

候老师回复我的帖子，好激动。候老师，“The molybdenum and sulfur atoms are relaxed in the strained unit cell into the minimum potential energy configuration both in and out of the plane”  这个面内和面外都弛豫是怎么做的呢？ a对应的比如3%应变下的a'是（1+3%）a吗？那这种情况下我能求出弹性模量吗？侯老师，非常感谢您  我刚开始做这个，有太多的不懂了
赞一下 回复此楼
3楼2014-12-01 10:32:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
引用回帖:
3楼: Originally posted by 骆驼忍着 at 2014-12-01 11:32:52
候老师回复我的帖子，好激动。候老师，“The molybdenum and sulfur atoms are relaxed in the strained unit cell into the minimum potential energy configuration both in and out of the plane”  这个面内和 ...

"这个面内和面外都弛豫是怎么做的呢？"
如果MoS2层放在xy平面的话，上面的一句话意思就是原子的x, y, z坐标允许被优化。

" a对应的比如3%应变下的a'是（1+3%）a吗？"
是的。

"那这种情况下我能求出弹性模量吗？"
通过计算一系列的a'，比如(1-8%)a, (1-4%)a, (1-2%) a, a  (1+2%)a , (1+4%)a, (1+8%)a等等，然后通过拟合。总能对应变量的一元二次方程拟合。。
多看看文献。可以找到很多类似的文献，它们对计算的简单理论公式和步骤都写的很详细。
如这个帖子里提到的文献：
http://muchong.com/bbs/viewthread.php?tid=8095223
赞一下(1人) 回复此楼
4楼2014-12-01 18:51:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

骆驼忍着

铜虫 (小有名气)
引用回帖:
4楼: Originally posted by valenhou001 at 2014-12-01 18:51:00
"这个面内和面外都弛豫是怎么做的呢？"
如果MoS2层放在xy平面的话，上面的一句话意思就是原子的x, y, z坐标允许被优化。

" a对应的比如3%应变下的a'是（1+3%）a吗？"
是的。

"那 ...

恩 谢谢候老师
回复此楼
5楼2014-12-01 19:01:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

骆驼忍着

铜虫 (小有名气)
引用回帖:
4楼: Originally posted by valenhou001 at 2014-12-01 18:51:00
"这个面内和面外都弛豫是怎么做的呢？"
如果MoS2层放在xy平面的话，上面的一句话意思就是原子的x, y, z坐标允许被优化。

" a对应的比如3%应变下的a'是（1+3%）a吗？"
是的。

"那 ...

候老师，您好。我下载了这篇文献研究了一下，下面是我的理解过程：单轴加压时，取F1={λ1 0；0 1}，此时根据η = 1/2(FTF − I)和
a'i= Fai，可以确定此时的应变张量矩阵，但我个人觉得这样的话咱取的3%就是λ1而不再是应变幅度δ了。我是想着能不能类比着您的教程上的块体的计算方法来计算二维材料的弹性系数，也就是取η=（δ，0,0），然后拟合总能对应变量的一元二次方程（其实你已经告诉我了，我只是在理清我的思路），但到这我发现我有两个问题又不懂了，1：我不知道对于二维材料这个总能对应变量的一元二次方程怎么推出来，其实您教程上的那个六角晶体的我就没看懂。2：类比块体的话，这个defector.x要怎么变呢？我想是不是把对应的体积改为面积？  我这么多问题，还是我的基础不好导致的，也耽误您时间了
赞一下 回复此楼
6楼2014-12-04 21:32:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wss08

新虫 (小有名气)

【答案】应助回帖

引用回帖:
2楼: Originally posted by valenhou001 at 2014-11-30 15:41:00
简单的方法：
如果平面内晶格常数是a,b，法向方向上的晶格常数为c（带真空层）
施加应变可以通过改变晶格常数a, b的值来实现。
比如沿着a方向上的单轴应变，拉伸应变的话，将a增加，对每一个增加后的值a'，固定该 ...

借楼主问题，请教一下候老师，最近我也在做二维材料应变的结构优化，出现了以下两个问题，想请教一下您：1，我的vasp是重新编译过的，未加应变时优化结构使用了OPTCELL这个文件，固定了c轴，如果是加双轴应变再优化的话，是不是改变a，b的值后直接用isif=2优化即可。2，在之前做加缺陷的二维材料结构优化时，也是采用了固定c轴，isif=3去优化，但是contcar结果中a，b值变化很大，a，b由原来的8.5，8.5分别变到了9.9和7.6，不知道是啥原因导致的，还麻烦候老师解惑，谢谢
INCAR如下：
SYSTEM = OPT
PREC = Accurate
ENCUT = 500
EDIFF = 1e-7
IBRION = 2
ISIF = 3
NSW = 500
GGA = PE
ISMEAR = 0
SIGMA = 0.05
POTIM = 0.5
EDIFFG = -1e-2
NPAR = 4
NELM = 100

OPTCELL如下：
100
110
000
赞一下 回复此楼
成功绝非偶然！
7楼2024-12-02 17:16:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 骆驼忍着 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭