计算化学公社

标题: DFT+U氧空位Ceo2掺Pt吸附丙酮加H结构 SCF不收敛 [打印本页]

作者
Author:
yamin    时间: 2019-11-25 10:21
标题: DFT+U氧空位Ceo2掺Pt吸附丙酮加H结构 SCF不收敛
最终SIGMA = 0.05, ALGO = A ,收敛。供参考。谢谢各位前辈!


问题:DFT+U,SCF不收敛,我不知道怎么办了,卡半个月了。拜托各位前辈帮我看看。文中按顺序贴出了我的INCAR,尝试过的方法以及结构信息。(再补充,不加U时,收敛)
1.  INCAR信息如下

# control
SYSTEM=ceo2-OPT
#SYMPREC=1E-10
LCHARG=.FALSE.
LWAVE=.FALSE.
NPAR=4
LREAL=A
LMAXMIN=6

#SCF
ALGO=F
EDIFF=1E-5
ISMEAR=0
SIGMA=0.1
PREC=NORMAL
NELM=150
NELMIN=6
ENCUT=400

#LVDW=.T.

LDAU=.True.
LDAUL=-1 3  -1 -1 -1
LDAUU=0  5  0  0  0

#OPT
EDIFFG=-0.05
NSW=200
IBRION=2   # 0: MD, 1: Newton, 2: CG, 7:frequency
ISIF=2     # 3 for bulk, 2 for surface and others
#IDIPOL=3   # 1: x, 2:y, 3, z, 4: all

ADDGRID=.TRUE.  # use when IBRION=1
AMIN = 0.01

AMIX     =   0.20
BMIX     =   0.00001
AMIX_MAG =   0.80
BMIX_MAG =   0.00001


2. 按照前辈的方法依次尝试

http://bbs.keinsci.com/thread-11126-1-2.html

(1)调整了很多构型,DFT能优化成功,DFT+U(u=5失败)。排除构型不合理原因。
(2)添加LMAXMIN=6,仍然不收敛
(3)没有读取波函数
(4)设的是ALGO=Normal。
(5)SIGMA=0.5是粗收敛,提高SIGMA=0.1后不收敛
(6)我的ISPIN=1,以下都添加了。
AMIX = 0.2
BMIX = 0.0001 #almost zero, but 0 will crash some versions

