滿梅玲 陸春海 陳文凱,* 李 奕 章永凡

(1福建省光催化重點實驗室—省部共建國家重點實驗室培育基地,福州大學化學系,福州350002; 2成都理工大學核技術與自動化工程學院,成都610059)

一氧化碳分子在Pt/t-ZrO2(101)表面的吸附性質

滿梅玲1陸春海2陳文凱1,*李 奕1章永凡1

(1福建省光催化重點實驗室—省部共建國家重點實驗室培育基地,福州大學化學系,福州350002;2成都理工大學核技術與自動化工程學院,成都610059)

運用廣義梯度密度泛函理論(GGA-PW91)結合周期平板模型方法,研究了CO分子在完整與Pt負載的四方ZrO2(101)表面的吸附行為.結果表明:表面第二層第二氧位和表面第二橋位分別為CO分子和Pt原子在完整ZrO2(101)表面的穩(wěn)定吸附位,且覆蓋度為0.25 ML(monolayer)時均為穩(wěn)定吸附構型,吸附能分別為56.2和352.7 kJ·mol-1.CO分子在負載表面的穩(wěn)定吸附模式為C-end吸附,吸附能為323.8 kJ·mol-1.考察了CO分子在負載表面吸附前后的振動頻率、態(tài)密度和軌道電荷布居分析,并與CO分子和Pt原子在ZrO2表面的結果進行比較.結果表明,C端吸附CO分子鍵長為0.1161 nm,與自由的和吸附在ZrO2表面后的CO相應值(0.1141和0.1136 nm)相比伸長.吸附后C―O鍵伸縮振動頻率為2018 cm-1,與自由CO分子相比發(fā)生紅移;吸附后CO帶部分正電荷,電子轉移以Pt 5d?CO 2π的π反饋機理占主導地位.

密度泛函理論;四方型二氧化鋯;一氧化碳;吸附;鉑

1 引言

一氧化碳在氫燃料電池中被稱為失活催化劑,1一氧化碳的氧化反應是研究氧化型催化劑活性最普遍的反應;2在聚合物電解質燃料電池中,即使很低濃度的CO也能夠在陽極產生CO中毒而引起很嚴重的破壞.3如何消除一氧化碳或將其轉化成無害物質已成為實驗4-9和理論10-22研究的熱點,而金屬氧化物負載貴金屬是非常好的消除CO的催化劑. CO在貴金屬表面上吸附的反應機制通常起關鍵作用的是活潑金屬和載體,23貴金屬如Pt、Pd、Rh等是三效催化劑中的有效成分;貴金屬在金屬氧化物上的負載能增強催化活性和減少催化劑成本,同時能夠將CO、HCs(hydrocarbons)和NOx等污染氣體轉化為CO2和N2.在眾多金屬氧化物中,ZrO2由于具有特定的物理化學性質,包括表面酸堿性、氧化還原性、良好的熱穩(wěn)定性、抗催化中毒性,使其成為一種很具吸引力的催化劑和催化劑載體.24-27

江凌等28報道了CO在過渡金屬表面上的吸附,結果表明吸附強度為Pd＞Pt＞Cu＞Ag＞Au.任云鵬等29用密度泛函理論(DFT)方法研究了CO分子在Pt低指數(shù)面的吸附,結果發(fā)現(xiàn),當CO的覆蓋率為0.25 ML時,其在Pt(100)面橋位的吸附能為203.6 kJ· mol-1;在Pt(100)面最大吸附能為228.7 kJ·mol-1;在Pt(111)面hcp位的吸附能為189.1 kJ·mol-1,與文獻30計算的201.7 kJ·mol-1相吻合.Yamagishi等31用DFT方法研究了CO分子在Pt(100)、Pt(410)和Pt(110)面的吸附,研究考慮了吸附的吸附能、吸附構型和振動頻率,結果表明,覆蓋度對CO分子在不同Pt表面上的吸附能影響并不明顯,且穩(wěn)定吸附能和頻率的理論預測值與實驗數(shù)據(jù)32,33相吻合.Anez等34用Gaussian 03軟件包ONIOM2方法35研究了CO和H2O分子在四方ZrO2(110)和(101)表面的吸附,結果表明CO分子在(110)表面的最佳吸附位為連接兩個Zr原子的橋位,吸附能為90.8 kJ·mol-1;在(101)表面的最佳吸附位為連接Zr和O原子的橋位,吸附能為97.9 kJ·mol-1.Dilara等36研究發(fā)現(xiàn)在300 K時Pt可以在ZrO2(100)面形成層狀結構,但當溫度超過700 K時,Pt層就會凝聚成顆粒狀,高分辨電子能量損失譜證實CO能在Pt-ZrO2(100)界面發(fā)生吸附. Zhu等37研究了CO在Pt/ZrO2上的吸附,結果表明,隨著燒結溫度的增加,CO的吸附量明顯減少.Graf等38報道了Pt/ZrO2催化水煤氣轉換反應中羥基的反應性,研究表明單配位的羥基在ZrO2載體上能與CO反應形成甲酸鹽,多配位的羥基甲酸鹽需在Pt存在時才會發(fā)生分解生成H2和CO2,說明甲酸鹽物種的反應性與Pt顆粒有關.同時,Pt/ZrO2已經(jīng)在許多領域有廣泛應用,如作為傳感器、固體氧化物燃料電池39和催化劑.40,41目前對于CO在Pt/ZrO2的吸附研究理論研究方面未見報道.

