拉曼光谱在新型毒品快速检测中的应用
隋成华1,2,*, 屠潜1, 鄢波1,2, 陈乃波1,2, 卢忠1,2, 吕斌1,2
1. 浙江工业大学理学院, 浙江 杭州 310023
2. 浙江工业大学生物与医学物理信息技术协同创新中心, 浙江 杭州 310023
*通讯联系人

作者简介: 隋成华, 1956年生, 浙江工业大学理学院应用物理系教授 e-mail: suich@zjut.edu.cn

摘要

新型毒品的日益泛滥及快速更新对执法部门现场快速检测提出了越来越高的要求。 以三种典型的新型毒品为例, 利用密度泛函理论中的B3LYP杂化泛函, 在6-31G基组下优化分子几何结构以及拉曼振动频率的计算, 并利用拉曼光谱仪进行了实验检测, 用以研究拉曼光谱技术在新型毒品快速检测中的应用价值。 结果表明: 各毒品样品的理论计算拉曼光谱与实验拉曼光谱基本吻合, 理论计算光谱可以为实验光谱特征峰的归属提供参考; 各毒品拉曼特征峰的峰位差异明显, 其中冰毒特征峰为837和1003 cm-1, K粉最具有鉴别价值的特征峰为463, 659和1 046 cm-1, 而麻古最明显的特征峰为556, 1 329和1 699 cm-1, 拉曼光谱可被用于毒品的鉴别和认定; 微量冰毒及K粉残留物的实验拉曼光谱与其常量时基本一致, 拉曼光谱可被用于毒品微量残留物的准确识别; 伪冰毒N-异丙基苄胺的拉曼特征峰853 cm-1与冰毒拉曼特征峰837 cm-1存在明显差异, 因此拉曼光谱可被有效地用于以上毒品的检测; 而聚丙烯材料制成的透明包装对自身存在强烈荧光干扰的麻古拉曼光谱有较大影响。

关键词: 新型毒品; 拉曼光谱; 快速检测; 密度泛函
中图分类号:O657.3 文献标志码:A
Application of Raman Spectroscopy in Rapid Detection of New-Type Drugs
SUI Cheng-hua1,2,*, TU Qian1, YAN Bo1,2, CHEN Nai-bo1,2, LU Zhong1,2, LÜ Bin1,2
1. College of Science, Zhejiang University of Technology, Hangzhou 310023, China
2. Collaborative Innovation Center for Information Technology in Biological and Medical Physics, Hangzhou 310023, China
Abstract

The increasing spread and rapid updating of new-type drugs have put forward higher and higher requirements for the rapid detection by law enforcement departments. In this paper, Raman spectroscopy in the rapid detection of new-type drugs was studied by theoretical calculation and experimental detection of three typical new-type drugs. The molecules of the three drugs were optimized and calculated by B3LYP hybrid functional and 6-31G basis set. The Raman spectrometer was used for the experimental detection. The results indicated that the theoretical Raman spectra tallies with the experimental Raman spectra, and can be used for providing a reference for the attribution of experimental spectral peaks. The Raman peak positions of the three drugs are significantly different. The characteristic Raman peaks of methamphetamine are located at 837 and 1 003 cm-1. Those of ketamine are located at 463, 659 and 1 046 cm-1. And the most obvious characteristic peaks of Ma Gu are located at 556, 1 329 and 1 699 cm-1. Clearly the Raman spectra can be used for the identification of drugs. The Raman spectra of methamphetamine and ketamine residues tally with those obtained from the constant samples, which just shows that drug residues can be accurately identified by Raman spectroscopy. The characteristic Raman peak of pseudomethamphetamine (N-benzylisopropylamine) 853 cm-1 is significantly different from that of methamphetamine 837 cm-1. Therefore the Raman spectra can also be used to identify the authenticity of new-type drugs. Moreover, transparent packaging has a great influence on the Raman spectrum of Ma Gu which has strong fluorescence interference.

Keyword: New-type drugs; Raman spectroscopy; Rapid detection; Density functional theory
引 言

人工合成的新型毒品因其合成方式多样、 制备简单、 品种繁多、 价格低廉、 隐蔽性较强等特点, 近年来呈现出日益泛滥的趋势, 且已蔓延至以青少年群体为主的社会各个阶层, 给社会稳定、 经济发展、 人民健康造成了严重危害, 同时也给缉毒部门就毒品和疑似毒品的现场识别和快速检测带来了极大的挑战。

