计算化学公社

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

[VASP] VASP分子动力学ICONST怎么直接设置原子坐标作为变量

[复制链接 Copy URL]

28

帖子

0

威望

553

eV
积分
581

Level 4 (黑子)

通过ICONST文件指定原子间的距离或者角度,作为collective variable,可以做各种增强采样方法研究反应,这个没问题。
但我目前单纯想要让某个原子沿着z轴方向走,看vasp的Wiki似乎可以设置坐标(flag=X, Y, and Z; cX, cY, and cZ)作为变量,但我试了下,在ICONST里设置:cZ 7 0,感觉像是把7号原子的z坐标作为变量了,但似乎不行,vasp跑了一个离子步就终止了,是不是我的理解有问题,不是这么设置的

5

帖子

0

威望

117

eV
积分
122

Level 2 能力者

2#
发表于 Post on 2023-8-14 14:44:44 | 只看该作者 Only view this author
您好,如果体系中有真空层或许可以参考文献J.Am.Chem. Soc. 2022, 144, 17140−17148,图S1中提到的方法。另外,请问,您是否已有设置坐标的思路。

28

帖子

0

威望

553

eV
积分
581

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2023-8-14 20:05:11 | 只看该作者 Only view this author
deng111 发表于 2023-8-14 14:44
您好,如果体系中有真空层或许可以参考文献J.Am.Chem. Soc. 2022, 144, 17140−17148,图S1中提到的方 ...

巧了,一开始就是学习的这篇文章里的方法
在POSCAR里明明固定了氦原子,设置好坐标一跑,氦原子也跟着动了,不明白什么原因

5

帖子

0

威望

117

eV
积分
122

Level 2 能力者

4#
发表于 Post on 2023-8-15 15:30:11 | 只看该作者 Only view this author
cX 61 0.62247 0 。我参照你的方法试了下,加不加最后的0似乎都可以跑,但是跑的方向反了

28

帖子

0

威望

553

eV
积分
581

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2023-8-16 17:07:37 | 只看该作者 Only view this author
deng111 发表于 2023-8-15 15:30
cX 61 0.62247 0 。我参照你的方法试了下,加不加最后的0似乎都可以跑,但是跑的方向反了

你好请问0.62247是什么参数啊

5

帖子

0

威望

117

eV
积分
122

Level 2 能力者

6#
发表于 Post on 2023-8-17 11:06:39 | 只看该作者 Only view this author
终态的分数坐标,加了看他会不会跑到目标坐标上,不过似乎是失败的

28

帖子

0

威望

553

eV
积分
581

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2023-10-16 20:33:11 | 只看该作者 Only view this author
更新,问题解决。最近想到这个问题又试了下,发现设置笛卡尔坐标cX、cY、cZ没用,提示没有这个flag,换成分数坐标X、Y、Z一切正常,结果符合预期!感觉是我的vasp版本问题?

2

帖子

0

威望

61

eV
积分
63

Level 2 能力者

8#
发表于 Post on 2024-10-8 11:11:57 | 只看该作者 Only view this author
yxzinc 发表于 2023-10-16 20:33
更新,问题解决。最近想到这个问题又试了下,发现设置笛卡尔坐标cX、cY、cZ没用,提示没有这个flag,换成分 ...

您好想问一下,具体是怎么解决的,我按照您说的试了一下发现还是不行,需要在计算中移动的原子一直不发生移动

197

帖子

1

威望

1108

eV
积分
1325

Level 4 (黑子)

9#
发表于 Post on 2024-10-8 13:50:42 | 只看该作者 Only view this author
VASP做这种限制性动力学太不好用了,所以我一般都是CP2K+PLUMED来解决

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 09:33 , Processed in 0.370348 second(s), 24 queries , Gzip On.

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