| 查看: 15085 | 回复: 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朵）
» 猜你喜欢
可靠消息，8月21号(明天)放榜
已经有3人回复
-
只有每年这种时候来逛逛小木虫
已经有20人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有125人回复
-
今天基金会出结果吗？20260819
已经有15人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何用VASP做第一性原理拉伸？
已经有7人回复
- 测试POSCAR的脚本不运行
已经有10人回复
- preparation of znse nanocrystals using water-in-oil microemulsions
已经有6人回复
- 【分享】CRC出版社Handbook of Advanced Ceramics Machining（已搜索无重复）
已经有152人回复
- 700 pieces of composition written by american
已经有25人回复
- [求助]硫酸浸渍三氧化钨的原理
已经有7人回复
- 怎么准备proof/certificate of PhD degree?
已经有3人回复
- 关于oscillator strength
已经有3人回复
20楼2013-10-28 16:32:02
★★★★★ 五星级,优秀推荐
|
感谢楼主的分享！这里有些疑惑，一直没有弄明白，希望楼主帮忙答疑解惑，谢谢了！！ 1、这个 python script 主要功能是什么？它的作用主要是产生一个加若干大小应变以后的POSCAR吗？我们在使用的时候，是否需要写个脚本，先对产生的POSCAR做relax，然后在做static scf，即需要准备两个INCAR（INCAR.relax 和INCAR.static）？ 2、“If you want to calculate ideal tensile strength, you should add \'FCELL(1,1)=0.0\' to constr_cell_relax.F of vasp code.” 我记得您以前有一个帖子（http://muchong.com/bbs/viewthread.php?tid=1928920&page=1）是说，先在拉伸方向建立原胞，然后设置拉伸方向为X轴；再在constr_cell_relax.F中加入： SAVY=FCELL(2,2) SAVZ=FCELL(3,3) FCELL=0 FCELL(2,2)=SAVY FCELL(3,3)=SAVZ 只弛豫y和z方向，这样即可计算tensile strength了。那现在使用python script代码是不是只要在constr_cell_relax.F中增加一个\'FCELL(1,1)=0.0\'即可了吗（如果要算ideal tensile strength）？ 3、在input.dat中，您提到 \"The fourth line is degree of rotation ”，这一行的作用是不是实现将加的应变作用旋转至特定的方向，如（110），那么的话，这个数值是怎么计算得到的？ |
13楼2012-02-16 20:35:02
|
本帖内容被屏蔽 |
77楼2020-08-10 10:09:49
15楼2012-03-07 20:18:14
72楼2018-04-10 20:45:10
2楼2011-05-20 12:10:29
好东西，值得收藏！ 5楼2011-05-20 14:19:51
7楼2011-06-26 23:29:59
8楼2011-06-27 22:51:07
9楼2011-06-28 21:51:35
10楼2011-06-29 13:14:54
11楼2011-07-12 17:01:54
14楼2012-03-07 19:05:01
16楼2013-09-02 22:39:04
17楼2013-09-07 16:48:59
18楼2013-10-15 20:12:27
19楼2013-10-16 02:37:33
27楼2013-11-21 12:09:41
35楼2015-01-06 09:27:47
36楼2015-01-14 19:09:01
38楼2015-03-06 11:12:38
40楼2015-04-20 09:59:18
41楼2015-04-20 10:50:16
45楼2015-10-23 13:37:54
简单回复
sars5183楼
2011-05-20 13:18
回复
五星好评 感谢分享~
2011-05-20 13:51
回复
五星好评 
kgdu6楼
2011-05-20 16:33
回复
五星好评 谢谢分享！！
cup.qzu12楼
2011-12-23 11:25
回复
五星好评
墨明棋妙陈21楼
2013-11-04 08:54
回复
顶一下，感谢分享！
墨明棋妙陈22楼
2013-11-04 09:01
回复
送红花一朵
我等我的爱23楼
2013-11-05 20:35
回复
五星好评 顶一下，感谢分享！
crante24楼
2013-11-07 09:55
回复
五星好评 顶一下，感谢分享！
了甲灰25楼
2013-11-20 09:08
回复
五星好评
dyc_200826楼
2013-11-20 16:48
回复
五星好评 顶一下，感谢分享！
漫天飘雪28楼
2013-12-07 22:27
回复
五星好评 顶一下，感谢分享！
haizhiwuyu29楼
2014-02-10 16:44
回复
五星好评 顶一下，感谢分享！
for_xili30楼
2014-03-02 09:48
回复
五星好评 
beibeilover31楼
2014-04-20 18:42
回复
五星好评 顶一下，感谢分享！
123左叉32楼
2014-10-17 17:21
回复
五星好评 顶一下，感谢分享！
曾犟33楼
2014-12-21 14:48
回复
三星好评 顶一下，感谢分享！
123左叉34楼
2015-01-06 08:46
回复
顶一下，感谢分享！
txcokokok37楼
2015-02-10 19:34
回复
五星好评 顶一下，感谢分享！
maomao791039楼
2015-03-24 09:47
回复
顶一下，感谢分享！
Honorjian42楼
2015-06-29 15:43
回复
五星好评 顶一下，感谢分享！
zhy小桥流水43楼
2015-09-12 17:58
回复
五星好评 顶一下，感谢分享！
Hysen44楼
2015-09-29 22:29
回复
五星好评 顶一下，感谢分享！
冰河世纪！46楼
2015-12-07 14:49
回复
五星好评 顶一下，感谢分享！
cheapskate47楼
2015-12-15 11:52
回复
五星好评 顶一下，感谢分享！
xieqingxing48楼
2016-03-14 21:57
回复
五星好评 顶一下，感谢分享！
lisu55549楼
2016-05-05 10:19
回复
五星好评 顶一下，感谢分享！
shxincui51850楼
2016-05-19 10:26
回复
五星好评 顶一下，感谢分享！