常用的毒品检测方法主要有化学法、 色谱法和光谱法等[1, 2, 3], 其中化学法简单快速, 但不易检测出微量或化学结构相近的毒品, 且结果受杂质干扰较大; 色谱法需对样品进行预处理, 程序较为复杂, 且需专业人员位于实验室操作, 无法做到快速检测。 相比之下, 拉曼光谱法因其具有实时、 快速、 准确、 非接触、 无损、 样品量小、 样品无需预处理等一系列优点, 在毒品检测研究领域备受青睐。 目前, 研究人员已利用拉曼光谱技术对苯丙胺类毒品, 例如甲基苯丙胺(俗称冰毒), 进行过检测并对这些毒品的拉曼光谱特征峰进行了区分和总结[4]; 也曾利用拉曼光谱对致幻剂类毒品, 例如氯胺酮(俗称K粉), 做过检测分析[5, 6]。 但现有的研究均围绕实验展开, 很少有对毒品进行过理论拉曼光谱的计算与分析, 且缺乏针对现场执法中常见的毒品痕迹、 疑似毒品、 毒品包装等的拉曼光谱应用研究。 而且, 利用拉曼光谱进行毒品认定往往需要有拉曼数据库, 目前国内大都依赖进口设备及其自带的数据库, 一旦缴获现场出现了未在数据库中的毒品或者毒品混合物, 则可能会导致漏检或误检。 例如, 对于冰毒与咖啡因混合片剂(俗称麻古)[7], 目前尚未见有通过拉曼光谱进行分析检测的相关文献报道。

为此, 本文通过拉曼频率计算(gaussian09), 利用密度泛函理论(DFT), 获得了三种新型毒品样品(冰毒、 K粉、 麻古)的理论拉曼光谱, 并将其与实验光谱进行了比对, 用以进行理论计算的可靠性验证以及特征峰的归属分析, 为今后扩充拉曼数据库或指认新毒品做准备。 在此基础上, 为拓展拉曼光谱在新型毒品快速检测中的应用, 本文结合执法部门现场发现和提取毒品时的特点, 分别以上述三种样品为例, 进行了微量毒品残留物、 疑似毒品(N-异丙基苄胺)[8]等情况下的拉曼光谱检测比对。 考虑到在毒品物证的检测过程中毒品往往会被透明包装包裹, 打开包装检测, 不仅会增加样品准备时间甚至会污染物证或使物证损失, 因此采用透过聚丙烯透明塑料包装对上述毒品样品进行了检测。 这些工作将为未来新型毒品的拉曼数据库完善、 执法部门现场认定以及新型拉曼光谱毒品检测仪研制提供有益的参考。

1 实验部分
1.1 样品来源

所有毒品样品均由浙江省公安厅提供, 为现场缴获所得, 其中冰毒为透明晶体状, K粉为白色粉末状, 麻古为红色片剂; 而非毒品样品N-异丙基苄胺(冰毒同分异构体)透明晶体则从市场直接购得。

1.2 仪器及参数

实验采用由上海如海光电科技有限公司生产的SEED3000型拉曼光谱仪, 其主要技术指标: 光谱采集范围为200~3 200 cm-1, 激光光源的激发波长为785 nm, 激光输出功率最大为500 mW, 积分时间范围为1 ms~20 s。

1.3 方法

取质量为1.00 g的常量冰毒、 K粉、 麻古、 N-异丙基苄胺样品各5份进行激光拉曼光谱测试, 其中积分时间均为10 s, 激光器功率为320 mW。 根据对光谱仪以往测试结果的分析, 为保证置信度大于95%[9], 对每份样品选择不同部位做10次以上的拉曼光谱采集, 并将经过最小二乘法处理后的结果作为本文的实验拉曼光谱谱图。 鉴于毒品缴获现场实时检测的需要, 又取5.0 mg微量的冰毒和K粉并采用上述统计方法进行了毒品残留物拉曼光谱分析, 而聚丙烯透明包装对检测结果可靠性的影响也通过拉曼光谱分析得到。

