计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First Principle) › 关于DeePMD-Kit支持CP2K的数据格式的一个疑问(已解决)

关于DeePMD-Kit支持CP2K的数据格式的一个疑问(已解决)

查看数: 4171 | 评论数: 13 | 收藏 Add to favorites 5

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

joeson

发布时间: 2022-7-23 15:01

正文摘要:

本帖最后由 joeson 于 2022-7-24 21:32 编辑 DeePMD-Kit数据转化支持的格式有一个是CP2K的aimd_output，请问这个格式如何得到？已解决：原来就是指需要AIMD的.xyz与包含受力的.log文件。在github上有人回 ...

回复 Reply

Uus/pMeC6H4-/キ 发表于 Post on 4 day ago

Weldingspock 发表于 2024-11-6 09:38
没啥区别，只不过frc-1.xyz里面的格式可能更规整，更好文本处理

今天试了一下，原来设&FORCE_EVAL/&PRINT/&FORCES输出的受力信息并非在我以为的另一个.log文件，而是就在CP2K的主输出文件.out里。在2024.3版上，这段信息一共6列，依次是原子序号（1开始）、原子种类（1开始，同元素原子同种类）、元素名以及X、Y、Z方向受力分量，但数据位数比-frc-1.xyz里的少。

还有另一个小问题，用于DeePMD-kit训练的原子坐标数据对PBC是怎么考虑的？CP2K虽然输入有定义PBC以及盒子的A、B、C三个矢量，但输出的-pos-1.xyz并未把原子卷到盒子内而是可能漂移出去。有没有必要在训练前先自行wrap处理坐标，让原子仅出现在box.raw定义的盒子内呢？

Weldingspock 发表于 Post on 2024-11-6 09:38:02

Uus/pMeC6H4-/キ发表于 2024-11-3 14:05
有两点不大明白：

1.github里面提到要设置&FORCE_EVAL/&PRINT/&FORCES这个字段为ON来输出受力信息到.l ...

没啥区别，只不过frc-1.xyz里面的格式可能更规整，更好文本处理

Uus/pMeC6H4-/キ 发表于 Post on 2024-11-3 14:05:44

已解决：原来就是指需要AIMD的.xyz与包含受力的.log文件。在github上有人回复是用这两个文件转化成一个dpmd的npy文件。https://github.com/deepmodeling/dpdata/issues/103
补充说明：CP2K跑MD时，PRINT_LEVEL要设置为MEDIUM。

有两点不大明白：

1.github里面提到要设置&FORCE_EVAL/&PRINT/&FORCES这个字段为ON来输出受力信息到.log文件，而看手册说明这个字段默认输出水平是HIGH。显式设置这个字段为ON的话，应该无论&GLOABL/&PRINT_LEVEL设置了什么输出水平都能产生输出的.log文件吧，那为什么补充说明又说要设置&GLOABL/&PRINT_LEVEL为MEDIUM呢，是不是还有什么重要信息是平时常用的输出水平LOW不会输出但DeePMD-kit需要用的？

2.平时用CP2K跑分子动力学的时候用&MOTION/&PRINT/&FORCES也能按一定步数间隔输出类似坐标.xyz格式的、包含各帧各原子受力分量信息的<proj>-frc-1.xyz文件，那&FORCE_EVAL/&PRINT/&FORCES输出受力信息的.log文件和<proj>-frc-1.xyz相比，有什么对DeePMD-kit很关键的区别吗？

whutzty 发表于 Post on 2024-10-11 17:01:05