Prettyswan

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2919
帖子: 145
在線: 59.7小時
蟲號: 958555

[交流] 【求助成功】結(jié)構(gòu)優(yōu)化之前，結(jié)果比較好；結(jié)構(gòu)優(yōu)化之后，基態(tài)與文獻報道不符

實驗和文獻中都證實了CaFe2As2低溫基態(tài)為AFM1（我們自己命名，幾種可能的此結(jié)構(gòu)如下）

我發(fā)現(xiàn)如果直接使用實驗值而不進行結(jié)構(gòu)優(yōu)化（ISIF=3），那么計算結(jié)果的趨勢符合文獻報道，跟用wine2k的趨勢也一樣。

引用回帖:

AFM1
  free  energy TOTEN  =    -61.279719 eV
AFM2
  free  energy TOTEN  =    -61.255735 eV
AFM3
  free  energy TOTEN  =    -61.189768 eV
AFM4
  free  energy TOTEN  =    -60.938984 eV
AFM5
  free  energy TOTEN  =    -60.879493 eV
AFM6
  free  energy TOTEN  =    -61.222014 eV
AFM7
  free  energy TOTEN  =    -61.169375 eV
FM
  free  energy TOTEN  =    -60.824646 eV
NM
  free  energy TOTEN  =    -61.00932 eV
End of all

如果對實驗值進行優(yōu)化，然后使用優(yōu)化后的結(jié)構(gòu)進行計算，不僅趨勢不一樣，而且基態(tài)也不是AFM1了：

引用回帖:

AFM1
  free  energy TOTEN  =    -61.301294 eV
AFM2
  free  energy TOTEN  =    -61.301168 eV
AFM3
  free  energy TOTEN  =    -61.222014 eV
AFM4
  free  energy TOTEN  =    -61.3327854 eV
AFM5
  free  energy TOTEN  =    -61.334349 eV
AFM6
  free  energy TOTEN  =    -61.222427 eV
AFM7
  free  energy TOTEN  =    -61.222014 eV
FM
  free  energy TOTEN  =    -61.332874 eV
NM
  free  energy TOTEN  =    -61.32873 eV
End of all

這是我在結(jié)構(gòu)優(yōu)化以及總能計算中用到的對FM的INCAR文件，其他的類似

引用回帖:

#######結(jié)構(gòu)優(yōu)化的INCAR#############
System = CaFe2As2 fm:cell
PREC = Accurate
ISPIN = 2
ISMEAR = 1
ISTART = 0
ICHARG = 2
IBRION = 2
SIGMA = 0.18
EDIFF = 1E-5
EDIFFG = -1E-4
GGA = PE
NSW = 100
ENCUT = 400
ISIF = 3 #Full relaxtion
POTIM = 0.2
ISYM = 0

LORBIT = 11
LREAL = .FALSE.
NPAR = 1
LPLANE = .TRUE.
LSCALU = .FALSE.
NSIM = 4
##########靜態(tài)計算的INCAR##############
System = CaFe2As2 FM:dos
PREC = Accurate
ISPIN = 2
ISMEAR = -5
ISTART = 0
ICHARG = 2
SIGMA = 0.18
EDIFF = 1E-5
EDIFFG = -1E-4
GGA = PE
ENCUT = 400
POTIM = 0.2
ISYM = 0

LORBIT = 11
LREAL= .FALSE.
NPAR = 1
LPLANE = .TRUE.
NSIM = 4
###########非自洽計算的INCAR###########
System = CaFe2As2 fm:dos
PREC = Accurate
ISPIN = 2
ISMEAR = -5
ISTART = 1
ICHARG = 11
SIGMA = 0.18
EDIFF = 1E-5
EDIFFG = -1E-4
GGA = PE
ENCUT = 400
POTIM = 0.2
ISYM = 0

LORBIT = 11
LREAL= .FALSE.
NPAR = 1
LPLANE = .TRUE.
NSIM = 4

是結(jié)構(gòu)優(yōu)化的INCAR哪里不對？

