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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】计算单个Fe原子，无法收敛，请教
52/1返回列表
查看: 2508  |  回复: 19
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

marklau

银虫 (小有名气)

[交流] 【求助】计算单个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发消息 送TA红花 TA的回帖

marklau

银虫 (小有名气)
引用回帖:
Originally posted by hlshi at 2010-05-05 15:48:58:
加上一个参数LMAXMIX=4 试一下

我将ISMEAR 换成1，就可以收敛了。谢谢大虾～～
CODE:
vasp.4.6.28 25Jul05 complex
POSCAR found :  1 types and    1 ions
LDA part: xc-table for Ceperly-Alder, Vosko type interpolation para-ferro
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: wrap around errors must be expected
FFT: planning ...           1
reading WAVECAR
WARNING: chargedensity file is incomplete
entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.228579425515E+02    0.22858E+02   -0.94363E+02    36   0.253E+02
DAV:   2     0.323294218380E+00   -0.22535E+02   -0.21933E+02    70   0.479E+01
DAV:   3    -0.177611603020E+01   -0.20994E+01   -0.13943E+01    38   0.172E+01
DAV:   4    -0.184327818011E+01   -0.67162E-01   -0.61485E-01    38   0.646E+00
DAV:   5    -0.184732976998E+01   -0.40516E-02   -0.37287E-02    44   0.807E-01    0.209E+01
DAV:   6    -0.302844171571E+01   -0.11811E+01   -0.67150E+00    36   0.219E+01    0.459E+00
DAV:   7    -0.324208515997E+01   -0.21364E+00   -0.44398E+00    36   0.115E+01    0.730E+00
DAV:   8    -0.317459475747E+01    0.67490E-01   -0.38011E-01    36   0.614E+00    0.497E+00
DAV:   9    -0.322182782026E+01   -0.47233E-01   -0.15701E+00    36   0.602E+00    0.748E+00
DAV:  10    -0.260313193399E+01    0.61870E+00   -0.39587E-02    36   0.208E+00    0.563E+00
DAV:  11    -0.272945384789E+01   -0.12632E+00   -0.10652E+00    36   0.586E+00    0.509E+00
DAV:  12    -0.315110478641E+01   -0.42165E+00   -0.27704E-01    36   0.293E+00    0.192E+00
DAV:  13    -0.320200415799E+01   -0.50899E-01   -0.61462E-02    36   0.209E+00    0.100E+00
DAV:  14    -0.323922376333E+01   -0.37220E-01   -0.13802E-02    36   0.813E-01    0.217E+00
DAV:  15    -0.321361156158E+01    0.25612E-01   -0.47182E-02    36   0.117E+00    0.410E+00
DAV:  16    -0.319325145724E+01    0.20360E-01   -0.26023E-01    36   0.186E+00    0.261E+00
DAV:  17    -0.326741792662E+01   -0.74166E-01   -0.93128E-02    36   0.107E+00    0.624E-01
DAV:  18    -0.325784243572E+01    0.95755E-02   -0.66986E-03    36   0.754E-01    0.178E+00
DAV:  19    -0.323449325840E+01    0.23349E-01   -0.25929E-03    36   0.485E-01    0.250E+00
DAV:  20    -0.320865773067E+01    0.25836E-01   -0.83411E-04    36   0.201E-01    0.248E+00
DAV:  21    -0.320828740813E+01    0.37032E-03   -0.30582E-04    45   0.891E-02    0.238E+00
DAV:  22    -0.320919420029E+01   -0.90679E-03   -0.16676E-05    27   0.345E-02    0.235E+00
DAV:  23    -0.321027147707E+01   -0.10773E-02   -0.24401E-05    29   0.382E-02    0.239E+00
DAV:  24    -0.320667821151E+01    0.35933E-02   -0.38039E-03    36   0.731E-01    0.898E-01
DAV:  25    -0.324283163783E+01   -0.36153E-01   -0.25141E-03    36   0.335E-01    0.775E-01
DAV:  26    -0.325147233929E+01   -0.86407E-02   -0.34432E-04    36   0.122E-01    0.759E-01
DAV:  27    -0.326560502424E+01   -0.14133E-01   -0.67164E-03    38   0.618E-01    0.219E+00
DAV:  28    -0.323477424056E+01    0.30831E-01   -0.18790E-03    36   0.308E-01    0.226E+00
DAV:  29    -0.323448768723E+01    0.28655E-03   -0.40494E-03    36   0.431E-01    0.328E+00
DAV:  30    -0.321090439600E+01    0.23583E-01   -0.18260E-03    36   0.261E-01    0.317E+00
DAV:  31    -0.320958193682E+01    0.13225E-02   -0.48577E-05    26   0.541E-02    0.310E+00
DAV:  32    -0.321086879179E+01   -0.12869E-02   -0.12477E-05    18   0.186E-02    0.312E+00
DAV:  33    -0.321632280987E+01   -0.54540E-02   -0.49685E-04    35   0.213E-01    0.366E+00
DAV:  34    -0.320690152758E+01    0.94213E-02   -0.36791E-03    36   0.311E-01    0.354E+00
DAV:  35    -0.321346015566E+01   -0.65586E-02   -0.11672E-02    52   0.414E-01    0.365E+00
DAV:  36    -0.320802785997E+01    0.54323E-02   -0.51965E-04    36   0.140E-01    0.352E+00
DAV:  37    -0.321491374965E+01   -0.68859E-02   -0.41715E-03    36   0.302E-01    0.359E+00
DAV:  38    -0.321360563769E+01    0.13081E-02   -0.32780E-03    52   0.225E-01    0.336E+00
DAV:  39    -0.321430774096E+01   -0.70210E-03   -0.37629E-03    36   0.304E-01    0.285E+00
DAV:  40    -0.322043138211E+01   -0.61236E-02   -0.53157E-03    36   0.460E-01    0.216E+00
DAV:  41    -0.322924632608E+01   -0.88149E-02   -0.56379E-03    36   0.678E-01    0.794E-01
DAV:  42    -0.326425783331E+01   -0.35012E-01   -0.57882E-03    36   0.602E-01    0.687E-01
DAV:  43    -0.327409304375E+01   -0.98352E-02   -0.38324E-03    36   0.449E-01    0.338E-01
DAV:  44    -0.327453412663E+01   -0.44108E-03   -0.12725E-03    36   0.198E-01    0.680E-01
DAV:  45    -0.327297514671E+01    0.15590E-02   -0.46248E-04    36   0.186E-01    0.633E-01
DAV:  46    -0.327239895959E+01    0.57619E-03   -0.16750E-04    27   0.658E-02    0.507E-01
DAV:  47    -0.327303232208E+01   -0.63336E-03   -0.85596E-04    28   0.122E-01    0.360E-01
DAV:  48    -0.327355242643E+01   -0.52010E-03   -0.10669E-03    36   0.136E-01    0.280E-01
DAV:  49    -0.327366076535E+01   -0.10834E-03   -0.13545E-03    36   0.125E-01    0.191E-01
DAV:  50    -0.327385292927E+01   -0.19216E-03   -0.41903E-04    27   0.812E-02    0.128E-01
DAV:  51    -0.327381373369E+01    0.39196E-04   -0.11043E-04    27   0.577E-02
   1 F= -.32738137E+01 E0= -.32621087E+01  d E =-.351151E-01  mag=     4.0000
writing wavefunctions
赞一下 回复此楼
15楼2010-05-05 16:36:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

marklau

银虫 (小有名气)
自己顶一下，老大们来帮小弟看看撒～～
赞一下 回复此楼
2楼2010-05-05 09:11:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

六198

金虫 (正式写手)
★
marklau(金币+1):谢谢参与
搂主既然设定了NELM＝20，自洽循环怎么会有60步阿？？
赞一下(1人) 回复此楼
穷理求真真贯纵横皆入咏致知格物物含妙理总堪寻
3楼2010-05-05 09:48:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zengxiangm

木虫 (小有名气)
★
marklau(金币+1):谢谢参与
K点为何取得这么小
赞一下(1人) 回复此楼
4楼2010-05-05 12:06:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