Understanding the separator pore size inhibition effect on lithium dendrite via phase-field simulations

Yjie Li,Geng Zhng,Bin Chen,Wei Zho,Liting Sh,D Wng,Ji Yu,Siqi Shi,,d,∗∗

a School of Materials Science and Engineering,Shanghai University,Shanghai 200444,China

b Physical Science and Engineering Division,King Abdullah University of Science and Technology,Thuwal 23955-6900,Saudi Arabia

c Materials Genome Institute,Shanghai University,Shanghai 200444,China

d Zhejiang Laboratory,Hangzhou 311100,China

Keywords:Lithium-ion batteries Pore size Ionic distribution Phase-field simulations Dendrite growth

ABSTRACT Dendrite growth in lithium-ion batteries may bring thermal run-away especially at high current densities,which remains the major bottleneck to implement safe and fast charging for portable electronic devices or electronical vehicles.Designing dendrite inhibition separators with proper pore size is considered to be one of the most promising strategies to guarantee the battery safety.However,due to the impossible observation of lithium-ion distribution under separator by experiments,the underlying dendrite inhibition mechanism is still not fully understood.Here,we apply the phase-field model,which takes the separator phase into account to construct the electrochemical system total free energy,to study the ion re-distribution behavior of porous separator and understand the pore size inhibition effect on lithium dendrite.The numerical results indicate that separator with smaller pore size is beneficial to smoother electrodeposition,since the lithium-ion concentration on the electrode surface is more uniform under denser separator pores,when their sizes is larger than the critical nucleus.The proposed model could capture the physicochemical process of electrodeposition under multiphase structures,so it could also be used to explore dendrite growth under composite electrodes and composite solid electrolytes.

The rapid development of portable electronic devices,electronical vehicles and large-scale energy storage devices calls for advanced lithium-ion batteries (LIBs) with high energy density,reliable safety,and long cycle life [1–3].Nevertheless,lithium dendrite growth may bring poor cycle stability and thermal run-away especially at high current densities,which remains the major bottleneck to implement safe and fast charging for LIBs [4–6].To meet the rising demand for battery safety,separators,acting as physical barriers between electrodes,have been considered to be important components in suppressing dendrite [7,8].Many experimental researches have revealed the relation between separator structure and dendrite growth,including pore size [9–11],pore uniformity [12–16],surface coating [17–23],and electrolyte wettability[9,11,24–27].For example,Gaoet al.[11]fabricated polydopaminecoated polyacrylonitrile separator through spin coating and electrospinning.The fiber diameter and surface roughness of polyacrylonitrile fibers are decreased by the polydopamine coatings;thus the separator pore size is decreased and the lithium-ion flux is homogenized.Compared with polyacrylonitrile separator,batteries with polydopamine-coated polyacrylonitrile separator have more stable cycle life and better rate capability.Qinet al.[14]designed a porous water-based separator with uniform pore distribution,and demonstrated that it has better dendrite inhibition ability than the separator with uneven porous structure.Surface coating of separator can not only improve battery thermal stability,but also suppress dendrite growth.Yanget al.[10]reported a mesoporous silica-coated aluminum oxide separator.The perpendicular nanochannels of mesoporous silica films can effectively regulate the ionic flux and lead to uniform lithium deposition.The Li-Li battery with silica-coated aluminum oxide separator exhibits excellent cycling performance at ultrahigh current densities(10 mA/cm2) for more than 1600 h of operation.Similarly,Shenet al.[19]applied electrospraying and polymer photopolymerization to synthesize a chitosan,polyethylene oxide,and poly(triethylene glycol dimethacrylate) coated CelgardR○separator.The surface coatings provide a lithiophilic environment,which enables uniform lithium deposition,facilitates homogeneous solid electrolyte interphase layer,and thus prevents dendrite growth.

Among the above solutions,regulating separator pore size is the easiest way with lowest cost,and almost without additional weight.Eexperimental studies have shown that separators with smaller pore size can suppress dendrite growth;this is presumably attributed to the re-distribution of lithium ions by the dense pores.However,the dynamic lithium-ion distributions inside batteries cannot be directly observed by experiments,and thus the mechanisms behind the lithium inhibition behavior are still not fully understood.Hence it is necessary to apply computational method,which is neither as difficult as theoretical study nor as time-consuming and laborious as experimental study,to investigate the evolution during electrodeposition [28].The phase-field model is a computational model which describes microstructure evolution of material systems by using functions of space and time[29,30].During 2003–2012,the phase-field model has been developed and applied in electrodeposition by Monroe and Newman[31],Guyeret al.[32,33],Okajimaet al.[34,35].While these models assume linear kinetics,which are not applicable for systems far from equilibrium.At 2012,Lianget al.[36]firstly reported a nonlinear phase-field model,which takes the Butler-Volmer reaction kinetics into account,to simulate the electrode-electrolyte interface evolution during highly nonequilibrium processes.Later,they applied the nonlinear phase-field model to describe the electrodeposition process in electrochemical systems [37].At 2015,Chenet al.[38]further developed the above model to investigate the dendritic patterns during the lithium electrodeposition process,and the model is first verified by comparison with the Nernst equation in one-dimensional equilibrium system.For the application of phase-field model in dendrite growth under separators,Janaet al.[39]applied phase-field model to explore the effect of pore size on dendrite growth,and identified four regimes of dendrite growth: suppression regime,permeable regime,penetration regime,and short circuit regime.However,this study ignored the lithium-ion concentration during electrodeposition,which is important to understand dendrite growth mechanism.

In this paper,the influence of separator pore size on dendrite growth is studied by the phase-field model using self-written MATLAB code.This extended model takes the separator phase into account to construct the total free energy of the electrochemical system.The simulation results are firstly compared with previous studies conducted by integrated software to check validity of the model and accuracy of the code.Then we study the ion redistribution behavior of porous separator and the pore size inhibition effect on lithium dendrite.

Generally,the simulations for electrochemical dendrite growth include the following steps: (1) constructing the total free energy of the electrochemical system,(2) numerically solving the evolution equations of all fields with MATLAB,and (3) visualizing the simulation results.The figures of the phase-fields,ion concentrations,and electric potential at different times can be obtained from visualization of the numerical solution.The flow chart of the code is shown in Fig.1.

We introduce two phase-field variablesξandφ(ranging from 0 to 1) and a concentration setci(i=Li,Li+and anionAm−) to describe the system.(ξ=1,φ=0)denotes the electrode,(ξ=0,φ=1)denotes the separator,and(ξ=0,φ=0)denotes the electrolyte.The total free energy of the system can be expressed as

whereΦis the electric potential,fch,fgradandfelecare the Helmholtz free energy density,gradient energy density,and electrostatic energy density,respectively.

Fig.1.Flow chart of phase-field simulation code.

The chemical free energy density is expressed as follows:

whereμ◦iis the reference chemical potential of componenti,Ris the molar gas constant,andTis the temperature.The last three terms of the right side represents interphase barrier potentials,withW01,W02andW12being constants.The gradient energy density is expressed as:

whereκ1is the gradient energy coefficients,δis the anisotropic strength,ωis the anisotropic mode,θis the angle between the interface normal vector and the reference axis.The electrostatic energy density is expressed as

whereFis Faraday’s constant andziis the valence.The electrochemical reaction rate can be expressed from the Butler-Volmer equation.The evolution equations forξandφcan be described as:

whereLis the interfacial mobility between the electrode and separator,Rηis the reaction constant,αand 1-αare the chargetransfer coefficients,c0is the initial concentration of the electrolyte,h(ξ)=ξ3(6ξ2-15ξ+10) is an interpolating function.The overpotential is defined asη=Φ–Φ0,Φis the electric potential of Li metal anode,Φ0is the equilibrium potential without electric current.The evolution equation of Li+concentration is

whereDeff=(h(ξ)+h(φ))DLi+[1 −h(ξ)−h(φ)]Deis the effective diffusion coefficient,DLi+is the diffusion coefficient of Li+in the electrolyte,Deis the diffusion coefficient of Li+in the electrode,andKis the accumulation constant.The evolution equation of electrostatic potential is described by the law of charge conservation:

whereσeff=h(ξ)σLi+[1-h(ξ)]σeis the effective electric conductivity,σLiandσeare the electric conductivity of electrode and electrolyte,respectively.

The numerical solutions of Eqs.5–8 are obtained by using the finite difference method on MATLAB.The simulation is conducted in a two-dimensional region with a size of 16×16 μm2,and the grid spacing is 0.08 μm.The parameters used are listed in Table 1.The adiabatic boundary condition is used for the four boundaries of the phase-field variable,and left and right boundaries of the concentration and electric potential.The top and bottom boundaries of concentration are set as 1 and 0,respectively.The top and bottom boundaries of electric potential are set as 0 V and −0.7 V,respectively.

Table 1 Parameters in the phase-field model.

