反應(yīng)pH值對(duì)微波法合成的Mo修飾的Pt/C催化劑結(jié)構(gòu)和乙醇電氧化催化性能的影響

馬 忠 原鮮霞 夏小蕓 杜 娟 李 琳 馬紫峰

(上海交通大學(xué)化學(xué)工程系,上海200240)

1 引言

質(zhì)子交換膜燃料電池因其極高的能量密度和低溫運(yùn)行特性等優(yōu)點(diǎn)近年來備受關(guān)注.1,2鑒于所用常規(guī)燃料——?dú)錃獾纳a(chǎn)、儲(chǔ)存和運(yùn)輸?shù)葐栴},液體醇類(如甲醇、乙醇等)成為替代燃料之一.其中,甲醇存在一定的毒性和不穩(wěn)定性,且乙醇具有更高的能量密度.3因此,直接乙醇燃料電池(DEFC)已成為低溫燃料電池的研究熱點(diǎn).4-6

Pt催化劑已被廣泛用于直接乙醇燃料電池中作為乙醇氧化的催化劑.3,7-9然而,該類催化劑盡管具備優(yōu)秀的催化活性,但其表面很容易吸附乙醇氧化過程中產(chǎn)生的中間產(chǎn)物類CO物種,從而發(fā)生中毒現(xiàn)象而失去活性.10同時(shí),昂貴的成本也是制約其大規(guī)模應(yīng)用的主要原因之一.因此,開發(fā)低成本、高效率的乙醇電氧化催化劑已成為DEFC方向的一個(gè)重要任務(wù).為此,多種金屬如Ir、Mo、Os、Ru、Sn、W和Zr等已被引入Pt催化劑形成二元、三元甚至多元的Pt基催化劑,并對(duì)乙醇的電催化氧化表現(xiàn)出良好的性能.10-15其中,Mo元素盡管沒有電催化乙醇的作用,但具備較好的催化一氧化碳氧化的作用,16,17Samjeské等18利用電化學(xué)質(zhì)譜發(fā)現(xiàn)Mo元素對(duì)弱吸附態(tài)的CO具有很好的氧化作用.

微波法具有快速、高效、均勻等優(yōu)點(diǎn),可用于制備顆粒細(xì)小、均勻的納米粒子,已被廣泛應(yīng)用于有機(jī)合成及無機(jī)材料制備等方面.19-21另?yè)?jù)文獻(xiàn)報(bào)道,前驅(qū)體溶液的pH值對(duì)所制得催化劑的粒徑有重要影響,從而影響到催化劑的催化活性.22,23Chu等24采用微波法制備了PtRu/C催化劑,發(fā)現(xiàn)隨著pH值的增加,催化劑顆粒的粒徑逐漸變小.Zheng等25的研究表明,pH值為7時(shí)利用微波法制得的Ru85Se15/MWCNTs催化劑顆粒分散最均勻.

因此,本文采用微波法,通過調(diào)節(jié)不同反應(yīng)pH值制備了一系列Mo修飾的Pt/C催化劑,對(duì)比考察了反應(yīng)pH值對(duì)所制得的催化劑的微觀結(jié)構(gòu)、形貌及乙醇電氧化催化性能的影響,并對(duì)研究結(jié)果進(jìn)行了討論.

2 實(shí)驗(yàn)部分

2.1 催化劑的制備

Mo修飾的Pt/C催化劑采用微波法制備,具體合成方法為:將濃度分別為0.0488和0.0386 mol·L-1的H2PtCl6(國(guó)藥,分析純)和Na2MoO4(國(guó)藥,分析純)的乙二醇溶液加入到Vulcan XC-72(Cabot,美國(guó))的乙二醇(國(guó)藥,分析純)漿液中,超聲混合均勻后利用HCl(國(guó)藥,分析純)和NaOH(國(guó)藥,分析純)調(diào)節(jié)混合液的pH值(本文選取1、8、12和14四個(gè)pH值作為考察對(duì)象),然后將之置于微波反應(yīng)器(Sanle WP650D)中,在功率為325 W、冷凝回流的狀態(tài)下反應(yīng)15 min,然后洗去氯離子,真空干燥后得產(chǎn)物.

本文所制備的Mo修飾的Pt/C催化劑中,Pt:Mo的摩爾比均為2:1,Pt的質(zhì)量分?jǐn)?shù)為20%(w).

