如何分开读取开壳层plot文件（.33）中的轨道系数？

lukas · 发表于 Post on 2025-1-25 16:32:07

本帖最后由 lukas 于 2025-1-28 11:23 编辑

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

NBO-OXY036.33 (8.02 KB, 下载次数 Times of downloads: 0)

sobereva · 发表于 Post on 2025-1-26 11:32:10

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

lukas · 发表于 Post on 2025-1-26 11:53:44

本帖最后由 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资源，但可能查找的不细心，没找到对文件格式具体说明的部分。

sobereva · 发表于 Post on 2025-1-26 13:46:43

找个简单体系如H2测试便知
NAO属于自然轨道一大类，不严格分占据和非占据轨道。轨道占据数都是非整数的。

lukas · 发表于 Post on 2025-1-28 11:22:01

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[NBO相关] 如何分开读取开壳层plot文件（.33）中的轨道系数？

浏览过的版块