施斌 袁荔 唐天宇 陸利敏 趙先豪 魏曉楠 唐延林

(貴州大學(xué)物理學(xué)院, 貴陽 550025)

特丁基對苯二酚是重要的食品抗氧化劑.理論上, 基于密度泛函理論, 采用B3LYP 泛函及6?311G(d,p)基組在氣相環(huán)境下優(yōu)化分子的結(jié)構(gòu)并進(jìn)行頻率計(jì)算.在此基礎(chǔ)上, 基于含時(shí)密度泛函理論, 選用SMD(solvation model based on density)溶劑模型, 利用B3LYP 泛函并結(jié)合def2?TZVP 基組計(jì)算分子在無水乙醇溶劑中的前50 個(gè)激發(fā)態(tài).再通過Multiwfn軟件對紅外光譜做振動分析并考察分子間相互作用對紅外光譜的影響, 對紫外光譜做分子軌道和電子空穴分析.實(shí)驗(yàn)上, 通過KBr 壓片法, 利用傅里葉紅外變換光譜儀測定樣品紅外光譜.采用液相法, 以乙醇為溶劑, 利用紫外可見分光光度計(jì)測定樣品紫外光譜.通過對比分析可知, 理論光譜與實(shí)驗(yàn)光譜總體吻合較好.紅外光譜各基團(tuán)的特征吸收峰都較為明顯且較好吻合, 特丁基對苯二酚二聚體存在氫鍵作用, 這使得O—H 鍵的強(qiáng)度被削弱, 導(dǎo)致吸收頻率降低并在3670—3070 cm—1 處出現(xiàn)一個(gè)寬峰.紫外光譜主要由基態(tài)躍遷至第1, 2, 6, 7 激發(fā)態(tài)形成, 最大吸收峰位于200 nm 以下, 為π→π*和σ→π*躍遷形成, 268.8 nm 和221.4 nm 處的吸收峰均為n→π*和π→π*躍遷形成.由電子空穴圖可知, 這4 個(gè)主要激發(fā)均為電子局域激發(fā).

1 引 言

食品抗氧化劑作為一種十分重要的食品添加劑,能夠有效阻止食品的氧化, 提高食品的安全性[1].特丁基對苯二酚(tert?butylhydroquinone, TBHQ)作為主要的人工合成抗氧化劑, 具有十分優(yōu)異的抗氧化能力, 被廣泛的運(yùn)用到動植物脂肪和富脂食品中[2?4].TBHQ 比食品抗氧化劑丁基羥基茴香醚(butylated hydroxyanisole, BHA)、二丁基羥基甲苯、沒食子酸丙酯的抗氧化能力強(qiáng)5—7 倍.然而,過量的TBHQ 可能具有致癌的風(fēng)險(xiǎn)[5].因此, 對食品抗氧化劑的檢測及結(jié)構(gòu)性質(zhì)的研究十分重要.

目前, 對TBHQ 的抗氧化性質(zhì)及機(jī)理的研究較多, 對光譜的研究較少[6].逯美紅等[7]基于密度泛函理論(density functional theory, DFT)和實(shí)驗(yàn)對TBHQ、BHA 的拉曼光譜進(jìn)行對比研究, 并對主要的拉曼特征峰的振動模式進(jìn)行了指認(rèn), 為TBHQ 的檢測提供了一定的依據(jù), 但該研究僅運(yùn)用拉曼光譜進(jìn)行分析且沒有對TBHQ 的分子間相互作用和激發(fā)性質(zhì)進(jìn)行研究.Rode 等[8]同樣基于DFT, 模擬手性對噻吩磺酰胺單體和二聚體的紅外、紫外光譜并與實(shí)驗(yàn)光譜進(jìn)行對比研究, 充分研究了其幾何結(jié)構(gòu)和理化性質(zhì).Snehalatha 等[9]對食品添加劑卡莫依辛進(jìn)行了光譜分析和DFT 計(jì)算, 通過對比理論和實(shí)驗(yàn)紅外光譜并進(jìn)行振動分析, 對了解分子結(jié)構(gòu)具有重要意義, 但是該研究并沒有詳細(xì)解釋氫鍵對紅外光譜可能造成的影響.以上基于量子化學(xué)方法, 通過對比實(shí)驗(yàn)和理論光譜, 對了解物質(zhì)的結(jié)構(gòu)、性質(zhì)及相關(guān)檢測十分重要.而目前關(guān)于TBHQ 的紅外、紫外光譜檢測和結(jié)構(gòu)性質(zhì)的研究卻很少.

