版塊導航: 正在加載中...

客戶端APP下載

論文輔導

調(diào)劑小程序

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

marklau

銀蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 247.8
散金: 10
帖子: 122
在線: 37.4小時
蟲號: 619105
注冊: 2008-10-07
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

[交流] 【求助】計算單個Fe原子，無法收斂，請教已有12人參與

新學，看了手冊上的關于單原子能量計算的地方，就找了幾個原子算算。
用Pd，C都能很好收斂，但Fe就是不行，請教怎么弄才行呢？

計算設置：
贗勢都用的PAW-GGA
INCAR: 手冊上的一點都沒有改

CODE:

ISYM = 0 

ISPIN = 2 

VOSKOWN = 1 

ISMEAR = 0 

SIGMA = 0.1 

AMIX = 0.2 

BMIX = 0.0001 

NELM = 20 

ICHARG = 1

KPOINTS

CODE:

Monkhorst PAck

0

Monkhorst Pack

1 1 1

0 0 0

POSCAR： C 用8 Fe 也用8 Pd 用10

CODE:

atom

1

8.0 0.0 0.0

0.0 8.0 0.0

0.0 0.0 8.0

cart

0.0 0.0 0.0

C的計算結果：

CODE:

 vasp.4.6.28 25Jul05 complex 

 POSCAR found :  1 types and    1 ions

 LDA part: xc-table for Ceperly-Alder, standard interpolation

 POSCAR, INCAR and KPOINTS ok, starting setup

 WARNING: wrap around errors must be expected

 FFT: planning ...           1

 reading WAVECAR

 entering main loop

       N       E                     dE             d eps       ncg     rms          rms(c)

DAV:   1     0.100381037641E+02    0.10038E+02   -0.79036E+02    12   0.256E+02

DAV:   2     0.126096636712E+01   -0.87771E+01   -0.87771E+01    24   0.203E+01

DAV:   3    -0.575116297636E-01   -0.13185E+01   -0.12479E+01    18   0.168E+01

DAV:   4    -0.149852958512E+00   -0.92341E-01   -0.73476E-01    12   0.288E+00

DAV:   5    -0.152634806381E+00   -0.27818E-02   -0.27706E-02    18   0.624E-01    0.652E-02

DAV:   6    -0.152161921161E+00    0.47289E-03   -0.44607E-04    12   0.184E-01    0.378E-02

DAV:   7    -0.152101278865E+00    0.60642E-04   -0.11383E-04    12   0.129E-01

   1 F= -.15210128E+00 E0= -.10232008E+00  d E =-.149344E+00

 writing wavefunctions

Fe的計算結果：

CODE:

 vasp.4.6.28 25Jul05 complex 

 POSCAR found :  1 types and    1 ions

 LDA part: xc-table for Ceperly-Alder, standard interpolation

 POSCAR, INCAR and KPOINTS ok, starting setup

 WARNING: wrap around errors must be expected

 FFT: planning ...           1

 reading WAVECAR

 entering main loop

       N       E                     dE             d eps       ncg     rms          rms(c)

DAV:   1     0.230243849226E+02    0.23024E+02   -0.11100E+03    16   0.315E+02

DAV:   2     0.430828794596E+00   -0.22594E+02   -0.22612E+02    32   0.609E+01

DAV:   3    -0.688701305186E+00   -0.11195E+01   -0.89686E+00    16   0.158E+01

DAV:   4    -0.802889341277E+00   -0.11419E+00   -0.86993E-01    24   0.698E+00

DAV:   5    -0.804345096997E+00   -0.14558E-02   -0.14526E-02    16   0.517E-01    0.635E+00

DAV:   6    -0.409616975908E+01   -0.32918E+01   -0.29338E+01    32   0.374E+01    0.135E+01

DAV:   7    -0.740117134359E+00    0.33561E+01   -0.46177E+00    16   0.292E+01    0.908E+00

