基于高阶耦合簇与语法驱动异构体生成的星际CnHN（n=2–5）含氮碳链电子结构与振转光谱基准研究

《Physical Chemistry Chemical Physics》：Electronic structure and spectroscopic properties of astrochemically relevant CnHN isomersOpen Access

　　

CnHN异构体家族包含氰基（–C≡N）与异氰基（–N≡C）官能团，是构成星际介质（ISM）中一类重要含氮碳链。研究人员采用SmilX语法驱动异构体生成结合高阶从头算方法，包括显相关CCSD(T)-F12与CCSD(T)-F12/RI（分辨率恒等近似），对n=2–5选定低能中性CnHN异构体进行量子化学表征。研究共刻画21个最低能单重态与三重态物种的相对能量、平衡几何、电荷分布、偶极矩、垂直电离势（VIP）、振动频率与旋转常数（Ae、Be、Ce）。除C4HN全局极小为环-链结构外，低能异构体多为含端基–C≡N的线型或准线型物种；并首次表征此前未报道的C5HN异构体5-5（三元HC3环连接–CCN取代基）。偶极矩在结构相关系列中随碳链增长而增大；VIP主要由局域成键环境与端基连接性决定而非链长。非谐振动处理较谐波频率更贴近实验值，准线型物种因大振幅运动偏差较大。计算旋转常数与实验值偏差为0.02%–2.15%。该数据集为未探测CnHN异构体提供光谱与电子性质参考，助益未来实验室与天文探测。

论文解读：星际C_nHN含氮碳链的电子结构与光谱性质研究

研究背景与问题提出

在星际介质（ISM）中，含氮碳链分子参与复杂有机分子（COMs）形成网络，并可能是氨基酸等前生物分子的前体。C_nHN家族涵盖氰基（–C≡N）、异氰基（–N≡C）、卡宾与亚胺-累积烯等多种结构 motif。尽管HCN、HNC、HCCN及系列氰多炔（cyanopolyynes）已在金牛座分子云TMC-1等源中被探测，但既往理论多聚焦单一家族（如腈类或异腈类），缺乏跨构型异构体的统一高阶基准比较。同时，高丰度异构体分布不能仅由热力学稳定性解释，需结合形成路径与光谱可探测性。因此，开展系统量子化学表征，提供可靠旋转常数、偶极矩与振动频率，对指导实验室与射电天文搜寻具有关键意义。

研究人员工作与核心结论

研究人员通过SmilX（基于TokenSMILES语法框架的分子生成器）枚举C_nHN（n=2–5）构型空间，经DFT预筛后选取21个低能单重/三重态物种，利用显相关波函数方法CCSD(T)-F12（n≤4）与CCSD(T)/def2-TZVP（n=5）计算几何与能量，并扩展计算电荷、偶极矩、垂直电离势与振转光谱。结果表明：奇碳数家族全局极小为线型闭壳层腈结构，偶碳数家族倾向弯曲卡宾或环状单重态（C₄HN全局极小为三元环链结构4-1）；首次理论表征环状C₅HN异构体5-5；偶极矩随同系链长增大，VIP受端基与局域成键主导；计算旋转常数与实验偏差仅0.02%–2.15%，证明协议可外推至未探测异构体。

主要关键技术方法

研究以SmilX语法驱动生成848个初始结构，ωB97X-D/def2-TZVP做DFT预优化与频次验证；对低能子集采用ORCA中的CCSD(T)-F12与CCSD(T)-F12/RI做几何重优化与能量评估；偶极矩由CCSD密度给出；垂直电离势用IP-EOM-CCSD/aug-cc-pVTZ；振动非谐性采用混合GVPT2（revDSD-PBEP86-D3结合B3LYP-D3 QFF）；多参考特征通过T₁、S诊断与CASSCF-NEVPT2交叉检验。无生物样本队列，均为理论构型集。

研究结果分节浓缩

3.1 电子结构与自旋态评估：通过T₁与S诊断、〈S²〉及CASSCF-NEVPT2分析，确认低能三重态2-1、2-5、4-2虽具轻微自旋污染但仍以单组态为主；NEVPT2复现CCSD(T)自旋序：C₂HN三重态更低，C₄HN单重态环状更稳。

3.2 结构参数与相对能量：21个异构体中17个开链、4个含三元环；奇n全局极小为聚炔线型，偶n取弯曲/环化避规价键张力；4-1比线型三重4-2低约3.2 kcal mol^?1；5-5为首次刻画环丙二烯基接–CCN结构。

3.3 电荷分布与偶极矩：CM5电荷显示N恒负、H恒正，极化强度由局域键型主导；亚胺-累积烯卡宾3-4、5-6电荷分离最强，偶极矩分别达6.22 D与8.58 D；同系中偶极矩随链长上升（如4-1至5-5由3.44升至4.07 D）。

3.4 垂直电离势：VIP跨度9.635–11.961 eV；最低见于5-6、4-3、4-2；前线π轨道离域与轨道混合比链长更能解释VIP趋势，卡宾位点局域化降低电离能。

3.5 旋转常数：平衡常数A_e/B_e/C_e与实验可比，线型分子B偏差<0.6%，4-1环结构全常数偏差<0.5%；准线型2-1的A常数因实验拟合受限与大幅振动偏离较大；GVPT2对刚体改善有限，对准线型大振幅运动不足。

3.6 谐波与非谐振动频率：非谐处理总体优于谐波；C–N伸缩频率与Mayer键级正相关（r=0.77）、与键长负相关；准线型与卡宾物种因低频大振幅运动致VPT2失效，需变分法。

讨论与结论翻译

讨论指出，相对能量仅提供热力学参照，星际丰度还受形成/破坏动力学控制；高能在位异构体（如5-6）的探测证明非平衡路径的重要性。旋转常数与偶极矩分量限定可观测跃迁，而形成效率决定实际可探测性。

结论部分：本研究结合语法驱动异构体生成与显相关CCSD(T)-F12，建立n=2–5中性C_nHN共21个物种的一致基准数据集。C₄HN全局极小为环状单重态，C₅HN新环状异构体5-5具光谱潜力。偶极矩随链长增，VIP由局域成键与端基连接主导。计算旋转常数与实验吻合，支持该协议可用于预测含氮碳链光谱特征，为实验室表征与星际搜寻提供关键参考。