Citation: XUE Ji-long, FANG Lei, LUO Wei, MENG Yue, CHEN Tao, XIA Sheng-jie, NI Zhe-ming. Density functional study of water gas shift reaction catalyzed by Cu-Pt-Au ternary alloy[J]. Journal of Fuel Chemistry and Technology, 2019, 47(6): 688-696.
Cu-Pt-Au三元合金催化水煤气变换反应的密度泛函研究
-
关键词:
- 密度泛函理论
- / 水煤气变换反应(WGSR)
- / 机理
- / 三元合金
English
Density functional study of water gas shift reaction catalyzed by Cu-Pt-Au ternary alloy
-
Key words:
- density functional theory
- / water gas shift reaction(WGSR)
- / mechanism
- / ternary alloy
-
氢能因其清洁、高效、安全和可再生等特点,被视为二十一世纪最具发展潜力的新能源。以往的H2都由天然气经重整反应得到[1],重整气中含有的CO易导致催化剂发生中毒现象,然而,水煤气变换反应(WGSR)不仅能够减低CO浓度,还可以提高重整气中H2的含量。如何高效的催化WGSR也成为了世界各国的研究热点。且WGSR在氨合成、质子膜燃料电池(PEMFC)、甲醇合成、煤气化等方面都有广泛的应用[2-4]。现今WGSR催化剂主要的研究热点是Pt-Au双金属催化剂,因其双功能机制和电子效应能够有效地提高催化活性和耐毒性,并且能克服CO的截面尺寸效应[5]而受到了广泛的关注。如Guo等[6]研究发现,Pt-Au合金相较于Pd/C和Pt/C具有更好的稳定性和反应活性,并且对CO表现出了较强的耐毒性。Thuy-Duong等[7]比较了不同含量的Au/Pt同时负载在CeO2和TiO2时对WGSR的催化活性,发现当Au/Pt负载量都为5%时,300 ℃以下CO的转换率接近90%,表明该催化剂对WGSR具有较佳的低温催化能力。Yu等[8]制备的Pt-Au/CeO2材料对CO转化反应具有较高的活性,在250 ℃下的转化率为78%,对WGSR的良好催化活性可能是由于负载Au和Pt之后能产生更多的活性位点。
为了提高催化剂活性及降低成本,考虑引入对WGSR也有较好催化活性的Cu[9, 10]以减少催化剂中Pt和Au的含量。相比二元合金,三元合金具有更多调整自身性质的可能性。三元合金常用共电沉积方法制备[11]。例如Wang等[5]采用两步法在玻碳电极(GCE)上制备了Cu-Pt-Au三元催化剂,且发现Cu-Pt-Au三元合金材料中存在合金化微结构,当Cu:Pt:Au物质的量比为1.82:1.11:1时,催化剂对甲醇电氧化具有较优的催化活性,表明Cu-Pt-Au三元合金催化剂具有用于直接甲醇燃料电池的巨大潜力。Sohn等[12]发现,Pt1-xAuxCu3/C对氧还原反应的活性为Pt/C材料的2.4倍,且Pt0.97Au0.3Cu3表现出较优的电催化活性、质量比活性和面积比活性。Brault等[13]利用密度泛函理论(DFT)对Pt70Pd15Au15和Pt50Pd25Au25催化剂生长进行了分子动力学的模拟,在构成团簇结构时,Au层会自发包覆在表面并最终形成Pt-Pd-Au的核壳结构,发现随着Au原子含量从15%升高到25%时,催化剂的氧还原反应活性也逐渐提高。Guo等[14]利用DFT研究了PtRuM(M=Fe、Ni、Cu、Mo、Sn)三元合金团簇上O2解离的反应机理,通过反应能垒表明三金属团簇的催化性能优于双金属团簇,对O2催化解离的效果更佳。Zhang等[15]通过一步化学还原法在石墨烯片(G)上合成了具有良好分散性的AuPtPd三元合金纳米粒子,发现AuPtPd/G对甲醇氧化反应具有优异的催化活性。Hong等[16]利用非均相及外延生长方法制备了树枝状的Au/Pt和Au/PtCu纳米线,Au/PtCu较Au/Pt而言对甲醇电氧化表现出更优的电催化活性。
以上研究表明,Cu-Pt-Au三元合金催化剂可能对甲醇氧化等反应具有较优的催化活性,但其对水煤气变换反应的催化活性的研究报道极少。为探讨不同掺杂量的三元催化剂对水煤气变化反应机理的影响,本研究设计了Cux-Pty-Au(111)(x, y=1、2、3)三元合金催化剂,利用DFT理论[17]对WGSR中的小分子吸附情况及反应机理进行研究,并且阐明WGSR三条机理在三元合金催化剂上的反应历程,以期能够为WGSR在实验上设计出更加高效的催化剂提供理论依据。
1. 计算方法和模型
本研究均在Dmol3模块下进行[18, 19],采用广义梯度积分和Perdew-Burke-Ernzerhof泛函相结合(GGA-PBE)的密度泛函理论方法。为了减小范德华力引起的吸附误差,本研究运用了DFT-D的方法。且催化剂金属内层电子选用有效核势赝势(ECP)代替,价电子波函数采用双数值基加轨道极化函数(DNP)展开。Brillouinzone积分的k点取样采用Monkhorst-Pack的自动生成方法,网格参数设为Medium,Methfessel-Paxton smearing设置为0.005 eV。结构优化以能量差异小于2×10-5 Ha、原子位移小于5.0×10-4 nm以及各个原子上力的收敛4.0×10-2 Ha/nm为判据。金属周期性平板结构的自旋极化现象对吸附构型及催化剂本身的性质影响很小,计算时选择自旋极化不受限制[20]。选择稳定吸附构型的反应物和产物分别作为各基元反应的初态(IS)和终态(FS),执行单个LST极大值计算并重复使用共轭梯度法优化和QST极大值进行过渡态(TS)搜索,得出其过渡态结构。精度设置为Medium,RMS convergence精度为0.1 Ha/nm。
根据谐波过渡态理论(HTST)计算每个基元反应的速率常数,遵循以下等式[21, 22]:
$k = \frac{{{k_{\text{b}}}T}}{h} $ ${\text{e}}\left( {\frac{{ - \Delta {G_{\text{a}}}}}{{RT}}} \right) $ ,公式里的k指代的是反应速率常数,kb是玻尔兹曼常数,T代表温度,h为普朗克常数,ΔGa是吉布斯反应自由能,R为气体常数。吉布斯自由能是振动能量,旋转及平动能量之和,并考虑零点校正(ZPE),并且焓变H和熵变S都是考虑平动、转动、振动的能量贡献,而吉布斯自由能可通过ΔGa=H-TS得出。基于课题组之前在M-Au(111)二元合金表面上的研究[23],且Wu等[24]考察了不同比例的Cu-Pt同时掺杂在Au基催化剂上时的情况,发现Cu:Pt比为1:1时,催化剂的能量最低既结构最为稳定。当Au、Cu和Pt分别负载在CeO2等载体上时,对WGSR都有较好的催化活性[7]。考虑计算的精确度及效率,本研究构建Cu:Pt=1:1的3×3周期性平板模拟Cux-Pty-Au(111)(x, y=0、1、2、3)面,图 1为已经优化好的Cu3-Pt3-Au(111)催化剂模型。
图 1
2. 结果与讨论
2.1 催化剂的性质
2.1.1 催化剂的稳定性
多元催化剂的结构稳定性和热力学稳定性对催化活性有较为显著的影响。原子尺寸因素ΔR可以很好的对催化剂结构稳定性进行衡量[25]。当ΔR≤ 0.15时,原子之间能够形成较为稳定的合金,且ΔR越小,引起合金的晶格畸变越小[26]。而多元催化剂的结合能ΔEb可以从热力学角度衡量催化剂的稳定性,结合能越大表明多元合金催化剂的热力学稳定性越好[27]。本研究定义原子尺寸因素[28]ΔR1=|RCu- RAu|/RAu,ΔR2=|RPt- RAu|/RAu,ΔR3=|RCu- RPt|/RPt。式中,RAu、RCu、RPt分别代表Au、Cu、Pt的原子半径,分别为0.179、0.157、0.183 nm。ΔR1、ΔR2、ΔR3分别表示Cu-Au、Pt-Au、Cu-Pt二元合金的原子尺寸因素。计算得ΔR1、ΔR2、ΔR3分别为0.123、0.022、0.142。由此可得,三种合金都可以稳定存在,且Pt-Au合金最为稳定。
本研究中定义结合能ΔEb=|Esurface- aEAu- bECu- cEPt|,其中, aEAu代表a个Au原子的能量,bECu代表b个Cu原子的能量,cEPt代表c个Pt原子的能量Esurface代表催化剂平面的总能量。结合能数据见表 1。当Cu、Pt分别掺杂在Au(111)面上时,结合能都逐渐增大,当掺杂量相同时(即x=y),ΔEb[Pty-Au(111)] > ΔEb[Cux-Au(111)]。其中, ΔEb[Cu3-Au(111)]为72.97 eV,而ΔEb[Pt3-Au(111)]为77.15 eV。因此, Pt-Au合金催化剂较Cu-Au催化剂具有更好的热力学稳定性。并且Cux-Pty-Au(111)面的结合能都大于二元金属表面,其中, Cu3-Pt3-Au的结合能最大,为77.99 eV。
表 1
Surface ΔEb/eV Surface ΔEb/eV Surface ΔEb/eV Surface ΔEb/eV Surface ΔEb/eV Cu1-Au 72.35 Pt3-Au 77.15 Cu1-Pt1-Au 74.05 Cu1-Pt2-Au 75.76 Cu1-Pt3-Au 77.52 Cu2-Au 72.70 Pt1-Au 73.61 Cu2-Pt1-Au 74.38 Cu2-Pt2-Au 76.07 Cu2-Pt3-Au 77.79 Cu3-Au 72.97 Pt2-Au 75.36 Cu3-Pt1-Au 74.62 Cu3-Pt2-Au 76.29 Cu3-Pt3-Au 77.99 2.1.2 态密度图分析
d带中心可以用于衡量过渡金属催化剂表面的电子活性。d带中心距离费米能级(X= 0 eV)越近,d带空穴越多,催化剂的电子活性越高[29],进而可以通过d带中心的大小来检测和验证不同金属表面的吸附能力。态密度分析见图 2,掺杂量不同的催化剂的d带中心见表 2。
图 2
表 2
Surface Surface Surface Surface Surface Cu1-Au -3.34 Pt1-Au -3.32 Cu1-Pt1-Au -3.26 Cu1-Pt2-Au -3.20 Cu1-Pt3-Au -3.14 Cu2-Au -3.29 Pt2-Au -3.25 Cu2-Pt1-Au -3.22 Cu2-Pt2-Au -3.15 Cu2-Pt3-Au -3.09 Cu3-Au -3.24 Pt3-Au -3.18 Cu3-Pt1-Au -3.17 Cu3-Pt2-Au -3.11 Cu3-Pt3-Au -3.05 由表 2可知,二元Pt-Au合金的d带中心明显都要比Cu-Au合金的d带中心更加靠近费米能级,具有更优的电子活性,与文献报道的一致[30]。在二元合金中,随着掺杂量的增加,合金催化剂的d带中心都逐渐向费米能级移动,其中, Cu3-Au(111)和Pt3-Au(111)面的d带中心值为-3.24、-3.18 eV。在三元合金中,当y值相同时,Cu3-Pty-Au(111)的d带中心比Cu1-Pty-Au(111)更加靠近费米能级0.09 eV,当x值相同时,Cux-Pt3-Au(111)面较Cux-Pt1-Au(111)面的d带中心更靠费米能级0.12 eV。由此可得,Pt掺杂量的增加更有利于提高Cu-Pt-Au三元合金催化剂的电子活性,其中, Cu3-Pt3-Au(111)的d带中心为-3.05 eV,电子活性最佳。
2.1.3 不同平板上CO和H2O的吸附
在15种催化剂表面上计算了CO和H2O的吸附能及键长。由表 3可得,Cu3-Pt3-Au(111)面上CO和H2O的最优吸附构型最为稳定,并且C-O及O-H之间的键长都符合实验文献所述[31]。CO分子在Cu3-Pt3-Au(111)面上的吸附最为稳定,吸附能为2.32 eV,相较于其他Cu-Au、Pt-Au等二元催化剂来说,CO在Cux-Au(111)面上的吸附最弱,但在引入了Pt元素后,对CO的吸附能明显增大,可能是由于C中的电子向过渡金属表面转移,使C-O键被削弱,从而加强了分子在表面的吸附。综上可知,对不同催化剂表面的结合能、稳定性分析、态密度图分析可知,Cu3-Pt3-Au(111)面在所研究的催化剂中结合较为稳定,且d带中心更加靠近费米能级,具有更优越的电子活性。故在以下的讨论中选取Cu3-Pt3-Au(111)来模拟三元合金催化剂表面,计算了CO及H2O小分子在表面吸附时的电荷分布,并探讨了水煤气变换反应在Cu3-Pt3-Au(111)面上的反应机理。
表 3
表 3 不同掺杂比例金属表面上的CO和H2O的吸附能及结构参数Table 3. Adsorption energy and structural parameters of CO and H2O on metal surfaces with different doping ratiosSurface-CO Site C-O
(nm)Eads
/eVSurface-H2O Site O-H1
(nm)O-H2
(nm)Bond angle Eads
/eVCu1-Au(111) Top-Cu 0.1150 -0.86 Cu1-Au(111) Top-Cu 0.0979 0.0979 104.144 -0.50 Cu2-Au(111) Bri-Cu 0.1174 -0.88 Cu2-Au(111) Top-Cu 0.0976 0.0978 104.462 -0.54 Cu3-Au(111) Bri-Cu 0.1171 -0.92 Cu3-Au(111) Bri-Cu 0.0995 0.1055 104.444 -0.55 Pt1- Au(111) Top-Pt 0.1156 -1.92 Pt1- Au(111) Top-Pt 0.0977 0.0977 104.552 -0.46 Pt2- Au(111) Bri-Pt 0.1185 -2.02 Pt2- Au(111) Top-Pt 0.0977 0.0976 104.446 -0.4 Pt3- Au(111) Bri-Pt 0.1184 -2.08 Pt3- Au(111) Bri-Pt 0.0971 0.0971 104.674 -0.18 Cu1-Pt1-Au(111) Top-Pt 0.1157 -1.97 Cu1-Pt1-Au(111) Top-Cu 0.0978 0.0977 104.570 -0.50 Cu1-Pt2-Au(111) Top-Pt 0.1157 -2.00 Cu1-Pt2-Au(111) Top-Cu 0.0976 0.0977 104.658 -0.53 Cu1-Pt3-Au(111) Top-Pt 0.1156 -1.86 Cu1-Pt3-Au(111) Top-Cu 0.0977 0.0982 103.876 -0.55 Cu2-Pt1-Au(111) Top-Pt 0.1157 -2.00 Cu2-Pt1-Au(111) Top-Cu 0.0978 0.0981 104.344 -0.49 Cu2-Pt2-Au(111) Top-Pt 0.1157 -2.00 Cu2-Pt2-Au(111) Top-Cu 0.0977 0.0983 104.054 -0.54 Cu2-Pt3-Au(111) Top-Pt 0.1158 -1.95 Cu2-Pt3-Au(111) Top-Cu 0.0978 0.0981 104.046 -0.50 Cu3-Pt1-Au(111) Fcc 0.1141 -2.01 Cu3-Pt1-Au(111) Cu-Pt-Bri 0.0958 0.1018 104.204 -0.54 Cu3-Pt2-Au(111) Top-Pt 0.1159 -2.02 Cu3-Pt2-Au(111) Top-Pt 0.0977 0.0985 102.304 -0.55 Cu3-Pt3-Au(111) Fcc 1.188 -2.32 Cu3-Pt3-Au(111) Cu-Pt-Bri 0.979 0.974 104.138 -0.54 2.2 反应机理
2.2.1 Cu3-Pt3-Au(111)表面的电荷布局分析
在吸附体系中,Mulliken电荷布居分析[32.33]可反映出物质中的原子与键在吸附前后的电子分布,电子转移数越多则说明小分子与金属表面的相互作用更强,吸附能更大。为了分析CO及H2O在Cu3-Pt3-Au(111)面上的吸附情况,对CO及H2O分子在表面的最优吸附构型进行了Mulliken电荷布居分析,具体见表 4。
表 4
Atom Free charge(e) Adsorption charge(e) CCO 0.096 0.611 OCO -0.096 -0.178 HH2O 0.245 0.266 HH2O 0.245 0.269 OH2O -0.490 -0.412 Total 0 0.556 吸附之前CO小分子显电中性,其中, O原子的电负性大于C原子,因此, 在吸附前的CO分子中,成键电子尤其是p电子倾向于O原子,所以O原子带负电荷,C原子带正电荷。CO的整体电荷从0.000 e分别增加为0.433 e,H2O的整体电荷从0.000 e增加为0.123 e,说明在吸附过程中,部分电子从CO、H2O分子转移到了金属表面,增强了CO及H2O在催化剂表面的吸附。
2.2.2 WGSR中小分子在Cu3-Pt3-Au(111)表面上的吸附
为了明确WGSR中各个小分子在催化剂表面上的吸附程度,本研究引入吸附能(Eads)来衡量Eads= E(A/surface) - EA - Esurface,其中, E(A/surface)为小分子吸附在催化剂表面上时体系的总能量,EA表示小分子的能量,Esurface表示Cu3-Pt3-Au(111)面的总能量。Eads的符号和数值大小可衡量发生吸附的可能性和吸附作用的强弱,其值越大表示吸附过程中放出热量越多,形成的吸附体系越稳定。水煤气变换反应中共有10个吸附物质,分别为:CO、H2O、OH、H、O、COOH、CO2、H2、HCO、O和CHO,在Cu3-Pt3-Au(111)面吸附位的单点吸附共有70种。分别将10个吸附物质在Cu3-Pt3-Au(111)面上吸附的结构进行优化计算,通过对Eads的比较,得到10个吸附物种在催化剂上的最优吸附构型。吸附能如表 5所示,10种吸附小分子在Cu3-Pt3-Au(111)面的Cu-Pt-Bri和Fcc处的吸附能都普遍较大。当CO以碳端吸附在Fcc的时候吸附能最大值为2.32 eV,H2O吸附在Cu-Pt-Bri时的吸附能最大,其值为0.54 eV,并且CO和H2O分子的键长都稍有拉长,其中, 碳氧键拉长了0.0047 nm,O-H键拉长了0.0009 nm。并且在Pt原子附近的吸附位置都明显较Cu原子附近的吸附位的吸附能更大,所得到的吸附构型更加稳定。
表 5
表 5 Cu3-Pt3-Au(111)表面上小分子的吸附能Table 5. Adsorption energy of different small molecules on the Cu3-Pt3-Au(111) surfaceMolecule Eads/eV Top-Cu Top-Pt Bri-Cu Bri-Pt Cu-Pt-Bri Hcp Fcc CO -2.12 -1.94 -1.94 -2.40 -2.32 -2.18 -2.32 CO2 -0.09 -0.11 -0.09 -0.02 -0.12 -0.11 -0.10 O -4.38 -4.39 -4.54 -4.38 -4.54 -4.29 -4.54 H -2.79 -2.49 -2.39 -3.03 -3.03 -2.88 -2.89 H2O -0.37 -0.39 -0.27 -0.39 -0.54 -0.66 -0.26 COOH -1.78 -2.68 -1.79 -2.71 -2.77 -1.54 -1.68 OH -2.50 -2.97 -2.98 -2.71 -2.98 -2.75 -2.99 H2 -0.06 -0.07 -0.05 -0.07 -0.06 -0.11 -0.08 HCOO -2.61 -2.66 -2.50 -2.66 -2.66 -2.45 -2.65 CHO -1.47 -2.63 -1.47 -2.63 -2.61 -2.75 -2.69 2.2.3 WGSR机理
水煤气变换反应可按氧化还原机理、羧基机理、甲酸机理三种机理进行[34],并计算了WGSR中各个基元反应的能垒及反应热,见表 6,具体反应路径见图 3。
表 6
表 6 Cu3-Pt3-Au(111)表面上各基元反应的活化能(Ea)、反应热(ΔE)和速率常数kTable 6. Energy barriers, reaction heat and k of primer reactions on Cu3-Pt3-Au(111)Mechanism Element reaction Cu3-Pt3-Au(111) Ea/eV ΔE/eV k/(cm3·mol-1·s-1) 1: H2O+*=H*+OH* 0.02 0.01 4.85×1014 Redox mechanism 2-a: OH*=O*+H* 1.24 0.49 6.35×1013 2-b: OH*+OH*=H2O*+O* 1.36 1.35 5.33×1012 3: CO*+O*=CO2*+* 4.84 -0.05 8.30×1012 Carboxyl mechanism 4: CO*+OH*=COOH*+* 3.15 0.48 1.99×1011 5-a: COOH*+*=H*+CO2* 0.19 -0.19 2.29×1012 5-b: COOH*+OH*=H2O*+CO2* 1.07 -0.20 1.67×1012 Formate mechanism 2-a: OH*=O*+H* 1.24 0.49 6.35×1013 6: CO*+H*=CHO* 0.51 0.05 2.85×1012 7: CHO*+O*=HCOO* 0.85 -0.02 4.14×1013 8: HCOO*+*=CO2*+H* 3.09 -0.24 7.03×1012 9: 2H*=H2*+2* 0.43 -0.08 6.70×1013 图 3
由图 3可知,在氧化还原机理中,水解离后产生的OH经过2-a继续解离。以O端竖直吸附在Cu-Pt-Bri位上的OH解离成为游离的O、H,并平行吸附在Cu-Pt-Bri位的两侧,相距1.386 nm,该步骤需跨越1.24 eV的能垒,吸收热量0.49 eV,这是O原子的来源之一。另外,OH可以在催化剂表面发生歧化反应产生O,具体反应为2-b,该步骤能垒为1.36 eV,需要吸收热量1.35 eV,最终产生的H2O吸附在Top-Cu位,O原子吸附在Cu-Pt-Bri位。之后O原子在逐渐靠向CO分子,形成CO2最终吸附在Cu-Pt-Bri位上,该步骤为氧化还原机理中的速控步骤,需要跨过4.84 eV大小的能垒,并放出0.05 eV的热量,产生的CO2随后脱附,两个活性吸附位上H靠近形成H2,最终得到H2。在羧基机理中,水解离产生的OH吸附在Cu-Bri位,随后OH逐渐向Pt-Bri位置的CO靠近,形成COOH中间体,最终以C端吸附在催化剂表面,该步骤需要跨过3.15 eV能垒,且为吸热反应,吸热0.48 eV。形成CO2存在两种可能的反应,在图 3的5-a中,中间体COOH继续自解离成CO2和H,CO2逐渐远离催化剂表面,其中, C原子距表面0.36145 nm,该反应需克服0.19 eV的能垒,并放出0.19 eV的热量。但在图 3的5-b中,COOH中的O-H断裂,H向吸附在Cu-Bri位置上的OH逐渐靠近,形成的H2O脱附,OCO中间体则转化为CO2,并且图 3的5-b中的反应产物都逐渐远离催化剂表面,最终脱附。
在甲酸机理中,CO与一个H原子相互靠近形成CHO中间产物并以C、O平行吸附在Fcc位,该基元反应需克服0.51 eV的能垒,并吸热0.05 eV。之后O向CHO中的C端移动,最终形成C-O键,构成中间产物HCOO,能垒为0.85 eV,该反应为放热反应,释放0.02 eV热量。在反应8中HCOO发生自解离,解离出的CO2在表面脱附,H则和其他活性位点上的H结合成为H2,进而脱附,该步需克服3.09 eV的能垒,并放热0.24 eV。具体的过渡态见图 4所示。各机理能量变化如图 5所示。
图 4
图 5
3. 结论
本研究计算了不同掺杂量的二元合金(Cu-Au, Pt-Au)及三元合金(Cu-Pt-Au)的性质,表明在二元合金中Pty-Au较Cux-Au具有更好的稳定性,当x=y时,ΔEb[Pty-Au(111)] > ΔEb[Cux-Au(111)]成立。Pty-Au(111)与Cux-Au(111)相比较而言其d带中心更靠近费米能级,表现更优的电子活性。当再引入第三金属构成Cu-Pt-Au催化剂时,稳定性有较大提高,Cu3-Pt3-Au(111)面最为稳定,结合能为77.99 eV,结合能大于二元合金催化剂,且Cu3-Pt3-Au(111)的d带中心最靠近费米能级,为-3.05 eV,电子活性最佳。
分析了15种不同金属平面对CO及H2O的吸附构型及吸附能,发现Cu3-Pt3-Au(111)对CO及H2O有较强的吸附,更有利于水煤气变换反应进行。
讨论了WGSR中三条机理在Cu3-Pt3-Au(111)上的反应历程,发现在三元合金表面水解离反应极易发生。而反应3、反应4和反应8则分别作为氧化还原机理、羧基机理和甲酸机理中的速控步骤,其中, 氧化还原机理中形成TS3需要克服的能垒值最大,为4.84 eV,而羧基机理和甲酸机理的能垒分别为3.15、3.09 eV。甲酸机理相对比其他两种机理而言更容易在Cu3-Pt3-Au(111)上发生反应,具体反应路径为:1→2-a→6→7→8→9。并发现存在CO或OCO中间体的基元反应能垒值较大,可能是Pt与CO的强吸附作用。
-
-
[1]
JAMES J S. Catalysis in the development of clean energy technologies[J]. Catal Today, 2005, 100(1): 171-180.
-
[2]
ROZOVSKII A Y, LIN G I. Fundamentals of methanol synthesis and decomposition[J]. Top Catal, 2003, 22(3): 137-150.
-
[3]
AMMAL S C, HEYDEN A. Water-gas shift catalysis at corner atoms of Pt clusters in contact with a TiO2 (110) support surface[J]. ACS Catal, 2014, 4(10): 3654-3662. doi: 10.1021/cs5009706
-
[4]
CHOI Y, STENGER H G. Fuel cell grade hydrogen from methanol on a commercial Cu/ZnO/Al2O3 catalyst[J]. Appl Catal B:Environ, 2002, 38(4): 259-269. doi: 10.1016/S0926-3373(02)00054-1
-
[5]
WANG C Q, REN F F, ZHAI C Y, ZHANG K, YANG B B, BIN D, WANG H W, YANG P, DU Y K. Au-Cu-Pt ternary catalyst fabricated by electrodeposition and galvanic replacement with superior methanol electrooxidation activity[J]. Rsc Adv, 2014, 4(11): 57600-57607.
-
[6]
GUO Z, LIU B, ZHANG Q H, DENG W P, WANG Y, YANG Y H. Recent advances in heterogeneous selective oxidation catalysis for sustainable chemistry[J]. Chem Soc Rev, 2014, 43(1): 3480-3524.
-
[7]
THUY-DUONG N, ASHLEIGH E B, JOSÉ A R, SANJAYA D S. Au and Pt nanoparticle supported catalysts tailored for H2 production:From models to powder catalysts[J]. Appl Catal A:Gen, 2016, 518(25): 18-47.
-
[8]
YU Q Q, CHEN W, LI Y, JIN Y, SUO Z H. The action of Pt in bimetallic Au-Pt/CeO2 catalyst for water-gas shift reaction[J]. Catal Today, 2010, 158(3): 324-328.
-
[9]
KONG G X, MA X J, LIU Q J, LI Y, LIU Z T. Structural stability, elastic and thermodynamic properties of Au-Cu alloys from first-principles calculations[J]. Physica B, 2018, 533(15): 58-62.
-
[10]
MUHAMMAD A S, AKHTAR H, MUHAMMAD S, ORTWIN L, ALEXANDRE A. DFT study of synergistic catalysis of the water-gas-shift reaction on Cu-Au bimetallic surfaces[J]. ChemCatChem, 2016, 8(6): 1208-1217. doi: 10.1002/cctc.201501312
-
[11]
XU Y, ZHANG B. Recent advances in porous Pt-based nanostructures:Synthesis and electrochemical applications[J]. Chem Soc Rev, 2014, 43(8): 2439-2450. doi: 10.1039/c3cs60351b
-
[12]
SOHN Y, JI B J, PIL K. Chemically dealloyed Pt-Au-Cu ternary electrocatalysts with enhanced stability in electrochemical oxygen reduction[J]. Res Chem Interm, 2018, 44(6): 3697-3712. doi: 10.1007/s11164-018-3375-3
-
[13]
BRAULT P, COUTANCEAU C, JENNINGS P C, VEGGE T, BERNDT J, CAILLARD A, BARANTON S, LANKIANG S. Molecular dynamics simulations of ternary PtxPdyAuz fuel cell nanocatalyst growth[J]. Int J Hydrogen Energy, 2016, 41(47): 22589-22597. doi: 10.1016/j.ijhydene.2016.08.035
-
[14]
GUO W L, LIAN X, XIAO P, LIU F L, YANG Y, ZHANG Y H, ZHANG X X. DFT studies on the interaction of PtxRuyMz(M=Fe, Ni, Cu, Mo, Sn, x+y+z=4, x ≥ 1, y ≥ 1) alloy clusters with O2[J]. Mol Phys, 2015, 113(8): 854-865. doi: 10.1080/00268976.2014.983573
-
[15]
ZHANG Y Z, GU Y, LIN S X, WEI J P, WANG Z H, WANG C M, DU Y L, YE W C. One-step synthesis of PtPdAu ternary alloy nanoparticles on graphene with superior methanol electrooxidation activity[J]. Electrochim Acta, 2011, 56(24): 8746-8751. doi: 10.1016/j.electacta.2011.07.094
-
[16]
HONG W, WANG J, WANG E. Dendritic Au/Pt and Au/PtCu nanowires with enhanced electrocatalytic activity for methanol electrooxidation[J]. Small, 2014, 10(16): 3262-3265. doi: 10.1002/smll.v10.16
-
[17]
GHOLIZADEH R, YU Y X. N2O + CO reaction over Si-and Se-doped graphenes:An ab initio DFT study[J]. Appl Surf Sci, 2015, 357(1): 1187-1195.
-
[18]
张田, 郭琛, 魏淑贤, 武中华, 韩兆翔, 鲁效庆. 甲硫醇在Co修饰MoS2团簇边缘位的脱硫机理研究[J]. 化学学报, 2018,76,(1): 62-67. ZHANG Tian, GUO Chen, WEI Shu-xian, WU Zhong-hua, HAN Zhao-xiang, LU Xiao-qing. Investigation on CH3SH desulfurization mechanism at the edge site of Co-doped MoS2 cluster[J]. Acta Chim Sin, 2018, 76(1): 62-67.
-
[19]
钱梦丹, 薛继龙, 夏盛杰, 倪哲明, 蒋军辉, 曹勇勇. Pd/Cu(111)双金属表面催化糠醛脱碳及加氢的反应机理[J]. 燃料化学学报, 2017,45,(1): 34-42. doi: 10.3969/j.issn.0253-2409.2017.01.006QIAN Meng-dan, XUE Ji-long, XIA Sheng-jie, NI Zhe-ming, JIANG Jun-hui, CAO Yong-yong. Decarbonylation and hydrogenation reaction of furfural on Pd/Cu(111) surface[J]. J Fuel Chem Technol, 2017, 45(1): 34-42. doi: 10.3969/j.issn.0253-2409.2017.01.006
-
[20]
GOVIND N, PETERSEN M, FITZGERAL G, KING-SMITH D, ANDZELM J. A generalized synchronous transit method for transition state location[J]. Comput Mater Sci, 2003, 28(2): 250-258. doi: 10.1016/S0927-0256(03)00111-3
-
[21]
VINEYARD G H. Frequency factors and isotope effects in solid state rate processes[J]. J Phys Chem Solids, 1957, 3(1): 121-127.
-
[22]
PISKORZ W, ZASADA F, STELMACHOWSKI P, DIWALD O, KOTARBA A, SOJKA Z. Computational and experimental investigations into N2O decomposition over MgO nanocrystals from thorough molecular mechanism to Ab initio microkinetics[J]. J Phys Chem C, 2011, 115(45): 22451-22460. doi: 10.1021/jp2070826
-
[23]
蒋军辉, 钱梦丹, 薛继龙, 夏盛杰, 倪哲明, 邵蒙蒙. In-Au(111)和Ir-Au(111)合金表面的性质及其对巴豆醛的吸附比较[J]. 物理化学学报, 2016,32,(12): 2932-2940. doi: 10.3866/PKU.WHXB201609302JIANG Jun-hui, QIAN Meng-dan, XUE Ji-long, XIA Sheng-jie, NI Zhe-ming, SHAO Meng-meng. Comparison of properties of In-Au(111) and Ir-Au(111) alloy surfaces, and their adsorption to crotonaldehyde[J]. Acta Phys-Chim Sin, 2016, 32(12): 2932-2940. doi: 10.3866/PKU.WHXB201609302
-
[24]
WU G H, SUN Y, WU X, CHEN R, WANG Y. Large scale structural optimization of trimetallic Cu-Au-Pt clusters up to 147 atoms[J]. Chem Phys Lett, 2017, 686(16): 103-110.
-
[25]
陶杰, 姚正军, 薛烽.材料科学基础[M].北京:化学工业出版社, 2006, 50-51.TAO Jie, YAO Zheng-jun, XUE Feng. Fundamentals of Material Science[M]. Beijing:Chemical Industry Press, 2006, 50-51.
-
[26]
SUSAN M A, RYAN L A, KOJI S, RYO K, KAZUKI K, NGUYEN H L, HIROSHI N, HIDEAKI K. First principles calculations of transition metal binary alloys:Phase stability and surface effects[J]. J Electron Mater, 2017, 46(6): 3776-3783. doi: 10.1007/s11664-017-5402-3
-
[27]
黄敏, 徐畅, 程龙玖. [BxAl13-x]-(x=0~13)二元团簇的密度泛函理论研究[J]. 化学学报, 2016,74,(9): 758-763. HUANG Min, XU Chang, CHENG Long-jiu. Density functional theory studies of the binary systems[BxAl13-x]-(x=0-13)[J]. Acta Chim Sin, 2016, 74(9): 758-763.
-
[28]
钱梦丹, 罗伟, 倪哲明, 夏盛杰, 薛继龙, 蒋军辉. Ru修饰前后Pd(111)面的性质及对糠醛吸附的比较研究[J]. 高等学校化学学报, 2017,38,(9): 1611-1618. QIAN Meng-dan, LUO Wei, NI Zhe-ming, XIA Sheng-jie, XUE Ji-long, JIANG Jun-hui. Comparative study on the properties and adsorption of furfural of Pd(111) surface before and after Ru modification[J]. Chem J Chin Univ, 2017, 38(9): 1611-1618.
-
[29]
CHANG F F, SHAN S Y, VALERI P, ZAKIYA S, AOLIN L, JONATHAN R, WU J F, LUO J, YU G, REN Y, ZHONG C J. Composition tunability and (111)-dominant facets of ultrathin platinum-gold alloy nanowires toward enhanced electrocatalysis[J]. J Am Chem Soc, 2016, 138(37): 12166-12175. doi: 10.1021/jacs.6b05187
-
[30]
PALOTAS K, BAKO I, BUGYI L. Structural, electronic and adsorption properties of Rh(111)/Mo(110) bimetallic catalyst:A DFT study[J]. Appl Surf Sci, 2016, 389(15): 1094-1103.
-
[31]
JALILI S, ISFAHANI-ZEINI A, HABIBPOUR R. DFT investigations on the interaction of oxygen reduction reaction intermediates with Au (100) and bimetallic Au/M (100) (M=Pt, Cu, and Fe) surfaces[J]. Comput Theor Chem, 2013, 4(1): 18-26.
-
[32]
HAMED A, MAJID V. Computational designing ultra-sensitive nano-composite based on boron doped and CuO decorated graphene to adsorb H2S and CO gaseous molecules[J]. Mater Res Express, 2017, 4(7): 1-7.
-
[33]
LIU R Q. Adsorption and dissociation of H2O on Au(111) surface:A DFT study[J]. Comput Theor Chem, 2013, 1019(1): 141-145.
-
[34]
AMIT A G, JAMES A D, MANOS M. On the mechanism of low-temperature water gas shift reaction on copper[J]. J Am Chem Soc, 2008, 130(4): 1402-1414. doi: 10.1021/ja0768237
-
[1]
-
表 1 不同掺杂比例催化剂的结合能
Table 1. Binding energy of metal surfaces with different doping amounts
Surface ΔEb/eV Surface ΔEb/eV Surface ΔEb/eV Surface ΔEb/eV Surface ΔEb/eV Cu1-Au 72.35 Pt3-Au 77.15 Cu1-Pt1-Au 74.05 Cu1-Pt2-Au 75.76 Cu1-Pt3-Au 77.52 Cu2-Au 72.70 Pt1-Au 73.61 Cu2-Pt1-Au 74.38 Cu2-Pt2-Au 76.07 Cu2-Pt3-Au 77.79 Cu3-Au 72.97 Pt2-Au 75.36 Cu3-Pt1-Au 74.62 Cu3-Pt2-Au 76.29 Cu3-Pt3-Au 77.99 表 2 不同掺杂比例Au(111)的d带中心
Table 2. d-band center of different doping ratio Au(111)
Surface Surface Surface Surface Surface Cu1-Au -3.34 Pt1-Au -3.32 Cu1-Pt1-Au -3.26 Cu1-Pt2-Au -3.20 Cu1-Pt3-Au -3.14 Cu2-Au -3.29 Pt2-Au -3.25 Cu2-Pt1-Au -3.22 Cu2-Pt2-Au -3.15 Cu2-Pt3-Au -3.09 Cu3-Au -3.24 Pt3-Au -3.18 Cu3-Pt1-Au -3.17 Cu3-Pt2-Au -3.11 Cu3-Pt3-Au -3.05 表 3 不同掺杂比例金属表面上的CO和H2O的吸附能及结构参数
Table 3. Adsorption energy and structural parameters of CO and H2O on metal surfaces with different doping ratios
Surface-CO Site C-O
(nm)Eads
/eVSurface-H2O Site O-H1
(nm)O-H2
(nm)Bond angle Eads
/eVCu1-Au(111) Top-Cu 0.1150 -0.86 Cu1-Au(111) Top-Cu 0.0979 0.0979 104.144 -0.50 Cu2-Au(111) Bri-Cu 0.1174 -0.88 Cu2-Au(111) Top-Cu 0.0976 0.0978 104.462 -0.54 Cu3-Au(111) Bri-Cu 0.1171 -0.92 Cu3-Au(111) Bri-Cu 0.0995 0.1055 104.444 -0.55 Pt1- Au(111) Top-Pt 0.1156 -1.92 Pt1- Au(111) Top-Pt 0.0977 0.0977 104.552 -0.46 Pt2- Au(111) Bri-Pt 0.1185 -2.02 Pt2- Au(111) Top-Pt 0.0977 0.0976 104.446 -0.4 Pt3- Au(111) Bri-Pt 0.1184 -2.08 Pt3- Au(111) Bri-Pt 0.0971 0.0971 104.674 -0.18 Cu1-Pt1-Au(111) Top-Pt 0.1157 -1.97 Cu1-Pt1-Au(111) Top-Cu 0.0978 0.0977 104.570 -0.50 Cu1-Pt2-Au(111) Top-Pt 0.1157 -2.00 Cu1-Pt2-Au(111) Top-Cu 0.0976 0.0977 104.658 -0.53 Cu1-Pt3-Au(111) Top-Pt 0.1156 -1.86 Cu1-Pt3-Au(111) Top-Cu 0.0977 0.0982 103.876 -0.55 Cu2-Pt1-Au(111) Top-Pt 0.1157 -2.00 Cu2-Pt1-Au(111) Top-Cu 0.0978 0.0981 104.344 -0.49 Cu2-Pt2-Au(111) Top-Pt 0.1157 -2.00 Cu2-Pt2-Au(111) Top-Cu 0.0977 0.0983 104.054 -0.54 Cu2-Pt3-Au(111) Top-Pt 0.1158 -1.95 Cu2-Pt3-Au(111) Top-Cu 0.0978 0.0981 104.046 -0.50 Cu3-Pt1-Au(111) Fcc 0.1141 -2.01 Cu3-Pt1-Au(111) Cu-Pt-Bri 0.0958 0.1018 104.204 -0.54 Cu3-Pt2-Au(111) Top-Pt 0.1159 -2.02 Cu3-Pt2-Au(111) Top-Pt 0.0977 0.0985 102.304 -0.55 Cu3-Pt3-Au(111) Fcc 1.188 -2.32 Cu3-Pt3-Au(111) Cu-Pt-Bri 0.979 0.974 104.138 -0.54 表 4 Mulliken电荷布局分析
Table 4. Mulliken charge population analysis
Atom Free charge(e) Adsorption charge(e) CCO 0.096 0.611 OCO -0.096 -0.178 HH2O 0.245 0.266 HH2O 0.245 0.269 OH2O -0.490 -0.412 Total 0 0.556 表 5 Cu3-Pt3-Au(111)表面上小分子的吸附能
Table 5. Adsorption energy of different small molecules on the Cu3-Pt3-Au(111) surface
Molecule Eads/eV Top-Cu Top-Pt Bri-Cu Bri-Pt Cu-Pt-Bri Hcp Fcc CO -2.12 -1.94 -1.94 -2.40 -2.32 -2.18 -2.32 CO2 -0.09 -0.11 -0.09 -0.02 -0.12 -0.11 -0.10 O -4.38 -4.39 -4.54 -4.38 -4.54 -4.29 -4.54 H -2.79 -2.49 -2.39 -3.03 -3.03 -2.88 -2.89 H2O -0.37 -0.39 -0.27 -0.39 -0.54 -0.66 -0.26 COOH -1.78 -2.68 -1.79 -2.71 -2.77 -1.54 -1.68 OH -2.50 -2.97 -2.98 -2.71 -2.98 -2.75 -2.99 H2 -0.06 -0.07 -0.05 -0.07 -0.06 -0.11 -0.08 HCOO -2.61 -2.66 -2.50 -2.66 -2.66 -2.45 -2.65 CHO -1.47 -2.63 -1.47 -2.63 -2.61 -2.75 -2.69 表 6 Cu3-Pt3-Au(111)表面上各基元反应的活化能(Ea)、反应热(ΔE)和速率常数k
Table 6. Energy barriers, reaction heat and k of primer reactions on Cu3-Pt3-Au(111)
Mechanism Element reaction Cu3-Pt3-Au(111) Ea/eV ΔE/eV k/(cm3·mol-1·s-1) 1: H2O+*=H*+OH* 0.02 0.01 4.85×1014 Redox mechanism 2-a: OH*=O*+H* 1.24 0.49 6.35×1013 2-b: OH*+OH*=H2O*+O* 1.36 1.35 5.33×1012 3: CO*+O*=CO2*+* 4.84 -0.05 8.30×1012 Carboxyl mechanism 4: CO*+OH*=COOH*+* 3.15 0.48 1.99×1011 5-a: COOH*+*=H*+CO2* 0.19 -0.19 2.29×1012 5-b: COOH*+OH*=H2O*+CO2* 1.07 -0.20 1.67×1012 Formate mechanism 2-a: OH*=O*+H* 1.24 0.49 6.35×1013 6: CO*+H*=CHO* 0.51 0.05 2.85×1012 7: CHO*+O*=HCOO* 0.85 -0.02 4.14×1013 8: HCOO*+*=CO2*+H* 3.09 -0.24 7.03×1012 9: 2H*=H2*+2* 0.43 -0.08 6.70×1013 -
扫一扫看文章
计量
- PDF下载量: 0
- 文章访问数: 0
- HTML全文浏览量: 0

下载:
下载: