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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » La掺杂BiFeO3时反铁磁性的设置问题
111/2返回列表12下一页
查看: 2191  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

future_wl

木虫 (著名写手)

[求助] La掺杂BiFeO3时反铁磁性的设置问题

我用R3c菱方表示的模型基础上建了2*2*2的超胞，然后用La取代了一个Bi，分子式变成了Bi15LaFe16O48
POSCAR 如下：
BiLaFeO3
1.0
       11.2735004425         0.0000000000         0.0000000000
        5.7346471374         9.7059587077         0.0000000000
        5.7346471374         3.2725638144         9.1376124079
    O   Bi   Fe   La
   48   15   16    1
Direct
     0.198650002         0.261649996         0.471150011
     0.471150011         0.198650002         0.261649996
     0.261649996         0.471150011         0.198650002
     0.011650000         0.448650002         0.221149996
     0.448650002         0.221149996         0.011650000
     0.221149996         0.011650000         0.448650002
     0.698650002         0.261649996         0.471150011
     0.971149981         0.198650002         0.261649996
     0.761650026         0.471150011         0.198650002
     0.511650026         0.448650002         0.221149996
     0.948650002         0.221149996         0.011650000
     0.721149981         0.011650000         0.448650002
     0.198650002         0.761650026         0.471150011
     0.471150011         0.698650002         0.261649996
     0.261649996         0.971149981         0.198650002
     0.011650000         0.948650002         0.221149996
     0.448650002         0.721149981         0.011650000
     0.221149996         0.511650026         0.448650002
     0.698650002         0.761650026         0.471150011
     0.971149981         0.698650002         0.261649996
     0.761650026         0.971149981         0.198650002
     0.511650026         0.948650002         0.221149996
     0.948650002         0.721149981         0.011650000
     0.721149981         0.511650026         0.448650002
     0.198650002         0.261649996         0.971149981
     0.471150011         0.198650002         0.761650026
     0.261649996         0.471150011         0.698650002
     0.011650000         0.448650002         0.721149981
     0.448650002         0.221149996         0.511650026
     0.221149996         0.011650000         0.948650002
     0.698650002         0.261649996         0.971149981
     0.971149981         0.198650002         0.761650026
     0.761650026         0.471150011         0.698650002
     0.511650026         0.448650002         0.721149981
     0.948650002         0.221149996         0.511650026
     0.721149981         0.011650000         0.948650002
     0.198650002         0.761650026         0.971149981
     0.471150011         0.698650002         0.761650026
     0.261649996         0.971149981         0.698650002
     0.011650000         0.948650002         0.721149981
     0.448650002         0.721149981         0.511650026
     0.221149996         0.511650026         0.948650002
     0.698650002         0.761650026         0.971149981
     0.971149981         0.698650002         0.761650026
     0.761650026         0.971149981         0.698650002
     0.511650026         0.948650002         0.721149981
     0.948650002         0.721149981         0.511650026
     0.721149981         0.511650026         0.948650002
     0.000000000         0.000000000         0.000000000
     0.250000000         0.250000000         0.250000000
     0.500000000         0.000000000         0.000000000
     0.750000000         0.250000000         0.250000000
     0.000000000         0.500000000         0.000000000
     0.500000000         0.500000000         0.000000000
     0.750000000         0.750000000         0.250000000
     0.000000000         0.000000000         0.500000000
     0.250000000         0.250000000         0.750000000
     0.500000000         0.000000000         0.500000000
     0.750000000         0.250000000         0.750000000
     0.000000000         0.500000000         0.500000000
     0.250000000         0.750000000         0.750000000
     0.500000000         0.500000000         0.500000000
     0.750000000         0.750000000         0.750000000
     0.110600002         0.110600002         0.110600002
     0.360599995         0.360599995         0.360599995
     0.610599995         0.110600002         0.110600002
     0.860599995         0.360599995         0.360599995
     0.110600002         0.610599995         0.110600002
     0.360599995         0.860599995         0.360599995
     0.610599995         0.610599995         0.110600002
     0.860599995         0.860599995         0.360599995
     0.110600002         0.110600002         0.610599995
     0.360599995         0.360599995         0.860599995
     0.610599995         0.110600002         0.610599995
     0.860599995         0.360599995         0.860599995
     0.110600002         0.610599995         0.610599995
     0.360599995         0.860599995         0.860599995
     0.610599995         0.610599995         0.610599995
     0.860599995         0.860599995         0.860599995
     0.250000000         0.750000000         0.250000000
