MOMAP或Dushin能否计算QM/MM模型的重组能分解

JefferSon

sob老师和各位老师同学好，

我最近在做光学相关的研究，看了帅志刚老师的一篇文章“Hydrogen Bonding-Induced Morphology Dependence of Long-Lived Organic Room-Temperature Phosphorescence: A Computational Study， https://pubs.acs.org/doi/10.1021/acs.jpclett.6b01156”，其中提及了不同聚集体环境（单晶，无规，多晶）的重组能计算，但是里面貌似没有提及他的重组能计算时是仅计算该聚集体下分子孤立结构的重组能还是说直接以QM/MM模型计算。同时，我还想请教一下大家是否成功使用过DUSHIN计算QM/MM模型的重组能分解，此处粘贴一位老师同学的帖子（http://bbs.keinsci.com/thread-27073-1-1.html），貌似无法实现。sob老师在帖子下方回复的冻结原子方法我之前试过，但是会有一个报错问题就是（-1）固定后是没有频率的，但貌似dushin会先根据原子数预估了比如说你应该要有10个频率，此时因为固定后频率比如说只有8个了，软件就会直接报错“错误: 数据不足于填充下三角"，导致矩阵没有办法形成，计算就会卡在输出s_matr.dat文件步骤，不知道大家有遇见这个问题和知道解决方法嘛。谢谢大家，祝大家工作顺利！