| 查看: 1352 | 回复: 0 | |||
gswylq木虫 (正式写手)
|
[求助]
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人回复
» 本主题相关商家推荐: (我也要在这里推广)