考虑到实际采集到的初始拉曼数据易受噪声和荧光等多种因素的干扰, 会产生一个倾斜背景, 使基线偏离, 信噪比下降, 影响数据的分析处理, 为此本文采用了自适应迭代重加权惩罚最小二乘法(air-PLS)算法进行背景扣除。 该方法通过迭代重加权惩罚最小二乘法逼近背景, 从原始信号中减掉拟合背景, 得到扣除背景的光谱信号[10]。 图1即为冰毒样品在背景扣除前后的拉曼光谱对比图, 其中红色谱线为未扣除背景的拉曼光谱, 黑色则为扣除背景后的拉曼光谱。 由图可见, 扣除背景并不影响各拉曼峰的峰形和峰位, 但可以有效滤除高频噪声, 便于检测分析。

图1 冰毒样品在背景扣除前后的拉曼光谱对比Fig.1 Raman spectra of methamphetamine before and after baseline removed

1.4 理论计算

各样品的理论拉曼光谱由Gaussian view 5.0进行分子构型, 通过Gaussian09软件, 应用DFT中的Becke-3-Lee-Yang-Parr (B3LYP)杂化泛函, 在6-31G基组下进行结构优化以及拉曼频率计算, 求得能量对坐标的二阶导数再除以原子质量, 获得振动频率。 考虑到由计算方法本身的误差所致的理论计算频率与实验值之间总存在差异, 因此往往需要通过校正因子来进行校正。 在本文中, 最终的计算频率通过乘以0.971的校正因子得到[11, 12]

2 结果与讨论
2.1 毒品的分子几何构型

本文中, 冰毒、 K粉和麻古(主要成分为咖啡因)这三种新型毒品的分子几何构型通过B3LYP杂化泛函方法和6-31G基组优化得到。 在优化过程中, 力的最大值(Maximum Force)、 力的均方根(RMS Force)、 最大位移(Maximum Displacement)和均方根位移(RMS Displacement)四个判据都已收敛, 且在此优化结构下计算得到振动频率都没有虚频存在, 表明优化后的分子几何结构是稳定的。

三种新型毒品优化后稳定的分子几何构型如图2所示。 其中, 冰毒分子苯环上的一个氢原子被异丙基和甲胺构成的支链取代, 形成单取代苯类化合物; K粉分子由环己酮上的两个氢原子分别被甲氨基与氯苯基取代形成; 而咖啡因分子则由黄嘌呤上氮原子连接的氢原子被甲基取代而形成。 对于分子结构中每个原子都标有数字编号, 如1C表示1号碳原子, 7H表示7号氢原子。

图2 冰毒(a), K粉(b)和麻古(c)优化后的分子结构Fig.2 Optimized structures of methamphetamine(a), ketamine(b)and magu(c)

2.2 毒品的理论和实验拉曼光谱比对

图3给出了三种毒品的理论和实验拉曼光谱图, 其中图3(a)为冰毒样品, 图3(b)为K粉样品, 图3(c)为麻古样品。 由图可见, 在三种毒品各自对应的拉曼光谱中, 理论与实验峰位基本一致, 表明理论计算结果与实验结果较为吻合, 且实验获得的特征峰位与文献报道相符[4, 6], 从而有效证明了本文理论与实验结果的可靠性。 相比实验光谱, 理论光谱上些许的差别可能源于理论计算过程中只计算了单个分子, 并没有涉及到分子间相互作用的影响。

图3 三种样品的理论与实验拉曼光谱图比对
(a): 冰毒; (b): K粉; (c): 麻古
Fig.3 Comparison on experimental and calculated Raman spectra of the samples
(a): Methamphetamine; (b): Ketamine; (c): Magu

结合理论计算结果和相关文献, 表1分别对实验测得的三种毒品各自的拉曼特征峰做了归属。 图3(a)冰毒样品实验拉曼光谱中存在4个明显振动峰: 837, 1 003, 1 180和1 209 cm-1。 实验与理论计算结果一致, 最强峰都位于1 003 cm-1附近, 归属为苯环的环呼吸振动。 此外, 837 cm-1为NH摇摆振动; 1 180 cm-1为C— N伸缩振动; 1 209 cm-1处归属为CH2, CH和NH摇摆振动的耦合[4]

表1 三种毒品理论计算与实验光谱振动频率对比及归属 Table 1 Comparison of the theoretical and experimental vibration frequencies of three drugs and assignment

