李敏雯,肖少軍,張欽發(fā),*,向　紅

(1.無(wú)限極(中國(guó))有限公司,廣東廣州 510623;2.華南農(nóng)業(yè)大學(xué)食品學(xué)院,廣東廣州 510640)

?

Scatchard-Hildebrand熱力學(xué)模型在預(yù)測(cè)PE中有害物質(zhì)遷移分配系數(shù)的應(yīng)用

李敏雯1,2,肖少軍2,張欽發(fā)2,*,向紅2

(1.無(wú)限極(中國(guó))有限公司,廣東廣州 510623;2.華南農(nóng)業(yè)大學(xué)食品學(xué)院,廣東廣州 510640)

Scatchard-Hildebrand 熱力學(xué)模型，PE，分配系數(shù)，預(yù)測(cè)

包裝中有害物質(zhì)向食品中遷移一直是人們關(guān)注的重點(diǎn),而遷移測(cè)試時(shí)間長(zhǎng)以及對(duì)設(shè)備儀器的要求高,遷移檢測(cè)方法又不具普適性,而利用遷移模型來(lái)預(yù)測(cè)塑料中有害物質(zhì)的遷移情況成為目前及未來(lái)的研究熱點(diǎn)[1]。遷移物分配系數(shù)客觀上決定了塑料材料中有害化學(xué)物進(jìn)入食品的最大遷移量,是衡量包材對(duì)食品安全造成影響程度的關(guān)鍵參數(shù)。目前對(duì)分配系數(shù)的估算模型有正規(guī)溶液理論和基團(tuán)貢獻(xiàn)法[2],但這兩種模型都需要大量計(jì)算,并且缺乏推廣性。Scatchard-Hildebrand統(tǒng)計(jì)熱力學(xué)模型是根據(jù)多相體系中內(nèi)能的平衡而建立的預(yù)測(cè)分配系數(shù)的模型,在金屬溶劑有機(jī)萃取分配系數(shù)的預(yù)測(cè)中得到廣泛的應(yīng)用,具有良好的預(yù)測(cè)效果,而食品包裝材料中有害物質(zhì)向食品中的遷移與金屬萃取在原理上是相近的[3-5]。研究表明[6]Scatchard-Hildebrand 熱力學(xué)模型在PE中抗氧劑向脂肪類食品中遷移的分配系數(shù)預(yù)測(cè)過(guò)程中,預(yù)測(cè)結(jié)果與實(shí)驗(yàn)結(jié)果雖存在差距,但模型對(duì)某些特定的遷移體系還是具有一定的適用性,這主要是由于模型只是考慮了塑料、遷移物和模擬液的相容性差異,忽略了實(shí)際遷移過(guò)程中諸如傳質(zhì)阻力、塑料溶脹等因素對(duì)分配系數(shù)造成的影響。因此,需要對(duì)Scatchard-Hildebrand 熱力學(xué)模型進(jìn)行實(shí)驗(yàn)數(shù)據(jù)修正,使其對(duì)塑料中有害物質(zhì)向食品遷移的分配系數(shù)預(yù)測(cè)更加準(zhǔn)確。

表1　物質(zhì)特征參數(shù)表Table 1　Material characteristic parameters list

本文利用Scatchard-Hildebrand熱力學(xué)模型對(duì)PE材料內(nèi)五種抗氧劑向不同脂肪類模擬液中遷移平衡分配系數(shù)進(jìn)行預(yù)測(cè),并與文獻(xiàn)中的實(shí)驗(yàn)測(cè)試結(jié)果進(jìn)行對(duì)比分析,對(duì)熱力學(xué)模型進(jìn)行修正,以期對(duì)食品安全包裝材料的生產(chǎn)和選擇提供指導(dǎo)和幫助。

1　材料與方法

1.1研究材料

表1中物質(zhì)的醇-水分配系數(shù)利用ChemOffice Ultra 2004軟件(美國(guó)劍橋公司)進(jìn)行分子模擬得到,其他各物質(zhì)特征參數(shù)數(shù)據(jù)則來(lái)源于文獻(xiàn)[7-9]。