AMIX = 0.2
BMIX = 0.0001 #almost zero, but 0 will crash some versions
AMIX MAG = 0.8
BMIX MAG = 0.0001 #almost zero, but 0 will crash some versions
(7)  CeO2是半导体,所以设置ISMEAR=0
(8)  检查体系是不是特殊的磁性排列,即MAGMOM设置是否合理。(关于这点,我的ISPIN=1,还需要考虑吗?如何检查体系的磁性排列?)
(9)  提高积分精度,PREC=Accurate 。(尝试不收敛)
(10) 提高格点精度,ADDGRID = .TRUE. 。(尝试不收敛)
(11) 先用1 1 1 K点计算收敛,再读取CHGCAR,用高K点计算。(现在是K=331,尝试G=111,不收敛)
(12) 尝试不同的ALGO 比如:ALGO=Conjugate。(测试中)
(13) 如果在结构优化或者MD过程中,某一步突然不收敛,使用MAXMIX = 50(不相干,未尝试)
(14) 尝试用更小的ENCUT或者更大的ENCUT的预收敛。(原先是400,调整过450,不收敛)
(15) 换更小的赝势或者更soft的赝势。(赝势中Ce有Ce Ce_3 Ce_GW Ce_h,四个文件,目前选用的是Ce <PAW_PBE Ce 23Dec2003>  。我不知道如何区分软硬,待查询

(16) 最后给出一个VASP官方教程里解决的不收敛的方法:
                       1.用ALGO=N (是否收敛N to 2,Y to 6)
                       2. ICHARG=12 (no charge update,非自洽计算,N to5, Y to 3)
                       3. ICHARG=2 AMIX=0.1 ; BMIX=0.01(N to 4, Y to 6)
                       4. increase BMIXBMIX=3.0 ; AMIN=0.01(N to 5, Y to 6)            
                       5. Bug report
                       6. Use this setting
(这里关于括号中的注释我不理解意思)
3. 结构信息
FM3M 110面  supercell:3*3   3层固定底部一层
(, 下载次数 Times of downloads: 46)


POSCAR

111\(4)                                 
   1.00000000000000     
    10.6602001190000006    0.0000000000000000    0.0000000000000000
     0.0000000000000000   11.3069000244000009    0.0000000000000000
     0.0000000000000000    0.0000000000000000   27.5379009247000006
   O    Ce   Pt   H    C
    36    17     1     8     3
Selective dynamics
Direct
  0.3775976917672311  0.9967906676681029  0.3264518635253797   T   T   T
  0.1211986434165990 -0.0000848686787089  0.3270287977703940   T   T   T
  0.1192760621045083  0.1656408098932839  0.3995894548253406   T   T   T
  0.3716852820342998  0.1525365254326861  0.4020499135214294   T   T   T
  0.3750000000000000  0.1666699980000033  0.2579399959999975   F   F   F
  0.1250000000000000  0.1666699980000033  0.2579399959999975   F   F   F
  0.8774694741001148  0.9992992022559726  0.3271248790005488   T   T   T
  0.6229683574426530  0.9976310514423520  0.3265325703341083   T   T   T
  0.6274086574454923  0.1537290895820219  0.4027789192146946   T   T   T
  0.8795833421614427  0.1662461431575759  0.3987606742540441   T   T   T
  0.8749999779999982  0.1666699980000033  0.2579399959999975   F   F   F
  0.6250000220000018  0.1666699980000033  0.2579399959999975   F   F   F
  0.3742843279251741  0.3401316930774476  0.3287722007531013   T   T   T
  0.1220642187488985  0.3316013640282435  0.3270442505743299   T   T   T
  0.0994840006880943  0.5001926254715579  0.3971099664542040   T   T   T
  0.3785524341893844  0.5162417493073541  0.4057347165646560   T   T   T
  0.3750000000000000  0.5000000000000000  0.2579399959999975   F   F   F
  0.1250000000000000  0.5000000000000000  0.2579399959999975   F   F   F
  0.8796340509623183  0.3356433321028414  0.3275343178164884   T   T   T
  0.6248963732519363  0.3436536937999709  0.3306064802096449   T   T   T
  0.8435817111533865  0.4995391402130468  0.4081039090254628   T   T   T
  0.8749999779999982  0.5000000000000000  0.2579399959999975   F   F   F
  0.6250000220000018  0.5000000000000000  0.2579399959999975   F   F   F
  0.3769137325538578  0.6685716524959775  0.3275622508566277   T   T   T
  0.1217327632120917  0.6685923685000844  0.3269431229156253   T   T   T
  0.1193260043664254  0.8334624044781700  0.3996939172452387   T   T   T
  0.3761815150953545  0.8333900772473088  0.4015075418891398   T   T   T
  0.3750000000000000  0.8333299600000004  0.2579399959999975   F   F   F
  0.1250000000000000  0.8333299600000004  0.2579399959999975   F   F   F
  0.8789188698510234  0.6632847182505904  0.3274594802077431   T   T   T
  0.6220521243423986  0.6585866059742553  0.3312907822427772   T   T   T
  0.6204822334121542  0.8343854027256173  0.3997747452275983   T   T   T
  0.8779765190728333  0.8322824386028105  0.3987521569336377   T   T   T
  0.8749999779999982  0.8333299600000004  0.2579399959999975   F   F   F
  0.6250000220000018  0.8333299600000004  0.2579399959999975   F   F   F
  0.6149774347258308  0.5569219941551554  0.4648991585575832   T   T   T
  0.0001758699549367 -0.0005592783655826  0.3945725910156291   T   T   T
  0.0000000000000000  0.0000000000000000  0.2579399959999975   F   F   F
  0.2494305707209341  0.1651485782496097  0.3269199099009110   T   T   T
  0.4995405547362476  0.9963082183921594  0.3955322205201541   T   T   T
  0.5000000000000000  0.0000000000000000  0.2579399959999975   F   F   F
  0.7505406817268330  0.1648289993144933  0.3273925529985168   T   T   T
  0.0011306858380165  0.3239030770546590  0.3945533385504572   T   T   T
  0.0000000000000000  0.3333300019999967  0.2579399959999975   F   F   F
  0.2467882201629610  0.5011816691804805  0.3265267677040936   T   T   T
  0.5000000000000000  0.3333300019999967  0.2579399959999975   F   F   F
  0.7526582697237455  0.4990302793178982  0.3211712717147434   T   T   T
  0.0003471108900767  0.6753034257090438  0.3947624183296475   T   T   T
  0.0000000000000000  0.6666700399999996  0.2579399959999975   F   F   F
  0.2488348441758422  0.8343305069859456  0.3257020373134448   T   T   T
  0.4946151045993857  0.6833324401168691  0.3962414126667447   T   T   T
  0.5000000000000000  0.6666700399999996  0.2579399959999975   F   F   F
  0.7500706456105846  0.8329205414590178  0.3267695608056639   T   T   T
  0.4994543399795751  0.2515879317709884  0.3691297266503954   T   T   T
  0.7696651224867754  0.5029053448962930  0.4291124766044110   T   T   T
  0.4065865104807400  0.4345279946662692  0.4105320223500483   T   T   T
  0.4778148470147996  0.5511223446689605  0.5619901718087674   T   T   T
  0.5837150147524066  0.6659344414187797  0.5436186222771064   T   T   T
  0.6370334935431059  0.5424160398015719  0.5792818760486264   T   T   T
  0.5188409476122946  0.3354389718683318  0.5181079531725636   T   T   T
  0.6776773804636974  0.3443715112884200  0.5352557416931272   T   T   T
  0.6359505385974850  0.3301423451606869  0.4719758407686335   T   T   T
  0.5741166411296662  0.5706055749721644  0.5498166661754910   T   T   T
  0.6018515266954427  0.5029977260090216  0.5039647639507168   T   T   T
  0.6107685616565711  0.3708952042421303  0.5068137870757166   T   T   T







作者
Author:
yamin    时间: 2019-11-25 10:35
以下是读取SIGMA=0.5时的WAVECAR,用SIGMA=0.3时跑的汇总文件。37-38步de由+00变为+01,我看着能量发散,自洽情况也不乐观,就把任务杀了。
37 F= -.47725522E+03 E0= -.47695581E+03  d E =-.885580E+00
curvature:  -0.00 expect dE=-0.223E-02 dE for cont linesearch -0.867E-04
ZBRENT: interpolating
opt :   0.0001  next Energy=  -473.348392 (dE= 0.302E+01)
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.478793951078E+03   -0.22822E+01   -0.25928E+00  1620   0.289E+00    0.382E+01
RMM:   2    -0.507263048971E+03   -0.28469E+02   -0.11899E+02  1719   0.190E+01    0.592E+01
RMM:   3    -0.497612217431E+03    0.96508E+01   -0.80148E+01  1584   0.219E+01    0.237E+01
RMM:   4    -0.486075957777E+03    0.11536E+02   -0.13463E+01  1432   0.109E+01    0.200E+01
RMM:   5    -0.482756926935E+03    0.33190E+01   -0.55835E+00  1530   0.557E+00    0.161E+01
RMM:   6    -0.479114569654E+03    0.36424E+01   -0.55208E+00  1471   0.526E+00    0.140E+01
RMM:   7    -0.476790388857E+03    0.23242E+01   -0.43832E+00  1582   0.440E+00    0.128E+01
RMM:   8    -0.474032758651E+03    0.27576E+01   -0.50460E+00  1531   0.431E+00    0.903E+00
RMM:   9    -0.473110209977E+03    0.92255E+00   -0.35100E+00  1667   0.348E+00    0.755E+00
RMM:  10    -0.472464291543E+03    0.64592E+00   -0.20089E+00  1670   0.235E+00    0.608E+00
RMM:  11    -0.472250408343E+03    0.21388E+00   -0.77492E-01  1646   0.155E+00    0.528E+00
RMM:  12    -0.471955122787E+03    0.29529E+00   -0.71154E-01  1569   0.171E+00    0.439E+00
RMM:  13    -0.471738839606E+03    0.21628E+00   -0.96531E-01  1550   0.196E+00    0.463E+00
RMM:  14    -0.471722632427E+03    0.16207E-01   -0.53718E-01  1681   0.830E-01    0.337E+00
RMM:  15    -0.471658097444E+03    0.64535E-01   -0.27024E-01  1518   0.919E-01    0.306E+00
RMM:  16    -0.471656904570E+03    0.11929E-02   -0.17959E-01  1716   0.510E-01    0.257E+00
RMM:  17    -0.471671777476E+03   -0.14873E-01   -0.11539E-01  1610   0.482E-01    0.241E+00
RMM:  18    -0.471643882610E+03    0.27895E-01   -0.14341E-01  1530   0.613E-01    0.235E+00
RMM:  19    -0.471632223507E+03    0.11659E-01   -0.25686E-01  1653   0.616E-01    0.177E+00
RMM:  20    -0.471688404117E+03   -0.56181E-01   -0.36893E-01  1686   0.459E-01    0.167E+00
RMM:  21    -0.471768922302E+03   -0.80518E-01   -0.54300E-01  1694   0.247E-01    0.220E+00
RMM:  22    -0.472076190847E+03   -0.30727E+00   -0.12331E+00  1513   0.517E-01    0.542E+00
RMM:  23    -0.472539509599E+03   -0.46332E+00   -0.22261E+00  1616   0.579E-01    0.875E+00
RMM:  24    -0.473112655866E+03   -0.57315E+00   -0.32797E+00  1580   0.799E-01    0.108E+01
RMM:  25    -0.473473365941E+03   -0.36071E+00   -0.34008E+00  1593   0.628E-01    0.123E+01
RMM:  26    -0.473739792179E+03   -0.26643E+00   -0.24943E+00  1573   0.691E-01    0.130E+01
RMM:  27    -0.474054053607E+03   -0.31426E+00   -0.16259E+00  1539   0.792E-01    0.134E+01
RMM:  28    -0.474153596291E+03   -0.99543E-01   -0.76797E-01  1466   0.742E-01    0.134E+01
RMM:  29    -0.474125799452E+03    0.27797E-01   -0.33817E-01  1513   0.660E-01    0.133E+01
RMM:  30    -0.473627207277E+03    0.49859E+00   -0.22692E-01  1332   0.877E-01    0.127E+01
RMM:  31    -0.473350666071E+03    0.27654E+00   -0.19202E-01  1525   0.632E-01    0.123E+01
RMM:  32    -0.473309818506E+03    0.40848E-01   -0.69633E-02  1524   0.349E-01    0.122E+01
RMM:  33    -0.473306607803E+03    0.32107E-02   -0.35654E-02  1438   0.249E-01    0.122E+01
RMM:  34    -0.473624839239E+03   -0.31823E+00   -0.48482E-02  1385   0.434E-01    0.125E+01
RMM:  35    -0.473614985980E+03    0.98533E-02   -0.30829E-02  1419   0.177E-01    0.125E+01
RMM:  36    -0.473702767506E+03   -0.87782E-01   -0.54942E-02  1449   0.213E-01    0.126E+01
RMM:  37    -0.473464411707E+03    0.23836E+00   -0.13017E-01  1301   0.386E-01    0.125E+01
RMM:  38    -0.473235903983E+03    0.22851E+00   -0.18945E-01  1490   0.646E-01    0.123E+01
RMM:  39    -0.474035892744E+03   -0.79999E+00   -0.16174E-01  1304   0.108E+00    0.129E+01
RMM:  40    -0.474532500614E+03   -0.49661E+00   -0.18173E-01  1495   0.749E-01    0.135E+01
RMM:  41    -0.474223577581E+03    0.30892E+00   -0.22521E-02  1294   0.391E-01    0.133E+01
RMM:  42    -0.473725500429E+03    0.49808E+00   -0.77815E-02  1336   0.651E-01    0.129E+01
RMM:  43    -0.473709113713E+03    0.16387E-01   -0.28689E-02  1502   0.195E-01    0.127E+01
RMM:  44    -0.473769764482E+03   -0.60651E-01   -0.65016E-03  1144   0.161E-01    0.127E+01
RMM:  45    -0.473269301876E+03    0.50046E+00   -0.58446E-02  1351   0.576E-01    0.124E+01
RMM:  46    -0.472512716471E+03    0.75659E+00   -0.22173E-01  1480   0.915E-01    0.118E+01
RMM:  47    -0.470479757447E+03    0.20330E+01   -0.14468E+00  1340   0.302E+00    0.636E+00
RMM:  48    -0.472474762100E+03   -0.19950E+01   -0.48715E-01  1376   0.125E+00    0.855E+00
RMM:  49    -0.473382405173E+03   -0.90764E+00   -0.16190E-01  1438   0.120E+00    0.893E+00
RMM:  50    -0.473375768633E+03    0.66365E-02   -0.82209E-02  1719   0.410E-01    0.896E+00
RMM:  51    -0.473353197201E+03    0.22571E-01   -0.13805E-02  1843   0.131E-01    0.891E+00
RMM:  52    -0.473145797969E+03    0.20740E+00   -0.36968E-02  1589   0.293E-01    0.787E+00
RMM:  53    -0.473145392685E+03    0.40528E-03   -0.91109E-03  1476   0.222E-01    0.792E+00
RMM:  54    -0.473162182602E+03   -0.16790E-01   -0.83174E-03  1732   0.115E-01    0.781E+00
RMM:  55    -0.473463988558E+03   -0.30181E+00   -0.86439E-02  1623   0.557E-01    0.828E+00
RMM:  56    -0.473578623992E+03   -0.11464E+00   -0.79517E-02  1646   0.367E-01    0.849E+00
RMM:  57    -0.473828615282E+03   -0.24999E+00   -0.63134E-02  1632   0.442E-01    0.897E+00
RMM:  58    -0.473971852222E+03   -0.14324E+00   -0.30628E-02  1653   0.251E-01    0.925E+00
RMM:  59    -0.473998144114E+03   -0.26292E-01   -0.90931E-03  1675   0.134E-01    0.928E+00
RMM:  60    -0.474063061808E+03   -0.64918E-01   -0.10784E-02  1607   0.134E-01    0.944E+00
RMM:  61    -0.474058471134E+03    0.45907E-02   -0.25173E-03  1466   0.833E-02    0.944E+00
RMM:  62    -0.473588565804E+03    0.46991E+00   -0.10403E-01  1534   0.693E-01    0.889E+00
RMM:  63    -0.473437422851E+03    0.15114E+00   -0.78559E-02  1616   0.445E-01    0.873E+00
RMM:  64    -0.473276594571E+03    0.16083E+00   -0.34660E-02  1543   0.381E-01    0.857E+00
RMM:  65    -0.473459576653E+03   -0.18298E+00   -0.35288E-02  1417   0.544E-01    0.842E+00
RMM:  66    -0.473906703480E+03   -0.44713E+00   -0.20838E-01  1556   0.967E-01    0.850E+00
RMM:  67    -0.468043908764E+03    0.58628E+01    0.47768E+00  1800   0.166E+01    0.203E+01
RMM:  68    -0.468118031671E+03   -0.74123E-01   -0.17777E+01  1580   0.110E+01    0.186E+01
RMM:  69    -0.470176659222E+03   -0.20586E+01   -0.14630E+01  1774   0.361E+00    0.265E+01
RMM:  70    -0.470109088226E+03    0.67571E-01   -0.26807E+01  1540   0.609E+00    0.215E+01
RMM:  71    -0.471662131475E+03   -0.15530E+01   -0.23393E+01  1657   0.331E+00    0.197E+01
RMM:  72    -0.470294317280E+03    0.13678E+01   -0.21169E+01  1347   0.842E+00    0.156E+01
RMM:  73    -0.473017442650E+03   -0.27231E+01   -0.31128E+01  1529   0.353E+00    0.151E+01
RMM:  74    -0.472890467696E+03    0.12697E+00   -0.16122E+00  1401   0.204E+00    0.144E+01
RMM:  75    -0.472991447648E+03   -0.10098E+00   -0.90727E-01  1554   0.724E-01    0.144E+01
RMM:  76    -0.473030508068E+03   -0.39060E-01   -0.20781E-01  1358   0.649E-01    0.144E+01
RMM:  77    -0.472975188892E+03    0.55319E-01   -0.47559E-01  1410   0.118E+00    0.143E+01
RMM:  78    -0.472967318089E+03    0.78708E-02   -0.34940E-01  1498   0.792E-01    0.140E+01
RMM:  79    -0.473007437170E+03   -0.40119E-01   -0.28850E-01  1468   0.676E-01    0.140E+01
RMM:  80    -0.473115609223E+03   -0.10817E+00   -0.25148E-01  1453   0.700E-01    0.140E+01
RMM:  81    -0.473097453366E+03    0.18156E-01   -0.21414E-01  1436   0.667E-01    0.138E+01
RMM:  82    -0.473061768134E+03    0.35685E-01   -0.37164E-01  1477   0.102E+00    0.138E+01
RMM:  83    -0.472974086856E+03    0.87681E-01   -0.22528E-01  1461   0.826E-01    0.135E+01
RMM:  84    -0.473051093501E+03   -0.77007E-01   -0.14788E-01  1520   0.527E-01    0.135E+01
RMM:  85    -0.473002321310E+03    0.48772E-01   -0.11180E-01  1443   0.573E-01    0.133E+01
RMM:  86    -0.472786820536E+03    0.21550E+00   -0.21350E-01  1515   0.756E-01    0.131E+01
RMM:  87    -0.472768049834E+03    0.18771E-01   -0.24683E-01  1384   0.100E+00    0.130E+01
RMM:  88    -0.472712580166E+03    0.55470E-01   -0.19001E-01  1537   0.593E-01    0.129E+01
RMM:  89    -0.472715014964E+03   -0.24348E-02   -0.74407E-02  1490   0.400E-01    0.129E+01
RMM:  90    -0.472727186435E+03   -0.12171E-01   -0.69669E-02  1473   0.316E-01    0.128E+01
RMM:  91    -0.472755650110E+03   -0.28464E-01   -0.55389E-02  1393   0.338E-01    0.128E+01
RMM:  92    -0.473335144101E+03   -0.57949E+00   -0.16214E+00  1395   0.253E+00    0.156E+01
RMM:  93    -0.476138101488E+03   -0.28030E+01   -0.48250E+00  1546   0.203E+00    0.326E+01
RMM:  94    -0.478364475280E+03   -0.22264E+01   -0.47308E+00  1387   0.256E+00    0.423E+01
RMM:  95    -0.478145896215E+03    0.21858E+00   -0.14477E+00  1396   0.193E+00    0.469E+01
RMM:  96    -0.478710247528E+03   -0.56435E+00   -0.51949E-01  1403   0.137E+00    0.493E+01
RMM:  97    -0.480017081958E+03   -0.13068E+01   -0.45196E-01  1477   0.121E+00    0.498E+01
RMM:  98    -0.479666796419E+03    0.35029E+00   -0.17149E-01  1475   0.794E-01    0.485E+01
RMM:  99    -0.477807770978E+03    0.18590E+01   -0.69721E-01  1370   0.183E+00    0.456E+01
RMM: 100    -0.478381297911E+03   -0.57353E+00   -0.39901E-01  1432   0.119E+00    0.454E+01
RMM: 101    -0.478487919647E+03   -0.10662E+00   -0.98770E-02  1538   0.585E-01    0.453E+01
RMM: 102    -0.467626232836E+03    0.10862E+02   -0.23642E+00  1403   0.332E+00    0.257E+01
RMM: 103    -0.468761008131E+03   -0.11348E+01   -0.14614E+00  1547   0.197E+00    0.573E+00
RMM: 104    -0.471457423884E+03   -0.26964E+01   -0.84932E-01  1457   0.239E+00    0.305E+00
RMM: 105    -0.472427790552E+03   -0.97037E+00   -0.18975E+00  1647   0.252E+00    0.123E+01
RMM: 106    -0.473993165283E+03   -0.15654E+01   -0.82174E-01  1371   0.165E+00    0.131E+01
RMM: 107    -0.473687180538E+03    0.30598E+00   -0.45205E-01  1501   0.130E+00    0.127E+01
RMM: 108    -0.473546764887E+03    0.14042E+00   -0.11319E-01  1551   0.477E-01    0.126E+01
RMM: 109    -0.473485623064E+03    0.61142E-01   -0.60987E-02  1476   0.461E-01    0.126E+01
RMM: 110    -0.473559206775E+03   -0.73584E-01   -0.66562E-02  1542   0.383E-01    0.127E+01
RMM: 111    -0.473601130744E+03   -0.41924E-01   -0.19535E-02  1557   0.186E-01    0.128E+01
RMM: 112    -0.471794422468E+03    0.18067E+01   -0.71637E-01  1406   0.193E+00    0.113E+01
RMM: 113    -0.470868907930E+03    0.92551E+00   -0.83567E-01  1535   0.140E+00    0.104E+01
RMM: 114    -0.470929557123E+03   -0.60649E-01   -0.10368E-01  1418   0.455E-01    0.103E+01
RMM: 115    -0.471113382613E+03   -0.18383E+00   -0.15623E-02  1458   0.242E-01    0.105E+01
RMM: 116    -0.471248699938E+03   -0.13532E+00   -0.14656E-02  1506   0.207E-01    0.106E+01
RMM: 117    -0.472146088975E+03   -0.89739E+00   -0.24649E-01  1420   0.113E+00    0.109E+01
RMM: 118    -0.472228799173E+03   -0.82710E-01   -0.91457E-02  1525   0.319E-01    0.111E+01
RMM: 119    -0.472140860568E+03    0.87939E-01   -0.39536E-02  1433   0.416E-01    0.111E+01
RMM: 120    -0.472096622125E+03    0.44238E-01   -0.36782E-02  1529   0.273E-01    0.111E+01
RMM: 121    -0.472149327345E+03   -0.52705E-01   -0.43699E-03  1315   0.185E-01    0.111E+01
RMM: 122    -0.472868887473E+03   -0.71956E+00   -0.20181E-01  1387   0.105E+00    0.116E+01
RMM: 123    -0.473463525443E+03   -0.59464E+00   -0.32326E-01  1535   0.931E-01    0.122E+01
RMM: 124    -0.472727017893E+03    0.73651E+00   -0.26683E-01  1397   0.120E+00    0.120E+01
RMM: 125    -0.472402997591E+03    0.32402E+00   -0.28440E-01  1536   0.756E-01    0.116E+01
RMM: 126    -0.472267089612E+03    0.13591E+00   -0.36130E-02  1518   0.310E-01    0.115E+01
RMM: 127    -0.472549347340E+03   -0.28226E+00   -0.25294E-02  1389   0.366E-01    0.117E+01
RMM: 128    -0.472566037139E+03   -0.16690E-01   -0.72947E-03  1439   0.104E-01    0.118E+01
RMM: 129    -0.472556022885E+03    0.10014E-01   -0.10724E-03  1217   0.573E-02    0.118E+01
RMM: 130    -0.472137183182E+03    0.41884E+00   -0.60618E-02  1418   0.538E-01    0.116E+01
RMM: 131    -0.472137659927E+03   -0.47675E-03   -0.16955E-02  1436   0.142E-01    0.115E+01
RMM: 132    -0.472769265852E+03   -0.63161E+00   -0.12285E-01  1419   0.786E-01    0.117E+01
RMM: 133    -0.472916250353E+03   -0.14698E+00   -0.76585E-02  1527   0.418E-01    0.120E+01
RMM: 134    -0.472928041246E+03   -0.11791E-01   -0.67483E-03  1544   0.131E-01    0.121E+01
RMM: 135    -0.472676186914E+03    0.25185E+00   -0.12750E-01  1504   0.660E-01    0.120E+01
RMM: 136    -0.472697851957E+03   -0.21665E-01   -0.22912E-02  1516   0.195E-01    0.120E+01
RMM: 137    -0.473894053131E+03   -0.11962E+01   -0.73608E-01  1437   0.181E+00    0.124E+01
RMM: 138    -0.473166922504E+03    0.72713E+00   -0.15321E-01  1398   0.927E-01    0.124E+01
RMM: 139    -0.472813058068E+03    0.35386E+00   -0.15834E-01  1526   0.641E-01    0.123E+01
RMM: 140    -0.472269523933E+03    0.54353E+00   -0.26784E-01  1508   0.931E-01    0.119E+01
RMM: 141    -0.471842079455E+03    0.42744E+00   -0.11585E-01  1535   0.565E-01    0.116E+01
RMM: 142    -0.471132394315E+03    0.70969E+00   -0.44276E-01  1525   0.127E+00    0.111E+01
RMM: 143    -0.470144609801E+03    0.98778E+00   -0.62467E-01  1476   0.157E+00    0.789E+00
RMM: 144    -0.470757289583E+03   -0.61268E+00   -0.27119E-01  1502   0.913E-01    0.461E+00
RMM: 145    -0.472174189585E+03   -0.14169E+01   -0.39117E-01  1582   0.159E+00    0.496E+00
RMM: 146    -0.471977057091E+03    0.19713E+00   -0.94269E-01  1702   0.174E+00    0.374E+00
RMM: 147    -0.471887099507E+03    0.89958E-01   -0.32495E-01  1697   0.965E-01    0.316E+00
RMM: 148    -0.471768767985E+03    0.11833E+00   -0.19803E-01  1626   0.784E-01    0.234E+00
RMM: 149    -0.471786270147E+03   -0.17502E-01   -0.73868E-02  1688   0.385E-01    0.256E+00
RMM: 150    -0.471816836392E+03   -0.30566E-01   -0.35021E-02  1565   0.358E-01    0.283E+00
RMM: 151    -0.471834389204E+03   -0.17553E-01   -0.26905E-02  1684   0.176E-01    0.377E+00
RMM: 152    -0.471840980690E+03   -0.65915E-02   -0.11321E-01  1504   0.625E-01    0.650E+00
RMM: 153    -0.472188174007E+03   -0.34719E+00   -0.21682E-01  1512   0.522E-01    0.111E+01
RMM: 154    -0.472793129488E+03   -0.60496E+00   -0.52502E-01  1307   0.886E-01    0.137E+01
RMM: 155    -0.472985159801E+03   -0.19203E+00   -0.34462E-01  1421   0.859E-01    0.139E+01
RMM: 156    -0.473032180244E+03   -0.47020E-01   -0.26226E-01  1522   0.632E-01    0.140E+01
RMM: 157    -0.473014112037E+03    0.18068E-01   -0.25001E-01  1421   0.654E-01    0.139E+01
RMM: 158    -0.473066678942E+03   -0.52567E-01   -0.24964E-01  1501   0.625E-01    0.139E+01
RMM: 159    -0.473128550310E+03   -0.61871E-01   -0.18455E-01  1367   0.588E-01    0.139E+01
RMM: 160    -0.473152522151E+03   -0.23972E-01   -0.20497E-01  1443   0.495E-01    0.139E+01
RMM: 161    -0.473191857392E+03   -0.39335E-01   -0.11591E-01  1280   0.515E-01    0.138E+01
RMM: 162    -0.473228388931E+03   -0.36532E-01   -0.14126E-01  1495   0.449E-01    0.138E+01
RMM: 163    -0.473255984859E+03   -0.27596E-01   -0.82176E-02  1256   0.444E-01    0.137E+01
RMM: 164    -0.473274057770E+03   -0.18073E-01   -0.10109E-01  1484   0.404E-01    0.136E+01
RMM: 165    -0.473306770613E+03   -0.32713E-01   -0.80567E-02  1379   0.532E-01    0.136E+01
RMM: 166    -0.473335576125E+03   -0.28806E-01   -0.10435E-01  1506   0.457E-01    0.135E+01
RMM: 167    -0.473304704936E+03    0.30871E-01   -0.79282E-02  1395   0.550E-01    0.134E+01
RMM: 168    -0.473310356004E+03   -0.56511E-02   -0.93147E-02  1520   0.385E-01    0.134E+01
RMM: 169    -0.473333928468E+03   -0.23572E-01   -0.59510E-02  1373   0.440E-01    0.133E+01
RMM: 170    -0.473273216351E+03    0.60712E-01   -0.10643E-01  1504   0.546E-01    0.133E+01
RMM: 171    -0.473260372785E+03    0.12844E-01   -0.58251E-02  1434   0.407E-01    0.132E+01
RMM: 172    -0.473258016223E+03    0.23566E-02   -0.91305E-02  1500   0.473E-01    0.130E+01
RMM: 173    -0.473250575081E+03    0.74411E-02   -0.52764E-02  1388   0.392E-01    0.130E+01
RMM: 174    -0.473233390606E+03    0.17184E-01   -0.70116E-02  1464   0.350E-01    0.129E+01
RMM: 175    -0.473287861429E+03   -0.54471E-01   -0.48319E-02  1332   0.410E-01    0.129E+01
RMM: 176    -0.473370423437E+03   -0.82562E-01   -0.70858E-02  1504   0.371E-01    0.129E+01
RMM: 177    -0.473372485650E+03   -0.20622E-02   -0.56634E-02  1425   0.447E-01    0.128E+01
RMM: 178    -0.473368287930E+03    0.41977E-02   -0.56073E-02  1383   0.280E-01    0.128E+01
RMM: 179    -0.473382180798E+03   -0.13893E-01   -0.36955E-02  1336   0.323E-01    0.127E+01
RMM: 180    -0.473374893275E+03    0.72875E-02   -0.45951E-02  1378   0.256E-01    0.127E+01
RMM: 181    -0.473498985426E+03   -0.12409E+00   -0.54712E-02  1433   0.463E-01    0.126E+01
RMM: 182    -0.473637455919E+03   -0.13847E+00   -0.79883E-02  1423   0.530E-01    0.127E+01
RMM: 183    -0.473669597762E+03   -0.32142E-01   -0.40614E-02  1458   0.314E-01    0.127E+01
RMM: 184    -0.473675588756E+03   -0.59910E-02   -0.33895E-02  1317   0.236E-01    0.127E+01
RMM: 185    -0.473608524768E+03    0.67064E-01   -0.25951E-02  1365   0.288E-01    0.127E+01
RMM: 186    -0.473556211404E+03    0.52313E-01   -0.29876E-02  1410   0.235E-01    0.126E+01
RMM: 187    -0.472707182570E+03    0.84903E+00   -0.13323E+00  1329   0.279E+00    0.123E+01
RMM: 188    -0.472636271877E+03    0.70911E-01   -0.62112E-01  1543   0.809E-01    0.122E+01
RMM: 189    -0.472693481755E+03   -0.57210E-01   -0.49571E-02  1446   0.429E-01    0.121E+01
RMM: 190    -0.472718235790E+03   -0.24754E-01   -0.13300E-01  1524   0.620E-01    0.121E+01
RMM: 191    -0.472740041307E+03   -0.21806E-01   -0.58830E-02  1414   0.544E-01    0.120E+01
RMM: 192    -0.473610322706E+03   -0.87028E+00   -0.42386E-01  1432   0.145E+00    0.125E+01
RMM: 193    -0.473516646410E+03    0.93676E-01   -0.13539E-01  1550   0.408E-01    0.127E+01
RMM: 194    -0.473450438326E+03    0.66208E-01   -0.21455E-02  1460   0.254E-01    0.126E+01
RMM: 195    -0.473312702409E+03    0.13774E+00   -0.63571E-02  1529   0.427E-01    0.125E+01
RMM: 196    -0.472960773611E+03    0.35193E+00   -0.10958E-01  1460   0.662E-01    0.124E+01
RMM: 197    -0.479522489694E+03   -0.65617E+01   -0.35100E+01  1535   0.106E+01    0.415E+01
RMM: 198    -0.474900371155E+03    0.46221E+01   -0.54812E+00  1295   0.644E+00    0.276E+01
RMM: 199    -0.472948424023E+03    0.19519E+01   -0.44641E+00  1528   0.400E+00    0.190E+01
RMM: 200    -0.471915897273E+03    0.10325E+01   -0.82965E+00  1551   0.584E+00    0.102E+01
RMM: 201    -0.472287064043E+03   -0.37117E+00   -0.37061E+00  1506   0.430E+00    0.642E+00
RMM: 202    -0.472050513353E+03    0.23655E+00   -0.21180E+00  1635   0.220E+00    0.472E+00
RMM: 203    -0.471889912675E+03    0.16060E+00   -0.64991E-01  1662   0.125E+00    0.387E+00
RMM: 204    -0.471807125421E+03    0.82787E-01   -0.16661E-01  1630   0.816E-01    0.323E+00
RMM: 205    -0.471797020394E+03    0.10105E-01   -0.74487E-02  1691   0.450E-01    0.293E+00
RMM: 206    -0.471770270106E+03    0.26750E-01   -0.52578E-02  1636   0.394E-01    0.259E+00
RMM: 207    -0.471768033692E+03    0.22364E-02   -0.35974E-02  1686   0.299E-01    0.240E+00
RMM: 208    -0.471759291356E+03    0.87423E-02   -0.32931E-02  1715   0.244E-01    0.227E+00
RMM: 209    -0.471733953055E+03    0.25338E-01   -0.37730E-02  1556   0.336E-01    0.200E+00
RMM: 210    -0.471735204030E+03   -0.12510E-02   -0.36887E-02  1697   0.238E-01    0.187E+00
RMM: 211    -0.471742615618E+03   -0.74116E-02   -0.15158E-02  1560   0.148E-01    0.205E+00
RMM: 212    -0.471780420348E+03   -0.37805E-01   -0.31787E-02  1646   0.310E-01    0.319E+00
RMM: 213    -0.471887218898E+03   -0.10680E+00   -0.44689E-02  1573   0.236E-01    0.699E+00
RMM: 214    -0.472599171865E+03   -0.71195E+00   -0.31115E-01  1404   0.654E-01    0.130E+01
RMM: 215    -0.473080226542E+03   -0.48105E+00   -0.58840E-01  1298   0.113E+00    0.137E+01
RMM: 216    -0.473105980478E+03   -0.25754E-01   -0.23059E-01  1471   0.682E-01    0.139E+01
RMM: 217    -0.473133067374E+03   -0.27087E-01   -0.15004E-01  1503   0.542E-01    0.141E+01
RMM: 218    -0.473186926442E+03   -0.53859E-01   -0.95639E-02  1477   0.411E-01    0.142E+01
RMM: 219    -0.473209478650E+03   -0.22552E-01   -0.62514E-02  1364   0.338E-01    0.143E+01
RMM: 220    -0.473215823504E+03   -0.63449E-02   -0.50109E-02  1444   0.280E-01    0.144E+01
RMM: 221    -0.473216006902E+03   -0.18340E-03   -0.25342E-02  1363   0.229E-01    0.145E+01
RMM: 222    -0.473162143178E+03    0.53864E-01   -0.40272E-02  1505   0.348E-01    0.144E+01
RMM: 223    -0.473143963027E+03    0.18180E-01   -0.15057E-02  1464   0.206E-01    0.144E+01
RMM: 224    -0.473130230569E+03    0.13732E-01   -0.32883E-02  1354   0.378E-01    0.144E+01
RMM: 225    -0.473056725234E+03    0.73505E-01   -0.33243E-02  1527   0.302E-01    0.144E+01
RMM: 226    -0.472982221489E+03    0.74504E-01   -0.19354E-02  1501   0.198E-01    0.144E+01
RMM: 227    -0.472923967097E+03    0.58254E-01   -0.15635E-02  1462   0.233E-01    0.144E+01
RMM: 228    -0.472759011708E+03    0.16496E+00   -0.38552E-02  1475   0.385E-01    0.144E+01
RMM: 229    -0.472673312861E+03    0.85699E-01   -0.24928E-02  1533   0.231E-01    0.143E+01
RMM: 230    -0.472608065226E+03    0.65248E-01   -0.13475E-02  1451   0.212E-01    0.143E+01
RMM: 231    -0.472537964031E+03    0.70101E-01   -0.23696E-02  1542   0.218E-01    0.143E+01
RMM: 232    -0.472553373873E+03   -0.15410E-01   -0.19879E-02  1396   0.299E-01    0.143E+01
RMM: 233    -0.472526116543E+03    0.27257E-01   -0.15123E-02  1502   0.165E-01    0.142E+01
RMM: 234    -0.472619119063E+03   -0.93003E-01   -0.15406E-02  1399   0.238E-01    0.142E+01
RMM: 235    -0.472772421858E+03   -0.15330E+00   -0.34292E-02  1510   0.274E-01    0.143E+01
RMM: 236    -0.472801802980E+03   -0.29381E-01   -0.16192E-02  1373   0.190E-01    0.143E+01
RMM: 237    -0.472704044086E+03    0.97759E-01   -0.39379E-02  1500   0.364E-01    0.142E+01
RMM: 238    -0.472711642377E+03   -0.75983E-02   -0.21661E-02  1386   0.180E-01    0.142E+01
RMM: 239    -0.472691276755E+03    0.20366E-01   -0.19851E-02  1410   0.221E-01    0.141E+01
RMM: 240    -0.472714869585E+03   -0.23593E-01   -0.31734E-02  1480   0.239E-01    0.141E+01
RMM: 241    -0.472757461153E+03   -0.42592E-01   -0.38356E-02  1489   0.335E-01    0.141E+01
RMM: 242    -0.472328987112E+03    0.42847E+00   -0.14933E-01  1395   0.876E-01    0.141E+01
RMM: 243    -0.472300271172E+03    0.28716E-01   -0.89386E-02  1542   0.353E-01    0.141E+01
RMM: 244    -0.472379166930E+03   -0.78896E-01   -0.34400E-02  1436   0.331E-01    0.140E+01
RMM: 245    -0.472322405604E+03    0.56761E-01   -0.48565E-02  1527   0.264E-01    0.140E+01
RMM: 246    -0.472309593532E+03    0.12812E-01   -0.29572E-02  1330   0.271E-01    0.140E+01
RMM: 247    -0.472335920299E+03   -0.26327E-01   -0.57679E-02  1527   0.332E-01    0.138E+01
RMM: 248    -0.472350286273E+03   -0.14366E-01   -0.32438E-02  1388   0.275E-01    0.137E+01
RMM: 249    -0.472293172381E+03    0.57114E-01   -0.49319E-02  1488   0.242E-01    0.137E+01
RMM: 250    -0.472305225913E+03   -0.12054E-01   -0.34435E-02  1169   0.278E-01    0.136E+01
RMM: 251    -0.472328183967E+03   -0.22958E-01   -0.51105E-02  1385   0.221E-01    0.135E+01
RMM: 252    -0.472190505279E+03    0.13768E+00   -0.52534E-02  1394   0.436E-01    0.134E+01
RMM: 253    -0.472189005663E+03    0.14996E-02   -0.61541E-02  1452   0.240E-01    0.133E+01
RMM: 254    -0.472173395925E+03    0.15610E-01   -0.43687E-02  1363   0.341E-01    0.133E+01
RMM: 255    -0.472064923399E+03    0.10847E+00   -0.73482E-02  1519   0.341E-01    0.132E+01
RMM: 256    -0.471981569051E+03    0.83354E-01   -0.51441E-02  1379   0.373E-01    0.131E+01
RMM: 257    -0.483068028124E+03   -0.11086E+02   -0.42218E+01  1514   0.125E+01    0.304E+01
RMM: 258    -0.477345937718E+03    0.57221E+01   -0.68867E+00  1292   0.810E+00    0.171E+01
RMM: 259    -0.474078013158E+03    0.32679E+01   -0.51638E+00  1406   0.543E+00    0.105E+01
RMM: 260    -0.473669421606E+03    0.40859E+00   -0.19525E+00  1523   0.262E+00    0.861E+00
RMM: 261    -0.473211535640E+03    0.45789E+00   -0.48857E-01  1660   0.135E+00    0.803E+00
RMM: 262    -0.472038748236E+03    0.11728E+01   -0.21708E+00  1613   0.305E+00    0.434E+00
RMM: 263    -0.471946783456E+03    0.91965E-01   -0.61011E-01  1679   0.106E+00    0.422E+00
RMM: 264    -0.471860594132E+03    0.86189E-01   -0.15484E-01  1698   0.607E-01    0.414E+00
RMM: 265    -0.471944085319E+03   -0.83491E-01   -0.59887E-02  1687   0.283E-01    0.705E+00
RMM: 266    -0.472533786296E+03   -0.58970E+00   -0.29050E-01  1502   0.763E-01    0.130E+01
RMM: 267    -0.473172995874E+03   -0.63921E+00   -0.69816E-01  1396   0.128E+00    0.140E+01
RMM: 268    -0.473241672930E+03   -0.68677E-01   -0.22863E-01  1486   0.656E-01    0.143E+01
RMM: 269    -0.473209260949E+03    0.32412E-01   -0.11258E-01  1496   0.488E-01    0.144E+01
RMM: 270    -0.473242134987E+03   -0.32874E-01   -0.74633E-02  1504   0.373E-01    0.145E+01
RMM: 271    -0.473353244440E+03   -0.11111E+00   -0.59515E-02  1430   0.411E-01    0.147E+01
RMM: 272    -0.473341520387E+03    0.11724E-01   -0.10562E-01  1460   0.614E-01    0.147E+01
RMM: 273    -0.472994480299E+03    0.34704E+00   -0.17163E-01  1478   0.757E-01    0.146E+01
RMM: 274    -0.472885004225E+03    0.10948E+00   -0.77916E-02  1532   0.340E-01    0.145E+01
RMM: 275    -0.472877461243E+03    0.75430E-02   -0.20159E-02  1451   0.229E-01    0.145E+01
RMM: 276    -0.472853492580E+03    0.23969E-01   -0.26032E-02  1518   0.185E-01    0.146E+01
RMM: 277    -0.473019433531E+03   -0.16594E+00   -0.42029E-02  1469   0.405E-01    0.147E+01
RMM: 278    -0.472941115093E+03    0.78318E-01   -0.21398E-02  1344   0.268E-01    0.147E+01
RMM: 279    -0.472832180532E+03    0.10893E+00   -0.24175E-02  1484   0.278E-01    0.146E+01
RMM: 280    -0.472533223514E+03    0.29896E+00   -0.49779E-02  1428   0.425E-01    0.144E+01
RMM: 281    -0.472270919778E+03    0.26230E+00   -0.60426E-02  1506   0.397E-01    0.142E+01
RMM: 282    -0.472251153246E+03    0.19767E-01   -0.28382E-02  1525   0.232E-01    0.141E+01
RMM: 283    -0.472181246568E+03    0.69907E-01   -0.28519E-02  1492   0.298E-01    0.141E+01
RMM: 284    -0.472175030144E+03    0.62164E-02   -0.23983E-02  1383   0.162E-01    0.140E+01
RMM: 285    -0.472143053319E+03    0.31977E-01   -0.18792E-02  1405   0.236E-01    0.140E+01
RMM: 286    -0.472139463543E+03    0.35898E-02   -0.25072E-02  1338   0.162E-01    0.140E+01
RMM: 287    -0.472244839080E+03   -0.10538E+00   -0.44075E-02  1498   0.401E-01    0.139E+01
RMM: 288    -0.472262862629E+03   -0.18024E-01   -0.31265E-02  1417   0.185E-01    0.139E+01
RMM: 289    -0.472269398755E+03   -0.65361E-02   -0.21618E-02  1217   0.237E-01    0.138E+01
RMM: 290    -0.472186274618E+03    0.83124E-01   -0.35755E-02  1458   0.245E-01    0.138E+01
RMM: 291    -0.472046334897E+03    0.13994E+00   -0.36954E-02  1393   0.367E-01    0.137E+01
RMM: 292    -0.472034120687E+03    0.12214E-01   -0.48030E-02  1481   0.301E-01    0.136E+01
RMM: 293    -0.472141499098E+03   -0.10738E+00   -0.30173E-02  1313   0.290E-01    0.136E+01
RMM: 294    -0.472084465248E+03    0.57034E-01   -0.41858E-02  1395   0.234E-01    0.135E+01
RMM: 295    -0.472132528294E+03   -0.48063E-01   -0.30814E-02  1244   0.275E-01    0.134E+01
RMM: 296    -0.472064390037E+03    0.68138E-01   -0.46736E-02  1454   0.252E-01    0.134E+01
RMM: 297    -0.472075551720E+03   -0.11162E-01   -0.34011E-02  1346   0.305E-01    0.133E+01
RMM: 298    -0.472090755355E+03   -0.15204E-01   -0.49241E-02  1479   0.239E-01    0.132E+01
RMM: 299    -0.472015972928E+03    0.74782E-01   -0.36607E-02  1344   0.321E-01    0.131E+01
RMM: 300    -0.471978450539E+03    0.37522E-01   -0.50913E-02  1482   0.246E-01    0.130E+01
RMM: 301    -0.471904741216E+03    0.73709E-01   -0.38319E-02  1339   0.328E-01    0.129E+01
RMM: 302    -0.470851123317E+03    0.10536E+01   -0.37611E-01  1385   0.136E+00    0.111E+01
RMM: 303    -0.470862611980E+03   -0.11489E-01   -0.16133E-01  1444   0.448E-01    0.989E+00
RMM: 304    -0.471068771312E+03   -0.20616E+00   -0.67162E-02  1376   0.586E-01    0.100E+01
RMM: 305    -0.471343987673E+03   -0.27522E+00   -0.90098E-02  1537   0.454E-01    0.113E+01
RMM: 306    -0.471316769199E+03    0.27218E-01   -0.62459E-02  1325   0.428E-01    0.118E+01
RMM: 307    -0.470798012435E+03    0.51876E+00   -0.28979E+00  1428   0.389E+00    0.892E+00
RMM: 308    -0.472368494202E+03   -0.15705E+01   -0.19073E+00  1370   0.439E+00    0.544E+00
RMM: 309    -0.472227182065E+03    0.14131E+00   -0.62207E-01  1434   0.173E+00    0.514E+00
RMM: 310    -0.472211177435E+03    0.16005E-01   -0.17297E-01  1683   0.581E-01    0.592E+00
RMM: 311    -0.472486054759E+03   -0.27488E+00   -0.20932E-01  1507   0.388E-01    0.103E+01
RMM: 312    -0.472700813238E+03   -0.21476E+00   -0.13521E+00  1498   0.186E+00    0.136E+01
RMM: 313    -0.472791884330E+03   -0.91071E-01   -0.74059E-01  1505   0.919E-01    0.137E+01
RMM: 314    -0.472825950172E+03   -0.34066E-01   -0.19049E-01  1527   0.420E-01    0.137E+01
RMM: 315    -0.472974309143E+03   -0.14836E+00   -0.19791E-01  1442   0.460E-01    0.136E+01
RMM: 316    -0.473090512443E+03   -0.11620E+00   -0.22194E-01  1355   0.593E-01    0.133E+01
RMM: 317    -0.472726554817E+03    0.36396E+00   -0.67736E-01  1457   0.152E+00    0.126E+01
RMM: 318    -0.472712271452E+03    0.14283E-01   -0.36679E-01  1462   0.104E+00    0.122E+01
RMM: 319    -0.472733234688E+03   -0.20963E-01   -0.19634E-01  1429   0.682E-01    0.121E+01
RMM: 320    -0.472828546524E+03   -0.95312E-01   -0.12688E-01  1498   0.546E-01    0.121E+01
RMM: 321    -0.472863454576E+03   -0.34908E-01   -0.54159E-02  1525   0.367E-01    0.121E+01
RMM: 322    -0.472995928008E+03   -0.13247E+00   -0.48507E-02  1513   0.357E-01    0.123E+01
RMM: 323    -0.472997047052E+03   -0.11190E-02   -0.11356E-02  1395   0.162E-01    0.123E+01
RMM: 324    -0.472849830797E+03    0.14722E+00   -0.18353E-02  1467   0.231E-01    0.122E+01
RMM: 325    -0.472909457275E+03   -0.59626E-01   -0.73553E-03  1293   0.135E-01    0.122E+01
RMM: 326    -0.472871001193E+03    0.38456E-01   -0.71717E-03  1402   0.123E-01    0.122E+01
RMM: 327    -0.472829687390E+03    0.41314E-01   -0.21398E-02  1445   0.299E-01    0.121E+01
RMM: 328    -0.472986033629E+03   -0.15635E+00   -0.44362E-02  1413   0.445E-01    0.122E+01
RMM: 329    -0.472977958698E+03    0.80749E-02   -0.35852E-02  1512   0.273E-01    0.122E+01
RMM: 330    -0.472973087727E+03    0.48710E-02   -0.71461E-03  1288   0.108E-01    0.122E+01
RMM: 331    -0.472920170771E+03    0.52917E-01   -0.55246E-03  1396   0.138E-01    0.121E+01
RMM: 332    -0.473123845966E+03   -0.20368E+00   -0.53602E-02  1434   0.491E-01    0.122E+01
RMM: 333    -0.473050167216E+03    0.73679E-01   -0.20549E-02  1517   0.207E-01    0.122E+01
RMM: 334    -0.473027611759E+03    0.22555E-01   -0.86309E-03  1513   0.155E-01    0.122E+01
RMM: 335    -0.472900291908E+03    0.12732E+00   -0.16429E-02  1425   0.238E-01    0.121E+01
RMM: 336    -0.472788662709E+03    0.11163E+00   -0.11986E-02  1468   0.187E-01    0.120E+01
RMM: 337    -0.472047503678E+03    0.74116E+00   -0.25677E-01  1387   0.119E+00    0.117E+01
RMM: 338    -0.471886642115E+03    0.16086E+00   -0.17074E-01  1543   0.496E-01    0.114E+01
RMM: 339    -0.471828578879E+03    0.58063E-01   -0.24378E-02  1420   0.304E-01    0.113E+01
RMM: 340    -0.471790964801E+03    0.37614E-01   -0.16849E-02  1558   0.146E-01    0.113E+01
RMM: 341    -0.471779922576E+03    0.11042E-01   -0.89595E-03  1448   0.106E-01    0.112E+01
RMM: 342    -0.472132724294E+03   -0.35280E+00   -0.72099E-02  1477   0.545E-01    0.114E+01
RMM: 343    -0.472133053205E+03   -0.32891E-03   -0.18381E-02  1398   0.156E-01    0.115E+01
RMM: 344    -0.472018182373E+03    0.11487E+00   -0.76007E-03  1366   0.158E-01    0.114E+01
RMM: 345    -0.472025876117E+03   -0.76937E-02   -0.64699E-03  1100   0.747E-02    0.114E+01
RMM: 346    -0.472046806397E+03   -0.20930E-01   -0.41760E-03  1042   0.929E-02    0.114E+01
RMM: 347    -0.471585209774E+03    0.46160E+00   -0.93584E-02  1332   0.707E-01    0.110E+01
RMM: 348    -0.471480549329E+03    0.10466E+00   -0.49036E-02  1521   0.315E-01    0.109E+01
RMM: 349    -0.471428136153E+03    0.52413E-01   -0.16835E-02  1530   0.219E-01    0.108E+01
RMM: 350    -0.471464822985E+03   -0.36687E-01   -0.17548E-02  1528   0.214E-01    0.109E+01
RMM: 351    -0.471514912249E+03   -0.50089E-01   -0.77040E-03  1492   0.122E-01    0.109E+01
RMM: 352    -0.471779171662E+03   -0.26426E+00   -0.39985E+00  1411   0.509E+00    0.119E+01
RMM: 353    -0.473801085113E+03   -0.20219E+01   -0.99946E-01  1407   0.311E+00    0.833E+00
RMM: 354    -0.472747294836E+03    0.10538E+01   -0.82128E-01  1460   0.232E+00    0.592E+00
RMM: 355    -0.472297809791E+03    0.44949E+00   -0.64513E-01  1677   0.146E+00    0.458E+00
RMM: 356    -0.472159535773E+03    0.13827E+00   -0.39784E-01  1674   0.109E+00    0.469E+00
RMM: 357    -0.471877626683E+03    0.28191E+00   -0.21899E+00  1527   0.311E+00    0.598E+00
RMM: 358    -0.471904911336E+03   -0.27285E-01   -0.12208E+00  1375   0.212E+00    0.623E+00
RMM: 359    -0.472060389569E+03   -0.15548E+00   -0.10793E+00  1659   0.822E-01    0.105E+01
RMM: 360    -0.472609390322E+03   -0.54900E+00   -0.14324E+00  1496   0.968E-01    0.131E+01
RMM: 361    -0.472731772920E+03   -0.12238E+00   -0.42822E-01  1466   0.643E-01    0.133E+01
RMM: 362    -0.472743426907E+03   -0.11654E-01   -0.20608E-01  1517   0.550E-01    0.133E+01
RMM: 363    -0.472781467623E+03   -0.38041E-01   -0.19482E-01  1546   0.434E-01    0.133E+01
RMM: 364    -0.472826066092E+03   -0.44598E-01   -0.22823E-01  1457   0.484E-01    0.133E+01
RMM: 365    -0.472887170516E+03   -0.61104E-01   -0.28589E-01  1483   0.550E-01    0.133E+01
RMM: 366    -0.472950805956E+03   -0.63635E-01   -0.29253E-01  1402   0.609E-01    0.133E+01
RMM: 367    -0.473064189241E+03   -0.11338E+00   -0.30840E-01  1505   0.687E-01    0.133E+01
RMM: 368    -0.473126861399E+03   -0.62672E-01   -0.22038E-01  1442   0.640E-01    0.133E+01
RMM: 369    -0.473137552982E+03   -0.10692E-01   -0.21527E-01  1470   0.596E-01    0.133E+01
RMM: 370    -0.473093405770E+03    0.44147E-01   -0.17493E-01  1425   0.528E-01    0.132E+01
RMM: 371    -0.473157926530E+03   -0.64521E-01   -0.14250E-01  1405   0.469E-01    0.131E+01
RMM: 372    -0.472912186281E+03    0.24574E+00   -0.24341E-01  1470   0.844E-01    0.129E+01
RMM: 373    -0.472972101513E+03   -0.59915E-01   -0.12402E-01  1531   0.444E-01    0.128E+01
RMM: 374    -0.473007867055E+03   -0.35766E-01   -0.10897E-01  1436   0.491E-01    0.128E+01
RMM: 375    -0.473071360584E+03   -0.63494E-01   -0.83733E-02  1430   0.371E-01    0.127E+01
RMM: 376    -0.473051433159E+03    0.19927E-01   -0.89927E-02  1474   0.441E-01    0.127E+01
RMM: 377    -0.472537767627E+03    0.51367E+00   -0.31600E-01  1507   0.101E+00    0.123E+01
RMM: 378    -0.472406459907E+03    0.13131E+00   -0.16694E-01  1539   0.522E-01    0.120E+01
RMM: 379    -0.472536931491E+03   -0.13047E+00   -0.61846E-02  1457   0.395E-01    0.120E+01
RMM: 380    -0.472518992392E+03    0.17939E-01   -0.77421E-02  1459   0.267E-01    0.120E+01
RMM: 381    -0.472530548634E+03   -0.11556E-01   -0.47207E-02  1234   0.336E-01    0.119E+01
RMM: 382    -0.472551727241E+03   -0.21179E-01   -0.87659E-02  1527   0.377E-01    0.119E+01
RMM: 383    -0.472282507297E+03    0.26922E+00   -0.69969E-02  1397   0.500E-01    0.117E+01
RMM: 384    -0.472183383652E+03    0.99124E-01   -0.77764E-02  1517   0.300E-01    0.116E+01
RMM: 385    -0.472171558268E+03    0.11825E-01   -0.49080E-02  1428   0.379E-01    0.116E+01
RMM: 386    -0.472126423448E+03    0.45135E-01   -0.64831E-02  1524   0.274E-01    0.116E+01
RMM: 387    -0.472615471694E+03   -0.48905E+00   -0.16124E-01  1364   0.905E-01    0.118E+01
RMM: 388    -0.472630168076E+03   -0.14696E-01   -0.90133E-02  1399   0.309E-01    0.119E+01
RMM: 389    -0.472620985254E+03    0.91828E-02   -0.66832E-02  1478   0.472E-01    0.118E+01
RMM: 390    -0.472694774151E+03   -0.73789E-01   -0.69886E-02  1519   0.292E-01    0.119E+01
RMM: 391    -0.472689696325E+03    0.50778E-02   -0.29236E-02  1191   0.271E-01    0.119E+01
RMM: 392    -0.472726929558E+03   -0.37233E-01   -0.80382E-02  1528   0.448E-01    0.119E+01
RMM: 393    -0.472959871773E+03   -0.23294E+00   -0.43303E-02  1370   0.420E-01    0.120E+01
RMM: 394    -0.473022187574E+03   -0.62316E-01   -0.46570E-02  1512   0.237E-01    0.121E+01
RMM: 395    -0.472805300001E+03    0.21689E+00   -0.32737E-02  1332   0.372E-01    0.120E+01
RMM: 396    -0.472612665014E+03    0.19263E+00   -0.54967E-02  1518   0.328E-01    0.118E+01
RMM: 397    -0.474859020516E+03   -0.22464E+01   -0.22247E+00  1378   0.328E+00    0.135E+01
RMM: 398    -0.475099282580E+03   -0.24026E+00   -0.84080E-01  1565   0.853E-01    0.144E+01
RMM: 399    -0.474830893083E+03    0.26839E+00   -0.67093E-02  1390   0.557E-01    0.142E+01
RMM: 400    -0.474595688841E+03    0.23520E+00   -0.16782E-01  1548   0.643E-01
  38 F= -.47459569E+03 E0= -.47458491E+03  d E =0.177395E+01
作者
Author:
yamin    时间: 2019-11-25 10:43
题中POSCAR第倒数八行不知道为什么我这里显示不出来,编辑时是完整的。完整数据为 0.5837150147524066  0.6659344414187797  0.5436186222771064   T   T   T

作者
Author:
Accelerator    时间: 2019-11-25 12:22
VASP一些特无论如何都弄不收敛的体系 构型又没有不合理的 可以尝试用不收敛的波函数多优化几步 有望走几步后SCF得以收敛
作者
Author:
yamin    时间: 2019-11-25 16:56
Accelerator 发表于 2019-11-25 12:22
VASP一些特无论如何都弄不收敛的体系 构型又没有不合理的 可以尝试用不收敛的波函数多优化几步 有望走几步 ...

您好,是指读取波函数再优化吗?
作者
Author:
啦啦黑还黑    时间: 2019-11-25 22:25
为啥用ISPIN=1呀,你这体系很可能自旋极化,doped金属原子和Ce4+之间有电子转移,打开自旋ISPIN=2
作者
Author:
yamin    时间: 2019-11-26 17:03
啦啦黑还黑 发表于 2019-11-25 22:25
为啥用ISPIN=1呀,你这体系很可能自旋极化,doped金属原子和Ce4+之间有电子转移,打开自旋ISPIN=2

好,我来试试,谢谢前辈!
作者
Author:
yamin    时间: 2019-11-26 17:16
啦啦黑还黑 发表于 2019-11-25 22:25
为啥用ISPIN=1呀,你这体系很可能自旋极化,doped金属原子和Ce4+之间有电子转移,打开自旋ISPIN=2

完美二氧化铈表面在吸附丙酮和2个H时,ISPIN=1,同样SCF不收敛,单独吸附丙酮收敛,单独吸附H收敛。看文献U值范围4或5,截断能范围400-500(看过一篇300的,找不到出处了),我打算尝试分别修改ISPIN=2,ENCUT=500,U=4试试。
作者
Author:
yamin    时间: 2019-11-27 16:16
啦啦黑还黑 发表于 2019-11-25 22:25
为啥用ISPIN=1呀,你这体系很可能自旋极化,doped金属原子和Ce4+之间有电子转移,打开自旋ISPIN=2

您好,我打开了ISPIN=2,还是不收敛,还能咋整?
作者
Author:
yamin    时间: 2019-11-28 14:43
目前老师给我改了SIGMA=0.05 ALGO=A,收敛中
作者
Author:
yamin    时间: 2019-11-29 01:07
Accelerator 发表于 2019-11-25 12:22
VASP一些特无论如何都弄不收敛的体系 构型又没有不合理的 可以尝试用不收敛的波函数多优化几步 有望走几步 ...

最终SIGMA = 0.05, ALGO = A ,收敛。谢谢前辈指导!
作者
Author:
风飞    时间: 2020-8-28 11:26
yamin 发表于 2019-11-25 10:35
以下是读取SIGMA=0.5时的WAVECAR,用SIGMA=0.3时跑的汇总文件。37-38步de由+00变为+01,我看着能量发散,自洽 ...

您好,请问您搭建CeO2(111)氧空位的步骤是?
我之前切出CeO2(111)面,然后进行优化,然后在111面删除一个O,然后对这个模型进行优化,结果怎么也不收敛。请问是不是我这样的构建的方法有问题
作者
Author:
Penson    时间: 2020-8-30 15:33
体系比较大时,用ALGO = V会更容易收敛,但是可靠性相对弱一些,不过我觉得基本靠谱
作者
Author:
yamin    时间: 2020-8-31 14:28
风飞 发表于 2020-8-28 11:26
您好,请问您搭建CeO2(111)氧空位的步骤是?
我之前切出CeO2(111)面,然后进行优化,然后在111面删除 ...

是这么做的,直接删掉了O原子。你试试调整参数
作者
Author:
风飞    时间: 2020-9-1 22:04
yamin 发表于 2020-8-31 14:28
是这么做的,直接删掉了O原子。你试试调整参数

你好,我是用QE进行计算的,请问您调整的这个参数:SIGMA = 0.05, ALGO = A   代表的什么意思呢?我看看能否借鉴到我的计算中

作者
Author:
yamin    时间: 2020-9-12 20:21
https://www.vasp.at/wiki/wiki/index.php/ALGO
ALGO是控制电子算法的
作者
Author:
syuusama    时间: 2024-8-31 11:29
风飞 发表于 2020-9-1 22:04
你好,我是用QE进行计算的,请问您调整的这个参数:SIGMA = 0.05, ALGO = A   代表的什么意思呢?我看看 ...

SIGMA为展宽大小,ISMEAR和SIGMA这两个参数确定弥散(smearing)的方法和展宽的参数,SIGMA = 0.05通常用于非金属和半导体;ALGO为电子迭代算符,不同的 ALGO 对应的是不同的算法,本质上都是在找势能面上的能量极小值,ALGO = All即IALGO = 58,可以看 VASP wiki 里的介绍。https://www.vasp.at/wiki/index.php/ALGO




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3