计算化学公社

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

[Multiwfn使用咨询] 键级曲线的研究请教

[复制链接 Copy URL]

10

帖子

0

威望

39

eV
积分
49

Level 2 能力者

本人按照“通过键级曲线和ELF-LOL-RDG等值面动画研究化学反应过程”教程做了键级曲线如图,我的反应是(戊烯、烯丙基甲醚、丙烯酸甲酯的)双键打开变成单键。遇到以下疑问,希望老师能解答:1.做出的图的双键键级没有达到2,单键键级低于1,这是正常的吗?
2.键级曲线在鞍点突变,并没有老师做得那么平滑上升,是哪里出错了吗?

K}5OZCCDZMXTL4TRMEBK.png (7.69 KB, 下载次数 Times of downloads: 54)

K}5OZCCDZMXTL4TRMEBK.png

6万

帖子

99

威望

5万

eV
积分
120080

管理员

公社社长

2#
发表于 Post on 2015-10-10 00:42:12 | 只看该作者 Only view this author
用的哪种键级?哪个计算级别?
涉及突变的几个点做一下波函数稳定性测试看是否稳定,也看看能量曲线是否有突跃。
北京科音自然科学研究中心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

10

帖子

0

威望

39

eV
积分
49

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2015-10-10 01:20:05 | 只看该作者 Only view this author
sobereva 发表于 2015-10-10 00:42
用的哪种键级?哪个计算级别?
涉及突变的几个点做一下波函数稳定性测试看是否稳定,也看看能量曲线是否有 ...

这是我IRC输入命令 :#p irc=(maxpoints=50,calcall,StepSize=3) ub3lyp/6-311++g(d,p) nosymm geom=connectivity

IRC的曲线是没有突变的,请问老师我可以参考什么资料看波函数稳定性测试?

3097

帖子

29

威望

1万

eV
积分
17094

Level 6 (一方通行)

4#
发表于 Post on 2015-10-10 02:08:42 | 只看该作者 Only view this author
elerun 发表于 2015-10-10 01:20
这是我IRC输入命令 :#p irc=(maxpoints=50,calcall,StepSize=3) ub3lyp/6-311++g(d,p) nosymm geom=conn ...

Gaussian 关键词 stable

6万

帖子

99

威望

5万

eV
积分
120080

管理员

公社社长

5#
发表于 Post on 2015-10-10 02:34:39 | 只看该作者 Only view this author
elerun 发表于 2015-10-10 01:20
这是我IRC输入命令 :#p irc=(maxpoints=50,calcall,StepSize=3) ub3lyp/6-311++g(d,p) nosymm geom=conn ...

涉及到断键,会产生对称破缺态,应该再加上guess(mix,always),这样通常才能让每一步都得到对称破缺态。
参考
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82
北京科音自然科学研究中心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

10

帖子

0

威望

39

eV
积分
49

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2015-10-12 20:18:56 | 只看该作者 Only view this author
sobereva 发表于 2015-10-10 00:42
用的哪种键级?哪个计算级别?
涉及突变的几个点做一下波函数稳定性测试看是否稳定,也看看能量曲线是否有 ...

sob老师,请问是用高斯命令#p stable=opt freq guess(mix,always) ub3lyp/6-311++g(d,p) geom=connectivity 算突变的几个点吗?

6万

帖子

99

威望

5万

eV
积分
120080

管理员

公社社长

7#
发表于 Post on 2015-10-12 21:15:02 | 只看该作者 Only view this author
elerun 发表于 2015-10-12 20:18
sob老师,请问是用高斯命令#p stable=opt freq guess(mix,always) ub3lyp/6-311++g(d,p) geom=connectivi ...

不是。计算方法和基组还写原来的,弄成个普通单点任务,然后写stable关键词,看输出文件中提示波函数是否稳定,比如# B3LYP/6-311++G(d,p) stable。如果显示不稳定,说明原先的计算级别不合适。这往往是因为键拉长过程中没有考虑闭壳层波函数变化成对称破波函数缺导致的,需要按我文中的方式得到对称破缺态。
北京科音自然科学研究中心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

10

帖子

0

威望

39

eV
积分
49

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2015-10-13 21:30:53 | 只看该作者 Only view this author
sobereva 发表于 2015-10-12 21:15
不是。计算方法和基组还写原来的,弄成个普通单点任务,然后写stable关键词,看输出文件中提示波函数是否 ...

谢谢老师

10

帖子

0

威望

39

eV
积分
49

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2015-10-14 10:37:13 | 只看该作者 Only view this author
sobereva 发表于 2015-10-12 21:15
不是。计算方法和基组还写原来的,弄成个普通单点任务,然后写stable关键词,看输出文件中提示波函数是否 ...

我按照老师的方法做了波函数稳定性测试,鞍点附近突变的点结果都显示 The wavefunction is stable under the perturbations considered.是不是代表那些突变的点技术上是没有错误的呢?

6万

帖子

99

威望

5万

eV
积分
120080

管理员

公社社长

10#
发表于 Post on 2015-10-14 12:45:59 | 只看该作者 Only view this author
elerun 发表于 2015-10-14 10:37
我按照老师的方法做了波函数稳定性测试,鞍点附近突变的点结果都显示 The wavefunction is stable under  ...

建议你先作个最基本的乙烷C-C断裂的图,没问题再做你自己的。
有些细节远程很难说清楚。
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2025-8-13 02:58 , Processed in 0.155152 second(s), 23 queries , Gzip On.

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