因此, 本文對TBHQ 的理論和實(shí)驗(yàn)紅外、紫外光譜進(jìn)行比較分析, 一方面基于DFT 對紅外光譜進(jìn)行振動分析并考察分子間弱相互作用對紅外光譜的影響, 另一方面對紫外光譜的主要激發(fā)態(tài)做分子軌道和電子空穴分析, 可以充分了解光譜特征、分子結(jié)構(gòu)及激發(fā)性質(zhì), 對改造、檢測相關(guān)食品抗氧化劑具有重要意義.

2 方 法

2.1 計(jì)算方法

本文利用Gaussian16 軟件[10]基于密度泛函理論, 對TBHQ 分子的紅外、紫外光譜進(jìn)行計(jì)算.首先, 利用B3LYP/6?311G(d, p)級別[11?12]在氣相環(huán)境下優(yōu)化分子的結(jié)構(gòu)并進(jìn)行頻率計(jì)算.在幾何優(yōu)化的每一步都精確地計(jì)算Hessian 矩陣, 優(yōu)化結(jié)果均收斂且頻率計(jì)算無虛頻, 相關(guān)結(jié)構(gòu)優(yōu)化參數(shù)如表1 所示.再利用Multiwfn 3.7[13]軟件對頻率進(jìn)行校正并繪制紅外光譜.由于能量的計(jì)算對基組的敏感性遠(yuǎn)高于結(jié)構(gòu)優(yōu)化和振動分析, 因此, 在結(jié)構(gòu)優(yōu)化的基礎(chǔ)上, 基于含時(shí)密度泛函理論(time?depen?dent density functional theory, TDDFT), 在B3L?YP/def2?TZVP 級別下[14]計(jì)算分子在無水乙醇溶劑中前50 個(gè)激發(fā)態(tài)并繪制紫外光譜和電子空穴圖[15,16].為了充分考慮溶劑效應(yīng)對結(jié)果的影響, 體現(xiàn)溶質(zhì)與溶劑間的靜電相互作用和非靜電相互作用, 激發(fā)能的計(jì)算都是在SMD 隱式溶劑模型[17]下進(jìn)行.

為了研究TBHQ 分子間的弱相互作用對紅外光譜的影響, 本文利用genmer 軟件結(jié)合molclus 1.9.6 軟件[18]對該分子的二聚體構(gòu)型進(jìn)行搜索, 并得到15 個(gè)構(gòu)型.再調(diào)用Gaussian 16 軟件, 在B3LYP?D3(BJ)/6?311G(d, p)級別對15 個(gè)二聚體構(gòu)型進(jìn)行結(jié)構(gòu)優(yōu)化, 并采用M062X?D3 泛函在6?311+G(2d, p)基組水平上計(jì)算它們的單點(diǎn)能[19?23].對TBHQ 二聚體的弱相互作用的研究, 是利用約化密度梯度(reduced density gradient, RDG)函數(shù)方法[24].RDG 函數(shù)是用來描述均勻電子分布偏差的實(shí)空間函數(shù), 表達(dá)式為

其中ρ 表示電子密度.sign (λ2)函數(shù)是電子密度Hessian 矩陣的第二大本征值λ2的符號.而弱相互作用的特點(diǎn)為RDG 函數(shù)值和ρ 函數(shù)值都較小, 可以通過在實(shí)空間相應(yīng)區(qū)域形成梯度等值面進(jìn)行表示, 再將sign(λ2)ρ 函數(shù)投影到RDG 等值面上, 則可以表現(xiàn)出弱相互作用的強(qiáng)度和類型[23].

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

