栗悟饭与龟波功 发表于 2024-5-4 15:12
社长您好！抱歉再次打扰，昨天电脑不在身边，不是很确定，今天仔细查看hyper_polar.f90后发现似乎beta || ...

读取的数值的有效位数已经足够多了，没有改进的余地，读更多位数也不会对研究有任何实际的影响

sobereva 发表于 2024-5-3 23:15
这部分功能的目的是处理（超）极化率的数据，偶极矩仅仅是顺带输出。如果就是要偶极矩的数据，直接从输出文 ...

社长您好！抱歉再次打扰，昨天电脑不在身边，不是很确定，今天仔细查看hyper_polar.f90后发现似乎beta ||的值也是根据前者读取的dipole计算出来的，尽管在dipole如此小的情况谈这些量都没什么意义了。但个人感觉有必要向您汇报一下这个情况，如果可行的话，能否用在计算beta ||其他模块的dipole代替？毕竟这个值确实有时候得和实验对比

谢谢社长！看来我体系对称性太高了，导致与dipole相关的β||也没有可信度了，我是在做不同参数泛函cpdft下结果与ccsd的比较时遇到了困难，β本身对的不错但一旦牵扯上dipole就完蛋了。多亏社长一语点醒梦中人，应该选择具有相对较强dipole的体系。

这部分功能的目的是处理（超）极化率的数据，偶极矩仅仅是顺带输出。如果就是要偶极矩的数据，直接从输出文件读，或者用Multiwfn的主功能300的主功能5直接基于波函数信息计算就行了。Multiwfn当前功能保留的小数位位数足矣满足一般需要，而且偶极矩再怎么小，a.u.为单位给出六位小数也足矣了，而且根本没有什么计算级别能把偶极矩计算到小数位第六位的精度。