计算化学公社

标题: VASP如何实现单轴拉伸 [打印本页]

作者
Author: 莫雨萱 时间: 2023-12-25 21:23
标题: VASP如何实现单轴拉伸
1.最近学习文献发现vasp可以对材料实现单轴拉伸模拟，最后得到应力与应变的关系，在网上看到有帖子分享了借助VASPKIT来实现，主要是通过VPKIT.in来实现应力的施加，但是我看很多文章施加应力有明确的方向，比如沿着Z轴对材料进行拉伸，但是VAPIT.in文件好像没办法明确指定应力施加方向？
2.另外，对于单轴拉伸我一直不确定自己理解的对不对：如果我需要在Z轴施加应力分析材料所发生的应变，是不是只改POSCAR中晶格常数c,而a,b保持不变哇？

作者
Author: kotori 时间: 2023-12-26 14:38
你的理解是对的，实际上如果只是做单轴应变，只需要改变POSCAR中对应方向的晶格参数就行了。
这里一般需要用到一个补丁https://github.com/Chengcheng-Xiao/VASP_OPT_AXIS重新编译VASP，因为VASP本身是不支持在结构优化的时候固定特定方向的基矢的，使用这个脚本重新编译VASP后在INCAR中使用IOPTCELL参数就可以了。

作者
Author: 莫雨萱 时间: 2023-12-27 17:13

kotori 发表于 2023-12-26 14:38
你的理解是对的，实际上如果只是做单轴应变，只需要改变POSCAR中对应方向的晶格参数就行了。
这里一般需要 ...

好的好的，确认了理解上没问题就是迈出了一大步，谢谢你的回答！
另外对于重新编译VASP来固定某一方向基矢的问题我看到网上大部分是对二维材料的优化，我的结构不是二维材料，同时在网上也看到有帖子说计模拟单轴拉伸时，不再进行晶格优化(ISIF=2)，是不是可以不用重新编译vasp添加IOPTCELL参数哇？

作者
Author: kotori 时间: 2023-12-27 17:27

莫雨萱发表于 2023-12-27 17:13
好的好的，确认了理解上没问题就是迈出了一大步，谢谢你的回答！
另外对于重新编译VASP来固定某一方向基 ...

跟是否计算二维材料关系不大，主要取决于你的需求。如果你不需要再继续优化晶格了，确实直接ISIF=2就可以了，不需要使用IOPTCELL了。

作者
Author: 莫雨萱 时间: 2023-12-27 19:22

kotori 发表于 2023-12-27 17:27
跟是否计算二维材料关系不大，主要取决于你的需求。如果你不需要再继续优化晶格了，确实直接ISIF=2就可以 ...

好的，谢谢谢谢，明白了。

作者
Author: zzaa 时间: 2024-4-2 11:25
两位大佬，问下做拉伸时候，获取的应力是直接在OUTCAR 获取in kb嘛

作者
Author: tangbo123 时间: 2024-10-18 21:37

zzaa 发表于 2024-4-2 11:25
两位大佬，问下做拉伸时候，获取的应力是直接在OUTCAR 获取in kb嘛

你好，朋友。请问你是做的什么材料啊？模型建立了多少原子？跑一个应变点大概用多久啊？

作者
Author: tangbo123 时间: 2024-10-18 21:38

莫雨萱发表于 2023-12-27 19:22
好的，谢谢谢谢，明白了。

你好，朋友。请问你是做的什么材料啊？模型建立了多少原子？跑一个应变点大概用多久啊？

欢迎光临计算化学公社 (http://bbs.keinsci.com/)