计算化学公社

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

[CP2K] AIMD中GHOST原子不固定

[复制链接 Copy URL]

44

帖子

0

威望

314

eV
积分
358

Level 3 能力者

各位老师同学们大家好,请问为什么我在CP2K任务中通过以下模块设GHOST原子来得到固定坐标参照:
&COORD
    X x y z
&END COORD

&KIND X
    ELEMENT H
    MASS 0
    GHOST T
    BASIS_SET  NONE
&END KIND

跑aimd的时候发现第零帧.xyz还在我要求的固定点中,但是在第一帧及以后,虚原子为什么就变成0 0 0坐标?

此外我还尝试了设置&MOTION中的&CONSTRAINS模块固定原子,结果仍然是一样:
#  &CONSTRAINT
#     &FIXED_ATOMS
#      COMPONENTS_TO_FIX XYZ        
#      LIST 924
#    &END FIXED_ATOMS
#  &END CONSTRAINT

541

帖子

5

威望

3445

eV
积分
4086

Level 6 (一方通行)

2#
发表于 Post on 2025-4-7 18:23:56 | 只看该作者 Only view this author
文中的
“此外我还尝试了设置&MOTION中的&CONSTRAINS模块固定原子,结果仍然是一样:
#  &CONSTRAINT
#     &FIXED_ATOMS
#      COMPONENTS_TO_FIX XYZ        
#      LIST 924
#    &END FIXED_ATOMS
#  &END CONSTRAINT”
这个部分为什么带着注释,实际跑的时候注释删掉了么?

44

帖子

0

威望

314

eV
积分
358

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2025-4-7 20:00:14 | 只看该作者 Only view this author
Stardust0831 发表于 2025-4-7 18:23
文中的
这个部分为什么带着注释,实际跑的时候注释删掉了么?

注释的意思是我尝试过的意思,我知道注释是不读入的意思

1412

帖子

6

威望

3150

eV
积分
4682

Level 6 (一方通行)

傻傻的木瓜

4#
发表于 Post on 2025-4-7 20:22:26 | 只看该作者 Only view this author
固定坐标的点应该不需要设为ghost atom才能用于定义其他集合变量吧,手册搜FIX_POINT可以看到很多colvar都支持定义这个代替原子几何中心GEO_CENTER的。有可能鬼原子的设计只适合计算单点能时做BSSE校正用,不完全兼容涉及原子坐标变化的任务,但是还得找到源代码相应部分才能确认。

另外,求不删此帖……不知道为啥楼主好多以前发的帖子都看不了了,感觉对辛苦码字回答问题、参与讨论的其他人不是很友好。
√546=23.36664289109

44

帖子

0

威望

314

eV
积分
358

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2025-4-7 20:37:19 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-4-7 20:22
固定坐标的点应该不需要设为ghost atom才能用于定义其他集合变量吧,手册搜FIX_POINT可以看到很多colvar都 ...

好的 谢谢回答我再尝试一下。删帖是因为之前的一些问题没有什么参考意义,都是因为一些不仔细问的。

44

帖子

0

威望

314

eV
积分
358

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2025-4-7 20:41:03 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-4-7 20:22
固定坐标的点应该不需要设为ghost atom才能用于定义其他集合变量吧,手册搜FIX_POINT可以看到很多colvar都 ...

此外,我之前尝试过定义COLVAR中配位数相关的&POINT/设为FIX_POINT,但貌似不起作用

44

帖子

0

威望

314

eV
积分
358

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2025-4-7 20:57:43 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-4-7 20:22
固定坐标的点应该不需要设为ghost atom才能用于定义其他集合变量吧,手册搜FIX_POINT可以看到很多colvar都 ...

例如我设置了GEO_CENTER,我应该如何在ATOMS_FROM后面引用所设置的GEO_CENTER呢,输出文件总是报错说引用错误。另外我会把有意义的帖子留下,之前的问题也并没有得到解决,所以就删帖了

342

帖子

0

威望

1490

eV
积分
1832

Level 5 (御坂)

8#
发表于 Post on 2026-7-10 18:57:58 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-4-7 20:22
固定坐标的点应该不需要设为ghost atom才能用于定义其他集合变量吧,手册搜FIX_POINT可以看到很多colvar都 ...

你好,我最近想做一些电场相关的计算。我看到很多文章在研究电场(不是外部电场,而是自发形成的电场)对某个共价键影响时,会计算环境原子(例如水中的水分子)产生的电场投影在该键上的分量。这时候为了避免成键原子自电场影响,会把成键原子改成GHOST。我在CP2K手册里面看到GHOST可以选择保留基组和不保留基组两种,想问下这种情况下要不要保留基组呢