2.2 催化劑的表征

催化劑的微觀形貌采用日本電子株式會(huì)社的JEOL JEM 2000 EX透射電子顯微鏡(TEM)進(jìn)行觀察.

催化劑的晶體結(jié)構(gòu)采用日本理學(xué)Rigaku D/MAX PC 2000型X射線衍射(XRD)儀器進(jìn)行表征,測(cè)試條件為:CuKa射線源,管電壓為35 kV,管電流為25 mA,掃描速率為6(°)·min-1.

催化劑的表面形態(tài)采用美國(guó)PHI公司的PHI 5000C ESCA System儀器進(jìn)行測(cè)試.測(cè)試條件為鋁/鎂靶,高壓14.0 kV,功率250 W,通能93.9 eV,真空優(yōu)于1.33×10-6Pa.采用美國(guó)RBD公司的RBD147數(shù)據(jù)采集卡和AugerScan 3.21軟件分別采集樣品的0-1200(1100)eV的全掃描譜,而后采集各元素相關(guān)軌道的窄掃描譜.以C 1s(284.5 eV)為基準(zhǔn)進(jìn)行結(jié)合能校正.

2.3 催化劑的電化學(xué)性能表征

催化劑對(duì)乙醇電氧化過程的催化性能用循環(huán)伏安曲線、計(jì)時(shí)電流法和電化學(xué)阻抗譜進(jìn)行測(cè)試.測(cè)試采用常規(guī)的三電極電解池,其中以覆有所制備催化劑的玻碳電極為工作電極、鉑絲為對(duì)電極、Hg/Hg2SO4電極為參比電極、0.5 mol·L-1H2SO4和 0.5 mol·L-1H2SO4+0.5 mol·L-1C2H5OH為電解液.工作電極的具體制備方法為:取5 mg催化劑,加入一定配比的水和5%(w)的Nafion乳液配成5 g·L-1的催化劑懸乳液,超聲分散10 min后移取10 μL涂于有效面積為0.096 cm2的玻碳電極表面,干燥后得到工作電極.循環(huán)伏安測(cè)試的電位掃描速率為50 mV·s-1;計(jì)時(shí)電流測(cè)試采用恒電位0.273 V(vsHg/Hg2SO4),時(shí)間為1000 s;交流阻抗測(cè)試采用美國(guó)阿美特克公司的Solartron SI1287電化學(xué)工作站進(jìn)行,頻率掃描范圍為100 kHz-1 Hz,交流電壓振幅為5 mV,測(cè)試電極電位為0.8 V.

3 結(jié)果與討論

圖1為不同反應(yīng)pH值所制備的Mo修飾Pt/C催化劑的TEM圖.圖2為TEM圖中催化劑組分顆粒對(duì)應(yīng)的粒徑分布圖.從圖1和圖2可以看出,當(dāng)pH值為1時(shí),負(fù)載在碳載體上的催化劑顆粒數(shù)目較少、粒徑較大且分布范圍較廣,說明強(qiáng)酸環(huán)境不利于活性組分在碳載體上的均勻負(fù)載;當(dāng)提高pH值為8時(shí),催化劑顆粒的數(shù)目有所增加,但顆粒尺寸仍舊較大且粒徑分布較廣;當(dāng)pH值進(jìn)一步提升為12時(shí),碳載體上催化劑組分的顆粒較細(xì)小,但局部有所團(tuán)聚,不利于催化劑活性的充分發(fā)揮；當(dāng)提高pH值為14時(shí),催化劑組分在碳載體上分布均勻、且粒徑細(xì)小.這些結(jié)果表明,反應(yīng)pH值對(duì)微波法合成的催化劑中催化活性組分的粒徑大小與分布具有可控性的調(diào)節(jié)作用,強(qiáng)堿環(huán)境有利于粒徑細(xì)小的催化活性組分在碳載體上的均勻分布.這可能與催化劑合成過程中乙二醇氧化的產(chǎn)物乙醇酸和乙醇酸鹽有關(guān).其中,乙醇酸鹽對(duì)金屬離子具備很強(qiáng)的穩(wěn)定作用,能夠減緩金屬在碳載體上的沉積,從而使得形成的金屬顆粒尺寸細(xì)小并均勻分散.23當(dāng)溶液的pH值較高時(shí),溶液中乙醇酸鹽的濃度也較高,因而有利于制得粒徑細(xì)小且均勻分布的金屬粒子.當(dāng)溶液的pH值較低時(shí),溶液中乙醇酸鹽的濃度較低,乙醇酸的濃度較高.而乙醇酸對(duì)金屬粒子的穩(wěn)定作用很微弱,因而制備的金屬粒子的粒徑較大、分布也不均勻.

