24小时热门版块排行榜    

查看: 13973  |  回复: 77
本帖产生 1 个 1ST强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[资源] calculating ideal strength of crystal using VASP

http://code.google.com/p/ideal-strength-vasp/downloads/list
Recently, I write a python script to calculate ideal strength of crystal combination with VASP. You can download it from above website.
the details is found from readme in wiki at that website.

Introduction
About method details, please see PRL 82,2713(1999)

Details
First you need to modify vasp code. If you want to calculate ideal tensile strength, you should add 'FCELL(1,1)=0.0' to constr_cell_relax.F of vasp code. And if you choose ideal shear strength, you should add 'FCELL(1,3)= 0.0' and 'FCELL(3,1) = 0.0' to constr_cell_relax.F of vasp code. Here, you need to recompile vasp.

second, you should prepare input.dat and this file as below:

POSCAR
0.02 #strain
20 #step
-45.0 -35.264390 0.0 # rotate Z, Y and X.
1 # 1 tensile, 2 shear
mpiexec -np 8 vasp.4.6

The first line is the name of POSCAR. The second line is strain of distortion. The third line is total step of distortion. The fourth line is degree of rotation. This is for calculating special orientation. For example, if you want to calculating ideal strength of Diamond along 100, you just set 0.0 0.0 0.0. If you want to calculate 110 orientation, you need set -45.0 0.0 0.0. If you want to calculate 111 orientation, you need set -45.0 -35.264390 0.0. The five line is to choose tensile of shear. The sixth line is execute command.



[ Last edited by gavinliu7390 on 2011-5-20 at 14:09 ]

由于最近发现很多人,管我要执行的例子,特此传一个,希望有助于大家。
如有问题,可以发email给我 lhy@calypso.cn[ Last edited by gavinliu7390 on 2013-11-4 at 01:42 ]

由于之前用的google code,但是google code虽然现在支持下载,但以后有可能关闭。所以我又提供了一个新的下载链接。
https://sourceforge.net/projects/ideal-strength-vasp/

[ Last edited by gavinliu7390 on 2015-3-14 at 08:31 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : diamond_strength.tar.bz2
  • 2013-11-04 01:40:57, 529.88 K

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenguo_whu

新虫 (初入文坛)


老师,您好。想问下计算一些参杂了异类原子体系的Tensile stress-strain,得到的数据作图,遇到曲线上下波动频繁的情况怎么办?
72楼2018-04-10 20:45:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 78 个回答
好东西,值得收藏!
5楼2011-05-20 14:19:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sjtuwyf

铜虫 (小有名气)


★★★★★ 五星级,优秀推荐

Very good! Thank you !
7楼2011-06-26 23:29:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao7910

金虫 (小有名气)


★★★ 三星级,支持鼓励

1。我测试剪切应力过程中在终端出现:hit a member that was already found in another star。不知哪里出了问题。
2。楼主如果用拉伸和剪切完的晶格直接进行计算会导致垂直于应变方向上也会有很大的残余应力。所以希望楼主能把程序改进下,对拉伸和剪切后能够对垂直于应变方向进行晶格和原子位置的驰豫,使得在垂直应变方向应力尽可能小。
我把运行中出的提示和我测试的文件打包给你,希望能告知原因何在?
8楼2011-06-27 22:51:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
sars5183楼
2011-05-20 13:18   回复  
五星好评  感谢分享~
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见