计算化学公社

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

[CP2K] 求助:用CP2K跑NVT,采用周期性边界条件时原子跑出盒子

[复制链接 Copy URL]

10

帖子

0

威望

139

eV
积分
149

Level 2 能力者

本帖最后由 we12306 于 2023-4-28 21:10 编辑

我在使用xTB跑石墨晶体中夹杂Xe原子的NVT时,发现明明是周期性边界条件,而越过边界的Xe原子没有到盒子另一侧,而是跑出了盒子是怎么回事?

微信截图_20230428210843.png (249.48 KB, 下载次数 Times of downloads: 11)

微信截图_20230428210843.png

graphite_1073K_Xe_0.01_01.inp

49.48 KB, 下载次数 Times of downloads: 15

197

帖子

1

威望

1108

eV
积分
1325

Level 4 (黑子)

2#
发表于 Post on 2023-4-29 15:21:25 | 只看该作者 Only view this author
请问一下这个。inp参数默认设置的是BOMD吗

10

帖子

0

威望

139

eV
积分
149

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2023-4-29 16:07:45 | 只看该作者 Only view this author
Weldingspock 发表于 2023-4-29 15:21
请问一下这个。inp参数默认设置的是BOMD吗

按照http://bbs.keinsci.com/thread-21981-1-1.html中的解释,基于xTB就是BOMD。请问你知道为啥这个原子会出去吗?这跟MM模拟中的周期性是不是不一样,还是我的设置有问题?

1061

帖子

0

威望

3271

eV
积分
4332

Level 6 (一方通行)

4#
发表于 Post on 2023-4-29 17:08:43 | 只看该作者 Only view this author
we12306 发表于 2023-4-29 16:07
按照http://bbs.keinsci.com/thread-21981-1-1.html中的解释,基于xTB就是BOMD。请问你知道为啥这个原子 ...

用 VMD 打开 .xyz 轨迹。

在 VMD 的主功能界面的菜单栏中找到 Extensions ——> Tk Console,然后输入
  1. pbc set {a, b, c, alpha, beta, gamma} -all   ##号后面不用输入,意义为 a, b, c 是盒子(晶格)长度,alpha, beta, gamma 是角度
  2. pbc box -all
  3. pbc wrap -all
复制代码

10

帖子

0

威望

139

eV
积分
149

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2023-4-29 17:17:56 | 只看该作者 Only view this author
乐平 发表于 2023-4-29 17:08
用 VMD 打开 .xyz 轨迹。

在 VMD 的主功能界面的菜单栏中找到 Extensions ——> Tk Console,然后输入 ...

非常感谢!!!我还要一个问题是请问如果不这样处理,直接将原始的XYZ坐标和受力输入到神经网络中是不是训练的会有问题?

1061

帖子

0

威望

3271

eV
积分
4332

Level 6 (一方通行)

6#
发表于 Post on 2023-4-29 17:27:56 | 只看该作者 Only view this author
we12306 发表于 2023-4-29 17:17
非常感谢!!!我还要一个问题是请问如果不这样处理,直接将原始的XYZ坐标和受力输入到神经网络中是不是 ...

轨迹本身就已经是满足周期性边界条件的,上面的设置只是将原子“显示”在盒子里而已。
不影响你后续的计算。

10

帖子

0

威望

139

eV
积分
149

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2023-4-29 17:56:24 | 只看该作者 Only view this author
乐平 发表于 2023-4-29 17:27
轨迹本身就已经是满足周期性边界条件的,上面的设置只是将原子“显示”在盒子里而已。
不影响你后续的计 ...

好的,非常感谢您的回复!!

本版积分规则 Credits rule

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

GMT+8, 2024-11-26 17:32 , Processed in 1.569048 second(s), 24 queries , Gzip On.

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