计算化学公社

 找回密码 Forget password
 注册 Register
Views: 12189|回复 Reply: 10
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] ECD结果拟合的不好

[复制链接 Copy URL]

36

帖子

0

威望

112

eV
积分
148

Level 2 能力者

本帖最后由 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

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

2#
发表于 Post on 2017-11-14 17:32:56 | 只看该作者 Only view this author
TDDFT任务垃圾关键词太多,scfcon=6还导致精度下降。加弥散根本没必要
改成
#p td(nstates=15) PBE1PBE/6-311G* scrf(solvent=methanol)

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

为了和实验对比较好,对理论谱的激发能、强度乘以一定的系数因子消除系统偏差是可以的
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

36

帖子

0

威望

112

eV
积分
148

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2017-11-14 18:53:17 | 只看该作者 Only view this author
sobereva 发表于 2017-11-14 17:32
TDDFT任务垃圾关键词太多,scfcon=6还导致精度下降。加弥散根本没必要
改成
#p td(nstates=15) PBE1PBE/6 ...

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

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

4#
发表于 Post on 2017-11-14 21:17:57 | 只看该作者 Only view this author
xuj 发表于 2017-11-14 18:53
谢谢sobereva老师的指点!以上ECD图都是加权以后的结果,另外请问是否能推荐一下ECD计算中基组和函数的选 ...


不是“函数”是“泛函”
小分子用PBE0,大共轭体系用CAM-B3LYP
基组用6-311G*足矣,用更大了没帮助
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

36

帖子

0

威望

112

eV
积分
148

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2017-11-15 20:49:31 | 只看该作者 Only view this author
sobereva 发表于 2017-11-14 21:17
不是“函数”是“泛函”
小分子用PBE0,大共轭体系用CAM-B3LYP
基组用6-311G*足矣,用更大了没帮助

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

1510749875(1).png (7.28 KB, 下载次数 Times of downloads: 49)

1510749875(1).png

1

帖子

0

威望

7

eV
积分
8

Level 1 能力者

6#
发表于 Post on 2017-11-16 22:17:04 | 只看该作者 Only view this author
如果计算值和实测值一直相差很大的话,还要考虑是不是结构定的有点问题。

36

帖子

0

威望

112

eV
积分
148

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2017-11-20 17:33:26 | 只看该作者 Only view this author
sawyer 发表于 2017-11-16 22:17
如果计算值和实测值一直相差很大的话,还要考虑是不是结构定的有点问题。

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

147

帖子

4

威望

1583

eV
积分
1810

Level 5 (御坂)

科研泥瓦匠

8#
发表于 Post on 2019-5-14 07:37:57 | 只看该作者 Only view this author
本帖最后由 欢乐多 于 2019-5-14 07:53 编辑
xuj 发表于 2017-11-20 17:33
谢谢您的回复,是需要进行紫外矫正

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

Q0JB6ROX109Z88QL_Z)F9TB.png (218.95 KB, 下载次数 Times of downloads: 44)

Q0JB6ROX109Z88QL_Z)F9TB.png
终日寻春不见春
芒鞋踏破岭头云
归来偶把梅花嗅
春在枝头已十分

147

帖子

4

威望

1583

eV
积分
1810

Level 5 (御坂)

科研泥瓦匠

9#
发表于 Post on 2019-5-22 12:20:22 | 只看该作者 Only view this author
本帖最后由 欢乐多 于 2019-5-22 12:22 编辑
sobereva 发表于 2017-11-14 21:17
不是“函数”是“泛函”
小分子用PBE0,大共轭体系用CAM-B3LYP
基组用6-311G*足矣,用更大了没帮助

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

n421.gif (45.93 KB, 下载次数 Times of downloads: 56)

n421.gif

N421R ECD5.gjf

3.08 KB, 下载次数 Times of downloads: 16

N421R ECD5.LOG

66.94 KB, 下载次数 Times of downloads: 8

N421R opt frq.gjf

3.07 KB, 下载次数 Times of downloads: 20

终日寻春不见春
芒鞋踏破岭头云
归来偶把梅花嗅
春在枝头已十分

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

10#
发表于 Post on 2019-5-22 12:33:44 | 只看该作者 Only view this author
欢乐多 发表于 2019-5-22 12:20
sob老师,这是自己摸索1个月得出的计算ECD结果,可是跟实验结果并不吻合,与发文章的结果距离甚是遥远, ...

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

泛函尝试PBE0,甭用B3LYP,仔细看
乱谈激发态的计算方法
http://sobereva.com/265http://bbs.keinsci.com/thread-415-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

147

帖子

4

威望

1583

eV
积分
1810

Level 5 (御坂)

科研泥瓦匠

11#
发表于 Post on 2019-5-22 12:40:56 | 只看该作者 Only view this author
sobereva 发表于 2019-5-22 12:33
优化时候连极化函数都没有等于搞笑,计算结果毫无意义
用aug-cc-pVDZ远不如用def-TZVP,或更好的def2-TZ ...

好的,老师,我再试试吧
终日寻春不见春
芒鞋踏破岭头云
归来偶把梅花嗅
春在枝头已十分

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-25 03:31 , Processed in 0.192510 second(s), 29 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list