pupu1616 发表于 2023-7-13 16:13 materials project的U是按照专门高通量做的处理,没办法一定保证可靠,如果你感兴趣可以看一下他们在网站上放的文章。特别复杂的磁性体系DFT+U甚至是杂化泛函都未必能解决很好。 |
卡开发发 发表于 2021-2-8 00:34 卡老师打扰了,请问materials project网站上给出的磁性和DOS数据是从哪里得来的呢,是否可靠呢?我算Ni9S8遇见了类似的问题,不加U,是无自旋的,和materials project一样。如果加4的U值,算出来是有自旋的。不知道怎么是合理的 |
八十番 发表于 2021-2-8 22:42 如果有条件尽可能是线性响应或其他物理的方式确定U。U对电子结构和势能面确实会有影响,但不见得是正向的。 |
卡开发发 发表于 2021-2-8 22:08 嗯嗯好的,不过请问,评判加的U值合不合理是看总能量还是看晶格常数和实验值接不接近呢?这个我没有什么经验 |
八十番 发表于 2021-2-8 21:06 自旋极化应该是可以不考虑了,U该考虑的情况下还是要考虑。 |
本帖最后由 卡开发发 于 2021-2-8 01:59 编辑 八十番 发表于 2021-2-7 23:38 接近的值我可能能帮你凑出来一个,但是没有实际意义,算出来搞不好还是非磁性的能量比较低,从Materials Project给的结果来看应该也是非磁性的。U增加引起磁矩突越的原因可能是U在增加的时候导致晶格变大,晶格变大再引起磁矩变得更大。 |
八十番 发表于 2021-2-8 00:08 你取的应该是对的。一般按照形成能低的取,我们认为是稳定结构,你可以参考(Energy Above Hull / Atom),但是也因你要研究的问题而定,有时候不一定要研究最稳定的结构。 |
卡开发发 发表于 2021-2-7 23:03 请问,在MP里搜索ni3s2会出来三种文件,具体用那种还有讲究么?我曾在知乎上看到说用mp后数字小的比较好 |
卡开发发 发表于 2021-2-7 23:03 这个是我的POSCAR文件,麻烦您过目一下 |
1.91 KB, 下载次数 Times of downloads: 6
本帖最后由 卡开发发 于 2021-2-7 23:11 编辑 八十番 发表于 2021-2-7 22:58 可以把结构发来看看,我不确定我测试的结构和你的结构是否相同。我在Materials Project上面找的mp-1220109是符合你的计量数,然后sigma=0.05,U取3的时候,算出来有磁矩。当然,只是为了测试,如果U值要取合理就还得额外讨论,另外POTIM可能不宜过大。 |
卡开发发 发表于 2021-2-7 14:15 我把sigma设置成0.05,k点增大,磁矩算出来还是基本上约等于0.想不明白为什么会这样,应该和低估磁矩关系不大吧,低估的话应该是减小,可我现在的情况是总磁矩总是约等于0 |
减小展宽增加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.