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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » La掺杂BiFeO3时反铁磁性的设置问题
51/1返回列表
查看: 2194  |  回复: 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的回帖

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的回帖
查看全部 11 个回答

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的回帖
查看全部 11 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学336分求调剂 +3 zbzihdhd 2026-03-15 3/150 2026-03-16 16:44 by 我的船我的海
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考研] 一志愿华中师范071000，325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 328求调剂 +3 5201314Lsy！ 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 一志愿安徽大学材料工程专硕313分，求调剂的学校 +8 Yu先生 2026-03-10 10/500 2026-03-14 01:04 by JourneyLucky
[考研] 招收0805（材料）调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 一志愿湖师大化学289求调剂 +6 XMCMM3.14159 2026-03-10 6/300 2026-03-14 00:28 by JourneyLucky
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 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
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
信息提示
关闭
请填处理意见
关闭