双光谱二维异步相关光谱表征分子间相互作用的可行性研究
李晓佩1,*, 张勇杰2, 薛丽贞1
1. 大连工业大学实验仪器中心, 辽宁 大连 116034
2. 大连工业大学纺织与材料工程学院, 辽宁 大连 116034
*通讯联系人

作者简介: 李晓佩, 1984年生, 大连工业大学实验仪器中心讲师 e-mail: li_xp@dlpu.edu.cn

摘要

双光谱二维异步相关光谱(2T2D-异步相关光谱)是一种利用两幅一维光谱创建异步相关光谱的新方法。 相对于最少需要三幅一维光谱的传统异步相关光谱, 2T2D-异步相关光谱可简化实验过程, 有利于样品昂贵体系的表征。 利用数学分析、 模拟及实际体系实验对2T2D-异步相关光谱应用于表征分子间相互作用的可行性进行了研究。 首先建立一个包含P和Q两种溶质的模拟体系, 设定P有光谱峰, Q没有光谱峰。 数学分析表明: (1) 不正确设置P和Q初始浓度, 可导致2T2D-异步相关光谱的强度恒为零; (2) 2T2D-异步相关光谱不能反映与分子间相互作用相关的峰强变化; 因此, 利用2T2D-异步相关光谱表征分子间相互作用可能得到错误的结论。 为将2T2D-异步相关光谱发展成为表征分子间相互作用的可靠方法, 首先对2T2D-异步相关光谱中P和Q初始浓度的设置方法进行了研究, 得出当P和Q初始浓度满足文中式(6)时, 可避免不正确设置P和Q初始浓度导致2T2D-异步相关光谱强度恒为零的情况; 在此基础上, 为解决2T2D-异步相关光谱不能反映与分子间相互作用相关的峰强变化问题, 通过向体系中加入具有独立光谱峰和适合浓度的虚拟物质S, 发展出带有辅助交叉峰的2T2D-异步相关光谱(ASAP-2T2D-异步相关光谱)。 模拟体系实验表明ASAP-2T2D-异步相关光谱可正确反映与分子间相互作用相关的峰宽、 峰位及峰强变化, 是表征分子间相互作用的可靠方法。 最后, 将ASAP-2T2D-异步相关光谱应用于表征苯并-15-冠醚-5(BC)与Li+间的相互作用, 实验结果表明ASAP-2T2D-异步相关光谱可同时反映BC特征峰的峰位及峰强变化, 进一步证实ASAP-2T2D-异步相关光谱可正确表征分子间相互作用。

关键词: 分子间相互作用; 2T2D-异步相关光谱; 辅助交叉峰
中图分类号:O657.3 文献标志码:A
Application of Two-Trace Two-Dimensional Asynchronous Correlation Spectroscopy in Characterizing Intermolecular Interactions: A Feasibility Study
LI Xiao-pei1,*, ZHANG Yong-jie2, XUE Li-zhen1
1. Instrumental Analysis Center, Dalian Polytechnic University, Dalian 116034, China
2. School of Textile and Material Engineering, Dalian Polytechnic University, Dalian 116034, China
*Corresponding author
Abstract

Two-trace two-dimensional (2T2D) asynchronous correlation spectroscopy is a new method of generating two-dimensional asynchronous correlation spectroscopy based on a pair of one-dimensional (1D) spectra. Compared with the conventional two-dimensional asynchronous correlation spectroscopy where at least three 1D spectra are needed, 2T2D asynchronous correlation spectroscopy makes the experiment easier, and is a better method for the expensive samples. In the present paper, the feasibility of using 2T2D asynchronous correlation spectroscopy to characterize intermolecular interactions was explored. Firstly, a model system containing two solutes P and Q were set up. P possessed a characteristic peak, and Q had no characteristic peaks. The mathematical analysis demonstrated that the intensity of 2T2D asynchronous correlation spectra was always zero when the initial concentrations of P and Q were set incorrectly. Furthermore, the variations of the absorptivity induced by intermolecular interactions could not be discerned by 2T2D-asynchronous correlation spectra. Therefore, incorrect results may be obtained when 2T2D asynchronous correlation spectroscopy is adopted to characterize intermolecular interactions. In order to develop 2T2D asynchronous correlation spectroscopy into a reliable method for characterizing intermolecular interactions, the setting method of the initial concentrations of P and Q in the 2T2D asynchronous correlation spectroscopy was first studied. The case that the intensity of 2T2D asynchronous correlation spectra was always zero induced by incorrectly setting the initial concentrations of P and Q could be avoided when the initial concentrations of P and Q satisfied the requirement of Eq.(6) in the text. On this basis, the 2T2D-asynchronous correlation spectroscopy with auxiliary cross peaks (ASAP-2T2D-asynchronous correlation spectroscopy) was developed by introducing a virtual substance S with an isolated peak and proper concentration into the research system to solve the problem that 2T2D-asynchronous correlation spectra could not reflect the variations of the absorptivity induced by intermolecular interactions. The results of computer simulation experiments demonstrated that the ASAP-2T2D-asynchronous correlation spectra could correctly reflect the variation of the peak position, bandwidth and absorptivity indicating that the ASAP-2T2D-asynchronous correlation spectroscopy was a reliable method for characterizing intermolecular interactions. Finally, the ASAP-2T2D-asynchronous correlation spectroscopy was adopted to characterize the intermolecular interaction between Li+ and benzo-15-crown-5 (BC). The results demonstrated that the variation of the peak position and absorptivity of the characteristic peak of BC could be reflected by the ASAP-2T2D-asynchronous correlation spectrum. These results further confirmed that the ASAP-2T2D-asynchronous correlation spectroscopy could correctly characterize intermolecular interactions.

Keyword: Intermolecular interactions; 2T2D-asynchronous correlation spectroscopy; Auxiliary cross peaks
引 言

分子间相互作用, 包括范德华力[1]、 氢键[2]、 卤键[3]、 静电作用[4]、 疏水作用[5]、 π — π 堆积[6]等, 是众多化学、 生物、 物理体系的重要研究内容。 光谱法是一种重要的分子间相互作用表征方法[7, 8, 9, 10]。 分子间相互作用的光谱分析基于光谱的变化, 然而众多实际体系中的分子间相互作用引起的一维光谱变化不明显, 甚至淹没在相互作用分子原有的谱峰中, 这使得光谱法在分子间相互作用研究中的应用受到了很大的限制。

二维相关光谱是对外部扰动下获得的一系列一维光谱进行相关分析得到的谱图。 由于扩展了一个维度, 二维相关光谱相对于一维光谱具有简化重叠峰和增强光谱分辨率的突出优势, 为提取淹没在一维光谱中的分子间相互作用信息提供了契机[11, 12]。 二维相关光谱分为同步相关光谱和异步相关光谱; 异步相关光谱没有自相关峰, 具有更为丰富的交叉峰峰型, 这两个特点决定异步相关光谱更适合应用于分子间相互作用的表征。 传统异步相关光谱的创建最少需要三幅一维光谱[13]。 最近, 二维相关光谱的创始人Noda教授提出一种利用两幅一维光谱创建异步相关光谱的新方法[式(1)], 即2T2D-异步相关光谱[13]。 2T2D-异步相关光谱可反映出两幅一维光谱的细微差别, 这种细微差别无法通过肉眼观察和传统的光谱比较法(比如加权光谱差减法和皮尔逊相关系数)鉴别开。 因此, 2T2D-异步相关光谱在质量控制和精密检测等方面有重要的应用前景。 鉴于2T2D-异步相关光谱与传统异步相关光谱有许多共性, 如能将2T2D-异步相关光谱应用于分子间相互作用的检测, 则有望进一步简化实验过程, 更有利于样品昂贵体系的表征。 本工作就2T2D-异步相关光谱应用于表征分子间相互作用的可行性进行了研究。

Ψ2T2D(x, y)=12[A1(x)A2(y)-A1(y)A2(x)](1)

其中A1(x), A1(y), A2(x)和A2(y)为两幅一维光谱在坐标xy处的强度。

1 实验部分
1.1 模拟体系

模拟体系包含P和Q两种溶质, 在分子间相互作用下, 一部分P转化为U, 一部分Q转化为V[式(2)]。

$P+O \mathop{}_{\rightleftarrows}^{K} U+V$ (2)

设定Q没有特征光谱峰, P有特征光谱峰, 且与溶剂光谱峰不发生重叠。 符合朗伯-比尔定律的模拟体系一维光谱表达式如式(3)所示。

A(ν)=fP(ν)cP(init)+[fU(ν)-fP(ν)]cU(eq)(3)

设定光程为1; ν 为波长; cP(init)cU(eq)分别为模拟溶液中P的初始浓度及U的平衡浓度; fP(ν )和fU(ν )为光谱项, 用高斯函数描述[式(4)]。

fjν=εjgjν=εje-(ln2)(ν-νj)2wj2(4)

其中j表示P, U或S。 ε j, ν jwj分别为摩尔吸光系数、 峰位和半高半宽。

1.2 实际体系

苯并-15-冠醚-5/氯化锂/甲醇体系, 具体实验过程参考文献[14]。 模拟及实际体系数据处理程序由MATLAB R2014a软件编写。

2 结果与讨论
2.1 2T2D-异步相关光谱表征分子间相互作用的可靠性研究

为正确表征分子间相互作用, 2T2D-异步相关光谱应同时具有两个性质: (1)P和Q间不存在相互作用时, 2T2D-异步相关光谱的强度为零; (2)P和Q间存在相互作用时, 2T2D-异步相关光谱的强度不为零, 且2T2D-异步相关光谱可反映与分子间相互作用相关的P光谱峰的峰宽、 峰位和峰强变化。 基于这两个性质, 利用数学分析和模拟体系实验对2T2D-异步相关光谱表征分子间相互作用的可靠性进行研究。

配置两个P/Q溶液, 基于式(1)和式(3), 可得式(5)。

Ψ2T2D(x, y)=F(x, y)c(5)

c=12(cP1(init)cU2(eq)-cP2(init)cU1(eq))F(x, y)=fP(x)fU(y)-fP(y)fU(x)

当P和Q间不存在相互作用时, cU1(init)cU2(init)为零, 则c=0; 因此2T2D-异步相关光谱的强度为零, 表明2T2D-异步相关光谱满足性质(1)。

当P和Q间存在相互作用时, 根据性质(2), 应同时满足c≠ 0和F(x, y)≠ 0。 首先对c进行分析; c为两项的差值, 不正确选择P和Q的初始浓度, 可导致c=0; 比如, 当 cP1(init)=0时, 对应的 cU1(init)=0, 因此c=0; 基于c的表达式可得: 当两个溶液中P和Q的初始浓度满足式(6)时可使c≠ 0。

cP1(init)0cP2(init)0cQ1(init)0cP1(init)0cP2(init)0cQ2(init)0(6)

F(x, y)进行分析, 将式(4)代入F(x, y)的表达式中可得

F(x, y)=εPεU[gP(x)gU(y)-gP(y)gU(x)](7)

当P和Q间的相互作用只导致峰强变化时, gP(x)=gU(x)和gP(y)=gU(y), 因此F(x, y)=0; 表明即使P和Q间存在相互作用, 2T2D-异步相关光谱的强度依然可能为零。

以下模拟实验, 进一步证实了以上数学分析结果。 P和U的光谱参数如表1所示; P和Q的初始浓度按照式(6)进行设置, 列于表2。 图1(A)给出P和Q间不存在相互作用及P和Q间的相互作用使P的峰宽、 峰位及峰强发生变化时, 表2中第一个溶液对映一维光谱。 由于分子间相互作用的强度弱(K=0.01)及分子间相互作用导致的谱峰变化较小(ν U-ν P=1; WU-WP=1; ε U-ε P=-0.02), 四幅一维光谱没有明显的区别, 基于一维光谱无法得到P和Q间的相互作用信息。 图1(B)给出图1(A)所列四种情况对应的2T2D-异步相关光谱。 当P和Q间不存在相互作用时, 2T2D-异步相关光谱的强度为零[图1(B(a))]。 当分子间相互作用使P的峰宽和峰位发生变化时, 2T2D-异步相关光谱都出现了明显的交叉峰[图1(B(b, c))], 且二者交叉峰峰型有明显差别, 表明2T2D-异步相关光谱可反映与分子间相互作用相关的P光谱峰的峰宽和峰位变化。 然而, 当分子间相互作用仅使P的峰强发生变化时, 2T2D-异步相关光谱的强度为零[图1(B(d))]。

表1 模拟体系中P和U的谱峰参数 Table 1 Peak parameters of the chemical species P and U in the model system
表2 模拟体系中P和Q的初始浓度 Table 2 Initial concentrations of the chemical species P and Q in the model system

图1 (A) P和Q间不存在相互作用(表1ν U=100, WU=20和ε U=1, a), P和Q间的相互作用使P的峰宽(ν U=100, WU=21和ε U=1, b)、 峰位(ν U=101, WU=20和ε U=1, c)及峰强(ν U=100, WU=20和ε U=0.98, d)发生变化时, 表2中第一个溶液对应的一维光谱; (B)图1(A)所列四种情况对应的2T2D-异步相关光谱Fig.1 (A) The one dimensional spectra of the first solution in Table 2 of the model system where no intermolecular interactions exist between P and Q (ν U=100, WU=20和ε U=1 in Table 1, a) and intermolecular interactions between P and Q cause thebandwidth (ν U=100, WU=21和ε U=1, b), peak position (ν U=101, WU=20和ε U=1, c) and absorptivity (ν U=100, WU=20和ε U=0.98, d) variation of the characteristic peak of P; (B) 2T2D-asynchronous correlation spectra in the conditions listed in Figure 1(A)

结果表明: 2T2D-异步相关光谱满足性质1, 但不满足性质2。 因此2T2D-异步相关光谱不能作为表征分子间相互作用的可靠方法。

2.2 带有辅助交叉峰的2T2D-异步相关光谱(ASAP-2T2D-异步相关光谱)

仅P有光谱峰时, 传统异步相关光谱也不能反映与分子间相互作用相关的峰强变化[14]。 通过向体系中引入具有独立光谱峰的虚拟物质S, 传统异步相关光谱解决了该问题。 基于虚拟物质的传统异步相关光谱称为带有辅助交叉峰的异步相关光谱(ASAP-异步相关光谱)[14]。 借助ASAP-异步相关光谱的思想, 将S引入到2T2D-异步相关光谱中; 基于虚拟物质的2T2D-异步相关光谱命名为ASAP-2T2D-异步相关光谱。 以下, 通过数学分析、 模拟及实际体系实验对ASAP-2T2D-异步相关光谱表征分子间相互作用的可靠性进行研究。

由于S的加入, 式(3)转变为式(8)

A(ν)=fP(ν)cP(init)+[fU(ν)-fP(ν)]cU(eq)+fS(ν)cS(8)

其中fS(ν )为光谱项, 用高斯函数描述[式(4)]; cS为S的浓度。

基于式(1)和式(8)可得

Ψ2T2D(x, y)=12i=13WiW1=[fP(x)fS(y)-fP(y)fS(x)]·[cP1(init)cS⁃2-cS1cP2(init)]W2=[fP(x)fU(y)-fP(y)fU(x)]·[cP1(init)cU2(eq)-cP2(init)cU1(eq)]W3=[fU(x)fS(y)+fP(y)fS(x)-fP(x)fS(y)-fU(y)fS(x)][c(eq)U1cS2-cS1cU2(eq)](9)

当P和Q间不存在相互作用时, cU1(init)cU2(init)为零, 则W2=0, W3=0。 这时ASAP-2T2D-异步相关光谱还有一个非零项W1, 因此P和Q间不存在相互作用时, ASAP-2T2D-异步相关光谱的强度并不为零, 表明ASAQ-2T2D-异步相关光谱中存在与分子间相互作用无关的假峰。 以下模拟体系实验表明, ASAP-2T2D-异步相关光谱的假峰对分子间相互作用的正确识别造成严重干扰。 模拟体系中, 第一和第二个溶液中S的浓度分别设置为7和3; S的峰位、 半高半宽和摩尔吸光系数分别设置为250, 20和1。 如图2(A)所示, 当P和Q间不存在相互作用及P和Q间的相互作用使P的峰宽、 峰位及峰强发生变化时, ASAP-2T2D-异步相关光谱都存在相同峰型的交叉峰, 表明ASAP-2T2D-异步相关光谱不能正确表征分子间相互作用。

图2 P和S的浓度设置不满足式(10)(A)和满足式(10)(B)时, (A)所列四种情况对应的ASAP-2T2D-异步相关光谱Fig.2 ASAP-2T2D-asynchronous correlation spectra in the conditions listed in (A) when the concentrations of P and S do not satisfy the requirement of Eq.(10) (A) and do satisfy the requirement of Eq.(10) (B)

为正确表征分子间相互作用, ASAP-2T2D-异步相关光谱中的W1项应该去除。 根据W1的表达式可得: 当P和S的浓度设置满足式10时W1=0。

cP1(init)cS2-cS1cP2(init)=0(10)

如图2(B)所示[根据式(10), 第一和第二个溶液中S的浓度分别设置为2和8], 当P和Q间不存在相互作用时, ASAP-2T2D-异步相关光谱的强度为零[图2B(a)]; 当P和Q间的相互作用使P的峰宽、 峰位及峰强发生变化时, ASAP-2T2D-异步相关光谱中均出现明显的交叉峰[图2B(b— d)]。 以图2B(b)为例, ASAP-2T2D-异步相关光谱可分为三个区: 一区的交叉峰即2T2D-异步相关光谱中的交叉峰, 不能反映与分子间相互作用相关的峰强变化[图2B(d)中一区的交叉峰强度为零]; 二区和三区的交叉峰为源于S的辅助交叉峰, 二者关于对角线呈反对称。 图2B(b— d)中辅助交叉峰的峰型各不相同, 表明按照式(10)设置S浓度时, ASAP-2T2D-异步相关光谱可反映P的峰位、 峰宽及峰强变化。 图2(B)证实ASAP-2T2D-异步相关光谱可同时满足性质1和2, 是表征分子间相互作用的可靠方法。

基于以上结论, ASAP-2T2D-异步相关光谱进一步应用于实际体系中分子间相互作用的表征。 在前期工作中, 配置了四个溶液, 利用ASAP-异步相关光谱对苯并-15-冠醚-5(BC)与Li+之间的相互作用进行了研究[14]。 从四个溶液中选取前两个进行ASAP-2T2D-异步相关光谱分析。 根据式(10), S的浓度分别设置为0.41和0.30 mol· L-1; S的峰位、 半高半宽和摩尔吸光系数分别设置为1450, 10和1。 如图3所示, BC/Li+/CH3OH体系对应的ASAP-2T2D-异步相关光谱的交叉峰峰型与ASAP-异步相关光谱的交叉峰峰型完全相同[14], 即一区中的交叉峰峰型表明Li+与BC之间的相互作用使BC特征峰的峰位发生了变化, 而二区和三区中的交叉峰左右两部分强度的绝对值不同(左右两部分强度的绝对值大小表现为等高线的圈数, 圈数越多, 绝对值越大), 表明Li+与BC之间的相互作用使BC特征峰的峰强也发生了变化。 以上结果进一步证实了ASAP-2T2D-异步相关光谱可正确表征分子间相互作用。

图3 BC和S浓度设置满足式(10)时, BC/Li+/CH3OH体系对应的ASAP-2T2D-异步相关光谱Fig.3 ASAP-2T2D-asynchronous correlation spectrum of BC/Li+/CH3OH system when the concentrations of BC and S satisfy the requirement of Eq.(10)

3 结 论

对2T2D-异步相关光谱应用于表征分子间相互作用的可行性进行了研究。 当P有光谱峰, Q没有光谱峰时, 数学分析表明P和Q的初始浓度需满足文中式(6)的要求, 否则2T2D-异步相关光谱的强度恒为零。 数学分析和模拟体系实验表明2T2D-异步相关光谱只能反映与分子间相互作用相关的峰宽和峰位的变化, 而不能反映峰强的变化。 针对此问题, 通过向体系中加入具有独立光谱峰和适合浓度的虚拟物质S, 发展出ASAP-2T2D-异步相关光谱。 模拟体系实验表明ASAP-2T2D-屿步相关光谱可正确反映与分子间相互作用相关的峰宽、 峰位及峰强变化, 是表征分子间相互作用的可靠方法。 实际体系实验的结果进一步证实了此结论。

参考文献
[1] Jia H Y, Huang Z J, Fei Z F, et al. Journal of Materials Chemistry B, 2017, 5: 8193. [本文引用:1]
[2] Mitsuhashi R, Suzuki T, Hosoya S, et al. Crystal Growth & Design, 2017, 17: 207. [本文引用:1]
[3] Kolář M H, Hobza P. Chemical Reviews, 2016, 116: 5155. [本文引用:1]
[4] Chen P Y, Zhang L, Zhu S G, et al. Journal of Molecular Structure, 2017, 1131: 250. [本文引用:1]
[5] Yim D B, Kim J E, Kim H I, et al. Small, 2018, 14: 1800026. [本文引用:1]
[6] Kang M, Zhang P C, Cui H G, et al. Macromolecules, 2016, 49: 994. [本文引用:1]
[7] WANG Yan, Da, GUO Ming, et al(王燕, 吕达, 郭明, ). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2018, 38(2): 494. [本文引用:1]
[8] Alanazi A M, Abdelhameed A S, Bakheit A H, et al. Journal of Molecular Liquids, 2017, 238: 3. [本文引用:1]
[9] Yang F L, Yang X, Wu R Z, et al. Physical Chemistry Chemical Physics, 2018, 20: 11386. [本文引用:1]
[10] HAN Guo-cheng, SU Xiao-rui, HOU Jia-ting, et al(韩国成, 苏晓瑞, 侯嘉婷, ). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2018, 38(12): 3958. [本文引用:1]
[11] Wu Y Q, Zhang L P, Jung Y M, et al. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2018, 189: 291. [本文引用:1]
[12] He A Q, Zeng X Z, Xu Y Z, et al. Journal of Physical Chemistry A, 2017, 121: 7524. [本文引用:1]
[13] Noda I. Journal of Molecular Structure, 2018, 1160: 471. [本文引用:2]
[14] Li X P, He A Q, H K, et al. RSC Advances, 2015, 5: 87739. [本文引用:5]