如何进行COSMO、SCRF溶剂化计算

软件版本为AMS2019.3及之前的版本,AMS2020以后的版本请参考链接:COSMO、SCRF溶剂化的计算

在使用的时候,例如结构优化、光谱计算、过渡态计算等等,只需要在原先不考虑溶剂化的基础上,增加如下所示的设置即可。

这是ADF的两种设置最简单的溶剂化,也是最常用的溶剂化方法。

COSMO溶剂化模型如上图所示,选好溶剂化理论以及溶剂即可。其中Solvent epsilon和Solvent radius不需要设置——除非在溶剂列表中没有该溶剂。

SCRF的溶剂化,需要设置Default radii为Noodelman。

溶剂化能

COSMO的溶剂化能,参考:ADF计算溶剂化自由能(Free Energy of Solvation)

SCRF的溶剂化能,在out文件的尾部,搜Net可以搜到,Net即表示溶剂化能,例如:

                   ANALYSIS (Energy in kcal/mol)
 Cycle     Cost       Gainr     GainP      Net       EleDip
 ------------------------------------------------------------
  0      0.00    -24.16      0.00    -24.16   -24.165    1.2979
  1      7.69    -41.82      0.00    -41.82   -34.126    2.2712
  2     13.94    -48.54      0.00    -48.54   -34.595    2.9494
  3     17.07    -51.57      0.00    -51.57   -34.500    3.4460
  4     18.73    -53.13      0.00    -53.13   -34.402    3.8248
  5     19.64    -53.98      0.00    -53.98   -34.342    4.1198
  6     20.15    -54.45      0.00    -54.45   -34.300    4.3551
  7     20.45    -54.73      0.00    -54.73   -34.279    4.5462

其中第6列数据Net即表示在SCRF的迭代过程中,溶剂化能的不同值,最后一次为-34.279kcal/mol。

注意

SCRF搜索Net会有多个,因为SCRF本身是一个自洽迭代的过程,因此会多次计算溶剂与溶液,直到达到平衡。在out文件中,也会对应的出现多次类似上诉内容,以最后一次出现的为准(最后一次出现的,才是收敛的结果)