|
|
本帖最后由 ldatea 于 2020-6-29 19:28 编辑
你说的二聚体,是指一个溶质分子+一个溶剂分子 吗?这种一般不叫二聚体,二聚体通常指同一种分子形成的复合物。
如果你还是用DMSO溶剂,隐式溶剂模型就可以了,DMSO是典型偶极溶剂,偶极矩较大,且正电中心在硫上,藏在内部,负电中心露在外面(氧),只考虑静电作用的时候,阳离子在DMSO中的溶剂化作用强,阴离子溶剂化作用弱。现在你的溶质也可以近似看成一个正电荷中心藏在内部,负电荷中心露在外面的偶极子。
如果用能和酮形成氢键的溶剂(比如水,位阻小的醇,胺,氯仿等),最最最少得带一个溶剂(体现氢键),然后对这个复合物用隐式溶剂模型。具体细节请看sob老师的博文(《谈谈隐式溶剂模型下溶解自由能和体系自由能的计算》)http://sobereva.com/327
我其实现在都不敢确定你这个体系到底会有几个构象,目前只是优化出两个极小值点而已。
我这边算出来的结果:
6-31g*(D3BJ)几何优化,然后用ORCA的wB97M-V/ma-def2-TZVP(RIJCOSX加速) 计算单点(这里用弥散函数的目的是为了和之后要算的阴离子体系用同一个基组,本身中性分子加弥散函数是不划算的,还有一个次要目的是估算后面算阴离子的时间,就是看wB97M-V/ma-def2-TZVP能不能算得动)
O-C-C-H二面角168.106的构象,wB97M-V/ma-def2-TZVP单点-733.841302485153Hartree,M05-2X/6-31g*气相单点-733.793625308Hartree,SMD溶剂模型下M05-2X/6-31g*单点-733.812249149Hartree
O-C-C-H二面角8.964的构象,wB97M-V/ma-def2-TZVP单点-733.835717944143Hartree,M05-2X/6-31g*气相单点-733.787252648Hartree,SMD溶剂模型下M05-2X/6-31g*单点-733.803927949Hartree第三个减第二个就是溶剂化自由能(状态是默认的298.15K,1atm,溶液标态1mol\L)。
|
|