四丁基氢氧化铵（TBAH）分子中TBA+和OH-的力场参数选择的问题？

模小助

拓扑处理这里有个容易踩坑的细节，拿出来讨论一下。

我首先是按照TOP程序中例2（ em=gd3bj scrf=(solvent=h2o)下优化了TBA+和OH-的结构；在关键词：B3LYP/def2TZVP em=gd3bj下计算了TBA+和OH-的气相下的单点能；

主要想确认：这种拓扑/力场处理是否合理，有没有更稳妥的设置方式？