由于在光信号存储、光通信、光计算、生物成像、信号处理和图像重建技术等领域应用广泛，非线性光学材料发展得非常迅速^[1-6]，如何获得高效的非线性光学材料成为研究热点。D-π-A型分子的非线性光学性质可通过改变给体、受体的能力以及改变共轭桥的性质进行调节^[7-9]，D代表电子给体(Donor)、π代表含有共轭π键的桥体、A代表电子受体(Acceptor)。这些D-π-A型分子在设计高效非线性光学(Nonlinear Optics, NLO)材料方面具有很大的潜力，因此，大量研究工作聚焦于如何构建D-π-A型分子并有效调节非线性光学响应。

分子开关是重要的纳米组件，由于其在通信、数字处理、金属检测等方面的潜在应用引起了研究人员极大的兴趣。通过外加刺激作用(如：电场、pH、温度、溶剂等)，分子体系由对称结构向不对称结构转变，引起分子结构的扭曲、电子结构的破坏，驱动电荷转移，从而产生“开关效应”^[10-11]，实现分子“开”与“关”状态的快速可逆转换，形成NLO分子开关。NLO开关要求分子必须具备2种NLO响应性质截然不同的形式，并且都很稳定。设计合成这类具有高非线性开关效应的有机分子已经引起了很多实验工作者的兴趣。理论计算方法也被用来研究非线性开关效应且有足够的能力解释其非线性光学性质，甚至可以在实验研究之前就设计出具有显著非线性光学开关效应的有机分子。例如通过给对称分子施加电场的方法，发现了一系列分子开关的非线性光学效应，所以高度对称的分子被认为是一类潜在的非线性开关材料^[12-15]。

高度对称环蕃分子的化学名称为“(1, 2, 4, 5)cyclophane”，具有高稳定性的π电子共轭体系，已经被SAKAMOTO等^[16]在实验上合成，由于其高度对称的结构使其成为潜在的NLO开关分子，但该分子在外电场作用下的非线性光学性质的研究尚未被报道。此外，给环蕃分子引入吸电子取代基和供电子取代基，从而得到不同的环蕃衍生物，由此猜想环蕃分子及其衍生物在一定强度范围的外电场作用下有可能展现出可逆的非线性光学开关效应。那么在多大的电场强度范围内，能设计出有效稳定的非线性光学开关？它具有多大的非线性光学响应强度？随着外电场强度的变化，分子结构和非线性光学性质如何变化？环蕃分子与其衍生物相比，哪种分子具有更优良的“开关效应”？为了探究解答这一系列问题，该研究使用基于密度泛函理论的量子化学计算方法研究了环蕃分子及其衍生物的非线性光学性质。

1.   计算方法

密度泛函理论(DFT)方法的计算准确度主要取决于交换相关泛函的选择，采用新发展的杂化型泛函M06-2X，该泛函被推荐用于主族热化学、非键作用力、价态或里德堡态的激发能、以及较大分子的电子结构的计算^[17-18]，具有较广泛的应用范围。M06-2X适合于本文研究对象分子中π-π相互作用的计算，

此外还采用了D3方法来校正色散相互作用。选择了较大的基组以尽可能的消除基组误差，所有计算均在M06-2X-D3/6-311++G(2df, 2pd)水平下进行。

假设当分子处在一个较弱的外电场时，它的总能量可以展开为：

其中，E⁰表示无外电场时分子的总能量，U_i(U_j或U_k)表示外电场U在i(j或k)方向上的分量，μ、α和β分别表示偶极距、极化率和超极化率。

对于1个分子，平均偶极距(μ₀)和极化率(α₀)的定义如下：

第一超极化率(β₀)定义为：

其中，

本文所有计算使用Gaussian 16程序包^[19]完成，光学响应性质通过求解耦合Hartree-Fock(CPHF)方程来获得，分子结构图通过其附带GaussView软件绘制。

2.   结果与讨论

2.1   分子几何结构

在M06-2X-D3/6-311++G(2df, 2pd)水平上进行几何结构优化和振动频率的计算，优化得到的分子构型在不同电场下均未出现虚频，证明计算所得结构均为分子势能面上的极小点稳定结构。环蕃分子(图 1A)是在实验上已被合成出来的分子，它包含2个平行苯环并且其中4个氢原子被烷基取代，构成了一个立体环状结构体系，具有堆叠的π-π弱相互作用。考虑到分子的空间位阻效应，本文选择了2个体积较小的取代基团，即吸电子基(-Cl)和给电子基(-NH₂)取代苯环上剩余的4个氢原子形成D-π-A型结构，并构造出具有中心反演对称性的环蕃分子衍生物(图 1B)。

