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

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

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

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

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

谢谢老师

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

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

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

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

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

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

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

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

用的哪种键级？哪个计算级别？
涉及突变的几个点做一下波函数稳定性测试看是否稳定，也看看能量曲线是否有突跃。