计算化学公社

标题: 求助贴-CuI溶剂化能量值问题 [打印本页]

作者
Author: shirely 时间: 2018-8-2 15:28
标题: 求助贴-CuI溶剂化能量值问题
各位大神好，最近用CuI做催化剂算一个体系，但是在做溶剂化的时候只有CuI的能量值变化很大，而其他结构能量值变化较小，从而导致基于零点的能量差值很大，一直找不到原因，希望能在这里得到解答，谢谢各位！对CuI单独优化的输入文件如下
%chk=CuI
%mem=800mb
%nprocshared=4
#p wb97xd/LANL2DZ opt freq=noraman

a3

0 1
Cu 0.00000000 0.00000000 -1.56565600
I 0.00000000 0.00000000 0.85668000

溶剂化下单点能计算的输入文件如下
%chk=CuI
%mem=800mb
%nprocshared=4
#p wb97xd/SDD scrf(SMD,solvent=Toluene) scf=tight geom=check guess=read

a3

0 1

作者
Author: Accelerator 时间: 2018-8-2 23:53
不要对I用LanL2DZ，精度太差而且可能会导致结构不合理（亲身经历）。不如对Cu和I就用def2系列
我不知道wB97xD算碘化亚铜怎么样，这种情况下不同泛函对Cu-I键的离子成分描述差异会很大
freq=noraman不用写
不知道碘化亚铜在溶液中的存在形式是什么，不知直接用单体作为模型合理性如何；如果没有加配体的话我觉得可以从碘化亚铜与底物的配合物作为能量的起点
此外我觉得单独的碘化亚铜单体溶剂化能算不准是可以理解的

作者
Author: sobereva 时间: 2018-8-3 00:22
怎么一会儿用lanl2dz一会用SDD？

"做溶剂化"指代不明，到底是溶剂下优化，还是算溶解自由能什么的得说请

noraman、scf=tight是完全多余，看
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331

作者
Author: shirely 时间: 2018-8-3 10:06

sobereva 发表于 2018-8-3 00:22
怎么一会儿用lanl2dz一会用SDD？

"做溶剂化"指代不明，到底是溶剂下优化，还是算溶解自由能什么的得说请 ...

谢谢，是溶剂下优化，优化用的LANL2DZ基组然后改用大基组进行溶剂化，其他机构能量变化相差不大，只有CuI变化很大。

作者
Author: shirely 时间: 2018-8-3 10:07

Accelerator 发表于 2018-8-2 23:53
不要对I用LanL2DZ，精度太差而且可能会导致结构不合理（亲身经历）。不如对Cu和I就用def2系列
我不知道wB9 ...

谢谢，按您的方法先试试

作者
Author: sobereva 时间: 2018-8-4 02:46

shirely 发表于 2018-8-3 10:06
谢谢，是溶剂下优化，优化用的LANL2DZ基组然后改用大基组进行溶剂化，其他机构能量变化相差不大，只有CuI ...

“进行溶剂化”依然指代不明
如果你的目的就是溶剂下优化结构，直接用SDD带着隐式溶剂模型做就行了，完全没必要先用lanl2DZ预处理一遍，不仅多余，还会因此显著增加整个过程的耗时
如果是比较气相和溶剂下优化的结构的差异，也应当都用SDD

欢迎光临计算化学公社 (http://bbs.keinsci.com/)