王程輝，胡永華*，王健，徐迎波，潘洋

1.安徽中煙工業(yè)有限責(zé)任公司技術(shù)中心，合肥市高新區(qū)天達路9號230088

2.中國科學(xué)技術(shù)大學(xué)國家同步輻射實驗室，合肥市合作化南路42號230029

新鮮的卷煙主流煙氣是一種高活性的氣溶膠，其物質(zhì)組成處在一個不斷變化的動態(tài)過程中［1-3］。傳統(tǒng)的卷煙煙氣分析方法由于要涉及到取樣、分離、純化以及衍生化等復(fù)雜的分析過程，會導(dǎo)致煙氣的陳化，而且在取樣和分析過程中還可能會引入人為的后生物質(zhì)，從而使其分析結(jié)果不能夠真實地反映吸煙者口腔中的煙氣狀況［4-5］。近年來，對卷煙煙氣進行在線分析研究日益受到了人們的關(guān)注和重視［5-6］。離子-分子反應(yīng)和真空紫外光電離兩種軟電離質(zhì)譜技術(shù)具有速度響應(yīng)快、高選擇性和適合多組分檢測等優(yōu)點，已成為卷煙煙氣在線和實時分析的有效手段［5，7-8］。這類軟電離質(zhì)譜技術(shù)除了依據(jù)待測物的分子量進行初步定性外，還需要結(jié)合電離能作進一步的定性分析，然而，目前有很多物質(zhì)在文獻中找不到精確的電離能數(shù)據(jù)，這給定性分析研究帶來了困難。

物質(zhì)的電離能可以通過實驗方法直接測量，也可以根據(jù)分子軌道理論計算得到。實驗測定雖然能得到可靠的結(jié)果，但是由于化合物數(shù)量眾多，直接對所有化合物的電離能進行測定，工作量大，難度較高；理論計算不受實驗條件的限制，是獲得物質(zhì)電離能的另一可行方法。目前，文獻報道中關(guān)于電離能的計算方法較多［9-11］，但這些方法主要適用于特定類型化合物電離能的計算，普適性不強。因此，探討能夠精確計算電離能的普適性方法并考察其計算結(jié)果在卷煙煙氣組分定性分析中的應(yīng)用，具有重要的理論研究和應(yīng)用價值。

1 材料與方法

1.1 材料與儀器

實驗所使用的材料為1種國產(chǎn)烤煙型卷煙；實驗儀器包括1個自制的卷煙煙氣采樣器和1臺自制的反射式飛行時間質(zhì)譜儀。

1.2 物質(zhì)分子的電離能計算方法

先用ChemOffice2002軟件繪制分子的初始結(jié)構(gòu)，并用分子力學(xué)方法MM2進行預(yù)優(yōu)化；然后在Gaussian 03軟件包中分別采用密度泛函方法［B3LYP/6-311++G（d，p）］和組合從頭算方法（G3B3）對分子離子和中性分子結(jié)構(gòu)作進一步優(yōu)化和能量計算，并且每一個最優(yōu)構(gòu)型都經(jīng)過頻率計算確認(rèn)不存在任何虛頻，是勢能面上的一個極小點。計算完成后，即可獲得分子離子和中性分子的能量，物質(zhì)分子的電離能為這兩種能量的差值［12］。

1.3 卷煙煙氣氣相化學(xué)成分的在線分析

圖1 卷煙煙氣氣相化學(xué)成分在線分析實驗裝置示意圖

卷煙煙氣氣相化學(xué)成分的在線分析實驗裝置示意圖見圖1。該裝置主要由1個卷煙煙氣采樣器和1臺反射式飛行時間質(zhì)譜儀組成。煙氣采樣器的主體部分是1個三端口注射器，3個端口分別是：通過閥a與卷煙夾持器相連接的煙氣吸入端口，過閥b與煙氣傳輸石英毛細(xì)管相連的煙氣輸出端口和活塞端口。反射式飛行時間質(zhì)譜儀的前端設(shè)置有同步輻射真空紫外光電離區(qū)，真空紫外光由800 MeV電子儲環(huán)的彎轉(zhuǎn)磁鐵產(chǎn)生。卷煙煙氣由質(zhì)譜內(nèi)真空負(fù)壓產(chǎn)生的吸力經(jīng)過石英毛細(xì)管引入電離區(qū)，并與真空紫外光交叉后被電離，離子在質(zhì)譜中飛行一段時間后被多通道探測器（MCP）檢測。

