计算化学公社

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

[CP2K] 用metadynamic (cp2k)求体系自由能的metadynlog文件内容问题

[复制链接 Copy URL]

1

帖子

0

威望

13

eV
积分
14

Level 1 能力者

跳转到指定楼层 Go to specific reply
楼主
请问用metadynamic(cp2k)求体系自由能,目前有metadynlog文件,2个CV。 col2 和col3 是两个cv。后面6个col是什么? 请问可以用这个文件求出来自由能吗?
      0.00000000      3.31399      5.44380      0.00000      0.00000      0.00000      0.00000      0.00000      0.00000
      0.50000000      3.31267      5.44584      0.00000      0.00000      0.00000      0.00000      0.00000      0.00000
      1.00000000      3.31030      5.44884      0.00000      0.00000      0.00000      0.00000      0.00000      0.00000
      1.50000000      3.30683      5.45288      0.00000      0.00000      0.00000      0.00000      0.00000      0.00000
      2.00000000      3.30219      5.45802      0.00000      0.00000      0.00000      0.00000      0.00000      0.00000
      2.50000000      3.29635      5.46428      0.00000      0.00000      0.00000      0.00000      0.00000      0.00000
      3.00000000      3.28929      5.47167      0.00000      0.00000      0.00000      0.00000      0.00000      0.00000
      3.50000000      3.28103      5.48010      0.00000      0.00000      0.00000      0.00000      0.00000      0.00000
      4.00000000      3.27160      5.48947      0.00000      0.00000      0.00000      0.00000      0.00000      0.00000
      4.50000000      3.26106      5.49961      0.00000      0.00000      0.00000      0.00000      0.00000      0.00000
      5.00000000      3.24951      5.51038      0.00000      0.00000      0.00000      0.00000      0.00000      0.00000
      5.50000000      3.23700      5.52158      0.00000      0.00000      0.00000      0.00000      0.00000      0.00000
      6.00000000      3.22360      5.53306      0.00000      0.00000      0.00000      0.00000      0.00000      0.00000


6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

2#
发表于 Post on 2020-3-30 20:30:13 | 只看该作者 Only view this author
如置顶的新社员必读贴和论坛首页的公告栏所示,求助帖必须在帖子标题明确体现出此帖内容是求助或提问,并反映出帖子具体内容,仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“aimd 自由能 metadynamic”改了,以后务必注意
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

186

帖子

1

威望

4409

eV
积分
4615

Level 6 (一方通行)

3#
发表于 Post on 2020-3-31 12:45:19 | 只看该作者 Only view this author
请参考 https://www.cp2k.org/exercises:2015_cecam_tutorial Second task 的最后一段:
Other quantities that can be monitored from the MOTION/FREE_ENERGY/METADYN/PRINT/COLVAR output, beside the instantaneous values of the three CVs (2nd, 3rd, and 4th col.) are : the instantaneous gradient of the bias potential computed with respect to CV (5th, 6th, 7th col.), the gradients wit respect to the CVs of wall potentials, if present, (8th,9th,10th col), the instantaneous value of the bias potential, and the instantaneous values of the wall potentials.

1

帖子

0

威望

63

eV
积分
64

Level 2 能力者

4#
发表于 Post on 2021-1-30 18:22:50 | 只看该作者 Only view this author
请问您体系的自由能求值问题解决了吗?我也遇到了这个问题。

13

帖子

0

威望

95

eV
积分
108

Level 2 能力者

5#
发表于 Post on 2025-8-1 16:33:56 | 只看该作者 Only view this author
请问解决了吗,我也想知道

本版积分规则 Credits rule

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

GMT+8, 2025-8-13 05:49 , Processed in 0.214553 second(s), 26 queries , Gzip On.

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