24小时热门版块排行榜

返回列表

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

marklau

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 247.8
散金: 10
帖子: 122
在线: 37.4小时
虫号: 619105
注册: 2008-10-07
性别: GG
专业: 凝聚态物性 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

请懂得其中缘由的给小弟讲解一下，我实在是迷糊啊。

回复此楼

» 猜你喜欢

孩子确诊有中度注意力缺陷已经有12人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
2026博士申请-功能高分子，水凝胶方向已经有6人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复
请问2026国家基金面上项目会启动申2停1吗已经有5人回复
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？已经有3人回复

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

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

marklau

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 247.8
散金: 10
帖子: 122
在线: 37.4小时
虫号: 619105
注册: 2008-10-07
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

Originally posted by 锐利的碎片 at 2010-05-05 13:04:38:
fe的磁矩呢？

加Fe的磁距也试过，不收敛。
只有加了LNONLLINEAR= .TURE. 才能收敛。

我改掉INCAR 使它更简单，加上非线性磁性计算，去掉ISMEAR 使其为默认值1。

CODE:

SYSTEM        = Fe single atom

ISTART        = 0

ISYM          = 0

LNONCOLLINEAR = .TRUE.

MAGMOM        = 0 0 3

VOSKOWN       = 1

LORBIT        = 11

KPOINTS和POSCAR POTCAR和上面一样，这样的设置是可以算的：其原子能量和BULK计算的能量差值也大约在5eV 左右，最后磁矩为4，难道Fe的单原子能量该这么计算？

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

 entering main loop

       N       E                     dE             d eps       ncg     rms          rms(c)

DAV:   1     0.256646941351E+02    0.25665E+02   -0.39633E+02    40   0.516E+01

DAV:   2     0.326313207207E+01   -0.22402E+02   -0.12133E+02    48   0.290E+01

DAV:   3    -0.314465270856E+01   -0.64078E+01   -0.53989E+01    64   0.228E+01

DAV:   4    -0.346624888441E+01   -0.32160E+00   -0.29638E+00    48   0.645E+00

DAV:   5    -0.347538439844E+01   -0.91355E-02   -0.89922E-02    40   0.774E-01    0.111E+01

DAV:   6    -0.342588678283E+01    0.49498E-01   -0.14823E+01    40   0.160E+01    0.198E+01

DAV:   7    -0.338974532228E+01    0.36141E-01   -0.21376E+01    40   0.125E+01    0.620E+00

DAV:   8    -0.329693557492E+01    0.92810E-01   -0.17043E-01    40   0.162E+00    0.549E+00

DAV:   9    -0.315637743274E+01    0.14056E+00   -0.14421E-01    40   0.141E+00    0.588E+00

DAV:  10    -0.324903074074E+01   -0.92653E-01   -0.31170E-02    40   0.674E-01    0.583E+00

DAV:  11    -0.327963051985E+01   -0.30600E-01   -0.67846E-02    40   0.103E+00    0.565E+00

DAV:  12    -0.318913147059E+01    0.90499E-01   -0.65430E-02    40   0.109E+00    0.547E+00

DAV:  13    -0.321556128242E+01   -0.26430E-01   -0.14070E-02    40   0.586E-01    0.514E+00

DAV:  14    -0.314643210543E+01    0.69129E-01   -0.13432E-01    40   0.164E+00    0.457E+00

DAV:  15    -0.316514297143E+01   -0.18711E-01   -0.61759E-03    40   0.488E-01    0.458E+00

DAV:  16    -0.312696103824E+01    0.38182E-01   -0.21595E-02    40   0.702E-01    0.406E+00

DAV:  17    -0.312710171656E+01   -0.14068E-03   -0.19803E-03    40   0.259E-01    0.403E+00

DAV:  18    -0.308253641974E+01    0.44565E-01   -0.18697E-02    40   0.627E-01    0.359E+00

DAV:  19    -0.307244883783E+01    0.10088E-01   -0.64998E-03    40   0.458E-01    0.359E+00