1412

帖子

6

威望

3150

eV
积分
4682

Level 6 (一方通行)

傻傻的木瓜

9#
发表于 Post on 2026-7-10 20:14:11 | 只看该作者 Only view this author
ljh123 发表于 2026-7-10 18:57
你好,我最近想做一些电场相关的计算。我看到很多文章在研究电场(不是外部电场,而是自发形成的电场)对 ...

说实话我没看懂什么叫“电场投影在该键上的分量”,化学键既有方向又有长度,电场往该方向投影以后在沿着化学键上每个点的大小未必相等,按哪个点算你说的“分量”?

另外,为什么不能直接把成键原子完全去掉,算一个只保留环境分子的体系的电场,再在后处理的时候用键的端点(成键原子)坐标来做数学运算?
√546=23.36664289109

342

帖子

0

威望

1490

eV
积分
1832

Level 5 (御坂)

10#
发表于 Post on 2026-7-10 20:59:40 | 只看该作者 Only view this author
本帖最后由 ljh123 于 2026-7-10 21:05 编辑
Uus/pMeC6H4-/キ 发表于 2026-7-10 20:14
说实话我没看懂什么叫“电场投影在该键上的分量”,化学键既有方向又有长度,电场往该方向投影以后在沿着 ...

https://pubs.acs.org/doi/10.1021/ja403917z
参考文献是这个您可以看看(我找了好几篇论文,大家都是这么处理),他在SI里也给了处理脚本。
按照我的理解,例如C=O键,就是把电场投影到C→O的单位矢量上,端点分别在羰基氧和碳上,然后取平均。至于为什么不直接删,主要是论文作者就是用ghost原子的。所以我想问问要不要保留基函数,如果不保留的话应该等价于直接删去?
(题外话,不管需不需要保留基函数我可能都会选择用ghost,因为CP2K中这样能最大程度保留我的原子序号一一对应,便于后处理)

1412

帖子

6

威望

3150

eV
积分
4682

Level 6 (一方通行)

傻傻的木瓜

11#
发表于 Post on 2026-7-12 11:45:16 | 只看该作者 Only view this author
ljh123 发表于 2026-7-10 20:59
https://pubs.acs.org/doi/10.1021/ja403917z
参考文献是这个您可以看看(我找了好几篇论文,大家都是这 ...

这篇文章用了GROMACS做经典力场分子动力学,看4.2节的介绍电场对体系的影响体现在原子电荷上面,直接拿溶剂原子电荷不为0和为0的情况下溶质原子 (C或O) 的总受力之差作为静电作用产生的力,除以原子电荷作为原子感受到的“电场” ["The subtraction of all nonelectrostatic forces (as determined by rerunning the trajectory with the charge-neutralized topology) from the total force results in a force exerted on a particular atom due only to electrostatic interactions (eq 2). The electrostatic force can be converted into an electrostatic field (eq 3) simply by dividing by the partial charge of the atom in question (either C or O of acetophenone or p-Ac-Phe)."] ,再平均两个原子的“电场”作为键的“电场”。说实话我认为这种处理方式不太能类推到量子化学计算上。

至于问题本身,鬼原子带基函数主要是用在BSSE计算里,处理原子中心基函数的不完备性误差对相互作用能的影响的。有时候中性分子获得过剩电子产生的的temporary anion state可能也需要额外基函数来表现弥散的电子分布。除此之外,如果鬼原子仅作为三维空间的坐标定位用、不应该对量子化学计算产生什么影响,我觉得没太多必要带基函数。
√546=23.36664289109

342

帖子

0

威望

1490

eV
积分
1832

Level 5 (御坂)

12#
发表于 Post on 2026-7-12 17:07:05 | 只看该作者 Only view this author
本帖最后由 ljh123 于 2026-7-12 17:16 编辑
Uus/pMeC6H4-/キ 发表于 2026-7-12 11:45
这篇文章用了GROMACS做经典力场分子动力学,看4.2节的介绍电场对体系的影响体现在原子电荷上面,直接拿溶 ...

我想的是,不用原子电荷来计算电场,用静电势的负梯度。在计算得到静电势的格点数据后,再计算某点某方向的数值梯度。
https://pubs.acs.org/doi/10.1021/ja2069592 例如这篇

本版积分规则 Credits rule

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

GMT+8, 2026-7-29 13:03 , Processed in 0.214001 second(s), 20 queries , Gzip On.

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