圖1 不同反應(yīng)pH值制得的Mo修飾的Pt/C催化劑的TEM圖Fig.1 TEM images of Mo-modified Pt/C catalysts synthesized at various pH values

圖2 不同反應(yīng)pH值制得的Mo修飾的Pt/C催化劑的粒徑分布圖Fig.2 Particle size distributions in the Mo-modified Pt/C catalysts synthesized at various pH values

圖3為不同反應(yīng)pH值制得的Mo修飾的Pt/C催化劑的XRD譜圖.可以看出,所有催化劑的XRD譜圖中在24.9°附近均有碳C(002)面的衍射峰,且都有明顯的Pt的衍射特征峰,但都沒有Mo的特征峰.這可能是由于催化劑中Mo元素以無定形型態(tài)存在或已嵌入到Pt的晶格中.26當(dāng)pH值為1時(shí),所制得的Mo修飾的Pt/C催化劑的Pt(111)的特征峰比較明顯,其他峰強(qiáng)度則較弱或不明顯.說明當(dāng)反應(yīng)溶液呈酸性時(shí),所制備的催化劑中Pt的結(jié)晶度較差.當(dāng)pH值為8、12或14時(shí),所制得的催化劑中面心立方晶格結(jié)構(gòu)的Pt的四個(gè)特征衍射峰Pt(111)、(200)、(220)和(311)均很明顯,說明堿性環(huán)境有利于制備結(jié)晶度完好的Pt顆粒.同時(shí),通過比較可以發(fā)現(xiàn),Mo修飾的Pt/C催化劑中Pt的晶面有所負(fù)移,說明摻雜的Mo嵌入到了Pt的晶格當(dāng)中并使其晶格參數(shù)增大.通常,Pt(111)、(200)和(311)晶面的衍射峰強(qiáng)度會(huì)受到載體碳材料特征峰的影響,而Pt(220)晶面特征峰不受載體干擾.27從圖中可以看出,隨著pH值的增加,Pt(220)晶面衍射峰的半峰寬依次增大,說明催化劑中活性組分的晶粒逐漸變小,這與TEM觀察到的結(jié)果是一致的.

圖3 不同反應(yīng)pH值制得的Mo修飾的Pt/C催化劑的XRD譜圖Fig.3 XRD patterns of Mo-modified Pt/C catalysts synthesized at various pH values pH:(a)1;(b)8;(c)12;(d)14

圖4為不同反應(yīng)pH值制得的Mo修飾的Pt/C催化劑中Pt元素4f軌道的XPS譜圖,其中71.0-71.5 eV(4f7/2軌道)和74.3-74.8 eV(4f5/2軌道)的峰對(duì)應(yīng)單質(zhì)金屬Pt,72.4-72.7 eV(4f7/2軌道)和75.5-75.9 eV(4f5/2軌道)對(duì)應(yīng)Pt(II),可能為PtO或Pt(OH)2.對(duì)XPS譜圖的擬合結(jié)果(見圖4)表明,各催化劑中的Pt主要以單質(zhì)金屬的形式存在,二價(jià)Pt只占很小的一部分,由此獲得的不同價(jià)態(tài)Pt的結(jié)合能列于表1.可以看出,隨著pH值的增加,無論是4f7/2還是4f5/2軌道,Pt(0)的結(jié)合能位置都逐漸正移,這可能是由于隨著反應(yīng)pH值的增加,所合成的催化劑中Mo元素和Pt元素之間的作用力逐漸增強(qiáng),28也進(jìn)一步證明了在所合成的催化劑中Mo已經(jīng)嵌入到Pt的晶格當(dāng)中.