图3(b)K粉样品的特征峰主要有598, 659, 1 046, 1 442和1 593 cm-1。 根据相关文献598 cm-1为C— Cl伸缩振动; 659 cm-1处特征峰由苯环上C— H变形振动引起; 1 046 cm-1为苯环呼吸振动; 1 442 cm-1为CH3变形振动与CH2剪切振动的耦合; 1 593 cm-1为C=C不对称伸缩振动[6]

实验检测的麻古样品外形为红色片剂, 主要成分为咖啡因。 图3(c)麻古样品的实验拉曼光谱中仅存在556, 1 329和1 699 cm-1三个特征峰。 由于目前尚未有文献通过拉曼光谱对麻古进行系统的分析, 因此本文结合理论计算对麻古实验拉曼光谱特征峰进行了归属。 根据理论计算结果: 556 cm-1为咖啡因分子呼吸振动; 1 329 cm-1为C— N伸缩和C— N— C反对称伸缩振动; 1 699 cm-1为C=O伸缩振动。

上述三种毒品拉曼光谱的分析结果显示, 各毒品样品的特征峰差异明显: 837和1 003 cm-1为冰毒的两个最明显特征峰; K粉最具有鉴别价值的特征峰则位于463, 659和1 046 cm-1; 而麻古的三个最明显特征峰分别为556, 1 329和1 699 cm-1。 因此, 利用拉曼光谱对不同类型的毒品进行现场鉴别和认定是切实可行的。 同时, 这种通过密度泛函进行拉曼光谱振动频率计算的方法, 以往常被用于对药物分子拉曼光谱的理论计算[13], 而在毒品拉曼光谱计算上的应用则鲜有报道。 因此, 面对新型毒品更新速度快、 种类繁多的现状, 这一拉曼光谱的理论计算方法或将为今后现场实时分析所需的新型毒品拉曼光谱数据库的建立和完善起到重要的作用。

2.3 微量和常量样品的实验拉曼光谱比对

毒品犯罪现场通常会留有毒品痕迹, 但剂量极低, 利用常用的检测手段较难实现现场认定。 拉曼光谱技术快速、 实时、 无损的优点在对毒品残留物的检测上则有着潜在的优势。 因此, 本文对5.0 mg微量的冰毒和K粉残留物均进行了检测, 并将检测得到的拉曼光谱与1.00 g常量样品拉曼光谱进行了对比。

图4(a)和(b)分别为微量和常量冰毒的拉曼光谱。 通过比较我们发现, 微量冰毒与常量冰毒样品的拉曼光谱图基本一致, 出峰位置均为746, 837, 1 003, 1 209, 1 454和1 602 cm-1。 尽管微量的实验样品在检测过程中荧光干扰显著增强, 但用算法扣除后, 其拉曼光谱基本与常量的冰毒拉曼光谱一致。

图4 微量冰毒(a)和常量冰毒(b)实验拉曼光谱对比Fig.4 Comparison of Raman spectra between trace methamphetamine (a) and constant methamphetamine (b)

图5(a)和(b)分别为微量和常量K粉样品的实验拉曼光谱。 通过比较我们发现, 图5(a)中微量K粉拉曼光谱中的特征峰数量明显少于图5(b)中常量K粉样品, 几乎观察不到598, 1 442和1 593 cm-1这三个峰, 但461, 659和1 050 cm-1三个最具鉴别价值的特征峰则依然存在, 即通过这三个特征峰我们仍可以识别出该样品为K粉。

图5 微量K粉(a)和常量K粉(b)实验拉曼光谱对比Fig.5 Comparison of Raman spectra between trace ketamine (a) and constant ketamine (b)

上述实验结果表明, 利用拉曼光谱可有效识别微量的冰毒和K粉残留物, 这将对禁毒一线的现场痕迹检测起到积极的作用。

2.4 聚丙烯透明包装对于拉曼检测的影响

在检测现场缴获毒品物证的过程中, 毒品常常会有透明塑料包装包裹, 打开包装取出样品检测不仅费时, 而且容易污染毒品样品, 甚至造成毒品物证损失。 因此, 采用透过聚丙烯材料制成的透明包装, 分别对其中的冰毒、 K粉和麻古样品进行了拉曼光谱检测。 为方便比较, 本文对聚丙烯透明塑料包装也单独做了拉曼光谱检测。 图6即为聚丙烯透明包装的拉曼光谱图, 其中双峰814和846 cm-1为其最具鉴别价值的特征峰, 归属为碳原子与氢原子之间摇摆振动[14]

