删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

烷基环己苯异硫氰酸液晶材料太赫兹波吸收

本站小编 Free考研考试/2021-12-29

摘要:基于密度泛函理论, 利用Gaussian09程序包, 通过Opt + Freq优化结构, 以B3LYP/6-311g基组对4-(trans-4-n-alkylcyclohexyl) isothiocyanatobenzenes (CHBT) 16种液晶分子在0.1—5.0 THz波段的吸收进行研究, 提出了通过考察分子转动惯量和质量重心偏移对太赫兹吸收的影响方法. 计算结果表明, 随着分子烷基链上碳原子数目增加, 分子的转动惯量和重心的偏移将对太赫兹波的吸收产生影响. 在0.1—5.0 THz波段, 3—7个碳原子的分子吸收相对较强. 以此为参照, 减少和增加烷基链上的碳原子数目都会降低分子对太赫兹波的吸收. 在0.3—3.0 THz范围内, 将计算结果与10种有实验数据的分子进行了对比. 结果发现, 低频波段计算结果与实验测量相比存在差异, 其中吸收峰位置的差异可能来源于氢键. 比较吸收强度的相对大小值, 发现实验测量与计算结果较为一致, 表明吸收强度来源于偶极子振动、转动吸收, 展现了计算模拟的积极意义. 研究结果可对相关分子的设计与合成提供有用的建议.
关键词: 液晶/
太赫兹/
烷基链/
高斯

English Abstract


