作者简介: 汪 源, 1984年生, 三亚学院理工学院教授 e-mail: wangyuan0155@163.com
相干反斯托克斯拉曼光谱(CARS) 和相干反斯托克斯超拉曼光谱(CAHRS)等高阶非线性光谱技术已经应用在动力学过程、 基因表达谱筛选、 高分辨率光谱分析等诸多领域。 但因其涉及到的高阶微观极化率张量元数量众多, 对其光谱信号进行定量分析非常困难。 文献通过理论分析将CARS和CAHRS的微观极化率张量元分解为拉曼微观极化率张量元微分 α'i'j'和超拉曼微观极化率张量元微分 β'i'j'k'的乘积, 将高阶光谱定量分析的困难简化为对低阶光谱的分析。 该研究处理了具有 C2 v对称性分子基团, 提出一种利用键极化加和模型及其实验修正的方法简化超拉曼微观极化率张量元微分 β'i'j'k'比值的方案。 首先利用键极化加和模型方法, 在对 C2 v分子基团局域模式假设的基础上, 进行单键的局域模式假设, 通过对单键伸缩振动进行简正振动模式的分析, 并与单键的偶极矩表达式进行比较, 得到单键的超拉曼微观极化率张量元微分 β'i'j'k'。 在单键 β'i'j'k'结果的基础上利用键极化加和模型, 对 C2 v分子基团进行对称性分析, 将基团中两个单键伸缩振动的 β'i'j'k'耦合, 得到 C2 v分子基团对称振动和反对称振动两种简正振动模式的超拉曼微观极化率张量元微分 β'i'j'k'的表达式。 在此基础上对 C2 v分子基团对称振动和反对称振动的 β'i'j'k'表达式进行理论推导, 结合对应振动模式下不同偏振组合的超拉曼退偏率 ρHR-SS和 ρHR-AS及HV偏振组合时实验测量的超拉曼光谱强度比, 得到超拉曼实验修正的 C2 v分子基团 β'i'j'k'之间的比值。 再结合文献中拉曼实验修正的拉曼微观极化率张量元微分 α'i'j'之间的比值, 即可得到 C2 v分子基团CARS和CAHRS的微观极化率张量元之间的比例关系, 为定量分析高阶非线性光谱信息提供理论基础。
Coherent Anti-Stokes Raman Spectroscopy (CARS) and Coherent Anti-Stokes Hyper-Raman Spectroscopy (CAHRS), as the high-order nonlinear spectroscopy techniques, have been used in many research studies fields such as dynamic process, gene expression spectrum screening, and high-resolution spectroscopy. However, it is difficult to make a quantitative analysis of the spectral signals that involved many high-order micro-polarizability tensors. It has reported that the CARS and CAHRS microscopic hyperpolarizability tensor elements can be decomposed into the product of the differentiation of Raman microscopic polarizability tensor α'i'j' and Hyper-Raman microscopic polarizability tensor β'i'j'k' which simplifies the difficulty of quantitative analysis of higher-order spectra to the analysis of lower-order spectra. In this paper, we use the bond additivity model (BAM) combined with the experimental correction to address the C2 v symmetry molecular and present the simplified scheme for differentiation of Hyper-Raman microscopic polarizability tensor elements β'i'j'k'. The differentiation of Hyper-Raman microscopic polarizability tensor elements β'i'j'k' of the C2 v molecular group is obtained and the expressions of β'i'j'k' for symmetric and antisymmetric vibrations of the C2 v molecular group also are deduced. According to the experimental measurement of the Hyper-Raman spectral intensity ratio of different polarization combinations of the ρHR-SS and ρHR-AS, along with HV polarization combinations of the corresponding vibration modes, the ratio between β'i'j'k' of C2 v molecular group modified by Hyper-Raman experiment is obtained. Finally, substituting the experimentally corrected differentiation of Raman microscopic polarizability tensor elements α'i'j' in the literature into the ratio above can obtain the proportional relationship between CARS the microscopic polarizability tensor elements and CAHRS of the C2 v molecular group. This method can provide the basis for the quantitative analysis of high-order nonlinear spectral.
近十几年来, 相干反斯托克斯拉曼光谱(Coherent Anti-Stokes Raman spectroscopy, CARS)及其表面增强方法[1, 2, 3, 4]以及相干反斯托克斯超拉曼光谱(Coherent Anti-Stokes hyper-Raman spectroscopy, CAHRS)[5, 6, 7]等高阶非线性光谱技术在动力学过程、 基因表达谱筛选、 高分辨率光谱技术等方面得到重要的应用。 然而上述研究仅仅对光谱信号进行定性分析, 缺乏对光谱信号进行定量分析。 这在某种程度上限制了其应用。 作为高阶非线性相干光学过程, CARS和CAHRS的光谱信号由分子基团的三、 四阶微观极化率张量元贡献[8], 对于n阶非线性相干光谱, 分子微观极化率张量元的数量为3n+1个, 例如CARS涉及81个分子微观极化率张量元。 这些数量众多的分子微观极化率张量元将导致对高阶非线性相干光学过程进行定量分析存在着较大的困难[8]。 加之作为高阶非线性相干光学过程, CARS和CAHRS的实验条件较为苛刻, 为了便于分析结果, 目前使用较多的实验构型为垂直入射式[2]。 但随着对高阶非线性光学过程理论的深入研究和高阶非线性光学技术的发展, 实验技术和实验方法将会逐渐完善, 非垂直入射的非共线实验构型的情况下, 高阶微观极化率张量元的获得对光谱信号定量分析尤为重要。
基于二阶非线性光学过程的和频振动光谱(sum frequency generation vibrational spectroscopy, SFG-VS)定量分析的方法[9, 10, 11], 在以往研究中提出了一种旨在简化CARS和CAHRS的高阶微观极化率张量元的方法[8, 12], 通过对分子基团对称性分析使独立非零的高阶微观极化率张量元数量减少, 由此将CARS微观极化率张量元β i'j'k'l'表示成两个拉曼微观极化率张量元微分α 'i'j'和α 'k'l'的乘积; 将CAHRS微观极化率张量元β i'j'k'l'm'表示成超拉曼微观极化率张量元微分β 'i'j'k'和拉曼微观极化率张量元微分α 'l'm'的乘积, 再利用α 'i'j'之间的比值及β 'i'j'k'之间的比值简化β i'j'k'l'和β i'j'k'l'm'之间的比值, 从而将求解高阶微观极化率张量元比值的问题简化为求解两种低阶微观极化率张量元微分比值的问题。 最后归结为如何求拉曼和超拉曼微观极化率张量元微分比值。
键极化加和模型方法(bond additivity model, BAM)在文献中[9, 13]通常用来求解拉曼微观极化率张量元微分比值。 Wu[14]曾经在讨论空气/甲醇界面的SFG-VS信号时发现, 将传统BAM得出的拉曼微观极化率张量元微分α 'i'j'之间比值带入SFG-VS定量分析方法得到的结果与实验测量不符。 为此其用实验修正的方法对BAM的结果进行了修正, 使其结果与实验测得的界面和频振动光谱能很好地吻合。 本工作利用键极化加和模型及其实验修正的方法, 提出了计算C2v对称性分子基团的超拉曼微观极化率张量元微分β 'i'j'k'的简化方案。
BAM方法是在分子基团局域模式假设的基础上, 进行单个键的局域模式假设。 该模型通过对称性分析, 把基团中单个键的伸缩振动所对应的超拉曼微观极化率张量元微分耦合起来, 得到对应分子基团各简正振动模式的超拉曼微观极化率张量元微分β 'i'j'k'的表达式。 在BAM方法中单键的讨论是整个分子基团讨论的基础, 因此首先讨论图1中具有C∞ v对称性的单键AB的β 'i'j'k'。
单键AB的折合质量GA1和对称振动模式的简正坐标QA1定义如式(1)和式(2)[14]
式(1)和式(2)中MA和MB为原子A和B的质量, Δ r为AB沿着键坐标系(ξ , η , ζ )中ζ 轴的振动矢量。 从表达式(1)和式(2)可得式(3)
对于AB, 分子坐标系(a, b, c)与键坐标系(ξ , η , ζ )完全重合, 因此振动矢量Δ r在分子坐标系中的投影分量分别在式(4)中表达
AB的超拉曼微观极化率张量元微分之间的关系如式(5)
其中
其中Eξ , Eη 和Eζ 分别为外界光场在键坐标系中的投影。 将式(6)—式(8)变换到分子坐标系, 得
将|Δ r|的式(3)代入式(9)—式(11)中, 并与分子坐标系中偶极矩表达式
将β i'j'k'对QA1求导, 可得到C∞ v对称性单键AB的β 'i'j'k'表达式
考虑到式(15)—式(17), 文献[12]中提及的C∞ v对称性分子基团β 'i'j'k'比值
把分子基团AB2中两个单键的β 'i'j'k'耦合起来, 即可得到AB2各简正振动模式的β 'i'j'k'表达式。
对于图2所示的具有C2v对称性的分子基团AB2, 简正模式的折合质量GA1和GB1如式(20)和式(21)定义[14]
其中τ 为AB2的键角。 简正坐标QA1和QB1的定义如式(22)和式(23)[14]
其中|Δ ri|为单键ABi(i=1, 2)沿着各自键坐标系(ξ i, η i, ζ i)中ζ i轴的振动矢量Δ ri的模。 从式(22)和式(23)可得式(24)和式(25)
振动矢量Δ r1和Δ r2在分子坐标系中的投影分量分别为式(26)和式(27)
按照分子基团偶极矩的定义, 单键ABi(i取1, 2)在键坐标系中偶极矩的三个分量表达式见式(28)—式(30)
键AB1和AB2除了取向以外完全相同, 所以两个单键皆具有相同的β 0, r1和r2; Eξ i, Eη i和Eζ i为外界光场在各自键坐标系中的投影, 从其表达式为(见图2)式(31)和式(32)
将两个单键的偶极矩投影到分子坐标系中并加和, 可得式(33)和式(34)
将式(33)—式(35)与分子坐标系中的偶极矩表达式
由|Δ r1|和|Δ r2|的式(24)和式(25)可知
将式(46)和式(47)代入β i'j'k'的式(36)—式(45), 并将所得的β i'j'k'对QA1求导, 得到C2v-A1振动模式β 'i'j'k'表达式如式(48)—式(52)
由式(48)—式(52)即可求得C2v-A1振动模式中β 'i'j'k'的比例关系
从式(53)—式(56)可知, 各比值满足
将β i'j'k'对QB1求导, 即可得C2v-B1振动模式β 'i'j'k'表达式及其比例关系。
通过对AB2分子基团进行入射光和出射光不同偏振组合的超拉曼光谱测量, 可以实现对AB2的β 'i'j'k'进行实验修正。 图3表示测量超拉曼退偏率的实验光路图, (x, y, z)为实验室坐标系, 入射光的偏振为平行(H)或垂直(V)于水平面(xy面)的方向, 出射光在与入射光成90° 的方向进行检测, 检测的偏振方向为垂直(V)方向。
通过选择入射光的偏振可以得到两个偏振组合的超拉曼光谱信号
式(57)中< β
式(58)中RIi为从分子坐标系到实验室坐标系中欧拉变换的矩阵元[8], β 'i'j'k'为超拉曼微观极化率张量元微分。
将τ , r1和r2表示的分子基团某振动模式的β 'i'j'k'代入式(57), 即可得到该分子基团相应振动模式的超拉曼退偏率ρ HR与r1和r2的函数关系。 同时测量对称振动模式和反对称振动模式的超拉曼退偏率ρ HR-SS和ρ HR-AS, 就可利用ρ HR-SS和ρ HR-AS求出r1和r2的值。 例如对于具有C2v对称类型的乙烯中CH2基团, 有键角τ =120° , MC=12, MH=1, 则ρ HR-SS和ρ HR-AS的表达式为
将实验测量的ρ HR-SS和ρ HR-AS代入式(59)和式(60)中即可得CH单键的r1和r2值。
与处理拉曼微观极化率张量元微分相同[14], 在处理不同对称振动模式的超拉曼微观极化率张量元微分比值时需要结合同一偏振组合的超拉曼光谱的实验数据对BAM方法进行进一步修正。 例如对于乙烯的CH2分子基团, HV偏振组合时对称振动模式和反对称振动模式的超拉曼光谱强度比值见式(61)
利用BAM得到的β 'i'j'k'的表达式可以表示任意两个β 'i'j'k'的比值, 例如β
结合式(61)和式(62)可得式(63)
通过结合超拉曼光谱实验测量求出r1和r2, 及对称振动模式和反对称振动模式的超拉曼光谱强度比, 就可得到实验修正的超拉曼微观极化率张量元微分β 'i'j'k'之间的比值, 再结合文献中实验修正的拉曼微观极化率张量元微分α 'i'j'之间的比值, 即可得到CARS和CAHRS的微观极化率张量元之间的比例关系。
(1)本研究以C2v对称性分子基团为例, 利用键极化加和模型方法, 通过对单键伸缩振动进行简正振动模式的分析, 得到单键的超拉曼微观极化率张量元微分β 'i'j'k'。
(2)在单键β 'i'j'k'分析的基础上, 通过分析C2v分子基团的对称性, 得到C2v分子基团对称振动和反对称振动两种简正振动模式的超拉曼微观极化率张量元微分β 'i'j'k'的表达式。
(3)推导了键极化加和模型得到的C2v分子基团β 'i'j'k'表达式, 结合测量得到的不同偏振组合下超拉曼光谱信号强度, 得到修正的C2v分子基团β 'i'j'k'之间的比值。
(4)由文献中实验修正的拉曼微观极化率张量元微分α 'i'j'之间的比值, 即可得到CARS和CAHRS的微观极化率张量元之间的比例关系, 为定量分析高阶非线性光谱信号提供理论基础。
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|