Kinetic Monte Carlo simulation of surface heterogeneity in graphite anodes for lithium-ion batteries: Passive layer formation

The properties and chemical composition of the solid-electrolyte-interface (SEI) layer have been a subject of intense research due to their importance in the safety, capacity fade, and cycle life of Li-ion secondary batteries. Kinetic Monte Carlo (KMC) simulation is applied to explore the formation of the passive SEI layer in the tangential direction of the lithium- ion intercalation in a graphite anode. The simulations are found to consistent with observations in the literature that the active surface coverage decreases with time slowly in the initial stages of the battery operation, and then decreases rapidly. The effects of operating parameters such as the exchange current density and temperature on the formation of the passive SEI layer are investigated. The active surface coverage at the end of each charging cycle was initially lower at higher temperature, but remained constant for more cycles. The temperature that optimizes the active surface in a lithium-ion battery at Cycle 1 can result in less active surface area for most of the battery life.

[1]  M. Kalos,et al.  Time evolution of a quenched binary alloy. IV. Computer simulation of a three-dimensional model system , 1974 .

[2]  A. B. Bortz,et al.  A new algorithm for Monte Carlo simulation of Ising spin systems , 1975 .

[3]  W. H. Weinberg,et al.  Theoretical foundations of dynamical Monte Carlo simulations , 1991 .

[4]  H. Hazan,et al.  First Principles , 2011, For the University.

[5]  Ralph E. White,et al.  Corrosion initiation time of steel reinforcement in a chloride environment—A one dimensional solution , 1997 .

[6]  Anton Van der Ven,et al.  Lithium Diffusion in Layered Li x CoO2 , 1999 .

[7]  Herbert L Case,et al.  An accelerated calendar and cycle life study of Li-ion cells. , 2001 .

[8]  Ralph E. White,et al.  Mathematical modeling of the capacity fade of Li-ion cells , 2003 .

[9]  Ralph E. White,et al.  Solvent Diffusion Model for Aging of Lithium-Ion Battery Cells , 2004 .

[10]  R. Braatz,et al.  Coarse-Grained Kinetic Monte Carlo Simulation of Copper Electrodeposition with Additives , 2004 .

[11]  Ganesan Nagasubramanian,et al.  Modeling capacity fade in lithium-ion cells , 2005 .

[12]  Ralph E. White,et al.  Analysis of capacity fade in a lithium ion battery , 2005 .

[13]  Anton Van der Ven,et al.  Thermodynamics of spinel LixTiO2 from first principles , 2005 .

[14]  Ralph E. White,et al.  Review of Models for Predicting the Cycling Performance of Lithium Ion Batteries , 2006 .

[15]  Timothy O. Drews,et al.  Stochastic Simulation of the Early Stages of Kinetically Limited Electrodeposition , 2006 .

[16]  R. Braatz,et al.  Monte Carlo Simulation of Kinetically Limited Electrodeposition on a Surface with Metal Seed Clusters , 2007 .

[17]  Anton Van der Ven,et al.  Nondilute diffusion from first principles: Li diffusion in Li x TiS 2 , 2008 .

[18]  Richard D. Braatz,et al.  A hybrid multiscale kinetic Monte Carlo method for simulation of copper electrodeposition , 2008, J. Comput. Phys..

[19]  Ralph E. White,et al.  Capacity fade analysis of a lithium ion cell , 2008 .

[20]  V. Subramanian,et al.  Mathematical Model Reformulation for Lithium-Ion Battery Simulations: Galvanostatic Boundary Conditions , 2009 .

[21]  Anton Van der Ven,et al.  Phase stability and nondilute Li diffusion in spinel Li 1 + x Ti 2 O 4 , 2010 .

[22]  Richard D. Braatz,et al.  Parameter Estimation and Capacity Fade Analysis of Lithium-Ion Batteries Using First-Principles-Based Efficient Reformulated Models , 2009 .