24小时热门版块排行榜

返回列表

wenhuakul

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 63
帖子: 5
在线: 9.7小时
虫号: 4135464
注册: 2015-10-11
专业: 物理学I

[求助] 小程序Murn.f拟合晶格常数结果出现NaN如何解决？

初学vasp，在用程序murn.f进行拟合晶格常数时，有的能拟合出来，有的拟合不出来，拟合不出来的情况是输出文件out.m中出现NaN，不知道是什么原因？下面是我的输入和输出文件：
输入文件in.m

2
0.25
5.20 5.70 40
11
5.20 -11.805731
5.25 -11.874634
5.30 -11.920540
5.35 -11.945297
5.40 -11.950657
5.45 -11.938253
5.50 -11.909634
5.55 -11.866236
5.60 -11.809424
5.65 -11.740470
5.70 -11.660566

输出文件out.m

least-squares fit of etot(v) by Murnaghan eq.

         UNIT OF LENGTH: ULA =    0.529177  ANGSTROEM
         UNIT OF ENERGY: IUNIT = 2:
                                       (ENERGIES INPUT IN ELECTRONVOLTS)

IN THE STRUCTURE CONSIDERED, THE UNIT CELL VOLUME = A**3 *    0.250000

min and max alat(bohr) to compute Etot for at  40 points =          5.2000000 and          5.7000000
      NUMBER OF DATA POINTS: N = 11
-------------------------------------------------------------------------------

   I    ALATT       VOLUME       ENERGY (EV)          ENERGY (RY)

   1    5.200000    35.152000    -11.8057310000    -11.8057310000
   2    5.250000    36.175781    -11.8746340000    -11.8746340000
   3    5.300000    37.219250    -11.9205400000    -11.9205400000
   4    5.350000    38.282594    -11.9452970000    -11.9452970000
   5    5.400000    39.366000    -11.9506570000    -11.9506570000
   6    5.450000    40.469656    -11.9382530000    -11.9382530000
   7    5.500000    41.593750    -11.9096340000    -11.9096340000
   8    5.550000    42.738469    -11.8662360000    -11.8662360000
   9    5.600000    43.904000    -11.8094240000    -11.8094240000
10    5.650000    45.090531    -11.7404700000    -11.7404700000
11    5.700000    46.298250    -11.6605660000    -11.6605660000
-------------------------------------------------------------------------------

   ALATT0 VOL0 (ULA**3)  B0 (MBAR)    B0PRIME    E0 (EV)    E0 (RY)

      5.4000    39.36600    1.00000    4.00000 -11.95066    0.00000
                                                            (=INITIAL GUESS)

MIN: 2.7000    4.92075    0.01000    0.01000
MAX: 8.1000 132.86025    10.00000    10.00000
-------------------------------------------------------------------------------

FIT BY MURNAGHAN EQUATION EQUATION OF STATE
-------------------------------------------------------------------------------

ITER ALATT0 VOL0 (ULA**3)  B0 (MBAR)    B0PRIME    E0 (EV) SUM OF SQUARES
                                                                        IERR

20    NaN       NaN       NaN       NaN       NaN       NaN
                                                                           10
40    NaN       NaN       NaN       NaN       NaN       NaN
                                                                           0
40    NaN       NaN       NaN       NaN       NaN       NaN
                                                                           0
-------------------------------------------------------------------------------

   ALATT0 VOL0 (ULA**3)  B0 (MBAR)    B0PRIME    E0 (EV)    E0 (RY)

alat=    NaN b0=    NaN b0p=    NaN E0=       NaN
         NaN       NaN       NaN       NaN       NaN       NaN
         NaN
-------------------------------------------------------------------------------

      alat(A)       energy(Ryd)          alat(bohr)       energy(H)
-------------------------------------------------------------------------------

      2.75172                NaN       5.20000                NaN
      2.75851                NaN       5.21282                NaN
      2.76529                NaN       5.22564                NaN
      2.77207                NaN       5.23846                NaN
      2.77886                NaN       5.25128                NaN
      2.78564                NaN       5.26410                NaN
      2.79243                NaN       5.27692                NaN
      2.79921                NaN       5.28974                NaN
      2.80600                NaN       5.30256                NaN
      2.81278                NaN       5.31538                NaN
      2.81956                NaN       5.32821                NaN
      2.82635                NaN       5.34103                NaN
      2.83313                NaN       5.35385                NaN
      2.83992                NaN       5.36667                NaN
      2.84670                NaN       5.37949                NaN
      2.85349                NaN       5.39231                NaN
      2.86027                NaN       5.40513                NaN
      2.86705                NaN       5.41795                NaN
      2.87384                NaN       5.43077                NaN
      2.88062                NaN       5.44359                NaN
      2.88741                NaN       5.45641                NaN
      2.89419                NaN       5.46923                NaN
      2.90098                NaN       5.48205                NaN
      2.90776                NaN       5.49487                NaN
      2.91454                NaN       5.50769                NaN
      2.92133                NaN       5.52051                NaN
      2.92811                NaN       5.53333                NaN
      2.93490                NaN       5.54615                NaN
      2.94168                NaN       5.55897                NaN
      2.94847                NaN       5.57179                NaN
      2.95525                NaN       5.58462                NaN
      2.96204                NaN       5.59744                NaN
      2.96882                NaN       5.61026                NaN
      2.97560                NaN       5.62308                NaN
      2.98239                NaN       5.63590                NaN
      2.98917                NaN       5.64872                NaN
      2.99596                NaN       5.66154                NaN
      3.00274                NaN       5.67436                NaN
      3.00953                NaN       5.68718                NaN
      3.01631                NaN       5.70000                NaN

回复此楼