24小时热门版块排行榜

>论坛更新日志 (4201)
>文献求助 (447)
>虫友互识 (361)
>导师招生 (357)
>考博 (195)
>博后之家 (176)
>硕博家园 (160)
>休闲灌水 (148)
>论文投稿 (129)
>基金申请 (92)
>考研 (69)
>招聘信息布告栏 (60)
>教师之家 (56)
>公派出国 (45)
>绿色求助(高悬赏) (43)
>找工作 (38)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

egg19900805

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 60.5
帖子: 67
在线: 71.4小时
虫号: 2442352
注册: 2013-04-30
专业: 结构陶瓷

[求助] 固定C軸鬆弛A B軸的方法

我用PAW的勢能算BiFeO3 將它relax後(ISIF=3)
發現體積(115.5)比�?驗值(126)少了約9%

侯老師vasp軟件第23頁的複雜情況
有說到調整晶格常數?碚业阶畹偷哪芰课恢�
因此想試試看固定C軸優化A B軸 ?硎贵w積更合理

所以我將POSCAR的
Fe2 Bi2 O6
1.00000000000000
   4.7056534361214437 -0.0473600491087546 2.7444919032575892    (a軸)
   1.5539706021362183 4.4419135148306212 2.7444919082264456    (b軸)
-0.0674598255581248 -0.0473599951737402 5.4470968385121745
a、b軸分別乘以0.98~1.05?砩炜s體積

INCAR如下
SYSTEM=BiFeO3_R3C  (Fe Bi O)
ISTART = 0 ; ICHARG = 2
ISPIN = 2                            # spinpolarized calculations
ENCUT = 450
PREC = Accurate
ISMEAR = -5 ; SIGMA = 0.2
#EDIFF = 1E-5
#EDIFFG = -1E-3
NPAR = 1
MAGMOM = 6 -6 0 0 0 0 0 0 0 0
NSW = 60
IBRION = 2
ISIF = 2
LDAU = .TRUE.          #Switches on the L(S)DA+U
LDAUTYPE = 2          #Dudarev’s approach to LSDA+U
LDAUL = 2 -1 -1
LDAUU = 0 0 0
LDAUJ = 0 0 0
#RWIGS = 1.164 1.635 0.82

但是得到最低的能量卻仍然是原本relax後之晶格常數
而不是�?驗值
請問我的方法有哪一部分出錯了嗎
感謝

擷取.PNG

回复此楼

» 猜你喜欢

到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有7人回复
申请2026年博士已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有6人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

1楼2013-05-13 18:03:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

egg19900805

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 60.5
帖子: 67
在线: 71.4小时
虫号: 2442352
注册: 2013-04-30
专业: 结构陶瓷

引用回帖:

2楼: Originally posted by ApolloYang at 2013-05-13 18:39:49
1. 优化后的晶格常数往往不等于实验值。
一半的用LDA会低一些  用GGA貌似会高点。
  DFT计算对应于T＝0K，而实验室不可能T=0K。

2. 晶格常数优化，有必要算磁性嘛，  有必要加U ？

謝謝指教!!第一點讓我恍然大悟原?頊囟纫灿杏绊�

好像先前在哪裡看到磁性跟自旋一開始就要加，因為也會改變晶格的位置

謝謝你的回應@@