cpmd metadynamic tutorial day4中的SI54不能完成

gswylq

请计算过的老师、同学指导一下：
最近我重复cpmd metadynamic tutorial day4中的SI54，META_CELL里面的inp2文件（如下输入文件1）不能完成计算，我的计算步骤为：
1. 先将下面输入文件里面&CPMD部分修改，执行波函数优化（MAXSTEP 200，去掉METASTEPNUM部分）；
2. 然后进行分子动力学模拟（&CPMD 部分设置为 MOLECULAR DYNAMICS  RESTART WAVEFUNCTIONS COORDINATES  CELL LATEST MAXSTEP=2000，去掉METASTEPNUM部分）；
3. 上一步分子动力学跑完后，将&CPMD 部分重新设置(MOLECULAR DYNAMICS   RESTART WAVEFUNCTIONS COORDINATES ACCUMULATORS VELOCITIES CELL LATEST  QUENCH BO  PARRINELLO-RAHMAN， MAXSTEP  60000), 并加入METASTEPNUM部分，结果发现分子动力学跑完一部分后停止了，说“LATEST文件读入错误”。

我想问一下，我这个计算步骤是不是有问题？烦请指导一下。
我发现第三步，如果不加入METASTEPNUM部分，分子动力学一直会计算，可是加入了METASTEPNUM部分后，只计算一会就报错。是不是第2步的MAXSTEP和第三步的MAXSTEP需要相同？都要设置为60000？


#输入文件1 ####
&CPMD
MOLECULAR DYNAMICS
RESTART WAVEFUNCTIONS COORDINATES ACCUMULATORS VELOCITIES CELL LATEST
QUENCH BO
PARRINELLO-RAHMAN
TEMPERATURE IONS
300.
TEMPCONTROL IONS
300. 30.
MAXSTEP
60000
TRAJECTORY XYZ  SAMPLE
50
STRESS TENSOR
10
CMASS
250000
EMASS
800
  SPLINE RANGE
  3.0
MAXCPUTIME
85000
&END

&SYSTEM
  CELL VECTORS
  -15.3878  0.00 15.3878  0.00 15.3878  15.3878  -15.3878  15.3878  0.00
  CUTOFF
  20.
  PRESSURE
  95.
&END

&DFT
   FUNCTIONAL LDA
&END

&ATOMS
*SI_MT_LDA  KLEINMAN-BYLANDER
LMAX=P
54
     0.00000      0.00000      0.00000
    -5.12927      5.12927      0.00000
   -10.25855     10.25855      0.00000
     0.00000      5.12927      5.12927
    -5.12927     10.25855      5.12927
   -10.25855     15.38782      5.12927
     0.00000     10.25855     10.25855
    -5.12927     15.38782     10.25855
   -10.25855     20.51710     10.25855
    -5.12927      0.00000      5.12927
   -10.25855      5.12927      5.12927
   -15.38782     10.25855      5.12927
    -5.12927      5.12927     10.25855
   -10.25855     10.25855     10.25855
   -15.38782     15.38782     10.25855
    -5.12927     10.25855     15.38782
   -10.25855     15.38782     15.38782
   -15.38782     20.51710     15.38782
   -10.25855      0.00000     10.25855
   -15.38782      5.12927     10.25855
   -20.51710     10.25855     10.25855
   -10.25855      5.12927     15.38782
   -15.38782     10.25855     15.38782
   -20.51710     15.38782     15.38782
   -10.25855     10.25855     20.51710
   -15.38782     15.38782     20.51710
   -20.51710     20.51710     20.51710
     2.56464      2.56464      2.56464
    -2.56464      7.69391      2.56464
    -7.69391     12.82319      2.56464
     2.56464      7.69391      7.69391
    -2.56464     12.82319      7.69391
    -7.69391     17.95246      7.69391
     2.56464     12.82319     12.82319
    -2.56464     17.95246     12.82319
    -7.69391     23.08174     12.82319
    -2.56464      2.56464      7.69391
    -7.69391      7.69391      7.69391
   -12.82319     12.82319      7.69391
    -2.56464      7.69391     12.82319
    -7.69391     12.82319     12.82319
   -12.82319     17.95246     12.82319
    -2.56464     12.82319     17.95246
    -7.69391     17.95246     17.95246
   -12.82319     23.08174     17.95246
    -7.69391      2.56464     12.82319
   -12.82319      7.69391     12.82319
   -17.95246     12.82319     12.82319
    -7.69391      7.69391     17.95246
   -12.82319     12.82319     17.95246
   -17.95246     17.95246     17.95246
    -7.69391     12.82319     23.08174
   -12.82319     17.95246     23.08174
   -17.95246     23.08174     23.08174


META DYNAMICS COLLECTIVE VARIABLE
DEFINE VARIABLES
3
CELLSIDE   1  SCF 1. KCV 1.0  MCV 100.0
CELLSIDE   2  SCF 1. KCV 1.0  MCV 100.0
CELLSIDE   3  SCF 1. KCV 1.0  MCV 100.0
END  DEFINE

METASTEPNUM
800
META_RESTART
1657
HILLS SHIFT  RCUT  2.0  1.0  = 0.3  0.005
TUNING HHEIGHT  = 0.0010 0.09
LAGRANGE  TEMPERATURE
500.00
LAGRANGE TEMPCONT
500.00 100.00
MAXSTEPNUM  INTERMETA
200
MINSTEPNUM  INTERMETA
70
METASTORE
20  1 1000
MOVEMENT CHECK
0.01
MAXKINEN
0.07
RESTRAIN VOLUME
0.0  1.1  2.0
END METADYNAMICSMETA

&END

多次拒绝斯佳丽

你好 能不能分享一下cpmd metadynamic 4-days tutorial，网上资源很多都被删除了，感激不尽