计算化学公社

标题: 求助：四氢萘氢转移过渡态和Beta键键能计算问题 [打印本页]

作者
Author: zhangxiance 时间: 2023-4-13 17:42
标题: 求助：四氢萘氢转移过渡态和Beta键键能计算问题
本帖最后由 zhangxiance 于 2023-4-13 22:06 编辑

计算四氢萘发生氢转移的能垒和四氢萘（带有正碳离子）Beta键键能，有几个问题不能确定，在此向大家请教。
首先，分别优化四氢萘和C10+，使用b3lyp/6-31++g(d,p) freq opt，得到最优化结构后，让两个产物摆在一起，得到IM1，具体结构为下图。之后使用 opt(calcfc,noeigen,TS,tight) freq
找到过渡态TS，最后找到中间IM1。但对于IM1，带有伯正碳离子的C10+总会被苯环吸引，最终连到苯环上的碳上。TS倒是有且只有一个负频，而且负频的振动是H-从四氢萘向C10+移动，但IM2的结构优化出来与TS差太多，不知道这个反应路径是否准确？能量计算使用自由能Sum of electronic and thermal Free Energies。
(, 下载次数 Times of downloads: 5) (, 下载次数 Times of downloads: 4) (, 下载次数 Times of downloads: 7)
之后计算四氢萘+（正碳离子在1位）Beta键断裂所需的能量，首先使用guess(mix,always) opt=modredundant pop=always nosymm将Beta键两端C分离，再opt freq 得到优化结构，但无论怎么优化，都得不到伯碳正碳离子，只能得到C-C-C形成一个三角形，这个结构能够理解，但不知道这个结果计算出的键能是否准确，四氢萘+和Beta键断裂后的结构如下图，两个结构使用的电荷为1，多重态为1。计算键能使用的能量为Sum of electronic and thermal Enthalpies。
(, 下载次数 Times of downloads: 3) (, 下载次数 Times of downloads: 4)
还有一个问题，我加了scf=noincore，但计算的时候还只是使用一个核，计算速度很慢。
#SBATCH --nodes=1
#SBATCH --ntasks-per-node=8
#SBATCH --mem=16G
请大家不吝赐教，再次感谢！！！

作者
Author: dnlx 时间: 2023-4-13 22:15
一、模拟你的体系的任何一个分子离子都不应该加弥散函数，尤其不应该给氢加，查社长关于弥散函数的文章；
二、碳正离子加成到四氢化萘上很有可能不存在能垒，所以优先于夺氢的反应。

作者
Author: zhangxiance 时间: 2023-4-13 22:33

dnlx 发表于 2023-4-13 22:15
一、模拟你的体系的任何一个分子离子都不应该加弥散函数，尤其不应该给氢加，查社长关于弥散函数的文章；
...

非常感谢，我去学习一下。

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