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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 体相和表面做应变
91/1返回列表
查看: 2009  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿吉林大学材料学硕321求调剂 +3 Ymlll 2026-03-18 5/250 2026-03-18 19:32 by Ymlll
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 9/450 2026-03-18 15:05 by stone_128
[考研] 314求调剂 +8 无懈可击的巨人 2026-03-12 8/400 2026-03-18 14:50 by haxia
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
信息提示
关闭
请填处理意见
关闭