分子超极化率作为铁电向列液晶二阶非线性光学筛选描述符

《Physical Chemistry Chemical Physics》：Molecular hyperpolarisability as a screening descriptor for second-order nonlinear optics in ferroelectric nematic liquid crystalsOpen Access

　　

铁电向列（NF）液晶兼具流动性与宏观极性有序性。尽管原型NF材料RM734展现出巨大的非线性光学系数，但大多数已知铁电向列材料在设计时并未考虑光学非线性。在此，研究人员评估了电子结构计算能否用于在已知极性液晶材料中识别有前景的非线性光学候选物。研究人员采用一系列DFT方法和基组计算了频率依赖的分子超极化率张量（βijk），并利用取向气体模型（oriented-gas model）结合经验〈P1〉和〈P3〉值将其转换为宏观d系数。计算结果与代表性材料可获得的实验d33、d15和d13/d31系数进行了基准对比。研究人员发现，绝对值强烈依赖于方法和体参数假设，而相对趋势对筛选更为有用。显式构象平均并未持续改善与实验的一致性。将性能最优的单构象方案应用于更广泛的极性液晶数据集，识别出具有增强预测非线性光学响应的候选材料，并揭示了基于供体-受体不对称性、共轭长度和连接基团选择的简单设计规则。

论文解读：分子超极化率作为铁电向列液晶二阶非线性光学筛选描述符

一、研究背景与科学问题

二阶非线性光学（NLO）现象源于材料极化对施加电场的二次响应，其宏观极化可展开为线性、二阶和三阶电极化率项。二阶项负责倍频（SHG）、自发参量下转换等过程，在光子学和量子光学中具有重要应用。二阶非线性光学过程要求材料具有非中心对称结构；在常规向列相中，指向矢是无极性的，+n≡?n，因此该相在平均意义上具有中心对称性。近年来发现的铁电向列（N_F）液晶相及类似的极性近晶相（如SmA_F、SmC_P）提供了一条在软流体介质中实现非中心对称组织的途径，因而成为二阶NLO材料的理想候选。

铁电液晶可通过小电场重新取向，且极性有序自发形成，无需外加极化电场，这使其相较于固态材料和聚合物体系具有显著优势。原型N_F材料RM734在N_F相中展现出巨大的非线性光学系数，但该材料（以及大多数其他N_F材料）并非作为NLO发色团设计。与极性液晶组织兼容且能最大化分子超极化率（β）的结构特征尚不明确，这限制了理性发现新型NLO优化的N_F材料。

二、研究内容与意义

研究人员旨在评估电子结构计算能否作为筛选工具，在已知极性液晶材料中识别有前景的NLO候选物。研究通过计算频率依赖的分子超极化率张量，结合取向气体模型将其转换为宏观d系数，并利用实验数据对多种DFT方法和基组进行基准验证。研究最终确定了最佳计算方案，应用于更广泛的极性液晶数据集，识别出具有增强NLO响应的候选材料，并揭示了基于供体-受体不对称性、共轭长度和连接基团选择的简单设计规则。该研究为理性设计新型NLO优化的铁电向列材料提供了计算筛选框架和分子设计指导。论文发表在《Physical Chemistry Chemical Physics》。

三、关键技术方法

研究人员采用多层次的电子结构计算方法。初始几何结构通过GFN2-xTB半经验紧束缚方法优化；构象通过ORCA 6.1中的GOAT算法生成。频率依赖超极化率张量在800 nm入射波长下，使用Gaussian 16软件包中的多种方法计算，涵盖Hartree–Fock、杂化GGA（B3LYP、BHandHLYP）、杂化meta-GGA（M06、M06-2X、M06-HF、TPSSH）和范围分离杂化（CAM-B3LYP、wB97X-D、LC-BLYP）等，配以6-31+G(d)、6-311++G(d,p)、aug-cc-pVDZ和aug-cc-pVTZ基组。宏观d系数通过取向气体模型计算，采用经验〈P₁〉和〈P₃〉序参数值，各向同性局域场修正采用折射率n_e=1.8、n_o=1.5，密度统一取1.3 g cm^?3。基准验证使用RM734、DIO和C1三种材料的实验数据。筛选数据集来源于文献报道的所有极性液晶材料。

