24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lijuanshu

铜虫 (正式写手)

应助: 1 (幼儿园)
金币: 330.2
帖子: 360
在线: 150.7小时
虫号: 1415189

[交流] 第一性原理拉伸，怎么出现两处断裂？

我研究元素偏聚在二元合金晶界处对该体系的影响，用第一性原理拉伸试验对含有晶界的体系进行拉伸，图1是不加应变的电荷密度，图2是应变为60%的电荷密度，从图二中可以明显的看出有两处断裂，并且晶界处并没有发生断裂，断裂是发生在晶界的附近，我的结果与文献上的差别很大，还希望与了解这方面的虫友们探讨下，共同解决问题，

图1CHGCAR.gif

图2CHGCAR.gif

回复此楼

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2013-04-25 16:35:52

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

acridine

木虫 (正式写手)

1ST强帖: 11
应助: 13 (小学生)
金币: 4901.6
帖子: 706
在线: 389.5小时
虫号: 305478

★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
lijuanshu: 金币+1 2013-05-09 16:02:32
franch: 金币+2, 谢谢回帖交流，， 2013-05-13 21:28:02
lijuanshu: 金币+1 2013-06-26 20:38:44

1.从你引用的句子看文章中的模型应该是没问题的,他的slab厚度够,而且他也做了收敛测试,足以证明模型的准确性.至于真空层,你之前做的表面,自然真空层要比较厚，以使得表面与镜像的距离足够远。而这里研究的是界面，即使表面受到影响，也不会影响到13个slab之外的界面的性质。我个人认为他这里不加真空层也可以的，我猜他可能是在界面掺杂Al原子后，slab厚度会改变，用真空层来起到释放slab应变的作用。你的情况应该不同吧，首先拉伸的时候有真空是不合理的，另外既然你用的是平面波法，加真空会无故引入额外的计算量，那又是何必呢？
2.我没使用过 VASP,所以对你说的软件不了解。不过实际上你的图应该是宽度相等，长度不同（是长度方向拉伸吧？）这个可能是画图时候设置的问题，也可能是数据点设置的问题。一般软件提供的工具都是方便我们查看结果，画图效果并不好。真正写作的时候都要把数据导出，再用专门的绘图软件画图
另外你的分析的角度也不太对。一般是研究断裂开始的状态，并不是等到完全断裂。你看你40%的结果，已经明显发生断键了，就没有必要再拉伸了。而且目测断裂的开端的应该就是界面原子。
不过我最终的建议还是你把模型以及拉伸过程想清楚设置明确之后再分析结果。对于你自己不确定是否可靠的结果，你分析得再清楚也是没有意义的