--> --> -->
液晶是既有液体的流动性, 又有晶体的各向异性物质. 大多数液晶是具有长棒状或园盘状分子形状的有机化合物. 由于其半流体性质, 在电场作用下, 液晶分子由于偶矩的响应会按电场方向取向, 导致分子原有的排列方式发生变化, 从而使液晶的光学性质也随之发生改变, 这种电光效应通常都难免伴随着对光吸收的改变[1,2]. 因此, 根据电光器件的性能要求, 液晶物质与光的相互作用机理研究显得十分重要[3,4].
太赫兹技术是目前光学领域最重要的新兴技术之一[5]. 在许多方面, 比如医学[6]、安全监测系统[7]、生命科学等领域都有应用[8,9]. 通常, 太赫兹波是指频率范围为0.1—10.0 THz之间的电磁波, 波长范围为0.03—3.00 mm, 介于微波与远红外波段之间. 由于太赫兹的频率较高, 所以具有良好的空间分辨率和时间分辨率[10,11].
在安全检查和医学检查方面, 由于X光的光子能量高, 对人体伤害大, 而低频率的红外光能量过低, 穿透力差, 不易被探测器接收, 同时, 红外探测器还易受环境的干扰[12]. 处于高频区域的太赫兹波克服了两者的缺点, 光子能量相比X光较低, 只有毫电子伏特, 不易破坏被检测物质[13]. 因此, 太赫兹波是安全检查和医学影像技术的理想波谱. 上述器件中, 液晶材料作为调制光波的必备材料, 其与太赫兹光的相互作用是关键因素[14,15]. 液晶分子的电光效应是建立在分子极性基础之上, 极性越强, 对THz波的响应与吸收都十分敏感. 如果对THz波吸收太强, 波能量的损失就越大. 另外, 液晶分子的振动和旋转也对波吸收有较大的贡献. 所以, 选择一种对THz波具有较小的吸收且敏感的电光响应材料是较为困难的事情. 红外光谱是材料分子吸收某些波长的红外光的谱线. 该波长的光能够引起分子振动能级和转动能级的跃迁, 根据分子内部原子间的相对振动和分子转动等信息可以确定物质的分子结构. 因此, 红外吸收光谱又称分子振动光谱或振转光谱[16].
通常将红外光谱按波长分为三个区域: 近红外区(0.75—2.50 μm)、中红外区(2.5—25.0 μm)和远红外区(25—300 μm). 中红外光谱属于分子的基频振动光谱; 远红外光谱则属于分子的转动光谱和某些基团的振动光谱. 太赫兹波谱是位于红外光谱上的中、远红外区的光波. 因此, 研究太赫兹波的问题, 事实上是研究中、远红外波的工作. 根据相关光谱文献[16], 关于物质对红外光谱吸收的原因基本上都有较为详细的结论. 例如, 有机化合物分子中的官能团特征吸收峰表现在中红外光的较高频段, 而在远红外波段, 产生吸收峰的主要原因是分子的振动与转动对波的吸收. 通常, 有机分子内低能电子跃迁、含氢原子团伸缩振动的倍频吸收位于近红外(大于120 THz)波段, 绝大多数有机化合物和无机离子的基频吸收带, 即分子内与单个化学键相关的伸缩和弯曲等振动模式出现在中红外区(6—120 THz), 而各种有机分子之间的氢键弱相互作用、偶极子的振动和转动跃迁以及晶体中声子的低频骨架振动所对应的吸收频率均位于远红外区(小于6 THz)的太赫兹波段[10].
由于THz波的产生对波源有较高的要求, 而实验研究有较多的不可控制因素, 因此, 关于液晶在THz波段的实验测量数据并不多, 关于THz吸收方面的物理和化学机理分析也较少.
由于理论模拟方法简单便利, 能够有效了解和研究物质分子的THz吸收特性, 工作效率高、成本低, 近些年在液晶THz吸收方面也有一些计算模拟工作[17,18]. 这些工作能够对液晶分子的性质进行设计与预测, 对材料的合成具有指导意义.
本研究参考了具有THz实验测量数据的液晶nCHBT(4-(trans-4-n-alkylcyclohexyl) isothiocy-anatobenzenes) (n = 0—15)系列分子[19], 通过计算得到分子的中远红外吸收光谱, 分析分子在0.1—5.0 THz范围内的波吸收状况, 讨论烷基链的碳原子数目与THz吸收的关联性.
本研究利用基于密度泛函理论的Gaussian09程序包[20]对分子进行模拟, 特别是对分子进行红外光谱方面的计算, 为分子合成与实验测量进行了有效的理论预测[21-26]. 关于液晶分子THz波吸收方面的模拟计算工作也有一定的研究.
图1(a)图1(b)展示了液晶分子5CB[18]和MBBA(N-(4-methoxybenzylidene)-4-butylaniline)[17]在THz波段的实验以及计算结果. 由于实验测量中液晶分子的吸收是由于氢键弱相互作用、偶极子的振动和转动跃迁引起的, 有些振动频率简并或十分接近不易分辨, 还有一些简正振动无红外活性, 使分子不产生红外吸收, 所以吸收谱难以观察到吸收峰. 另一方面, 由于实验误差和噪音的影响, 低频率附近的测量结果总是难以表征. 而理论模拟采用的是单分子模型, 没有考虑分子之间的相互作用, 因此没有体现氢键的作用, 吸收峰值体现了偶极子振动和转动的吸收, 吸收峰值比较稀疏. 但是可以看出计算结果和实验测量的吸收趋势十分相似, 说明计算模拟得到偶极子振动和转动的吸收值与实验测量相比是准确的. 通常, 在较高频段, 液晶分子产生的吸收绝大部分为官能团产生的特征吸收峰, 因而计算结果与实验测量符合得较好. 因此, 使用Gaussian计算模拟结果是可靠的, 可以详细地展示分子对波吸收的复杂性. 这些计算结果为本工作的可靠性提供支撑[27-29].
图 1 (a) 5CB计算值和实验值对比; (b) MBBA计算值和实验值对比
Figure1. (a) Comparison of calculations and experimental results about 5CB; (b) comparison of calculations and experimental results about MBBA.

本文所研究的16种液晶分子, 均是在CHBT分子的基本结构基础上, 通过改变烷基链碳原子的数目n构成不同的液晶分子. 其结构如图2所示.
图 2 16个棒状液晶分子的化学结构 (n = 0—15)
Figure2. Chemical structure of 16 liquid crystal molecules (n = 0?15)

在使用Gaussian09程序包过程中, 以Opt + Freq计算类型, 使用B3LYP/6G311g基组优化了16种液晶分子, 得到了红外吸收光谱等数值[20]. 通过这些数值, 对影响分子吸收的因素进行了讨论. 讨论总体分为两部分:第一部分对0.1—5.0 THz频率范围内分子的吸收进行了分析;第二部分针对Chodorow等[19]的实验数据, 比较了0.3—3.0 THz波段范围内3—12 CHBT液晶分子的吸收光谱.
本文讨论16种向列相液晶分子的THz吸收状况, 根据文献, 只有n = 3—12这10种分子具有实验数据, 其他6种分子暂未发现实验数据. 但是, 作为参考, 我们也对其进行了计算讨论.
2
3.1.0.1—5.0 THz波段的波吸收分析
-->图3展示了0.1—5.0 THz波段16个分子的吸收状况. 可以看出, 所有分子的吸收都有一个较高的峰值, 该吸收峰的分布十分相似, 大约分布在0.4—0.9 THz之间. 但分子峰值的大小有所有不同. 当n = 3, 4, 5, 6, 7时, 分子的吸收峰值较高. 0.1—5.0 THz属于远红外波段, 对于单个分子模型而言, 分子的振动与转动是对波吸收的主要因素. 因此, 这些吸收峰差异是由分子自身振动与转动差异带来的[10,16].
图 3 0—5 THz波段分子吸收光谱
Figure3. Molecular absorption spectroscopy at 0?5 THz.

另一方面, 红外吸收谱的强度是分子振动跃迁概率的表征, 与分子振动时偶极矩的变化相关. 因此, 与基团固有偶极矩有关, 极性越强, 振动时偶极矩变化越大, 吸收谱峰值越高; 分子中的原子对称性好, 振动时偶极矩变化就小, 吸收谱峰值越弱. 液晶分子在低太赫兹波段的吸收同时又是分子转动的一个重要体现. 分子的转动取决于分子内原子质量分布和转动惯量. 因此, 为了分析0.1—5.0 THz波段吸收峰产生的原因, 我们计算了这些分子的偶极矩和转动惯量.
16个分子转动惯量和偶极矩的计算结果随烷基链长度的变化展示在图4图5中. 从图4图5可以看出, 分子的转动惯量随碳原子数目的增加而增大. 偶极矩只有n = 3, 5候有些异常, 可能是烷基链碳氢的增加对分子重心从苯环到环己烷的迁移产生的电荷中心的差异造成的. 其他分子的偶极矩整体差异不大, 但是微小差别显示出偶极矩的奇偶效应. 这说明分子的电荷分布随烷基链的增加变化较小, 显示烷基链对分子的极性贡献不大.
图 4 16个分子在其长轴的转动惯量值的变化趋势
Figure4. Trend of the rotational inertia of 16 molecules on their long axis.

图 5 16个分子在其长轴的总偶极矩值的变化趋势
Figure5. Trend of the total dipole moment of 16 molecules on their long axis.

影响分子振动和转动的因素主要考虑分子的重心和转动惯量, 因此, 需要讨论分子内的原子质量分布(如表1)和转动惯量对吸收的影响. 我们把整个分子分为NCS-Benzene-Cyclohexane-CnH2n+1几个部分. 根据各部分的原子质量之和, 确定分子的重心位置. 表1图6展示了分子质量分布和重心位置的标定. 以n = 1为例, 经过计算发现, 分子的质量中心在靠近环己烷一端的苯环位置.
图 6 1CHBT的分子质量分布和重心位置标定
Figure6. Calibration of mass distribution and center of gravity for 1CHBT.

nCHBTNCS/
g·mol–1
Benzene/
g·mol–1
NCS+ Benzene/
g·mol–1
Cyclohexane/
g·mol–1
NCS+ Benzene +
Cyclohexane/g· mol–1
CnH2n+1/
g·mol–1
Cyclohexane + CnH2n+1/
g·mol–1
0CHBT587613482216183
1CHBT5876134822161597
2CHBT58761348221629111
3CHBT58761348221643125
4CHBT58761348221657139
5CHBT58761348221671153
6CHBT58761348221685167
7CHBT58761348221699181
8CHBT587613482216113195
9CHBT587613482216127209
10CHBT587613482216141223
11CHBT587613482216155237
12CHBT587613482216169251
13CHBT587613482216183265
14CHBT587613482216197279
15CHBT587613482216211293


表116个分子各个组成部分的原子质量
Table1.Atom mass of each component of 16 molecules.

下面具体分析所有分子的质量分布和重心位置与吸收的关系. 从表1可以看出, 当n ≤ 3时, 异硫氰酸脂加苯环的质量大于环己烷加烷基链的质量, 分子的重心在苯环上. 因为苯环的结构坚硬, 所以分子短而且不容易发生弯曲, 振动吸收较弱. 但当n = 3时, 分子的重心十分接近苯环与环己烷相连的C—C键, 导致该C—C键容易弯曲, 分子的振动变得容易, 从而增加3CHBT对光的振动吸收, 所以3CHBT峰值较高, 而(0—2)CHBT峰值较低. 当4 ≤ n ≤ 7时, 异硫氰酸脂加苯环的质量小于环己烷加烷基链的质量, 分子的重心迁移在环己烷上. 因为环己烷的结构柔软, 导致环己烷的C—C—H与连接苯环的C—C—C形成的面角扭曲或者弯曲振动变得容易, 从而增加4CHBT, 5CHBT, 6CHBT和7CHBT对THz的吸收, 所以它们的峰值较高. 当n ≥ 8时, 虽然分子的重心继续迁移, 结构稳定性较差. 但是, 随着烷基链上碳原子的数目不断增加, 分子逐渐变长, 质量逐渐增大, 转动惯量增大变得显著, 相对于同样的温度和分子内能, 分子的转动速度变低, 振动与转动状态变得稳定, 对THz波的吸收变得较弱, 所以吸收峰值较低.
2
3.2.计算与实验结果对比
-->Chodorow等[19]通过实验测量了0.3—3.0 THz波段(n = 3—12)CHBT分子的吸收光谱. 图7对比了10组实验数据与计算结果. 由于实验测量分为o光和e光, 根据相关理论, 我们将o光与e光吸收数据加和绘制为一条曲线. 由于实验测量和计算结果吸收的单位不同, 实验测量的单位是吸收系数, 而计算结果只是吸收的相对数值大小, 实验结果和计算结果的吸收峰值高低的绝对值比较没有意义, 只需要观察不同分子相对峰值的高低. 因此, 为了方便对比和观察, 我们将计算结果的数值进行倍数放大.
图 7 0.3—3.0 THz波段(3—12)CHBT吸收光谱的实验值与计算值对比
Figure7. Comparison of calculations and experimental results about (3?12)CHBT in 0.3?3.0 THz.

通过图7的数据对比, 可以看出, 就像其他计算得到的结果一样[17,18], 在0.3—0.8 THz之间计算得到的强吸收峰, 实验测量都没有得到. 实验测量没有发现明显的吸收峰, 吸收只是随着频率的增加而增加, 但是在1.7 THz之后, 液晶开始出现吸收, 大部分材料都有平坦的吸收峰. 而计算值在0.5 THz处有特别明显的吸收峰, 在之后又有一些较低的吸收峰.
实验测量在低频段峰值很少, 在这个频段, 有机分子之间的氢键(N—H)弱相互作用在吸收中占据主要成分, 而不同液晶材料的分子结构与分子间的氢键是相似的. 偶极子的振动-转动跃迁以及晶体中声子的低频骨架振动所对应的吸收频率占据次要成分, 因此吸收的峰值位置差异难以体现, 而吸收的强度差别体现在偶极子跃迁. 换句话说, 实验测量更取决于氢键和声子的贡献, 由于液晶分子有些振动频率简并或十分接近不易分辨, 还有一些简正振动无红外活性, 使分子不产生红外吸收, 所以从吸收谱难以观察到吸收峰. 另一方面, 理论模拟采用的是在300 K设置温度下的单分子模型, 没有考虑分子之间的相互作用, 因此氢键的作用没有体现. 计算得到的峰值取决于偶极子的振动跃迁, 这个跃迁吸收只对吸收强度有贡献, 而吸收峰值比较稀疏.
为了说明上述现象和便于观察, 图8展示了计算结果在0.5 THz之后的较低峰值与实验值的对比. 可以发现, 对于实验测量中吸收强度较大的((n = 3—7)CHBT)的分子, 计算结果在相应的位置也表现出强的吸收, 表明该强度来源于偶极子振动吸收的贡献(图8(b)中标出). 相反, 其他几个分子具有较低吸收峰. 也就是说, 计算结果的峰值体现的是偶极子振动和转动吸收, 对吸收强度的贡献结果与实验测量体现出来的规律是一致的.
图 8 计算结果与实验测量的吸收峰值分析 (a)实验测量; (b) 计算在1.2—3.0 THz的吸收峰值
Figure8. Analysis of the absorption peaks between the calculation and the experimental data: (a) The experimental measurements; (b) the calculation of absorption peaks in 1.2?3.0 THz.

