计算化学公社

标题: Ni3S2打开自旋后显示磁矩为零，加较大的U值后出现磁矩但与实验值相差较远 [打印本页]

作者
Author: 八十番 时间: 2021-2-6 21:38
标题: Ni3S2打开自旋后显示磁矩为零，加较大的U值后出现磁矩但与实验值相差较远
是这样的，我首先对Ni3S2（POSCAR是Ni9S6）进行了结构优化，没有打开自旋，没有加U。测出来的晶格常数与实验值很接近。因为有Ni，我结构优化以后便在incar里加入ispin=2，magmom设为1，再次优化，结果oszicar显示总磁矩为零。我便做了次静态运算，发现每个原子算出来磁矩为零。然后我在incar里加U，U值试了试，1.5、2、2.5、3、4、5、5.5、6、6.5、7、7.5、8、9。发现U值小于6的时候算出来总磁矩依然是零。当U值大于6的时候出现了磁矩，值是11左右，算下来一个Ni的磁矩有一点多，远大于参考磁矩0.66.并且晶格常数和实验值偏离较大，快达到了1埃。目前就很疑惑，有Ni的话不是应该加自旋和U么。为什么加了自旋测不出磁矩，加的U必须大于一定值才能出磁矩，而且磁矩还偏离参考值这么多，难道这种物质不需要加自旋和U？图片是，我把所有参数加进去的INCAR

作者
Author: 卡开发发 时间: 2021-2-7 14:15
减小展宽增加k点试试，展宽可能会导致低估磁矩。

作者
Author: 八十番 时间: 2021-2-7 22:58

卡开发发发表于 2021-2-7 14:15
减小展宽增加k点试试，展宽可能会导致低估磁矩。

我把sigma设置成0.05，k点增大，磁矩算出来还是基本上约等于0.想不明白为什么会这样，应该和低估磁矩关系不大吧，低估的话应该是减小，可我现在的情况是总磁矩总是约等于0

作者
Author: 卡开发发 时间: 2021-2-7 23:03
本帖最后由卡开发发于 2021-2-7 23:11 编辑

八十番发表于 2021-2-7 22:58
我把sigma设置成0.05，k点增大，磁矩算出来还是基本上约等于0.想不明白为什么会这样，应该和低估磁矩关系 ...

可以把结构发来看看，我不确定我测试的结构和你的结构是否相同。我在Materials Project上面找的mp-1220109是符合你的计量数，然后sigma=0.05，U取3的时候，算出来有磁矩。当然，只是为了测试，如果U值要取合理就还得额外讨论，另外POTIM可能不宜过大。

作者
Author: 八十番 时间: 2021-2-7 23:38

卡开发发发表于 2021-2-7 23:03
可以把结构发来看看，我不确定我测试的结构和你的结构是否相同。我在Materials Project上面找的mp-122010 ...

这个是我的POSCAR文件，麻烦您过目一下

作者
Author: 八十番 时间: 2021-2-8 00:08

卡开发发发表于 2021-2-7 23:03
可以把结构发来看看，我不确定我测试的结构和你的结构是否相同。我在Materials Project上面找的mp-122010 ...

请问，在MP里搜索ni3s2会出来三种文件，具体用那种还有讲究么？我曾在知乎上看到说用mp后数字小的比较好

作者
Author: 卡开发发 时间: 2021-2-8 00:32

八十番发表于 2021-2-8 00:08
请问，在MP里搜索ni3s2会出来三种文件，具体用那种还有讲究么？我曾在知乎上看到说用mp后数字小的比较好

你取的应该是对的。一般按照形成能低的取，我们认为是稳定结构，你可以参考（Energy Above Hull / Atom），但是也因你要研究的问题而定，有时候不一定要研究最稳定的结构。

作者
Author: 卡开发发 时间: 2021-2-8 00:34
本帖最后由卡开发发于 2021-2-8 01:59 编辑

八十番发表于 2021-2-7 23:38
这个是我的POSCAR文件，麻烦您过目一下

接近的值我可能能帮你凑出来一个，但是没有实际意义，算出来搞不好还是非磁性的能量比较低，从Materials Project给的结果来看应该也是非磁性的。U增加引起磁矩突越的原因可能是U在增加的时候导致晶格变大，晶格变大再引起磁矩变得更大。

作者
Author: 八十番 时间: 2021-2-8 21:06

卡开发发发表于 2021-2-8 00:34
接近的值我可能能帮你凑出来一个，但是没有实际意义，算出来搞不好还是非磁性的能量比较低，从Materials ...

也就是说我这个东西从计算结果来看不一定非得考虑磁性和加U么。我今天又尝试改参数算了下，ISMEAR改了改，加了LORBIT跟LASPH试了试。算出来还是0.。
然后做了次静态计算，incar是附图，算出来总磁矩和各个原子磁矩都是零。。。所以说我这个物质是不是可以放弃考虑磁矩以及加U了

作者
Author: 卡开发发 时间: 2021-2-8 22:08

八十番发表于 2021-2-8 21:06
也就是说我这个东西从计算结果来看不一定非得考虑磁性和加U么。我今天又尝试改参数算了下，ISMEAR改了改 ...

自旋极化应该是可以不考虑了，U该考虑的情况下还是要考虑。

作者
Author: 八十番 时间: 2021-2-8 22:42

卡开发发发表于 2021-2-8 22:08
自旋极化应该是可以不考虑了，U该考虑的情况下还是要考虑。

嗯嗯好的，不过请问，评判加的U值合不合理是看总能量还是看晶格常数和实验值接不接近呢？这个我没有什么经验

作者
Author: 卡开发发 时间: 2021-2-8 22:54

八十番发表于 2021-2-8 22:42
嗯嗯好的，不过请问，评判加的U值合不合理是看总能量还是看晶格常数和实验值接不接近呢？这个我没有什么 ...

如果有条件尽可能是线性响应或其他物理的方式确定U。U对电子结构和势能面确实会有影响，但不见得是正向的。

作者
Author: pupu1616 时间: 2023-7-13 16:13

卡开发发发表于 2021-2-8 00:34
接近的值我可能能帮你凑出来一个，但是没有实际意义，算出来搞不好还是非磁性的能量比较低，从Materials ...

卡老师打扰了，请问materials project网站上给出的磁性和DOS数据是从哪里得来的呢，是否可靠呢？我算Ni9S8遇见了类似的问题，不加U，是无自旋的，和materials project一样。如果加4的U值，算出来是有自旋的。不知道怎么是合理的

作者
Author: 卡开发发 时间: 2023-7-14 00:58

pupu1616 发表于 2023-7-13 16:13
卡老师打扰了，请问materials project网站上给出的磁性和DOS数据是从哪里得来的呢，是否可靠呢？我算Ni9S ...

materials project的U是按照专门高通量做的处理，没办法一定保证可靠，如果你感兴趣可以看一下他们在网站上放的文章。特别复杂的磁性体系DFT+U甚至是杂化泛函都未必能解决很好。

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