计算化学公社

标题: CP2K计算偶极矩的问题 [打印本页]

作者
Author: 1739140886 时间: 2022-5-3 15:57
标题: CP2K计算偶极矩的问题
最近在计算只包含前两周期元素的纯有机分子体系（原子数500左右），想对其进行偶极矩计算，之前一直在用Gaussian09，感觉计算效率不高，想转用CP2K试试（基组def2-SVPD就可以），想了解一下如何输出偶极矩信息呢？谢谢大家解答！！

作者
Author: Eudaimonia 时间: 2022-5-3 16:26
纯有机体系不含周期性可能ORCA会更好一点？不确定

作者
Author: 1739140886 时间: 2022-5-3 16:36

Eudaimonia 发表于 2022-5-3 16:26
纯有机体系不含周期性可能ORCA会更好一点？不确定

目前也是在尝试，组里比较建议用CP2K，就想跑一下试试，就不太懂得如何输出偶极矩信息

作者
Author: Eudaimonia 时间: 2022-5-3 19:33

1739140886 发表于 2022-5-3 16:36
目前也是在尝试，组里比较建议用CP2K，就想跑一下试试，就不太懂得如何输出偶极矩信息

我在之前的一次计算中，本来目的是计算电荷布居，把Multiwfn支持的所有类型都print出来了，里面输出了一个后缀为.voronoi的文件，里面数据如下所示
# Atom    Z             Radius       Position(X)       Position(Y)       Position(Z)    Voronoi_Volume             Z(eff)             Charge          Dipole(X)          Dipole(Y)          Dipole(Z)    ChargeCenter(X)    ChargeCenter(Y)    ChargeCenter(Z)    Quadrupole(XX)    Quadrupole(XY)    Quadrupole(XZ)    Quadrupole(YX)    Quadrupole(YY)    Quadrupole(YZ)    Quadrupole(ZX)    Quadrupole(ZY)    Quadrupole(ZZ
可能是你需要的偶极矩
在&DFT的&PRINT内加入
   &VORONOI
      VORONOI_RADII Covalent
   &END VORONOI
即可

你可以尝试一下看看是不是你需要的偶极矩

作者
Author: sobereva 时间: 2022-5-3 19:34
没必要尝试CP2K
对此体系ORCA是首选

作者
Author: 1739140886 时间: 2022-5-3 20:11

Eudaimonia 发表于 2022-5-3 19:33
我在之前的一次计算中，本来目的是计算电荷布居，把Multiwfn支持的所有类型都print出来了，里面输出了一 ...

感谢感谢，我去试试！！

作者
Author: 1739140886 时间: 2022-5-3 20:12

sobereva 发表于 2022-5-3 19:34
没必要尝试CP2K
对此体系ORCA是首选

感谢sob老师，那我也去了解下ORCA！！

作者
Author: qczgzly 时间: 2023-12-17 16:24
试试在&DFT的&PRINT内加入
&MOMENTS
MAX_MOMENT 1 #1是偶极矩，2是四级矩，依此类推，周期性体系最大只能设置1
PERIODIC Y #如果不是周期性体系就改成N
&END MOMENTS

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