计算化学公社

标题: 如何分开读取开壳层plot文件（.33）中的轨道系数？ [打印本页]

作者
Author: lukas 时间: 2025-1-25 16:32
标题: 如何分开读取开壳层plot文件（.33）中的轨道系数？
本帖最后由 lukas 于 2025-1-28 11:23 编辑

各位好，我希望读取33文件中轨道的展开系数信息构建，但33文件中不同自旋不分开记录，如何分别读取？两部分的信息的分界点在何处，对于不同自旋多重度时候如何确定？ (, 下载次数 Times of downloads: 0)

作者
Author: sobereva 时间: 2025-1-26 11:32
问题要描述完整，我假定你说的是NBO plot文件
看Multiwfn源代码文件fileIO.f90中读入NBO plot文件的子程序subroutine read31，特别是if (chartemp(1:11)==" ALPHA SPIN".or.chartemp(1:11)==" alpha spin") then后面的部分

作者
Author: lukas 时间: 2025-1-26 11:53
本帖最后由 lukas 于 2025-1-26 11:55 编辑

sobereva 发表于 2025-1-26 11:32
问题要描述完整，我假定你说的是NBO plot文件
看Multiwfn源代码文件fileIO.f90中读入NBO plot文件的子程序 ...

是NBO的，但发现33文件似乎没有alpha和beta的索引，折中的方法是先使用脚本对半添加beta spin的分隔，再进行您说的读取，但我不确定是否33文件中的排列是如此，先生成α再生成β的原子轨道展开系数？此外还想请教一下，这些轨道是包含空轨道吗？我读过您整理的NBO资源，但可能查找的不细心，没找到对文件格式具体说明的部分。

作者
Author: sobereva 时间: 2025-1-26 13:46
找个简单体系如H2测试便知
NAO属于自然轨道一大类，不严格分占据和非占据轨道。轨道占据数都是非整数的。

作者
Author: lukas 时间: 2025-1-28 11:22

sobereva 发表于 2025-1-26 13:46
找个简单体系如H2测试便知
NAO属于自然轨道一大类，不严格分占据和非占据轨道。轨道占据数都是非整数的。

我采用了典型的氧气三重态作为计算案例，得到了l607.exe产生的33文件，设置的uhf/3-21g，其中基组为18个，但在33中也只有18个轨道按理来说开壳层轨道是基组的两倍，这令我很困惑，此外33文件最后出现了几行我不了解的数据，麻烦老师能否给我解答。我先附上33文件。

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