Gaussian用ONIOM做QM/MM参数怎么补？

铭律凉

我想要优化在一大堆相同的分子中两个分子（Dimer）被激发时的结构
为了节约时间，选择对除了Dimer以外的分子选用MM Amber，对Dimer用RODFT计算激发态（认为这个五重激发态在所有五重态中能量最低）
设置完ONIOM以及方法后，提示我缺少一些参数

缺的都是键角信息

我看过了，需要自己添加，缺的太多了我不太可能手动补，我倒是可以写一个脚本自己去把缺的键角补了，但是力常数我不知道该区哪里找
我看到有例子TAO package可以快速补，但是使用了以后parmlookup一直在重复
Use of uninitialized value &angle1atm1 in uc at ………………
Use of uninitialized value &angle1atm2 in uc at ………………
Use of uninitialized value &angle1atm3 in uc at ………………

始终不给我输出
求助各位希望能找到解决办法  谢谢

之后检查发现TAO package需要你写入原子坐标信息内带有电荷信息，我就统一随便设了一个值都为0.00，然后parmlook给我输出了
HrmBnd1 CM CM CM 0.00000 0.00000
HrmBnd1 CM CA CA 0.00000 0.00000
HrmBnd1 CM CA HA 0.00000 0.00000
不知道能不能用这个结果加上去?

铭律凉

sobereva 发表于 2021-10-24 03:40
有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必 ...

非常不好意思，已经改了，请麻烦删除之前的楼层与这条回复

sobereva

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。