您的当前位置:首页正文

锂离子电池硅薄膜电极充电膨胀的有限元仿真及其实验验证

2021-11-17 来源:品趣旅游知识分享网
锂离子电池硅薄膜电极充电膨胀的有限元仿真及其实验验证

季家磊;朱孔军;刘鹏程;钱国明;王裕;刘劲松

【摘 要】硅基负极锂离子电池材料因其具有高的理论容量(4200m Ah/g)而成为最有希望的高容量负极材料之一.但硅负极在充放电过程中的体积效应,将引起电极材料粉化以及循环性能变差.为解决上述问题,将硅与惰性过渡金属材料复合,过渡金属充当体积效应的缓冲层.本文利用有限元软件abaqus对比了三种不同的硅薄膜材料(Si/Si-M n/Si-Zr).通过磁控溅射方法制备了上述三种硅薄膜材料,并对其进行了SEM、XRD、循环性能等测试,实验得出的结论与仿真结果一致,加入的过渡金属材料有利于缓解体积效应,且Mn材料的缓解效应更强.%Si is the most promising anode material for high energy lithium ion battteries because of its high specific capacity(4200mAh·g -1).But its undesirable volume enpansion results in me-chanical degration and capacity reduction.It is a promising way to combine Si and inert metal to relieve the expansion during Li+insertion/extraction.In this article,use Abaqus to compare three different Si thin films(Si,Si-Mn,Si-Zr).Si thin film was deposited on Cu foil by magne-tron supttering for use as lithium ion battery anode

material.The electrochemical performance of Si film was investgated by cyclic voltammetry and constant current charge/discharge test.The results consistent with simulation.The use of metal material is useful for the electronical per-formance and Zr is more useful than M n. 【期刊名称】《电池工业》 【年(卷),期】2017(021)006

【总页数】9页(P19-27)

【关键词】锂离子电池;硅薄膜负极;惰性金属;Abaqus;磁控溅射 【作 者】季家磊;朱孔军;刘鹏程;钱国明;王裕;刘劲松

【作者单位】南京航空航天大学航空宇航学院,江苏南京 210016;南京航空航天大学航空宇航学院,江苏南京 210016;广州大学机械与电子工程学院,广东广州 510006;南京航空航天大学航空宇航学院,江苏南京 210016;南京航空航天大学材料科学与工程学院,江苏南京 210016;南京航空航天大学材料科学与工程学院,江苏南京 210016 【正文语种】中 文 【中图分类】TM910.2

1 Introduction

Due to their advantages of high energy density, long life, low toxicity and environmental friendliness, lithium-ion batteries(LIBs) have become the most promising and widely applied rechargeable batteries.[1] LIBs have been widely used in portable electronics such as mobile phone, digital camera, DV, laptop, and (hybrid) electrical vehicles. The theoretical capacity of commercial graphite (used as anode) is only 372mAh/g[2], and can not meet the increasing demands for lithium-ion batteries with high energy density and long cycling life. In recent years, the development of new high capacity anode material has attracted significant interest. It is well known that some elements can electrochemical react with Li with high capacity.

Some alloying elements with high theoretical capacities, such as Si, Sn, Ge Al[3-6], and conversion electrodes such as NiO, and Co3O4[7-8],have been studied extensively. Among these material, Si has high theoretical capacity,4100mAh/g, ten times of graphite[9]. However, Si shows a massive volume expansion/contraction during Li+ insertion/extraction, larger than 300% after fullly lithium insertion[10]. This causes the pulverization of Si particle and loose contacts between Si particles and current collector, which will further result in mechanical in mechanical instability and poor cyclability[11-13]. To solve such problems, combine Si and inert metal materials which can relieve the huge volume change of Si thin films during lithiation and delithiation. Researchers have made attempts to improve the electrochemical performance of Si thin films as anode material, among which, the introduction of a secondary material is an effective way[14-16].

In this study, choose secondary materials which have good conductivity and ductility and act as buffer to alleviate the particle pulverization. Use Abaqus to compare the stree and strain in three different thin films(Si, Si-Mn, Si-Zr) during Li+ insertion/extraction and analyse the role of the inert metal. Then, fabricate above Si thin films by magnetron supttering. The electrochemical performance of Si film was investgated by cyclic

voltammetry and constant current charge/discharge test. The experimental results consist with the simulation. The use of metal material is useful for the cycling performance and Mn is more useful. 2 Finite Element Model

