多孔與平面硅電極界面的電化學行為

呂京美　程璇

(1北京理工大學珠海學院化工與材料學院，廣東珠海519088；2廈門大學材料學院材料科學與工程系，福建廈門361005；3福建省特種先進材料重點實驗室，福建廈門361005)

多孔與平面硅電極界面的電化學行為

呂京美1程璇2,3,*

(1北京理工大學珠海學院化工與材料學院，廣東珠海519088；2廈門大學材料學院材料科學與工程系，福建廈門361005；3福建省特種先進材料重點實驗室，福建廈門361005)

研究了重摻雜n-型單晶硅(CSi)在氫氟酸體系中生成多孔硅(PSi)的電化學行為，根據線性極化曲線，選取不同的電流密度，采用恒電流陽極極化法，制備了一系列多孔硅層。利用掃描電子顯微鏡對其進行了表面和斷面形貌的表征，通過線性掃描極化技術和計時電位法，比較了單晶硅電極和多孔硅電極的電化學行為，分析了多孔硅形成前后的塔菲爾曲線和計時電位曲線，給出了多孔硅形成過程中的重要電化學參數，如腐蝕電流、開路電位、塔菲爾斜率等。并對其進行深入分析，根據實驗結果，提出了單晶硅電極/電解質界面和多孔硅電極/電解質界面的結構模型，并利用該模型討論了兩種電極界面的電化學特性。

多孔硅；硅電極界面；腐蝕速率；開路電位；極化曲線；界面電化學

Based on our former work11,12,this paper w as carried out to systemically investigate theelectrochemicalbehaviorsof CSiand PSi formedwith differentm icrostructures by anodic polarization in HFand ethanolsolutions.The potential－currentdensity(E－j) and potential－time(E－t)curvesweremeasured during the PSi formation.Morphology and surface species analysis,the electrochem ical responses of CSi electrode/electrolyte and PSi electrode/electrolyte interfaces are compared,and further discussed based on the interfacialmodelsproposed.

2 Experimental

The electrochemicalmeasurementswere done using a flat cell made of polytetrafluoroethy lene(PTFE)w ith a three-electrode system in 40%HF:ethanol(1:2 by volume)solutions on flat n(111)CSi(~25mΩ?cm)as described previously12.Theworking, counter,and reference electrodeswere CSi,Pt foil,and saturated calomel electrode(SCE),respectively.The AutoLab electrochem icalworking station(Model PGSTAT 30from Metrohm AG, Switzerland)wasused to record theelectrochemicaldata.The PSi layerswere fabricated by constantly applying different current densities(jA)based on the polarization curves for 10m in.The polarization curvewasmeasured from－0.335to 5.665V,while Tafel curvesin thepotential rangeof±250mV versus open circuit potential(EOCP)before(CSielectrode)and after the PSi formation (PSielectrode).Themorphologiesof the PSi layerswereobserved by using the field em ission scanning electronm icroscope(SEM) (LEO 1530,Oxford instrument,Germany)w ithout surface coating.A ll the reagentsused in ourworkweremostly analytical reagent,except the HFsolutionwhichwasguarantee reagent.

Fig.1 Typicalpolarization curveof highly doped n(111)-Si(CSi) in 40%HF:ethanol(1:2 by volume)solutionsat50m V?s－1The insetshows thepolarization behaviorofCSiat thepotentials<1.0V.

3 Results and discussion

A typical polarization curve of heavily doped n(111)-Si(referred to CSi)electrode in HF-ethanol solutions is presented in Fig.1.It is evident from Fig.1 that the current density initially increased very slow ly w ith the increase of potential,and then increased rapidly at the potentials>0.5V.Minor current oscillationswereobserved at the potentials>3.0V due perhaps to the hydrogen bubbles produced during the anodic polarization.No characteristic peak corresponding to PSi formation was found in Fig.1.

