# Modeling the Lithium Ion/Electrode Battery Interface Using Fick’s Second Law of Diffusion, the Laplace Transform, Charge Transfer Functions, and a [4, 4] Padé Approximant

Authors: John H. Summerfield; Charles N. Curtis
Journal: International Journal of Electrochemistry
 (2015)

---

## Abstract

This work investigates a one-dimensional model for the solid-state diffusion in a LiC6/LiMnO2 rechargeable cell. This cell is used in hybrid electric vehicles. In this environment the cell experiences low frequency electrical pulses that degrade the electrodes. The model’s starting point is Fick’s second law of diffusion. The Laplace transform is used to move from time as the independent variable to frequency as the independent variable. To better understand the effect of frequency changes on the cell, a transfer function is constructed. The transfer function is a transcendental function so a Padé approximant is found to better describe the model at the origin. Consider ∂c(r,t)/∂t=D∂2c(r)/∂2r+(2/r)(∂c(r)/∂r).


---

## Body

## 1. Introduction

The Li-ion batteries in hybrid electric vehicles experience electric current pulses during their use. These transient electric pulses stress the crystal structure of the electrodes and cause lithium salt crystal growth, both of which inhibit easy ion diffusion. Also the lithium ion cells are attached in series of 72 cells. This added variable increases the chance of varying pulses. These stressors lead to shortened battery life. In order to build better batteries these pulses need to be better understood. The electric pulses have low frequency, between 0.01 Hz and two Hz [1].The purpose of this work is to present a relatively simple mathematical model of the Li+ ions in the electrolyte of a LiC6/LiMnO2 battery diffusing into the electrodes. Lithium ion battery electrolyte conductivity has been explored [2]. Software to help students better understand the chemistry of batteries is available as well [3]. The Li+ battery electrochemical reactions have also been addressed [4]. Most recently Li+ ion movement in the electrolyte has been investigated [5].The model presented here is appropriate for Analytical Chemistry [6], Instrumental Analysis [7], Physical Chemistry [8], or Electrochemistry [9] classes. In Analytical Chemistry it could be presented as an example of stationary electrode chemistry. In Instrumental Analysis it could be presented as part of the introduction to disk electrode chemistry. In Physical Chemistry it could be presented as a diffusion example. In an Electrochemistry class it could be explored as an example of electrode chemistry or as an aspect of battery systems management. In all courses the model is presented as an in-class lecture. Also it could be extended into an independent study project. This topic will be discussed at the end of the paper.Transport properties have been known for over a hundred years yet they continue to be relevant [10]. Specifically batteries rely on the movement of ions due to the difference in voltage from electrode to electrode. This motion is termed ion migration. The ions also move from spot to spot in the battery due to differences in their concentrations. This motion is characterized as solid-state diffusion when the ions move in and out of the electrodes.Mathematical modeling of migration and diffusion within a lithium ion battery is quite involved. The voltages and the solution concentrations must be simulated for the surface of the electrodes and for the electrolyte between them. These equations are made more complex because the kinetics and transport parameters that describe the surfaces and the highly concentrated electrolyte are nonlinear [11].Presently four differential equations are used to simulate the ion/electrode diffusion of a typical lithium ion battery that consists of a LiMnO2 cathode and a LiC6 anode, separated by the electrolyte, LiPF6 [12]. There is great interest in moving away from a graphite anode to a better charge carrying semiconductor and replacing the manganese in the cathode with less toxic silicon or sulfur. Both changes are made ultimately to increase electric current [13–16].This paper relies on one differential equation to model the Li-ion diffusion into the solid-state electrode. The equation is a form of Fick’s second law of diffusion [17]. Diffusion is stressed because during the cell’s operation a mean electric field is created at the electrode’s surface so diffusion plays an important role. Also, simplification of solid-state diffusion as a whole is a current topic [18].


## 2. Li-Ion/Electrode Model

There are three simple atomic models. The atom can be viewed as a plane, a cylinder, or a sphere. The lithium ion is best described as a sphere. The lithium ion’s concentration as a function of time is best described by the Laplacian in spherical coordinates:(1)∂cr,θ,ϕ,t∂t=1r2∂∂rr2∂cr,θ,ϕ,t∂r+1r2sin⁡θ∂∂θsin⁡θ∂cr,θ,ϕ,t∂θ+1r2sin2⁡θ∂2cr,θ,ϕ,t∂ϕ2,where r is the distance from the origin—the center of the lithium ion—θ is the polar angle measured down from the North Pole, and ϕ is the azimuthal angle.For the case of one-dimensional solid-state diffusion (1) becomes(2)∂cr,t∂t=1r2∂∂rr2∂cr,t∂r,where D is the solid-state diffusion coefficient. The value of D is 2.0 × 10−12 cm2 s−1 for the LiC6 electrode and is 3.7 × 10−12 cm2 s−1 for the LiMnO2 electrode [12].To simplify, the derivatives on the right side of (2) are solved. This yields(3)∂cr,t∂t=D∂2cr∂r2+2r∂cr∂r.This partial differential equation has two boundary values. At r=0, the center of the lithium ion, ∂cr,t/∂r=0. At r=R, the radius of the lithium ion, 1 × 10−4cm, ∂cr,t/∂r=-jt/DaF, where for the negative LiC6 electrode(4)-jt=itAL.In (4), A is the total electrode plate area, L is the electrode length, and i(t) is the current flowing through the electrode [19]. Returning to the original boundary condition, D is the solid-state diffusion coefficient. The parameter a is the electrode volume fraction, 0.580 for the LiC6 electrode and 0.500 for the LiMnO2 electrode. The Faraday constant is F, 96,487 C/mol [12].


