应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 救命啊！VASP计算磁性出错
151/2返回列表12下一页
查看: 3365  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

junfeiyan

新虫 (小有名气)

[求助] 救命啊！VASP计算磁性出错 已有2人参与

小弟在复现一篇关于PTO缺陷诱导磁性的文章（Vacancy-driven ferromagnetism in ferroelectric PbTiO3），我是先优化PTO单胞，然后组建超胞，2*2*2，主要参数和这篇APL是一样的，附我的INCAR
Name of system
   SYSTEM    =     PbTiO3

Start parameter
   ISTART    =     0
   INIWAV    =     1
   ICHARG    =     2

Write flags
   NWRITE    =     2
   LWAVE     =     T
   LCHARG    =     T

Parallelization
   NPAR      =     4
   LPLANE    =     T
   LSCALU    =     F
   LSCALAPACK =    F

Electronic Relaxation 1
   ENCUT     =     500.00  eV
   EDIFF     =     1.0E-06  eV
   PREC      =     High
   NELM      =     200
   NELMIN    =     2

Electronic Relaxation 2
   ALGO      =     Normal
   LDIAG     =     T
   LREAL     =     .FALSE.

Ionic Relaxation
   #NSW       =     150
   #IBRION    =     2
   #ISIF      =     3
   #EDIFFG    =     -0.01
   #POTIM     =     0.5
   #ISYM      =     2
   #IWAVPR    =     1
   #LCORR     =     T

DOS related values
   ISMEAR    =     0
   SIGMA     =     0.025
   #NEDOS     =     3000
#   NBANDS    =     3958

Magnetic parameter
   ISPIN     =     2
   LORBIT    =     11
Exchange-correlation functional
   GGA       =     PE
   VOSKOWN = 1
   MAGMOM = 8*1 7*1 24*1
我分别用ISIF=2(GGA),ISIF=3(GGA),ISIF=3(GGA+U)去优化晶胞，然后自洽计算磁性，但结果和文献的差距很大
比如，我用和文献一样的GGA，那么我用ISIF=3计算的话，当初始磁矩为2或3，5的时候，总磁矩为0.88，文献里Tii缺陷诱导的磁矩为2.223，我的结果小很多，但是，自旋电荷密度图像是差不多的然后我把每个原子的初始磁矩设为9，结果很离谱了，磁矩为3.13，而且缺陷附近的原子磁矩大，远的磁矩小，明显不符合事实。ISIF=2的话，和ISIF的情况类似，只是磁矩偏小，才0.3，加U也很离谱，我试着把初始磁矩设为0.5和9两个极端，我发现结果一样，磁矩为3.1，但是结果肯定是错误的，因为对比文献，磁矩该大的地方小，该小的地方大。希望大神出来解救小弟，我已经捣鼓一天了，试了N多种情况，求解答，不胜感激。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-05-24 22:42:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

白菜胡萝卜

木虫 (正式写手)
★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-05-27 09:17:38
引用回帖:
3楼: Originally posted by junfeiyan at 2014-05-25 08:22:54
你是在帮我顶贴还是套取金币呢...

我昨天在手机上回复的，所以全是乱码。我是在说，你多做几遍静态计算看看磁矩会不会降下去。
赞一下(1人) 回复此楼
4楼2014-05-25 22:29:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

junfeiyan

新虫 (小有名气)
引用回帖:
5楼: Originally posted by 零下1℃ at 2014-05-26 08:28:00
PbTiO3？？？？初始磁矩MAGMOM = 8*1 7*1 24*1?????LZ你是不是把O原子也设初始磁矩啦？这错误是不是太明显啦，氧原子可是没有磁性的啊

我在研究的是缺陷诱导磁性，Ti缺陷会诱导周围氧原子的磁性，此外晚生想问下，VASP磁矩不是从大到小的一个优化过程吗，如果没有磁矩的原子给了一定得初始磁矩，难道最终不会吧磁矩优化为0吗
赞一下(1人) 回复此楼
6楼2014-05-26 13:41:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

junfeiyan

