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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 体相和表面做应变
91/1返回列表
查看: 1936  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

菲顶顶

铜虫 (正式写手)

[求助] 体相和表面做应变 已有1人参与

想请问大家，用vasp计算分别计算体相和表面情况的单轴和双轴应变的时候的方法是一样的吗？
我的想法是直接改变晶格常数不就行了，无论体相还是表面，如果是单轴的就改变一个方向的晶格常数，双轴的同时改变两个方向的晶格常数，然后优化。ISIF是怎么取的？体相取3？表面取2？
可是我在查一些资料的时候看到大家说的什么拉伸哪个方向就固定该方向不优化之类之类的。表示不懂，有做过应力应变的大神们可否来解答一下。
回复此楼

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

应变

» 猜你喜欢

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

stractor

金虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
菲顶顶: 金币+10, ★★★很有帮助 2015-01-31 22:24:33
liliangfang: 金币+2, 谢谢交流 2015-02-01 09:49:00
“ISIF是怎么取的？体相取3？表面取2？”
不能直接取3和2，那样弛豫后又回到平衡状态，应变消失。
听说有人改写了VASP弛豫的代码，可以实现你的要求。
还有一种方式就是EOS拟合，比如，我在xy平面施加双轴应变，让z方向逐渐改变一点晶格常数(一些列dz)，计算在该应变下，不同z的改变对于的体系总能，然后通过E~dz找到能量最低的.
赞一下 回复此楼
2楼2015-01-31 21:57:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
引用回帖:
2楼: Originally posted by stractor at 2015-01-31 21:57:39
“ISIF是怎么取的？体相取3？表面取2？”
不能直接取3和2，那样弛豫后又回到平衡状态，应变消失。
听说有人改写了VASP弛豫的代码，可以实现你的要求。
还有一种方式就是EOS拟合，比如，我在xy平面施加双轴应变， ...

这样确实可以的，我先试试，那ISIF到底该取什么呢，你的意思就是不优化直接计算自洽能量是吧。
赞一下 回复此楼
3楼2015-01-31 22:24:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
菲顶顶: 金币+10, ★★★很有帮助 2015-02-01 14:08:19
引用回帖:
3楼: Originally posted by 菲顶顶 at 2015-01-31 22:24:26
这样确实可以的，我先试试，那ISIF到底该取什么呢，你的意思就是不优化直接计算自洽能量是吧。...

不用优化了，直接电子自洽计算。然后通过总能来确定。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
4楼2015-02-01 11:21:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
引用回帖:
4楼: Originally posted by stractor at 2015-02-01 11:21:33
不用优化了，直接电子自洽计算。然后通过总能来确定。
...

好的谢谢。
我还有个问题想和您讨论下。这种方法对双轴应变适用，那如果是单轴应变呢，需要只固定一个方向，优化另外两个方向，用这种方法，是同时改变另外两个方向的大小，这样就很难确定同时改变多少。

我自己想到一个方法就是可以通过固定POSCAR的xy坐标和优化z坐标，应该也是可以的吧。这样就可以同时优化多个方向，只是不知道可不可行。
赞一下 回复此楼
5楼2015-02-01 14:08:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
"我自己想到一个方法就是可以通过固定POSCAR的xy坐标和优化z坐标，应该也是可以的吧。这样就可以同时优化多个方向，只是不知道可不可行"
你是说xy平面加应变吧，感觉行不通，因为你用ISIF=3的话，x，y的长度变了，如果用ISIF=2，你的施加应变的胞的Z方向的晶格长度又如何取呢

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
6楼2015-02-02 01:49:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
引用回帖:
6楼: Originally posted by stractor at 2015-02-02 01:49:40
"我自己想到一个方法就是可以通过固定POSCAR的xy坐标和优化z坐标，应该也是可以的吧。这样就可以同时优化多个方向，只是不知道可不可行"
你是说xy平面加应变吧，感觉行不通，因为你用ISIF=3的话，x，y的 ...

好的 谢谢 我还听人说可以在想加应变的方向改变晶格常数，然后用ISIF=4，自己测试过了，感觉结果是不对的。非常感谢你
赞一下 回复此楼
7楼2015-02-02 15:14:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redsnowolf

银虫 (小有名气)
vasp算应变的时候，变形是通过修改POSCAR自己加的，并且应变量要逐步增大，算大变形的时候不能把应变量一步加到位。结构优化前要修改constr_cell_relax.F文件，将施加应变的方向固定不优化，然后ISIF该怎么取就怎么取。例如只在x方向上有变形，可以在constr_cell_relax.F里加上FCELL(1,1)=0，把x方向上的力置零，重新编译，再优化时x方向的晶格常数就不变了。其它方向的、单双轴变形的情况类似，具体的在constr_cell_relax.F文件里有说明。
赞一下(1人) 回复此楼
8楼2015-02-04 14:15:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
引用回帖:
8楼: Originally posted by redsnowolf at 2015-02-04 14:15:28
vasp算应变的时候，变形是通过修改POSCAR自己加的，并且应变量要逐步增大，算大变形的时候不能把应变量一步加到位。结构优化前要修改constr_cell_relax.F文件，将施加应变的方向固定不优化，然后ISIF该怎么取就怎么 ...

谢谢
回复此楼
9楼2015-02-04 21:38:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 菲顶顶 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