【求助】计算单个Fe原子，无法收敛，请教 - 第一性原理 - Vasp&MedeA

版块导航: 正在加载中...

客户端APP下载

论文辅导

申博辅导

登录注册

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

24小时热门版块排行榜

返回列表

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

[ Last edited by marklau on 2010-5-5 at 15:53 ]

赞一下

回复此楼

11楼2010-05-05 15:16:01

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

hlshi

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 39.2
散金: 11
红花: 1
帖子: 360
在线: 89.1小时
虫号: 467913
注册: 2007-11-27
性别: GG
专业: 凝聚态物性 II ：电子结构

★ ★
marklau(金币+1):谢谢参与
zzy870720z(金币+1):感谢建议 2010-05-10 08:01:48

加上一个参数LMAXMIX=4 试一下

引用回帖:

Originally posted by marklau at 2010-05-04 22:30:11:
新学，看了手册上的关于单原子能量计算的地方，就找了几个原子算算。
用Pd，C都能很好收敛，但Fe就是不行，请教怎么弄才行呢？

计算设置：
赝势都用的PAW-GGA
INCAR: 手册上的一点都没有改
[code]
ISYM = ...

赞一下(2人)

回复此楼

12楼2010-05-05 15:48:58

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

marklau

银虫 (小有名气)

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

引用回帖:

Originally posted by hlshi at 2010-05-05 15:48:58:
加上一个参数LMAXMIX=4 试一下

好的，我试试看～～

赞一下

回复此楼

13楼2010-05-05 15:55:27

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

marklau

银虫 (小有名气)

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

引用回帖:

Originally posted by hlshi at 2010-05-05 15:48:58:
加上一个参数LMAXMIX=4 试一下

果然是高手，一出手就是这么拉风。

这下连NELM＝20都好使了算到20步就“嘎”，停了。

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.228549245353E+02    0.22855E+02   -0.94363E+02    36   0.253E+02

DAV:   2     0.295334820715E+00   -0.22560E+02   -0.21934E+02    70   0.479E+01

DAV:   3    -0.182618532218E+01   -0.21215E+01   -0.14135E+01    38   0.173E+01

DAV:   4    -0.189912810496E+01   -0.72943E-01   -0.65897E-01    38   0.658E+00

DAV:   5    -0.190337083148E+01   -0.42427E-02   -0.40454E-02    44   0.824E-01    0.209E+01

DAV:   6    -0.309485403745E+01   -0.11915E+01   -0.66578E+00    36   0.219E+01    0.438E+00

DAV:   7    -0.325813761933E+01   -0.16328E+00   -0.39052E+00    36   0.114E+01    0.709E+00

DAV:   8    -0.317855284913E+01    0.79585E-01   -0.38526E-01    36   0.618E+00    0.467E+00

DAV:   9    -0.292741746749E+01    0.25114E+00   -0.11151E+00    36   0.786E+00    0.655E+00

DAV:  10    -0.319609665261E+01   -0.26868E+00   -0.40323E-01    36   0.651E+00    0.492E+00

DAV:  11    -0.301861264412E+01    0.17748E+00   -0.15888E-01    36   0.390E+00    0.391E+00

DAV:  12    -0.298525732405E+01    0.33355E-01   -0.13364E-02    36   0.683E-01    0.360E+00

DAV:  13    -0.295361280139E+01    0.31645E-01   -0.74603E-01    36   0.407E+00    0.491E+00

DAV:  14    -0.307903200415E+01   -0.12542E+00   -0.11658E+00    36   0.569E+00    0.413E+00

DAV:  15    -0.308863987758E+01   -0.96079E-02   -0.98570E-03    36   0.517E-01    0.411E+00

DAV:  16    -0.308889814071E+01   -0.25826E-03   -0.84297E-04    36   0.226E-01    0.409E+00

DAV:  17    -0.308509568675E+01    0.38025E-02   -0.25643E-04    27   0.108E-01    0.408E+00

DAV:  18    -0.305561209154E+01    0.29484E-01   -0.21902E-02    36   0.132E+00    0.356E+00

DAV:  19    -0.304792993586E+01    0.76822E-02   -0.21612E-03    36   0.306E-01    0.355E+00

DAV:  20    -0.304243348797E+01    0.54964E-02   -0.28878E-03    36   0.420E-01

   1 F= -.30424335E+01 E0= -.29961322E+01  d E =-.926025E-01  mag=     3.9931

 writing wavefunctions

去掉NELM，跑到60步还是没有收敛，高手，还有什么地方不对么？指点指点

赞一下

回复此楼

14楼2010-05-05 16:12:06

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

marklau

银虫 (小有名气)

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

引用回帖:

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