关于高斯中计算酶催化磺胺甲恶唑的电荷和多重度的设定问题

春秋

向老师请教自旋多重度的问题：
一个有机污染物（磺胺甲恶唑）电荷中性单重态，一个酶活性位点（P450酶）电荷中性二重态，我想研究酶低自旋状态催化该有机污染物的反应能垒，也就是把污染物和酶放在一个计算体系中计算，此时体系的净电荷数仍然是0，但自旋多重度应设为几合适呢？

sobereva

看具体怎么催化的。如果此过程Fe部分配位方式没怎么变的话，自旋多重度还是2。如果有显著变化，可以同时考虑在4重态势能面上搜索反应过程的驻点，并和二重态情况的进行比较

wzkchem5

另外注意，因为这个体系里可能有互相呈反铁磁耦合的自旋，只把自旋多重度设对是不够的，仍然可能收敛到错误的电子态（通俗来讲就是自旋向上和向下的单电子数都是对的，但是单电子分布位置错误）。所以算完以后必须检查自旋密度分布等，确认收敛到了想要的电子态，如果收敛到了错误的电子态，需要通过交换轨道等方式收敛到正确的电子态

春秋

非常感谢老师们的回复！我补充了文献中的反应过程机理，供老师们参考。
sob老师，Fe部分配位方式应该发生了改变，您的意思是计算四重态的极小点能量去和二重态的极小点进行比较吗？我还有点纠结计算体系在酶的基础上还添加了污染物（自旋多重度为1），如果我要分别计算酶在高自旋（四重态）和低自旋（低自旋）状态下的体系能量变化，是否直接选取自旋多重度为4和2计算就可以呢？需要考虑这个自旋多重度为1的污染物吗？
wxk老师，您应该是提醒开壳层和闭壳层的选择，非常感谢，请问正确的电子态有统一标准还是要通过文献进行判断呢？