DAV:  20    -0.309031874113E+01   -0.17870E-01   -0.13321E-02    40   0.355E-01    0.287E+00

DAV:  21    -0.309386068450E+01   -0.35419E-02   -0.53594E-03    40   0.368E-01    0.271E+00

DAV:  22    -0.314831841382E+01   -0.54458E-01   -0.30603E-02    40   0.837E-01    0.254E+00

DAV:  23    -0.317311771316E+01   -0.24799E-01   -0.22394E-02    40   0.528E-01    0.191E+00

DAV:  24    -0.319382706095E+01   -0.20709E-01   -0.36579E-02    48   0.664E-01    0.280E+00

DAV:  25    -0.308897637713E+01    0.10485E+00   -0.90647E-02    40   0.142E+00    0.403E+00

DAV:  26    -0.318460828028E+01   -0.95632E-01   -0.31788E-02    40   0.823E-01    0.357E+00

DAV:  27    -0.319106981052E+01   -0.64615E-02   -0.21742E-02    40   0.736E-01    0.258E+00

DAV:  28    -0.322147704565E+01   -0.30407E-01   -0.27897E-02    40   0.813E-01    0.419E+00

DAV:  29    -0.323060786405E+01   -0.91308E-02   -0.16676E-02    40   0.629E-01    0.413E+00

DAV:  30    -0.325465605110E+01   -0.24048E-01   -0.49151E-02    40   0.104E+00    0.462E+00

DAV:  31    -0.318388308436E+01    0.70773E-01   -0.54351E-02    40   0.888E-01    0.388E+00

DAV:  32    -0.317862991610E+01    0.52532E-02   -0.91390E-02    40   0.969E-01    0.907E-01

DAV:  33    -0.329614729376E+01   -0.11752E+00   -0.28569E-02    40   0.653E-01    0.180E+00

DAV:  34    -0.328384359774E+01    0.12304E-01   -0.62473E-02    40   0.767E-01    0.124E+00

DAV:  35    -0.330161680537E+01   -0.17773E-01   -0.15300E-02    40   0.551E-01    0.653E-01

DAV:  36    -0.330753664342E+01   -0.59198E-02   -0.39928E-03    40   0.270E-01    0.736E-01

DAV:  37    -0.330759624964E+01   -0.59606E-04   -0.11276E-03    40   0.120E-01    0.823E-01

DAV:  38    -0.330729244196E+01    0.30381E-03   -0.52779E-04    40   0.106E-01    0.853E-01

DAV:  39    -0.330348527106E+01    0.38072E-02   -0.30225E-03    40   0.249E-01    0.450E-01

DAV:  40    -0.331145234543E+01   -0.79671E-02   -0.45436E-03    40   0.320E-01    0.158E+00

DAV:  41    -0.330686482915E+01    0.45875E-02   -0.38753E-03    40   0.231E-01    0.169E+00

DAV:  42    -0.329753989421E+01    0.93249E-02   -0.55994E-03    40   0.343E-01    0.776E-01

DAV:  43    -0.330948772079E+01   -0.11948E-01   -0.43688E-03    40   0.260E-01    0.969E-01

DAV:  44    -0.331038980037E+01   -0.90208E-03   -0.93586E-04    40   0.145E-01    0.106E+00

DAV:  45    -0.330885261332E+01    0.15372E-02   -0.61868E-04    40   0.920E-02    0.852E-01

DAV:  46    -0.330870312790E+01    0.14949E-03   -0.41874E-04    40   0.891E-02    0.610E-01

DAV:  47    -0.330963051006E+01   -0.92738E-03   -0.48464E-04    40   0.736E-02    0.643E-01

DAV:  48    -0.330963060850E+01   -0.98439E-07   -0.11608E-05    32   0.172E-02

   1 F= -.33096306E+01 E0= -.32872876E+01  d E =-.670292E-01  mag=    -0.0026     0.0020     4.0006

 writing wavefunctions