## 3. Laplace Transform

The Laplace transform converts the independent variable of an equation from time to frequency. All other independent variables are unaffected. The Laplace transform is a set of rules that are applied depending on whether the function is a polynomial, exponential derivative, and so forth [20]. For example, consider the boundary condition at r=R, ∂cr,t/∂r=-jt/DaF. The Laplace transform is (5)∂Cr,s∂r=-JsDaF,where the capital letter indicates that the function has undergone the transform. The independent variable is now s, frequency, rather than time.Turning to (3), a time derivative is transformed by multiplication by frequency plus a term that includes initial conditions. This change converts ∂cr,t/∂t to sCr,s-cr,0. The second term is zero since if the frequency is 0 Hz then there is no current pulse and no change in concentration. The transformation of the entire equation can be written as(6)∂2cr∂r2+2r∂cr∂r-sDCr,s=0.The solution to this differential equation is provided in the appendix. The final result atr=R is (7)CR,s=-RJssinh⁡Rs/DDaFRs/Dcosh⁡Rs/D-sinh⁡Rs/D.


## 4. Transfer Function

The interest is in the surface concentration of the ions on the electrode so a transfer function fromJ(s) to C(R,s) is sought. A transfer function is the ratio of output to input for a system. The roots of the numerator of a transfer function are the zeros of the function. Zeros will block the frequency transmission. The roots of the numerator are the singularities. Frequency goes to infinity at these points.The transfer function shows that the Laplace transform of the output is the product of the transfer function of the system and the transform of the input. In the transform domain the action of a linear system on the input is simply a multiplication with the transfer function. The transfer function is a natural generalization of the concept of gain—the increase in frequencies in the operating frequency range—of a system [21].The model yields the transfer function(8)CR,sJs=-Rsinh⁡Rs/DDaFRs/Dcosh⁡Rs/D-sinh⁡Rs/D.This transfer function is shown in Figure 1 for the anode and cathode.Figure 1
Transfer functions for a one-dimensional model of the solid-state diffusion in a typical electric car battery cell. Both transfer functions are moving away from the singularity ats=0. On discharge, as the electron output of the anode moves toward zero, the input into the cathode increases, making its transfer function decrease.


## 5. Padé Approximation

Polynomials—from a Maclaurin series—are often used in science, rather than a transcendental function, to speed computing time. Nonetheless the polynomials may oscillate which may increase error depending on step size. To help remedy this problem mathematicians have turned to a richer class of functions—the ratio of two polynomials.This ratio of polynomials is computed as a Padé approximant. The Padé approximation to the model’s transfer function is(9)Ps=∑m=0Nbmsm1+∑n=1Nansn=numsdens.The denominator begins at one so the solutions are normalized. The coefficients are found from the Maclaurin series. More specifically,(10)dens=∑k=02N+1cksk-nums=0,where ck’s are from the Maclaurin series [22].The transfer function has a singularity ats=0 so q2CR,s/Js is expanded, where q=Rs/D, to cancel the singularity [23]. The result for the anode is(11)q2CR,sJs=-RDaF-3-q25+q4175-2q67875+37q83031875.For N=2, with the variable change, (10) becomes (12)-3-b0+-3a1-15-b1q+1175-b2-a15-3a2q2+a1175-27875-a25q3+373031875-2a17875+a2175q4+37a13031875-2a27875q5+37a23031875q6=0.The q3 and q4 terms provide two equations for the two unknowns, a1 and a2. The value of a1 is 3/55. The value for a2 is 1/3465. With these two values, b0, b1, and b2 are found using the first three terms—the constant term, the q coefficient, and the q2 coefficient. The values are −3, −4/11, and −1/165, respectively.Substituting these values into (9) yields a Padé approximant,(13)q2CR,sJs=-RDaF-3-4q2/11-q4/1651+3q2/55+q4/3465.As a check, this is the expression returned byMathematica. This function is graphed for each electrode in Figure 2.Figure 2
The transfer function has been linearized. The singularity at the origin has been removed. Physically, at higher frequency more lithium ions leave the electrodes.This Padé approximant was used because it is the lowest order ratio that provides linearization of the modified charge transfer function. The even powers relate to the powers series forsinh⁡(x) and cosh⁡(x).


## 6. Conclusion

Degradation of Li+ ion battery electrodes is a current topic. The lithium ions move in and out of the electrodes distorting the electrodes’ crystal structure. This leads to a decrease in capacitance and thus cell life. There are also side reactions that form lithium salts on the electrodes and other reactions that dissolve the electrodes. Low frequency electrical pulses also disrupt the electrodes [24].This work describes a one-dimensional model of the lithium ions’ solid-state diffusion. It incorporates the low frequency electrical pulses that occur in the lithium ion cells of hybrid electric vehicles. Simple expressions for solid-state diffusion at low frequencies may provide guidance for models that incorporate the higher frequencies, 500 kHz to 500 MHz, used to charge the cells [1].In this work a transfer function was determined analytically in terms of hyperbolic functions. These nonpolynomial functions make it difficult to produce a standard transfer function in a simple variable. Instead a Padé approximant was created to produce a linear transfer function.The same type of analysis can be done for a nickel-metal hydride rechargeable cell. In this cell the ions are more cylindrical than spherical so cylindrical coordinates are relied on. For completeness, a planar molecule would be modeled best using Cartesian coordinates. Both differential equations are in the literature [25].



---

*Source: 496905-2015-06-17.xml*