计算化学公社

标题: CuO晶胞的优化问题 [打印本页]

作者
Author: yiranfengbai 时间: 2018-12-18 10:17
标题: CuO晶胞的优化问题
RT，MS中优化后的CuO晶胞常数与实验值相差很大。我尝试在primitive cell之后设定Cu原子的spin为+1，但是优化后的晶胞常数还是不好。请问各位，你们平时在优化CuO这类晶胞时是怎么处理的呢？
另外，我在优化primitive cell时，没有勾选use sysmmetry ,优化之后find sysmmetry就不能回到原来的对称性了，这怎么办呢？ (, 下载次数 Times of downloads: 71) (, 下载次数 Times of downloads: 73)

作者
Author: yiranfengbai 时间: 2018-12-18 10:18
还有个问题，在优化表面，比如（111）面时，还要勾选use sysmmetry吗？

作者
Author: xylz6188 时间: 2018-12-18 11:15
方法不同，得到的结果肯定不同，你得首先采用不同的方法计算初晶胞晶体结构，从晶胞参数、键长键角与实验数据比较，选择最接近的方法。或者参考文献，选取方法

作者
Author: get-it 时间: 2018-12-18 11:22
关于第二个问题，可以把tolerance调大

作者
Author: yiranfengbai 时间: 2018-12-18 14:44

xylz6188 发表于 2018-12-18 11:15
方法不同，得到的结果肯定不同，你得首先采用不同的方法计算初晶胞晶体结构，从晶胞参数、键长键角与实验数 ...

我时按照其他计算论文中的参数设置的，基组，smearing，K点什么的都一样，但是计算结果就是和别人不一样。我感觉关键应该时CuO的反铁磁性，就和Fe2O3一样，但是就是不知道怎么设置。以前我计算Fe2O3时也出现这种情况，后来我设置了自旋磁性，计算结果就比较好了。所以我感觉应该也时这个问题

作者
Author: yiranfengbai 时间: 2018-12-18 14:46

get-it 发表于 2018-12-18 11:22
关于第二个问题，可以把tolerance调大

tolerance=0.5的时候也是一样，还是P1,引入这个对称性之后，再Find sysmmetry，会发现对称性变成了P42/MMC。好奇怪

作者
Author: xylz6188 时间: 2018-12-19 08:08

yiranfengbai 发表于 2018-12-18 14:44
我时按照其他计算论文中的参数设置的，基组，smearing，K点什么的都一样，但是计算结果就是和别人不一样 ...

加自旋和不加那肯定区别很大的

作者
Author: yiranfengbai 时间: 2018-12-19 15:33

xylz6188 发表于 2018-12-19 08:08
加自旋和不加那肯定区别很大的

谢谢您的回复。我试着在primitive cell中选取Cu原子，然后都赋予spin+2，但是计算完之后晶胞常数还是很差，跟别人的计算文献相差很多。不知道您有没有时间测试一下，帮我试试是什么问题呢？

作者
Author: xylz6188 时间: 2018-12-19 15:48

yiranfengbai 发表于 2018-12-19 15:33
谢谢您的回复。我试着在primitive cell中选取Cu原子，然后都赋予spin+2，但是计算完之后晶胞常数还是很差 ...

你把你做的结果和实验值放出来看看，到底差别多大，这样我也无法判断啊

作者
Author: yiranfengbai 时间: 2018-12-19 15:55
(, 下载次数 Times of downloads: 66) (, 下载次数 Times of downloads: 62)

作者
Author: yiranfengbai 时间: 2018-12-19 15:57

xylz6188 发表于 2018-12-19 15:48
你把你做的结果和实验值放出来看看，到底差别多大，这样我也无法判断啊

实验值：a=4.684 , b= 3.423 , c=5.129 , β=99.540 。和我的计算相差很多。可以说，CuO是目前最让我头疼的一个晶胞，其他晶胞的优化还没有像这个这么夸张。

作者
Author: yiranfengbai 时间: 2018-12-19 15:59
我已经做了不下30次测试，包括K-points，Energy cutoff等参数，都不是很好。

作者
Author: xylz6188 时间: 2018-12-19 20:37

yiranfengbai 发表于 2018-12-19 15:57
实验值：a=4.684 , b= 3.423 , c=5.129 , β=99.540 。和我的计算相差很多。可以说，CuO是目前最让我头疼 ...

从新把所有设置复原（Tools——Settings Organizer ——reset），然后好好看看文献，重新设置参数，计算一遍。如果还是同样的情况，请换一台电脑计算一下看看。

作者
Author: yiranfengbai 时间: 2018-12-20 09:41

xylz6188 发表于 2018-12-19 20:37
从新把所有设置复原（Tools——Settings Organizer ——reset），然后好好看看文献，重新设置参数，计算 ...

按照您的方法重新试了，结果还是一样。另外找了台电脑，用同样的参数计算了一下，还是不行。难受