1.2實(shí)驗(yàn)方法

按照首先選取經(jīng)典成熟的測(cè)試方法[10]得出的數(shù)據(jù),其次盡量選擇同一實(shí)驗(yàn)室(或課題組)測(cè)試得出數(shù)據(jù)的篩選原則,從大量的文獻(xiàn)中搜集并選取了文獻(xiàn)[11-14]中塑料PE內(nèi)五種抗氧劑,即BHA、BHT、Irganox1076、Irganox1010以及Irgafos168在不同溫度下向脂肪模擬液異辛烷、正己烷及95%乙醇中遷移平衡時(shí)的實(shí)驗(yàn)測(cè)試結(jié)果,其中文獻(xiàn)[13-15]中分配系數(shù)數(shù)據(jù)Fexp被直接引用,文獻(xiàn)[11-12]中的結(jié)果圖形通過(guò)圖像數(shù)據(jù)化工具軟件GetData Graph Digitizer 2.24[16]取得遷移結(jié)果數(shù)據(jù),再通過(guò)已知條件轉(zhuǎn)化計(jì)算得到分配系數(shù)實(shí)驗(yàn)值Fexp,同時(shí)將表1中的相關(guān)物質(zhì)參數(shù)代入Scatchard-Hildebrand熱力學(xué)模型公式如下[6]:

式(1)

其中,Fsim為遷移物的模擬分配系數(shù),V遷移物為遷移物的摩爾體積/(cm3/mol),δ聚合物、δ模擬液、δ遷移物分別為聚合物、模擬液和遷移物的溶解度參數(shù)((J/cm3)0.5),φ模擬液、φ聚合物分別為模擬液在模擬液/遷移物混合體系中、聚合物在聚合物/遷移物混合物體系中所占的體積分?jǐn)?shù),因?yàn)檫w移物在各混合體系中所占體積都非常小,所以φ模擬液、φ聚合物一般都設(shè)為1,R為氣體常數(shù)(8.314 J·mol-1·K-1),T為溫度(K)。

計(jì)算得出分配系數(shù)模擬值Fsim,由于模型是在諸多假設(shè)下成立的[6],與實(shí)際遷移結(jié)果有一定的差距,需要對(duì)模型進(jìn)行修正,利用實(shí)驗(yàn)值對(duì)模擬值進(jìn)行修正。為方便計(jì)算,將實(shí)驗(yàn)值Fexp取對(duì)數(shù),然后將模擬值對(duì)數(shù)與實(shí)驗(yàn)值對(duì)數(shù)相減得到修正系數(shù)Ks(即Ks=ln(Fsim/Fexp)),研究修正系數(shù)Ks與體系溫度(T)、遷移物分子量(M)、遷移物醇-水分配系數(shù)(logP)以及模擬液極性Ps間的數(shù)學(xué)關(guān)系,最后通過(guò)軟件SPSS16.0(SPSS公司)利用多元線性回歸分析所有參數(shù)(T、M、logP、Ps)與修正系數(shù)Ks的關(guān)系,得到最終修正系數(shù)方程,即:Ks=ln(Fsim/Fexp)=f(T,M,logP,Ps),其中f(T,M,logP,Ps)表示與溫度、遷移物分子量和醇-水分配系數(shù)以及模擬液極性相關(guān)的函數(shù)。那么,遷移平衡時(shí)分配系數(shù)就可以通過(guò)修正后的如下模型公式(2)進(jìn)行預(yù)測(cè)。

式(2)

2　結(jié)果與分析

2.1計(jì)算結(jié)果

PE中有害物質(zhì)向食品模擬物中遷移的實(shí)驗(yàn)測(cè)定值Fexp和采用Scatchard-Hildebrand熱力學(xué)模型計(jì)算值Fsim以及模型修正系數(shù)Ks列于表2。