新虫 (小有名气)
引用回帖:
10楼: Originally posted by 白菜胡萝卜 at 2014-05-27 13:14:43
O原子设置初始磁性是没有问题的。首先，O原子自身是P轨道价电子，自身不会形成磁性；其次，在形成化合物之后，3d电子或者4f电子会与O原子的p电子杂化成键，这样会使得O元素的p电子性质改变，能不能产生磁性就不是单 ...

前辈，我经过与一个师兄的讨论，他跟我说，按照洪特定则和泡立不相容原理的电子排布，去给定初始磁矩，比如Pb电子排布4f145d106s26p2，那么初始磁矩为2，氧本身按电子排布没有磁性，不设也没关系，加入诱导出来的话，最终也会表现出来，我刚接触磁性，即得以为前辈的话，是vasp磁矩是从大到小优化的，给给的大了会优化到合理的值，你给笑了，可能就出错，我现在也不知道，初始磁矩该怎么给，是由特定的规则，还是我们要不断的尝试不同的大小，以得出不同的结果去对比。
赞一下(1人) 回复此楼
11楼2014-05-27 13:54:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

白菜胡萝卜

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
liliangfang: 出错 2014-05-25 15:06:47
你多???几次?????计算看看。估计就趋于一致了

[ ????С??????? ]
赞一下 回复此楼
2楼2014-05-25 00:59:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

junfeiyan

新虫 (小有名气)
???????:
2?: Originally posted by ??????? at 2014-05-25 00:59:38
你多???几次?????计算看看。估计就趋于一致了

????????????????????????
赞一下 回复此楼
3楼2014-05-25 08:22:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

零下1℃

铁虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
franch: 金币+2, 谢谢回帖交流。。。 2014-05-27 09:17:33
PbTiO3？？？？初始磁矩MAGMOM = 8*1 7*1 24*1?????LZ你是不是把O原子也设初始磁矩啦？这错误是不是太明显啦，氧原子可是没有磁性的啊
赞一下 回复此楼
5楼2014-05-26 08:28:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kjsun

木虫 (正式写手)
★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-05-27 09:17:52
VASP 计算磁性体系有缺陷，通过调整MAGMOM能解决一部分问题，但是也不能完全保证好用。
赞一下 回复此楼
7楼2014-05-26 13:58:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

零下1℃

铁虫 (正式写手)
★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-05-27 09:17:57
引用回帖:
6楼: Originally posted by junfeiyan at 2014-05-26 13:41:48
我在研究的是缺陷诱导磁性，Ti缺陷会诱导周围氧原子的磁性，此外晚生想问下，VASP磁矩不是从大到小的一个优化过程吗，如果没有磁矩的原子给了一定得初始磁矩，难道最终不会吧磁矩优化为0吗...

你研究的方向理论我不懂，不能妄下定论。。但是做了有一段时间的磁性计算，还是第一次听说没有磁性的原子也可以赋予初始磁矩。个人观点，这样做不是很合适。
赞一下 回复此楼
8楼2014-05-26 14:56:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

junfeiyan

新虫 (小有名气)
引用回帖:
8楼: Originally posted by 零下1℃ at 2014-05-26 14:56:54
你研究的方向理论我不懂，不能妄下定论。。但是做了有一段时间的磁性计算，还是第一次听说没有磁性的原子也可以赋予初始磁矩。个人观点，这样做不是很合适。...

多谢前辈，可能确实是我的理论和经验不足，我也是刚接触磁性计算，多谢前辈指导。
赞一下 回复此楼
9楼2014-05-26 15:20:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

白菜胡萝卜

木虫 (正式写手)

【答案】应助回帖

引用回帖:
8楼: Originally posted by 零下1℃ at 2014-05-26 14:56:54
你研究的方向理论我不懂，不能妄下定论。。但是做了有一段时间的磁性计算，还是第一次听说没有磁性的原子也可以赋予初始磁矩。个人观点，这样做不是很合适。...

O原子设置初始磁性是没有问题的。首先，O原子自身是P轨道价电子，自身不会形成磁性；其次，在形成化合物之后，3d电子或者4f电子会与O原子的p电子杂化成键，这样会使得O元素的p电子性质改变，能不能产生磁性就不是单纯分析p轨道电子能不能产生磁性那么简单了。不然还要第一性原理计算干嘛。
赞一下 回复此楼
10楼2014-05-27 13:14:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 junfeiyan 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