本文通過實(shí)驗(yàn), 采用KBr 壓片法對TBHQ 紅外光譜進(jìn)行測定.實(shí)驗(yàn)前, 將TBHQ 和KBr 粉末烘干并按1∶100 比例充分研磨, 再壓制成均勻透明薄片, 利用Bruker Daltonics 公司生產(chǎn)的傅里葉紅外變換光譜儀測定.紫外光譜的測定是采用液相法, 以乙醇為溶劑, 利用MAPADA 公司生產(chǎn)的UV?6100S 型紫外可見分光光度計(jì)測定.實(shí)驗(yàn)采用純度為99%的食品抗氧化劑TBHQ 均為浙江一諾生物科技有限公司生產(chǎn).

3 結(jié)果與分析

3.1 幾何結(jié)構(gòu)

食品抗氧化劑TBHQ 為白色粉末, 化學(xué)分子式為C10H14O2.利用B3LYP 泛函, 在6?311G(d, p)基組水平上將其幾何構(gòu)型優(yōu)化至局域能量極小值的穩(wěn)定結(jié)構(gòu), 如圖1 所示, 其主要鍵長、鍵角和二面角列于表2.

3.2 紅外光譜

在結(jié)構(gòu)優(yōu)化的相同級別下, 對TBHQ 分子的頻率進(jìn)行計(jì)算, 無虛頻.再利用Multiwfn 3.7 軟件對所有頻率進(jìn)行校正, 校正因子為0.967, 并繪制氣相下4000—400 cm—1的紅外光譜.實(shí)驗(yàn)上, 利用KBr 壓片法測定紅外光譜.其中, 計(jì)算光譜為圖中的紅色譜線, 對應(yīng)于左邊Y 軸的摩爾吸收系數(shù),實(shí)驗(yàn)光譜為圖中的藍(lán)色譜線, 對應(yīng)于右邊Y 軸的透過率.TBHQ 的紅外光譜如圖2 所示, 其特征吸收峰的振動分析如表3 所示.

表1 TBHQ 分子幾何結(jié)構(gòu)優(yōu)化參數(shù)Table 1.Optimized geometry structure parameters of TBHQ.

圖1 TBHQ 分子幾何構(gòu)型Fig.1.Geometric structure of TBHQ molecule.

表2 TBHQ 分子主要鍵長、鍵角、二面角Table 2.The main bond length, bond angle, and dihedral angle of TBHQ.

圖2 TBHQ 紅外光譜Fig.2.IR spectra of TBHQ.

通過對TBHQ 理論光譜和實(shí)驗(yàn)光譜的對比分析可知特征吸收峰整體吻合較好.位于3001.16 cm—1處苯環(huán)上的C—H 伸縮振動, 1589.30, 1485.15 cm—1處苯環(huán)的環(huán)伸縮振動和771.51 cm—1處苯環(huán)上C—H 面外彎曲振動的吸收峰為芳烴的特征吸收峰; 2962.58 cm—1處C—H 伸縮振動和1442.72 cm—1處C—H 彎曲振動為烷烴的特征吸收峰; 3288.54 cm—1處O—H 伸縮振動和1309.63 cm—1處C—O 伸縮振動為酚的特征吸收峰.2358.88 cm—1處的小吸收峰可能是由于樣品中含有碳氮三鍵的雜質(zhì)引起.而理論光譜和實(shí)驗(yàn)光譜在3000 cm—1以上的區(qū)域存在較大差異, 可能是由于分子間氫鍵的作用造成.