本文采用密度泛函理論結合周期平板模型方法,模擬了Pt負載前后CO分子的吸附行為.對不同吸附構型進行構型優(yōu)化以及電荷、頻率和態(tài)密度的分析,目的主要是探討CO分子吸附前后吸附性質的變化,為實驗研究提供一定的理論依據(jù).

2 模型和計算方法

四方ZrO2晶型相當于立方螢石型結構沿著c軸伸長而變形的晶體結構,空間群為P42/nmc.由于四方ZrO2最穩(wěn)定表面為(101)表面,42因此本文計算選取t-ZrO2(101)-2×2九層平板厚度來討論CO分子在Pt/ZrO2表面的吸附,如圖1所示.相鄰兩層平板間的真空層厚度設為1 nm,以確保平板間的分子間相互作用力足夠小.

通過對蚊蟲的標記、釋放、回捕研究發(fā)現(xiàn)，由于蚊蟲不具有對宿主和孽生地位置的記憶能力，所以蚊蟲在初期是通過隨機擴散搜尋吸血目標；一旦蚊蟲飛行到足以感應目標所在位置的距離之后，蚊蟲通過定向飛行接近宿主并完成吸血活動［10］。本研究以蚊籠中飼養(yǎng)的斯氏按蚊為研究對象，旨在通過對比不同條件下的按蚊吸血率，討論影響按蚊吸血行為中定向飛行階段的諸多因素。

本文采用廣義梯度密度泛函理論和周期平板模型方法模擬CO分子的吸附行為,所有計算均由Dmol3軟件包43-45實現(xiàn).廣義梯度近似下的PW91泛函46被用來描述交換相關能,Pt和Zr原子內層電子由有效核電勢(ECP)代替,價電子波函數(shù)采用雙數(shù)值基加極化函數(shù)(DNP)展開,非金屬原子C和O采用全電子基組.Brillouin zone積分的Monkhorst-Pack網(wǎng)絡參數(shù)設為2×3×1.優(yōu)化時固定下面六層,表面三層放開.

（1）紅柳泉油田下干柴溝組下段圈閉主要為巖性圈閉，類型有砂巖上傾尖滅巖性圈閉、物性封閉巖性圈閉（成巖圈閉）、砂巖透鏡體三種。

對于Pt原子在ZrO2(101)表面的吸附同樣考慮如圖1所示的八種吸附位.構型優(yōu)化后,除在Osub2吸附位外,其他吸附位構型均轉化為同一種構型,吸附能在331.8-352.7 kJ·mol-1之間.表3只列出了Osub2位和bri2位的構型參數(shù)、吸附能及電荷.由表可見,Pt原子帶部分負電荷,說明電子由底物向Pt原子轉移,且bri2位所帶電荷更多.吸附后Zr原子均得電子,O原子均失電子,說明O原子轉移電子給Pt原子,Pt原子又將部分電子轉移給Zr原子.表4給出了吸附前后Pt原子的軌道布居,固體Pt的價電子組成為6s15d9.吸附后在Osub2和bri2位分別變成6s1.0545d8.8226p0.276和6s1.1705d8.8846p0.142.可見,吸附后兩種吸附位的Pt原子d軌道分別失去0.178e和0.116e,而s軌道分別得到0.054e和0.170e,p軌道分別得到0.276e和0.142e,凈結果為Pt原子分別得到0.152e和0.196e.

圖1 t-ZrO2(101)-2×2超晶胞平板模型Fig.1 Side view and top view for the slab model of t-ZrO2(101)-2×2 supercellLight blue and red balls represent Zr and O atoms,respectively. Osurf1and Osurf2represent two different O sites of top surface;Osub1and Osub2represent two different O sites of sub-surface;Zraand Zrb represent two different Zr sites;bri1 and bri2 represent two different bridge sites.

表1 CO分子在ZrO2(101)表面吸附的幾何構型參數(shù)和吸附能Table 1 Geometrical parameters and adsorption energies for CO adsorbed on ZrO2(101)surface

3.3.1 吸附構型、吸附能和振動頻率

3 結果與討論

3.1 CO在ZrO2(101)面的吸附研究

3.1.1 不同吸附位的吸附活性

對CO分子在t-ZrO2(101)表面的吸附,本文選取如圖1所示的Osurf1、Osurf2(表層兩種不同O位),Osub1、Osub2(次表層兩種不同的O位),Zra、Zrb(a,b指代兩種不同Zr位),bri1和bri2(1,2指代兩種不同橋位)共八種吸附位,分別考慮了CO分子以C-end、O-end和parallel(以下簡稱para)吸附模式吸附于表面,不同構型優(yōu)化后,計算得到的幾何構型參數(shù)及吸附能列于表1.

基于以上分析,以下僅考慮bri2位上的吸附情況.對Pt原子的覆蓋度我們采用和CO分子一樣的方法,即定義0.25、0.5和1 ML分別指在ZrO2(101)-2×2超晶胞、1×2超晶胞和1×1超晶胞上分別吸附一個Pt原子時的覆蓋度.計算結果列于表5.如圖1俯視圖中的Zrb位有4個等價吸附位,Zrb位和Zr4之間的原子距離為0.3531 nm,而Zrb位和Zr6之間的原子距離為0.6060 nm.由于前者距離較小會使得Pt原子易與鄰近單胞的Pt原子發(fā)生相互作用而形成二聚體結構.

具體選擇哪一種接班模式，家族企業(yè)需要充分考慮內、外部因素及其對能力和忠誠兩個維度的影響。企業(yè)外部因素短期是無法去改變的，企業(yè)只能去適應。決定家族企業(yè)接班模式選擇的主要是內部因素。[7]企業(yè)在不同的發(fā)展階段，企業(yè)內部影響因素有很大的變化，從而出現(xiàn)不同階段選擇不同的接班模式，[8]見表2。

圖2給出了Pt原子在ZrO2(101)表面吸附前后的態(tài)密度圖,以利于進一步分析Pt原子的作用機理.從圖中可以看到,與吸附前相比,吸附后兩種吸附位的s、p軌道均向低能級方向移動,位于Fermi能級之下,說明兩種軌道均得電子.d軌道發(fā)生離域化,說明Pt原子的d軌道與底物存在較強的相互作用.與前面的分析相吻合,綜合說明bri2位吸附更穩(wěn)定.同時吸附能352.7 kJ·mol-1與Eichler57的計算值(363.8 kJ·mol-1)比較接近.

3.1.2 不同覆蓋度下CO分子的吸附活性比較

上述分析知道,CO分子以C-end在Osub2位吸附時最穩(wěn)定,因此以下僅討論C-end在Osub2位上的吸附.不同覆蓋下吸附活性不同,Yang等49在c-ZrO2(110)表面吸附CO分子時考慮了不同的覆蓋度的影響,計算結果表明,并不是覆蓋度越大吸附活性越大,而是在覆蓋度為0.25 ML時更穩(wěn)定.因此本文也采用同樣方法討論覆蓋度的影響.覆蓋度通過計算吸附CO分子數(shù)與表面Zr原子數(shù)的比值得到,即定義0.25、0.5和1 ML分別指在ZrO2(101)-2×2超晶胞、1×2超晶胞和1×1超晶胞上分別吸附一個CO分子時的覆蓋度.計算結果列于表2.需要指出的是,如圖1俯視圖中的Zra位有4個等價吸附位,Zra位和Zr1之間的原子距離為0.3633 nm,而Zra位和Zr2之間的原子距離為0.6402 nm.由于前者距離較小使得CO分子會與鄰近單胞的CO發(fā)生相互作用,從而導致吸附能降低.50-53

從表5可以看到,隨著覆蓋度的增加,吸附能相差不大,并且是逐漸減小的,鍵長相差也不大,因此三種覆蓋度下Pt原子在覆蓋度為0.25 ML時吸附構型更穩(wěn)定,也說明Pt趨向于分散而不是聚集.以下對于CO分子的吸附采用0.25 ML覆蓋度來討論.

3.2 Pt原子在ZrO2(101)表面的吸附

3.2.1 幾何構型參數(shù)和布居分析

吸附能定義為吸附前后各物質總能量的變化,吸附能為正值表示較強的吸附作用.

表2 CO分子在ZrO2(101)表面吸附的幾何構型參數(shù)和吸附能Table 2 Geometrical parameters and adsorption energies for CO adsorbed on ZrO2(101)surface

表3 Pt及與Pt成鍵的Zr、O原子在ZrO2(101)表面的幾何構型參數(shù)、吸附能和Mulliken電荷Table 3 Geometrical parameters,adsorption energies,and Mulliken population for Pt and Zr, O bonded with Pt on ZrO2(101)surface

表4 Pt原子吸附前后的軌道電子布居Table 4 Orbital population for atoms of Pt before and after adsorption

3.2.2 態(tài)密度分析

由表1可知,O-end吸附均為物理吸附,CO分子C-end吸附在Osub2、Zra、Zrb、bri1和bri2位及para吸附在Zra、Zrb、bri1和bri2位為化學吸附.吸附能越大,說明吸附質與襯底Zr原子的吸附鍵越強.CO以C-end吸附在Osub2位最穩(wěn)定,吸附能為56.2 kJ· mol-1.計算得到自由CO分子的鍵長為0.1141 nm,與實驗值(0.1128 nm)47相接近,吸附后C―O鍵與自由的CO分子相比,均有不同程度的縮短.吸附后C―O鍵的伸縮振動頻率為2172.8 cm-1,與文獻48中利用紅外光譜得到的CO分子吸附在四方ZrO2晶體上的振動頻率(2183 cm-1)很接近.

3.2.3 不同覆蓋度下Pt原子吸附活性的比較

圖2 Pt原子吸附前及吸附后在Osub2和bri2位的分態(tài)密度圖Fig.2 Partial density of states of Pt atom before and after adsorption on Osub2and bri2 sites

(5)白堊系。測區(qū)白堊系零星分布于西南部的柳州市大、小鵝山、涼水鐵路橋下、進德、魯比、咸甫、良泗、柳城鳳山地區(qū)以及洛埠附近等地，面積約10.96 km2。根據(jù)孢粉組合、巖石組合特征和與鄰區(qū)相對比，劃分下統(tǒng)和上統(tǒng)。測區(qū)內出露不連續(xù)，保存不完整，與下伏地層呈角度不整合接觸。下統(tǒng)稱永福(K1Y)群，上統(tǒng)稱西垌(K2x)組。

由表2可以看到,隨著覆蓋度的增加吸附能逐漸減小,54-56C―Zr和C―O鍵伸長.三種覆蓋度中覆蓋度為0.25 ML時CO分子的吸附構型為穩(wěn)定構型,這與我們選擇ZrO2(101)-2×2超晶胞是一致的.

3.3 CO分子在Pt/ZrO2(101)表面的吸附

其中Eadsorbate和Esubstrate分別表示吸附前吸附質和底物的能量,Eadsorptionsystem表示吸附后體系的總能量.

以上計算已經(jīng)得到Pt原子在t-ZrO2(101)-2×2表面的最有利吸附位為bri2位,因此我們選取此表面研究CO分子的吸附.吸附分別考慮C端、O端單齒配位和CO分子雙齒配位吸附共3種情況.構型優(yōu)化后,C端和雙齒配位吸附構型轉化為同一種構型. C-end和O-end兩種吸附構型的吸附能分別為323.8和77.9 kJ·mol-1.Anez等34用Gaussian 03軟件包ONIOM2方法35研究CO分子在四方ZrO2(101)表面的吸附,得到CO分子在連接Zr與O的橋位上的吸附能為97.9 kJ·mol-1,與我們的計算值(323.8 kJ· mol-1)相比較,說明Pt的負載更有利于C-end的吸附;Dulaurent等58利用紅外光譜研究了CO分子在Pt/ZrO2面的Pt位上的吸附,結果表明,當覆蓋度θ=0時的吸附能E0=195 kJ·mol-1和θ=1時的吸附能E1= 96 kJ·mol-1相比差別較大,這主要是由于文獻中Pt的分散度較低(5%)造成的.一般而言,CO分子都是采用C原子與金屬相連方式,O-end吸附模式不穩(wěn)定,因此本文只對C-end吸附進行研究.為清晰起見,構型圖僅選取表面三層.圖3為最穩(wěn)定吸附構型的側視圖.從圖3可以看出,構型優(yōu)化后CO分子發(fā)生傾斜,但CO分子與Pt原子之間的鍵角仍為180°.

表5 Pt原子在ZrO2(101)表面吸附的幾何構型參數(shù)和吸附能Table 5 Geometrical parameters and adsorption energies for Pt adsorbed on ZrO2(101)surface

由圖3可知,吸附后CO分子的鍵長為0.1161 nm,與自由CO分子相比鍵伸長,C―Pt鍵長為0.1838 nm,與文獻59報道的CO在Pt(111)面的吸附中,C―O鍵長為0.1140 nm,C―Pt鍵長為0.1850 nm相吻合.C端吸附時C―O鍵的振動頻率為2018 cm-1,Dulaurent等58利用紅外光譜研究CO吸附在Pt/ZrO2面Pt位上的頻率為2068 cm-1,我們的計算結果與其接近.吸附后C―O鍵的伸縮振動頻率與自由分子的伸縮振動頻率(2112 cm-1)相比發(fā)生紅移,與自由分子的實驗值(2138 cm-1)47相差不大.

3.3.2 Mulliken電荷布居和軌道布居分析

算法綜合比較了NLMS和FDNLMS算法的性能,比較指標包括處理時間以及自適應濾波器失調系數(shù)和回聲返回損耗增益值:

CO分子的電子組態(tài)為[(1σ)2(2σ)2(3σ)2(4σ)2(1π)4(5σ)2(2π)0].60表6列出了Pt原子吸附前后CO分子、Pt原子及與Pt、C成鍵的Zr、O(為了與CO分子中的O原子相區(qū)別,以O?表示)原子的Mulliken電荷布居.可見,與自由狀態(tài)下相比,吸附后CO分子均失去電子,Pt原子得到電子,底物O原子和Zr原子均得電子.這也可以由表7的C、O、Pt及與Pt成鍵的Zr、O原子吸附前后的軌道電子布居來解釋:負載Pt后,C原子的s軌道失去0.614e,p軌道得到0.198e;O原子的s、p軌道分別得到0.008e和0.052e;Pt原子的s、p軌道分別得到0.172e和0.198e,d軌道失去0.174e.底物Zr原子的s、p軌道分別失去0.002e和0.030e,d軌道得到0.110e;O原子s軌道失去0.004e, p軌道得到0.108e.吸附后CO分子帶部分正電荷,與表6結果相吻合.計算結果表明Pt原子負載后更有利于CO分子在ZrO2表面的吸附.

圖3 最穩(wěn)定吸附構型的側視圖Fig.3 Side view of the most stable adsorbed structuresGray balls represent C atoms.bond length in nm

圖4 吸附前后CO分子,表面Pt、Zr、O原子的態(tài)密度圖Fig.4 Density of states for CO,Pt,Zr,O atoms before and after CO adsorption

3.3.3 態(tài)密度分析

依照填埋場的服務范圍，對2017年及以后理論垃圾填埋量進行預測。隨著城市人口的穩(wěn)定和收運點數(shù)量的穩(wěn)定，城鄉(xiāng)環(huán)衛(wèi)一體化建設的基本實現(xiàn)，依據(jù)經(jīng)驗，生活垃圾產生量增長率按1%計。依據(jù)填埋場的運行預估，預測結果見表2。

對于自由CO分子,它的前線軌道成分主要來自于C原子的2s和2p、O原子的2s和2p軌道的貢獻,CO的5σ軌道成分來自于C原子的2s和2p的貢獻,而2π軌道成分則來自于C原子的2p和O原子的2p的貢獻.61由圖4可見,自由CO分子的最高占據(jù)軌道(HOMO)是位于Fermi能級處的5σ軌道,最低空軌道(LUMO)是2π軌道,且5σ和1π軌道之間存在一定的能級差,而一般認為1π軌道是定域化的,不參與成鍵.吸附后CO分子的4σ、5σ和2π軌道均向低能級方向移動,5σ軌道遠離Fermi能級,而且可以清楚看到,C-end吸附時,CO的5σ軌道與1π軌道很接近,CO的5σ軌道與Pt的6s軌道發(fā)生重疊,2π軌道與Pt的5d軌道發(fā)生重疊.這說明,C-end吸附時同時存在兩種轉移機理,即CO 5σ?Pt 6s的σ成鍵機理和Pt 5d?CO 2π的π反饋機理.從以上討論可以看出,吸附的CO分子帶部分正電荷,也說明CO 5σ?Pt 6s的σ成鍵機理占主要地位,但是CO的5σ軌道為非鍵軌道,對于鍵長和頻率的影響不大.因此Pt 5d?CO 2π的π反饋機理就占主導地位,導致C―O鍵伸長,相應的頻率發(fā)生紅移.與Pt原子成鍵的Zr、O原子吸附后軌道稍微向低能級發(fā)生偏移,這些結果也與表6的電荷布居分析結果相符.

三是要保證公正性。名師評選必須在推薦、審核、評審和審批四個環(huán)節(jié)嚴格把關，杜絕暗箱操作和名師“不名”的情況發(fā)生，確保評選結果的公正性、權威性和示范性，確保評選出的教學名師質量。

一定時期內傳統(tǒng)產業(yè)對一國經(jīng)濟發(fā)展起著支撐作用，發(fā)達國家也不例外.下面以日本、瑞典為例分析傳統(tǒng)產業(yè)發(fā)展的國際經(jīng)驗.

表6 CO分子及各原子和Pt以及與Pt成鍵的Zr、O(O?)原子的平均Mulliken電荷布居Table 6 Mulliken population for CO,atoms of CO,Pt and Zr,O bonded with Pt

表7 C、O、Pt及與Pt成鍵的Zr、O原子吸附前后的軌道電子布居Table 7 Orbital population for atoms of C,O,Pt and Zr, O bonded with Pt atom before and after adsorption

4 結論

采用密度泛函理論GGA-PW91方法模擬了CO分子在Pt/四方ZrO2(101)面最穩(wěn)定位的吸附行為.對CO分子和Pt原子在完整ZrO2(101)及CO分子在負載表面的吸附分別進行了構型優(yōu)化,并分析了它們的電子結構性質.研究表明:CO分子和Pt原子的穩(wěn)定吸附位分別為Osub2和bri2位,吸附能分別為56.2和352.7 kJ·mol-1.而CO分子在負載表面吸附時,C端吸附最穩(wěn)定,吸附能為323.8 kJ·mol-1,且伸縮振動頻率分析結果表明,C端吸附的頻率發(fā)生紅移,鍵伸長;態(tài)密度和布居電荷分析結果表明:CO吸附過程中帶部分正電荷,電子轉移以Pt 5d?CO 2π的π反饋機理占主導地位.通過對CO分子在負載型穩(wěn)定ZrO2面上的吸附研究為實驗中消除CO或將其轉化為無害物質提供理論依據(jù).

(1) Zou,H.B.;Chen,S.Z.;Liu,Z.L.;Lin,W.M.Int.J.Hydrog. Energy 2009,34,9324.

(2) Gong,X.Q.;Liu,Z.P.;Raval,R.;Hu,P.J.Am.Chem.Soc. 2004,126,8.

(3) Igarashi,H.;Fujino,T.;Watanabe,M.J.Electroanal.Chem. 1995,391,119.

(4) Dekkers,M.A.P.;Lippits,M.J.;Nieuwenhuys,B.E.Catal. Lett.1998,56(4),195.

(5) Rupprechter,G.;Dellwig,T.;Unterhalt,H.J.Phys.Chem.B 2001,105(18),3797.

(6)Hadjiivanov,K.;Kn?zinger,H.;Mihaylov,M.J.Phys.Chem.B 2002,106(10),2618.

(7) Shaikhutdinov,S.K.;Meyer,R.;Naschitzki,M.Catal.Lett. 2003,86(4),211.

(8)Chen,J.G.;Xiang,H.W.;Dong,Q.N.;Wang,X.Z.;Sun,Y.H. Acta Phys.-Chim.Sin.2001,17(2),161. [陳建剛,相宏偉,董慶年,王秀芝,孫予罕.物理化學學報,2001,17(2),161.]

(9)Chen,Y.Q.;Gong,M.C.;Zhu,X.H.;Ming,H.;Tang,S.H. Chem.J.Chin.Univ.2002,23(4),715.[陳耀強,龔茂初,祝小紅,明 虹,唐水花.高等學?；瘜W學報,2002,23(4),715.]

(10) Jia,Y.X.;Guo,X.Y.Acta Phys.-Chim.Sin.2005,21(3),306. [賈玉香,郭向云.物理化學學報,2005,21(3),306.]

(11) Li,P.Chem.Res.Appl.2005,17(6),741.[李 平.化學研究與應用,2005,17(6),741.]

(12)Nygren,M.A.;Pettersson,L.G.M.J.Chem.Phys.1996,105 (20),9339.

(13) Koper,M.T.;Shubina,T.E.;Van Santen,R.A.J.Phys.Chem. B 2002,106(3),686.

(14) Bredow,T.J.Phys.Chem.B 2002,106(28),7053.

(15) Orita,H.;Inada,Y.J.Phys.Chem.B 2005,109(47),22469.

(16)Zhou,J.;Li,Z.H.;Wang,W.N.;Fan,K.N.J.Phys.Chem.A 2006,110(22),7167.

(17) Dholabhai,P.P.;Ray,A.K.Phys.Scr.2007,75,506.

(18) Seebauer,E.G.;Kong,A.C.F.;Schmidt,L.D.J.Vac.Sci. Technol.A 1987,5(4-1),464.

(19)Wang,Y.Z.;Wu,R.F.;Zhao.Y.X.Catal.Today 2010,158 (3-4),470.

(20)Wasileski,S.A.;Koper,M.T.M.;Weaver,M.J.J.Phys.Chem. B 2001,105(17),3518.

(21) Doll,K.Surf.Sci.2004,573(3),464.

(22) Olsen,R.A.;Philipsen,P.H.T.;Baerends,E.J.J.Chem.Phys. 2003,119(8),4522.

(23) Panagiotopoulou,P.;Christodoulakis,A.;Kondarides,D.I.; Boghosian,S.J.Catal.2006,240(2),114.

(24) Srinivasan,R.;Davis,B.H.Catal.Lett.1992,14(2),165.

(25)Dang,Z.;Anderson,B.G.;Amenomiya,Y.;Morrow,B.A. J.Phys.Chem.1995,99(39),14437.

(26)Tanabe,K.Mater.Chem.Phys.1985,13(3-4),347.

(27)Tanabe,K.;Misono,M.;Ono,Y.;Hattori,H.Stud.Surf.Sci. Catal.1989,51,1.

(28) Jiang,L.;Wang,G.C.;Guan,N.J.;Wu,Y.;Cai,Z.S.;Pan,Y. M.;Zhao,X.Z.;Huang,W.;Li,Y.W.;Sun,Y.H.;Zhong,B. Acta Phys.-Chim.Sin.2003,19(5),393.[江 凌,王貴昌,關乃佳,吳 楊,蔡遵生,潘蔭明,趙學莊,黃 偉,李永旺,孫予罕,鐘 炳.物理化學學報,2003,19(5),393.]

(29) Ren,Y.P.;Lu,Y.X.;Lou,Q.Acta Phys.-Chim.Sin.2007,23 (11),1728.[任云鵬,魯玉祥,婁 琦.物理化學學報,2007, 23(11),1728.]

(30) Lyncha,M.;Hu,P.Surf.Sci.2000,458(1-3),1.

(31) Yamagishi,S.;Fujimoto,T.;Inada,Y.;Orita,H.J.Phys.Chem. B 2005,109,8899.

(32) Martin,R.;Gardner,P.;Bradshaw,A.M.Surf.Sci.1995,342, 69.

(33) Banholzer,W.F.;Parise,R.E.;Masel,R.I.Surf.Sci.1985,155 (2-3),653.

(34) Anez,R.;Sierraalta,A.;Martorell,G.J.Mol.Struct.2009,900, 59.

(35)Maseras,F.;Morokuma,K.J.Comp.Chem.1995,16(9),1170.

(36) Dilara,P.A.;Vohs,J.M.J.Phys.Chem.1995,99(47),17259.

(37) Zhu,X.;Xie,Y.;Liu,C.J.;Zhang,Y.P.J.Mol.Catal.A 2008, 282(1-2),67.

(38) Graf,P.O.;de Vlieger,D.J.M.;Mojet,B.L.;Lefferts,L. J.Catal.2009,262(2),181.

(39) Mitterdorfer,A.;Gauckler,L.J.Solid State Ionics 1999,117 (3-4),187.

(40) Bitter,J.H.;Seshan,K.;Lercher,J.A.J.Catal.1997,171(1), 279.

(41) Bitter,J.H.;Seshan,K.;Lercher,J.A.Top.Catal.2000,10 (3-4),295.

(42) Hofmann,A.;Clark,S.J.;Oppel,M.;Hahndorfk,I.Phys. Chem.Chem.Phys.2002,4,3500.

(43) Delly,B.J.Chem.Phys.1990,92(1),508.

(44) Delly,B.J.Chem.Phys.2000,113(18),7756.

(45) Delly,B.J.Phys.Chem.1996,100,6107.

(46) Perdew,J.P.;Chevary,J.A.;Vosko,S.H.;Jackson,K.A.; Pederson,M.R.;Singh,D.J.;Fiolhais,C.Phys.Rev.B 1992, 46(11),6671.

(47)Lide,D.R.CRC Handbook of Chemistry and Physics;CRC Press:Boca Raton,2003;pp 9-20.

(48)Ma,Z.Y.;Yang,C.;Wei,W.;Li,W.H.;Sun,Y.H.J.Mol. Catal.A:Chem.2005,227,119.

(49)Yang,Y.L.;Chen,W.K.;Guo,X.;Li,Y.;Zhang,Y.F.Chin.J. Struct.Chem.2010,29(7),1021.

(50)Nibbelke,R.H.;Campman,M.A.J.;Hoebink,J.H.B.J.; Marin,G.B.J.Catal.1997,171(2),358.

(51) Herz,R.K.;Marin,S.P.J.Catal.1980,65,281.

(52) Oh,H.;Fisher,G.B.;Carpenter,J.E.;Goodman,D.W. J.Catal.1986,100,360.

(53)Yeo,Y.Y.;Vattuone,L.;King,D.A.J.Chem.Phys.1997,106, 392.

(54) Chafk,T.;Dulaurent,O.;Gass,J.L.;Bianchi,D.J.Catal.1998, 179,503.

(55) Dulaurent,O.;Chandes,K.;Bouly,C.;Bianchi,D.J.Catal. 1999,188,237.

(56) Dulaurent,O.;Bianchi,D.Appl.Catal.A 2000,196,271.

(57) Eichler,A.Phys.Rev.B 2003,68,205408.

(58) Dulaurent,O.;Bianchi,D.Appl.Catal.A 2001,207(1-2),211.

(59) Bleakley,K.;Hu,P.J.Am.Chem.Soc.1999,121,7644.

(60) Pan,D.A.;Zhao,C.D.;Zheng,Z.X.Structures of Matters,2nd ed.;Higher Education Press:Beijing,1989;p 200 [潘道皚,趙成大,鄭載興.物質結構.第2版.北京:高等教育出版社, 1989:200]

(61) Sun,B.Z.;Chen,W.K.;Zheng,J.D.;Lu,C.H.Appl.Surf.Sci. 2008,255(5),3141.

August 25,2011;Revised:September 16,2011;Published on Web:September 30,2011.

Adsorption Properties of CO Molecules on Pt/t-ZrO2(101)Surfaces

MAN Mei-Ling1LU Chun-Hai2CHEN Wen-Kai1,*LI Yi1ZHANG Yong-Fan1
(1State Key Laboratory Breeding Base of Photocatalysis,Department of Chemistry,Fuzhou University,Fuzhou 350002,P.R.China;2College of Nuclear Technology and Automation Engineering,Chengdu University of Technology,Chengdu 610059,P.R.China)

We studied the adsorption of CO molecules on perfect and Pt-adsorbed t-ZrO2(101)surfaces using a periodic slab model by PW91 of generalized gradient approximation(GGA)within the framework of density functional theory.The results indicated that the second sub-surface oxygen and the second bridge sites are the most stable adsorbed sites for CO and Pt on the ZrO2(101)surface,respectively.When the coverage is 0.25 ML(monolayer)the most stable models were obtained with adsorption energies of 56.2 and 352.7 kJ·mol-1.The most stable model of CO adsorbed on the Pt/t-ZrO2(101)surface was obtained by C-end adsorption with an adsorption energy of 323.8 kJ·mol-1.We considered vibrational frequency calculations,density of states and the Mulliken population of the adsorption systems before and after adsorption and the