计算化学公社

 找回密码 Forget password
 注册 Register

Ni3S2打开自旋后显示磁矩为零,加较大的U值后出现磁矩但与实验值相差较远

查看数: 11173 | 评论数: 13 | 收藏 Add to favorites 2
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-2-6 21:38

正文摘要:

是这样的,我首先对Ni3S2(POSCAR是Ni9S6)进行了结构优化,没有打开自旋,没有加U。测出来的晶格常数与实验值很接近。因为有Ni,我结构优化以后便在incar里加入ispin=2,magmom设为1,再次优化,结果oszicar显示总 ...

回复 Reply

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

materials project的U是按照专门高通量做的处理,没办法一定保证可靠,如果你感兴趣可以看一下他们在网站上放的文章。特别复杂的磁性体系DFT+U甚至是杂化泛函都未必能解决很好。
pupu1616 发表于 Post on 2023-7-13 16:13:10
卡开发发 发表于 2021-2-8 00:34
接近的值我可能能帮你凑出来一个,但是没有实际意义,算出来搞不好还是非磁性的能量比较低,从Materials  ...

卡老师打扰了,请问materials project网站上给出的磁性和DOS数据是从哪里得来的呢,是否可靠呢?我算Ni9S8遇见了类似的问题,不加U,是无自旋的,和materials project一样。如果加4的U值,算出来是有自旋的。不知道怎么是合理的
卡开发发 发表于 Post on 2021-2-8 22:54:46
八十番 发表于 2021-2-8 22:42
嗯嗯好的,不过请问,评判加的U值合不合理是看总能量还是看晶格常数和实验值接不接近呢?这个我没有什么 ...

如果有条件尽可能是线性响应或其他物理的方式确定U。U对电子结构和势能面确实会有影响,但不见得是正向的。
八十番 发表于 Post on 2021-2-8 22:42:53
卡开发发 发表于 2021-2-8 22:08
自旋极化应该是可以不考虑了,U该考虑的情况下还是要考虑。

嗯嗯好的,不过请问,评判加的U值合不合理是看总能量还是看晶格常数和实验值接不接近呢?这个我没有什么经验
卡开发发 发表于 Post on 2021-2-8 22:08:17
八十番 发表于 2021-2-8 21:06
也就是说我这个东西从计算结果来看不一定非得考虑磁性和加U么。我今天又尝试改参数算了下,ISMEAR改了改 ...

自旋极化应该是可以不考虑了,U该考虑的情况下还是要考虑。
卡开发发 发表于 Post on 2021-2-8 00:34:06
本帖最后由 卡开发发 于 2021-2-8 01:59 编辑
八十番 发表于 2021-2-7 23:38
这个是我的POSCAR文件,麻烦您过目一下

接近的值我可能能帮你凑出来一个,但是没有实际意义,算出来搞不好还是非磁性的能量比较低,从Materials Project给的结果来看应该也是非磁性的。U增加引起磁矩突越的原因可能是U在增加的时候导致晶格变大,晶格变大再引起磁矩变得更大。
卡开发发 发表于 Post on 2021-2-8 00:32:04
八十番 发表于 2021-2-8 00:08
请问,在MP里搜索ni3s2会出来三种文件,具体用那种还有讲究么?我曾在知乎上看到说用mp后数字小的比较好

你取的应该是对的。一般按照形成能低的取,我们认为是稳定结构,你可以参考(Energy Above Hull / Atom),但是也因你要研究的问题而定,有时候不一定要研究最稳定的结构。
八十番 发表于 Post on 2021-2-8 00:08:11
卡开发发 发表于 2021-2-7 23:03
可以把结构发来看看,我不确定我测试的结构和你的结构是否相同。我在Materials Project上面找的mp-122010 ...

请问,在MP里搜索ni3s2会出来三种文件,具体用那种还有讲究么?我曾在知乎上看到说用mp后数字小的比较好
八十番 发表于 Post on 2021-2-7 23:38:27
卡开发发 发表于 2021-2-7 23:03
可以把结构发来看看,我不确定我测试的结构和你的结构是否相同。我在Materials Project上面找的mp-122010 ...

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

POSCAR.txt

1.91 KB, 下载次数 Times of downloads: 6

卡开发发 发表于 Post on 2021-2-7 23:03:44
本帖最后由 卡开发发 于 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可能不宜过大。
八十番 发表于 Post on 2021-2-7 22:58:18
卡开发发 发表于 2021-2-7 14:15
减小展宽增加k点试试,展宽可能会导致低估磁矩。

我把sigma设置成0.05,k点增大,磁矩算出来还是基本上约等于0.想不明白为什么会这样,应该和低估磁矩关系不大吧,低估的话应该是减小,可我现在的情况是总磁矩总是约等于0
卡开发发 发表于 Post on 2021-2-7 14:15:18
减小展宽增加k点试试,展宽可能会导致低估磁矩。

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 08:36 , Processed in 0.173573 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list