求助~460个原子有机体系计算HOMO和LUMO能级需要怎么计算

LipengWu

请教一下各位有过相关计算经验的老师们，我有个~460个原子（C225H225N5O5)的有机分子需要计算HOMO和LUMO能级，想问一下需要怎么计算，或者是否可以截取部分片段约100原子来计算呢？感谢各位老师

sgwzq

如果你确实需要计算HOMO和LUMO，你可以把这个问题发在量子化学区。

只给分子式不够，你需要描述清楚具体是什么样的分子，有没有对称性，算HOMO和LUMO干什么。

这么多C，看起来像是芳香体系。如果必须算的话，几何优化大概是要用半经验方法了（比如加色散校正的PM7或者xTB等等），正常的DFT方法算不动。

另外，DFT的HOMO-LUMO 和 Hartree-Fock的HOMO-LUMO又不一样。对于DFT方法，Koopman's Theorem不一定成立。详见https://roaldhoffmann.com/sites/default/files/fromd6/444s_0.pdf

sobereva

普通学术问题不要发到大杂烩版块。这次给你移动了，以后注意！

注意把下文看了，搞清楚HOMO和LUMO的意义，然后交代清楚算HOMO、LUMO目的何在，并且把体系特征充分描述清楚，否则不可能准确回答你的问题。很多情况一些新人说要计算HOMO、LUMO，但是他们实际真正需要的并不是HOMO和LUMO
正确地认识分子的能隙(gap)、HOMO和LUMO
http://sobereva.com/543（http://bbs.keinsci.com/thread-16758-1-1.html）


并且把你当前的计算资源说清楚。仔细看下文末尾部分了解怎么正确提问此类问题
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚、完整
http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html）

yaonaite

我本地的配置是8375C双路服务器，虚拟机中运行的是物理核心32个，内存60G，研究对象是多环芳烃
几何优化b3lyp/6-31g**能算动，需要几天的样子，很少到一周
单点加到3zeta级别也能算动，不加弥散函数时间不长，不到10小时吧，加上弥散时间会长很多乃至不收敛，但如果能收敛，1-2天也能完

如果有可以合理简化的结构，可以简化，比如对研究的性质没啥影响的叔丁基改成甲基，甚至直接删掉之类的简化

江月待何人

至少需要说明是什么样的体系，或者给出结构式，才能判断简化哪些原子；
一般思路是先预优化，再做结构优化，最后算高精度单点能，仅需要数据的话，看输出文件即可，如果需要等值面图，使用GaussView或者Multiwfn生成，如果还是一头雾水，建议好好搜索相关博文看看，直接参与量子化学培训班是最快最省时的入门方式