| 查看: 2791 | 回复: 7 | ||||||
| 【奖励】 本帖被评价4次，作者浔荆增加金币 3.2 个 | ||||||
| 本帖产生 1 个 1ST强帖 ，点击这里进行查看 | ||||||
[资源]
用VASP固定应力计算的python脚本（20171016更新 增加了对切应力的支持）
|
||||||
|
----------------------------------------- 原帖： 脚本目的：用VASP实现指定应力张量下材料的DFT模拟，目前只支持主应变，尚未考虑切向应力 由于VASP程序本身不支持固定应力计算，本脚本通过广义胡克定律反推获得指定应力所需的应变，相应的调整POSCAR，不断重复迭代该过程直至获得所需应力状态。 附件是对bcc钨的单轴拉伸算例，脚本头部注释行有使用说明。 若发现任何bug，请致信jie.hou2@mail.mcgill.ca ----------------------------------------- 20171016更新： 最近有算切应力的需求，还是完善了对切应力的支持。附件1.0中是bcc钨的一个算例，随便设了一个应力张量： -100 10 -50 10 100 100 -50 100 50 测试正常收敛， 另附简单说明 Instructions: 1. Input the pressure (stress) tensor in the script 2. Input estimated elastic modulus in the script 3. Set ISIF=2 in the INCAR 4. Copy POSCAR to poscar.0 5. Configure mpiexec commands in the script based on your linux environment 6. Use 'qsub $YourJobPath/fixpressure.py' to submit your job Notes: This script calculates stress tensor by using vasp program, and adjusts the POSCAR iteratively until the stress tensor converges to the value you want. The POSCAR adjustment is estimated according to generalized Hooke's law based on elastic modulus you inputted. The iteration algorithm can tolerate certain errors in the modulus, but bad modulus input may lead to a slow convergence. In order to obtain an accurate stress result, please use a high ECUT. See Vasp manual for details. You can set ISTART=1 and LWAVE=.T. in the INCAR file, this can accelerate the calculation. However, this may also raise convergence problems when shear stresses are applied. Do your own tests. Shoule you find any bug, please contact jie.hou2@mail.mcgill.ca.[ Last edited by 浔荆 on 2017-10-17 at 08:51 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : fixpressure.rar
- 附件 2 : fix-stress1.0.rar
2017-09-30 05:25:27, 89.32 K
2017-10-17 02:19:47, 127.43 K
» 收录本帖的淘帖专辑推荐
 vasp+wannier90 | dft_vasp | 方法不错 | vasp |
» 猜你喜欢
让我中一个面上吧！
已经有5人回复
-
科研孤儿太难了
已经有12人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
-
重要消息，中午系统在维护
已经有12人回复
-
今天基金会出结果吗？20260819
已经有15人回复
» 本主题相关商家推荐: (我也要在这里推广)
3楼2017-09-30 10:10:20
4楼2017-09-30 10:12:01
5楼2017-10-02 11:59:49
简单回复
漫天飘雪2楼
2017-09-30 08:43
回复
五星好评 
雁儿霏霏6楼
2017-10-03 10:54
回复
五星好评 顶一下，感谢分享！
2017-10-04 17:08
回复
五星好评 顶一下，感谢分享！
浔荆8楼
2017-10-17 03:08
回复
40