Li+ insertion will result in a distorted lattice, volumetric expansion, mechanical stresss occures because of the constraint of Cu substrate. The size and stiffness of the substrate(Cu foil) is much lower than Si thin film, the deformation of Cu foil is then much lower than Si thin film and can be neglected, we assume the substrate to be rigid. Cracking and interface debonding are not considered, body force and inertia effects are neglected. Mechaniacl deformaion is thought to be quasi-static because it is much slower than diffusion process.

An axisymetric finite element model under a cylindrical polar coordinate system(r,θ,z) is used in Abaqus. Si thin film is assumed to be homogeneous and isotropic and be firmly bonded to the rigid substate. Because

mechanical stress and diffusion process influence each other, fully coupled thermal-mechanical transient analysis procedure is used. First-order elements are used for the highly nonliner problem, finite element size is set to 1% of the height of Si thin film and fine mesh is used due to stress concentration. To improve convergence of the nonliner problem, liner search algorithm and maximum 5 interations are used.

There is no diffusion-stress aanalysis in Abaqus,use the method proposed by Prussin[17] as convention. Mechanical response under concetration loading is analogous to that under temperature loading, stress caused by diffusion is analogous to thermal stress.

Extending the 1D relation given by Prussin[17] to 3D, the constitutive equation for diffudsion-induced deformation of an elastic solid can be expressed

(1)

Fig.1 Structure of thin film

where εij(i,j=1,2,3) are componts of strain tensor; σij(i,j=1,2,3) are componts of stress tensor; c(mol m3) is concentration of diffusion componts; Ω is partial molar volume representing volume expansion caused by diffusion of Li+; E is elastic modulus; υ is posson’s ratio. Stress caused by diffusion is analogous to that caused by temperature gradient, Ω/3 plays the same role as thermal expansion coefficient in thermal stress analysis.

2.1 Structure and Material

The model of anode is based on the 2016-type cell which is used to be tested later. The anode of 2016-type cell is wafer thin film The thickness and radius of the Si thin film is D and R, the thickness of transition metal material is d. According to 2016-type cell, R is set to be 6μm. Bourderau S[18] fabricated the Si thin film with 1.2μm which had bad cycling performance, while thin film with 275nm[19] had better cycling performance, thus D is set to be 500nm. Transition metal just works as buffer layer and not participate in LI+ insertion/extraction, d is smaller than Si and set to be 200nm(Fig. 1.).

Based on volumes of lithiated silicon at different Li-Si alloy phases[20], and linear relations between Li fraction and elastic constants[21- 22], dependence of elastic constants on concentration c (fmol mm-3) is expressed.

E=E0+k1 c,υ=υ0+k2 c. (2)

Where E0=130Gpa, V0=0.22[23]. k1= -0.13Gpa.μm3fmol-1,k2=-0.00047μm3fmol-1(minus k1、k2 represents the soften of Si electrode during lithium intercalation.)

The choice of transition metal must have good ductility, it acts as the buffer to alleviate the huge expansion, at the same time, it doesn’t act with Li+. Metal choosed here is Mn and Zr. 2.2 Boundary Condition

As metioned previous,the structure of the electrode is wafer type and symmetry, also the electrode is surrounded by invariant Li-ion

concentration,the electrod can be treated as a symmetric finite model, and for simplify, we choose a section for analysis. The initial condition is when t=0, c=0. (3)

In potentiostatic operation,the electrode surfaces are surrounded by an invariant Li-ion concentration, cs, so the concentration of Li-ion on the top surface and edge surface is fixed. when 0when 0Cu foil is rigid substrate and doesn’t take part in Li+ diffusion,

(6)

Under the cylindrical polar coordinate system, the structure, boundary conditions and loading conditions are all axisymmetric. when t>0,r=0, ur=0. (7)

Volume change consistsin stress because of the constraint of substrate. Because the film is firmly adherent to the substrate, there is no lateral displacements occures on the interface. when t>0,z=0,ur=uθ=uz=0. (8)

There is no mechanical loading applied on the top surface and side sur face.

when t>0,在z=h处,σz=0. (9)

when t>0, r=R, z>0, σr=0. (10)

3 Simulation Results

3.1 Concentratin, Displacement and Stress Fields

The concentration field before fully insertion is showed in Fig. 2a. Due to the edge diffusion, concention is dependent on radial coordinate. For the central region of the electrode, concentration is dependent on axial displacement.

The displacement and stress field after fully insertion is showed in Fig. 2b-e.

Fig. 2d-e. shows the expansion caused by lithium-ion insertion includes radial extension and bending.

