计算化学公社

标题: CP2K中，非正交模型周期性边界问题 [打印本页]

作者
Author: 桂の猿 时间: 2024-3-26 21:43
标题: CP2K中，非正交模型周期性边界问题
各位老师好

我最经在做晶体的吸附计算，计算使用cp2k，晶体模型如图，是六方晶系，我想使用周期性边界条件，我按过去的正交模型的Cell给出的边界条件&CELL
   A    9.33600000    0.00000000    0.00000000
   B -4.66800000    8.08521300    0.00000000
   C    0.00000000    0.00000000 14.94600000
   PERIODIC XYZ
但是，这样的话，周期是不是不对？是否会默认成xyz方向的周期性，而不是晶体侧面的方向？请问应该如何修正或者有类似例子吗？谢谢

作者
Author: sobereva 时间: 2024-3-26 22:20
明显当前的晶胞信息不对。想想按照当前晶胞平移复制成超胞变成什么样子自然就知道了

PERIODIC XYZ指的是三维周期性，不是指的各个晶格矢非得平行于各个笛卡尔轴

作者
Author: 乐平 时间: 2024-3-26 22:54
本帖最后由乐平于 2024-3-26 23:12 编辑

其实可以偷懒

查一下 CP2K 手册关于 SUBSYS 小节的里的 CELL 和 TOPOLOGY，你会发现，CP2K 可以直接读 cif 文件里的晶胞参数，以及晶体坐标信息。你不需要从六方转换到正交。

CELL
https://manual.cp2k.org/trunk/CP ... LL.CELL_FILE_FORMAT

TOPOLOGY
https://manual.cp2k.org/trunk/CP ... UBSYS/TOPOLOGY.html

下面是例子：

&CELL
PERIODIC XYZ
CELL_FILE_FORMAT CIF
CELL_FILE_NAME YOUR_FILE_NAME
&END CELL
&TOPOLOGY
COORD_FILE_FORMAT CIF
COORD_FILE_NAME YOUR_FILE_NAME
&END TOPOLOGY

复制代码

也就是说，你只需要把你的晶体结构的 .cif 文件放在 CP2K 输入文件（.inp 文件）的同一个文件夹内，把晶体结构名替换掉上述代码里的 YOUR_FILE_NAME 就可以了。CP2K 就能自动读取晶胞参数和晶体结构信息。

这样做还有另外的好处，那就是，对于一系列计算可以保持计算参数完全相同，只用改结构文件名。可以很大程度避免粗心造成计算参数设置偏差（虽然 Multiwfn 设置 CP2K 输入文件很好用，但是 CP2K 的参数实在有些多，如果不关闭 Multiwfn 一次输入产生很多输入文件可以保证参数一致，但是关闭 Multiwfn 后，上次设置的参数会清零，万一不小心漏掉了某个设置，检查起来有点小头疼）。

作者
Author: 桂の猿 时间: 2024-3-27 21:35

sobereva 发表于 2024-3-26 22:20
明显当前的晶胞信息不对。想想按照当前晶胞平移复制成超胞变成什么样子自然就知道了

PERIODIC XYZ指的是 ...

谢谢sob老师，我是利用您开发的Multiwfn先将GaN的cif建立了一个4*4*3的超胞，如果输出xyz文件，然后用ovito看，就是如图所示，如果我输出cif文件，则是正确的。
如果我希望利用 PERIODIC XYZ ，是不是应该将现在的文件正交化才可以？谢谢

作者
Author: 桂の猿 时间: 2024-3-27 21:39

乐平发表于 2024-3-26 22:54
其实可以偷懒

查一下 CP2K 手册关于 SUBSYS 小节的里的 CELL 和 TOPOLOGY，你会发现，CP2K 可以直接读 c ...

谢谢乐平老师，我是利用sob开发的Multiwfn先将GaN的cif建立了一个4*4*3的超胞，如果输出xyz文件，然后用ovito看，就是如图所示，如果我输出cif文件，则是正确的。

如果我用您的方法，使用超胞cif文件，那么这里的周期性边界是不是也是有问题的？谢谢

作者
Author: sobereva 时间: 2024-3-28 00:16

桂の猿发表于 2024-3-27 21:35
谢谢sob老师，我是利用您开发的Multiwfn先将GaN的cif建立了一个4*4*3的超胞，如果输出xyz文件，然后用ovi ...

我从来不用OVITO
xyz文件没盒子信息，任何可视化程序都不可能还原出原本的盒子
跟正交化毫无关系，前面我都强调了

PERIODIC XYZ指的是三维周期性，不是指的各个晶格矢非得平行于各个笛卡尔轴

欢迎光临计算化学公社 (http://bbs.keinsci.com/)