图6 聚丙烯透明包装拉曼光谱Fig.6 Raman spectra of polypropylene transparent plastic packaging

图7(a)和(b)分别为有聚丙烯透明包装和无聚丙烯透明包装时检测冰毒样品得到的拉曼光谱图。 通过对比我们发现: 经透明包装后检测到的冰毒拉曼光谱在1 006 cm-1左右为最强峰, 这与无透明包装时直接检测的结果基本一致; 同时, 800 cm-1之前区域以及1 200 cm-1之后区域的特征峰也基本吻合。 最显著的区别主要在原837 cm-1处的特征峰。 显然, 在加了一层聚丙烯透明包装后该峰被810和845 cm-1两个峰所取代, 而在图6聚丙烯透明包装的拉曼光谱中814和846 cm-1恰好为聚丙烯最主要的特征峰。 由此可知, 透过聚丙烯透明包装鉴别冰毒时, 聚丙烯材料虽会对冰毒样品的拉曼光谱产生轻微的干扰, 但并不影响对冰毒的识别。

图7 透明包装包裹的冰毒样品(a)和无包装包裹的冰毒(b)拉曼光谱对比Fig.7 Raman spectra of methamphetamine measured directly through a plastic (polypropylene) packaging (a) and Raman spectra of methamphetamine without the plastic packaging (b)

相比之下, 图8中透过透明包装检测得到的K粉拉曼光谱不仅有461, 659, 1 048, 1 596和1 725 cm-1这几个K粉的特征峰, 还夹杂着813, 842, 1 151和1 460 cm-1几个聚丙烯透明包装的特征峰。 由此可见, 在透过透明包装检测K粉得到的拉曼光谱中, 我们可以同时识别出K粉和聚丙烯两种物质。

图8 透明包装包裹的K粉(a)和无包装包裹的K粉(b)拉曼光谱对比Fig.8 Raman spectra of ketamine measured directly through a plastic(polypropylene) packaging (a)and Raman spectra of ketamine without the plastic packaging(b)

不同于上述两种毒品, 图9中透过透明包装检测得到的麻古拉曼光谱与聚丙烯透明包装自身的拉曼光谱基本一致, 说明实际检测到的是聚丙烯包装而不是麻古。 究其原因可能是由于麻古自身的荧光干扰太大, 聚丙烯相对于麻古更容易被检测到, 麻古的拉曼峰被聚丙烯透明包装拉曼峰所掩盖。 如图10所示, 在对麻古样品进行拉曼光谱检测时, 其本身带有强烈的荧光背景, 因此检测得到的拉曼信号很微弱。 背景扣除前的麻古拉曼光谱几乎观测不到特征峰, 扣除背景后才有557, 1 325和1 699 cm-1三个明显特征峰出现。 由此可见, 透明包装对麻古的拉曼光谱检测具有重要的影响。

图9 透明包装包裹的麻古样品(a)和透明包装(b)拉曼光谱对比Fig.9 Raman spectra of magu measured directly through a plastic(polypropylene) packaging (a)and Raman spectra of the plastic (polypropylene) packaging(b)

图10 麻古样品在背景扣除前后的拉曼光谱对比Fig.10 Raman spectra of magu before and after baseline removed

2.5 利用拉曼光谱鉴别毒品真伪

N-异丙基苄胺为冰毒的同分异构体, 性状与冰毒极其相似, 无论是缉毒执法者或是滥用者都很难通过肉眼辨别。 目前它是化工和药物的中间体, 属于非管制化学品, 无任何兴奋作用, 被毒贩广泛用作冰毒的掺假剂, 甚至于伪装成冰毒出售, 以牟取高额的利润[8]。 图11为冰毒与N-异丙基苄胺的实验拉曼光谱对比图, 其中冰毒在837 cm-1处为NH摇摆振动, 而该峰在N-异丙基苄胺的拉曼光谱中则偏移到了853 cm-1附近。 此外, N-异丙基苄胺拉曼光谱中1 418 cm-1处为CH3摇摆振动, 冰毒在该波数位置并没有特征峰出现, 冰毒CH3摇摆振动位于1 454 cm-1左右。 上述差异表明, 利用拉曼光谱可鉴别毒品的真伪, 这对今后完善毒品拉曼数据库和进行毒品现场分析都具有重要意义。

