计算化学公社

标题: ECD结果拟合的不好 [打印本页]

作者
Author: xuj 时间: 2017-11-14 16:16
标题: ECD结果拟合的不好
本帖最后由 xuj 于 2017-11-16 12:35 编辑

请教各位老师，我最近计算的了几个化合物的ECD，其中有两个化合物ECD计算结果和实测拟合的不好，图片和相应的程序如下，请帮忙分析一下，万分感激！图一的命令：

%nprocshared=12
%mem=10GB
%chk=cs2A.chk
# opt freq b3lyp/6-31+g(d,p) scrf=(solvent=methanol) geom=connectivity

--link1--
%nprocshared=12
%mem=10GB
%chk=cs2A.chk
# td=(nstate=15) b3lyp/6-31+g(d,p) geom=allcheck

图二的命令：
%nprocshared=12
%chk=CS1831ECD2.chk
%rwf=CS1831ECD2.rwf
%nosave
%mem=6320MB
#p td=(nstates=15,full,singlet) b3lyp/aug-cc-pvdz scrf=(solvent=methanol)
geom=connectivity scfcon=6 scfcyc=300

作者
Author: sobereva 时间: 2017-11-14 17:32
TDDFT任务垃圾关键词太多，scfcon=6还导致精度下降。加弥散根本没必要
改成
#p td(nstates=15) PBE1PBE/6-311G* scrf(solvent=methanol)

注意考虑构象的影响，参见
使用Multiwfn绘制构象权重平均的光谱
http://sobereva.com/383

为了和实验对比较好，对理论谱的激发能、强度乘以一定的系数因子消除系统偏差是可以的

作者
Author: xuj 时间: 2017-11-14 18:53

sobereva 发表于 2017-11-14 17:32
TDDFT任务垃圾关键词太多，scfcon=6还导致精度下降。加弥散根本没必要
改成
#p td(nstates=15) PBE1PBE/6 ...

谢谢sobereva老师的指点！以上ECD图都是加权以后的结果，另外请问是否能推荐一下ECD计算中基组和函数的选择方面的参考书或文献？非常感谢您！

作者
Author: sobereva 时间: 2017-11-14 21:17

xuj 发表于 2017-11-14 18:53
谢谢sobereva老师的指点！以上ECD图都是加权以后的结果，另外请问是否能推荐一下ECD计算中基组和函数的选 ...

不是“函数”是“泛函”
小分子用PBE0，大共轭体系用CAM-B3LYP
基组用6-311G*足矣，用更大了没帮助

作者
Author: xuj 时间: 2017-11-15 20:49

sobereva 发表于 2017-11-14 21:17
不是“函数”是“泛函”
小分子用PBE0，大共轭体系用CAM-B3LYP
基组用6-311G*足矣，用更大了没帮助

谢谢sobereva老师，我按照您的推荐更改了基组和泛函对图二对应的化合物重新计算了一下，结果如图，麻烦您再帮忙分析一下！非常感谢您！

作者
Author: sawyer 时间: 2017-11-16 22:17
如果计算值和实测值一直相差很大的话，还要考虑是不是结构定的有点问题。

作者
Author: xuj 时间: 2017-11-20 17:33

sawyer 发表于 2017-11-16 22:17
如果计算值和实测值一直相差很大的话，还要考虑是不是结构定的有点问题。

谢谢您的回复，是需要进行紫外矫正

作者
Author: 欢乐多 时间: 2019-5-14 07:37
本帖最后由欢乐多于 2019-5-14 07:53 编辑

xuj 发表于 2017-11-20 17:33
谢谢您的回复，是需要进行紫外矫正

层主，您好，最近在学习ECD，计算结果拟合，X轴差别较大，就像您说的需要紫外校正，如何进行紫外矫正，是图上选择ECD曲线》右键》proprties 进行修改吗？其有标准没有

作者
Author: 欢乐多 时间: 2019-5-22 12:20
本帖最后由欢乐多于 2019-5-22 12:22 编辑

sobereva 发表于 2017-11-14 21:17
不是“函数”是“泛函”
小分子用PBE0，大共轭体系用CAM-B3LYP
基组用6-311G*足矣，用更大了没帮助

sob老师，这是自己摸索1个月得出的计算ECD结果，可是跟实验结果并不吻合，与发文章的结果距离甚是遥远，不知道问题出在哪了，希望老师能指点一二。

作者
Author: sobereva 时间: 2019-5-22 12:33

欢乐多发表于 2019-5-22 12:20
sob老师，这是自己摸索1个月得出的计算ECD结果，可是跟实验结果并不吻合，与发文章的结果距离甚是遥远， ...

优化时候连极化函数都没有等于搞笑，计算结果毫无意义
用aug-cc-pVDZ远不如用def-TZVP，或更好的def2-TZVP。加没多大用的弥散前干嘛不提升到3-zeta？仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

泛函尝试PBE0，甭用B3LYP，仔细看
乱谈激发态的计算方法
http://sobereva.com/265（http://bbs.keinsci.com/thread-415-1-1.html）

作者
Author: 欢乐多 时间: 2019-5-22 12:40

sobereva 发表于 2019-5-22 12:33
优化时候连极化函数都没有等于搞笑，计算结果毫无意义
用aug-cc-pVDZ远不如用def-TZVP，或更好的def2-TZ ...

好的，老师，我再试试吧

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