计算化学公社

标题: 键级曲线的研究请教 [打印本页]

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

作者
Author: sobereva 时间: 2015-10-10 00:42
用的哪种键级？哪个计算级别？
涉及突变的几个点做一下波函数稳定性测试看是否稳定，也看看能量曲线是否有突跃。

作者
Author: elerun 时间: 2015-10-10 01:20

sobereva 发表于 2015-10-10 00:42
用的哪种键级？哪个计算级别？
涉及突变的几个点做一下波函数稳定性测试看是否稳定，也看看能量曲线是否有 ...

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

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

作者
Author: liyuanhe211 时间: 2015-10-10 02:08

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

作者
Author: sobereva 时间: 2015-10-10 02:34

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

作者
Author: elerun 时间: 2015-10-12 20:18

sobereva 发表于 2015-10-10 00:42
用的哪种键级？哪个计算级别？
涉及突变的几个点做一下波函数稳定性测试看是否稳定，也看看能量曲线是否有 ...

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

作者
Author: sobereva 时间: 2015-10-12 21:15

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。如果显示不稳定，说明原先的计算级别不合适。这往往是因为键拉长过程中没有考虑闭壳层波函数变化成对称破波函数缺导致的，需要按我文中的方式得到对称破缺态。

作者
Author: elerun 时间: 2015-10-13 21:30

sobereva 发表于 2015-10-12 21:15
不是。计算方法和基组还写原来的，弄成个普通单点任务，然后写stable关键词，看输出文件中提示波函数是否 ...

谢谢老师

作者
Author: elerun 时间: 2015-10-14 10:37

sobereva 发表于 2015-10-12 21:15
不是。计算方法和基组还写原来的，弄成个普通单点任务，然后写stable关键词，看输出文件中提示波函数是否 ...

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

作者
Author: sobereva 时间: 2015-10-14 12:45

elerun 发表于 2015-10-14 10:37
我按照老师的方法做了波函数稳定性测试，鞍点附近突变的点结果都显示 The wavefunction is stable under ...

建议你先作个最基本的乙烷C-C断裂的图，没问题再做你自己的。
有些细节远程很难说清楚。

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