计算化学公社

标题: 求助过渡态搜索中遇到的IRC结果不对和报错问题 [打印本页]

作者
Author: jdg 时间: 2026-3-30 14:21
标题: 求助过渡态搜索中遇到的IRC结果不对和报错问题
目的：研究聚酰亚胺分子断键并形成新化学键的过程，在研究2号和3号化学键时遇到了问题 (, 下载次数 Times of downloads: 0)
问题：①针对2号化学键，断键后的自由基可进行优化且结果没有虚频，因此以拉伸键长为条件进行了柔性扫描，发现能量最高点是另一化学键断裂造成的，我以扫描结果中能量最高点左侧的几个模型进行过渡态优化，结果因有多个虚频而报错，是否说明这一过程中不存在过渡态或者我应如何修改结构以贴近过渡态 (, 下载次数 Times of downloads: 0)
②针对3号化学键，断键后的自由基存在虚频，依旧以拉伸键长为条件进行了柔性扫描，我以结构中的能量最高点做了过渡态优化，只有一个虚频，但是IRC的结果只往两边各走了一个点，我已尝试在过渡态优化时加入tight关键词，以及在IRC计算中改变Maxpoints和Stepsize，结果没有改善 (, 下载次数 Times of downloads: 0)

作者
Author: Uus/pMeC6H4-/キ 时间: 2026-3-30 14:37
有什么理由（含中间体表征与有出处的文献的合理计算）认为：
1. 实验中聚合物的这些化学键都能断裂；
2. 每根化学键都对应了“只有这根化学键断裂而其他化学键不改变”的反应（特别是考虑到涉及环状结构时“键解离能”定义的模糊性）；且
3. 反应在所用计算级别和/或客观现实中都有过渡态（而不是无势垒反应的情形）？

是否有类似http://sobereva.com/82的对自旋极化单重态的恰当考虑？

作者
Author: jdg 时间: 2026-3-30 15:42
1.对于化学键的断裂，我是在做了键级分析和单点能计算（算解离能）后，选部分解离能低的键来进行计算，认为在实际反应中有较大可能断裂
2.就扫描结果来看，当我改变目标键的键长后，会伴随有其他键的断裂或重新生成问题，但我考虑的是目标键的断裂是否存在过渡态，我也觉得扫描时单一改变键长这一参数是不准确的，但我也没有更好的思路
3.目前的计算我使用的是2-zeta基组，因为我看之前的博文

http://sobereva.com/460

已经够用，所以没有考虑计算级别的问题。并且我也在思考这个问题，该如何判断反应是无势垒还是过渡态初猜结构有问题导致的错误
4.我在做扫描的时候考虑过开壳层单重态计算，上面的结果图就是，过渡态优化时没有考虑，还是三重态，我会进行尝试

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