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, 给个红包，谢谢回帖
franch: 金币+2, 谢谢回帖交流，， 2013-05-13 21:27:39

引用回帖:

9楼: Originally posted by lijuanshu at 2013-04-27 14:30:20
非常感谢你的回复。
1.如果初始结构不加真空层的话，由于slab模型时周期性的，那个他们之间就会有相互作用啊，我一开始就比较纠结到底加不加真空层，现在更加纠结了。
2.我计算了很多有关拉伸的性质，应变从1%到 ...

1 肯定是不加真空的。不加真空时，由于周期性，你建立的超胞中是有两个晶界的（在你加真空的那个地方又出现了一个晶界）。事实上这么建模不仅是拉伸得需要，也模拟晶界的一般方法，我估计你可能刚接触界面的模拟，所以不太熟。而加真空是不合适的，引入了表面必然会影响体系的性质。
2. 我就是这个意思，应变是要逐渐加的。你通过改变z轴来作为初始结构，这在小应变下是完全可以的，但在大应变下这个初始模型就不合理的。我举个极端的例子吧，你把一根钢丝和一根橡皮筋串联绑在一起，拉伸，那么结果必然是橡皮筋先断裂。逐渐加应变就能模拟出这一过程。如果你直接加60%应变，就好比初始条件就是让钢丝和橡皮筋都被拉长了60%，那么得到的结果会是：橡皮筋没事，钢丝断了。这就足以说明初始结构无意义，会导致十分荒谬的计算结果。
另外60% 这个应变似乎太大了吧，你看看文献的应变范围取到多少