求助利用AICD分析开壳层pi体系

Exciplex

请问sob老师和各位大神，在使用AICD (http://sobereva.com/294) 分析开壳体系的时候，输入前线的pi轨道的编号之后，电流密度总比选择闭壳的减少一半。有的时候只是看上去瘦，但有的时候电流密度只有一半，严重影响分析。请问这是因为什么，要如何解决？

sobereva

确保轨道序号输对了，alpha和beta都得输入
而且先拿个闭壳层体系对比基于RKS和UKS的结果，检验是否相符。
如果确信不是序号的问题，可能本来结果就那样；如果怀疑是程序bug，发邮件问作者

Exciplex

sobereva 发表于 2023-2-19 00:56
确保轨道序号输对了，alpha和beta都得输入
而且先拿个闭壳层体系对比基于RKS和UKS的结果，检验是否相符。
...

这边alpha和beta输入的格式是什么呢？好像输入xxxa,xxxb或者用multiwfn输出的轨道编号都没法正确进行计算。。

sobereva

Exciplex 发表于 2023-2-19 15:27
这边alpha和beta输入的格式是什么呢？好像输入xxxa,xxxb或者用multiwfn输出的轨道编号都没法正确进行计算 ...

我的博文里难道没明确说么？

对于非限制性开壳层（U）形式的计算，IOp手册里没写清楚怎么指定要考虑的轨道。我这里根据我读源代码以及实际测试做出的理解介绍一下。比如计算三重态水分子，原本一共有6个占据的alpha轨道和4个占据的beta轨道，在你设置要考虑的轨道时，记住beta占据轨道的序号是排在alpha占据轨道后面的，而不用管alpha和beta的空轨道。比如，如果你要让原先的那10个占据轨道都被考虑，输入文件末尾就写1-10。如果你只想考虑alpha占据轨道当中后三个、beta占据轨道中第2个，就应该写4-6 8。之所以是8，是因为6+2=8，即原本的alpha占据轨道数目加上要考虑的这个beta的序号。空轨道不用管，比如你把11、12等序号超过原先占据轨道总数的轨道也纳入的话，不会影响结果。