CP2K计算偶极矩的问题

1739140886 · 发表于 Post on 2022-5-3 15:57:14

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

Eudaimonia · 发表于 Post on 2022-5-3 16:26:14

纯有机体系不含周期性可能ORCA会更好一点？不确定

1739140886 · 发表于 Post on 2022-5-3 16:36:27

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

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

Eudaimonia · 发表于 Post on 2022-5-3 19:33:59

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
即可

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

sobereva · 发表于 Post on 2022-5-3 19:34:19

没必要尝试CP2K
对此体系ORCA是首选

1739140886 · 发表于 Post on 2022-5-3 20:11:27

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

感谢感谢，我去试试！！

1739140886 · 发表于 Post on 2022-5-3 20:12:18

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

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

qczgzly · 发表于 Post on 2023-12-17 16:24:16

试试在&DFT的&PRINT内加入
&MOMENTS
MAX_MOMENT 1 #1是偶极矩，2是四级矩，依此类推，周期性体系最大只能设置1
PERIODIC Y #如果不是周期性体系就改成N
&END MOMENTS

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

[CP2K] CP2K计算偶极矩的问题