24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2555  |  回复: 32

葬春1993

新虫 (小有名气)


[交流] 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埃。
这样的模型可不可靠呢?周围没有人做这方面的,我也是新手,所以感觉很虚~~希望大家针对我建模的过程给出指导和不足!!谢谢!!

VASP计算矿物层状结构的尝试
白云母.jpg


VASP计算矿物层状结构的尝试-1
白云母 Image.jpg


VASP计算矿物层状结构的尝试-2
原胞.jpg
VASP计算矿物层状结构的尝试-3
slab模型.jpg
回复此楼

» 猜你喜欢

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

查看全部散金贴

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

wzbhit

至尊木虫 (著名写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
不太清楚你构建这样表面(slab)模型的目的是什么  就是研究什么性质?很难判断你的模型对不对

发自小木虫IOS客户端
27楼2016-06-10 22:36:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

葬春1993

新虫 (小有名气)


引用回帖:
27楼: Originally posted by wzbhit at 2016-06-10 22:36:28
不太清楚你构建这样表面(slab)模型的目的是什么  就是研究什么性质?很难判断你的模型对不对

我是想做白云母的掺杂,所以想先建立和白云母的结构模型,因为我用的VASP,不能用占有率,所以我想先用铝掺杂来建立白云母模型,再做其他原子掺杂。不知道行不行?

发自小木虫Android客户端
28楼2016-06-11 07:32:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
28楼: Originally posted by 葬春1993 at 2016-06-11 07:32:44
我是想做白云母的掺杂,所以想先建立和白云母的结构模型,因为我用的VASP,不能用占有率,所以我想先用铝掺杂来建立白云母模型,再做其他原子掺杂。不知道行不行?
...

你问题中提到了slab,大家都认为你做表面计算。
但你又说掺杂,一般大家都做体相结构掺杂。
所以 现在的问题是  你想做表面研究 还是体相研究?
如果体相,那你的slab模型肯定不对。
如果是表面掺杂,你的模型在一定程度上make sense。(还有很多问题)
BTW,800eV 太大了。你的测试肯定不对

» 本帖已获得的红花(最新10朵)

29楼2016-06-11 21:33:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

葬春1993

新虫 (小有名气)


送红花一朵
引用回帖:
29楼: Originally posted by emilyoyang at 2016-06-11 21:33:49
你问题中提到了slab,大家都认为你做表面计算。
但你又说掺杂,一般大家都做体相结构掺杂。
所以 现在的问题是  你想做表面研究 还是体相研究?
如果体相,那你的slab模型肯定不对。
如果是表面掺杂,你的模型 ...

谢谢你的回复!!我本来准备做体相掺杂,但是掺杂前的白云母结构中,原本就有一部分硅被铝取代,一个晶胞中总共四层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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
30楼: Originally posted by 葬春1993 at 2016-06-12 10:51:24
谢谢你的回复!!我本来准备做体相掺杂,但是掺杂前的白云母结构中,原本就有一部分硅被铝取代,一个晶胞中总共四层16个硅,有四个被铝取代,每层有四个取代位,组合太多,所以白云母的基础模型很难建立,我是想只 ...

大概明白你的意思了。其实你现在问题的本质就是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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)


引用回帖:
30楼: Originally posted by 葬春1993 at 2016-06-12 10:51:24
谢谢你的回复!!我本来准备做体相掺杂,但是掺杂前的白云母结构中,原本就有一部分硅被铝取代,一个晶胞中总共四层16个硅,有四个被铝取代,每层有四个取代位,组合太多,所以白云母的基础模型很难建立,我是想只 ...

另外 收敛性测试  是total energy per atom. 我没仔细算。但你的测试肯定有问题
不可能这么大。
我还是建议你先学会如何处理disorder的结构 然后再做bulk掺杂。很简单的
有问题站内信给我

» 本帖已获得的红花(最新10朵)

32楼2016-06-13 01:50:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

葬春1993

新虫 (小有名气)


送红花一朵
引用回帖:
32楼: Originally posted by emilyoyang at 2016-06-13 01:50:04
另外 收敛性测试  是total energy per atom. 我没仔细算。但你的测试肯定有问题
不可能这么大。
我还是建议你先学会如何处理disorder的结构 然后再做bulk掺杂。很简单的
有问题站内信给我...

非常非常感谢!!我会好好看!!

发自小木虫Android客户端
33楼2016-06-13 08:31:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
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): 谢谢参与
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客户端
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   回复  
葬春1993(金币+1): 谢谢参与
发自小木虫Android客户端
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): 谢谢参与
2016-06-10 22:02   回复  
葬春1993(金币+1): 谢谢参与
2016-06-10 22:08   回复  
葬春1993(金币+1): 谢谢参与
2016-06-10 22:09   回复  
葬春1993(金币+1): 谢谢参与
lvluohua26楼
2016-06-10 22:12   回复  
葬春1993(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 葬春1993 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +4 光伏爱好者 2026-05-17 7/350 2026-05-24 10:20 by 光伏爱好者
[基金申请] 西安交大新媒学院副院长用撤稿论文结题 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:16 by kudofaye
[教师之家] 论文撤稿了 +3 bjvtcliu 2026-05-24 5/250 2026-05-24 10:06 by Equinoxhua
[教师之家] 某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路! +4 zju2000 2026-05-21 5/250 2026-05-24 09:35 by songwz
[考博] 26/27申博自荐 10+4 ZXW0202 2026-05-22 9/450 2026-05-24 08:47 by bjvtcliu
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 hvkbtfonbv 2026-05-23 3/150 2026-05-24 08:01 by 9ps9vgkqva
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 pmo95bazuy 2026-05-23 7/350 2026-05-24 06:35 by fpo5ljpv91
[基金申请] 河北省自然科学基金 +5 Peterchao 2026-05-18 7/350 2026-05-23 16:57 by articlefan
[基金申请] 揭秘青基评审内幕:几个A才能顺利中标 +3 国自然国社科中 2026-05-23 4/200 2026-05-23 15:37 by 2000zf36392
[基金申请] 青B发送上会通知了吗 +5 chemBioBro 2026-05-22 7/350 2026-05-23 12:35 by zhuifengzhy
[考博] 博士申请 +3 焦晓明 2026-05-21 3/150 2026-05-23 11:26 by mlc840311
[文学芳草园] 献血感触 +7 呀呀好傻 2026-05-19 13/650 2026-05-21 20:15 by 呀呀好傻
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +14 1234567wang 2026-05-17 16/800 2026-05-21 17:58 by 脆脆的饼干
[基金申请] 国自然评分 +4 无名者登山 2026-05-20 5/250 2026-05-21 16:35 by swuq
[基金申请] 国自然上会要求 +7 无名者登山 2026-05-18 11/550 2026-05-21 15:50 by draco1987
[基金申请] 提交了我也来说说感想 +9 fummck 2026-05-20 10/500 2026-05-21 14:17 by draco1987
[基金申请] 评审有感 +15 popular289 2026-05-18 26/1300 2026-05-21 10:35 by 西葫芦炒鸡蛋
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
信息提示
请填处理意见