24小时热门版块排行榜

返回列表

gswylq

木虫 (正式写手)

应助: 10 (幼儿园)
金币: 2546.6
散金: 120
红花: 1
帖子: 782
在线: 730.1小时
虫号: 659377
注册: 2008-11-21
专业: 金属材料的凝固与结晶学

[求助] cpmd metadynamic tutorial day4中的SI54不能完成

请计算过的老师、同学指导一下：
最近我重复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@oxox6085@smutao

回复此楼

» 猜你喜欢

有没有人能给点建议已经有5人回复
假如你的研究生提出不合理要求已经有12人回复
实验室接单子已经有7人回复
全日制（定向）博士已经有5人回复
萌生出自己或许不适合搞科研的想法，现在跑or等等看？已经有4人回复
Materials Today Chemistry审稿周期已经有4人回复
参与限项已经有3人回复
对氯苯硼酸纯化已经有3人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复

» 本主题相关商家推荐: (我也要在这里推广)

1楼 2018-03-05 15:37:11

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

相关版块跳转我要订阅楼主 gswylq 的主题更新

返回列表