HOCH2CP的理论结构和光谱表征：乙醇腈的磷类似物

《Journal of Molecular Spectroscopy》：Theoretical structural and spectroscopic characterization of HOCH2CP: A phosphorus analogue of glycolonitrile

　　

由于含磷物种在星际介质（ISM）中的搜索因其在预生物化学（prebiotic chemistry）中的重要性而在近年来成为一个引人注目的课题。继在星际介质中检测到乙醇腈分子（HOCH₂CN）之后，研究人员采用高级量子化学方法来预测HOCH₂CP（乙醇腈的含磷类似物）的gauche和anti构象体的结构和光谱性质。具体而言，研究人员使用复合方案（composite scheme）并在CCSD(T)（耦合簇单双激发并包含三重激发微扰）水平上包含非谐校正（anharmonic corrections），详细表征了它们的平衡几何、相对稳定性、转动常数、振转参数（rovibrational parameters）和红外光谱（infrared spectra）。研究人员发现gauche构象体比anti构象体稳定1.22 kcal mol^?1，且上述构象体的互变势垒（interconversion barrier）仅为0.35 kcal mol^?1（约175 K），表明两种构象体可能共存于稠密冷分子云（dense cold molecular clouds）。此外，研究人员还研究了可能导致其形成的可能反应。

磷是生命必需元素，在生物体内丰度较高（如细菌中P/H约10^?3），但宇宙中磷的丰度极低（P/H=2.8×10^?7），主要来源于超新星爆发。尽管已检测到多种含磷星际分子（如PN、PO、HCP、CP等），磷化学在星际介质（ISM）中仍理解不足。2019年，预生物分子乙醇腈（HOCH₂CN）在太阳型原恒星IRAS 16293–2422 B中被探测到，其形成被认为与冰晶颗粒表面反应及气相反应（如OH+CH₂CN路径）相关。乙醇腈的等价磷类似物羟基甲基磷乙炔（HOCH₂CP）尚未被检测，但作为含磷物种，其结构和光谱数据对星际搜索具有重要意义。为此，研究人员采用高精度量子化学方法，系统预测了HOCH₂CP两种构象体（gauche和anti）的平衡几何、相对稳定性、转动常数、振转参数、红外光谱及形成路径。论文发表在《Journal of Molecular Spectroscopy》。

研究人员主要采用以下关键技术方法（不超过250字）：
- 密度泛函理论（DFT）计算：使用双杂化B2PLYP-D3泛函，结合aug-cc-pVTZ和aug-cc-pVQZ基组，表征构象体及其互变势垒。
- 高级从头算：在CCSD(T)理论水平上（包括非谐校正），采用复合方案（composite scheme）和aug-cc-pCVTZ等基组，获取高精度平衡几何、转动常数、离心畸变常数及振动光谱参数。
- 红外光谱预测：基于谐波和非谐分析，计算红外谱图及绝对强度。
- 形成路径研究：考虑中性-中性气相反应（如P+H₂COCH₂、OH+CH₂CPH等），采用CCSD(T)水平计算反应势垒和放热性。

研究结果：
- **Conformational landscape, structure and stability**：通过DFT和CCSD(T)计算确认HOCH₂CP存在两种构象体——gauche（两个对映体，HOCC二面角约±50°）和anti（Cs对称，HOCC二面角180°）。gauche构象体比anti稳定1.22 kcal mol^?1（CCSD(T)-F12/cc-pCVTZ水平），互变势垒仅0.35 kcal mol^?1（约175 K），表明在稠密冷分子云（如T<10 K）中两种构象体可共存。gauche构象体的平衡几何中P≡C键长约1.545 ?，C–O键长约1.397 ?。
- **Rotational spectroscopy**：研究人员预测了gauche和anti构象体的转动常数（A、B、C）及所有六阶离心畸变常数（如Δ_J、Δ_JK、δ_J等），并给出偶极矩分量（gauche构象体μ_a=1.21 D, μ_b=1.52 D, μ_c=1.18 D；anti构象体μ_a=3.09 D, μ_b=1.28 D, μ_c=0）。这些数据可直接用于指导实验室转动光谱探测。
- **Vibrational spectroscopy**：通过非谐计算（CCSD(T)/aug-cc-pVTZ水平）得到两种构象体的基频振动模式和红外谱图。gauche构象体最显著的红外峰位于P≡C伸缩振动区域（约1300 cm^?1，强度约200 km mol^?1）和O–H伸缩振动区域（约3500 cm^?1，强度约65 km mol^?1）；anti构象体类似但峰位偏移。这些光谱特征可用于实验室红外鉴定。
- **Formation routes**：研究人员探索了多条可能的气相反应路径产生HOCH₂CP，包括P+H₂COCH₂、OH+CH₂CPH、CH₂OH+CP等。其中最有利路径为OH+CH₂CPH（CH₂CPH为磷乙炔衍生物），该反应在CCSD(T)水平计算为放热（约?50 kcal mol^?1）且无显著势垒，认为可能在星际条件下有效。

在讨论部分，研究人员将HOCH₂CP与氮类似物HOCH₂CN进行了比较：两者构象稳定性相似（gauche均比anti稳定），但HOCH₂CP的互变势垒更低（0.35 vs. ~2 kcal mol^?1），表明HOCH₂CP的构象化更易进行。偶极矩方面，HOCH₂CP的anti构象体具有较大的a-型偶极矩（3.09 D），利于转动光谱的a-型跃迁检测。此外，HOCH₂CP的P≡C伸缩振动频率（~1300 cm^?1）显著低于HOCH₂CN的C≡N伸缩（~2240 cm^?1），反映了P≡C键的较弱强度。结合形成路径分析，HOCH₂CP可能通过类似的OH+CH₂PH反应生成，但其在稠密分子云中的丰度取决于前体物种（如CP、HCP）的可用性。

研究结论（翻译）：
在这项工作中，研究人员研究了乙醇腈（glycolonitrile）的含磷类似物羟基甲基磷乙炔（HOCH₂CP），详细描述了其构象体的平衡几何、相对稳定性、转动常数、振转参数和红外光谱。最初，研究人员采用双杂化B2PLYP-D3泛函在DFT水平表征构象体及其互变势垒。随后，研究人员进行了高级从头算，在CCSD(T)水平上使用复合方案（composite scheme），并广泛考虑基组扩展。研究人员确认gauche构象体是能量最低的构象形式，比anti构象体稳定1.22 kcal mol^?1，且两构象体间互变势垒仅为0.35 kcal mol^?1（≈175 K），表明gauche和anti构象体可能共存于稠密冷分子云中。研究人员提供了高精度转动光谱常数（包括转动常数、离心畸变常数和偶极矩分量），以及非谐校正的红外光谱数据。最后，研究人员提出了一个可产生HOCH₂CP的新气相反应路径，即OH+CH₂CPH，该反应在星际条件下在能量上可行且放热。本报告的数据预计将促进HOCH₂CP在实验室中的检测，并为未来在星际介质中搜寻这一新含磷物种提供钥匙。