The radial displacement is concentrated at the edge of the top surface and the maximum radial displacement occures at the edge on the top surface, also there is little radial displacement in the central region of the film. The maximum axial displacement occures at the center of the top surface. Axial displacement in the central region can be regarded to be

independent of radial coordinate. Due to the fixed constraint of the rigid substrate, negative axial displacement is possible near the edge on the interface. A dome-like morphology is formed due to the axial and lateral expansions.

Fig.2 Concentration, displacemt and stress fields (a.concentration field, b.stress field, c, d, e. displacement field in equilibrium state) Fig.2b shows the stress caused by lithium-ion insertion mainly occures at the center of the top surface and the edge of the bottom surface.

3.2 ComparasionofDisplacement/StressFieldsinDifferentSi-MThinFilms The displacement after fully insertion is showed in Fig.3. The distribution of displacement in different Si-M electrode is similar. No matter the total displacement or vertional/radial displacement declines while the metal is used.

Same conclusion can be achieved in the stress fields (Fig. 4). The maximun of von mises of Si-M thin film is less than Si thin film. Both the results of displacement and stress comparasion fields reveals that use of metal is beneficial to the Si anode to experience less destroy during insertion.

Fig.3 Comparasion of displacement fields (a.total, b.radial, c.axial displacement)

Fig. 4 Comparasion of stress fields(a.Si, b.Si-Mn,c.Si-Zr) 4 Experimental Results 4.1 Experimental

Si thin films were prepared in an PVD 75 multi-target magnetron

sputtering system(KJLC, Co.). The samples were deposited on both Si wafer for thickness measurement and Cu foil for electrochemical measurements. The target was N-type monocrystalline Si with 2 inch diameters and 99.999% purity, Mn with 99.9% purity and Zr with 99.5% purity. The target-substrate distance of the sputtering system was set to be 50mm. After the base pressure reached 8.3×10-4 Pa, Ar (99.999%) was introduced into the chamber. The working pressure was kept at 8mtorr. Si thin films were deposited using a constant radio frequency power supply of 100W.The film thickness was controlled by deposition time.The amount of deposited Si was calculated assuming a density of 2.33g·cm-3 for the Si thin film. The morphology and accurate thickness of Si thin films were measured by the field emission scanning electron microscopy (FSEM, SIGMA, Germany). The phase structure of was analyzed by X-ray diffraction (XRD, Bruker D8 Advance, Germany).

To evaluate the electrochemical properties of the Si thin film anode, 2025-type half-cells were assembled in an argon-filled glove box with H2O and O2 concentrations of less than 1ppm. A lithium metal foil was used as a counter electrode, and Celgard2400 was used as a separator. The

electrolyte solution was 1.0 M LiPF6 in EC/DEC (1∶1 vol/vol). Cyclic voltammetry measurements were performed using an electrochemical workstation (Princeton PARSTAT MC) at a scan rate of 0.01 mV in the potential range 0V~1.5V.Galvanostatic charge/discharge measurement was carried out using a Land battery test system (LAND CT2001A) with the cut off potentials being 0V versus Li/Li+ for discharge and 1.5V versus Li/Li+ for charge.

4.2 Results and Discussion

Fig.5 SEM images(a.cross-ssection, subface of b.Si, c. Si-Mn, d. Si-Zr) The cross-sectional SEM image of Si thin film deposited on a Si wafer is presented in Fig. 5a. The thickness of the dense Si、Mn、Zr can be observed,and the corresponding growth rate can be calculated,finally actual operating time was obtained accoring to the target thickness(Tab 1). Table1

SputteringParameterTargetSiMnZrPower(W)100100100Time(min)606030Thickness(nm)320200170Rate(nm/min)5.33.35.7TargetThickness(nm)500200200ActualTime(min)946035

Fig.6. shows the XRD pattern of Si thin film deposited on Cu foil. All the diffraction peaks are attributed to the Cu foil, and no peak of Si appears, especially the typical peak for crystal Si at 28°. This indicates that the Si thin film is amorphous.

Fig.6 XRD patterns of Si thin films

The L+ insertion/extraction reactions of Si thin film were studied by cyclic voltammetry. For all of the three thin films, three cyclic voltammetric

curves of the Si thin film are shown in Fig. 7. In the first scanning cycle, there is a cathodic peak at 0.32V, which disappears from the second cycle. This cathodic peak is attributed to the formation of a solid electrolyte interphase (SEI) layer due to decomposition of electrolyte on the film surface. Two cathodic peaks at 0.20V and 0.05V, as well as two anodic peaks at 0.50V and 0.33V, are observed on all three cyclic voltammograms; these are ascribed to the electrochemical reactions of Li+ insertion and extraction in the Si thin film. The slight difference in the intensity and the potential for each peak can be attributed to the kinetic effect involved in the cyclic voltammetry measurement.

Fig.7 Cyclic voltammetry plots (scanning rate 0.1Mv/s,potential range 0V~1.5V, a. Si; b.Si-Mn; c.Si-Zr)

Fig.8.shows the first three times of the discharge/charge curves. The first discharge capacity of the Si, Si-Mn, Si-Zr thin film is 2045.0mAh·g-1, 2203.1mAh·g-1, 2505.0mAh·g-1, and initial coulombic efficiency is 101.76%, 103.98%, 102.89%. The first and second reversible capacity of the Si-Mn thin film is 1900.3mAh·g-1 and 1976.0mAh g-1,for Si-Zr, 1997.0mAh·g-1, 2054.7mAh·g-1,which is much larger than that of a graphite

anode(1662.9mAh·g-1, 1692.3mAh· g-1,respectively). The irreversible capacity is attributed to the formation of a SEI layer in the first cycle. In evidence, a SEI-formation voltage plateau is observed near 0.32V, which disappears in the second cycle. This observation is also in good agreement with CV results.

Fig.8 Discharge/charge curves (a. Si; b.Si-Mn; c.Si-Zr)

Cycling performance of the Si thin films are shown in Fig. 9 a-c. The first reversible capacity 100mA/g for Si, Si/Mn, Si/Zr is 1692.3mAh/g, 1830.8mAh/g, 1955.6mAh/g respectively, and 71.2%, 83.9%, 88.2% capacity remained after 50 cycles.The introduce of transition metal can enhance both the first reversible capacity and the capacity retention effectively, which proves that the metal can improve the cycling performance of Si thin films. Rate performance of the Si thin films are shown in Fig. 9 d-f.

Fig.9 Eletronical performance of Si thin films (Cycling performance of a.Si, b. Si-Mn, c. Si-Zr, rate performance of d. Si, e. Si-Mn, f. Si-Zr)

To further evaluate the performance of Si thin films, the rate capability measurements (Fig. 9d-f) at the quickly increased current density from 0.1A/g to 1A/g were carried out. For Si thin films, the discharge capacity of 2045.3mAh/g, 1413.2mAh/g, 1128.8mAh/g, 919.5mAh/g, 732.7mAh/g can be obtained at 0.1A/g, 0.2A/g, 0.3A/g, 0.5A/g, 1.0A/g, 0.1A/g. For Si-Mn, 2203.1mAh/g, 1718.9mAh/g, 1535.7mAh/g, 1329.6mAh/g, 1044.3mAh/g can be obtained, and for Si/Zr, 2505.0mAh/g, 1859.4mAh/g, 1661.2mAh/g, 1500.6mAh/g, 1117.8mAh/g can be delievered. Although suffering from the rapid change of the current density, the cell can still exhibit a stable cycling at each current. Importantly, 80% of the first reversible capacity can be remained for Si thin film when the current density is turned back to 1A/g, and for Si-Mn, Si-Zr, 89.0% and 92.9% can be remained.It proved that use of metal is beneficial to the electrocal performance again. 5 Conclusion

In this article, Si and inert metal is combined to relieve the expansion during Li+ insertion/ extraction. Use Abaqus to compare three different Si thin films (Si, Si-Mn, Si-Zr).we found that the use of inert metal reduces the displacement and stress induced during the Li+ insertion. Also, Si-M thin film used as anode material.was deposited by magnetron supttering The morphology of the Si-M thin films are similar, and XRD results reveals that the structure of Si thin films is amorphous. The electrochemical

performance of Si thin films consistents with the simulation, use of metal can relieve the expansion and result in better cycling and rate performance. Among Mn and Zr, Mn is more useful. References:

【相关文献】

[1] Goodenough J B, Park K S, The Li-ion Rechargeable Battery: A Perspective[J],American Chemical Society, Journal,2013. 135(4):1167-76.

[2] Wachtler M, Besenhard J O, Winter M, Tin and tin-based intermetallics as new anode materials for lithium-ion cells[J], Journal of Power Sources, 2001 , 94 (2) :189-193. [3] Obrovac M N, Krause L J, Reversible Cycling of Crystalline Silicon Powder[J],Journal of the Electrochemical Society, 2007,154 (2) :A103-A108.

