关于Multiwfn的一些疑惑与建议

栗悟饭与龟波功

社长，各位老师同学你们好！

首先非常感谢社长无私地向我们分享了Multiwfn这样一款非常好的软件，解决了我工作中的许多问题。

但是，在使用Multiwfn的hyperpolarization的功能时（Version 3.8(dev), release date: 2024-Jan-16），我发现有一点地方可能需要一些改进。我发现在hyper_polar.f90中关于dipole moment的读取是定位了“Dipole moment (field-independent basis, Debye)”这一行的内容，正常情况下（假设standard orientation与input orientation一致）应该是与“Electric dipole moment (input orientation):”的输出是一致的，但实际上后者能显示的精度是更高的。这在遇到一些具有较小dipole的体系的情况时尤为重要，比如附件图片所示，由于“Dipole moment (field-independent basis, Debye)”那一行的精度过低导致输出的结果被截断了。



这是Multiwfn的输出，也就是“Dipole moment (field-independent basis, Debye)”除以au2debye，可以看到在这种情况下其实很大一部分的值被截断了。




这是我个人使用中遇到的一点小问题，向社长反映一下，希望Multiwfn越来越好！

sobereva

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

栗悟饭与龟波功

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

栗悟饭与龟波功

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

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

sobereva

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

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