| 查看: 2555 | 回复: 32 | |||
[交流]
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埃。 这样的模型可不可靠呢?周围没有人做这方面的,我也是新手,所以感觉很虚~~希望大家针对我建模的过程给出指导和不足!!谢谢!! 白云母.jpg 白云母 Image.jpg 原胞.jpg slab模型.jpg |
» 猜你喜欢
云南大学材料与能源学院解琳课题组钙钛矿博士招生
已经有7人回复
西安交大新媒学院副院长用撤稿论文结题
已经有5人回复
论文撤稿了
已经有5人回复
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路!
已经有5人回复
26/27申博自荐
已经有9人回复
售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
河北省自然科学基金
已经有7人回复
揭秘青基评审内幕:几个A才能顺利中标
已经有4人回复
青B发送上会通知了吗
已经有7人回复
» 抢金币啦!回帖就可以得到:
坐标北京,诚征女友
+2/444
中科院先进技术研究院集成电路先进封装博士后招聘
+1/287
征结婚对象,男找女,坐标【北京】【北京】【北京】
+5/265
诚邀津门师者共研
+1/85
海南医科大学肖华教授(国家级青年人才)招收博士生1名---报名从速
+4/40
【王宁利院士团队招聘】眼科学博士3人
+1/33
安徽大学俞洋洋课题组招收化学方向博士生1人
+1/33
江西理工大学稀土学院急招博士生(2026年9月入学)2名,稀土光功能材料方向,非诚勿扰
+1/31
安徽农业大学招聘师资博士后以及人才引进教师
+1/30
紧急招收2026年秋季入学博士生1名(湘潭大学 固体废弃物低碳利用湖南省工程研究中心)
+1/15
深圳大学物理与光电工程学院/深圳先进光源研究院招收2026级联合培养博士生
+1/13
上海理工大学-赵斌教授课题组招收申请考核制博士【新能源材料】
+1/11
紧急招收2026年秋季入学博士生1名(河北工大/北京科技大学联合 增材制造/生物材料)
+1/6
海南大学博士研究生第二批招生 脑空间信息学团队(图像计算方向)
+1/5
重庆新桥医院临床研究助理招聘
+1/5
北航国新院季梦奇副教授招聘博士(智能感知、立体视觉等)
+1/4
深圳理工大学核酸药物递送与干细胞工程课题组科研助理教授和博士后招聘启事
+1/4
肿瘤免疫课题组招聘-科研助理(1年左右推为博士生)
+1/3
首都医科大学药学院2026年申请考核制博士研究生招生/有机化学/药物化学 三年学制
+1/2
电力全国重点实验室双一流A类长江学者团队26年全日制学术博士1名(补招)
+1/1
27楼2016-06-10 22:36:28
|
我是想做白云母的掺杂,所以想先建立和白云母的结构模型,因为我用的VASP,不能用占有率,所以我想先用铝掺杂来建立白云母模型,再做其他原子掺杂。不知道行不行? 发自小木虫Android客户端 |
28楼2016-06-11 07:32:44
29楼2016-06-11 21:33:49
送红花一朵 |
谢谢你的回复!!我本来准备做体相掺杂,但是掺杂前的白云母结构中,原本就有一部分硅被铝取代,一个晶胞中总共四层16个硅,有四个被铝取代,每层有四个取代位,组合太多,所以白云母的基础模型很难建立,我是想只取晶胞中的一层,并做表面那层硅被铝取代,这样总共就只有4个硅原子,即使扩胞后也很少,同时,在硅铝取代后,结构电荷不平衡,考虑到原始白云母结构中层间有钾离子平衡电荷,所以我就想通过slab模型来做钾离子吸附,从而靠近原始白云母结构。这个slab模型被我当作理想的白云母结构来进行其他原子的取代掺杂!!我知道这样的模样感觉不伦不类,但是没有更好的办法了~~ 对于截断能测试,我是以原胞无硅铝取代去钾,共40个原子进行测试,测试结果是这样的: 400 -294.426485 450 -294.015464 500 -293.903562 550 -293.928305 600 -293.997559 650 -294.060051 700 -294.108222 750 -294.131795 800 -294.146404 850 -294.149370 900 -294.149391 到800后相差<0.01,所以选的800。 你看有没有问题!! 发自小木虫Android客户端 |
30楼2016-06-12 10:51:24
★
小木虫: 金币+0.5, 给个红包,谢谢回帖
小木虫: 金币+0.5, 给个红包,谢谢回帖
|
大概明白你的意思了。其实你现在问题的本质就是ordering,因为白云母是个disorder的structure。 关于disorder结构如何order 有很多算法(or package)可以帮你搞定。手动的排列组合 是很不靠谱的一种做法。 你可以看一下这篇文章 https://jcheminf.springeropen.co ... 6/s13321-016-0129-3 这里面提到了 目前关于disorder结构处理的所有方法。 另外,推荐你可以学习一下pymatgen(http://pymatgen.org/) 他可以很方便的处理disorder的结构(OrderDisorderedStructureTransformation, EnumerateStructureTransformation) |
31楼2016-06-13 01:48:34
32楼2016-06-13 01:50:04
33楼2016-06-13 08:31:30
简单回复
ieem2楼
2016-06-07 10:48
回复
葬春1993(金币+1): 谢谢参与
祝福
2016-06-07 10:52
回复
葬春1993(金币+1): 谢谢参与
祝福祝福祝福祝福祝福祝福祝福祝福
tzynew4楼
2016-06-07 17:35
回复
葬春1993(金币+1): 谢谢参与
2016-06-07 18:04
回复
葬春1993(金币+1): 谢谢参与
2016-06-07 18:13
回复
葬春1993(金币+1): 谢谢参与
2016-06-07 18:44
回复
葬春1993(金币+1): 谢谢参与
2016-06-07 19:06
回复
葬春1993(金币+1): 谢谢参与
寻找cock9楼
2016-06-07 19:07
回复
葬春1993(金币+1): 谢谢参与
祝福 发自小木虫IOS客户端
skyish10楼
2016-06-07 19:14
回复
葬春1993(金币+1): 谢谢参与
一 发自小木虫IOS客户端
挑战0111楼
2016-06-07 19:33
回复
葬春1993(金币+1): 谢谢参与
赞 发自小木虫Android客户端
BorjiginSMY12楼
2016-06-07 19:34
回复
葬春1993(金币+1): 谢谢参与
好 发自小木虫IOS客户端
nono200913楼
2016-06-09 16:38
回复
葬春1993(金币+1): 谢谢参与
shwkfd14楼
2016-06-09 16:49
回复
葬春1993(金币+1): 谢谢参与
祝福 发自小木虫IOS客户端
小蜗牛_15楼
2016-06-09 16:57
回复
jixiangbao16楼
2016-06-09 18:38
回复
葬春1993(金币+1): 谢谢参与
fdzjm17楼
2016-06-09 19:55
回复
葬春1993(金币+1): 谢谢参与
祝福 发自小木虫IOS客户端
shenrenren18楼
2016-06-09 22:04
回复
葬春1993(金币+1): 谢谢参与
caledonian19楼
2016-06-10 03:30
回复
葬春1993(金币+1): 谢谢参与
Jameskelly20楼
2016-06-10 04:05
回复
葬春1993(金币+1): 谢谢参与
祝福
gorgan21楼
2016-06-10 09:39
回复
葬春1993(金币+1): 谢谢参与
祝福 [ 发自手机版 http://muchong.com/3g ]
psylhh22楼
2016-06-10 21:58
回复
葬春1993(金币+1): 谢谢参与
yangdw12345923楼
2016-06-10 22:02
回复
葬春1993(金币+1): 谢谢参与
yangdw12345624楼
2016-06-10 22:08
回复
葬春1993(金币+1): 谢谢参与
sz_kebaoyuan25楼
2016-06-10 22:09
回复
葬春1993(金币+1): 谢谢参与
lvluohua26楼
2016-06-10 22:12
回复
葬春1993(金币+1): 谢谢参与













回复此楼