表2表明,PE中不同的遷移物向不同的食品模擬物中在不同的溫度條件下遷移的分配系數(shù)的實(shí)驗(yàn)值Fexp與Scatchard-Hildebrand熱力學(xué)模型計(jì)算值Fsim比值都不同,特別是隨著溫度、遷移物的種類以及模擬物的極性等不同,模擬值與實(shí)驗(yàn)值比值,即Ks值有明顯的差異,這可能是因?yàn)橛?jì)算值是在諸多假設(shè)條件下得出的,而實(shí)際遷移中分配系數(shù)受諸如傳質(zhì)阻力、塑料溶脹等實(shí)際存在因素的影響[6],研究Ks值與這些參數(shù)的關(guān)系,實(shí)現(xiàn)對(duì)Scatchard-Hildebrand熱力學(xué)模型的修正,從而使該熱力學(xué)模型可以進(jìn)行準(zhǔn)確的預(yù)測(cè)。

表2　PE材料中遷移物分配系數(shù)實(shí)驗(yàn)值Fexp與模擬值Fsim及修正系數(shù)KsTable 2　Partition coefficient experimental values Fexp, simulated values Fsim and correction coefficients Ks of migrants in PE material

2.2單因素實(shí)驗(yàn)結(jié)果與分析

2.2.1溫度對(duì)Ks的影響根據(jù)表2文獻(xiàn)[11-12]中的實(shí)驗(yàn)測(cè)試數(shù)據(jù),探討其他條件相同時(shí)溫度對(duì)修正系數(shù)Ks的影響,結(jié)果如圖1。

圖1　PE中BHT、BHA、Irganox1010、 Irganox1076向異辛烷遷移時(shí)溫度對(duì)Ks的影響Fig.1　Effect of temperature on Ks of BHT， BHA，Irganox1010 and Irganox1076 in PE migrate to isooctane

圖1顯示,修正系數(shù)Ks隨溫度升高都逐漸減小,Irganox1010相關(guān)系數(shù)只有0.875,其中原因一方面可能是軟件采集實(shí)驗(yàn)數(shù)據(jù)時(shí)出現(xiàn)的誤差,另一方面可能是Irganox1010分子結(jié)構(gòu)過(guò)于龐大,溫度造成其分子在聚合物間運(yùn)動(dòng)相比其他遷移物要困難[12],而其他遷移物的相關(guān)系數(shù)都達(dá)到0.97以上,說(shuō)明溫度與Ks有較好的線性關(guān)系,而Ks與溫度呈負(fù)相關(guān)是由于無(wú)論實(shí)驗(yàn)分配系數(shù)還是模型預(yù)測(cè)分配系數(shù),都會(huì)因溫度升高而降低,但是所采用的分配系數(shù)預(yù)測(cè)模型并未考慮諸如分子間吸引力、摩擦力等阻力對(duì)實(shí)際遷移造成的影響[6]。

2.2.2遷移物分子量對(duì)Ks的影響根據(jù)表1和表2,探討其他條件相同時(shí)遷移物分子量對(duì)修正系數(shù)Ks的影響,其結(jié)果如圖2。

圖2　PE中BHT、BHA、Irganox1010、Irganox1076 在相同溫度下向異辛烷遷移時(shí)遷移物分子量對(duì)Ks的影響Fig.2　Effect of migration molecular weight on Ks of BHT，BHA，Irganox1010 and Irganox1076 in PE migrate to isooctane

由圖2可以看出,Ks也隨著遷移物分子量增大而增大,線性相關(guān)系數(shù)達(dá)到0.994以上,說(shuō)明遷移物分子量與修正系數(shù)Ks間存在良好的線性關(guān)系。這是由于所建分配系數(shù)預(yù)測(cè)模型只是從熱力學(xué)角度考慮了遷移物與聚合物、食品(模擬液)間的相溶性,一般遷移物分子量越大,其空間結(jié)構(gòu)也就越龐大復(fù)雜,在實(shí)際遷移過(guò)程中所受的阻力越大,也越容易滯留于塑料材料中導(dǎo)致最終分配系數(shù)越大[12]。

2.2.3遷移物醇-水分配系數(shù)對(duì)Ks的影響醇-水分配系數(shù)logP反映的是物質(zhì)在油水兩相中的分配情況,logP越大,該物質(zhì)越親油,越小則水溶性越好。根據(jù)表1和表2探討其他條件相同時(shí)遷移物醇-水分配系數(shù)對(duì)修正系數(shù)Ks的影響,其結(jié)果如圖3。