图  1  环蕃分子及其衍生物几何结构

Figure  1.  The geometrical structure of cyclophane and its derivative

下载: 全尺寸图片 幻灯片

2.2   外电场强度范围的确定

沿着环蕃及其衍生物分子的z轴方向(图 2)固定电场方向，当外电场强度U与分子前线轨道能级差E_g相近时，分子会产生共振并产生不稳定的结构以及发散的NLO，因此需要避免该情况的发生。本文计算了不同外场强度下分子的E_g并将结果绘制成图 2B。当外电场强度U=14.40×10⁹ V/m时，经单位换算为0.76 eV，对应的环蕃分子能级差E_g为0.96 eV，而当U=14.90×10⁹ V/m时，经单位换算为0.79 eV，其对应的E_g为0.74 eV，因此可确定产生共振的U范围应在14.40×10⁹~14.90×10⁹ V/m。此外，计算了第一超极化率沿z分量β_z随U的变化(表 1、表 2)，发现当U>14.50×10⁹ V/m后，β_z变化非常大，可认为分子在该电场下已经发生共振，所以将环蕃分子外电场强度的研究范围确定为0~14.40×10⁹ V/m。此外，对于环蕃衍生物分子，当U>13.90×10⁹ V/m后，分子结构严重扭曲，不再具有研究意义，因此将其研究电场强度范围确定为0~13.90×10⁹ V/m。

图  2  环蕃分子及其衍生物的前线轨道能级差E_g随外电场强度的变化

Figure  2.  The E_g of cyclophane and its derivative with respect to the external electric field

下载: 全尺寸图片 幻灯片

表  1  环蕃分子的计算结果

Table  1.  The calculation results of cyclophane

电场强度/(×109 V·m-1)	μ0/D	α0/(×10-25esu)	β0/(×10-30esu)	βz/(×10-30esu)
0	0	371.9	0	0
2.571	2.618	371.9	2.765	2.765
5.142	5.237	373.4	5.901	5.901
7.713	7.906	376.4	10.09	10.09
10.28	10.63	379.4	18.06	18.06
12.86	13.45	391.2	85.67	85.67
13.88	14.67	425.3	972.8	972.8
13.99	14.79	438.6	1 397	1 397
14.09	14.92	457.9	2 030	2 030
14.19	15.07	486.1	3 042	3 042
14.30	15.20	532.0	5 204	5 201
14.40	15.35	628.3	14 977	14 328
14.45	15.46	869.9	412 204	110 827
14.50	15.53	172 6	1 865 960	719 489
注：为了计算方便，偶极矩(μ0)的单位采用德拜(1 D=3.334×10-30 C·m)，极化率(α0)、第一超极化率(β0)和第一超极化率z分量(βz)的单位esu表示“静电单位”。下表同。

下载: 导出CSV 

| 显示表格

表  2  环蕃衍生物分子的计算结果

Table  2.  The calculation results of cyclophane derivative

电场强度/(×109 V·m-1)	μ0/D	α0/(×10-25esu)	β0/(×10-30esu)	βz/(×10-30esu)
0	0	435.7	0	0
2.571	3.127	435.7	3.663	3.637
5.142	8.262	435.7	7.231	7.136
7.713	11.41	438.6	12.54	12.47
10.28	14.62	444.6	24.23	24.21
12.86	18.02	477.2	834.5	807.5
13.37	18.73	480.1	969.9	954.5
13.88	32.00	613.5	4 784	4 498

下载: 导出CSV 

| 显示表格

2.3   非线性光学性质

根据表 1和表 2中环蕃分子及其衍生物的特征NLO响应数据，绘制β_z随外电场强度的变化(图 3)。从表 1可知，第一超极化率β₀的贡献主要来自其z分量β_z，因为外电场的方向为垂直2个苯环平面的方向(即z轴方向)。环蕃分子的β_z随外电场的变化表现为：在低电场范围内，β_z缓慢增大，在高电场范围内，β_z迅速增大，表明外电场大小对环蕃分子的非线性光学响应有很大影响，它会大幅度破坏分子的对称性(分子点群由C_2h变为C₂)，并使电子结构发生改变，驱动π→π^*电荷转移，根据含时密度泛函理论，d→π和π→π^*的电荷转移对β₀有很大贡献，从而产生非线性光学响应。

图  3  环蕃分子及其衍生物的β_z随外电场强度的变化

Figure  3.  The change of β_z of cyclophane and its derivative with the strength of external electric field

下载: 全尺寸图片 幻灯片