To further investigatew hether the PSi layerswere formed on the CSi surface during the anodic polarization,four current densities,1,10,20,and 30mA?cm－2,were selected and applied to CSielectrodesw ith the etching time being fixed at10m in.The responded potentials(E)weremeasured asa function of time(t) and the E－t curves are given in Fig.2(a).The responding potentialsdecreased very sharply at the beginning of currentdensity applications,which is usually caused by the double layer discharge,and then reached a pseudo-steady state.However,the time to reach a pseudo-steady state(tps)variedw ith themagnitude of current density applied(jA).It required nearly 2 s to attain stable potential values at30mA?cm－2,while approximately 50s at1 mA?cm－2.A shorter tpsata higher jAgenerally indicatesa faster electrochem ical reaction at an electrode/electrolyte interface.A plotof the stabilized potentials(ES)against jAisshown in Fig.2(b). A linear increase in ESwith an increaseof jAwas clearly observed, show ing the rate-determ inate step being hole concentration in space charge layer(SCL)instead of electrolyte concentration.A larger jAvalue im pliesa higher hole concentration in n-type CSi, accordingly,causes a larger redox potential(ES).In addition, apparentoscillationswere observed at jA≥10mA?cm－2,which m ightbe caused by the formation and evolution of the hydrogen gasoccurred during the polarization process.

The surfacemorphologies of the CSi sam ples obtained after different jAapplications are presented in Fig.3.A sponge-like structurew ith pore diameters ranging 10－20nm was observed and the porosity increased slightly w ith an increase of jAwhich led tomoreuniformed pore distributions.The thicknessof a PSi layer (LPSi)could be roughly determined by the cross-sectional SEMimages(asshow in Fig.4(a,b)),and isplotted against jAin Fig.4(c).Itseemed that the LPSiincreased linearlywith jA,implying thatthe grow th rate of PSi layer strongly depended on themagnitude of jA.

Fig.2 Chronoam peromogram s recorded during the PSi formation(a)and variation of stabilized potentialw ith d ifferent app lied cu rrent densities(b)

In order to obtain kinetic information of charge-transfer reactions occurred on anodic and catholic sitesw ithout introducing any remarked changes into the Si surface,Tafel curves were measured with the Sielectrode before and after theapplications of jAin HF-ethanol solution in the range of±250mV(vs EOCP). The results are shown in Fig.5(a).It is obvious that the Tafel curvesobtainedw ith Sielectrode before jA(representing CSi)were well reproduciblewith the identical EOCPof－0.3V(dashed lines in Fig.5(a)),certified a relatively stable and repeatableelectrode surface state of CSi.How ever,the Tafel curvesobtained after jA(representing PSi)shifted tomore negative potentialsand larger currentdensitiesas jAincreased(solid lines in Fig.5(a)),suggesting that the electrochem ical reactions took place at the electrode surface,leading to PSiformation.

Fig.3 Sur facem orphologies of PSim icrostructures form ed on n(111)-Siw ith the app lied cur rent densitiesof (a)1m A?cm－2,(b)10m A?cm－2,(c)20m A?cm－2,and(d)30m A?cm－2for 10min

Fig.4Cross-sectionalSEMimagesof PSim icrostructures formed on n(111)-Siw ith theapplied currentdensitiesof(a)10m A?cm－2, (b)20m A?cm－2;(c)theobserved th icknessof PSi layers LPSifabricated atdifferent jAvalues for 10m in asa function ofapp lied currentdensity

When the electrode/electrolyte interface is controlled by the activation,the kinetics can be described by the Bulter-Volmer equation.When theoverpotential(η)is largeenough,the Bulter-Volmerequation can be simplified to

where a and b are Tafel constant and Tafel slope,respectively,and can be calculated by the following equations for anodic polarization.

where R is thegas constant,T is the temperature,α?is the transfer coefficient formultiple electron electrode reaction,which was usually used to estimate the influence of the potential for the electrode reaction activation energy,F is the Faraday constant,and jcorris the corrosion current density.The Tafel slopes includingcathode branch(bc)and anodic branch(ba)were determ ined by fitting the linear regions in Tafel plot,as illustrated in Fig.5(b). The jcorrcan be expressed in terms of the corrosion rate of Siaccording to

