计算化学公社

标题: Multiwfn_3.8_dev 产生 CP2K 输入文件时，设置周期性之后，晶胞参数发生变动？ [打印本页]

作者
Author: lightman 时间: 2023-3-12 23:25
标题: Multiwfn_3.8_dev 产生 CP2K 输入文件时，设置周期性之后，晶胞参数发生变动？
sob老师您好！我在使用 Multiwfn_3.8_dev_bin_Win64 产生 CP2K 输入文件时发现如下问题：

尝试学习一个您之前列举的CP2K例子，即在NaCl晶体表面吸附一个一氧化碳CO分子的结构优化的任务时，在选择周期性XY方向（-7 Set direction(s) of applying periodic boundary condition）后晶胞参数Z 发生了变动（由25变为了25.38778436）。下面附了GaussView产生的结构模型文件（NaCl_CO.gjf）以及Multiwfn_3.8_dev_bin_Win64 产生的CP2K输入文件：在周期性XY方向重复的CP2K输入文件（NaCl_CO_XY.inp）和在周期性XYZ方向重复的CP2K输入文件（NaCl_CO_XYZ.inp）。

sob老师您看看是怎么回事啊？谢谢啦

(, 下载次数 Times of downloads: 10) (, 下载次数 Times of downloads: 12)

作者
Author: sobereva 时间: 2023-3-12 23:34
Multiwfn对二维周期性情况会根据当前用的Poisson solver自动确定合适的真空区厚度，避免浪费又避免不够

(, 下载次数 Times of downloads: 9)

仅对于三维周期性情况才会直接用输入文件里记录的盒子信息

作者
Author: lightman 时间: 2023-3-12 23:43

sobereva 发表于 2023-3-12 23:34
Multiwfn对二维周期性情况会根据当前用的Poisson solver自动确定合适的真空区厚度，避免浪费又避免不够

...

原来是这样啊，我明白了，谢谢sob老师

作者
Author: pikachuupup 时间: 2025-8-29 17:47

sobereva 发表于 2023-3-12 23:34
Multiwfn对二维周期性情况会根据当前用的Poisson solver自动确定合适的真空区厚度，避免浪费又避免不够

...

请问老师，如果想让c方向的长度为一固定值，方便与之前算的结构进行对照，那么我可以直接在Multiwfn输出的inp上手动更改c的值吗，手动改完c后的POISSON字段需要更改吗？
&POISSON
PERIODIC XY
PSOLVER MT
&END POISSON

作者
Author: sobereva 时间: 2025-8-29 18:01

pikachuupup 发表于 2025-8-29 17:47
请问老师，如果想让c方向的长度为一固定值，方便与之前算的结构进行对照，那么我可以直接在Multiwfn输出 ...

可以

不改。注意c可以相对于Multiwfn默认的改大了，但如果改得比默认的更小，会不符合MT的要求

作者
Author: pikachuupup 时间: 2025-8-30 16:04

sobereva 发表于 2025-8-29 18:01
可以

不改。注意c可以相对于Multiwfn默认的改大了，但如果改得比默认的更小，会不符合MT的要求

明白了，谢谢老师~

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