气相色谱法测定及预测有机磷阻燃剂的蒸气压

王庆芝, 赵洪霞, 王 琰, 谢 晴, 陈景文, 全 燮

(工业生态与环境工程教育部重点实验室, 大连理工大学环境学院, 辽宁 大连 116023)

研究论文

气相色谱法测定及预测有机磷阻燃剂的蒸气压

王庆芝, 赵洪霞*, 王 琰, 谢 晴, 陈景文, 全 燮

(工业生态与环境工程教育部重点实验室, 大连理工大学环境学院, 辽宁 大连 116023)

气相色谱;蒸气压;定量构效关系;有机磷阻燃剂

有机磷阻燃剂(OPFRs)由于其出色的阻燃性被广泛应用到各种商业产品中,例如电子设备、涂料、软垫家具、建筑材料、纺织品等[1-4]。近年来,随着世界各国逐渐禁止对溴系阻燃剂的生产和使用,OPFRs被认为是常用溴代阻燃剂多溴联苯醚(PBDEs)最为合适的替代品。OPFRs的产量和销量也逐年上升[5],仅在2007年,中国对OPFRs的生产使用量就超过了7万吨[6], 2011年中国磷系阻燃剂的生产量达到10万吨,并且预估中国市场的需求量每年将以15%增加[7]。2011年全球有机磷酸酯的消费量达到了50万吨,2015年全球消费量达到68万吨[4,7]。与PBDEs相似,由于OPFRs主要通过物理方式添加于材料中,而不是和高分子材料化学键合,因此很容易通过挥发、磨损等方式释放到环境介质[8,9]。目前,OPFRs在灰尘[10]、沉积质[11]、土壤[12]、空气(北极)[13]、海洋表层水[14]和生物样品[13,15]中均有检出。有研究表明长期暴露于存在OPFRs的环境中能导致人体激素分泌紊乱、生殖功能下降及生长发育不良等结果[16-22]。由于OPFRs使用量大并具有潜在的环境危害性,因此越来越多的环境学者开始关注这类新型有机污染物。

为了更好地了解OPFRs在环境介质中的迁移情况,准确地获取相关的环境理化参数是非常必要的。蒸气压(P，Pa)作为描述有机化合物环境行为的一种重要理化参数,不仅能够用来表征化合物从液体或固体中逃逸的趋势以及颗粒物对该化合物的吸附效应,而且还可以进一步确定该持久性有机污染物的环境迁移和归趋行为。对于低蒸气压的化合物，一般可采用饱和气流法测定其P值；但对于OPFRs,由于其是半挥发或难挥发性物质,因而采用饱和气流法测定时,气相中的绝对浓度较低,只有用高灵敏度的检测手段才能准确测定其P值。尽管OPFRs作为一类新型有机污染物备受关注,然而这类化合物的P值数据较为缺乏。

本文以p,p′-二氯二苯三氯乙烷(p,p′-DDT)为参考化合物，7种多氯联苯(PCBs)为校正化合物,采用气相色谱法(GC)对14种OPFRs在不同温度(T, K)条件下的PGC值进行了测定。本文探讨了温度与logPGC的相关性,同时研究了分子结构(如摩尔体积、取代模式等)对logPGC的影响,并且通过定量构效关系(QSPR)模型探索logPGC与OPFRs相对保留时间(RRT)的关联,采用了最小二乘法建立logPGC与RRT的线性关系模型,从而预测其他没有标准品的OPFRs的logPGC值。本研究将为探究该类化合物的环境迁移行为和定量评估其环境风险提供一定的理论依据。

1 理论部分

根据Bidleman[23]和Wania等[24]提出的方法,利用待测化合物(i)在色谱柱中的保留时间(tR,i，min)与相同T下某一参考化合物(ref)在色谱柱中的保留时间tR,ref的关系,根据公式(1)对不同T下的ln (tR,i/tR,ref)和lnPref进行线性回归分析,采用最小二乘法,即可获得公式(2)中的蒸发热(ΔvapH，kJ/mol)系数ΔvapHi/ΔvapHref和常数项C。通过公式(2),根据参考化合物的蒸气压值,就可以对未知待测化合物的蒸气压进行测定。

ln (tR,i/tR,ref)=[1-(ΔvapHi/ΔvapHref)] lnPref-C

(1)

lnPi=(ΔvapHi/ΔvapHref) lnPref+C

(2)

由于利用GC测定的蒸气压值与饱和气流法测定的蒸气压值并不总是相同的,因此为了保证GC测定蒸气压的准确性,本文选用7种已知饱和蒸气压的PCBs为校正化合物,对GC的测定值和饱和气流法的测定值进行了校正,并根据实验测定的PGC值来预测OPFRs化合物的PL值(文献校正值)。校正曲线如下所示:

logPL=1.067 logPGC+0.065 4，r=0.999

(3)

2 实验部分

2.1仪器与试剂

HP-6890N气相色谱仪(美国Agilent公司),包括电子捕获检测器(ECD)、火焰离子化检测器(FID)和分流/不分流进样系统;Maxi MixⅡ涡旋混合仪(美国Thermo Scientific公司); MS-105万分之一天平(瑞士Mettler公司); Milli-Q纯水机(美国Millipore公司)。

烷基磷酸酯包括磷酸三正丁酯(tri-n-butyl phosphate, TBP)、磷酸三正戊酯(tri-n-pentyl phosphate, TPeP)、磷酸三己基酯(trihexyl phosphate, THP)、磷酸三(丁氧基乙基)酯(tris(2-butoxyethyl)phosphate, TBEP)和磷酸三异辛酯(tris (2-ethylhexyl) phosphate, TEHP)标准品,卤代烷基磷酸酯包括磷酸三(2-氯乙基)酯(tris (2-chloroethyl) phosphate, TCEP)、磷酸三(2-氯丙基)酯(tris (2-chloroisopropyl) phosphate, TCPP)、磷酸三(1,3-二氯异丙基)酯(tris (1,3-dichloro-2-propyl) phosphate, TDCP)和磷酸三(2,3-二溴丙基)酯(tris (2,3-dibromopropyl) phosphate, TRIS)标准品,芳基磷酸酯包括磷酸甲酚二苯酯(cresyl diphenyl phosphate, CDPP)、磷酸三苯酯(triphenyl phosphate, TPP)、三苯基氧化膦(triphenylphosphine oxide, TPPO)、2-乙基己基二苯基磷酸酯(2-ethylhexyl diphenyl phosphate, EHDP)和磷酸三甲苯酯(tricresyl phosphate, TCP)标准品，参考化合物p,p′-DDT，内标物包括d21-TPrP和d15-TPP，以上物质均购自美国AccuStandard公司,纯度均大于95.0%。7种校正化合物PCB28、PCB52、PCB101、PCB138、PCB153、PCB169和PCB180标准品均购自美国Cambridge Isotope Laboratories公司。

2.2分析条件

色谱柱:DB-1MS毛细管色谱柱(1.0 m×0.25 mm, 0.25 μm,美国Agilent公司);进样口温度:260 ℃;分流比:10∶1;载气:N2;柱流量:3.5 mL/min;程序升温:初始温度60 ℃,保持5.0 min,以15 ℃/min升温至290 ℃,保持2.0 min; FID温度:260 ℃; ECD温度:300 ℃;氢气流速:75 mL/min;空气流速:100 mL/min。

预测模型的建立：14种OPFRs采用DB-5MS、DB-1、DB-17MS和DB-XLB毛细管色谱柱(30 m×0.25 mm, 0.25 μm,美国Agilent公司)进行保留时间的测定分析。程序升温条件为初始温度为40 ℃,保持4.0 min;以15 ℃/min升温至190 ℃;以10 ℃/min升温至310 ℃,保持8.0 min。炉温是由热电偶温度计校正完成的。

2.3数据处理

在363.15～473.15 K温度范围内测定14种OPFRs化合物和p,p′-DDT在5～8个恒温条件下的保留时间，并对其进行数据分析。

本文选用p,p′-DDT作为参考化合物,其蒸气压PL与温度的关系如公式(4)所示:

logPL=12.48-4 699.5/T，r=0.998

(4)

根据公式(4)的斜率可计算出p,p′-DDT的蒸发热ΔvapHref。其中，斜率=ΔvapHref/2.303R(R=8.314 J/(mol5K))。因此可计算出p,p′-DDT的ΔvapHref为89.9 kJ/mol。

