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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » （BiFeO3）GGA和LDA算出来的V0都比实验值小约10~20个立方埃？请各位高手帮我检查一下
53/1返回列表
查看: 3970  |  回复: 27
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

donggaomin

金虫 (正式写手)

[求助] （BiFeO3）GGA和LDA算出来的V0都比实验值小约10~20个立方埃？请各位高手帮我检查一下

我在优化R3c菱形结构的BiFeO3，按“PHYSICAL REVIEW B 74, 224412 2006”（见附件）中提到的理论或实验值建模并优化，文献中参数如下：用GGA，Bi:5d(10)6s(2)6p(3); Fe: 3p(6)3d(6)4s(2); O: 2s(2)2(p(4)     [圆括号中的数表示是上标]


用MS建模：用的是“161  R3c     rhombohedral " 及上图上的参数一一输入，得：

然后“save file without calculate"，从“3D Atomistic (3).cell”文档中取出以下参数作为POSCAR：
%BLOCK LATTICE_CART
       4.606442169033813       1.656627159127671       2.914118391381455
       0.000000000000000       4.895275579680102       2.914118391381455
       0.000000000000000       0.000000000000000       5.697000000000000
%ENDBLOCK LATTICE_CART

%BLOCK POSITIONS_FRAC
  O   0.5342000126838684   0.9356999993324280   0.3865000009536743
  O   0.3865000009536744   0.5342000126838684   0.9356999993324279
  O   0.9356999993324280   0.3865000009536744   0.5342000126838684
  O   1.4356999993324280   1.0342000126838682   0.8865000009536744
  O   1.0342000126838684   0.8865000009536743   1.4356999993324278
  O   0.8865000009536743   1.4356999993324280   1.0342000126838684
Fe   0.2231999933719636   0.2231999933719636   0.2231999933719636
Fe   0.7231999933719636   0.7231999933719636   0.7231999933719635
Bi   0.0000000000000000   0.0000000000000000   0.0000000000000000
Bi   0.5000000000000000   0.5000000000000000   0.5000000000000000

我得到的POSCAR：
BiFeO3_R3c
1
4.606442169033813       1.656627159127671       2.914118391381455
0.000000000000000       4.895275579680102       2.914118391381455
0.000000000000000       0.000000000000000       5.697000000000000
6 2 2
Direct
0.5342000126838684   0.9356999993324280   0.3865000009536743
0.3865000009536744   0.5342000126838684   0.9356999993324279
0.9356999993324280   0.3865000009536744   0.5342000126838684
1.4356999993324280   1.0342000126838682   0.8865000009536744
1.0342000126838684   0.8865000009536743   1.4356999993324278
0.8865000009536743   1.4356999993324280   1.0342000126838684
0.2231999933719636   0.2231999933719636   0.2231999933719636
0.7231999933719636   0.7231999933719636   0.7231999933719635
0.0000000000000000   0.0000000000000000   0.0000000000000000
0.5000000000000000   0.5000000000000000   0.5000000000000000

POTCAR: 我选(GGA)：O    Fe_pv     Bi_d   >> POTCAR   

KPOINTS:
BFO
0
MP
5 5 5
0 0 0

INCAR:
ISTART = 0
ICHARG = 2
ENCUT = 875
ISMEAR = 0; SIGMA = 0.1  
NSW = 200
IBRION = 2
ISIF = 3
EDIFF = 1E-5; EDIFFG = -1E-2
LWAVE = .FALSE.
LCHARG =.FALSE.
PREC = High

结果是：V0=114.57立方埃，而文献中的理论值为128.48，实验值为：124.60；
我改赝势：LDA，则V0=105.53，
                  PBE，则V0=114.40，
                  LSDA+U，则V0=105.72  （U值为对Fe的d电子作修正，此处U=4）
                 PBEsol+U,则V0=109.00    （U=4）
我不用ISIF=3的一次优化，而改用侯柱峰老师的有两个参量要优化的方法优化（我改成六角来优化），再用BM方程拟合得V0，PBE ， V0=117.52
                                                                    PBEsol+U, V0=111.24

总之，无论我怎么改用GGA或LDA都没有用，均远小于理论及实验值，这是为什么啊？请高手指点！
谢谢！
回复此楼

» 猜你喜欢
好好学习
1楼 2011-05-16 21:11:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

byin

铜虫 (正式写手)
引用回帖:
Originally posted by donggaomin at 2011-05-22 20:15:17:
解决了，加了SPIN = 2 体积就上去了，再调一下磁矩更好

我也是加了spin的，但不知道为什么是115.

可不可以把你的4个输入文件发我邮箱，我参考一下。
byin87@gmail.com


多谢多谢！
赞一下 回复此楼
23楼2011-05-23 14:22:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 28 个回答

identation

金虫 (著名写手)
★ ★ ★
bingmou(金币+1): 感谢交流！ 2011-05-16 21:54:40
franch(金币+2): 谢谢回帖交流 2011-05-16 22:22:40
我没计算过BiFeO3，这个我听一个老师讲过，几何优化的时候variable cell优化容易出问题，最好是手动优化晶格常数，也就是设定不同的晶格常数，优化原子为止比较总能
赞一下(2人) 回复此楼
2楼2011-05-16 21:42:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

donggaomin

金虫 (正式写手)
引用回帖:
Originally posted by identation at 2011-05-16 21:42:25:
我没计算过BiFeO3，这个我听一个老师讲过，几何优化的时候variable cell优化容易出问题，最好是手动优化晶格常数，也就是设定不同的晶格常数，优化原子为止比较总能

谢谢！
其实我已经这样做了，是按侯柱峰老师优化复杂结构的方法（有两个以上的参量要优化), 结果没多大变化，问题可能不是出在这里
赞一下 回复此楼
好好学习
3楼2011-05-16 21:57:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuhuimeiwoo

金虫 (小有名气)

【答案】应助回帖

★
donggaomin(金币+3): 谢谢您的回答，还希望您继续看看！ 2011-05-16 22:17:09
franch(金币+1): 谢谢回帖交流 2011-05-16 22:22:15
我想知道的是你们用的是同一个软件计算的吗？
赞一下(1人) 回复此楼
paper我要好多好的的paper
4楼2011-05-16 22:14:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 28 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +36 医学老男孩 2026-08-13 85/4250 2026-08-17 21:23 by lcmdlut
[基金申请] 感觉是下周放榜了 +5 angus9576 2026-08-17 9/450 2026-08-17 20:57 by ssxclkj
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 时间戳变了，能看出什么问题？ +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
关闭
请填处理意见
关闭