计算化学公社

标题: 换方法分子结构破坏 [打印本页]

作者
Author: lei234 时间: 2016-4-30 06:34
标题: 换方法分子结构破坏
图1是用M06-2X/6-31+G**优化成功的分子，换B3LYP/6-31G**优化后不能收敛，结构变成图2，其中一个H2O掉了下来，试了读结构1的chk优化；直接opt优化；减小步长优化iop(1/8=3)；最后都会变成图2的结构，请问还有什么办法吗？主要想用B3LYP/6-31G**优化后读HOMO-LUMO gap,读Alpha occ. eigenvalues和Alpha virt. eigenvalues的值，是不是只有优化后才有这个值，只算单点应该没有这个值吧？如果我用图1的结构，用B3LYP/6-31G*做pop=full,这样应该也不行吧？读出的轨道应该还是用M06-2X/6-31+G**算的吧？

各位大神帮帮我吧~~

(, 下载次数 Times of downloads: 24) (, 下载次数 Times of downloads: 30)

作者
Author: sobereva 时间: 2016-4-30 06:57
本身B3LYP描述氢键就不如M06-2X好，没必要再用B3LYP优化。而且那个水本身位置和朝向也不好，没法和上面的氧形成氢键，掉下来也很正常。

只是读取轨道能直接用B3LYP算单点就完了，根本用不着优化，pop=full都可以不写。

作者
Author: lei234 时间: 2016-5-3 01:16

sobereva 发表于 2016-4-30 06:57
本身B3LYP描述氢键就不如M06-2X好，没必要再用B3LYP优化。而且那个水本身位置和朝向也不好，没法和上面的氧 ...

那要加溶剂模型吗？这个反应是在溶液中进行的，但好像一般算HOMO-LUMO band gap都在气相下的。

作者
Author: sobereva 时间: 2016-5-3 04:33

lei234 发表于 2016-5-3 01:16
那要加溶剂模型吗？这个反应是在溶液中进行的，但好像一般算HOMO-LUMO band gap都在气相下的。

加

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