[4] Graetz J, Ahn C C, Yazami R, and Fultz B, Nanocrystalline and Thin Film Germanium Electrodes with High Lithium Capacity and High Rate Capabilities[J], 2004, 151 (5): A698-A702.

[5] Wolfenstine J,Foster D,Read J, Behl W K, and Luecke W, Experimental confirmation of the model for microcracking during lithium charging in single-phase alloys[J], Journal of Power Sources, 2000 , 87 (1-2) :1-3.

[6] Liu Y, Hudak N S, Huber D L, Limmer S J, Sullivan J P, and Huang J Y, In situ

transmission electron microscopy observation of pulverization of aluminum nanowires

and evolution of the thin surface Al2O3 layers during lithiation-delithiation cycles[J], Nano Letters, 2011, 11 (10) :4188.

[7] Wang Y, Qin Q Z, A Nanocrystalline NiO Thin-Film Electrode Prepared by Pulsed Laser Ablation for Li-Ion Batteries[J], Journal of the Electrochemical Society, 2002,149 (7):A873-A878.

[8] Fu Z W, Wang Y, Zhang Y, and Qin Q Z, Electrochemical reaction of nanocrystalline Co3O4, thin film with lithium[J], Solid State Ionics, 2004,170:105-109.

[9] Huggins R A, Advanced batteries: Materials science aspects[M], SpringerBerlin, 2009. [10] Lee S J, Lee J K, Chung S H, Lee H Y, Lee S M, and Baik H K, Stress effect on cycle properties of the silicon thin-film anode[J],Journal of Power Sources, 2001, 97: 191-193. [11] Winter M, Besenhard J O, ChemInform Abstract: Electrochemical Lithiation of Tin and Tin‐Based Intermetallics and Composites[J], Electrochimica Acta, 1999. 45:31-50. [12] Yoshio M, Tsumura T, Dimov N, Electrochemical behaviors of silicon based anode material[J], Journal of Power Sources, 2006 , 153 (2) :375-379.

[13] Wang D Y, Wu X D, Wang Z X, and Chen L Q, Cracking causing cyclic instability of lifepo 4, cathode material[J], Journal of Power Sources, 2005, 140 (1) :125-128. [14] Datta M K, Maranchi J, Chung S J, Epur R, Kadakia K, and Jampani P, Amorphous silicon-carbon based nano-scale thin film anode materials for lithium ion batteries[J], Electrochimica Acta, 2011, 56 :4717-4723.

[15] Zhou Y N, Li W J, Chen H J, Liu C, Zhang L, and Fu Z, Nanostructured nisi thin films as a new anode material for lithium ion batteries[J], Electrochemistry Communications, 2011,13 (6) :546-549.

[16] Imai Y, Watanabe A, Energetics of compounds related to Mg 2 Si as an anode material for lithium-ion batteries using first principle calculations[J], Journal of Alloys & Compounds, 2011, 509 (30) :7877-7880.

[17] Prussin S, Generation and Distribution of Dislocations by Solute Diffusion[J], Journal of Applied Physics, 1961 ,32(10):1876-1881.

[18] Bourderau S, Brousse T, Schleich D M, Amorphous silicon as a possible anode material for Li-ion batteries[J], Journal of Power Sources, 1999 ,81-82 (9) :233-236. [19] Chen L B, Xie J Y, Yu H C, and Wang T H, An amorphous si thin film anode with high capacity and long cycling life for lithium ion batteries[J], Journal of Applied Electrochemistry, 2009 , 39 (8) :1157-1162.

[20] He Y, Yu X, Li G, Wang R, Li H, and Wang Y, Shape evolution of patterned amorphous and polycrystalline silicon microarray thin film electrodes caused by lithium insertion and extraction[J], Journal of Power Sources, 2012 , 216 (11) :131-138.

[21] Shenoy V B, Johari P, Qi Y, Elastic softening of amorphous and crystalline Li-Si Phases with increasing Li concentration: A first-principles study[J], Journal of Power Sources, 2010,

195 (19) :6825-6830.

[22] Zeng Z, Liu N, Zeng Q, Ding Y, Qu S, and Cui Y, Elastic moduli of polycrystalline li15si4, produced in lithium ion batteries[J], Journal of Power Sources, (2013), 253 :431-431. [23] Hopcroft M A, Nix W D, Kenny T W, What is the Young's Modulus of Silicon[J], Journal of Microelectromechanical Systems, 2010 ,19(2):229-238.

因篇幅问题不能全部显示,请点此查看更多更全内容