图 1 363.15～473.15 K下部分OPFRs的ln (tR,i/tR,ref)和ln Pref的线性关系Fig. 1 Plots of ln (tR,i/tR,ref) versus ln Pref for some OPFRs over the temperature range of 363.15-473.15 K OPFRs: organophosphate flame retardants; TDCP: tris (1,3-dichloro-2-propyl) phosphate; CDPP: cresyl diphenyl phosphate; TPP: triphenyl phosphate; TPPO: triphenylphosphine oxide; EHDP: 2-ethylhexyl diphenyl phosphate; TBEP: tris (2-butoxyethyl)phosphate; TEHP: tris (2-ethylhexyl) phosphate; TCP: tricresyl phosphate; TRIS: tris (2,3-dibromopropyl) phosphate; tR,i: retention time of the analyte; tR,ref: retention time of the reference compound; Pref: vapor pressure of the reference compound.

3 结果与讨论

3.1气相色谱法测定PGC

根据Bidleman[23]和Wania等[24]提出的气相色谱法测定化合物蒸气压的方法,以p,p′-DDT为参考化合物,测定了14种OPFRs在不同温度条件下的蒸气压。图1显示了ln (tR,i/tR,ref)和lnPref之间的线性关系,两者之间的线性关系良好,相关系数

均>0.98,标准偏差均<0.1。通过ln (tR,i/tR,ref)和lnPref的线性回归方程,可以得到公式(2)中的ΔvapHi/ΔvapHref系数和常数项C。根据公式(2)和不同温度下p,p′-DDT的蒸气压值,进而计算出14种OPFRs在不同温度下的蒸气压值。表1列出了在温度为298.15 K时采用GC测定的14种OPFRs的蒸气压值。此外,根据公式(2)中已获得的ΔvapHi/ΔvapHref值,可计算出14种OPFRs的蒸发热。由表1可知,在温度为298.15 K时14种OPFRs的蒸气压值范围为6.72×10-7Pa(TRIS)～5.62×10-2Pa(TBP),蒸发热ΔvapHi的范围为74.1 kJ/mol(TBP)～142.2 kJ/mol(THP),大部分OPFRs的蒸发热高于90 kJ/mol。这与之前报道的PBDEs[25]、四环以上的多环芳烃(PAHs)[26,27]、PCBs[26,28]和某些有机氯农药[26]的蒸发热值的范围基本一致。

3.2温度与logPGC的相关性

为了更好地了解温度对OPFR化合物logPGC的影响,先计算出14种OPFRs在温度范围为278.15～308.15 K下的PL值,再将上述各OPFRs的logPL与1/T根据公式(5)进行相关性分析。公式(5)方程如下所示:

logPL=B/T+A

(5)

图2显示的是各等度温度条件下,14种OPFRs的logPL和1 000/T的线性关系,并可计算出公式(5)中的截距A和斜率B(见表1)。PL是作为描述有机污染化合物从液相到气相的一个重要的理化参数,也是持久性有机污染物的基础特性,对于确定持久性有机污染物在环境中的归趋行为具有重要作用[29]。本文中TBP在温度为298.15 K时的PGC值为5.62×10-2Pa,明显高于TPeP(2.05×10-3Pa)和THP(1.58×10-5Pa),说明TBP比TPeP和THP更容易挥发到外界环境中,并进行长距离的大气传输。PL通常会随着温度的改变而不断变化,由图2可推算出,TPeP在温度为308.15 K时的蒸气压值为6.85×10-3Pa,在温度为278.15 K时的蒸气压值仅为1.42×10-4Pa; TCPP在温度为308.15 K时的蒸气压值为3.85×10-3Pa,其在温度为278.15 K时的蒸气压值也仅为5.42×10-5Pa。以上结果说明,随着环境温度的不断升高,TPeP和TCPP将会更加趋向于从土壤和气溶胶基质中挥发到外界环境中去。对于其他OPFRs化合物的分析也观察到了类似的结果。

3.3分子结构对logPGC的影响采用Le Bas的方法[30]计算得到本研究中的14种OPFRs的摩尔体积(VM, cm3/mol)。图3显示了14种OPFRs的logPGC与VM的线性关系,计算得出的线性方程如下所示:

表 1 14种OPFRs的PGC、PL和ΔvapHiTable 1 Vapor pressures PGC, PL and enthalpies of vaporization ΔvapHi of the 14 OPFRs

PGC: mean±SD (n=3);C: constant term in equation (2);B: slope in equation (5);A: intercept in equation (5); TBP: tri-n-butyl phosphate; TCEP: 2-chloroethyl phosphate; TCPP: tris (2-chloroisopropyl) phosphate; TPeP: tri-n-pentyl phosphate; THP: trihexyl phosphate.

