| 查看: 2246 | 回复: 8 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
Sparks__Fly新虫 (小有名气)
|
[求助]
计算FeCr中掺杂He的形成能（初学者还请见谅） 已有2人参与
|
|
|
先用MS建模2X2X2的Fe原子模型，再把某个原子替换成Cr原子，之后在不同的间隙位置上掺杂He原子，计算形成能，从而判断He原子可能出现的位置。 第一步先是结构优化：INCAR： SYSTEM=onestep-FeCrHe-1 ENCUT=250（是否合理？） ISTART=0;ICHARG=2 ISMEAR=1;SIGMA=0.2 NSW=3（是否过小？）;IBRION=2 ISIF=3 POTIM=0.2 EDIFF=1E-6;EDIFFG=-1E-3 PREC=M POSCAR选择为Selective dynamics，仅仅固定He原子的位置（是否正确？） K点的选择需要多少才能保证计算的大致准确性？ 之后将CONTCAR复制成POSCAR，再进行静态计算：INCAR： SYSTEM=FeCrHe ENCUT=250 ISTART=0;ICHARG=2 ISMEAR=-5 PREC=Accurate K点的选择需要多少才能保证计算的大致准确性？（并不需要太精确，因为是个人电脑进行计算，并且还是初学） 请各位大牛帮忙看一下，设置是否合理。第一次接触vasp这个软件，课题组也没有任何人模拟相关软件，全靠自己在小木虫上摸索。 可能你们的几分钟，就能给予初学者最大的帮助。你们的帮助我也将牢记，并传递给未来的初学者，谢谢！麻烦了！仅有的金币奉上 |
回复此楼» 猜你喜欢
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有15人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
梦想
已经有8人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
zhangfrank
金虫 (正式写手)
- 应助: 178 (高中生)
- 金币: 660.2
- 散金: 226
- 红花: 53
- 帖子: 891
- 在线: 344.5小时
- 虫号: 2679423
- 注册: 2013-09-25
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2016-12-30 12:54:34
温杨树
木虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 1950.4
- 帖子: 36
- 在线: 720.8小时
- 虫号: 2517772
- 注册: 2013-06-22
- 性别: GG
- 专业: 金属材料的微观结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
Sparks__Fly(蛋蛋小童鞋代发): 金币+3, 谢谢参与，小木虫鼓励积极回答问题者 2016-12-29 16:42:08
Sparks__Fly: 金币+5, ★★★很有帮助 2016-12-30 11:32:47
感谢参与，应助指数 +1
Sparks__Fly(蛋蛋小童鞋代发): 金币+3, 谢谢参与，小木虫鼓励积极回答问题者 2016-12-29 16:42:08
Sparks__Fly: 金币+5, ★★★很有帮助 2016-12-30 11:32:47
|
在结构优化之前应该先进行收敛性测试，小木虫论坛上有解释收敛性测试的帖子。ENCUT参数和KPOINTS文件中k点的个数肯定要测试。NSW=3确实太小了，手册上的说明默认值是60，EDIFF设置得太小了，可以参考别人做的文献上的EDIFF值。 在确定INCAR文件的参数之前应该仔细阅读VASP的官方手册，里面的说明更全面详细。 |
2楼2016-12-29 14:30:53
温杨树
木虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 1950.4
- 帖子: 36
- 在线: 720.8小时
- 虫号: 2517772
- 注册: 2013-06-22
- 性别: GG
- 专业: 金属材料的微观结构
3楼2016-12-29 14:31:04
zhangfrank
金虫 (正式写手)
- 应助: 178 (高中生)
- 金币: 660.2
- 散金: 226
- 红花: 53
- 帖子: 891
- 在线: 344.5小时
- 虫号: 2679423
- 注册: 2013-09-25
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
Sparks__Fly(蛋蛋小童鞋代发): 金币+2, 谢谢交流 2016-12-29 21:01:15
Sparks__Fly: 金币+7, ★★★★★最佳答案 2016-12-30 11:37:46
感谢参与，应助指数 +1
Sparks__Fly(蛋蛋小童鞋代发): 金币+2, 谢谢交流 2016-12-29 21:01:15
Sparks__Fly: 金币+7, ★★★★★最佳答案 2016-12-30 11:37:46
|
楼上讲的比较清楚了，是需要进行收敛性测试会更好。如果经验足够，也可以直接进行计算。 第一步先是结构优化：INCAR： SYSTEM=onestep-FeCrHe-1 ENCUT=250（是否合理？） ####因为设置了PREC，其实这个参数可以省略。也可以看下PREC=M时候的ENCUT推荐值是什么。 ISTART=0;ICHARG=2 ISMEAR=1;SIGMA=0.2 NSW=3（是否过小？）;IBRION=2 ###=3肯定太小了。要么不设置，用默认的60，要么设置300等，让它自己弛豫完。如果到了300 ####步还没收敛说明你的结构很可能就不能收敛了 ISIF=3 ###在你的体系中需要考虑优化晶格尝试么？ POTIM=0.2 EDIFF=1E-6;EDIFFG=-1E-3 ###EDIFFG如果取负数表示用力收敛标准，那这个数值太小了。 PREC=M POSCAR选择为Selective dynamics，仅仅固定He原子的位置（是否正确？）###既然你想看He在哪个位置稳定，我觉的你不能固定任何原子的位置吧？ K点的选择需要多少才能保证计算的大致准确性？###立方晶系的话，假设晶格常数a，建议kmesh>=30/a 之后将CONTCAR复制成POSCAR，再进行静态计算：INCAR： SYSTEM=FeCrHe ENCUT=250 ### 用默认的或者提高 ISTART=0;ICHARG=2 ISMEAR=-5 PREC=Accurate K点的选择需要多少才能保证计算的大致准确性？（并不需要太精确，因为是个人电脑进行计算，并且还是初学）###立方晶系的话，假设晶格常数a，建议kmesh>=30/a |
4楼2016-12-29 20:57:25
30