測量時，先將劍橋濾片放入卷煙夾持器中，并在夾持器的前端插入煙支，此時閥a開啟，閥b關(guān)閉。然后，抽吸卷煙，主流煙氣中的粒相物通過劍橋濾片攔截，氣相物進入采樣器。一次抽吸完畢后，立即關(guān)閉閥a，打開閥b，與此同時，啟動光電離質(zhì)譜進行檢測。質(zhì)譜檢測采用兩種掃描模式：①固定光子的能量進行掃描，可以獲得待測物質(zhì)的光電離質(zhì)譜圖；②掃描光子的波長，測量在不同波長下的光電離質(zhì)譜，通過積分質(zhì)譜峰的面積，可以得到不同質(zhì)量的離子隨波長（即能量）變化的曲線，該曲線稱之為光電離效率譜曲線。

2 結(jié)果與討論

2.1 兩種方法計算卷煙煙氣化學(xué)成分電離能的可靠性評估

為了評估B3LYP/6-311++G（d，p）方法和G3B3方法的電離能計算結(jié)果的可靠性，選取中國煙草科學(xué)數(shù)據(jù)庫——卷煙煙氣分析數(shù)據(jù)分庫［13］中分子量在100 Da以內(nèi)且文獻［14］中給出了精確電離能數(shù)據(jù)的70種煙氣化學(xué)組分作為驗證集進行計算和分析。表1中列出了這70種物質(zhì)電離能的文獻值、計算值以及二者之間的偏差。

表1 70種卷煙煙氣化學(xué)成分電離能的文獻值、計算值及其偏差（eV）

(續(xù)表1)

由表1可以看出，這70種驗證集組分包括了單質(zhì)以及烴類、醇類、酚類、醚類、醛類、酮類和雜環(huán)類化合物等，基本上覆蓋了常見煙氣化學(xué)成分的類別，表明所選的物質(zhì)具有較好的代表性。雖然方差分析結(jié)果顯示這兩組電離能的計算值與文獻值之間均沒有表現(xiàn)出顯著的差異性，但是與文獻中給出的精確電離能數(shù)值相比較，通過B3LYP/6-311++G（d，p）方法計算得到的電離能數(shù)值中，除了少數(shù)物質(zhì)電離能的偏差較小外，大多數(shù)電離能的偏差均大于0.10 eV，有的偏差甚至高于1.00 eV；而采用G3B3方法計算得到的電離能中，絕大多數(shù)物質(zhì)的電離能均能夠與文獻值很好的吻合。這說明G3B3方法計算得到的電離能的精確度要高于B3LYP/6-311++G（d，p）方法。

圖2為所有驗證集合電離能的理論值和文獻值的線性擬合關(guān)系。由圖2（a）中可知，B3LYP/6-311++G（d，p）方法的擬合相關(guān)系數(shù)為0.9868，標(biāo)準(zhǔn)偏差（SD）為0.20 eV，這進一步表明該方法在預(yù)測不同類型化學(xué)組分的電離能上存在較大的偏差。而從圖2（b）可以看出，G3B3方法計算的電離能與文獻值線性擬合的相關(guān)系數(shù)R2為0.9988，SD僅為0.06 eV，說明G3B3方法可以精確預(yù)測不同類型物質(zhì)的電離能。部分文獻中未給出準(zhǔn)確電離能的煙氣化學(xué)成分，電離能的G3B3方法計算結(jié)果如表2所示。根據(jù)文獻中給定的精確電離能以及表2中數(shù)值，可以對采用同步輻射光電離質(zhì)譜技術(shù)在線分析卷煙煙氣時的一些卷煙煙氣組分進行定性分析。

2.2 電離能計算結(jié)果在卷煙煙氣組分定性分析中的應(yīng)用

圖3是卷煙煙氣氣相物在光子能量為10.50 eV和采樣10 s的條件下得到的光電離質(zhì)譜圖。從圖3可以看出，大多數(shù)強度較高的卷煙煙氣氣相組分的分子量都在100 Da范圍內(nèi)。很多物質(zhì)都具有同分異構(gòu)體，其數(shù)量隨分子量的增加而迅速增多，由于同分異構(gòu)體的質(zhì)量數(shù)相同，故在質(zhì)譜上不能準(zhǔn)確分辨出每一個質(zhì)譜峰對應(yīng)的產(chǎn)物。為此，可以通過掃描光子能量得到每個化學(xué)組分的光電離效率譜曲線，并結(jié)合文獻報道的電離能數(shù)值和G3B3方法計算得到的電離能數(shù)值來確定產(chǎn)物的準(zhǔn)確歸屬，包括分辨出同分異構(gòu)體。

