24小时热门版块排行榜    

查看: 13972  |  回复: 77
【奖励】 本帖被评价58次,作者gavinliu7390增加金币 45.8
本帖产生 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

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

精华网帖收集 仿真建模与计算 研究生资料集 材料计算模拟实用技巧
VASP 第一性原理计算-tg 第一性原理计算经验 第一原理资料贴
计算-vasp 程序语言箱 第一性原理计算辅助工具 第一性原理和电化学
MS资料 分子动力学 Photochemistry 有效质量
计算化学-辅助工具 vasp 各种各样的软件

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
好东西,值得收藏!
5楼2011-05-20 14:19:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 78 个回答

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的回帖
引用回帖:
Originally posted by maomao7910 at 2011-06-27 22:51:07:
1。我测试剪切应力过程中在终端出现:hit a member that was already found in another star。不知哪里出了问题。
2。楼主如果用拉伸和剪切完的晶格直接进行计算会导致垂直于应变方向上也会有很大的残余应力。所 ...

这不是我程序的问题, 是你没有改好vasp代码。
因为vasp必须改源码,才能使它选择性优化。
除非你用其它软件 (e.g. abinit).
可以看看上面的图, 计算是没有任何问题的。

[ Last edited by gavinliu7390 on 2011-6-28 at 21:53 ]
9楼2011-06-28 21:51:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
sars5183楼
2011-05-20 13:18   回复  
五星好评  感谢分享~
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见