计算化学公社

标题: 怎么比较方便的得到波函数分析文件？ [打印本页]

作者
Author: 路人癸 时间: 2025-10-5 19:04
标题: 怎么比较方便的得到波函数分析文件？
求助各位大佬，老师，就我已知的方法是用G16做几何结构优化的时候，得到的chk文件转变为fch，但是用高斯做几何结构优化，分子数多的话算不了，那我就得不到完整体系的fch文件，求

作者
Author: wal 时间: 2025-10-5 19:08
换xtb

作者
Author: zjxitcc 时间: 2025-10-5 19:39
本帖最后由 zjxitcc 于 2025-10-5 19:43 编辑

你可能想说的是“原子数多的话算不了”，不是“分子数多的话算不了”。150原子以下的常规体系用Gaussian的DFT做结构优化和频率分析不是难题；一两百个原子的常规体系用ORCA的DFT配合RI加速，也不是难题。500原子以上的有机、生物大分子，适合采用线性标度分块方法（例如GEBF, MFCC, MBE）结合DFT进行相应的计算，算一两千原子没问题。如果是上千原子的常规体系，想算得快+精度要求不高，可以使用半经验量子化学方法，例如xTB程序里的GFN2-xTB方法。注意，xTB是程序，它里面实现了很多方法例如GFN-FF, GFN0-xTB, GFN1-xTB, GFN2-xTB等等，程序与方法不能混淆。

如果研究的体系是周期性的（石墨烯、MOF、固体、金属氧化物等等），应当采用CP2K/QE/VASP等周期性的第一性原理程序进行计算，不适合挖个孤立团簇用Gaussian/ORCA硬算。

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-10-5 19:42
服务器/超算机时及资源多少，分子总数多少、特征如何（基态还是激发态、有无双自由基或多自由基等复杂电子结构、有无过渡金属原子、每个分子总共多少原子等），已经有立体笛卡尔坐标还是要自行从头生成三维结构（包括从字符串或非立体结构文件生成，以及从晶体结构提取等），希望在什么级别优化，要拿波函数文件做什么分析？按http://sobereva.com/620把所有这些细节都讲清楚才能提供建议，另外写shell脚本等做批处理是免不了的。

作者
Author: sobereva 时间: 2025-10-6 11:45
“方便”用词不当，看你的意思，你只是想低耗时得到大体系的波函数文件用于Multiwfn分析

没有具体信息没法准确回答。如果坚持DFT计算，不重要的部分可以截去或用小基组。如果之前一直用的Gaussian，单点计算部分可以改用ORCA靠RIJ（纯泛函）、RIJCOSX（杂化泛函）显著节约耗时。如果显著放宽对波函数质量的要求，用xtb做PTB计算（比GFN2-xTB的波函数质量更好），产生的molden文件用Multiwfn也能分析。

鉴于很多人对PTB不了解，顺带给出北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）里详细讲GFN系列方法的理论知识和xtb程序部分的幻灯片：
(, 下载次数 Times of downloads: 3)

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