24小时热门版块排行榜    

查看: 1129  |  回复: 17

whqiaoling

银虫 (初入文坛)

[求助] CASTEP 已有3人参与

用CASTEP对Mg3N2的表面模型进行优化,构建了六个表面模型(001)、(010)、(101)面,每个面又分为Mg终止和N终止两种终止情况,这六个模型的计算参数都一样,Mg终止的模型很快就优化成功了,可是N终止的一直不成功,请大神指点一下是什么原因?该怎么解决?
主要参数设置如下
CASTEP
1.JPG


CASTEP-1
2.JPG


CASTEP-2
3.JPG


CASTEP-3
4.JPG


CASTEP-4
5.JPG


CASTEP-5
6.JPG
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

初之

木虫 (正式写手)

顶一下,我也遇到过这种情况,二氧化铈分别以Ce和O终结,结果只有以Ce终结的可以优化成功。
念念不忘,必有回响。
2楼2015-03-04 10:48:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

不知模型是否都符合计量比,另外,不同暴露面(上下表面都算)对于体系而言前线电子密度不同,适当调节density mixing可以利于收敛。
不一定挂在论坛,计算问题问题欢迎留言。
6楼2015-03-05 12:56:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

@Ellete-blue

新虫 (小有名气)

【答案】应助回帖

引用回帖:
9楼: Originally posted by whqiaoling at 2015-03-06 15:43:28
模型不满足化学计量比,是不是需要删除部分原子,一般删除哪些原子?...

建模的时候控制层的厚度就可以做到电中性,不需要删除原子,通常让厚度保持晶胞边长的整数倍,体系自然是电中性的。

[ 发自手机版 http://muchong.com/3g ]
10楼2015-03-06 16:48:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
9楼: Originally posted by whqiaoling at 2015-03-06 15:43:28
模型不满足化学计量比,是不是需要删除部分原子,一般删除哪些原子?...

楼上是正解,按照整数的厚度即可。

如果这样还是很难收敛,那么你得考虑计算的参数及精度问题,如果这样还不行,就得考虑一些数值技巧了。
不一定挂在论坛,计算问题问题欢迎留言。
12楼2015-03-06 23:13:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

@Ellete-blue

新虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
whqiaoling: 金币+30, ★★★很有帮助 2015-03-06 15:18:51
不知道你建的表面的具体结构,猜测可能在N终结的情况下,体系不是电中性的。你可以计算你最终材料的化学配比,通过减少电荷来平衡体系的电荷。同时增加表面的厚度,来降低电荷不平衡的影响。

[ 发自手机版 http://muchong.com/3g ]
3楼2015-03-04 19:42:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

@Ellete-blue

新虫 (小有名气)

【答案】应助回帖

刚才想到,最好的方法是通过调节层的厚度确保电荷是中性的。这个问题肯那个跟计算参数的设置没关系,关键是模型

[ 发自手机版 http://muchong.com/3g ]
4楼2015-03-04 19:47:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NCEPUSCR

铁虫 (小有名气)

好奇的问一下,为什么SCF设置为1000呢?系统一般默认的是50,在Dmol3中在尝试smearing之前,貌似也就提高到100。。。(小弟新手,问题可能太幼稚了,请多包涵)
5楼2015-03-04 22:33:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whqiaoling

银虫 (初入文坛)

引用回帖:
5楼: Originally posted by NCEPUSCR at 2015-03-04 22:33:38
好奇的问一下,为什么SCF设置为1000呢?系统一般默认的是50,在Dmol3中在尝试smearing之前,貌似也就提高到100。。。(小弟新手,问题可能太幼稚了,请多包涵)

因为以前计算的时候遇到过这种情况:SCF循环结束但是体系未达到基态。所以每次计算我就把SCF设置的很大,但貌似不需要这么大。
7楼2015-03-06 15:17:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whqiaoling

银虫 (初入文坛)

引用回帖:
3楼: Originally posted by @Ellete-blue at 2015-03-04 19:42:20
不知道你建的表面的具体结构,猜测可能在N终结的情况下,体系不是电中性的。你可以计算你最终材料的化学配比,通过减少电荷来平衡体系的电荷。同时增加表面的厚度,来降低电荷不平衡的影响。
...

您好,保持体系电中性就是建好表面模型之后删除几个原子吗?一般删除那些原子?
8楼2015-03-06 15:20:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whqiaoling

银虫 (初入文坛)

引用回帖:
6楼: Originally posted by 卡开发发 at 2015-03-05 12:56:46
不知模型是否都符合计量比,另外,不同暴露面(上下表面都算)对于体系而言前线电子密度不同,适当调节density mixing可以利于收敛。

模型不满足化学计量比,是不是需要删除部分原子,一般删除哪些原子?
9楼2015-03-06 15:43:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 whqiaoling 的主题更新
信息提示
请填处理意见