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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » La掺杂BiFeO3时反铁磁性的设置问题
111/2返回列表12下一页
查看: 2213  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 281求调剂（0805） +3 烟汐忆海 2026-03-16 8/400 2026-03-18 00:47 by 烟汐忆海
[考研] 材料专硕306英一数二 +8 z1z2z3879 2026-03-16 10/500 2026-03-18 00:01 by 星空星月
[基金申请] 被我言中：新模板不强调格式了，假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 材料与化工求调剂 +6 为学666 2026-03-16 6/300 2026-03-17 20:15 by peike
[考研] 328求调剂，英语六级551，有科研经历 +3 生物工程调剂 2026-03-16 8/400 2026-03-17 19:03 by Wangjingyue
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 工科材料085601 279求调剂 +3 困于星晨 2026-03-17 3/150 2026-03-17 14:08 by ms629
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] ［0860］321分求调剂，ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] （081700）化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
信息提示
关闭
请填处理意见
关闭