[ Last edited by Prettyswan on 2011-3-24 at 10:12 ]

該體系存在collapse phase transition，導(dǎo)致（使用贗勢的）理論計算跟實驗值不符。

[ Last edited by Prettyswan on 2011-4-19 at 10:36 ]

回復(fù)此樓

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-03-14 11:35:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cenwanglai

榮譽版主 (知名作家)

Prettyswan(金幣+5): 2011-04-26 09:50:29

占位學習！

回復(fù)此樓

2樓2011-03-15 14:07:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

尹珺加油

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1123.9
帖子: 290
在線: 60.3小時
蟲號: 1165595

Prettyswan(金幣+5): 2011-04-26 09:50:34

學習學習！順便幫頂！

贊一下

回復(fù)此樓

3樓2011-03-19 19:57:03

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

vimtex

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 443.2
帖子: 34
在線: 17.5小時
蟲號: 481120

★
ellsaking(金幣+1): 能詳細說明下如何加入U值嗎 2011-04-05 18:55:44
Prettyswan(金幣+10): 2011-04-26 09:50:20

引用回帖:

Originally posted by Prettyswan at 2011-03-14 11:35:32:
實驗和文獻中都證實了CaFe2As2低溫基態(tài)為AFM1（我們自己命名，幾種可能的此結(jié)構(gòu)如下）

但是我計算態(tài)密度的總能卻是：
AFM1：free energy TOTEN = -61 ...

能量差得不多
遇到過類似的情況

你的這個例子，可能可以通過GGA+U調(diào)整一下，可以先試試，反正體系不復(fù)雜

贊一下(1人)

回復(fù)此樓

4樓2011-03-20 05:08:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Prettyswan

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 2919
帖子: 145
在線: 59.7小時
蟲號: 958555

引用回帖:

Originally posted by vimtex at 2011-03-20 05:08:33:
能量差得不多
遇到過類似的情況

你的這個例子，可能可以通過GGA+U調(diào)整一下，可以先試試，反正體系不復(fù)雜

能量差的不多，這就是問題所在，應(yīng)該是1、2、3、7為一組，4、5為一組，明顯分開才是。而能量基態(tài)也應(yīng)該是AFM1才對。

贊一下

回復(fù)此樓

6樓2011-03-20 10:28:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

franch

版主 (著名寫手)

Prettyswan(金幣+5): 2011-04-26 09:50:53

這個問題我也遇到過，，，有點頭疼。。

贊一下

回復(fù)此樓

7樓2011-04-05 16:59:06

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

簡單回復(fù)

lexiaoyao5樓

2011-03-20 07:17 回復(fù)

Prettyswan(金幣+5): 2011-04-26 09:50:41

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 本人考085602 化學工程專碩 +19	不知道叫什么！ 2026-03-15	21/1050	2026-03-20 20:48 by zhukairuo
[考研] 一志愿吉林大學材料學碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	3/150	2026-03-20 11:58 by 呱呱呱呱叫
[考研] 一志愿中國海洋大學，生物學，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 0703化學調(diào)劑 +10	妮妮ninicgb 2026-03-15	14/700	2026-03-19 22:59 by 學員8dgXkO
[考研] 生物學調(diào)劑招人！��！ +3	山海天嵐 2026-03-17	4/200	2026-03-19 21:34 by 怎么釋懷
[考研] 0703化學 305求調(diào)劑 +4	FY_yy 2026-03-14	4/200	2026-03-19 05:54 by anny19840123
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無懈可擊111
[考研] 354求調(diào)劑 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 297求調(diào)劑 +8	戲精丹丹丹 2026-03-17	8/400	2026-03-18 14:30 by laoshidan
[考研] 302求調(diào)劑 +10	呼呼呼。。。。 2026-03-17	10/500	2026-03-18 12:45 by Linda Hu
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 考研化學學碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 277材料科學與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運氣yunqi
[考研] 0703化學調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 288求調(diào)劑 +4	奇點0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky