计算化学公社

 找回密码 Forget password
 注册 Register
Views: 2455|回复 Reply: 5
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] cp2k-7.1做MetaD问题求助

[复制链接 Copy URL]

24

帖子

0

威望

233

eV
积分
257

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
  各位老师,最近我在用cp2k-7.1跑MetaD,我定义了一个cv值,为Si-H键和H-O键键长的差值,跑了一段得到了cp2k-COLVAR.metadynLog和cp2k-HILLS.metadynLog两个输出文件。目前的一些问题是:(1).不会将输出数据转化为自由能随cv的变化文件,看了一些帖子,在cp2k/src/metadyn_tools里面有附带的处理文件,不过我发现7.1的版本对应的文件名不一样而且我调用不了,共graph.F、graph_methods.F、graph_utils.F和PACKAGE四个文件,没明白是用不了还是我不会使用的问题,加了权限了也不能用。(2).我看了一些文献中,出现了局部极小值的自由能对应的位点,如何清楚地定位在这各能量下对应的结构信息呢。(3).因为是初次尝试,尚不清楚定义CV的位置是否合理,还想了解一些定义的技巧。
自学太难了,还请各位老师能够不吝赐教。

24

帖子

0

威望

233

eV
积分
257

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2022-6-9 08:22:50 | 只看该作者 Only view this author
yjb 发表于 2022-6-8 21:57
可以用9.1,社长有些怎么编译,很方便

好的  我也准备换个版本试试

173

帖子

0

威望

1800

eV
积分
1973

Level 5 (御坂)

5#
发表于 Post on 2022-6-8 21:57:10 | 只看该作者 Only view this author
予拥12138 发表于 2022-6-8 19:47
报错内容如下:(1). cp2k-7.1/src/metadyn_tools/graph.F: line 2: FES:: command not found
(2). cp2k- ...

可以用9.1,社长有些怎么编译,很方便
独立之精神 自由之思想

24

帖子

0

威望

233

eV
积分
257

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2022-6-8 19:47:46 | 只看该作者 Only view this author
yjb 发表于 2022-6-8 18:34
1)  graph.psmp -ndim 1 -ndw 1 -file METADYN-1_10000.restart  -cp2k 可以试着把graph.psmp换成graph.F等 ...

报错内容如下:(1). cp2k-7.1/src/metadyn_tools/graph.F: line 2: FES:: command not found
(2). cp2k-7.1/src/metadyn_tools/graph.F: line 3: syntax error near unexpected token `C'
(3) cp2k-7.1/src/metadyn_tools/graph.F: line 3: `!   Copyright (C) 2002,2003,2004,2005,2006,2007,2008,2009,2010,2011,2012,2013
不能理解,哎。

24

帖子

0

威望

233

eV
积分
257

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2022-6-8 19:21:32 | 只看该作者 Only view this author
yjb 发表于 2022-6-8 18:34
1)  graph.psmp -ndim 1 -ndw 1 -file METADYN-1_10000.restart  -cp2k 可以试着把graph.psmp换成graph.F等 ...

谢谢您,您的帖子我看过,但是我给与权限并运行graph.F、graph_methods.F都不成功。另外,关于计算速度的问题,我是232个原子的体系,稳定状态下,1个离子步跑1min,,时间是我跑普通AIMD的两倍,不知道这个速度是否有改善的可能。

173

帖子

0

威望

1800

eV
积分
1973

Level 5 (御坂)

2#
发表于 Post on 2022-6-8 18:34:58 | 只看该作者 Only view this author
1)  graph.psmp -ndim 1 -ndw 1 -file METADYN-1_10000.restart  -cp2k 可以试着把graph.psmp换成graph.F等试试。 2)graph.psmp -ndim 1 -ndw 1  -file METADYN-1_10000.restart -find-minima -cp2k 会生成全部的极小值。3)去看Review。
独立之精神 自由之思想

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-20 04:00 , Processed in 0.184232 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list