四、研究结果

全局最小能量几何上的计算

对RM734、DIO和C1的基态几何进行单构象计算，并与实验值对比。HF方法对RM734的d₃₃、d₃₁和d₁₅均低估，对DIO则高估。B3LYP在所有分子上表现普遍较差而被排除。BHandHLYP系统性地倾向于产生较大的计算系数。杂化meta-GGA方法表现出对精确交换成分的明显依赖：M06给出最大系数并显著高估实验值，M06-HF给出最小系数且对RM734和C1的d₃₃估计较好，M06-2X介于两者之间。范围分离杂化方法普遍高估RM734的d₃₃，LC-BLYP在该家族中低估。aug-cc-pVDZ与aug-cc-pVTZ基组对计算结果影响不大。

构象系综上的计算

研究人员对RM734和DIO进行构象系综分析。结果显示，将β作为构象系综的加权平均值并未显著改善实验与计算NLO系数之间的相关性。RM734的构象系综导致d₁₅和d₃₁一致增加；DIO在基态和构象预测之间存在较大差异，但整体上构象平均并未带来系统性改进。

方法选择

通过计算d₃₃比值（以RM734为参考）并定义偏差和离散度两个评价指标，研究人员发现所有测试方法均一致高估d₃₃值。范围分离杂化方法一致产生较低的离散度，表明其最有可能成功复现实验比值。基于偏差-离散度分析，研究人员选择wB97X-D/aug-cc-pVDZ方法组合，因其具有最低的离散度，表明在基准数据中能一致地计算比值。

N_F数据集筛选

将wB97X-D/aug-cc-pVDZ方案应用于文献报道的极性液晶数据集。计算的d₃₃系数范围从0到49.7 pm V^?1，RM734计算值为29.7 pm V^?1（实验值22.5 pm V^?1）。通过UMAP投影和化学结构分析，研究人员识别出具有增强NLO响应的候选材料。

设计规则

研究发现，最大的β_xxx和d₃₃值一致出现在具有"推拉"电子体系的系统中。RM734具有2,4-二甲氧基苯甲酸酯作为供体、4-硝基苯基作为受体的推拉体系。使用N,N-二甲氨基（NF1）或甲硫基（S8）等更强供体可产生更大的β_xxx和d系数。2,5-二取代1,3-二氧杂环己烷是NLO性质的"死重"，仅在与其他官能团（如炔基）组合时才产生显著响应。三氧杂双环[2.2.2]辛烷仅在邻近苯乙炔基序时有效。连接基团方面，炔基显著优于酯基，而酯基、硫酯和直接苯-苯键相当。二氟亚甲氧基（Q桥）因其共轭中断而必然导致较小的β值。侧向氟取代虽能增强极性有序的热稳定性，但会降低计算NLO系数。扩展π共轭体系（增加芳香环或使用共轭连接基）与增大NLO系数相关。

五、讨论与结论

宏观NLO系数不仅取决于分子超极化率，还依赖于材料特定参数，包括奇秩序参数、分子数密度和折射率。这些体性质无法从孤立分子的电子结构计算中确定，且预计在不同材料间存在差异。使用RM734的代表性序参数和折射率值归一化了分子超极化率对计算非线性系数的贡献，为排名分子性能提供了一致基础，因此计算系数应解释为相对筛选指标。β_xxx与d₃₃之间存在线性关系，计算本质上是β_xxx的线性重缩放。

研究结论表明，具有大二阶NLO系数的N_F材料的设计规则相当简单：沿分子长轴由扩展、平面、共轭通路连接的强供体-受体不对称性。促进极性组织但中断共轭的基序（如二氟亚甲氧基）在这方面必然不利。该研究为铁电向列液晶材料的NLO性能预测和理性设计提供了有效的计算筛选框架。