DAV:   8    -0.652382166324E+00    0.87735E-01   -0.66724E+00    24   0.116E+01    0.625E+00

DAV:   9    -0.567736501516E+00    0.84646E-01   -0.10096E-01    16   0.137E+00    0.638E+00

DAV:  10    -0.475271277207E+00    0.92465E-01   -0.35441E-01    16   0.235E+00    0.568E+00

DAV:  11    -0.427238039158E+00    0.48033E-01   -0.20409E+00    16   0.297E+00    0.354E+00

DAV:  12    -0.573242811738E+00   -0.14600E+00   -0.33645E+00    16   0.338E+00    0.667E+00

DAV:  13    -0.581861511128E+00   -0.86187E-02   -0.86676E+00    16   0.516E+00    0.535E+00

DAV:  14    -0.811900353632E+00   -0.23004E+00   -0.73093E+00    16   0.552E+00    0.749E+00

DAV:  15    -0.974625231316E+00   -0.16272E+00   -0.16806E+01    16   0.549E+00    0.738E+00

DAV:  16    -0.915465786785E+00    0.59159E-01   -0.10641E+01    16   0.613E+00    0.723E+00

DAV:  17    -0.911197857256E+00    0.42679E-02   -0.16244E+01    16   0.568E+00    0.610E+00

DAV:  18    -0.844729299645E+00    0.66469E-01   -0.85582E+00    16   0.592E+00    0.694E+00

DAV:  19    -0.784040923533E+00    0.60688E-01   -0.13437E+01    16   0.541E+00    0.552E+00

DAV:  20    -0.758741612250E+00    0.25299E-01   -0.61771E+00    16   0.561E+00    0.664E+00

DAV:  21    -0.748458287164E+00    0.10283E-01   -0.12732E+01    16   0.520E+00    0.554E+00

DAV:  22    -0.774337121701E+00   -0.25879E-01   -0.65242E+00    16   0.566E+00    0.667E+00

DAV:  23    -0.814292862819E+00   -0.39956E-01   -0.14168E+01    16   0.525E+00    0.583E+00

DAV:  24    -0.818383056273E+00   -0.40902E-02   -0.77047E+00    16   0.582E+00    0.684E+00

DAV:  25    -0.855046073298E+00   -0.36663E-01   -0.14954E+01    16   0.540E+00    0.595E+00

DAV:  26    -0.827299653904E+00    0.27746E-01   -0.79869E+00    16   0.586E+00    0.683E+00

DAV:  27    -0.845403316473E+00   -0.18104E-01   -0.14740E+01    16   0.543E+00    0.584E+00

DAV:  28    -0.814167268451E+00    0.31236E-01   -0.76097E+00    16   0.582E+00    0.677E+00

DAV:  29    -0.831445991281E+00   -0.17279E-01   -0.14426E+01    16   0.539E+00    0.578E+00

DAV:  30    -0.806626560545E+00    0.24819E-01   -0.73934E+00    16   0.580E+00    0.674E+00

DAV:  31    -0.833127894980E+00   -0.26501E-01   -0.14441E+01    16   0.537E+00    0.582E+00

DAV:  32    -0.809186965030E+00    0.23941E-01   -0.74732E+00    16   0.581E+00    0.675E+00

DAV:  33    -0.841207002500E+00   -0.32020E-01   -0.14587E+01    16   0.538E+00    0.586E+00

DAV:  34    -0.812300294878E+00    0.28907E-01   -0.75757E+00    16   0.583E+00    0.676E+00

DAV:  35    -0.844501077548E+00   -0.32201E-01   -0.14635E+01    16   0.540E+00    0.586E+00

DAV:  36    -0.811633691714E+00    0.32867E-01   -0.75676E+00    16   0.583E+00    0.675E+00

DAV:  37    -0.843754872755E+00   -0.32121E-01   -0.14601E+01    16   0.540E+00    0.586E+00

