版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

gswylq

木蟲 (正式寫手)

應助: 10 (幼兒園)
金幣: 2546.6
散金: 120
紅花: 1
帖子: 782
在線: 730.1小時
蟲號: 659377
注冊: 2008-11-21
專業(yè): 金屬材料的凝固與結(jié)晶學

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

請計算過的老師、同學指導一下：
最近我重復cpmd metadynamic tutorial day4中的SI54，META_CELL里面的inp2文件（如下輸入文件1）不能計算，我的計算步驟為：
1. 先將下面輸入文件里面&CPMD部分修改，執(zhí)行波函數(shù)優(yōu)化（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部分，結(jié)果發(fā)現(xiàn)分子動力學跑完一部分后停止了，說“LATEST文件讀入錯誤”。

我想問一下，我這個計算步驟是不是有問題？煩請指導一下。
我發(fā)現(xiàn)第三步，如果不加入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

回復此樓

» 猜你喜歡

北科281學碩材料求調(diào)劑已經(jīng)有13人回復
289求調(diào)劑已經(jīng)有11人回復
網(wǎng)絡空間安全0839招調(diào)劑已經(jīng)有5人回復
材料277求調(diào)劑已經(jīng)有4人回復
274求調(diào)劑已經(jīng)有7人回復
085600 材料與化工 329分求調(diào)劑已經(jīng)有3人回復
302求調(diào)劑已經(jīng)有3人回復
312求調(diào)劑已經(jīng)有5人回復
329求調(diào)劑已經(jīng)有3人回復
【2026考研調(diào)劑】制藥工程 284分求相關專業(yè)調(diào)劑名額已經(jīng)有8人回復

» 本主題相關商家推薦: (我也要在這里推廣)

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉(zhuǎn) 我要訂閱樓主 gswylq 的主題更新