计算化学公社

标题: 求助利用AICD分析开壳层pi体系 [打印本页]
作者
Author: Exciplex    时间: 2023-2-18 23:53
标题: 求助利用AICD分析开壳层pi体系
请问sob老师和各位大神，在使用AICD (http://sobereva.com/294) 分析开壳体系的时候，输入前线的pi轨道的编号之后，电流密度总比选择闭壳的减少一半。有的时候只是看上去瘦，但有的时候电流密度只有一半，严重影响分析。请问这是因为什么，要如何解决？
作者
Author: sobereva    时间: 2023-2-19 00:56
确保轨道序号输对了，alpha和beta都得输入
而且先拿个闭壳层体系对比基于RKS和UKS的结果，检验是否相符。
如果确信不是序号的问题，可能本来结果就那样；如果怀疑是程序bug，发邮件问作者
作者
Author: Exciplex    时间: 2023-2-19 15:27
sobereva 发表于 2023-2-19 00:56
确保轨道序号输对了，alpha和beta都得输入
而且先拿个闭壳层体系对比基于RKS和UKS的结果，检验是否相符。
...

这边alpha和beta输入的格式是什么呢？好像输入xxxa,xxxb或者用multiwfn输出的轨道编号都没法正确进行计算。。
作者
Author: sobereva    时间: 2023-2-20 02:25
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等序号超过原先占据轨道总数的轨道也纳入的话，不会影响结果。





欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3