為了探究TBHQ 分子之間是否形成氫鍵[25],首先利用genmer 軟件結(jié)合molclus1.9.6 軟件對該分子的二聚體構(gòu)型進(jìn)行搜索, 并得到15 個(gè)構(gòu)型,再調(diào)用Gaussian16 軟件, 利用B3LYP?D3(BJ)泛函, 結(jié)合6?311G(d, p)基組對15 個(gè)二聚體構(gòu)型進(jìn)行結(jié)構(gòu)優(yōu)化, 并在M062X?D3/6?311+G(2d, p)下計(jì)算它們的單點(diǎn)能[23].15 個(gè)二聚體能量如圖3 所示.再通過Multiwfn 3.7 軟件結(jié)合VMD 軟件[26]繪制出TBHQ 二聚體的RDG 函數(shù)等值面圖[24,27].通過分析可知, 二聚體3, 8, 10, 13 和15 都存在分子間氫鍵的作用, 它們的RDG 函數(shù)等值面圖如圖4 所示.圖中, 等值面為藍(lán)色區(qū)域表示氫鍵、強(qiáng)鹵鍵等強(qiáng)吸引作用, 綠色區(qū)域表示范德瓦耳斯作用, 紅色區(qū)域表示在環(huán)、籠中出現(xiàn)較強(qiáng)的位阻效應(yīng)等強(qiáng)互斥作用[23].

表3 TBHQ 分子紅外光譜振動分析Table 3.Vibration analysis of IR spectra of TBHQ.

圖3 TBHQ 二聚體的單點(diǎn)能Fig.3.Conformational energy of the dimer.

二聚體8 和10 的能量較低, 為—1079.7784744和—1079.7784757 Hartree, 且二聚體8 和10的總均 方 根 位 移(root mean square displacement/deviation)為0.00004117 nm, 它們的結(jié)構(gòu)十分相似, 極可能為TBHQ 二聚體的穩(wěn)定結(jié)構(gòu), 因此, 做出二聚體10 的紅外光譜如圖5 所示.而二聚體3 和13 的結(jié)構(gòu)也較為相似.通過對圖4 進(jìn)行分析可知, TBHQ 二聚體RDG 函數(shù)等值面在兩個(gè)羥基之間形成藍(lán)色等值面, 即形成了分子間氫鍵.通過對比單分子的紅外光譜可知, 二聚體的紅外光譜在3563.85 cm—1處出現(xiàn)一個(gè)吸收峰, 由振動分析可知, 其主要是由于形成分子間氫鍵的兩個(gè)酚羥基的O—H 伸縮振動形成, 而3709.24 cm—1處的吸收峰為未形成分子間氫鍵的酚羥基的O—H 伸縮振動形成.所以, 分子間氫鍵的形成削弱了O—H 鍵的強(qiáng)度, 導(dǎo)致吸收頻率降低.不同的分子間氫鍵使O—H 鍵的削弱程度不同, 在疊加效果下使得O—H 伸縮振動在紅外光譜3670—3070 cm—1處出現(xiàn)一個(gè)寬峰.而理論光譜是模擬氣相環(huán)境下單分子的紅外光譜, 并不存在分子間氫鍵的作用, 導(dǎo)致理論光譜和實(shí)驗(yàn)光譜在3000 cm—1以上的區(qū)域存在較大差異.

圖4 TBHQ 二聚體RDG 函數(shù)等值面圖 (a) 二聚體3; (b) 二聚體8; (c) 二聚體10; (d) 二聚體13; (e) 二聚體15Fig.4.RDG function isosurface map of dimer: (a) Dimer 3; (b) dimer 8; (c) dimer 10; (d) dimer 13; (e) dimer 15.

圖5 二聚體10 的紅外光譜Fig.5.IR spectrum of dimer 10.

3.3 紫外光譜

基 于TDDFT, 在B3LYP/def2?TZVP 級 別下, 利用SMD 溶劑模型, 計(jì)算分子在無水乙醇溶劑中前50 個(gè)激發(fā)態(tài)并繪制紫外光譜, 再利用Mul?tiwfn 3.7 軟件分析主要激發(fā)態(tài)的分子軌道特征和電子?空穴圖[16].實(shí)驗(yàn)上, 利用液相法將樣品溶于乙醇溶液并測定紫外光譜.為了更好的將理論光譜和實(shí)驗(yàn)光譜進(jìn)行對比, 本文僅關(guān)注200—340 nm 波段振子強(qiáng)度大于0.05 的激發(fā)態(tài), TBHQ 分子的紫外光譜如圖6 和7 所示, 激發(fā)性質(zhì)如表4 所示.

通過對TBHQ 紫外光譜的對比分析可知, 理論和實(shí)驗(yàn)光譜整體吻合較好.理論光譜最大吸收峰位于200 nm 以下的真空紫外區(qū), 200 nm 以上有兩個(gè)小的吸收峰, 分別位于221.4 nm 和268.8 nm,主要由基態(tài)躍遷至1, 2, 6, 7 激發(fā)態(tài)形成.其中, 電子由第45 軌道(最高占據(jù)軌道)激發(fā)到第46 軌道(最低空軌道)對S0—S1 激發(fā)的貢獻(xiàn)較大; 電子由第45 軌道激發(fā)到第47 軌道對S0—S2 激發(fā)的貢獻(xiàn)較大; 電子由第44 軌道激發(fā)到第46 軌道對S0—S6激發(fā)的貢獻(xiàn)較大; 電子由第44 軌道激發(fā)到第47軌道對S0—S7 激發(fā)的貢獻(xiàn)較大.通過繪制以上4 個(gè)軌道的波函數(shù)等值面圖, 對激發(fā)特征進(jìn)行分析,如圖8所示.其中, 軌道波函數(shù)正相位用紅色等值面表示, 負(fù)相位用藍(lán)色等值面表示, 等值面大小未標(biāo)注時(shí)皆為0.05.通過分析可知, 第44 軌道的苯基上具有π軌道特征, 丁基上的C—C 鍵具有σ 軌道特征; 第45 軌道具有π 和n 軌道特征, 當(dāng)?shù)戎得娲笮?.18時(shí)n 軌道特征更明顯; 第46 軌道表現(xiàn)出π*軌道特征; 第47 軌道表現(xiàn)出n 和π*軌道特征, 當(dāng)?shù)戎得娲笮?.14 時(shí)π*軌道特征更明顯.因此, 268.8 nm處的吸收峰為n→π*和π→π*躍遷形成;221.4 nm處的吸收峰為n→π*和π→π*躍遷形成;200 nm 以下的吸收峰為π→π*和σ→π*躍遷形成.由實(shí)驗(yàn)光譜可知, 其最大吸收峰同樣位于200 nm 以下, 在229.1 nm 和293.7 nm 處有吸收峰, 而靠近200 nm附近出現(xiàn)了尖峰, 主要是由于該波段接近真空紫外區(qū), 受制于儀器自身測量范圍和精度的影響所造成的誤差.通過對上述光譜的分析可知, 實(shí)驗(yàn)光譜相對理論光譜出現(xiàn)輕微紅移, 一方面可能是由于乙醇溶劑的極性造成了π→π*躍遷譜帶向長波方向位移, 另一方面, 可能是由于理論光譜是模擬單分子在溶劑下的光譜情況, 而實(shí)驗(yàn)光譜是較多分子在溶劑下形成的光譜.

圖6 TBHQ 分子理論紫外光譜Fig.6.Theoretical UV spectrum of TBHQ.

圖7 TBHQ 分子實(shí)驗(yàn)紫外光譜Fig.7.Experimental UV spectrum of TBHQ.

表4 TBHQ 分子的激發(fā)性質(zhì)Table 4.Excited state properties of TBHQ.

圖8 分子軌道波函數(shù)等值面圖 (a) 第44 軌道; (b) 第45 軌道; (c) 第45 軌道(等值面為0.18); (d) 第46 軌道; (e) 第47 軌道;(f) 第47 軌道(等值面為0.14)Fig.8.Isosurface of molecular orbital wavefunction: (a) The 44th molecular orbital; (b) the 45th molecular orbital; (c) the 45th mo?lecular orbital (isosurface = 0.18); (d) the 46th molecular orbital; (e) the 47th molecular orbital; (f) the 47 th molecular orbital(isosurface = 0.14).

圖9 TBHQ 的電子空穴圖 (a) S0—S1; (b) S0—S2; (c) S0—S6; (d) S0—S7Fig.9.Electron?hole distribution diagram of TBHQ: (a) S0—S1; (b) S0—S2; (c) S0—S6; (d) S0—S7.

通過對TBHQ 紫外光譜的分析可知, 在200—340 nm 波段, 其主要由基態(tài)躍遷至第1, 2, 6, 7 激發(fā)態(tài)形成, 因此對這4 個(gè)激發(fā)做電子空穴分析.通過Multiwfn3.7 結(jié)合VMD 軟件, 分別做出它們的電子空穴圖, 如圖9 所示, 其等值面值均為0.002.其中, 綠色區(qū)域?yàn)殡娮拥戎得? 藍(lán)色區(qū)域?yàn)榭昭ǖ戎得?表5 為相關(guān)激發(fā)的電子空穴相關(guān)參數(shù).其中,Sr衡量電子和空穴的重疊程度, 最大值1 說明完全重疊, 最小值0 說明完全不重疊; D 指數(shù)衡量電子與空穴質(zhì)心之間的距離; H 指數(shù)是衡量電子與空穴的總體平均分布廣度; t 指數(shù)是衡量電子與空穴的分離程度, t ＞ 0 說明電子與空穴分離較為充分, t ＜ 0說明電子與空穴沒有充分分離; 空穴離域指數(shù)(hole delocalization index, HDI)和電子離域指數(shù)(electron delocalization index, EDI)分別衡量空穴和電子的離域情況, 它們的數(shù)值越小, 離域程度越高, 分布的均勻程度越大[23].通過結(jié)合圖形和相關(guān)參數(shù)分析可知, 這4 個(gè)激發(fā)的電子和空穴主要分布于苯環(huán)和酚羥基上, Sr指數(shù)都較大, 尤其是S0—S2 和S0—S6 激發(fā)的電子和空穴幾乎完全重疊;而D 指數(shù)較小, 說明電子與空穴質(zhì)心很接近; t 指數(shù)都為負(fù)值, 說明電子和空穴沒有充分分離; H 指數(shù)都較為接近, 說明電子和空穴的總體分布廣度都差不多, 基本都覆蓋于體系中的苯環(huán)和酚羥基上;這4 個(gè)激發(fā)的HDI 和EDI 值并不是很大且都較為接近, 說明空穴和電子的離域特征明顯, 離域于整個(gè)苯環(huán)上.通過以上分析可知, 這4 個(gè)激發(fā)為電子局域激發(fā).

表5 電子空穴相關(guān)參數(shù)Table 5.Electron?hole distribution related parameters of TBHQ.

4 結(jié) 論

本文通過理論和實(shí)驗(yàn)相結(jié)合的方式, 對比研究了食品抗氧化劑TBHQ 的紅外光譜和紫外光譜,基于DFT 和TDDFT 得到的理論光譜與實(shí)驗(yàn)光譜總體吻合較好.對紅外光譜的振動分析可知, 各基團(tuán)的特征吸收峰都較為明顯且較好吻合, 而實(shí)驗(yàn)光譜在3670—3070 cm—1處出現(xiàn)一個(gè)寬峰是由于分子間氫鍵的O—H 伸縮振動形成.紫外光譜最大吸收峰位于200 nm 以下, 為π→π*和σ→π*躍遷形成, 268.8 nm 和221.4 nm 處的吸收峰均為n→π*和π→π*躍遷形成.由電子?空穴分析可知, 基態(tài)躍遷至第1, 2, 6, 7 激發(fā)態(tài)均為電子局域激發(fā).

通過對TBHQ 光譜的對比研究, 可以充分了解紅外、紫外光譜的形成機(jī)理, 并充分了解TBHQ的分子結(jié)構(gòu)和激發(fā)性質(zhì).這對了解、改造以及檢測食品抗氧化劑TBHQ 具有重要意義.