作者
Author: xylz6188 时间: 2018-12-20 10:35

yiranfengbai 发表于 2018-12-20 09:41
按照您的方法重新试了，结果还是一样。另外找了台电脑，用同样的参数计算了一下，还是不行。难受

我看你的版本没有校正??? 你试试GGA-PW91加上色散校正试试

作者
Author: yiranfengbai 时间: 2018-12-20 11:16

xylz6188 发表于 2018-12-20 10:35
我看你的版本没有校正??? 你试试GGA-PW91加上色散校正试试

我们这边只有MS5.0的版权，好像不能用

作者
Author: xylz6188 时间: 2018-12-20 12:15

yiranfengbai 发表于 2018-12-20 11:16
我们这边只有MS5.0的版权，好像不能用

难怪，至少用5.5的吧，5.5版本是MS最经典的版本，到时候发文章的时候你也不用写得这么详细的版权吧，我觉得你这是软件本身的问题，你换个软件试一试。MS2017真的很狗屎，算个100原子以内的结构都跑不动，我算一个高岭石跑一个月，我电脑还是40核的，我现在都准备放弃MS转VASP了，它真的支配做一个看图软件。
给你最新的：链接: https://pan.baidu.com/s/17jQqnRkf9dtLJee8WDcBXQ 提取码: szu1
以及5.5版本的：链接: https://pan.baidu.com/s/1EUwT7Nr31rrUt4513qGPuQ 提取码: aq4b

作者
Author: 1234bubble 时间: 2018-12-20 13:39

yiranfengbai 发表于 2018-12-20 11:16
我们这边只有MS5.0的版权，好像不能用

如果不想换软件或者提高MS版本的话，建议不要试了，没希望。
若干年前参加CASTEP培训，我已经问过李明宪（台湾淡江大学物理系那个），他都搞不定

作者
Author: yiranfengbai 时间: 2018-12-20 15:38

xylz6188 发表于 2018-12-20 12:15
难怪，至少用5.5的吧，5.5版本是MS最经典的版本，到时候发文章的时候你也不用写得这么详细的版权吧，我觉 ...

十分感谢，那我试试5.5的，希望结果稍微好点

作者
Author: 1234bubble 时间: 2018-12-20 16:57

yiranfengbai 发表于 2018-12-20 15:38
十分感谢，那我试试5.5的，希望结果稍微好点

应该是不用试了。7.0版本没戏，7.0之前的也应该没戏

作者
Author: yiranfengbai 时间: 2018-12-20 19:44

1234bubble 发表于 2018-12-20 16:57
应该是不用试了。7.0版本没戏，7.0之前的也应该没戏

我发现金属氧化物中，好像就CuO难搞，其他的都还好。我主要是想研究CuO(111)面上吸附小团簇，做一下反应路径，我感觉没必要太纠结这个晶胞常数吧，将就着用Dmol3优化之后的晶胞就行了吧？

作者
Author: 1234bubble 时间: 2018-12-20 20:03

yiranfengbai 发表于 2018-12-20 19:44
我发现金属氧化物中，好像就CuO难搞，其他的都还好。我主要是想研究CuO(111)面上吸附小团簇，做一下反应 ...

以我的惨痛经历友情提示：不要这样做，一是以后结构很难收敛；二是结果太不准，没法比较。
这个体系最好用VASP之类的软件算

作者
Author: yiranfengbai 时间: 2018-12-20 20:10

1234bubble 发表于 2018-12-20 20:03
以我的惨痛经历友情提示：不要这样做，一是以后结构很难收敛；二是结果太不准，没法比较。
这个体系最好 ...

还好有您的提醒

，要不然我可要费大功夫了，那我还是用MS来搞其他体系吧，如NiO，MnO2这些。，

作者
Author: 1234bubble 时间: 2018-12-20 20:12

yiranfengbai 发表于 2018-12-20 20:10
还好有您的提醒，要不然我可要费大功夫了，那我还是用MS来搞其他体系吧，如NiO，MnO2这些。，

我是折腾一年多以后，非常心痛的放弃了

作者
Author: 1234bubble 时间: 2018-12-20 20:14

yiranfengbai 发表于 2018-12-20 20:10
还好有您的提醒，要不然我可要费大功夫了，那我还是用MS来搞其他体系吧，如NiO，MnO2这些。，

先用MS做其他体系，空闲时学学VASP，慢慢就会喜欢上VASP，并且把MS降级成绘图软件了

作者
Author: yiranfengbai 时间: 2018-12-20 20:21

1234bubble 发表于 2018-12-20 20:12
我是折腾一年多以后，非常心痛的放弃了

不过我觉得这1年的折腾还是挺有用的，毕竟也是我们自己的宝贵经验，想想也挺好玩的，哈哈

欢迎光临计算化学公社 (http://ccc.keinsci.com/)