DAV:  38    -0.809561258453E+00    0.34194E-01   -0.75176E+00    16   0.583E+00    0.674E+00

DAV:  39    -0.843400156208E+00   -0.33839E-01   -0.14577E+01    16   0.540E+00    0.586E+00

DAV:  40    -0.808495376940E+00    0.34905E-01   -0.74986E+00    16   0.583E+00    0.673E+00

DAV:  41    -0.844593949703E+00   -0.36099E-01   -0.14587E+01    16   0.540E+00    0.587E+00

DAV:  42    -0.808271591148E+00    0.36322E-01   -0.75059E+00    16   0.583E+00    0.673E+00

DAV:  43    -0.846001177321E+00   -0.37730E-01   -0.14601E+01    16   0.541E+00    0.588E+00

DAV:  44    -0.807929499882E+00    0.38072E-01   -0.75102E+00    16   0.583E+00    0.673E+00

DAV:  45    -0.846813631700E+00   -0.38884E-01   -0.14604E+01    16   0.541E+00    0.588E+00

DAV:  46    -0.807243572952E+00    0.39570E-01   -0.75043E+00    16   0.583E+00    0.672E+00

DAV:  47    -0.847297262534E+00   -0.40054E-01   -0.14601E+01    16   0.541E+00    0.589E+00

DAV:  48    -0.806507994835E+00    0.40789E-01   -0.74972E+00    16   0.584E+00    0.672E+00

DAV:  49    -0.847843379039E+00   -0.41335E-01   -0.14600E+01    16   0.542E+00    0.589E+00

DAV:  50    -0.805896505900E+00    0.41947E-01   -0.74940E+00    16   0.584E+00    0.671E+00

DAV:  51    -0.848938057707E+00   -0.43042E-01   -0.14607E+01    16   0.541E+00    0.589E+00

DAV:  52    -0.805401982771E+00    0.43536E-01   -0.74793E+00    16   0.583E+00    0.670E+00

DAV:  53    -0.846814008885E+00   -0.41412E-01   -0.14559E+01    16   0.542E+00    0.589E+00

DAV:  54    -0.804441488733E+00    0.42373E-01   -0.74536E+00    16   0.583E+00    0.670E+00

DAV:  55    -0.847741306154E+00   -0.43300E-01   -0.14574E+01    16   0.542E+00    0.589E+00

DAV:  56    -0.804995110492E+00    0.42746E-01   -0.74829E+00    16   0.584E+00    0.670E+00

DAV:  57    -0.848624031592E+00   -0.43629E-01   -0.14592E+01    16   0.543E+00    0.590E+00

DAV:  58    -0.806145226421E+00    0.42479E-01   -0.75222E+00    16   0.584E+00    0.671E+00

DAV:  59    -0.850438298920E+00   -0.44293E-01   -0.14625E+01    16   0.544E+00    0.591E+00

DAV:  60    -0.806757200440E+00    0.43681E-01   -0.75430E+00    16   0.584E+00

   1 F= -.80675720E+00 E0= -.79179809E+00  d E =-.448773E-01

 writing wavefunctions

Pd的計算結果：

CODE:

 vasp.4.6.28 25Jul05 complex 

 POSCAR found :  1 types and    1 ions

 LDA part: xc-table for Ceperly-Alder, standard interpolation

 POSCAR, INCAR and KPOINTS ok, starting setup

 WARNING: wrap around errors must be expected

 FFT: planning ...           1

 reading WAVECAR

 entering main loop

       N       E                     dE             d eps       ncg     rms          rms(c)

DAV:   1     0.628549838390E+02    0.62855E+02   -0.12860E+03    18   0.337E+02

DAV:   2     0.170291659933E+02   -0.45826E+02   -0.45379E+02    34   0.365E+01

DAV:   3    -0.104839770083E+01   -0.18078E+02   -0.18048E+02    35   0.621E+01

DAV:   4    -0.227506623423E+01   -0.12267E+01   -0.12727E+01    18   0.219E+01

DAV:   5    -0.229607774769E+01   -0.21012E-01   -0.20918E-01    20   0.188E+00    0.466E+00

DAV:   6    -0.148948556717E+01    0.80659E+00   -0.39923E+00    18   0.210E+01    0.104E+00

DAV:   7    -0.145823769572E+01    0.31248E-01   -0.21584E-01    18   0.218E+00    0.609E-01

DAV:   8    -0.143330551272E+01    0.24932E-01   -0.22237E-02    18   0.143E+00    0.754E-02

DAV:   9    -0.143328837479E+01    0.17138E-04   -0.95896E-04    27   0.247E-01

   1 F= -.14332884E+01 E0= -.14528045E+01  d E =0.585485E-01

 writing wavefunctions

請懂得其中緣由的給小弟講解一下，我實在是迷糊啊。

回復此樓

1樓 2010-05-04 22:30:11

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

marklau

銀蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 247.8
散金: 10
帖子: 122
在線: 37.4小時
蟲號: 619105
注冊: 2008-10-07
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

自己頂一下，老大們來幫小弟看看撒～～

贊一下

回復此樓

2樓2010-05-05 09:11:09

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

六198

金蟲 (正式寫手)

應助: 0 (幼兒園)
金幣: 1542.4
散金: 20
紅花: 2
帖子: 697
在線: 92小時
蟲號: 470500
注冊: 2007-11-29
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

★
marklau(金幣+1):謝謝參與

摟主既然設定了NELM＝20，自洽循環(huán)怎么會有60步阿？？

贊一下(1人)

回復此樓

窮理求真真貫縱橫皆入詠致知格物物含妙理總堪尋

3樓2010-05-05 09:48:32

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

zengxiangm

木蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 1981.2
帖子: 221
在線: 96.9小時
蟲號: 676900
注冊: 2008-12-18
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

★
marklau(金幣+1):謝謝參與

K點為何取得這么小

贊一下(1人)

回復此樓

4樓2010-05-05 12:06:43

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

fzx2008

榮譽版主 (著名寫手)

1ST強帖: 9
應助: 272 (大學生)
貴賓: 1.045
金幣: 21814.1
散金: 710
紅花: 33
沙發(fā): 1
帖子: 1808
在線: 1601小時
蟲號: 621341
注冊: 2008-10-09
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構
管轄: 第一性原理

★
marklau(金幣+1):謝謝參與

引用回帖:

Originally posted by zengxiangm at 2010-05-05 12:06:43:
K點為何取得這么小

很顯然，計算的是單個原子放在一個大盒子里。只需取一個1個k點！

贊一下(1人)

回復此樓

5樓2010-05-05 12:14:04

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

luxia0210

鐵桿木蟲 (著名寫手)

應助: 1 (幼兒園)
金幣: 6476.2
紅花: 6
帖子: 1497
在線: 329.7小時
蟲號: 877190
注冊: 2009-10-19
專業(yè): 材料物理化學

★ ★ ★
marklau(金幣+1):謝謝參與
hometownlove(金幣+2):xiexie 2010-05-05 16:01:57

我覺得問題一：你的K點太小了，不妨取6x6x6的mesh,然后你的晶格長度取的明顯偏小，F(xiàn)e的d 電子會和鄰近的Fe相互作用，你計算的記過其實是一個晶格的能量，而不是Fe原子的能量，不妨去15x15x15，還有就是沒有+U，導致在計算Fe的d電子能量石發(fā)散的，我覺得這是最關鍵的，你找一下文獻看看單個Fe原子的U值，或者你可以自己測試一下。

贊一下(2人)

回復此樓

潛心百年，為求一��！

6樓2010-05-05 12:22:36

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

銳利的碎片

木蟲 (正式寫手)

star watcher

1ST強帖: 1
應助: 136 (高中生)
金幣: 3637.1
散金: 252
紅花: 22
帖子: 988
在線: 1224.9小時
蟲號: 961933
注冊: 2010-03-05
專業(yè): 凝聚態(tài)物性 II ：電子結構

★
marklau(金幣+1):謝謝參與

fe的磁矩呢？

回復此樓

7樓2010-05-05 13:04:38

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

marklau

銀蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 247.8
散金: 10
帖子: 122
在線: 37.4小時
蟲號: 619105
注冊: 2008-10-07
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

引用回帖:

Originally posted by 六198 at 2010-05-05 09:48:32:
摟主既然設定了NELM＝20，自洽循環(huán)怎么會有60步阿？？

我也不清楚它就是走了這么多……

感覺好像沒有作用一般，我也試過設大，它還是60步后就停止了

贊一下

回復此樓

8樓2010-05-05 15:08:04

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

marklau

銀蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 247.8
散金: 10
帖子: 122
在線: 37.4小時
蟲號: 619105
注冊: 2008-10-07
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

引用回帖:

Originally posted by zengxiangm at 2010-05-05 12:06:43:
K點為何取得這么小

如5樓所說，呵呵

回復此樓

9樓2010-05-05 15:09:09

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

marklau

銀蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 247.8
散金: 10
帖子: 122
在線: 37.4小時
蟲號: 619105
注冊: 2008-10-07
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

引用回帖:

Originally posted by luxia0210 at 2010-05-05 12:22:36:
我覺得問題一：你的K點太小了，不妨取6x6x6的mesh,然后你的晶格長度取的明顯偏小，F(xiàn)e的d 電子會和鄰近的Fe相互作用，你計算的記過其實是一個晶格的能量，而不是Fe原子的能量，不妨去15x15x15，還有就是沒有+U，導 ...

謝謝回復，我試試看。

贊一下

回復此樓

10樓2010-05-05 15:13:55

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

相關版塊跳轉我要訂閱樓主 marklau 的主題更新

返回列表

普通表情龍兔虎貓高級回復 (可上傳附件)

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

[考研] 材料與化工求調(diào)劑 +6	為學666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	8/400	2026-03-17 19:03 by Wangjingyue
[考研] 268求調(diào)劑 +8	一定有學上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 304求調(diào)劑 +8	小熊joy 2026-03-14	8/400	2026-03-17 17:29 by ruiyingmiao
[考研] 一志愿南京大學，080500材料科學與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 藥學383 求調(diào)劑 +3	藥學chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 286求調(diào)劑 +3	lemonzzn 2026-03-16	5/250	2026-03-16 20:43 by lemonzzn
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 一志愿華中師范071000，325求調(diào)劑 +6	RuitingC 2026-03-12	6/300	2026-03-16 14:50 by 可淡不可忘
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7	mysdl 2026-03-11	9/450	2026-03-13 22:43 by JourneyLucky
[考研] 308求調(diào)劑 +5	是Lupa啊 2026-03-11	5/250	2026-03-13 22:13 by JourneyLucky
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學習一志愿哈工大 +4	林yaxin 2026-03-12	4/200	2026-03-13 22:04 by 星空星月
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 301求調(diào)劑 +6	Liyouyumairs 2026-03-11	6/300	2026-03-13 20:11 by JourneyLucky
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 考研調(diào)劑 +4	芬達46 2026-03-12	4/200	2026-03-13 16:04 by ruiyingmiao
[考研] 材料專碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 0817化學工程與技術考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考研] 070303一志愿西北大學學碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu

24小時熱門版塊排行榜

[交流] 【求助】計算單個Fe原子，無法收斂，請教 已有12人參與

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

[交流] 【求助】計算單個Fe原子，無法收斂，請教已有12人參與

» 本主題相關價值貼推薦，對您同樣有幫助: