CuO晶胞的优化问题

yiranfengbai

RT，MS中优化后的CuO晶胞常数与实验值相差很大。我尝试在primitive cell之后设定Cu原子的spin为+1，但是优化后的晶胞常数还是不好。请问各位，你们平时在优化CuO这类晶胞时是怎么处理的呢？
另外，我在优化primitive cell时，没有勾选use sysmmetry ,优化之后find sysmmetry就不能回到原来的对称性了，这怎么办呢？

yiranfengbai

还有个问题，在优化表面，比如（111）面时，还要勾选use sysmmetry吗？

xylz6188

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

get-it

关于第二个问题，可以把tolerance调大

yiranfengbai

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

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

yiranfengbai

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

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

xylz6188

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

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

yiranfengbai

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

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

xylz6188

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

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

yiranfengbai

yiranfengbai

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

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

yiranfengbai

我已经做了不下30次测试，包括K-points，Energy cutoff等参数，都不是很好。

xylz6188

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

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

yiranfengbai

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

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

xylz6188

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

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