我想用Gaussian 09计算离子液体[BIm][Tf2N]分别在气相、考虑溶剂效应下的阴阳离子间的相互作用能以及相互作用吉布斯自由能。计算步骤如下:
1. 在气相下分别优化[BIm][Tf2N]、[BIm]+、[Tf2N]- 的结构并做频率分析确定无虚频,这个时候获得三者在气相的E_gas和Ggas;
所用关键词: # opt freq b3lyp/cc-pvdz empiricaldispersion=gd3 temperature=313.15
2. 根据第1步获得的气相优化结构,用smd模型,溶剂选Methanol,考虑溶剂效应分别对[BIm][Tf2N]、[BIm]+、[Tf2N]- 进行优化并做频率分析确定无虚频,这个时候获得三者在考虑溶剂效应下的E_sol和 G_sol;
所用关键词: # opt freq b3lyp/cc-pvdz scrf=(smd,solvent=methanol) empiricaldispersion=gd3 temperature=313.15
3. 相互作用能和相互作用吉布斯自由能的计算:
3.1 相互作用能deltaE:
deltaE_gas = E_gas([BIm][Tf2N]) – E_gas([BIm]+) – E_gas([Tf2N]-) = -100.8 kcal/mol
deltaE_sol = E_sol([BIm][Tf2N]) – E_sol([BIm]+) – E_sol([Tf2N]-) = -17.1 kcal/mol
3.2 相互作用能deltaG:
deltaG_gas = G_gas([BIm][Tf2N]) – G_gas([BIm]+) – G_gas([Tf2N]-) = -85.9 kcal/mol
deltaG_sol = G_sol([BIm][Tf2N]) – G_sol([BIm]+) – G_sol([Tf2N]-) = -3.1 kcal/mol
请问:
1. 这样的计算过程是否可以?
2. 计算出来的气相和考虑溶剂效应下的相互作用能和相互作用吉布斯自由能相差这么大是否合理?
最后优化得到的结构见附件。