请问您已经找到这个旋转能垒了吗?我目前也遇到这个问题,能否讨论指教下,谢谢! |
sobereva 发表于 2020-9-5 02:46 好的,我再深入学习下,谢谢您 |
心辉 发表于 2020-9-4 09:59 我说的绝对不行是指用B3LYP算有机体系能量问题在精度上如今绝对不行 在这里我已经说得很明确了,不要试图自寻死路 坚持使用B3LYP算有机体系的人的下场 http://bbs.keinsci.com/thread-12773-1-1.html |
心辉 发表于 2020-9-4 09:59 绝对不行可能也不至于,但是性价比并非最优是肯定的,因为单点能计算(即使加大了基组)的时间远小于几何优化和频率计算的时间,所以单点能最好不仅加大基组,也换用更准确的泛函,或干脆用波函数方法。 之所以大家普遍用B3LYP/6-311+G(d,p)之类算单点,是因为十几年前市面上还没有双杂化泛函,杂化泛函也没有什么比B3LYP好的,DLPNO-CCSD(T)更没有的情况下,没有多少方法能比B3LYP/6-311+G(d,p)性价比更高。但是现在来看,这种算法已经过时了。像你的体系,如果R1、R2各不超过20个原子的话,完全可以DLPNO-CCSD(T)/def2-TZVP,精度比B3LYP-D3(BJ)/6-311+G(d,p)能高好几倍,单点能计算时间相比B3LYP增加10倍左右,但即使这样也不见得比小基组结构优化和频率计算的时间多少。 |
447951397 发表于 2020-9-4 10:50 好的。我再研究研究 |
谢谢 |
心辉 发表于 2020-9-4 09:59 可以 |
心辉 发表于 2020-9-4 09:59 你怕是对社长有什么误解,社长的意思是优化和振动分析时用级别低的基组没问题。而之后计算能量差时都用用高档一点的基组相对准确一些,不同基组下能量相减是无意义的。 |
sobereva 发表于 2020-9-3 21:58 我有点迷惑,社长推荐我B3LYP-D3(BJ)优化过渡态还行,而算能量绝对不行。 我有点没理解,我之前计算是在B3LYP-D3(BJ)/6-31G*下优化和振动并算到自由能校正项,然后在6-311G+**下算到单电子能,两者结合得到最后的自由能。 现在社长说能量绝对不行,我理解的是B3LYP-D3(BJ)/6-31G*下优化和振动并算到自由能校正项,然后用M06-2X/6-311G+**算单电子能,然后两者结合。这种方式可以么? |
sobereva 发表于 2020-9-3 21:58 谢谢您,受教了 |
wzkchem5 发表于 2020-9-3 20:08 好的谢谢您 |
扫描只适合确定过渡态初猜的位置,真正的旋转势垒应当用过渡态与极小点能量求差来得到。B3LYP-D3(BJ)优化过渡态还行,而算能量绝对不行,仔细看 坚持使用B3LYP算有机体系的人的下场 http://bbs.keinsci.com/thread-12773-1-1.html 2 算能量肯定得用SMD模型表现溶剂效应。优化的话对这个体系无所谓 3 仔细看 谈谈如何通过势垒判断反应是否容易发生 http://sobereva.com/506(http://bbs.keinsci.com/thread-14313-1-1.html) |
你这个扫描曲线有很大的hysteresis(能量稳步爬升然后突然掉下去),扫出来的最高点不一定能很好地代表过渡态。你可以试试倒着扫,然后你就会发现最高点的位置变了。 出现hysteresis一般是因为扫描的时候有些变化很大的自由度你没有考虑到,比如你这个体系,可能从你给的结构旋转到另一个enantiomer是分两步的,第一步左上边苯环翻到纸面下面去,但是R1暂时不翻上来,也就是萘和吲哚都是弯曲结构;第二步R1再翻上来。如果存在这种现象的话,那么什么基组大小之类的都是次要的,得先把反应途径搞得定性正确再说 |
喵星大佬 发表于 2020-9-3 19:12 这个我之前有考虑过计算过渡态,但是有的时候会有两种情况:1.结构不能稳定在能量最高的点,而是滑落到稳定结构了。2.计算出来的结构算频率,得不到虚频或者是虚频的振动方向不是旋转的方向。请教下,一般这种情况要怎么处理? |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-27 07:38 , Processed in 0.179835 second(s), 26 queries , Gzip On.