一、研究背景

金属卟啉催化的C−H键羟基化反应作为细胞色素P450酶催化羟基化的仿生小分子版本，备受关注。高价铁氧（Fe⁴⁺=O）卟啉-阳离子自由基物种（Compound I）被公认为该类反应中的亲电氧化剂，其催化乙苯C−H键羟基化被认为遵循攫H/氧回弹（HAA/OR）机制，但关于自由基对的性质及其对立体化学的影响，现有研究结论存在分歧，需通过分子动力学模拟深入探究。

二、研究核心问题

1. 铁卟啉催化乙苯羟基化反应的关键步骤及决速步是哪一个过程？
2. 不同自旋态（双重态、四重态、六重态）能量表面对反应路径和产物形成有何影响？
3. 自由基对的动态行为（回弹、分离、重排）如何影响反应的立体选择性？
4. 溶剂环境对反应轨迹和自由基对寿命有怎样的调控作用？

三、研究方法

1. 密度泛函理论（DFT）计算：采用Gaussian 16程序，通过（U）B3LYP/6-31G*-LANL2DZ水平计算自由能曲线、过渡态结构及热力学校正值，结合SMD溶剂化模型模拟隐式溶剂环境。
2. 准经典分子动力学（MD）模拟：通过Singleton的Progdyn代码与Gaussian 16对接，在298.15 K下从优化过渡态出发，以1 fs步长传播轨迹，分析键形成/断裂阈值及无重穿事件的有效轨迹。

四、研究结果

1、反应机制与决速步

- 铁氧物种抽取乙苯的氢原子是反应的决速步，生成铁羟基物种与苄基自由基的自由基对，随后自由基对重组形成羟基化产物。

- 双重态能量表面上，氧回弹步骤无势垒；四重态/六重态等高自旋能量表面上，氧回弹步骤存在明显势垒（2.8-3.9 kcal/mol）。

2、不同自旋态的反应活性

- 低自旋（双重态）主导反应：气相中45%的反应轨迹直接生成羟基化产物，隐式溶剂模型中该比例提升至56%，反应兼具动态协同性与分步性（C−H键断裂与C−O键形成的时间间隔为41-619 fs）。

- 高自旋（四重态/六重态）贡献极小：98%-100%的反应轨迹导致自由基对分离，仅2%以下轨迹生成羟基化产物，属于能量分步过程。

3、自由基对动态行为与立体选择性

- 自由基存在两种路径：一是快速回弹形成羟基化产物，二是分离后可能发生苄基自由基构型翻转（生成立体异构产物）或β-氢抽取（形成苯乙烯副产物）。

- 自旋污染影响：双重态能量表面存在一定自旋污染，导致自由基对分离轨迹比例被高估，且会侵蚀反应的立体选择性，这是仿生金属卟啉催化乙苯羟基化对映选择性较低的重要原因。

4、溶剂环境的影响

隐式溶剂（二氯甲烷）相比气相环境，不仅提高了氧回弹轨迹的比例，还延长了自由基对的平均寿命，使C−H键断裂与C−O键形成的平均时间间隔从99 fs增至156 fs。

五、结论

铁卟啉催化乙苯羟基化反应中，双重态能量表面的反应路径比高自旋态更有利，其反应机制兼具动态协同与分步特征；氢原子抽取是决速步，氧回弹无势垒，而自由基对的分离与构型翻转会影响产物立体选择性；溶剂环境可调控反应轨迹分布与自由基对寿命。高自旋态虽能形成自由基对，但对羟基化产物的贡献可忽略不计。

文案&审核：小H

文章链接：https://doi.org/10.1021/jacs.3c03773