使用Multiwfn非常便利地创建CP2K程序的输入文件

sobereva · 发表于 Post on 2022-3-8 08:59:11

今日更新的Multiwfn的CP2K输入文件创建功能新支持了产生BSSE任务的输入文件
任务类型里选BSSE，按提示输入即可

captain · 发表于 Post on 2022-4-28 16:30:19

本帖最后由 captain 于 2022-4-28 16:32 编辑

请问大神，我借助Multiwfn产生的CP2K输入文件如下：
&CELL
ABC 26.000 26.000 41.000
PERIODIC NONE #Direction of applied PBC (geometry aspect)
&END CELL
想进行周期性计算，而这里的cell比我的体系的尺寸大一些（15*15*30），需要再修改一下ABC后面的数值吗？
另外，PERIODIC 后面的NONE(以及下面的 &POISSON PERIODIC NONE #Direction(s) of PBC for calculating electrostatics)应该也相应地根据情况修改为如XY、XYZ等吧？

sobereva · 发表于 Post on 2022-4-29 08:38:32

captain 发表于 2022-4-28 16:30
请问大神，我借助Multiwfn产生的CP2K输入文件如下：
&CELL
ABC 26.000 26.000 41.000

如果你载入的输入文件里本身就有盒子信息，Multiwfn产生的CP2K输入文件的盒子信息将和输入文件一致，并且默认的PERIODIC就是实际的周期性
如果输入文件没有盒子信息，程序会认为这是个孤立体系，自动往各个方向扩展一个距离来定义盒子，且默认PERIODIC NONE

captain · 发表于 Post on 2022-4-29 08:44:22

sobereva 发表于 2022-4-29 08:38
如果你载入的输入文件里本身就有盒子信息，Multiwfn产生的CP2K输入文件的盒子信息将和输入文件一致，并且 ...

明白了，多谢大神指点！

西麦西麦 · 发表于 Post on 2022-5-8 17:18:13

请问sob老师，用multiwfn生成cp2k输入文件的时候，比如CuOx，怎么调校U值和enforce_occupation的参数，另外现在的cp2k的版本k点和+U的功能是否仍然冲突？有没有调节的方法？

annaqz · 发表于 Post on 2022-5-16 09:23:53

cp2k优化晶胞参数可以固定 xy 方向只优化 z 方向吗？需要在输入文件中怎么设置呢？

sobereva · 发表于 Post on 2022-5-17 02:10:51

西麦西麦发表于 2022-5-8 17:18
请问sob老师，用multiwfn生成cp2k输入文件的时候，比如CuOx，怎么调校U值和enforce_occupation的参数，另外 ...

只能自己手动输入，调校不是Multiwfn的事

仍然冲突

sobereva · 发表于 Post on 2022-5-17 02:12:38

annaqz 发表于 2022-5-16 09:23
cp2k优化晶胞参数可以固定 xy 方向只优化 z 方向吗？需要在输入文件中怎么设置呢？

用Multiwfn产生结构优化且开启原子位置约束的任务的输入文件，然后手动修改输入文件里COMPONENTS_TO_FIX后面的内容

annaqz · 发表于 Post on 2022-5-17 09:30:53

sobereva 发表于 2022-5-17 02:12
用Multiwfn产生结构优化且开启原子位置约束的任务的输入文件，然后手动修改输入文件里COMPONENTS_TO_FIX ...

好的，那我试试。谢谢社长！

sobereva · 发表于 Post on 2022-5-18 18:48:41

今日更新的Multiwfn保存CP2K文件时，如果用的是MOLOPT系列基组，或者DZVP-GTH等老的结合GTH的赝势基组，在得到的CP2K输入文件中会看到基组名后面已经带了诸如-q5这样的价电子数。这样便于用户使用Multiwfn基于CP2K的molden文件进行分析之前在molden文件里插入[Nval]字段设置元素的价电子数（此字段里的内容应当对应-q后面的数字）

wth1219 · 发表于 Post on 2022-7-14 12:56:36

卢老师好，最新版本CP2K 2022.1删除了&XC_FUNCTIONAL中的&LibXC，将交换相关泛函作为独立的模块，可能需要您稍微调整一下Multiwfn中涉及到的输入文件模板了。
（吐槽：CP2K此前一向重视向下兼容的，怎么感觉也开始和ORCA学了……）

例如SCAN泛函原先的输入格式为：
   &XC_FUNCTIONAL
      &LIBXC
      FUNCTIONAL MGGA_X_SCAN
      &END LIBXC
      &LIBXC
      FUNCTIONAL MGGA_C_SCAN
      &END LIBXC
   &END XC_FUNCTIONAL

使用2022.1版本的输入格式为：
   &XC_FUNCTIONAL
      &MGGA_X_SCAN
      &END MGGA_X_SCAN
      &MGGA_C_SCAN
      &END MGGA_C_SCAN
   &END XC_FUNCTIONAL

sobereva · 发表于 Post on 2022-7-15 10:15:26

wth1219 发表于 2022-7-14 12:56
卢老师好，最新版本CP2K 2022.1删除了&XC_FUNCTIONAL中的&LibXC，将交换相关泛函作为独立的模块，可能需要 ...

这种交替的时候很麻烦，改成兼容2022.1，1000%又一堆用<=9.1的人问怎么输入文件跑不了

wth1219 · 发表于 Post on 2022-7-15 11:06:36

sobereva 发表于 2022-7-15 10:15
这种交替的时候很麻烦，改成兼容2022.1，1000%又一堆用

方才查了一下CP2K的Reference Manual，最近三个版本（8.2、9.1、2022.1）的&XC_FUNCTIONAL部分均补入了大量先前必须依靠&LibXC定义的泛函（包括SCAN、B97M-V。），上面那种新的定义方式是适用的。因此最简单粗暴的解决方案是直接说明部分输入文件仅适用于8.2及之后版本。正常应该像您生成ORCA输入文件那样做一个版本选择，但是正如您所说估计还是会有人直接忽略了。

sobereva · 发表于 Post on 2022-7-16 10:12:13

wth1219 发表于 2022-7-15 11:06
方才查了一下CP2K的Reference Manual，最近三个版本（8.2、9.1、2022.1）的&XC_FUNCTIONAL部分均补入了大 ...

我已更新了官网上Multiwfn最新版本，兼容了此写法。
<=8.2版本的用户被我放弃了，他们早晚会灭绝的

wudiazhu · 发表于 Post on 2022-7-21 16:29:03

今天发现用Multiwfn生成cp2k的输入文件时把我结构文件里的Pb原子识别成了P原子，但我的Cl原子是正常的没有识别成C原子，所以我在想可能是个小bug，我用的Multiwfn 3.8 win版，任务是跑动力学。另外，初学cp2k，想问几个小问题：
1. 用cp2k跑AIMD是否需要保证体系为电中性？
2. AIMD怎么对重原子用赝势基组，对非重原子用全电子基组？我在Multiwfn设置基组那里选了DZVP-MOLOPT-SR-GTH，全部设置完跑任务出现报错说“The requested basis set <CFIT3> for element <PB> was not found in the basis set files <BASIS_MOLOPT><BASIS_ADMM>”，但我看cp2k-9.1/data/BASIS_MOLOPT里面是有Pb元素的赝势定义的，不知道我哪里出了问题？
3. 我没太搞懂cp2k输入文件里的基组定义部分为什么每个元素有2个BASIS SET？其中的一个BASIS SET能否删去？POTENTIAL行是指什么？
问的可能有点基础，先谢谢各位老师了！

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 使用Multiwfn非常便利地创建CP2K程序的输入文件

评分 Rate

评分 Rate