基于键极化加和模型及其实验修正方法的超拉曼微观极化率张量元微分的简化方案- C2 v对称性
汪源1, 张贞2,*, 郭源2,3
1.三亚学院理工学院, 海南 三亚 572022
2.北京分子科学国家研究中心, 中国科学院化学研究所, 北京 100010
3.中国科学院大学化学科学学院, 北京 100049
*通讯作者 e-mail: zhangz@iccas.ac.cn

作者简介: 汪 源, 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的微观极化率张量元之间的比例关系, 为定量分析高阶非线性光谱信息提供理论基础。

关键词: 键极化加和模型; 超拉曼光谱; 微观极化率张量元微分; C2 v对称性
中图分类号:O657.37 文献标志码:A
A Simplified Method of Microscopic Polarizability Tensor Differential of Hyper-Raman Spectroscopy Based on Experimental Correction Bond Additivity Model- C2 v Symmetry
WANG Yuan1, ZHANG Zhen2,*, GUO Yuan2,3
1. Institute of Technology, University of Sanya, Sanya 572022, China
2. Beijing National Laboratory for Molecular Sciences (BNLMS), Institute of Chemistry, Chinese Academy of Sciences, Beijing 100010, China
3. School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
*Corresponding author
Abstract

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.

Keyword: Bond additivity model; Hyper-Raman spectroscopy; Microscopic polarizability tensor differential; C2 v symmetry
引言

近十几年来, 相干反斯托克斯拉曼光谱(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'的简化方案。

1 BAM求解单键的β 'i'j'k'

BAM方法是在分子基团局域模式假设的基础上, 进行单个键的局域模式假设。 该模型通过对称性分析, 把基团中单个键的伸缩振动所对应的超拉曼微观极化率张量元微分耦合起来, 得到对应分子基团各简正振动模式的超拉曼微观极化率张量元微分β 'i'j'k'的表达式。 在BAM方法中单键的讨论是整个分子基团讨论的基础, 因此首先讨论图1中具有Cv对称性的单键ABβ 'i'j'k'

图1 C∞ v对称性分子基团的分子坐标系Fig.1 Molecular coordinate system of C∞ v symmetric molecular group

单键AB的折合质量GA1和对称振动模式的简正坐标QA1定义如式(1)和式(2)[14]

GA1=1/MA+1/MB(1)

QA1=|Δr|/GA1(2)

式(1)和式(2)中MAMB为原子A和B的质量, Δ r为AB沿着键坐标系(ξ , η , ζ )中ζ 轴的振动矢量。 从表达式(1)和式(2)可得式(3)

|Δr|=QA1GA1(3)

对于AB, 分子坐标系(a, b, c)与键坐标系(ξ , η , ζ )完全重合, 因此振动矢量Δ r在分子坐标系中的投影分量分别在式(4)中表达

(Δr)a=0, (Δr)b=0, (Δr)c=|Δr|(4)

AB的超拉曼微观极化率张量元微分之间的关系如式(5)

β'ζξξ=β'ζηη=r1β'ζζζ, β'ξξζ=β'ξζξ=β'ηηζ=β'ηζη=r2β'ζζζ(5)

其中 β'ξηζ=βξηζ/Δζ, Δζ为AB沿着ζ 轴的振动矢量Δ r的模|Δ r|; r1r2为AB的超拉曼微观极化率张量元微分的比值。 为简单计, 令β 0=∂ β ζ ζ ζ /∂ Δ ζ 。 按照分子基团偶极矩的定义, 可以得到键坐标系中单键AB偶极矩的三个分量表达式

μξ=β0(r2EζEξ+r2EξEζ)Δζ=2β0r2EξEζΔζ(6)

μη=β0(r2EζEη+r2EηEζ)Δζ=2β0r2EηEζΔζ(7)

μζ=β0(r1EξEξ+r1EηEη+EζEζ)Δζ(8)

其中Eξ , Eη Eζ 分别为外界光场在键坐标系中的投影。 将式(6)—式(8)变换到分子坐标系, 得

μa=β0(r2EcEa+r2EaEc)|Δr|=2β0r2EaEc|Δr|(9)

μb=β0(r2EcEb+r2EbEc)|Δr|=2β0r2EbEc|Δr|(10)

μc=β0(r1EaEa+r1EbEb+EcEc)|Δr|(11)

将|Δ r|的式(3)代入式(9)—式(11)中, 并与分子坐标系中偶极矩表达式 μi'=j', k'=a, b, cβi'j'k'Ej'Ek'比较, 可得超拉曼微观极化率张量元β i'j'k'表达式如式(12)—式(14)

βaac=βaca=βbbc=βbcb=r2β0QA1GA1(12)

βcaa=βcbb=r1β0QA1GA1(13)

βccc=β0QA1GA1(14)

β i'j'k'QA1求导, 可得到Cv对称性单键AB的β 'i'j'k'表达式

β'aacSS=β'acaSS=β'bbcSS=β'bcbSS=r2β0GA1(15)

β'caaSS=β'cbbSS=r1β0GA1(16)

β'cccSS=β0GA1(17)

考虑到式(15)—式(17), 文献[12]中提及的Cv对称性分子基团β 'i'j'k'比值 rSS, 1HRrSS, 2HR表达式可化简如式(18)和式(19)

rSS, 1HR=β'caa/β'ccc=r1(18)

rSS, 2HR=β'aac/β'ccc=r2(19)

2 BAM求解分子基团AB2β 'i'j'k'

把分子基团AB2中两个单键的β 'i'j'k'耦合起来, 即可得到AB2各简正振动模式的β 'i'j'k'表达式。

对于图2所示的具有C2v对称性的分子基团AB2, 简正模式的折合质量GA1GB1如式(20)和式(21)定义[14]

GA1=(1+cosτ)/MA+1/MB(20)

GB1=(1-cosτ)/MA+1/MB(21)

图2 C2v对称性分子基团的分子坐标系Fig.2 Molecular coordinate system of C2v symmetric molecular group

其中τ 为AB2的键角。 简正坐标QA1QB1的定义如式(22)和式(23)[14]

QA1=(|Δr1|+|Δr2|)/2GA1(22)

QB1=(|Δr1|-|Δr2|)/2GB1(23)

其中|Δ ri|为单键ABi(i=1, 2)沿着各自键坐标系(ξ i, η i, ζ i)中ζ i轴的振动矢量Δ ri的模。 从式(22)和式(23)可得式(24)和式(25)

|Δr1|=12(QA12GA1+QB12GB1)(24)

|Δr2|=12(QA12GA1-QB12GB1)(25)

振动矢量Δ r1和Δ r2在分子坐标系中的投影分量分别为式(26)和式(27)

(Δr1)a=|Δr1|sinτ2, (Δr1)b=0, (Δr1)c=|Δr1|cosτ2(26)

(Δr2)a=-|Δr2|sinτ2, (Δr2)b=0, (Δr2)c=|Δr2|cosτ2(27)

按照分子基团偶极矩的定义, 单键ABi(i取1, 2)在键坐标系中偶极矩的三个分量表达式见式(28)—式(30)

μξi=β0(r2EζiEξi+r2EξiEζi)|Δri|=2β0r2EξiEζi|Δri|(28)

μηi=β0(r2EζiEηi+r2EηiEζi)|Δri|=2β0r2EηiEζi|Δri|(29)

μζi=β0(r1EξiEξi+r1EηiEηi+EζiEζi)|Δri|(30)

键AB1和AB2除了取向以外完全相同, 所以两个单键皆具有相同的β 0, r1r2; Eξ i, Eη iEζ i为外界光场在各自键坐标系中的投影, 从其表达式为(见图2)式(31)和式(32)

Eξ1=-Eacosτ2+Ecsinτ2, Eη1=-Eb, Eζ1=Easinτ2+Eccosτ2(31)

Eξ2=Eacosτ2+Ecsinτ2, Eη2=Eb, Eζ2=-Easinτ2+Eccosτ2(32)

将两个单键的偶极矩投影到分子坐标系中并加和, 可得式(33)和式(34)

μa=-μξ1cosτ2+μζ1sinτ2+μξ2cosτ2-μζ2sinτ2(33)

μb=-μη1+μη22(34)

μc=μξ1sinτ2+μζ1cosτ2+μξ2sinτ2+μζ2cosτ2(35)

将式(33)—式(35)与分子坐标系中的偶极矩表达式 μi'=j', k'=a, b, cβi'j'k'Ej'Ek'进行比较, 可得式(36)—式(45)

βcbb=β0r1(|Δr1|+|Δr2|)cosτ2(36)

βabb=β0r1(|Δr1|-|Δr2|)sinτ2(37)

βbbc=βbcb=β0r2(|Δr1|+|Δr2|)cosτ2(38)

βbba=βbab=β0r2(|Δr1|-|Δr2|)sinτ2(39)

βcaa=β0(|Δr1|+|Δr2|)cosτ2·sin2τ2+r1cos2τ2-2r2sin2τ2(40)

βccc=β0(|Δr1|+|Δr2|)cosτ2·cos2τ2+r1sin2τ2+2r2sin2τ2(41)

βacc=β0(|Δr1|-|Δr2|)sinτ2·cos2τ2+r1sin2τ2-2r2cos2τ2(42)

βaaa=β0(|Δr1|-|Δr2|)sinτ2·sin2τ2+r1cos2τ2+2r2cos2τ2(43)

βaac=βaca=β0(|Δr1|+|Δr2|)cosτ2·sin2τ2-r1sin2τ2+r2cosτ(44)

βcca=βcac=β0(|Δr1|-|Δr2|)sinτ2·cos2τ2-r1cos2τ2-r2cosτ(45)

由|Δ r1|和|Δ r2|的式(24)和式(25)可知

|Δr1|+|Δr2|=Qa2Ga(46)

|Δr1|-|Δr2|=Qb2Gb(47)

将式(46)和式(47)代入β i'j'k'的式(36)—式(45), 并将所得的β i'j'k'QA1求导, 得到C2v-A1振动模式β 'i'j'k'表达式如式(48)—式(52)

β'cbbSS=β02GA1r1cosτ2(48)

β'bbcSS=β'bcbSS=β02GA1r2cosτ2(49)

β'caaSS=β02GA1cosτ2sin2τ2+r1cos2τ2-2r2sin2τ2(50)

β'aacSS=β'acaSS=β02GA1cosτ2sin2τ2-r1sin2τ2+r2cosτ(51)

β'cccSS=β02GA1cosτ2cos2τ2+r1sin2τ2+2r2sin2τ2(52)

由式(48)—式(52)即可求得C2v-A1振动模式中β 'i'j'k'的比例关系

RSS⁃1HR=β'acaSSβ'cccSS=1-r1+r2cot2τ2-1cot2τ2+r1+2r2(53)

RSS⁃2HR=β'caaSSβ'cccSS=1+r1cot2τ2-2r2cot2τ2+r1+2r2(54)

RSS⁃3HR=β'bcbSSβ'cccSS=r2csc2τ2cot2τ2+r1+2r2(55)

RSS⁃4HR=β'cbbSSβ'cccSS=r1csc2τ2cot2τ2+r1+2r2(56)

从式(53)—式(56)可知, 各比值满足 RSS⁃1HR+ RSS⁃4HR= RSS⁃2HR+ RSS⁃3HR的关系, 且各比值只与AB2的键角τ 及单键AB的超拉曼微观极化率张量元微分比值r1r2有关, 而r1r2可用实验测量超拉曼退偏率的方法得到[14]

β i'j'k'QB1求导, 即可得C2v-B1振动模式β 'i'j'k'表达式及其比例关系。

3 AB2分子基团β 'i'j'k'的实验修正

通过对AB2分子基团进行入射光和出射光不同偏振组合的超拉曼光谱测量, 可以实现对AB2β 'i'j'k'进行实验修正。 图3表示测量超拉曼退偏率的实验光路图, (x, y, z)为实验室坐标系, 入射光的偏振为平行(H)或垂直(V)于水平面(xy面)的方向, 出射光在与入射光成90° 的方向进行检测, 检测的偏振方向为垂直(V)方向。

图3 超拉曼光谱实验光路图Fig.3 Experimental light road of Hyper-Raman spectroscopy

通过选择入射光的偏振可以得到两个偏振组合的超拉曼光谱信号 IHRHVIHRVV。 分子基团AB2的超拉曼退偏率 ρHR[14]定义为

ρHR=IHRHV/IHRVV=< β'zyy2> /< β'zzz2> (57)

式(57)中< β 'IJK2> (I, J, K=x, y, z)的表达式[14]

< β'IJK2> =i, j, k, i', j', k'β'ijkβ'i'j'k'·02π02π0πRIiRJjRKkRIi'RJj'RKk'sinθdθdΨdϕ02π02π0πsinθdθdΨdϕ(58)

式(58)中RIi为从分子坐标系到实验室坐标系中欧拉变换的矩阵元[8], β 'i'j'k'为超拉曼微观极化率张量元微分。

τ , r1r2表示的分子基团某振动模式的β 'i'j'k'代入式(57), 即可得到该分子基团相应振动模式的超拉曼退偏率ρ HRr1r2的函数关系。 同时测量对称振动模式和反对称振动模式的超拉曼退偏率ρ HR-SSρ HR-AS, 就可利用ρ HR-SSρ HR-AS求出r1r2的值。 例如对于具有C2v对称类型的乙烯中CH2基团, 有键角τ =120° , MC=12, MH=1, 则ρ HR-SSρ HR-AS的表达式为

ρHR⁃SS=ISS, HRHVISS, HRVV=24+8r1-96r2-16r1r2+108r12+96r2239+6r1+12r2+100r1r2+25r12+100r22(59)

ρHR⁃AS=IAS, HRHVIAS, HRVV=8+40r1-32r2-80r1r2+92r12+32r2255+62r1+116r2+100r1r2+23r12+108r22(60)

将实验测量的ρ HR-SSρ HR-AS代入式(59)和式(60)中即可得CH单键的r1r2值。

与处理拉曼微观极化率张量元微分相同[14], 在处理不同对称振动模式的超拉曼微观极化率张量元微分比值时需要结合同一偏振组合的超拉曼光谱的实验数据对BAM方法进行进一步修正。 例如对于乙烯的CH2分子基团, HV偏振组合时对称振动模式和反对称振动模式的超拉曼光谱强度比值见式(61)

ISS, HRHVIAS, HRHV=24+8r1-96r2-16r1r2+108r12+96r228+40r1-32r2-80r1r2+92r12+32r22GA1GE(61)

利用BAM得到的β 'i'j'k'的表达式可以表示任意两个β 'i'j'k'的比值, 例如β 'cccSSβ 'aaaAS的比值为

β'cccSSβ'aaaAS=331+3r1+6r23+r1+2r2GA1GE(62)

结合式(61)和式(62)可得式(63)

β'cccSSβ'aaaAS=ISS, HRHVIAS, HRHV138+40r1-32r2-80r1r2+92r12+32r2224+8r1-96r2-16r1r2+108r12+96r22·1+3r1+6r23+r1+2r2(63)

通过结合超拉曼光谱实验测量求出r1r2, 及对称振动模式和反对称振动模式的超拉曼光谱强度比, 就可得到实验修正的超拉曼微观极化率张量元微分β 'i'j'k'之间的比值, 再结合文献中实验修正的拉曼微观极化率张量元微分α 'i'j'之间的比值, 即可得到CARS和CAHRS的微观极化率张量元之间的比例关系。

4 结论

(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] Wallance M P, Leonardo N F, Vita S, et al. Appl. Spec. , 2020, 74(7): 751. [本文引用:1]
[2] Sinjab F, Hashimoto K, Zhao X Q, et al. Opt. Lett. , 2020, 45(6): 1515. [本文引用:2]
[3] Liu G Y, Narducci F A, Malinovskaya S A. J. Mod. Opt. , 2020, 67(1): 21. [本文引用:1]
[4] PENG Ya-jing, SUN Shuang, SONG Yun-fei, et al(彭亚晶, 孙爽, 宋云飞, ). Acta Phys. Sin. (物理学报), 2018, 67(2): 152. [本文引用:1]
[5] Fumero G, Batignani G, Dorfman K E, et al. Chem. Phys. Chem. , 2015, 16(16): 3438. [本文引用:1]
[6] Fumero G, Batignani G, Dorfman K E, et al. Chem. Phys. Chem. , 2015, 16(16): 3533. [本文引用:1]
[7] Zhang Z D, Bennett K, Chernyak V, et al. J. Phys. Chem. Lett. , 2017, 8(14): 3387. [本文引用:]
[8] Wang Y, Cui Z F, Wang H F. Chin. J. Chem. Phys. , 2007, (20): 449. [本文引用:]
[9] Wang H F, Gan W, Lu R, et al. Int. Rev. Phys. Chem. , 2005, (24): 191. [本文引用:]
[10] Wang H F. Prog. Surf. Sci. , 2016, 91: 155. [本文引用:]
[11] Wang H F, Velarde L, Gan W, et al. Annu. Rev. Phys. Chem. , 2015, 66: 189. [本文引用:]
[12] WANG Yuan, ZHANG Zhen, GUO Yuan(汪源, 张贞, 郭源). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2019, 39(2): 465. [本文引用:]
[13] Zhang D, Gutow J, Eisenthal K B. J. Phys. Chem. , 1994, 98: 13729. [本文引用:]
[14] Wu H, Zhang W K, Gan W, et al. J. Chem. Phys. , 2006, 125: 133203. [本文引用:]