请问大家：
1. 在进行计算的时候，这个体系还需要和BiFeO3本征体系一样设置为反铁磁的吗？
2. 如果设置反铁磁，我在INCAR 中设置 MAGMOM= 48*0   15*0   8*-6   8*6   0是否可以，需不需要在POSCAR 中将原子数写为如下形式：
    O   Bi   Fe   Fe   La
   48   15   8   8     1
3. POTCAR中Fe的赝势文件，用一个还是两个

期待大家的帮助，谢谢！！！

附上结构优化的INCAR :
SYSTEM=BiLaFeO3
ENCUT=800
ISMEAR=0
IBRION=2
ISIF=3
NSW=400
EDIFF=1E-6
EDIFFG=1E-5
EDIFFG=-0.01
ISTART=0
ICHARG=2
ISPIN=2
PREC=A
SIGMA=0.1
MAGMOM= 48*0 15*0 8*-6 8*6 0
LREAL=AUTO
LWAVE=.FALSE
LCHARG=.FALSE
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
未来就是现在
1楼 2013-04-10 11:15:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

identation

金虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
future_wl: 金币+20, ★★★很有帮助 2013-05-01 10:13:12
可以按照G型反铁磁设置初始自旋：先把POSCAR导入VESTA中查看，保证面内和面外相邻Fe之间初始磁矩反号；另外，可以考虑GGA+U计算，U典型值是4.0eV，我们组常用3.8eV；也可以不加U，但要保证计算出来的DOS有GAP
赞一下(1人) 回复此楼
2楼2013-05-01 09:56:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

future_wl

木虫 (著名写手)
引用回帖:
2楼: Originally posted by identation at 2013-05-01 09:56:11
可以按照G型反铁磁设置初始自旋：先把POSCAR导入VESTA中查看，保证面内和面外相邻Fe之间初始磁矩反号；另外，可以考虑GGA+U计算，U典型值是4.0eV，我们组常用3.8eV；也可以不加U，但要保证计算出来的DOS有GAP

感谢专家指导。
另外我想问一下：对于这种掺杂BiFeO3的体系，它的自发极化如何求呢？
我知道纯相时自发极化是通过berry phase算出顺电相和铁电相的极化强度，然后相减得到。掺杂体系的顺电相和铁电相如何确定呢？
期待您的回答！
赞一下 回复此楼
未来就是现在
3楼2013-05-01 10:13:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by future_wl at 2013-05-01 10:13:03
感谢专家指导。
另外我想问一下：对于这种掺杂BiFeO3的体系，它的自发极化如何求呢？
我知道纯相时自发极化是通过berry phase算出顺电相和铁电相的极化强度，然后相减得到。掺杂体系的顺电相和铁电相如何确定呢？ ...

也可以用berry phase计算极化强度，方法和纯相差不多
纯相也不是顺电铁电直接减，是从顺电相除法一点点一共原子，然后极化对原子位移作图，从图上得出铁电相极化强度
赞一下 回复此楼
4楼2013-05-02 01:41:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
引用回帖:
4楼: Originally posted by identation at 2013-05-02 01:41:49
也可以用berry phase计算极化强度，方法和纯相差不多
纯相也不是顺电铁电直接减，是从顺电相除法一点点一共原子，然后极化对原子位移作图，从图上得出铁电相极化强度...

感谢回复，您所说的“极化对原子位移作图”是不是原子从R-3c相到R3c相的原子位移？在掺杂体系中顺电相和铁电相也能用R-3c和R3c吗？
谢谢！
赞一下 回复此楼
未来就是现在
5楼2013-05-02 08:03:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)