In order to verify the feasibility of the electrochemical phasefield model,we conducted a simulation without separator and compared the results with the previous work.As shown in Fig.2a,the blue area represents the electrolyte.The upper and lower boundary conditions are set as the Dirichlet boundary conditions,and three initial semicircular nuclei are set on the electrode.The evolutions of phase-field variable (ξ),Li+concentration(cLi+) and electrical potential (Φ) from 0 to 2 s are shown in Fig.2.Due to the tetragonal crystal structure of Li (BCC),four equal branches will appear.Since the bottom margin of the simulation region represents the solid electrode phase,there are only three branches growing towards the electrolyte phase in electrodeposition (Fig.2b) [40].The concentration/electric potential gradient at the interface is the electrodeposition driving force.

From Figs.2b–d,it can be observed that the growth rate at the dendrite tip is much faster than that at the bottom,since the diffusion path of lithium ions is shorter for the tips [41].When the secondary dendrite get in contact with each other,the length of the secondary dendrite between two nucleation sites is much less than that of the outer secondary dendrite due to the competition between adjacent dendrite.This is consistent with the results in previous work through COMSOL or MOOSE [38].The ionic concentration and potential in electrolyte is also in good agreement with the work of Chenet al.[38,41,42].These simulation results show validity of the proposed electrochemical phase-field model.

Fig.2.Evolution of (a–d) phase-field variable,(e–h) Li+ concentration,and (i–l)electric potential without separator.Scale bar: 2 μm.

Fig.3.Phase-field variable at 5 s under different pore sizes: (a) 0.21 μm,(b)0.42 μm,(c) 0.84 μm,(d) 1.05 μm,(e) 2.1 μm,(f) 6 μm,and (g) maximum dendrite length versus pore sizes.Scale bar: 2 μm.

Next,the effect of pore size on Li dendrite growth is explored by six separator pore sizes (0.21,0.42,0.84,1.05,2.1 and 6 μm).At 5 s,the Li dendrite morphology under different pore sizes is shown in Figs.3a–f.The line chart of maximum dendrite length with respect to pore size is plotted in Fig.3g.It can be observed that with the increase of pore size,the maximum height of dendrite increases gradually.From Figs.3e and f,when the pore size is too large,the dendrite can easily puncture the separator;thus inducing a short circuit of the battery.However,the pore size is not better if too small.By comparing Figs.3a and b,it is seen that the dendrite length under pore size of 0.21 μm is higher than that of 0.42 μm,because dendrite can go through the separator if the pore size is close to the critical nuclei [39].

To further analyze the mechanism behind dendrite growth,we plotted the lithium-ion concentration under different pore sizes from 0 to 5 s,as shown in Fig.4.The initial lithium-ion concentrations under different pore sizes are all the same.At 1 s,lithium ions are consumed due to electrochemical reaction and deposited on the electrode,resulting in the low lithium-ion concentration near the electrode.Concentration difference drives lithium-ion migration from the electrolyte to the anode.Lithium ions can pass through the separator pores but cannot pass through the polymer matrix of the separators.For maintaining the same porosity,when the pore size increases,the corresponding polymer matrix area(green strip area) also increases.Therefore,if the matrix area is too large,the lithium-ion concentration under it will be reduced obviously,resulting in uneven distribution of lithium-ions,as shown in Figs.4n and q.The uneven distribution of lithium ions aggravates the tip effect,which in turn increases ionic heterogeneity and promotes dendrite growth.

Fig.4.Lithium-ion concentration at 0,1 and 5 s under different pore sizes: (a–c)0.21 μm,(d–f) 0.42 μm,(g–i) 0.84 μm,(j–l) 1.05 μm,(m–o) 2.1 μm,and (p–r) 6 μm.Scale bar: 2 μm.

In this paper,an electrochemical phase-field model and its corresponding MATLAB code are developed to explore the electrodeposition behavior under separators with different pore size.The results indicate that separators with smaller pore size are beneficial to smoother electrodeposition.While when the pore size decreased close to the critical nuclei,the dendrite growth will be promoted instead.This simulation provides significant strategies to material scientists for developing separators with optimum dendrite inhibition effect.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Nos.52102280,U2030206,11874254,51622207),Shanghai Pujiang Program (No.2019PJD016),Foundation of China Academy of Engineering Physics-Key Laboratory of Neutron Physics(No.2019BB07),and Scientific Research Project of Zhijiang Laboratory (No.2021PE0AC02).It was also supported by funding from King Abdullah University of Science and Technology (KAUST).