求助gaussian十六极矩以上如何输出

cuihongqiang · 发表于 Post on 2022-7-6 21:40:27

使用高斯09能够输出多极矩最高阶到16极矩，请问如何可以输出16极矩以上的多极矩，比如32极矩和64极矩等

wzkchem5 · 发表于 Post on 2022-7-6 21:53:22

估计是不行的。不过这么高阶的多极矩一般也没太大用吧？很多场合这种超高阶多极矩应该都可以用原子电荷或静电势代替

sobereva · 发表于 Post on 2022-7-7 03:32:00

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“gaussian09 十六极矩以上输出 ”改了，以后务必注意

Gaussian输出不了，改改Multiwfn代码可以容易地实现，稍微有举一反三的能力就能做到。具体来说是修改integral.f90里的subroutine domultipoleint和subroutine genMultipoleprim，将能够计算的多极矩积分矩阵扩展到更高阶（参考里面原本代码依葫芦画瓢即可，很容易），然后修改otherfunc3.f90里的subroutine calc_multipole（对应的功能的介绍看Multiwfn手册3.300.5节），将输出的信息扩展到十六极矩以上。

cuihongqiang · 发表于 Post on 2022-7-7 09:31:44

wzkchem5 发表于 2022-7-6 21:53
估计是不行的。不过这么高阶的多极矩一般也没太大用吧？很多场合这种超高阶多极矩应该都可以用原子电荷或静 ...

谢谢

cuihongqiang · 发表于 Post on 2022-7-7 09:32:16

sobereva 发表于 2022-7-7 03:32
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此 ...

谢谢

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

[Gaussian/gview] 求助gaussian十六极矩以上如何输出