圖2 電離能（IE）的文獻值和兩種計算方法理論值的線性擬合關(guān)系

表2 部分未知電離能的煙氣化學(xué)成分的G3B3計算值（eV）

圖3 卷煙煙氣氣相物在光子能量為10.50 eV和采樣10 s時的光電離質(zhì)譜圖

由于同步輻射真空紫外光源具有波長連續(xù)可調(diào)的優(yōu)越性，而且其電離方式是一個單光子過程，能夠以軟電離的方式對物質(zhì)進行電離，幾乎只產(chǎn)生待測物的母體離子，而不產(chǎn)生碎片，所以當(dāng)光子的能量小于某物質(zhì)的電離能時，沒有信號產(chǎn)生，表現(xiàn)在光電離效率譜曲線上離子峰強度為0，一旦光子的能量大于這種物質(zhì)的電離能，離子信號會迅速上升，其“拐點”處所對應(yīng)的光子能量即為該物質(zhì)的電離能，如果同一個質(zhì)量數(shù)中包括多個電離能不一樣的物質(zhì)，光電離效率譜圖上就會出現(xiàn)幾個拐點，因為每種物質(zhì)都對應(yīng)著特定的電離能，據(jù)此便可以對待測物質(zhì)進行準(zhǔn)確的定性分析［15］。

圖4是卷煙煙氣氣相物中質(zhì)量數(shù)分別為30，54，60和87的光電效率譜曲線，其掃描波長范圍為8.00～11.00 eV。從圖4（a）中可以看到質(zhì)量數(shù)為30 Da的物質(zhì)在位于9.26和10.87 eV處出現(xiàn)了兩個拐點（也稱為“出現(xiàn)勢”），這與文獻［14］中所報道的一氧化氮（電離能：9.2642±0.00002 eV）和甲醛（電離能：10.88±0.01 eV）的電離能是相對應(yīng)的。由于本實驗測定的是光子能量為10.50 eV時的光電離質(zhì)譜圖（見圖3），在此能量下甲醛不會被電離。因此，圖3中質(zhì)量數(shù)為30 Da所對應(yīng)的物質(zhì)是一氧化氮。圖4（b）顯示質(zhì)量數(shù)為54 Da的物質(zhì)有3個拐點，對應(yīng)光子的能量分別為9.08，9.25和10.18 eV，其中9.08和10.18 eV分別與1，3-丁二烯[電離能：（9.072±0.007）eV]和1-丁炔[電離能：（10.18±0.01）eV]的電離能相對應(yīng)；而9.25 eV在文獻［14］給出的準(zhǔn)確電離能中找不出與之相對應(yīng)的電離能數(shù)值，但是該數(shù)值與采用G3B3方法計算得到的1，2-丁二烯的電離能（9.25 eV）完全吻合。所以，圖3中質(zhì)量數(shù)為54 Da的物質(zhì)對應(yīng)著1，3-丁二烯，1，2-丁二烯和1-丁炔3種組分。從圖4（c）中能夠發(fā)現(xiàn)質(zhì)量數(shù)為60 Da的物質(zhì)在掃描光電子能量為9.05 eV處存在一個出現(xiàn)勢，然而這一數(shù)值與中國煙草科學(xué)數(shù)據(jù)庫——卷煙煙氣分析數(shù)據(jù)分庫［13］中給出的所有組分的電離能都不吻合，這說明此種煙氣組分在該數(shù)據(jù)庫中是不存在的；另外，在文獻［14］給出的質(zhì)量數(shù)為60 Da的準(zhǔn)確電離能中也找不到與9.05 eV相對應(yīng)的數(shù)據(jù)。為此，采用G3B3方法對質(zhì)量數(shù)為60 Da的未知的或沒有經(jīng)過準(zhǔn)確評價的電離能進行了計算，結(jié)果發(fā)現(xiàn)9.05 eV這一數(shù)值與環(huán)硫乙烷電離能的計算結(jié)果（9.06 eV）很好地吻合。故圖3中質(zhì)量數(shù)為60 Da的組分為環(huán)硫乙烷。另外，從圖4（d）中能夠看出，質(zhì)量數(shù)為87 Da的物質(zhì)在光電子能量為9.42 eV時出現(xiàn)了一個拐點，這與采用G3B3方法計算得到的丙酮酰胺的電離能（9.45 eV）相一致，所以圖3中質(zhì)量數(shù)為87 Da的組分為丙酮酰胺。采用本方法可以對其他質(zhì)量數(shù)所對應(yīng)的物質(zhì)進行鑒定，卷煙煙氣氣相物中一些質(zhì)量數(shù)所對應(yīng)組分的定性分析結(jié)果見表3。

