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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问怎么用VASP对二维结构进行加压呢 新人不知道如何入手 求助求助
71/1返回列表
查看: 2899  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +8 Tide man 2026-08-11 8/400 2026-08-12 17:30 by Vivilian
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] FileCode能看出啥？ +8 要乐观耀哥 2026-08-10 21/1050 2026-08-12 16:08 by 云上清扬
[基金申请] filecode +13 documentary 2026-08-10 15/750 2026-08-12 15:56 by 云上清扬
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 35/1750 2026-08-12 14:38 by sdfapple719
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 是这周出结果还是下周出结果？ +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
信息提示
关闭
请填处理意见
关闭