| 查看: 3473 | 回复: 14 | ||
[求助]
救命啊！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多种情况，求解答，不胜感激。 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有7人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有8人回复
-
博士延得我，科研能力直往上蹿
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有3人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 计算磁性多层膜的软件，ATK和VASP哪个更好用？
已经有4人回复
- VASP计算时，TOTAL-FORCE这么高是什么原因，求大神
已经有8人回复
- 感觉VASP计算缺陷存在很大问题
已经有17人回复
- 关于VASP中MAGMON的磁矩设置问题
已经有6人回复
- vasp计算中的疑惑
已经有8人回复
- VASP求助，正交结构的，磁性材料，怎样计算磁各向异性常数？？
已经有4人回复
- vasp计算 磁性时，MAGMOM的设置
已经有9人回复
- 求CASTEP 和 VASP 计算磁性问题
已经有4人回复
- VASP静态自洽计算，能量升高7eV，然后局域磁矩消失，请问是什么问题？
已经有8人回复
- vasp中计算磁矩时初始磁矩的设置参数是哪些呀，对应的各种枚举
已经有10人回复
- vasp磁性计算问题
已经有8人回复
- vasp计算磁性的问题
已经有5人回复
- 为什么VASP算出来Gd的磁矩总是偏小？
已经有14人回复
- 大家交流一下用VASP计算非共线磁性的方法和结果分析
已经有12人回复
- 关于vasp中ISMEAR参数的几个问题
已经有6人回复
- vasp 那些参数设置可以影响到磁性计算的磁矩
已经有12人回复
- vasp计算磁矩的困惑
已经有10人回复
- 【求助】vasp的运行时候报错：Error reading item 'IMAGES' from file INCAR.
已经有13人回复
- 【求助】有人用vasp算过zigzag型石墨烯纳米带的磁性么
已经有6人回复
- 【求助】xp下p4vasp打不开vasprun.xml
已经有5人回复
- 【求助】VASP怎么考虑自旋极化？
已经有14人回复
- 【求助】磁性计算
已经有16人回复
- 【求助】用vasp计算磁矩的问题
已经有13人回复
- 【求助成功】如何用VASP计算反铁磁
已经有14人回复
白菜胡萝卜 
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 3033.6
- 散金: 1231
- 红花: 19
- 帖子: 861
- 在线: 340.9小时
- 虫号: 1944808
- 注册: 2012-08-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2014-05-25 22:29:49
6楼2014-05-26 13:41:48
11楼2014-05-27 13:54:57
白菜胡萝卜
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 3033.6
- 散金: 1231
- 红花: 19
- 帖子: 861
- 在线: 340.9小时
- 虫号: 1944808
- 注册: 2012-08-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2014-05-25 00:59:38
3楼2014-05-25 08:22:54
5楼2014-05-26 08:28:00
kjsun
木虫 (正式写手)
- 应助: 34 (小学生)
- 金币: 5557.5
- 散金: 100
- 红花: 5
- 帖子: 434
- 在线: 469.8小时
- 虫号: 285070
- 注册: 2006-10-14
- 性别: GG
- 专业: 催化化学
7楼2014-05-26 13:58:19
8楼2014-05-26 14:56:54
9楼2014-05-26 15:20:55
白菜胡萝卜
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 3033.6
- 散金: 1231
- 红花: 19
- 帖子: 861
- 在线: 340.9小时
- 虫号: 1944808
- 注册: 2012-08-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2014-05-27 13:14:43
10