图 2 278.15～308.15 K下14种OPFRs的log PGC和1000/T的线性关系Fig. 2 Plots of log PGC versus 1000/T for the 14 OPFRs over the temperature range of 278.15-308.15 K

卤代烷基磷酸酯:

logPGC=-0.028 3VM+6.664，r=0.906

(6)

图 3 14种OPFRs的log PGC与VM的线性关系Fig. 3 Plots of log PGC versus VM for the 14 OPFRsVM: Le Bas molar volume.

芳基磷酸酯:

logPGC=-0.006 4VM-1.592，r=0.566

(7)

烷基磷酸酯:

logPGC=-0.014 0VM+3.002，r=0.881

(8)

从公式(6)～(8)中可以看出,摩尔体积对OPFRs化合物的logPGC有一定的影响。除芳基磷酸酯的logPGC和VM之间基本无线性关系外,其他OPFRs随着摩尔体积的增大,化合物的logPGC均呈下降趋势;当摩尔体积大于300 cm3/mol时,烷基磷酸酯的logPGC高于卤代烷基磷酸酯。

此外,对于烷基磷酸酯来说,TBP的logPGC为-1.25,明显高于THP (-4.80),此结果说明随着碳链长度的增加,该类化合物的logPGC呈现下降趋势,且每增加一个-CH2-基团,相应化合物的logPGC平均下降约0.5个数量级;对于卤代烷基磷酸酯来说,该类化合物的logPGC随着卤素原子取代数目的增加而降低,例如,TCEP的logPGC为-1.61,明显高于TDCP (-3.77)。以上结果表明,取代模式的不同是影响OPFR化合物logPGC的另一个重要因素。

3.4logPGC-RRT预测模型

OPFRs化合物很多,而标准样品较少,且价格昂贵,这就进一步阻碍了综合评估其环境行为和迁移能力。为了有效预测该类化合物的logPGC值,从根本上解决因标准样品昂贵、短缺等问题所产生的对OPFRs研究工作的阻碍，本文通过对14种OPFRs在4种不同极性色谱柱(DB-5MS, DB-1, DB-17MS和DB-XLB)上的相对保留时间(RRT)与蒸气压之间的相关性分析,建立了基于logPGC与RRT之间的QSPR模型,并将模型用于OPFRs化合物蒸气压的预测。

本研究分别对14种OPFRs在4种不同极性色谱柱上的保留时间进行了准确地测定分析,获得的RRT是由每种OPFRs组分的保留时间除以d21-TPrP和d15-TPP保留时间之和得到的。由于色散力、诱导力、偶极力等作用力的影响,14种OPFRs在4种不同极性色谱柱上的保留时间不同,结果见表2。采用逐步回归的方法分别考察了14种OPFRs的PGC与4种色谱柱上的相对保留时间参数之间的相关性,在4个相对保留时间参数中,发现DB-5MS柱和极性稍弱于DB-5MS柱、色散力相对占优势的DB-1柱的RRT(DB-1)、RRT(DB-5MS)与logPGC存在很好的线性关系。因此,建立了在298.15 K时预测logPGC的模型:

logPGC=5.34-7.91 RRT(DB-1)-7.32 RRT(DB-5MS)

(9)

VIPRRT(DB-1)=1.004, VIPRRT(DB-5MS)=0.996 ﹤0.01 SE=0.128, n=10 表 2 14种OPFRs在4种色谱柱上获得的相对保留时间Table 2 Relative retention times (RRTs) of the 14 OPFRs on the four different gas chromatographic columns

图 4 14种OPFRs的log PGC实验值和预测值的相关性Fig. 4 Plot of the experimental log PGC versus the predicted log PGC for the 14 OPFRs

4 结论

本研究采用GC测定了14种OPFRs在几个不同温度条件下的蒸气压,并根据这些蒸气压值,分析出14种OPFRs的logPGC随着温度的增加而减小。此外,研究了分子结构与logPGC的相关性,结果表明分子取代模式、摩尔体积与OPFRs的logPGC具有很好的相关性。通过逐步回归和变量优化,建立了OPFRs的RRT与logPGC之间的QSPR模型,模型对OPFRs的logPGC的预测值与实验值存在很好的线性关系,对外部样本具有良好的预测能力。

[1] Marklund A, Andersson B, Haglund P. Chemosphere, 2003, 53(9): 1137

