逐步量子化学优化揭示多并苯层状鱼骨状堆积与多晶型的结构起源

《Communications Chemistry》：Stepwise quantum chemical optimization reveals the origin of layered herringbone packing and polymorphism in polyacenes

　　

多并苯（polyacenes）是典型的π共轭体系（π-conjugated systems），在有机半导体（organic semiconductors, OSCs）中具有重要应用。该系列表现出四种不同的层状堆积基元（layered packing motifs），显著影响半导体性质，但其结构起源仍不清楚。在此，研究人员表明，对具有最小且合理变量集的结构模型进行逐步量子化学优化（stepwise quantum-chemical optimization），能够系统重现所有多晶型（polymorphs）。该模型从由分子对组成的单分子层出发，其中短轴在滑移对称性（glide symmetry）下被任意定位和取向，而长轴保持平行且不发生位移。通过引入长轴倾斜（long-axis inclinations）、扭转（torsions）和层间位移（interlayer displacements），研究人员获得了局部稳定结构，重现了所有多晶型的几何结构和相对稳定性。势能图（potential energy map）分析揭示了各形式特有的稳定机制。这些发现阐明了多并苯多晶型的起源，并为分子堆积如何影响电荷传输（charge transport）提供了见解。该工作建立了一个预测性第一性原理（first-principles）框架，用于理解和设计功能性晶态分子材料（functional crystalline molecular materials）。

多并苯（polyacenes；线性稠合芳香烃，n=2至6）是典型扩展π共轭体系，在有机半导体（OSCs）中应用重要，尤其并五苯（pentacene, n=5）在有机场效应晶体管（OFETs）中表现突出。其晶体具有层状结晶性，层状鱼骨状堆积（layered herringbone, LHB；相邻平面棒状分子形成T形和滑移平行接触）是常见基元。实验上多并苯晶体分为四型：Type I为长轴近垂直层面的简单LHB；Type II为长轴倾斜的LHB；Type III为相邻分子长轴扭转但保留滑移对称性（glide symmetry）；Type IV为Type II与III之间的中间相。问题在于LHB堆积稳定性和多晶型（polymorphs）起源不清楚；成对相互作用分析与实验多晶型相对稳定性连接不直接；量子化学虽可计算相互作用，但难解释为何某堆积有利；从头晶体结构预测维度高、局部极小多且缺机理解释。为此，研究人员以多并苯为模型，阐明稳定机制及随稠合环数变化的堆积差异。研究人员采用基于色散校正密度泛函理论（DFT-D）的逐步优化：从滑移对称下单分子层分子对开始，优化参数得到规则R形式；再依次引入长轴倾斜、扭转和层间位移，重现四型并解释相对稳定性；势能图揭示各形式稳定机制，并计算转移积分（transfer integrals）联系电荷传输。该工作建立预测性第一性原理框架，连接分解相互作用与堆积/多晶型，为设计功能晶态分子材料提供依据。

关键技术方法：以van der Waals（vdW）接触模型做初始几何搜索；用DFT-D（B3LYP/6-311G**、Grimme D3、counterpoise校正基组重叠误差BSSE）计算分子间相互作用能；用SAPT(DFT)（PBE0-GRAC，Psi4，aug-cc-pVDZ）分解静电、诱导、色散和交换排斥；在B3LYP/6-31G*下由最高占据分子轨道（HOMO）重叠计算转移积分；用自动工作流和hill-climbing算法逐步优化参数（α、a、b、φ_incl、θ_incl、θ_twist、φ'_incl、θ'_incl、x、y、z），并模拟粉末衍射。实验晶体结构取自Cambridge Structural Database（CSD；CCDC 1216818、1103069、114446、1844642、114447、665900、619982、847357）。

结果部分：
规则层内分子排列的优化（Optimization of regular intralayer molecular arrangements）：通过以分子对为重复单元、滑移对称构建单层，优化α、a、b并用E_intra(8)评价，发现α≈25°/65°的HB最稳定，2α≈51°接近实验θ_HB≈53°，称R形式，对应Type-I层内堆积；π堆积（π-stacked, PS）和交叉网格（crosshatch, CH）为亚稳，PS因T形接触少而不如HB稳定；SAPT(DFT)表明色散和静电共同稳定HB，色散主导，交换排斥平衡距离；E_intra(8)/n随n增加，诱导项非线性增加。

通过均匀长轴倾斜的优化（Optimization by uniform long-axis inclination）：在R形式α=25°下同时滑移长轴，以θ_incl、φ_incl扫描E_intra(6)，得到R全局最小、N局部最小（θ_incl≈25.4°、φ_incl≈47.6°，无滑移对称）、G局部最小（θ_incl≈32.6°、φ_incl=0°/180°，保留滑移对称）。N对应Type II/IV，G对应Type III；短多并苯N不稳定，解释Type III见于萘/蒽；长多并苯θ_incl随n变化导致Type II/IV竞争。