环蕃衍生物分子的β_z在低电场范围内缓慢增大，对比二者在外电场下β_z的变化趋势可看出，环蕃及其衍生物分子在同一电场下的β_z非常接近，而且二者在低电场区间的变化趋势也基本一致。但在高电场范围内，环蕃衍生物分子稳定性较差，其外部结构趋向散开，因而环蕃分子比其衍生物的稳定性更强，具有更显著的NLO效应。

β₀的增加与“开”有关，β₀的减少与“关”有关。为了获得开关效应，β_z在“开”与“关”状态下的差别必须足够大。NLO响应性质的计算结果表明：环蕃分子及其衍生物分子在一定电场范围下可以实现分子的“开”与“关”，且二者之间可以快速可逆转换，“开”状态下的NLO响应强度大，“关”状态的NLO响应强度为零，分子结构的稳定性较高，可作为新型高性能NLO开关的候选材料。

2.4   几何构型分析

使用苯环之间的C—C平均键长L来测量苯环的面间距(图 4A)，L随外电场强度U的变化如图 4B所示。当U=0~14.4×10⁹ V/m时，环蕃分子的L随U的增加而略微增大(变化幅度仅在4 pm内)，而当U=14.4×10⁹~15.0×10⁹ V/m时，L由增大转变为减小，原因是分子性质在共振前后的趋势相反。相较于环蕃分子，其衍生物分子的结构参数受外电场影响的程度较大。由于当U=13.9×10⁹ V/m时分子结构散开，因此将其衍生物的电场研究范围确定为0~13.9×10⁹ V/m，在该电场范围内环蕃衍生物的结构参数随U的变化具体表现为：L随U增加而逐渐增大(变化幅度在6 pm内)，在U=4.1×10⁹~4.6×10⁹ V/m的平均键长L仅发生接近0.6 pm的突变。

图  4  苯环间C—C平均键长(L)随外电场强度(U)的变化

Figure  4.  The change of the average bond length (L) of the C—C bond with the strength of the external electric field (U)

下载: 全尺寸图片 幻灯片

检查分子在各个电场下的几何结构后发现，当U=4.6×10⁹ V/m时，环蕃衍生物分子中-NH₂取代基的空间朝向发生了翻转(图 4C)，因此，当U=4.6×10⁹ V/m时，分子的几何结构未发生较大变化，L突变的原因是分子由原来的准中心对称结构变为不对称，并导致了E_g的不连续变化(图 2B)和NLO的不连续变化(图 3B)。总体来看，相较于环蕃分子而言，其衍生物分子的结构参数受外电场影响的程度略大，但变化幅度都很小，分子仍能稳定存在，在外电场作用下仍可以产生较大的二阶NLO光学响应。环蕃及其衍生物分子在施加一定范围的外电场后，构型由原来的中心对称结构变成了不对称，即对称结构被破坏，使得分子在某一电场内具有大的二阶非线性光学响应，可以成为非线性光学开关的“开”状态，无外场下的中心对称结构的NLO响应强度为0，从而形成一个较大差值的“关”状态，因此可以作为非线性光学开关材料。

3.   结论

在密度泛函M06-2X-D3/6-311++G(2df, 2pd) 理论水平下，研究了环蕃分子及其衍生物的几何结构及其在外电场作用下的非线性光学性质。研究表明：环蕃分子及其衍生物在适合的外电场强度范围内均可以成为稳定的非线性光学开关。当无外电场时，二者均为中心对称结构，第一超极化率为零，当外电场强度增大但不发生共振破坏分子结构时，二者均具备较大的第一超极化率(8.64×10^-30 esu)。在一定电场强度范围，环蕃及其衍生物分子可实现从中心对称到非中心对称的可逆变换，分子结构会发生轻微扭曲，但总体保持稳定，受外电场的影响很小。分析表明：环蕃及其衍生物分子在一定电场强度范围内都表现出较大的NLO响应，外电场可以驱动分子中的π→π^*电荷转移，而根据含时密度泛函理论，d→π和π→π^*的电荷转移对第一超极化率有很大贡献，从而产生了“开关效应”，并实现“开”与“关”状态之间的快速可逆转换，而且“开”状态的NLO响应很大，“关”状态的NLO响应为零。在较高的外电场条件下，环蕃分子比其衍生物分子的稳定性能更强，具有更显著的开关效应。D-π-A结构并没有显著增加环蕃分子衍生物的NLO响应强度，原因可能是环蕃分子独特的双层π电子结构比单一的π共轭体系更稳定且不易极化，导致D-π-A结构不能起增强环蕃分子NLO响应强度的作用，当然也可能是取代基不合适或取代位置不合适。基于环蕃及其衍生物分子的高稳定性和优异的非线性光学性质，该分子在适宜范围的外电场强度作用下将是新型高性能NLO开关的候选材料。