关于QM计算选取含金属蛋白初始结构与HOMO−LUMO gap的关系的询问

18217265596

The QM calculations are conducted with the proposed charge states and final charges are chosen, for which the optimized structure most closely shows integer orbital occupations and a maximized HOMO−LUMO gap.
作者在优化一系列含金属蛋白，用 Schrödinger’s PROPKA pKaprediction tool 预测了ph7下的质子化状态和电荷，然后基于 maximized HOMO−LUMO gap选取了合适的输入结构
请问这点的原理大致是什么？

wzkchem5

这个的关键在于轨道是整数占据，HOMO-LUMO gap是决定占据数的，而直接决定选哪个结构的是占据数而不是HOMO-LUMO gap。
如果靠人来决定一个蛋白的电荷数，一般就是把每个残基以及各种辅因子的电荷加起来。但是如果让程序自动做的话，程序不一定知道每种辅因子是多少电荷。所以就要每个电荷数都跑一个SCF计算，因为实际的蛋白大概率是单参考态体系，所以轨道占据数接近整数，因此挑那个让轨道占据数最接近整数的电荷数，就最有可能是对的