计算化学公社

标题: 求解H2S与S2反应的过渡态的计算 [打印本页]

作者
Author: sylvester 时间: 2022-10-7 15:04
标题: 求解H2S与S2反应的过渡态的计算
在cp2k中跑了H2S与S2的反应，反应过程如下图，我现在用高斯寻找过渡态，就是第一步反应的过渡态，首先优化了S2与H2S的结构，均利用B3LYP/6-311G++**,优化后发现s2的双键变为单键，把优化后的H2S与S2拖入高斯，根据跑出的产物结构，在寻找过渡态，首先寻找出的过渡态没有虚频，且成一条直线，相距很远，后来寻找出的过渡态虚频很小（只有-10左右），而且寻找过渡态过程中能量一直在降低，试过了很多初猜结构，最后都是这样的，而且S2与H2S的距离也相距很远。请问这是出现了什么问题。

作者
Author: wzkchem5 时间: 2022-10-7 15:25
首先一个问题，你为什么觉得你的计算结果有问题？违背了什么实验结果吗？
其次，你有没有考虑到S2分子的基态是三重态

作者
Author: chands 时间: 2022-10-7 16:31
本帖最后由 chands 于 2022-10-7 17:18 编辑

你这个反应是气相反应还是溶液相反应，溶液中形成多硫化物很常见，我不知道气相中是否有类似反应。
还有，你明明是用Gaussian算的，CP2K的输入文件一匹布那么长。
(莫非你研究火山口的反应？）

作者
Author: sylvester 时间: 2022-10-8 10:18

wzkchem5 发表于 2022-10-7 15:25
首先一个问题，你为什么觉得你的计算结果有问题？违背了什么实验结果吗？
其次，你有没有考虑到S2分子的基 ...

因为虚频很小，而且过渡态的能量跟优化后的能量相差不大。我是考虑了S2的基态的三重态的。所以在高斯的输入中是0 3

作者
Author: sylvester 时间: 2022-10-8 10:19

chands 发表于 2022-10-7 16:31
你这个反应是气相反应还是溶液相反应，溶液中形成多硫化物很常见，我不知道气相中是否有类似反应。
还有， ...

你好，我这个就是在cp2k中跑AIMD，然后判断出这是化学反应，才来找过渡态的，这个是气相反应。

作者
Author: wzkchem5 时间: 2022-10-8 15:23

sylvester 发表于 2022-10-8 03:18
因为虚频很小，而且过渡态的能量跟优化后的能量相差不大。我是考虑了S2的基态的三重态的。所以在高斯的输 ...

虚频小、过渡态不高都未必是有问题的，有的过渡态本来就是这样的，除非实验表明这一步能垒不可能这么低，或者KIE实验表明这一步的虚频不可能这么小，那种情况下再去考虑计算可能算错了的问题

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