圖5為不同反應(yīng)pH值制得的Mo修飾的Pt/C催化劑在0.5 mol·L-1H2SO4溶液中的循環(huán)伏安曲線,由此計(jì)算得到的各催化劑中Pt的有效電化學(xué)比表面積(ECSA)列于表2.可以看出,隨著反應(yīng)pH值的增大,催化劑的有效電化學(xué)比表面積明顯增加,pH值為14時(shí)達(dá)到最大值(41.91 m2·g-1).結(jié)合前文對(duì)TEM結(jié)果的討論,這可能是因?yàn)榇呋瘎┑念w粒隨pH值增加而減小導(dǎo)致的.

圖4 不同反應(yīng)pH值制得的Mo修飾的Pt/C催化劑中Pt的XPS譜圖Fig.4 XPS spectra of Pt in Mo-modified Pt/C catalysts synthesized at various pH values

表1 不同反應(yīng)pH值制得的Mo修飾的Pt/C催化劑中不同價(jià)態(tài)Pt的結(jié)合能Table 1 Binding energy of Pt species in Mo-modified Pt/C catalysts synthesized at various pH values

圖5 不同反應(yīng)pH值制得的Mo修飾的Pt/C催化劑在0.5 mol·L-1硫酸溶液中的循環(huán)伏安曲線Fig.5 CV curves of Mo-modified Pt/C catalysts synthesized at various pH values in 0.5 mol·L-1 H2SO4solution

表2 不同反應(yīng)pH值制得的Mo修飾的Pt/C催化劑的乙醇電氧化催化性能Table 2 Electro-catalytic performance of Mo-modified Pt/C catalysts synthesized at various pH values towards ethanol oxidation

圖6 不同反應(yīng)pH值制備的Mo修飾的Pt/C在0.5 mol·L-1硫酸+0.5 mol·L-1乙醇溶液中的循環(huán)伏安曲線Fig.6 CV curves of Mo-modified Pt/C catalysts synthesized at various pH values in 0.5 mol·L-1H2SO4+0.5 mol·L-1C2H5OH solution

圖6為不同反應(yīng)pH值制得Mo修飾的Pt/C催化劑在0.5 mol·L-1H2SO4+0.5 mol·L-1C2H5OH溶液中的循環(huán)伏安曲線,其中各催化劑對(duì)乙醇電催化氧化反應(yīng)的起始電位(onset電位,取反應(yīng)電流達(dá)到2 mA·cm-2時(shí)的電位)和峰值電流密度列于表2.可以看出,隨著反應(yīng)pH的提高,所得催化劑對(duì)乙醇電催化氧化反應(yīng)的起始電位逐漸降低、峰值電流密度逐漸增大,當(dāng)pH值提高到14時(shí),其乙醇氧化起始電位和峰值電流密度分別達(dá)到為0.503 V和32.79 mA·cm-2,這說明溶液pH值的提高有利于制備乙醇氧化催化活性更高的Mo修飾的Pt/C催化劑.

不同反應(yīng)pH值制得的Mo修飾的Pt/C催化劑在0.5 mol·L-1H2SO4+0.5 mol·L-1C2H5OH溶液中的計(jì)時(shí)電流曲線如圖7所示.在本文所研究的時(shí)間范圍內(nèi),各催化劑所表現(xiàn)出來的乙醇電催化氧化電流隨反應(yīng)pH值的增加而增大.在1000 s時(shí),反應(yīng)pH值為1、8、12、14所制得的催化劑的電流密度分別為0.94、3.16、2.06和7.07 mA·cm-2.這些結(jié)果說明反應(yīng)pH值為14時(shí)制得的催化劑對(duì)乙醇的電催化氧化不僅具有最高的催化活性,而且具有最好的穩(wěn)定性.

圖7 不同反應(yīng)pH值制得的Mo修飾的Pt/C催化劑在0.5 mol·L-1硫酸+0.5 mol·L-1乙醇溶液中的計(jì)時(shí)電流曲線Fig.7 Chronoamperometry curves of Mo-modified Pt/C catalysts synthesized at various pH values in 0.5 mol·L-1H2SO4+0.5 mol·L-1C2H5OH solution

圖8 不同反應(yīng)pH值制得的Mo修飾的Pt/C催化劑在0.5 mol·L-1硫酸+0.5 mol·L-1乙醇溶液中的電化學(xué)阻抗譜Fig.8 Electrochemical impedance spectroscopy of Mo-modified Pt/C catalysts synthesized at various pH values in 0.5 mol·L-1H2SO4+0.5 mol·L-1C2H5OH solution

