三维荧光结合AWRCQLD测量化妆品中三种添加剂
王玉田1, 张艳1,*, 商凤凯1, 张靖卓2, 张慧1, 孙洋洋1, 王选瑞1, 王书涛1
1. 燕山大学河北省测试计量技术及仪器重点实验室, 河北 秦皇岛 066004
2. 天津做票君机器人科技有限公司, 天津 300450
*通讯联系人 e-mail: 2209479401@qq.com

作者简介: 王玉田, 1952年生, 燕山大学电气工程学院教授 e-mail: y.t.wang@163.com

摘要

没食子酸(GAa), 学名为3,4,5-三羟基苯甲酸(分子式为C7H6O5), 通常以水合物的形式存在, 作为一种重要的有机原料, 广泛的存在于植物中。 有研究证明GAa具有抗氧化、 抗炎、 抗肿瘤、 抗病毒、 抗突变等多种作用。 因此GAa常作为抗氧化剂添加于化妆品中。 对羟基苯甲酸(p-HA), 分子式为C7H6O3, 其中的R基为甲基、 乙基、 丙基、 丁基或庚烷基时分别称为对羟基苯甲酸乙酯、 对羟基苯甲酸丙酯、 对羟基苯甲酸丁酯和对羟基苯甲酸庚酯。 p-HA酯类的抗菌性强、 毒性低、 抑菌作用不受pH影响, 因此常添加于化妆品及药物中用作防腐剂。 间苯二酚(RE)又称1, 3苯二酚或间二苯酚(分子式为C6H6O2)。 RE具有杀菌作用, 可作为防腐剂添加于化妆品中。 以没食子酸(GAa)、 对羟基苯甲酸(p-HA)和间苯二酚(RE)三种化妆品常用添加剂为目标分析物, 通过引入第四维—溶剂, 构建四维荧光光谱数据, 使用甲醇(光谱级)、 乙醇(光谱级)、 超纯水分别获得三组实验样本, 三组样本的配置方法与加入药品量相同。 使用FS920稳态荧光光谱仪对样本进行检测, 设置激发波长为210~330 nm, 间隔4 nm记录一个数据; 发射波长为280~480 nm, 间隔2 nm记录一个数据。 初始发射波长总是滞后激发波长10 nm, 由此可消除一级瑞利散射的干扰。 随后使用空白扣除法对初始荧光数据进行预处理, 去除了溶剂的拉曼散射。 最后, 采用核一致诊断法确定待测样本的组分数为3, 使用交替加权残差约束四线性分解(alternating weighted residual constrained quadratic decomposition, AWRCQLD)算法对预处理后的三维荧光光谱数据进行分解。 结果表明, AWRCQLD算法分解得到GAa、 p-HA和RE的激发、 发射光谱图与目标光谱几乎重叠, 能实现光谱重叠严重的GAa、 p-HA和RE的快速定性和定量分析。

关键词: 三维荧光光谱; 交替加权残差约束四线性分解; 没食子酸; 对羟基苯甲酸; 间苯二酚
中图分类号:O657.3 文献标志码:A
Three-Dimensional Fluorescence Combined with AWRCQLD to Measure Three Additives in Cosmetics
WANG Yu-tian1, ZHANG Yan1,*, SHANG Feng-kai1, ZHANG Jing-zhuo2, ZHANG Hui1, SUN Yang-yang1, WANG Xuan-rui1, WANG Shu-tao1
1. Measurement Technology and Instrument Key Lab of Hebei Province, Yanshan University, Qinhuangdao 066004, China
2. Tianjin ZuoPiaoJun Robot Technology Co., Ltd., Tianjin 300450, China
*Corresponding author
Abstract

