24小时热门版块排行榜

返回列表

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

葬春1993

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1028
帖子: 52
在线: 77.2小时
虫号: 3270590

[交流] VASP计算矿物层状结构的尝试

本人研究白云母2：1层状结构，结构如下：

，效果图：http://muchongimg.xmcimg.com/data/b ... _1465266329_433.jpg
其中硅氧四面体中的硅有1/4被铝取代，使2：1层带负点，所以层间用阳离子平衡电荷。
在建立模型前，考虑到我的计算资源有限只能使用16个核，单胞有80多个原子，所以我采用原胞进行计算，结构如下：http://muchongimg.xmcimg.com/data/b ... _1465266457_935.jpg
建模时，由于有非等价的硅铝取代，即需要平衡电荷的引入，所以我考虑用slab模型进行尝试。
首先，需要建立无硅铝取代的结构模型，将层间阳离子去掉；
然后，进行cut-off和k-points收敛性测试，发现截断能在800eV才趋于稳定，这是为什么呢？
接着建立slab模型，我只取了一个2:1层状结构，扩成2x2x1的超胞（大概80多个原子），http://muchongimg.xmcimg.com/data/b ... _1465267498_920.jpg，考虑到c轴方向的非周期性，我人为的固定了底面的硅氧四面体层，只改变真空层厚度，进行测试，发现真空层厚度在25埃才趋稳，这样我的基础模型就建好了。
最后，就是进行铝取代硅，加入钾离子为电荷平衡离子，结构优化后个别键长比实际的白云母结构短0.05埃。
这样的模型可不可靠呢？周围没有人做这方面的，我也是新手，所以感觉很虚~~希望大家针对我建模的过程给出指导和不足！！谢谢！！