求助Gaussian中使用内置溶剂定义SMD混合溶剂的方式？

曹丽佳

各位老师，前辈们：
       本人请教Gaussian中内置溶剂定义混合溶剂的方式。溶剂：体积比为1:9=H2O:DMSO。问题如下：
       sobereva老师在http://bbs.keinsci.com/thread-3345-1-1.html 博文中提到对于几种内置溶剂的混合情况，就把SMD溶剂参数相应地按照体积比混合来定义新溶剂即可。那么学生本人定义的溶剂参数，例如：eps的计算数值=eps(H2O)*0.1+eps(DMSO)*0.9，其它的epsinf、HBondAcidity等参数类似如此计算，关键词的书写：scrf=(read,SMD,solvent=DMSO)，不知道这样理解是否正确？

        非常感谢各位老师、前辈的不吝赐教。

路人甲2

同学，解决了吗，我也想知道这个怎么搞

cokie

路人甲2 发表于 2026-5-21 10:41
同学，解决了吗，我也想知道这个怎么搞

sob 老师在博文里 1.5.3 已经写的很明确了：
1. 想自己定义SMD极性和非极性全部参数怎么做（以及关键词怎么写）；
2. 想只定义极性部分，怎么做（以及关键词怎么写）
3. 混合溶剂怎么做
4. 混合溶剂中包含水怎么做。

1.5.3 自定义溶剂及混合溶剂（同时考虑极性和非极性部分贡献时）
在SMD溶剂模型下计算时如果也要考虑溶剂的非极性部分，那么必须把溶剂的完整的SMD参数全都进行定义，即写scrf(read,SMD,solvent=generic)，然后在输入文件末尾空一行写比如
eps=11.5
epsinf=2.0449
HBondAcidity=0.229
HBondBasicity=0.265
SurfaceTensionAtInterface=61.24
CarbonAromaticity=0.12
ElectronegativeHalogenicity=0.24

其中SurfaceTensionAtInterface是表面张力，单位是cal/mol/A^2。CarbonAromaticity是芳香度，是芳环上的碳原子数占所有非氢原子数的比例，ElectronegativeHalogenicity是卤素度，是卤原子占所有非氢原子数的比例。HBondAcidity和HBondBasicity分别是Abraham酸度和碱度，程序内置溶剂的这俩参数来自于Abraham的文章，对于一种新的溶剂这是没法查到的。所以，简单来说，对于一种全新溶剂，是没有严格办法考虑非极性部分参数的。在笔者来看，碰上这种情况时的相对最佳的做法是写比如scrf(read,SMD,solvent=AAA)，然后照常定义eps和epsinf，这里AAA是与你当前要用的溶剂特征最接近的自带的溶剂，这样非极性部分的参数就会借用AAA溶剂的凑合着。虽然这样不严格，但总比完全不考虑非极性部分强一点。

通过自定义SMD溶剂模型参数做SMD计算的实例见《通过SMD溶剂模型描述离子液体溶剂环境的方法》（http://sobereva.com/431）。

如果是几种内置溶剂混合的情况，就把SMD溶剂参数相应地按照体积比混合来定义新溶剂即可。内置溶剂的SMD参数在明尼苏达溶剂描述符数据库http://comp.chem.umn.edu/solvation/mnsddb.pdf里可以查到。

注意水在SMD中是特殊的，没有直接对应的非极性参数，有人问如果溶剂是水和其它溶剂分子混合构成的怎么去定义混合溶剂。在我来看没有绝对理想的办法去定义这种混合溶剂。对于SMD算这种混合溶剂下的溶解自由能，我认为一个可以接受的做法是在水中和在另一个纯溶剂中分别计算溶解自由能，然后按照水和其它溶剂分子的体积比来权重得到混合溶剂下的溶解自由能。

如果你不是用的SMD，而是其它溶剂模型结合dis cav rep选项来考虑溶剂非极性部分的贡献，那更没法严格地自定义新溶剂或者定义混合溶剂。