运行NMR=CSGT相关问题

Reolli

按照基于Gaussian的NMR=CSGT任务得到各个轨道对NICS贡献的方法 - 思想家公社的门口：量子化学·分子模拟·二次元教程进行pi轨道对NICS(0)ZZ的贡献，在指定pi轨道的时候报错，提示
Range of integers=   655 to    666 is invalid.  Limit is     1 to   324.
Error in AddMOL.
Error termination via Lnk1e in /share/home/u22310164/install_app//g16/l1002.exe at Tue May 20 21:30:55 2025.
计算的体系是单重态开壳层，pi轨道经查看为61-72,79-81,655-666,673-675，想知道为啥提示限度为324，下面是multiwfn查看轨道的数据。
Range of alpha orbitals:    1 -  594      Range of beta orbitals:  595 - 1188
Note: Orbital    81 is alpha-HOMO, energy:   -0.268755 a.u.   -7.313184 eV
       Orbital   675 is beta-HOMO, energy:    -0.268755 a.u.   -7.313185 eV
       Orbital    82 is alpha-LUMO, energy:   -0.064333 a.u.   -1.750590 eV
       Orbital   676 is beta-LUMO, energy:    -0.064333 a.u.   -1.750585 eV
       HOMO-LUMO gap of alpha orbitals:    0.204422 a.u.    5.562594 eV
       HOMO-LUMO gap of beta orbitals:     0.204422 a.u.    5.562601 eV

sobereva

上传输入输出文件

Reolli

sobereva 发表于 2025-5-21 05:28
上传输入输出文件

老师您好，这是我的输入输出文件。

sobereva

当前情况下，Gaussian里的轨道序号顺序是：[alpha占据轨道][beta占据轨道]
alpha和beta都是81个电子，因此占据轨道是1-162，你指定的序号不应超过这个范围

http://sobereva.com/294里说了

对于非限制性开壳层（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等序号超过原先占据轨道总数的轨道也纳入的话，不会影响结果。

Reolli

sobereva 发表于 2025-5-22 04:31
当前情况下，Gaussian里的轨道序号顺序是：[alpha占据轨道]
alpha和beta都是81个电子，因此占据轨道是1-16 ...

哦哦，懂了老师，感谢答疑解惑