圖8為不同反應(yīng)pH值制得的Mo修飾的Pt/C催化劑的電化學(xué)阻抗譜圖.可以看出,乙醇在反應(yīng)pH值為14時(shí)所制備的Mo修飾的Pt/C催化劑上的電化學(xué)反應(yīng)電阻最小,說明在該催化劑上乙醇的電氧化反應(yīng)具備最快的動(dòng)力學(xué)過程.當(dāng)pH值依次降低到12、8和1時(shí),在所制備的催化劑上,乙醇的電化學(xué)反應(yīng)電阻依次增大,且pH值為1時(shí)制得的Mo修飾的Pt/C催化劑上的電化學(xué)反應(yīng)電阻遠(yuǎn)大于其它條件下制備的催化劑,說明乙醇在Mo修飾的Pt/C催化劑上的電催化氧化反應(yīng)隨催化劑制備過程中反應(yīng)pH值的增加而變得更容易,催化劑的活性更高.

綜上所述,當(dāng)采用微波法制備Mo修飾的Pt/C催化劑時(shí),增加反應(yīng)的pH值有利于提高所制得催化劑的乙醇電氧化催化活性和穩(wěn)定性.這不僅與催化劑中Pt元素和Mo元素之間存在的協(xié)同作用相關(guān),28,29而且催化劑中活性組分Pt顆粒的尺寸以及分散狀態(tài)也是可能的原因之一,細(xì)小Pt顆粒在載體表面的均勻分布有利于催化劑性能的改善.

4 結(jié)論

反應(yīng)pH值對(duì)微波法合成的Mo修飾的Pt/C催化劑作為乙醇電氧化反應(yīng)催化劑的微結(jié)構(gòu)、形貌和催化性能有重要的影響.反應(yīng)pH值的提高有利于形成結(jié)晶良好的、小顆粒均勻分散的催化劑,從而改善催化劑的乙醇電氧化催化性能.在本文的研究范圍內(nèi)，當(dāng)反應(yīng)pH值為14時(shí),所制備的Mo修飾的Pt/C催化劑顆粒細(xì)小,且均勻分布在碳載體的表面,表現(xiàn)出最好的乙醇電氧化催化活性和穩(wěn)定性.

(1) Wee,J.H.Renew.Sust.Energ.Rev.2007,11,1720.doi:10.1016/j.rser.2006.01.005

(2) Peighambardoust,S.J.;Rowshanzamir,S.;Amjadi,M.Int.J.Hydrog.Energy2010,35,9349.doi:10.1016/j.ijhydene.2010.05.017

(3) Zhou,W.;Zhou,Z.;Song,S.;Li,W.;Sun,G.;Tsiakaras,P.;Xin,Q.Appl.Catal.B-Environ.2003,46,273.doi:10.1016/S0926-3373(03)00218-2

(4) Garcia,A.;Linares,J.;Chatenet,M.;Ticianelli,E.Electrocatalysis2014,5,41.doi:10.1007/s12678-013-0162-1

(5) Kamarudin,M.Z.F.;Kamarudin,S.K.;Masdar,M.S.;Daud,W.R.W.Int.J.Hydrog.Energy2013,38,9438.doi:10.1016/j.ijhydene.2012.07.059

(6) Song,S.;Tsiakaras,P.Appl.Catal.B-Environ2006,63,187.doi:10.1016/j.apcatb.2005.09.018

(7) Rightmire,R.A.;Rowland,R.L.;Boos,D.L.;Beals,D.L.J.Electrochem.Soc.1964,111,242.doi:10.1149/1.2426092

(8) Lai,S.C.S.;Koper,M.T.M.Phys.Chem.Chem.Phys.2009,11,10446.doi:10.1039/b913170a

(9) Gomes,J.F.;Bergamaski,K.;Pinto,M.F.S.;Miranda,P.B.J.Catal.2013,302,67.doi:10.1016/j.jcat.2013.02.024

(10) Antolini,E.J.Power Sources2007,170,1.doi:10.1016/j.jpowsour.2007.04.009

(11)Almeida,T.S.;Palma,L.M.;Morais,C.;Kokoh,K.B.;De Andrade,A.R.J.Electrochem.Soc.2013,160,F965.

(12) Linares,J.;Zignani,S.;Rocha,T.;Gonzalez,E.J.Appl.Electrochem.2013,43,147.doi:10.1007/s10800-012-0496-z

(13)Antolini,E.ChemSusChem2013,6,966.doi:10.1002/cssc.v6.6

(14) Hou,Z.;Yi,B.;Yu,H.;Lin,Z.;Zhang,H.J.Power Sources2003,123,116.doi:10.1016/S0378-7753(03)00515-9

(15) Tayal,J.;Rawat,B.;Basu,S.Int.J.Hydrog.Energy2011,36,14884.doi:10.1016/j.ijhydene.2011.03.035

(16) Santiago,E.I.;Camara,G.A.;Ticianelli,E.A.Electrochim.Acta2003,48,3527.doi:10.1016/S0013-4686(03)00474-2

(17) Ioroi,T.;Fujiwara,N.;Siroma,Z.;Yasuda,K.;Miyazaki,Y.Electrochem.Commun.2002,4,442.doi:10.1016/S1388-2481(02)00341-7

(18) Samjeské,G.;Wang,H.;L?ffler,T.;Baltruschat,H.Electrochim.Acta2002,47,3681.doi:10.1016/S0013-4686(02)00338-9

(19) Zeng,X.;Yuan,X.X.;Xia,X.Y.;Du,J.;Zhang,H.J.;Ma,Z.F.J.Inorg.Mater.2010,25,359.[曾 鑫,原鮮霞,夏曉蕓,杜 娟,張慧娟,馬紫峰.無機(jī)材料學(xué)報(bào),2010,25,359.]doi:10.3724/SP.J.1077.2010.00359

(20) Du,J.;Yuan,X.X.;Yu,J.H.;Ma,Z.F.Chin.J.Power Sources2007,31,873.[杜 娟,原鮮霞,余江虹,馬紫峰.電源技術(shù),2007,31,873.]

(21)Wu,X.;Jiang,Q.Z.;Ma,Z.F.;Shangguan,W.F.Chin.J.Inorg.Chem.2006,22,341.[吳 省,蔣淇忠,馬紫峰,上官文峰.無機(jī)化學(xué)學(xué)報(bào),2006,22,341.]

(22)Wang,X.Z.;Zheng,J.S.;Fu,R.;Ma,J.X.Acta Phys.-Chim.Sin.2011,27,85.[王喜照,鄭俊生,符 蓉,馬建新.物理化學(xué)學(xué)報(bào),2011,27,85.]doi:10.3866/PKU.WHXB20110111

(23)Bock,C.;Paquet,C.;Couillard,M.;Botton,G.A.;MacDougall,B.R.J.Am.Chem.Soc.2004,126,8028.doi:10.1021/ja0495819

(24)Chu,Y.Y.;Wang,Z.B.;Jiang,Z.Z.;Gu,D.M.;Yin,G.P.Fuel Cells2010,10,914.doi:10.1002/fuce.v10.6

(25) Zheng,Q.;Cheng,X.;Jao,T.C.;Weng,F.B.;Su,A.;Chiang,Y.C.Int.J.Hydrog.Energy2011,36,14599.doi:10.1016/j.ijhydene.2011.08.038

(26) Li,L.;Xu,B.Q.Acta Phys.-Chim.Sin.2005,21,1132.[李 莉,徐柏慶.物理化學(xué)學(xué)報(bào),2005,21,1132.]doi:10.3866/PKU.WHXB20051014

(27)Zhou,W.J.;Li,W,Z.;Zhou,Z.H.;Song,S.Q.;Wei,Z.B.;Sun,G.Q.;Tsiakaras,P.;Xin,Q.Chem.J.Chin.Univ.2003,24,858.[周衛(wèi)江,李文震,周振華,宋樹芹,魏昭彬,孫公權(quán),Tsiakaras,P.,辛 勤.高等學(xué)?；瘜W(xué)學(xué)報(bào),2003,24,858.]

(28) Yan,Z.;Wang,H.;Zhang,M.;Jiang,Z.;Jiang,T.;Xie,J.Electrochim.Acta2013,95,218.doi:10.1016/j.electacta.2013.02.031

(29)Cui,G.;Shen,P.K.;Meng,H.;Zhao,J.;Wu,G.J.Power Sources2011,196,6125.doi:10.1016/j.jpowsour.2011.03.042