计算化学公社

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

[CP2K] CP2K联用plumed报错:无法传递电荷信息

[复制链接 Copy URL]

44

帖子

0

威望

309

eV
积分
353

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
各位老师大家好,
######最近我使用CP2K联用plumed进行元动力学采样,想尝试用偶极矩作为CV,但是发现plumed出现如下报错(无法将电荷信息传递给plumed):
PLUMED: ERROR in input to action DIPOLE with label d : charges were not passed to plumed
PLUMED:
PLUMED: The above exception was the direct cause of the following exception:
PLUMED:
PLUMED: (core/PlumedMain.cpp:972) void PLMD::PlumedMain::justCalculate()
PLUMED: An error happened while calculating d

######以下是我的plumed.dat:

#定义原子
c: GROUP ATOMS=2
f: GROUP ATOMS=12
h: GROUP ATOMS=726-728,405-407,843-845,801-803,468-470,396-398

#定义配位数
ccf: COORDINATION GROUPA=c GROUPB=f R_0=0.18 NN=8 MM=14
d: DIPOLE GROUP=h
INCULDE FILE=HIRSHFELD.CHG
# Metadynamics on multiple dihedral angles
METAD ...
  LABEL=metad
  ARG=ccf,d
  SIGMA=0.1,0.1
  HEIGHT=0.5
  PACE=20
  BIASFACTOR=10
  TEMP=300
  GRID_MIN=0,-1
  GRID_MAX=1,1
  GRID_BIN=50,50
...
# 输出内容
PRINT STRIDE=5 ARG=ccf,d,metad.bias FILE=COLVAR

######以下是CP2K输入文件输出电荷的参数:
    &PRINT
      &HIRSHFELD
        FILENAME =HIRSHFELD.CHG
        SHAPE_FUNCTION DENSITY
        &EACH
          MD 1
        &END EACH        
      &END HIRSHFELD   
    &END PRINT

本版积分规则 Credits rule

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

GMT+8, 2025-8-14 18:19 , Processed in 0.192559 second(s), 20 queries , Gzip On.

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