图11 冰毒与N-异丙基苄胺的拉曼光谱图Fig.11 Raman spectra of methamphetamine and N-benzylisopropylamine

3 结 论

采用拉曼光谱仪检测了冰毒、 k粉、 麻古这三种常见的新型毒品, 并利用密度泛函理论计算了拉曼光谱, 得到的结果与实验数据相符。 拉曼频率的理论计算可以为今后新型毒品实验拉曼光谱特征峰的归属提供参考。 同时实验还发现, 无论是微量毒品残留物, 还是化学结构与冰毒相似的N-异丙基苄胺, 未来都可以借助拉曼光谱技术进行有效识别和区分。 不过, 聚丙烯材料制成的透明包装则易对自身荧光干扰较强的麻古拉曼检测产生显著的影响。 这些研究结果将有助于我们未来进一步展开新型手持、 非接触、 快速与无损型拉曼光谱毒品检测仪的研制工作。

The authors have declared that no competing interests exist.

参考文献
[1] QIU Hui, LE Sheng, WU Feng-bo(邱晖, 乐胜, 吴逢博). Guangdong Chemical Industry(广东化工), 2016, 43(19): 178. [本文引用:1]
[2] LU Xiao-han, YANG Tian(陆笑寒, 杨甜). Studies of Trace Elements and Health(微量元素与健康), 2014, 31(03): 58. [本文引用:1]
[3] ZHAI Wan-feng, ZHANG Chun-shui, GAO Li-sheng(翟晚枫, 张春水, 高利生). Chinese Journal of Drug Abuse Prevention and Treatment(中国药物滥用防治杂志), 2014, 20(5): 270. [本文引用:1]
[4] Triplett J S, Hatfield J A, Kaeff T L. Journal of Forensic Sciences, 2013, 58(6): 1607. [本文引用:3]
[5] WANG Qing-hua, WANG Rao, WANG Yan-jun(王清华, 王尧, 王燕军). Forensic Science and Technology(刑事技术), 2014, (5): 33. [本文引用:1]
[6] Yang Y, Li Z Y, Yamaguchi K. Nanoscale, 2012, 4(8): 2663. [本文引用:3]
[7] MIAO Cui-ying, ZHOU Lu-zhi, QIAN Zun-lei(苗翠英, 周鲁智, 钱尊磊). Journal of People’s Public Security University of China·Science and Technology(中国人民公安大学学报·自然科学版), 2016, 22(3): 34. [本文引用:1]
[8] DI Yu-min, ZHANG Kai, XIAO Nan(邸玉敏, 张凯, 肖楠). Forensic Science and Technology(刑事技术), 2012, (2): 63. [本文引用:2]
[9] XIAO Ming-yao(肖明耀). Error Theory and its Applications(误差理论与应用). Beijing: Metrology Publishing House(北京: 计量出版社), 1985. [本文引用:1]
[10] ZHANG Xun, CHEN Sheng, WU Bo-shi, et al(张逊, 陈胜, 吴博士, ). Chinese Journal of Analytical Chemistry(分析化学), 2016, 44(12): 1846. [本文引用:1]
[11] SU Ya-qiong, WU De-yin, TIAN Zhong-qun(苏亚琼, 吴德印, 田中群). Acta Physico-Chimica Sinica(物理化学学报), 2014, 30(11): 1993. [本文引用:1]
[12] HUANG Shuang-gen, WU Yan, HU Jian-ping, et al(黄双根, 吴燕, 胡建平, ). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2017, 37(1): 135. [本文引用:1]
[13] WANG Lei, LI Lai-cai(王磊, 李来才). Journal of Sichuan Normal University·Natural Science(四川师范大学学报·自然科学版), 2013, 36(01): 92. [本文引用:1]
[14] DONG Kun, RAO Zhi-fan, YANG Xiao-yun, et al(董鹍, 饶之帆, 杨晓云, ). China Plastics Industry(中国塑料工业), 2011, 39(6): 67. [本文引用:1]