请教分子中氨基在298.15K下的氢解离能的计算

老少年

各位老师好：计算新手求教，我现在需要计算防老剂分子氨基的氢解离能（用来比较他们的防老化性能），在计算的时候遇到了以下问题
基组都是用的 b3lyp/6-31g*     

1.因为是要计算298.15K下的氢解离能，不知道可不可以使用频率计算完以后的自由能。

2.我先是把去掉氢以后的片段进行优化再进行频率计算，但是这样得出来的片段他的建发生了改变，这样我不知道计算完的自由能（和其他的能量）还能不能用。

3.我把去掉氢的片段直接进行频率计算，那样得出来的四个收敛判断标准都是NO，并且还出现虚频。
  
      优化完片段以后键发生了改变         

wzkchem5

1. 解离能和解离自由能是两码事。应该先确定你需要的是解离能、解离焓还是解离自由能，然后从频率计算的输出文件里读相应的校正量，加到大基组的单点能结果上面。
2. 仔细阅读http://sobereva.com/414
3. 去掉氢必须重新优化才能算频率

另外还有一个问题，做这个计算的时候必须加色散校正，必须加溶剂模型，且必须做构象搜索，否则结果不可信