where d N/d t represents the changesof Siatom numbers in aunit area and a unit time(atom?cm－2?s－1),and NAis the Avogadro′s constant(6.022×1023mol－1).The important kinetic parameters obtained are summarized in Table1.Evidently,the values of bcand bafor both CSiand PSielectrodeswere approximately the same(ba≈bc)upon each jA,show ing theexcellentsymmetry in the anodic and catholic branches.The relative bigger valuesof bcand ba(0.17－0.38V?dec－1)observed both for CSiand PSielectrodes indicted that the corrosion rateof Siwasnotcontrolled by the rate of charge transfer in SCL alone,which had a typicalslop of 0.06V?dec－1for a heavily doped Si13.Thus,the steeper slopswere caused by thepotentialdrop in Helmholtz layer.The typicalvalues of theα?averaged 0.30for CSi,however,reduced apparently to 0.16for PSi.The Ecorrbecamemore negative,w hile jcorrlowered w ith the increasesof jA.

Fig.5(a)Polarization curvesobtained before and after the PSiformation;(b)a Tafelplot for typical n(111)-SielectrodeThe Tafelslopes including cathodebranch(bc)and anodic branch(ba)(insets)weredeterm inedby fitting the linear regions in Tafelplot,as illustrated in Fig.5(b).

Table1 Fitting parametersof Tafel curvesduring the PSi formation

The differences in the kinetic parametersobtained w ith CSiand PSielectrodessuggested theapparentalterationsat theelectrode/ electrolyte interface.The abridgedmodel for n-type CSiand PSi electrode/electrolyte interface is proposed and schematically illustrated in Fig.6.As can be seen in Fig.6(a),fora CSielectrode, when equilibrium between electrodeand electrolyte is reached,the SCL of CSielectrode is in the state of electron depletion w ith few holes left,thus,negative ionsare adsorbed on the interfacewhich results in negative EOCP.While fora PSielectrode,thenano-sized particles observed in Fig.2 result in a higher potential barrier, holesmove to the pore bottom but cannot reach the pore wall according to the quantum confinementmodel14,which would change the charge carriers distribution both in SCL and Helmholtz layeras depicted in Fig.6(b).Comparingwith the CSielectrode, the formation of PSi layer isequivalent to the extended thickness of Helmholtz layer,the interfacial thickness increased with larger jA,in particular,from 3.2μm at5mA?cm－2to27.4μm at30mA?cm－2,12which would result inmore electrons adsorbed on the interface and thusmore negative EOCP(Fig.6(c)).The potential drop at the Helmholtz layerwould be increased when the nanostructures formedw ithin the PSi layersby changing jAfrom 1 to 30mA?cm－2.

In order to quantitatively characterize the difference between the PSiand CSielectrodes,the follow ing parametersare defined:

where N is the numberof Siatoms,t is the time,and EOCPis the open circuit potential.Themagnitude ofΔ(d N/d t)indicates the corrosion rate of Si,whileΔEOCPthe degree of potential shiftafterthe PSi formation.The calculatedΔ(d N/d t)andΔEOCPare plotted against ESin Fig.7.The differences betw een CSi and PSimay reflect indirectly the different structures of electrode and the variations of surface state.Evidently,Δ(d N/d t)increased linearly w ith ES,meaning that the electrochem ical reaction rate ata PSi/ electrolyte interface is faster than thatataCSi/electrolyte interface and closely related to the linearly increased thicknessof PSi layer (Fig.4)because of the larger active surface areas resulted.Furthermore,the amounts of hydrogenous surface species were thought to accelerate the corrosion rate of PSi in solution by inducing lotsof surface state,which can be taken as theelectronic energy level for the exchange of charge carriers.During the PSi formation by jA,all Si―H and Si―O species including Si3SiH, Si2SiH2,Si(O)SiH3,O2SiH2,O3SiH were identified from the FTIR spectra12.The absorption area under the FTIR spectrum(S)can be obtained by integration of the curve.Therefore,the S values represent sem i-quantitatively the amounts of those species.As evident in Fig.7,similar to the corrosion rateΔ(d N/d t),the S also increased linearly w ith ES.On the other hand,theΔEOCPvalues becamemore negative and strongly dependentupon ES,which can be fittedwellwith the Boltzmann function shown below.

