Uus/pMeC6H4-/キ 发表于 2025-6-11 08:35 大致明白了，但还是有些不懂的地方，我还得多学习！ 感谢您的耐心讲解和宝贵时间！ |
本帖最后由 Uus/pMeC6H4-/キ 于 2025-6-11 08:39 编辑 ipoy 发表于 2025-6-11 07:43 我没有说“O-Li-O不算弱相互作用”的意思，但讨论的重点是环硅氧烷的反应的话这种因素相对次要（除非要比较如不同碱金属离子对反应的影响之类，有这种前提得及早讲明）。注意当溶剂极性较大时，带与不带隐式溶剂模型所得结论可能差别很大，特别是涉及离子的反应在气相和溶液相会有定性不同。 社长在http://sobereva.com/413已讲解了什么是弱相互作用。在凝聚相中不存在暴露的单原子金属离子，周围必定有一圈配位原子（包括主要配体，也包括显式溶剂分子和部分阴离子等等）将第一配位层饱和，比如同样计算镧系的http://sobereva.com/581例子里除了冠醚还有几个水分子。此时金属和周围原子至少有很显著的吸引作用（甚至是配位键这样强到能算化学键的作用），我觉得再称其为“弱相互作用”有点别扭。不考虑所有直接接触金属离子的配体是大忌，会得到定性错误的结论，这点在论坛里讨论吸附能/结合能的很多帖子里都有说明。 闭壳层/开壳层之别与化学键/弱相互作用之别是两个不同的层面的问题，没有必然相关性。 |
ipoy 发表于 2025-6-11 06:45 当前讨论的不是闭壳层主族物种成键断键的化学反应么，为啥会往弱相互作用去想呢？真是弱相互作用情况的话B3LYP-D3(BJ)也有比较好的精度。 另外，如果反应是在溶液环境中发生的，记得加隐式溶剂模型。 |
Uus/pMeC6H4-/キ 发表于 2025-6-10 17:31 感谢您的细心指导！！ int这个结构我确实没有找到，我之前打算从ts1出发找int结构，看来不太合理。我会根据您的建议先找int然后往两边扩，这个是我之前没想到的。 我也觉得问题在泛函选择上，根据我原来的经验（有可能都不对），用明尼苏达算弱相互作用力比较合理，所以我第一选择就是M06-2X，我会根据您的建议用B3lyp试一下。 再次感谢！ |
ipoy 发表于 2025-6-10 16:57 两条路径的TS1和TS'不能看成同一结构，二者主要区别在环内硅氧键是否断开。考虑下方路径的话，我建议先从五配位硅中间体开始，优化极小点找到INT结构，然后分别只拉长环外硅氧键优化TS1过渡态、只拉长环内硅氧键优化TS2过渡态，再对两个过渡态跑IRC并优化端点得到R与P。 当前file1_TS1(SiMe2).log的振动虚频67.93 i cm-1与file3_TS1(SiCH2CH2CF3).log的振动虚频78.92 i cm-1看起来都不太对劲，有可能并不对应真实的反应过程。在积分格点（注意不叫晶格）不够精细时，M06-2X等明尼苏达泛函就是容易出这种毛病，所以拿它做几何优化和IRC等描述势能面的任务不是最佳选择，可以换成B3LYP（结合D3(BJ)色散校正）等。 |
Uus/pMeC6H4-/キ 发表于 2025-6-10 16:17 太感谢您的认真回答了，非常感动，而且非常清晰。 我理解您上述的意思，我承认从R到P不是一个基元反应，可以看到其实TS1和TS’是一个结构，我在做二甲基取代的环硅氧烷过渡态搜索的时候发现了TS1（file1.gjf, file.log)，之前确实是按照从R到TS’到P猜测的，但是做了IRC发现TS’(file2.gjf, file2.log与TS1结构相同)只连接了R和int，并不能到达P。所以我认为从R到TS1到P这个过程不存在（帖子中的结构式我简化了）。 当我将二甲基硅中的一个甲基变成三氟丙基的时候，在搜索TS1发现只有一个虚频，能够正常终止，但是震动方向明显不对（file3.gjf, file3.log是三氟丙基中的亚甲基的伸缩振动），因此我认为三氟丙基的结构影响了过渡态的能量（使其降低），因此更难搜索到了（不知道我的理解对不对）。如果咱们的结构猜测正确的话，三氟丙基取代的TS1一定存在，但是不知道用什么方法能够把他计算出来。对于TS2我认为可能不存在，由int可能直接重排成P（不知道我的理解对不对，请指正） 
file.zip
(1.26 MB, 下载次数 Times of downloads: 1)
|
本帖最后由 Uus/pMeC6H4-/キ 于 2025-6-10 16:18 编辑 ipoy 发表于 2025-6-10 15:14 我怀疑我的意思没传达到位。按我的表达重新绘制简化的反应途径如下： 
其中R是反应物，P是产物，INT是中间体，均为势能面极小点；TS1、TS2、TS'是过渡态，均为势能面一阶鞍点。INT中的五配位硅用绿色圈出。我想表达的意思是，从R到P不是经过上方只有TS'一个过渡态的单步基元反应，而是经过下方有TS1加成过渡态、INT1中间体、TS2消除过渡态的两步基元反应。与碳的SN2取代反应是成键断键协同机理不同，半径更大的硅的取代反应更可能是先加成后消除的分步机理，这在Arrow Pushing in Inorganic Chemistry等无机化学书上有提到。对TS1结构跑IRC连通的两端是R和INT，而对TS2结构跑IRC连通的两端是INT和P。 请参照此图的标记将问题描述得更清晰，并将相关输入输出文件打包上传为帖子附件，说明结构与文件对应关系、问题出在哪个文件。 |
Uus/pMeC6H4-/キ 发表于 2025-6-10 14:47 您说的很对，是通过一步过渡态成键，得到一个硅中间体，然后这个中间体再变成现性结构（这里面可以通过IRC看出来），IRC在往right方向走的时候是图中所示的结构，当往另一方向走的时候不能得到线性结构，而后一个类似楼主说的5配位结构。 但是我的问题还存在，就是在搜索第一步过渡态的时候，当甲基用三氟丙基取代的时候就不能成功找到震动方向合理的过渡态。 |
| 该过程涉及“环外新的硅氧键生成”和“环内旧的硅氧键断裂”，确认二者是同一个基元反应里发生的吗？有可能先是硅醇的氧与环硅氧烷的硅经一步过渡态成键，得到一个五配位硅中间体（极小点），然后再经过第二个过渡态开环，得到最终开链产物。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-13 17:11 , Processed in 0.449974 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites