计算化学公社

标题: 运行NMR=CSGT相关问题 [打印本页]

作者
Author: Reolli 时间: 2025-5-21 00:00
标题: 运行NMR=CSGT相关问题
按照基于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

作者
Author: sobereva 时间: 2025-5-21 05:28
上传输入输出文件

作者
Author: Reolli 时间: 2025-5-21 09:43

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

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

作者
Author: sobereva 时间: 2025-5-22 04:31
当前情况下，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等序号超过原先占据轨道总数的轨道也纳入的话，不会影响结果。

作者
Author: Reolli 时间: 2025-5-22 09:35

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

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

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