24小时热门版块排行榜

>基金申请 (4662)
>论坛更新日志 (4454)
>休闲灌水 (1363)
>虫友互识 (274)
>论文投稿 (257)
>导师招生 (247)
>文献求助 (179)
>有奖起名 (165)
>硕博家园 (97)
>健康生活 (88)
>考博 (87)
>留学DIY (37)
>考研 (33)
>公派出国 (29)
>绿色求助(高悬赏) (28)
>博后之家 (27)

返回列表

gswylq

木虫 (正式写手)

应助: 10 (幼儿园)
金币: 2547.4
散金: 120
红花: 1
帖子: 798
在线: 731小时
虫号: 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

回复此楼

» 猜你喜欢

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

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

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

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

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] 为什么国自然不能直接公布 +3	bjdxyxy 2026-08-26	3/150	2026-08-26 12:48 by 小丑鱼爱读书
[基金申请] 科研孤儿太难了 +21	我4大白菜 2026-08-20	22/1100	2026-08-26 11:29 by nlgza
[基金申请] 怎么查啊 +4	huang1991js 2026-08-26	4/200	2026-08-26 11:26 by panjy.cn
[基金申请] 国合里面能看到了 +7	一怀馨秋 2026-08-26	7/350	2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4	yanglien 2026-08-26	5/250	2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6	wittyboy 2026-08-26	14/700	2026-08-26 10:55 by wittyboy
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3	鲸鱼融金_浙江_� 2026-08-22	3/150	2026-08-26 09:01 by zzuzxg
[基金申请] 出来了 +8	trojank 2026-08-26	8/400	2026-08-26 09:00 by gxc
[基金申请] 国合现在查不到了吗？ +10	chengyan1220 2026-08-24	20/1000	2026-08-26 08:57 by peasantsprig
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +24	木水思豆 2026-08-25	27/1350	2026-08-26 08:26 by yanchen918
[基金申请] 有没有大神帮我看看基金代码 28+4	1234567wang 2026-08-24	10/500	2026-08-25 19:15 by lfy8008
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9	yuleib84 2026-08-25	10/500	2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8	scalable 2026-08-24	8/400	2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9	图啦图啦 2026-08-24	10/500	2026-08-25 11:59 by 南海小哥
[基金申请] 我面上完蛋了 +13	且听虎啸 2026-08-20	14/700	2026-08-25 09:10 by mrkang
[基金申请] 建议基金发布提前给出明确的时间点 +13	kulium 2026-08-21	16/800	2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧！ +13	大萍1987 2026-08-20	16/800	2026-08-24 10:23 by 太傻了
[基金申请] 时间戳又变了 +13	wuchongjun 2026-08-20	19/950	2026-08-21 17:21 by 紫杉醇
[基金申请] 今天放榜没戏了吧 +9	yuleib84 2026-08-19	11/550	2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4	longfie172 2026-08-20	4/200	2026-08-21 08:58 by mark mao