计算化学公社

标题: 求助：非限制性开窍层如何选取活性空间 [打印本页]

作者
Author: Lapheal 时间: 2024-5-2 19:23
标题: 求助：非限制性开窍层如何选取活性空间
本帖最后由 Lapheal 于 2024-5-2 19:34 编辑

使用UHF计算得到了UNO轨道，如图所示，α和β电子有各自的轨道。
如果想选取轨道放入活性空间，如选取25，26号轨道放入活性空间，请问需要将120，121号轨道交换放入活性空间吗？

作者
Author: hebrewsnabla 时间: 2024-5-2 19:38
UNO不分alpha beta，只有一套。这个图片显示的大概不是UNO。

作者
Author: zjxitcc 时间: 2024-5-2 21:19
你这不是UNO。如果你有UHF的fch文件，可以用MOKIT的uno()函数产生UNO的fch文件，只有1列轨道，并且UNO没有轨道能量的概念，有轨道占据数。

进一步地，可以不产生UNO，而使用automr自动完成多组态/多参考态计算，例如读取一个合理的UHF fch文件，gjf文件如下

%mem=4GB
%nprocshared=2
#p CASSCF/cc-pVDZ
mokit{ist=1,readuhf='h2o_uhf.fch'}

复制代码

automr会先产生UNO，然后产生局域配对UNO（作为GVB初始轨道），进行GVB计算获得GVB轨道（作为CAS初始轨道），进行CASSCF计算。最后得到一个_CASSCF_NO.fch文件，内含CASSCF自然轨道及其轨道占据数，可用GaussView/Multiwfn打开、可视化。以水分子为例，会自动算出CASSCF(4,4)的结果。

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