遗憾的是, 文献中的实验数据仅仅只有0.3—3.0 THz的结果, 不像5CB分子那样具有较宽的波段[18], 因此, 难以对计算值的吸收峰和实验值的趋势进行更宽域的对应. 就像前边讨论的那样, 液晶材料5CB的理论计算也是在较高的频率才能够与实验结果在吸收峰位置有重合, 在低频段, 实验测量也只是得到了线性增加的类似结果, 也说明了实验测量可能是来源于氢键的贡献. MBBA的测量中, TANAKA的数据让人难以置信, 我们已在前边的工作中做过陈述[17]. NISHIZAWA的测量数据也像5CB一样显示出了上述特性, 特别是温度对波吸收的影响也不是很大, 但是零度的测量明显具有较多的吸收峰[17]. 在本研究中, 我们尝试了使用其他计算方法和关键词, 但是得到的结果与本工作使用的参数相差不大, 规律相同. 因此, 实验测量的起始阶段, 是否是仪器设备无法精确探测, 我们不得而知, 只能寄希望于未来有更多的实验数据进行比对.
利用理论模拟的方法研究了16种向列相液晶分子在0.1—5.0 THz波段的吸收情况. 结果发现, 随着分子烷基链上碳原子数目的增加, 分子对太赫兹波的吸收会带来影响. 通过分析分子结构发现, 分子烷基链的长度决定分子质量、重心、转动惯量, 这些不同的力学参数导致分子的转动与振动差异, 最终带来了对太赫兹波吸收的差异. 当碳原子数目小于3或大于7, 分子对太赫兹的吸收较低, 其他分子吸收较高. 通过比较实验测量与计算结果的吸收峰位置, 结果发现, 虽然低频波段计算结果无法体现氢键对吸收的贡献, 但是实验测量的吸收强度的差异可能来源于偶极子振动、转动的吸收, 这个差异能够在计算结果中显示出来, 体现出了计算模拟的积极意义. 期待未来的实验测量结果, 对计算的方法和关键词以及计算温度等参数加以修正.
相关话题/计算 实验 测量 质量 光谱

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 水合肼还原的氧化石墨烯吸附NO<sub>2</sub>的实验研究
    摘要:还原氧化石墨烯由于独特的原子结构,作为气体检测领域有潜力的候选者引起了研究者们的广泛兴趣.本文采用水合肼作为还原剂来制备还原氧化石墨烯,并以此作为叉指电极的气体敏感层,研究了其对NO2气体的响应特性.结果表明,水合肼还原的氧化石墨烯可以实现在室温下对浓度为1—40ppm(1ppm=10–6)的 ...
    本站小编 Free考研考试 2021-12-29
  • 基于1556 nm光纤激光器频率分裂效应的应力测量
    摘要:光学元件在红外波段的应力-光学常数是众多光学系统关心的问题之一.本文提出一种基于1556.16nm掺铒光纤激光频率分裂效应的光学玻璃内应力致双折射测量方法.选择平面介质膜腔镜和光纤光栅(FBG)构成线形半开放式谐振腔,并分析了光纤自身弯曲引入谐振腔内的双折射.将待测光学玻璃附带力传感结构放置在 ...
    本站小编 Free考研考试 2021-12-29
  • 三氨基三硝基苯基高聚物粘结炸药热力学性质的理论计算研究
    摘要:高聚物粘结炸药(PBX)的热力学性质是用于炸药结构响应、安全性评估、数值模拟分析等的重要参数.由于PBX结构的多尺度特性,完全采取实验方法精细表征这些参数存在巨大的挑战.本文运用第一性原理和分子动力学计算的方法,系统研究了三氨基三硝基苯(TATB)基高聚物粘结炸药的热力学参数和界面热传导性质. ...
    本站小编 Free考研考试 2021-12-29
  • 一种K分布强湍流下的测量设备无关量子密钥分发方案
    摘要:研究了K分布强湍流下自由空间测量设备无关量子密钥分发协议模型,采用阈值后选择方法来减少大气湍流对密钥生成率的影响,对比分析了使用阈值后选择方法前后协议的密钥率和湍流强度之间的关系.仿真结果表明,使用阈值后选择方法可以有效地提高协议的密钥生成率,尤其是在高损耗和强湍流区域,而且其最佳阈值与湍流强 ...
    本站小编 Free考研考试 2021-12-29
  • 钫原子磁偶极超精细结构常数及其同位素的磁偶极矩的理论计算
    摘要:应用基于B样条基组的相对论耦合簇理论方法,计算了212Fr原子的nS(n=7—12),nP(n=7—12)和nD(n=6—11)态的磁偶极超精细结构常数.与精确实验值的比较说明这套理论方法能精确计算出磁偶极超精细结构常数,其中7P态的磁偶极超精细常数的理论值与实验值之间的差异小于1%.在忽略场 ...
    本站小编 Free考研考试 2021-12-29
  • 中国散裂中子源反角白光中子束流参数的初步测量
    摘要:中国散裂中子源(CSNS)已于2018年5月建设完工,随后进行了试运行.其中的反角白光中子束线(Back-n)可用于中子核数据测量、中子物理研究和核技术应用等多方面的实验.本文报道对该中子束的品质参数测量实验过程以及最终实验结果.实验主要采用中子飞行时间法,利用235U,238U裂变室和6Li ...
    本站小编 Free考研考试 2021-12-29
  • 大质量样品厚度对中子活化定量分析的影响和修正
    摘要:滑石粉的主要化学成分是Mg3[Si4010](OH)2,通过检测面粉中24Mg和28Si原子核的含量,可以计算出面粉中滑石粉的含量.利用中子活化分析方法测量面粉中24Mg和28Si原子核的含量问题时,被测样品内部中子通量和能量随厚度的变化以及γ射线自吸收效应会对测量结果有较大的影响.利用MCN ...
    本站小编 Free考研考试 2021-12-29
  • 基于共心球透镜的多尺度广域高分辨率计算成像系统设计
    摘要:针对实时广域高分辨率成像需求,充分利用具有对称结构的多层共心球透镜视场大且各轴外视场成像效果一致性好的特点,设计基于共心球透镜的多尺度广域高分辨率计算成像系统.该系统基于计算成像原理,通过构建像差优化函数获得光学系统设计参数,结合球形分布的次级相机阵列进行全局性优化,提高系统性能的同时有效简化 ...
    本站小编 Free考研考试 2021-12-29
  • W波段分布作用速调管的设计和实验研究
    摘要:基于运动学理论、感应电流定理和电荷守恒定律,分析了分布作用谐振腔的渡越时间效应,推导了各个谐振腔工作于${\text{π}}$模的电子注与微波之间的能量转换系数、电子负载电导和电子负载电纳,计算结果显示采用分布作用谐振腔有利于提高速调管的工作效率.利用三维电磁仿真软件,设计了一款工作于W波段的 ...
    本站小编 Free考研考试 2021-12-29
  • 基于可调谐激光吸收光谱技术的二硫化碳中红外光谱参数测量
    摘要:采用$4.6\;{\text{μ}}{m{m}}$附近的量子级联激光器作为光源,搭建了一套二硫化碳(CS2)吸收光谱测量系统,结合可调谐二极管激光吸收光谱技术,对光谱范围为2178.99—2180.79cm–1的CS2吸收光谱展开了深入研究,重点测量了2180.5—2180.74cm–1的 ...
    本站小编 Free考研考试 2021-12-29