圖3　PE中BHT、BHA、Irganox1010、Irganox1076 在相同溫度下向異辛烷遷移平衡時(shí) 遷移物醇-水分配系數(shù)對(duì)Ks的影響Fig.3　Effect of migration alcohol/water partition coefficient on Ks of BHT，BHA，Irganox1010 and Irganox1076 in PE migrate to isooctane

結(jié)合表1,由圖3可見(jiàn),除了抗氧劑1010(logP=10.571)外,Ks值基本隨遷移物的logP值增大而呈線性增長(zhǎng),這主要是由于模擬液會(huì)對(duì)PE產(chǎn)生溶脹作用,根據(jù)相似相容原理,其在溶脹過(guò)程中對(duì)親油的遷移物分子有更好的溶解效果,從而實(shí)際遷移中越親油的遷移物分子越易遷出[17],實(shí)際遷移實(shí)驗(yàn)中分配系數(shù)Fexp越小,并且抗氧劑1010因其呈幾何對(duì)稱樹枝狀的分子結(jié)構(gòu)使其實(shí)際分配系數(shù)相比更加小[18]。

2.2.4模擬液極性對(duì)Ks的影響根據(jù)表1和表2,探討其他條件相同時(shí)模擬液極性對(duì)修正系數(shù)Ks的影響,其結(jié)果如圖4。

圖4　模擬液極性值對(duì)Ks的影響Fig.4　Effect of simulation liquid polarity on Ks

由圖4可知,無(wú)論何種遷移物,Ks隨模擬液極性對(duì)數(shù)值增大而線性增大,由于遷移物大多都是極性物質(zhì),由相似相容原理,隨著模擬液極性增大,遷移物在實(shí)際遷移過(guò)程中就越易溶解,那么實(shí)驗(yàn)分配系數(shù)就相比更小,模擬分配系數(shù)與實(shí)驗(yàn)分配系數(shù)的比值Ks就會(huì)越大。

2.3多元線性回歸分析法篩選Ks的最優(yōu)參數(shù)及驗(yàn)證

為研究以上各種單因素對(duì)Ks影響的顯著性及其綜合影響效果,根據(jù)表1中物質(zhì)參數(shù)及表2中所有的實(shí)驗(yàn)結(jié)果,在軟件SPSS中利用多元線性回歸(MLR)對(duì)溫度(T)、遷移物分子量(M)、遷移物醇-水分配系數(shù)(logP)、模擬液極性(Ps)與修正系數(shù)Ks進(jìn)行多元線性回歸,并采用逐步回歸(SLR)篩選出三個(gè)最優(yōu)參數(shù)(M、logP、Ps),得到如下修正方程:

Ks=7.289+0.021·M+1.317·lnPs-0.901·logP,其中N=33;調(diào)整R2=0.902;F=99.122;sig=0.000。

多元線性回歸分析結(jié)果表明:M、logP、lnPs和校正值Ks具有顯著的相關(guān)性,而T與Ks不具有明顯的相關(guān)性。那么,根據(jù)公式(2)可得修正后的分配系數(shù)預(yù)測(cè)模型如下:

式(3)

圖5　PE中抗氧劑向脂肪模擬液遷移分配系數(shù)實(shí)驗(yàn)值Fexp與預(yù)測(cè)值對(duì)比圖Fig.5　Contrast of experimental partition coefficient value Fexpand predicted one of antioxidants migrate from PE to fat simulants

3　結(jié)論

[1]Helmroth E,Rijk R,Dekker M,et al. Predictive modeling of migration from packaging materials into food products for regulatory purposes[J]. Trends in Food Science & Technology,2002,13:102-109.

[2]皮林格,巴納主編. 范家起,張玉霞譯. 食品用塑料包裝材料——阻隔功能、傳質(zhì)、品質(zhì)保證和立法[M]. 北京:化學(xué)工業(yè)出版社,2004.

