怎么通过自然轨道占据数知道哪些轨道是价层电子轨道？

Xx_DerEnTh

zjxitcc 发表于 2025-5-11 15:16
可以让MOKIT不做UHF计算。可以使用UCCSD自然轨道（注意必须是空间自然轨道，只有1列轨道；而不是自旋自然 ...

老师您好，想再请教一下Gaussian如何产生空间自然轨道？我尝试了用UCCSD+density+naturalorbitals两步来得到自然轨道，发现得到的轨道并不符合要求。看了sob老师的博文http://sobereva.com/403也尝试了先计算UCCSD+density，然后生成fchk文件传入Multiwfn，用功能200的子功能16，输入CC和1，得到new.mwfn文件，再把new.mwfn文件转化成了fch，得到的轨道好像也不符合要求，搜索Alpha O下面的全是0。

zjxitcc

Xx_DerEnTh 发表于 2025-5-11 18:50
老师您好，想再请教一下Gaussian如何产生空间自然轨道？我尝试了用UCCSD+density+naturalorbitals两步来 ...

Gaussian本身不支持产生UCCSD空间自然轨道。不过，可以用MOKIT调用Gaussian做UCCSD计算，计算完成后当前目录下自动就有xx_CC_NO.fch文件，内含UCCSD空间自然轨道（这是MOKIT读取Gaussian的UCCSD弛豫密度然后算的），以单重态O2为例，gjf文件为

1. %mem=20GB
   
2. %nprocshared=10
   
3. #p CCSD/def2TZVP
   
4. 
   
5. mokit{NoRI,CC_prog=Gaussian,Relaxed_DM,NO}
   
6. 
   
7. 0 1
   
8. O   0.0   0.0   0.0
   
9. O   0.0   0.0   1.1616
   

复制代码

Gaussian目前不支持UCCSD非弛豫密度，所以用关键词Relaxed_DM指定了计算弛豫密度。提交任务，即运行

1. autosr O2.gjf >O2.out 2>&1

复制代码

大概1 min结束，当前目录下有O2_uhf_CC_NO.fch文件，可以用GaussView/Multiwfn打开可视化，发现HONO和LUNO占据数均为1.01，是完美双自由基、开壳层单重态。