The negative shift in EOCPfor n-type CSi is caused by the electrochem ical reaction at CSielectrode in HF solution which depletes holes and leads to accumulation of negative charges,as illustrated in Fig.6(b).The depletion of holes in the pore wall resulted in large amounts of negative charges which w ould be increased with the thicker PSi layer.However for the PSi fabricated at high jA,the accumulation of negative charges at the interface became slow er because of the linear increase of the reaction ratew ith ES.Themost reasonable interpretation is faster consumption of holesundera larger jAwhich lessensaccumulation of thenegative chargeand less influencesΔEOCP.

Fig.6Proposedmodel for n-type Sibefore(a)and after(b)the PSi formation at theelectrode/electrolyte inter face w ith a schem atic sketch of reaction ratesat the in ter face(c)tand for the exchange rates for holesand electrons,respectively.

Fig.7 Variations of param etersw ith responding poten tials determ ined from Fig.2(b)

4Conclusions

The PSi layerswith nanosized poreswere formed in a broaden rangebetween 1 and 30mA?cm－2on heavily doped n(111)-Si.The reaction rate occurred at the Si/electrolyte interface was determ ined by the concentration of holes.The thickness of PSi layer, the difference in reaction rate between CSiand PSi,Δ(d N/d t),as wellas the totalamountsof Si―H and Si―O species increased linearly w ith jA.The variation inΔEOCPfollow ed the Boltzmann equation.Compared to the flatCSielectrode,the formation of PSi layers increased the thicknessof Helmholtz layerand changed the potential distribution,accordingly,the electrochemical reaction at the PSi/electrolyte interface slightly changed.

References

(1)Zhang,X.G.Electrochemistry ofSilicon and ItsOxide; KluwerAcademic/Plenum:New York,2001;pp 5－98.

(2)Lemann,V.Electrochemistry ofSilicon:Instrumentation, Science,Materialsand Applications;Wiley-VCH:New York, 2002;pp 10－103.

(3)Yang,H.B.;Hu,M.;Liang,J.R.;Zhang,X.R.;Liu,Z.G.Acta Phys.-Chim.Sin.2008,24,1017.[楊海波,胡明,梁繼然,張緒瑞,劉志剛.物理化學學報,2008,24,1017.]doi: 10.3866/PKU.WHXB20080618

(4)Tian,J.;Zheng,D.;Zhang,X.G.;Zhang,B.H.;Xia,B.J.; Yang,H.Acta Phys.-Chim.Sin.2007,23,68.[田娟,鄭丹,張熙貴,張寶宏,夏保佳,楊輝.物理化學學報,2007, 23,68.]doi:10.3866/PKU.WHXB20070114

(5)Parkhutik,V.P.;Glinenko,L.K.;Labunov,V.A.Surf.Technol. 1983,20,265.doi:10.1016/0376-4583(83)90009-2

(6)Sailor,M.J.Fundamentalsof PorousSilicon Preparation.In Porous Silicon in Practice;W iley-VCH Verlag GmbH&Co. KGaA:Weinheim,Germany,2011;pp 1－42.

(7)Bertagna,V.;Plougonven,C.;Rouelle,F.;Chem la,M. J.Electroanal.Chem.1997,422,115.doi:10.1016/S0022-0728(96)04902-9

(8)deMierry,P.;Etcheberry,A.;Rizk,R.;Etchegoin,P.; Aucouturier,M.J.Electrochem.Soc.1994,141,1539.doi: 10.1149/1.2054959

(9)A llongue,P.;Henry de Villeneuve,C.;Pinsard,L.;Bernard,M. C.Appl.Phys.Lett.1995,67,941.doi:10.1063/1.114702