[3]李以圭. 金屬溶劑萃取熱力學(xué)[M]. 北京:清華大學(xué)出版社,1988:145-152.

[4]Gillet G,Vitrac O,Desobry S. Prediction of partition coefficients of plastic additives between packaging materials and food simulants[J]. Industrial & Engineering Chemistry Research,2010,49(16):7263-7280.

[5]Piringer OG,Baner AL. Plastic packaging:interactions with food and pharmaceuticals[M]. Second,Completely Revised Edition. Weinheim:Wiley-VCH,2008. 504.

[6]肖少軍,張欽發(fā),向紅,等. 熱力學(xué)數(shù)學(xué)模型對(duì)塑料中化學(xué)物向食品模擬液中遷移平衡時(shí)分配系數(shù)的預(yù)測(cè)研究[J].食品安全質(zhì)量檢測(cè)學(xué)報(bào),2014,5(6):1746-1750.

[7]項(xiàng)文娟,郭小然,韓靜.溶解度參數(shù)在固體分散體中的應(yīng)用[J].今日藥學(xué),2010,20(12)：5-8.

[8]Hansen CM. Hansen solubility parameters-A User’s Handbook[M]. New York:CRC Press,2000.

[9]Tehrany E A,Brandsch J,Mercea P,et al. Simple method to calculate partition coefficient of migrant in food simulant/polymer system[J]. Journal of Food Engineering,2006,77:135-139.

[10]李揮,墨偉,范斌,等.食品包裝材料遷移實(shí)驗(yàn)方法和食品模擬物選擇[J].塑料助劑,2008,2:13-17.

[11]林勤保,李小梅,宋歡,等.聚乙烯塑料包裝材料中4中抗氧化劑向脂肪食品模擬物遷移的研究[J].分析科學(xué)學(xué)報(bào),2010,26(6):631-635.

[12]李小梅.塑料與紙質(zhì)食品包裝材料中抗氧化劑的測(cè)定及遷移研究[D].大同:山西大學(xué),2010.

[13]胡長(zhǎng)鷹,黃肖紅.HDPE膜中抗氧劑1076向替代脂肪類食品模擬物遷移的研究[J].食品工業(yè)科技,2009，30(2):244-247.

[14]黃肖紅.HDPE膜中抗氧化劑1076向食品或食品模擬物的遷移研究[D].廣州:暨南大學(xué),2008.

[15]趙威威.塑料包裝材料中化學(xué)物向脂肪類食品的遷移研究[D].無(wú)錫:江南大學(xué),2007.

[16]http://getdata-graph-digitizer.com/registration.php

[17]Ozaki A,Gruner A,Stoermer A,et al. Correlation between partition coefficientspolymer/food simulant,KP,Fand octanol/water LogPOW-a new approach in support of migration modeling and compliance testing[J]. Deutsche Lebensmittel-Rundschau,2010,106(4):203-208.

[18]邸曉賀.星型抗氧劑的設(shè)計(jì)合成與性能研究[D].大慶：大慶石油學(xué)院,2007:18-20.

[19]黃秀玲,王志偉,李明.影響PE淋膜紙中化學(xué)遷移物分配系數(shù)的因素[J]. 高分子材料科學(xué)與工程，2010,26(8):23-26.

Application of predicting migration partition coefficient of harmful migrants in PE by Scatchard-Hildebrand thermodynamic modeling

LI Min-wen1，2，XIAO Shao-jun2，ZHANG Qin-fa2，*，XIANG Hong2

(1.Infinitus(China) Co.，LTD，Guangzhou 510623，China; 2.Food Science College of South China Agricultural University，Guangzhou 510642，China)

2015-07-13

李敏雯(1973-)，女，在讀碩士研究生，主要從事日化和保健品新包裝開(kāi)發(fā)研究，E-mail：Cathy.Lee@infinitus-int.com。

張欽發(fā)(1963-)，男，博士，教授，研究方向:包裝新材料與新技術(shù)，E-mail：zqfzgn@163.com。

國(guó)家自然科學(xué)基金項(xiàng)目(31171689)。

10.13386/j.issn1002-0306.2016.05.009