[2] Stapleton H M, Dodder N G, Offenberg J H, et al. Environ Sci Technol, 2005, 39(4): 925

[3] Shaw S D, Blum A, Weber R, et al. Rev Environ Health, 2010, 25(4): 261

[4] van der Veen I, de Boer J. Chemosphere, 2012, 88(10): 1119

[5] Weil E D, Levchik S V. J Fire Sci, 2008, 26(3): 243

[6] Wang X W, Liu J F, Yin Y G. Progress in Chemistry, 2010, 22(10): 1983

王晓伟, 刘景富, 阴永光. 化学进展, 2010, 22(10): 1983

[7] Ou Y X. Chemical Industry and Engineering Process, 2011, 30(1): 210

欧育湘. 化工进展, 2011, 30(1): 210

[8] Cristale J, Quintana J, Chaler R, et al. J Chromatogr A, 2012, 1241: 1

[9] Wang X W, He Y Q, Lin L, et al. Sci Total Environ, 2014, 470/471: 263

[10] Bacaloni A, Cavaliere C, Foglia P, et al. Rapid Commun Mass Spectr, 2007, 21(7): 1123

[11] Cao S, Zeng X, Song H, et al. Environ Toxicol Chem, 2012, 31(7): 1478

[12] Mihajlovic I, Miloradov M V, Fries E. Environ Sci Technol, 2011, 45(6): 2264

[13] Moeller A, Sturm R, Xie Z, et al. Environ Sci Technol, 2012, 46(6): 3127

[14] Bollmann U E, Moeler A, Xie Z, et al. Water Res, 2012, 46(2): 531

[15] Kim J, Isobe T, Chang K, et al. Environ Pollut, 2011, 159(12): 3653

[16] Farhat A, Crump D, Chiu S, et al. Toxicol Sci, 2013, 134(1): 92

[17] Li H, Su G Y, Zou M, et al. Environ Sci Technol, 2015, 49(21): 12975

[18] Liu C S, Su G Y, Giesy J P, et al. Sci Rep, 2016, 6: 19045

[19] Liu X S, Ji K, Jo A, et al. Aquat Toxicol, 2013, 134/135: 104

[20] McGee S P, Cooper E M, Stapleton H M, et al. Environ Health Perspect, 2012, 120(11): 1585

[21] Meeker J D, Stapleton H M. Environ Health Perspect, 2010, 118(3): 318

[22] Wang J Z, Hou Y Q, Zhang J S, et al. J Hazard Mater, 2013, 263(2): 778

[23] Bidleman T F. Anal Chem, 1984, 56: 2490

[24] Wania F, Lei Y D, Harner T. Anal Chem, 2002, 74(14): 3476

[25] Wong A, Lei Y D, Alaee M, et al. J Chem Eng Data, 2001, 46(2): 239

[26] Hinckley D A, Bidleman T F, Foreman W T. J Chem Eng Data, 1990, 35: 232

[27] Lei Y D, Chankalal R, Chan A, et al. J Chem Eng Data, 2002, 47(4): 801

[28] Puri S, Chickos J S, Welsh W J. Anal Chem, 2001, 73(7): 1480

[29] Rittfeldt L. Anal Chem, 2001, 73(11): 2405

[30] Reid R C, Prausnitz J M, Polling B E. The Properties of Gases and Liquids. 4th ed. New York: McGraw-Hill Irwin Press, 1987

Determination and prediction for vapor pressures of organophosphate flame retardants by gas chromatography

WANG Qingzhi, ZHAO Hongxia*, WANG Yan, XIE Qing, CHEN Jingwen, QUAN Xie

(Key Laboratory of Industrial Ecology and Environmental Engineering of Ministry of Education, School of Environmental Science and Technology, Dalian University of Technology, Dalian 116023, China)

gas chromatography (GC); vapor pressure (P); quantitative structure-property relationship (QSPR); organophosphate flame retardants (OPFRs)

10.3724/SP.J.1123.2017.05015

2017-05-18

.E-mail:hxzhao@dlut.edu.cn.

国家自然科学基金(21677023);城市水资源与水环境国家重点实验室开放基金项目(QAK201606).

O658

A

1000-8713(2017)09-1008-06

Foundation item: National Natural Science Foundation of China (No. 21677023); State Key Laboratory Open Foundation of Urban Water Resource and Environment (No. QAK201606).

*