离子液体[BIm][Tf2N] pKa值的计算

我想计算[BIm][Tf2N]（结构见图1）这种离子液体的pKa值，主要关注图1中红圈所标注的H，我拟定了如下的计算步骤：下面的说明[BIm][Tf2N]以AH表示，[BIm][Tf2N]解离H以后的结构（见图2）以A-表示

1、采用smd模型，溶剂为水，优化AH、A-的结构，并确定无虚频，关键词: opt freq b3lyp/6-31G* scrf(smd,solvent=water)  empiricaldispersion=gd3；

2、取1中优化得到的结构做气相下的频率分析，来获得AH、A-气相下的自由能（Ggas(AH)、Ggas(A-)）以及单点能（Egas(AH)、Egas(A-)），关键词：freq m052x/6-31G* empiricaldispersion=gd3 ;
3、取1中优化得到的结构，用smd模型，溶剂为水计算AH、A-气相下的单点能Esol(AH)、Esol(A-)，关键词：m052x/6-31G* scrf(smd,solvent=water) empiricaldispersion=gd3 ;
4、deltaG(AH)=Esol(AH)-Egas(AH), deltaG(A-)=Esol(A-)-Egas(A-), Ggas(H+)=-6.28 kcal/mol, deltaGsol(H+)=-265.9 kcal/mol, deltaG(1atm->1M)=1.89 kcal/mol;
5、deltaGsol=Ggas(A-)+Ggas(H+)-Ggas(AH)+deltaG(A-)+deltaGsol(H+)-deltaG(AH)+deltaG(1atm->1M);
6、pKa=deltaGsol/RT;

请问：
(1) 这整个计算流程是否可行？
(2) 第2、3步的计算都是在第1步的结构上进行，这样是否可以？
(3) 第2步进行频率计算的时候可能会出现虚频的情况，那么这一步计算出来的自由能是否可用？
(4) 每一步的方法/基组组合是否合理？有没有更好的方法/基组组合推荐？

谢谢。