| 是这样的,即使Cl把H抽走了,由于这是个孤立体系,可能还是HCl的H指向丢了H的C的时候能量更低,可能做不到HCl直接分开出去。或者说真正的反应的生成物得是溶剂化后的HCl。要描述好这种反应,或许得上混合溶剂模型(隐式+显式),甚至配个小溶剂球跑metadynamic了(你CV可以是两个键长的线性组合、设一个CV就行),ORCA就可以实现这种带边界约束和增强采样的分子动力学模拟,详见手册https://www.faccts.de/docs/orca/ ... n.html#metadynamics |
| 参与人数Participants 1 | eV +5 | 收起 理由Reason |
|---|---|---|
|
| + 5 |
| 处理Cl离子或者Cl自由基,最好加几个显式溶剂分子,将其用卤键或氢键固定住... |
dbxhtw 发表于 2026-3-17 20:35 那样的话你可以参考4#的建议,试试更换扫描变量 根据经验这个攫氢一般是C-H-Cl为直线,你也可以试试冻结这个键角 如果确实难找到你也可以试试用ORCA的NEB-TS |
|
(1)键长不仅可以从长往短扫,也可以从短往长扫。 (2)4-5(C-H)、12-5(Cl-H)这两个键长都可以扫扫看。 |
| 参与人数Participants 1 | eV +5 | 收起 理由Reason |
|---|---|---|
|
| + 5 | 我很赞同 |
KazusaT 发表于 2026-3-17 18:23 这两个点都试了,振动方向都不是正确的 |
|
能量骤变感觉主要是因为小分子的构型有很大的改变,能明显看出来C1-C4的构型不一样了 你可以试试用这两个结构做opt=TS看能不能找到 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-4-13 19:00 , Processed in 0.304343 second(s), 26 queries , Gzip On.