(10)Parkhutik,V.;Andrade Ibarra,E.Mat.Sci.Eng.B 1999,58, 95.doi:10.1016/S0921-5107(98)00297-9

(11)Lu,J.M.;Cheng,X.Adv.Mater.Res.2013,650,34.doi: 10.4028/www.scientific.net/AMR.650

(12)Lu,J.M.;Cheng,X.;Zhang,Y.;Zhu,X.F.Electrochim.Acta 2010,55,5084.doi:10.1016/j.electacta.2010.03.100

(13)Zhang,X.G.;Collins,S.D.;Smith,R.L.J.Electrochem.Soc. 1989,136,1561.doi:10.1149/1.2096961

(14)Lehmann,V.;G?sele,U.Adv.Mater.1992,4,114.doi: 10.1002/adma.19920040212

Electrochemical Behavior of Porous and Flat Silicon Electrode Interfaces

LüJing-Mei1CHENG Xuan2,3,*
(1College ofChemical Engineering and Materials Science,Beijing Institute ofTechnology,Zhuhai519088,Guangdong Province, P.R.China;2DepartmentofMaterials Science and Engineering,Co llege ofMaterials,Xiamen University,Xiamen 361005, Fujian Province,P.R.China;3Fujian Key Laboratory ofAdvanced Materials,Xiamen 361005,Fujian Province,P.R.China)

The e lectrochem ical responses of heavily doped n-type single-crysta lsilicon(CSi)during the formation ofporous silicon(PSi)layers in hydrofluoric acid-based electrolyteswere investigated.A series of PSi layers were fabricated by constantly app lying differentanodic currentdensities,which w ere selected according to the linear polarization curve.The surface and cross-sectionalmorphologies of the PSilayerswere studied by scanning electronm icroscopy.The electrochem ica lbehaviorofCSiand PSie lectrodeswas compared by linear sweep po larization and chronopotentiometry techniques.The importantparameters associated with the electrochem ical reactions at the electrodeswere evaluated by analyzing the Tafelplots and chronopotentiograms obtained before and after the PSi formation.Structuralmodels of the CSielectrode/electrolyte and PSielectrode/ electrolyte interfaceswere suggested based on the experimentaldata.Accordingly,the interfacialcharacteristics of CSiand PSie lectrodeswere discussed.

Porous silicon;Silicon electrode interface;Corrosion rate;Open circuitpotential; Pola riza tion curve;Inte rface e lectrochem istry

1 Introduction

Poroussilicon(PSi)w ith differentm icrostructures can beobtained from flat signal crystal silicon(CSi)electrode in hydrofluoric acid(HF)contained solutions w ith the electrochem icalreactions taking placeat the interfaceof Sielectrode/electrolyte1－4. For a typical sem iconductor,an interface of electrode and electrolyte includes space charge layer in solid electrode side and Helmholtz layer in solution side.However,the distribution and transfer of charge carriers(electrons and holes)at the interface could become complex and different typesof potential-time curves were obtained during the PSi formation5,6.The kinetics of electrochemical corrosion of CSielectrode in dilute HFsolutionshave been previously studied by exam ining open circuit potential, polarization resistance,corrosion current density,and reaction rate7.However,tillnow,mostworks focused on the interface of flat Si electrode/electrolyte.Few have been reported at the PSi electrode/electrolyte interfacewhich should be important for the application and the formation mechanism of the PSi.It was pointed out that PSielectrode possesses larger surfacearea leading tomoreadsorbed specieswhich would change the distribution of carriers such as silicon-hydrogen species8－10.Meanwhile,the formation of nanopore structure would change the charge distribution in space charge layerof Sielectrode.

October9,2015;Revised:December22,2015;Published onWeb:December24,2015.

O646

10.3866/PKU.WHXB201512242

*Corresponding author.Email:xcheng@xmu.edu.cn;Tel:+86-592-2185599;Fax:+86-592-2183937.

The projectwas supported by the National Natural Science Foundation of China(11372263).

國家自然科學基金(11372263)資助項目?Editorialofficeof Acta Physico-Chim ica Sinica