3 結(jié)論

圖4 卷煙煙氣氣相物中質(zhì)量數(shù)為30，54，60和87組分的光電效率譜

表3 卷煙煙氣氣相物中一些組分的定性分析結(jié)果

B3LYP/6-311++G（d，p）方法在預(yù)測電離能上存在較大的系統(tǒng)誤差；G3B3方法能精確計算不同類型卷煙煙氣化學(xué)組分的電離能，其計算結(jié)果可用于卷煙煙氣的定性分析中，這為進一步采用軟電離質(zhì)譜法，尤其是同步輻射光電離質(zhì)譜技術(shù)在線和實時分析卷煙煙氣化學(xué)成奠定了基礎(chǔ)。

［1］McRae D D.The physical and chemical nature of tobacco smoke［J］.Rec Adv Tob Sci，1990，16：233-323.

［2］Cueto R，Pryor W.Cigarette smoke chemistry：conversion of nitricoxidetonitrogendioxideandreactionsofnitrogen oxides with other smoke components as studied by Fourier transforminfraredspectroscopy［J］.VibSpectrosc，1994，7（1）：97-111.

［3］Ceschini P，Lafaye A.Evolution of the gas-bapour phase and the total particulate matter of cigarette smoke in a single puff［J］.Beitr Tabakforsch Int，1976，8（6）：378-381.

［4］Baker R R，Proctor C J.A smoke odyssey［C］//CORESTA Smoke-Technology Meeting，Xi'an，China，2001.

［5］Adam T，Mitschke S，Baker R R.Investigation of tobacco pyrolysis gases and puff-by-puff resolved cigarette smoke by single photon ionisation（SPI）-time-of-flight mass spectrometry（TOFMS）［J］.Beitr Tabakforsch Int，2009，23（4）：203-226.

［6］Plunkett S，Parrish M E，Shafer K H.Time-resolved analysis of cigarette combustion gases using a dual infrared tunable diode laser system［J］.Vib Spectrosc，2001，27（1）：53-63.

［7］LiuC，F(xiàn)engS，HeemstJ，etal.Newinsightsintothe formation of volatile compounds in mainstream cigarette smoke［J］.Anal Bioanal Chem，2010，396（5）：1817-1830.

［8］Mitschke S，Adam T，Streibel T，et al.Application of time-offlightmassspectrometrywithlaser-basedphotoionization methodsfortime-resolvedon-lineanalysisofmainstream cigarette smoke［J］.Anal Chem，2005，77（8）：2288-2296.

［9］曾榮今，曹晨忠，袁華.分子拓?fù)鋵W(xué)方法估算多環(huán)芳香烴類化合物的電離能［J］.化學(xué)物理學(xué)報，2005，18（3）：372-378.

［10］姜賽紅，聶長明，林達，等.誘導(dǎo)效應(yīng)描述符與胺、醇、醚類化合物電離能的關(guān)系［J］.分子科學(xué)學(xué)報，2007，23（6）：411-415.

［11］Joantéguy S，Pfister-Guillouzo G，Chermette H.Assessment of density functional methods for the calculation of ionization potentials of unsaturated molecules［J］.J Phys Chem A，1999，103（18）：3505-3511.

［12］Namazian M，Coote M L.G3 calculations of the proton affinity and ionization energy of dimethyl methylphosphonate［J］.J Chem Thermodyn，2008，40（7）：1116-1119.

［13］中國煙草科技信息中心.中國煙草科學(xué)技術(shù)數(shù)據(jù)庫［DB/OL］//http：//sjk.ztri.com.cn/tobacco/sjkList.jsp

［14］Linstrom P J，Mallard W G.NIST Chemistry Webbook［C/OL］//National Institute of Standards and Technology，Gaithersburg，MD，2008.http：//webbook.nist.gov.

［15］齊飛.燃燒：一個不息的話題——同步輻射單光子電離技術(shù)在燃燒研究中的應(yīng)用［J］.物理，2006，35（1）：1-6.