耦合簇理论计算的β-丙内酯衍生物的振动拉曼光学活性

《The Journal of Physical Chemistry A》：The Vibrational Raman Optical Activity of β-Propiolactone Derivatives Computed with Coupled-Cluster Theory

　　

振动拉曼光学活性（VROA）是一种信息丰富的技术，可用于探索手性分子的分子立体化学。模拟VROA光谱所需的分子性质可使用耦合簇（CC）理论计算。研究人员首次以原点无关的方式报道了使用CC理论计算的VROA光谱，并将该方法应用于β-丙内酯（氧杂环丁烷-2-酮）的若干衍生物——一类具有天然存在衍生物且呈现多种立体化学的张力四元杂环。光谱在若干方面进行了理论基准测试：首先，比较使用CC2理论计算的VROA光谱以评估大量单电子基组的性能；其次，比较CC2与CCSD计算的光谱以衡量成本较低的CC2方法的有效性；第三，比较基于不同几何结构与振动Hessian矩阵由CC2计算的VROA光谱；最后，使用B3LYP/aug-cc-pVTZ计算的几何结构与Hessian矩阵，以CC2生成了甲基、乙基、丙基取代β-丙内酯的VROA光谱。对于3位（邻近羰基）单取代结构，≈1130 cm–1处出现强峰，主要对应环与烷基取代基间的C–C伸缩；对于4位单取代结构，≈1115 cm–1处出现类似峰，主要对应环与取代基间的C–C伸缩；对于(3S,4R)立体化学双取代物种，涉及环呼吸运动及C3与C4上氢原子同侧不对称摇摆的振动模式在≈1130 cm–1产生大负信号。整体光谱受3位构型主导，对4位较不敏感。结果表明VROA可有效用于探测3位立体化学，但双取代物种中未识别出表征4位立体化学的明确VROA信号。

研究背景方面，振动拉曼光学活性（VROA）因能在远共振近似下通过电偶极极化率α_ij、电偶极/磁偶极极化率G′_ij及电偶极/电四极极化率A_ijk的几何导数提供丰富立体结构信息，成为表征手性中心构型的有力手段。此前VROA光谱多依赖密度泛函理论（DFT）计算，而耦合簇（CC）理论虽具更高精度却鲜见原点无关的系统基准研究。β-丙内酯（氧杂环丁烷-2-酮）衍生物作为高张力四元杂环，具生物活性与合成价值，其3位与4位立体化学的判别长期缺乏可靠光谱标记。研究人员因此在《The Journal of Physical Chemistry A》发表工作，首次以原点无关的CC线性响应理论计算此类分子的VROA光谱，系统评估基组、CC2与CCSD层级、几何/Hessian来源的影响，并探寻可迁移的构型特征峰。

关键技术方法上，研究人员开发Python模块Roa衔接Psi4与CFOUR：几何优化与Hessian由B3LYP/aug-cc-pVTZ或CCSD(T)/cc-pVTZ给出，手性张量由频率依赖CC2及CCSD线性响应（采用修正速度规范电偶极算符，故原点无关）在Psi4中计算，张量对简正模导数通过核位移有限差分获得；基组涵盖cc-pVDZ至aug-cc-pVQZ及LPol、ORP等系列；光谱以532 nm入射、SCP背向180°模式按Δ(180°)公式绘出；谱间相似度用SNO、DNO、IDF、RADF、MJAC五种度量评估；样本为(S)-甲基环氧乙烷及3位/4位单取代与(3S,4R)双取代甲基/乙基/丙基β-丙内酯。

研究结果部分，“电偶极/电四极项”小节指出，CC计算中β(A)²项在多数模式与β(G′)²同号，背向180°实验 neglect 该四极项不致严重误差，但Δ(0°)或Δ(90°)_z配置下其符号反转关键。“(S)-甲基环氧乙烷的单电子基组”小节表明，固定几何/Hessian时各基组与aug-cc-pVTZ参考谱DNO均>0.99，asp-cc-pVDZ已可用；全自洽计算时aug-cc-pVDZ因频率复现优而优于ORP与Sadlej小基组，LPol大基组遇线性依赖。“CC2与CCSD比较”及“几何/Hessian影响”小节显示，CC2/aug-cc-pVDZ与CCSD/aug-cc-pVDZ在(S)-甲基环氧乙烷上形状吻合实验气相谱，B3LYP几何与CCSD(T)几何差异仅微调峰位；β-丙内酯体系中CC2可替代CCSD定性寻峰。“取代β-丙内酯VROA特征”小节总结：3位单取代在≈1130 cm^–1有强C–C伸缩VROA峰，4位单取代移至≈1115 cm^–1，(3S,4R)双取代在≈1130 cm^–1因环呼吸+C3/C4氢同侧摇摆呈大负峰；整体谱形由3位构型主宰，4位无稳健标记。

讨论与结论翻译：研究人员判定VROA光谱整体受3位构型主导且对4位不敏感，CC2配合B3LYP/aug-cc-pVTZ几何与Hessian可低成本给出定性正确谱形；asp-cc-pVDZ级基组已足用于张量导数，但全CC流程需aug-cc-pVDZ级基组以保频率；β-丙内酯环刚性使构象扰动有限，故≈1130 cm^–1附近3位特征峰有望推广至更大天然类似物以判识立体化学，而4位构型仍需其他手段。该工作为CC理论原点无关VROA计算建立基准，并提供手性杂环光谱判据。