【答案】应助回帖

引用回帖:
5楼: Originally posted by future_wl at 2013-05-02 08:03:35
感谢回复，您所说的“极化对原子位移作图”是不是原子从R-3c相到R3c相的原子位移？在掺杂体系中顺电相和铁电相也能用R-3c和R3c吗？
谢谢！...

这个你要找参考相的，参考的中心对称相取决于你采用的顺电相。因为berry phase的计算结果中海油arbitrary phase，这导致即使中心对称相，计算的结果可能也不为零，所以你需要扣除这个arbitrary phase。事实上，berry phase计算的极化矢量(Px,Py,Pz)是一个周期矢量,对于晶格常数那个数为a，b，c的正交结构，其周期为(Tx,Ty,Tz),其中Tx,Ty,Tz分别为把一个电子沿着x，y，z方向移动晶格常数a，b和c造成的极化强度变化，几分别为ea/abc，eb/abc。ec/abc，即约等于e/bc，e/ac，e/ab大小，你要扣除这个周期，得到确切的极化强度！
赞一下 回复此楼
6楼2013-05-02 09:48:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
引用回帖:
6楼: Originally posted by identation at 2013-05-02 09:48:08
这个你要找参考相的，参考的中心对称相取决于你采用的顺电相。因为berry phase的计算结果中海油arbitrary phase，这导致即使中心对称相，计算的结果可能也不为零，所以你需要扣除这个arbitrary phase。事实上，ber ...

有点复杂，呵呵
请教一下，自发极化强度是不是用下面公式算：
如果我用DFPT计算出玻恩有效电荷，公式里的原子位移量用铁电相坐标减去顺电相坐标，带入公式计算是否就可以？
多谢！

1111111.jpg
赞一下 回复此楼
未来就是现在
7楼2013-05-02 17:38:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)

【答案】应助回帖

引用回帖:
7楼: Originally posted by future_wl at 2013-05-02 17:38:22
有点复杂，呵呵
请教一下，自发极化强度是不是用下面公式算：
如果我用DFPT计算出玻恩有效电荷，公式里的原子位移量用铁电相坐标减去顺电相坐标，带入公式计算是否就可以？
多谢！

1111111.jpg
...

可以用波恩有效电荷计算，先得到波恩有效电荷张量，然后找出原子位移，带入公式中直接计算即可
赞一下 回复此楼
8楼2013-05-03 02:08:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
引用回帖:
8楼: Originally posted by identation at 2013-05-03 02:08:58
可以用波恩有效电荷计算，先得到波恩有效电荷张量，然后找出原子位移，带入公式中直接计算即可...

这样就不用berry phase 了吧，那和berry phase比起来精度上如何呢？
赞一下 回复此楼
未来就是现在
9楼2013-05-03 08:34:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
future_wl: 金币+10, ★有帮助, 多谢指导 2013-05-03 19:10:21
引用回帖:
9楼: Originally posted by future_wl at 2013-05-03 08:34:26
这样就不用berry phase 了吧，那和berry phase比起来精度上如何呢？...

其实是一致的，因为born charge也可以通过berry phase获得
赞一下 回复此楼
10楼2013-05-03 10:38:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 future_wl 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂 +8 调剂的考研学生 2026-03-09 8/400 2026-03-15 22:14 by Winj1e
[考研] 0703化学调剂，求各位老师收留 +7 秋有木北 2026-03-14 7/350 2026-03-15 17:30 by 小物理化学
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 313分生物学求调剂 +6 Yyt杨1 2026-03-09 8/400 2026-03-14 03:00 by JourneyLucky
[考研] 306求调剂 +4 唐薏薏 2026-03-09 4/200 2026-03-14 01:19 by JourneyLucky
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 复试调剂 +9 Copy267 2026-03-10 9/450 2026-03-13 23:45 by userper
[考研] 0703，333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] 285化工学硕求调剂（081700） +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 285求调剂 +4 ytter 2026-03-12 4/200 2026-03-13 14:48 by jxchenghu
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
信息提示
关闭
请填处理意见
关闭