怎么比较方便的得到波函数分析文件？

路人癸

求助各位大佬，老师，就我已知的方法是用G16做几何结构优化的时候，得到的chk文件转变为fch，但是用高斯做几何结构优化，分子数多的话算不了，那我就得不到完整体系的fch文件，求

wal

换xtb

zjxitcc

你可能想说的是“原子数多的话算不了”，不是“分子数多的话算不了”。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硬算。

Uus/pMeC6H4-/キ

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

sobereva

“方便”用词不当，看你的意思，你只是想低耗时得到大体系的波函数文件用于Multiwfn分析

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

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