Gallic acid (GAa), known as 3, 4, 5-trihydroxybenzoic acid (C7H605), usually in the form of hydrates, as an important organic raw material, is widely found in plants. Studies have shown that GAa has many effects such as anti-oxidation, anti-inflammatory, anti-tumor, anti-viral and anti-mutation. Therefore, GAa is often added to cosmetics as an antioxidant. P-Hydroxybenzoic acid (p-HA), whose molecular formula is C7H603, wherein the R group is methyl, ethyl, propyl, butyl or heptyl, respectively, referred to as ethyl p-hydroxybenzoate, Propyl p-hydroxybenzoate, butyl p-hydroxybenzoate and heptyl p-hydroxybenzoate. The p-HA ester has strong antibacterial property, low toxicity, and antibacterial action against pH, so it is often used as a preservative in cosmetics and medicines. Resorcinol (RE) is also known as 1,3 benzenediol or m-diphenol (formula C6H602). RE has a bactericidal action and can be added to cosmetics as a preservative. In this paper, three kinds of cosmetic additives such as gallic acid (GAa), p-hydroxybenzoic acid (p-HA) and resorcinol (RE) were used as target analytes, and four-dimensional fluorescence spectra were constructed by introducing a fourth dimension solvent. Three sets of experimental samples were obtained for methanol (spectral level), ethanol (spectral level), and ultrapure water, and the configuration of the three groups was the same as that of the added drug. Using the FS920 steady-state fluorescence spectrometer (spectral wavelength response range of 200~900 nm, liquid nitrogen refrigeration range of 77~320 K, excitation source power is 450 W, signal-to-noise ratio is 6 000:1) to test the sample, and setting the excitation wavelength to 210~330 nm, a piece of data was recorded at intervals of 4 nm; the emission wavelength was 280~480 nm, and a piece of data was recorded at intervals of 2 nm. The initial emission wavelength always lagged the excitation wavelength by 10 nm, thereby eliminating the interference of the first-order Rayleigh scattering. The initial fluorescence data were then pretreated using a blank subtraction method to remove Raman scattering of the solvent. Finally, the nuclear consistent diagnosis method was used to determine the number of components of the sample to be tested to be 3, and the three-dimensional fluorescence spectrum data after pretreatment were decomposed using the alternating weighted residual constrained quadratic decomposition (AWRCQLD) algorithm. The results showed that the AWRCQLD algorithm decomposes the excitation and emission spectra of GAa, p-HA and RE almost overlapping with the target spectrum, and can achieve rapid qualitative and quantitative analysis of GAa, p-HA and RE with severe spectral overlap. Using the AWRCQLD algorithm to decompose the samples, the average recoveries of GAa, p-HA and RE were 98.3%, 98.4% and 98.1%, respectively, and the root mean square error (RMSEP) was 0.081, 0.111 and 0.001 μg·mL-1. Three-dimensional fluorescence combined with the AWRCQLD algorithm enables rapid detection of GAa, p-HA and RE in cosmetics.

Keyword: Three-dimensional fluorescence spectroscopy; Alternating weighted residual constrained quadratic decomposition; Gallic acid; p-hydroxybenzoic acid; Resorcinol
引言

没食子酸(Gallic acid, GAa)又称五倍子酸, 是一种存在于自然界的多酚类化合物, 作为一种重要的有机原料, 没食子酸广泛的应用于食品、 医疗、 燃料等行业[1]。 同时GAa具有抗菌抗病毒的特性[2], 以GAa为原料合成的酯类化合物可用作抗氧化剂[3]。 对羟基苯甲酸(p-Hydroxybenzoic acid, p-HA)酯类的抗菌性强、 毒性低与间苯二酚(Resorcinol, RE)均可作为防腐剂添加于化妆品中。

目前检测GAa, p-HA和RE的方法主要有高效液相色谱法(HPLC)、 色谱分析法、 毛细管电泳法等。 但部分检测方法限于实验室且不易小型化。 荧光检测法绿色环保, 操作简便, 成本低, 具有较好应用前景。 国内袁志鹰等[4]使用HPLC测量百合中的对香豆酸和没食子酸, 该方法操作简单, 测量精度高, 但分析成本高, 液相色谱仪价格及日常维护费用贵。 国外Garcia Freitas等[5]使用固定在聚酯树脂上结合Cu3(PO4)的碳纤维复合电极对食物中的丁基羟基茴香醚(BHA)、 和2, 6-二叔丁基对甲酚(BHT)进行了检测。 该方法灵敏度高, 检测范围宽, 但选择性差。

本文通过引入第四维— 溶剂, 构建了四维光谱数据, 采用三维荧光结合交替加权残差约束四线性分解(AWRCQLD), 对光谱严重重叠的GAa, p-HA和RE进行测定。

1 理 论
1.1 四线性模型

对于四维数据矩阵Xq, 其数学表示方法类似于三线性成分模型[6, 7], 用标量表示为

xijkl=n=1Nainbjnckndln+eijkli=1, 2, , I, j=1, 2, J, k=1, 2, , K, l=1, 2, , L(1)

