计算化学公社

标题: 求助gaussian十六极矩以上如何输出 [打印本页]

作者
Author: cuihongqiang 时间: 2022-7-6 21:40
标题: 求助gaussian十六极矩以上如何输出
使用高斯09能够输出多极矩最高阶到16极矩，请问如何可以输出16极矩以上的多极矩，比如32极矩和64极矩等

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

作者
Author: sobereva 时间: 2022-7-7 03:32
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，仔细看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节），将输出的信息扩展到十六极矩以上。

作者
Author: cuihongqiang 时间: 2022-7-7 09:31

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

谢谢

作者
Author: cuihongqiang 时间: 2022-7-7 09:32

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

谢谢

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