层间堆叠排列的优化（Optimization of interlayer stacking arrangements）：通过搜索层间横向位移(x,y)和层间距，计算晶格体积V(x,y)与E_inter(7)。R形式最优在(±1.8 ?,0)附近，接近Type-I但对称性有差异；G形式（萘）最优(1.9 ?,0)，后续称twisted G；N形式（并五苯）有两个局部极小N1(0,1.4 ?)和N2(1.4,-0.7 ?)，分别对应Type II和Type IV，能量差小；N1/N2相对稳定性随θ_incl和n变化，解释Type II多见于n=5,6，Type IV见于n=4,5。

滑移对称性下的长轴扭转（Long-axis twist under glide symmetry）：通过G形式单层引入θ_twist和Δz_T，保持滑移对称并最小化E_intra(6)；萘θ_twist≈14°、蒽≈3°能量下降，长多并苯无下降；进一步优化E_int(near)得萘13°、蒽9°，接近实验；扭转增强额外CH/π接触和层间稳定，解释Type III只见短多并苯。

非均匀长轴倾斜的进一步优化（Further optimization with non-uniform long-axis inclination）：通过保持反演对称，仅倾斜中心分子θ'_incl、φ'_incl，对并四苯N型单层得局部最小θ'_incl≈1.6°、φ'_incl≈252°；优化E_int(near)发现Type-II无明显稳定，Type-IV在θ'_incl≈2.5°明显稳定；并四苯Type-IV较Type-II稳定，并五苯/并六苯Type-II仍更稳定，符合实验。

与实验结构的比较（Comparison with experimental structures）：通过比较优化结构与实验结构参数及模拟粉末衍射，总体衍射角和相对峰强一致；Type III（萘、蒽）和Type II（并五苯、并六苯）吻合好；Type IV和Type I高角度有差异，可能因面内晶格常数略小及0 K计算与有限温实验差异；Type IV并五苯为高温相，受热运动/热膨胀影响。

实际并五苯中R形式的偏离（Deviation of the R-form in actual pentacene）：通过引入长轴微倾斜并保持空间反演对称优化，发现E_intra(6)增大、E_inter(7)略降但E_int(near)均增大，Type-I附近无局部能量最小；说明R形式亚稳，实验Type-I受基底/界面、薄膜生长和测量条件影响，只能在表面受限条件实现。

多晶型能量竞争与关键结构特征（Polymorph energy competition and key structural features）：通过考虑更远分子组E_int(G2)至E_int(G5)并外推E_int(all)，发现长程相互作用使Type II、III、IV可低于Type-I；最稳定为萘/蒽Type III、并四苯Type IV、并五苯/并六苯Type II，与实验多晶型倾向一致。总结：Type-I为R形式单层最稳定但层间弱、长程后相对不稳，实验结构因基底/界面偏离；Type-II为N形式、T形接触错位一环、层间强，利于大n；Type-III为G形式、T形接触均匀错位一环并扭转，利于小n；Type-IV为N2型、长轴由平行到略偏移，层间增强，利于并四苯。

转移积分的计算（Calculation of transfer integrals）：通过HOMO重叠计算各优化结构T形和滑移平行接触的转移积分。并五苯Type-I最大，体现二维电荷传输；Type-II/IV较低；T形接触无长轴位移时最大，位移半环显著降低，整环恢复；Type-II因θ_incl接近一个苯环错位而保持较大；说明分子/晶体设计需控制长轴错位以优化电荷传输。

讨论部分总结：滑移对称性减少独立位移和自由度，利于构建分子对并解释高对称短多并苯；构型参数选择关键，先二维排列得R，再均匀倾斜得G/N，扭转和非均匀倾斜复现Type III/IV；方法目前适用于刚性棒状π共轭分子，未来可扩展至其他多环芳烃、氮杂并苯、氟化并苯及含柔性侧链/互锁层状体系；Type-I的R形式单层最稳定但层间弱，长程作用后不如其他型，单晶未确认，仅薄膜观察到，受基底/界面影响。研究结论（译意）：逐步量子化学优化以最小且合理的变量集系统重现多并苯所有四种实验多晶型及其随环数的相对稳定性趋势；层内排列先定、层间后优化的策略有效；长轴倾斜产生G和N局部稳定形式，滑移对称性、扭转和非均匀倾斜分别解释Type III和Type IV；长程相互作用对相对多晶型稳定性重要；该框架为理解层状分子晶体结构起源提供物理透明路径。