其中, xijklXq中元素; ain, bjn, ckndln分别是矩阵中A(I× N), B(J× N), C(K× N), D(L× N)中的一个元素; eijkl是四维残差阵Eq的一个元素; N为四维数据的组分数[8]

1.2 交替加权残差约束四线性分解(AWRCQLD)

AWRCQLD是由付海燕[9]提出的, 基于四线性成分模型的不完全扩展矩阵形式, 得到四个新的目标函数

σ(A)=l=1LX...l-Adiag(d(l))(CB)TF2+λAl=1L(X...l((CBT)+-Adiag(d(l)))WDF2(2)

σ(B)=i=1IXi...-Bdiag(a(i))(DC)TF2+λBi=1I(Xi...((DCT)+-Bdiag(a(i)))WAF2(3)

σ(C)=j=1JX.j..-Cdiag(b(j))(AD)TF2+λCj=1J(X.j..((ADT)+-Cdiag(b(j)))WBF2(4)

σ(D)=k=1KX..k.-Ddiag(c(k))(BA)TF2+λDk=1K(X..k.((BAT)+-Ddiag(c(k)))WCF2(5)

其中权重矩阵WA=diag(sqrt(1./diagm(ATA))), WB=diag(sqrt(1./diagm(BTB))), WC=diag(sqrt(1./diagm(CTC))), WD=diag(sqrt(1./diagm(DTD)))。 λ A, λ B, λ Cλ D为约束系数, 用来平衡四不完全扩展矩阵形式的损失函数和加权残差函数这两部分的特性。 由交替最小二乘原理, 交替优化式(2)— 式(5)中的目标函数即得到A, B, C, D四个等式

A=(l=1LX...l(CB)+λA((CB)T)+WDWD)diag(d(l))·(l=1Ldiag(d(l))((CTC)* (BTB)+λAWDWD)diag(d(l)))+(6)

B=(i=1IXi...(DC)+λB((DC)T)+WAWA)diag(a(i))·(i=1Idiag(a(i))((DTD)* (CTC)+λBWAWA)diag(a(i)))+(7)

C=(j=1JX.j..(AD)+λC((AD)T)+WBWB)diag(b(j))·(j=1Jdiag(b(j))((ATA)* (DTD)+λCWBWB)diag(b(j)))+(8)

D=(k=1KX..k.(BA)+λD((BA)T)+WCWC)diag(c(k))·(k=1Kdiag(c(k))((BTB)* (ATA)+λDWCWC)diag(c(k)))+(9)

同时, 算法AWRCQLD与四维平行因子(four-way PARAFAC)相比, AWRCQLD具有对待测组分数不敏感且收敛速度快的特点; AWRCQLD与交替惩罚四线性分解(APQLD)相比, AWRCQLD具有抗噪声能力强的优点。

2 实验部分
2.1 仪器与参数

实验使用的仪器为Edinburgh Instruments公司的FS920稳态荧光光谱仪, 选用GAa, p-HA和RE为待测样品, 设置激发波长为210~330 nm, 间隔4 nm记录一个数据; 发射波长为280~480 nm, 间隔2 nm记录一个数据。 初始发射波长总是滞后激发波长10 nm, 由此可消除一级瑞利散射的干扰。

2.2 材料与样品

实验所用的待测样品GAa, p-HA和RE为纯度大于99.5%的标准样品; 溶剂为甲醇(光谱级)。

溶液配制: 称取GAa标准样品0.005 g, 称取p-HA和RE标准样品各0.01 g, 用甲醇溶液分别溶解于50 mL的容量瓶中, 获得浓度为100 μ g· mL-1的一级储备液, 并保存于低温(4 ℃)且避光的环境中。 实验时取5 mL GAa, 10 mL p-HA和10 mL RE的一级储备液, 用甲醇溶液稀释至50 mL, 得到浓度为10 μ g· mL-1的标准溶液。 取适量的标准溶液, 用甲醇溶液稀释成不同浓度的8组校正样本和7组待测样本。 其中, 三组样本中加入一定量未经处理的液态化妆品。 所配置溶液的具体浓度见表1

表1 样品配置浓度(μ g· mL-1) Table 1 Samples preparation concentration
2.3 构建四维光谱

分别以甲醇、 乙醇和超纯水作为溶剂配置样本, 其中3组溶液的配置方法与步骤相同, 最终得到一个四维数据矩阵X=58× 69× 3× 15。 其中58和69为激发和发射波长数, 3为使用的溶剂种类数, 15为样本数。

由图1可以看出p-HA在三种不同溶剂中的荧光峰位置大致相同, 由图2可以看出相同浓度时, p-HA溶解于甲醇溶剂中的荧光强度最大, 超纯水中的荧光强度最小, 这可能与p-HA在水中溶解度小有关。

图1 浓度为4 μ g· mL-1的p-HA分别溶解在三种溶剂中的光谱图
(a), (d): 甲醇溶剂; (b), (e): 乙醇溶剂; (c), (f): 超纯水
Fig.1 Spectra of p-HA at a concentration of 4 μ g· mL-1 dissolved in three solvents
a), (d): ; (b), (e): ; (c), (f): 补英文

图2 p-HA在三种不同溶剂中的荧光峰强度对比Fig.2 Comparison of fluorescence peak intensities of p-HA in three different solvents

3 结果与讨论
3.1 光谱数据预处理

由于大多数溶剂都会引起较弱的Raman散射, 而Raman散射的会干扰待测物本身的荧光光谱, 使用空白扣除法可以有效去除Raman散射。

3.2 AWRCQLD分析

采用核一致诊断法来确定待测样本的组分数, 得到组分数为3。 采用AWRCQLD算法对待测样本数据进行分解, 得到激发、 发射光谱分解图, 由图3可以看出预测荧光光谱与目标荧光光谱高度重合, 算法得到的目标浓度与相对浓度的相关系数均大于r=0.99, AWRCQLD算法的预测效果令人满意。

图3 激发与发射光谱分解图Fig.3 Decomposition of excitation and emission spectrum

AWRCQLD算法的预测浓度如表2所示, 三种物质的回收率为93.3%~106.0%, GAa, p-HA和RE的平均回收率分别为98.3%, 98.4%和98.1%。 GAa, p-HA和RE的预测均方根误差(RMSEP)分别为0.081, 0.011 1和0.001 μ g· mL-1。 由表2可以看出AWRCQLD预测浓度效果较好。

表2 使用AWRCQLD算法测得GAa, p-HA和RE的浓度及回收率 Table 2 Determination of the concentration and recovery of GAa, p-HA and GC using the AWRCQLD algorithm
4 结 论

结果表明构建四维光谱数据且利用AWRCQLD算法对荧光数据进行分解的效果较好, 由于三阶校正算法AWRCQLD抗噪声的能力较强, 因此本文没有对样本数据进行去噪处理。 可以看出使用AWRCQLD算法能够实现化妆品中GAa, p-HA和RE的快速、 准确检测。

参考文献
[1] WANG Zhi-de(王之德). Natural Gas Chemical Industry(天然气化工: C1化学与化工), 1995, (1): 39. [本文引用:1]
[2] LI Xiao-ling, CUI Lan, ZHU De-qiu, et al(李肖玲, 崔岚, 祝德秋, ). China Pharmacist(中国药师), 2004, (10): 767. [本文引用:1]
[3] CHANG Lian-ju, ZHANG Zong-he, HUANG Jia-ling, et al(常连举, 张宗和, 黄嘉玲, ). Biomass Chemical Engineering(生物质化学工程), 2010, (4): 48. [本文引用:1]
[4] YUAN Zhi-ying, LIU Xiang-dan, PEI Gang, et al(袁志鹰, 刘湘丹, 裴刚, ). Chinese Journal of Information on Traditional Chinese Medicine(中国中医药信息杂志), 2018, (5): 82. [本文引用:1]
[5] Garcia Freitas K H, Fatibello-Filho O. Talanta, 2010, 81: 1102. [本文引用:1]
[6] Harshman R A, UCLA. Working Papers in Phonetics, 1970, 16(1): 1. [本文引用:1]
[7] Carroll J D, Chang J J. Psychometrika, 1970, 35(3): 283. [本文引用:1]
[8] LI Shu-fang(李淑芳). Structural Analysis and Multidimensional Calibration of Second Order Correction Algorithm for Quantitative Analysis of Complex Dynamic System Drugs( 二阶校正算法结构剖析和多维校正用于复杂动态体系药物定量分析研究). Hunan University(湖南大学), 2010: 62. [本文引用:1]
[9] FU Hai-yan(付海燕). Chemical Pattern Recognition and Multidimensional Calibration and Its Application in Complex System Analysis( 化学模式识别和多维校正方法及其在复杂体系分析中的应用研究). Hunan University(湖南大学), 2010: 74. [本文引用:1]