Phase Field Crystal Modeling of Coherent Interfaces in Lithium Batteries
Authors: Renuka Balakrishna Ananya , Carter W. Craig 
Published: February 2020

The aim is to study the role of material crystallography in the formation of coherent interfaces in lithium ferrous phosphate (LFP) electrodes. In a typical LFP electrode, lithiated crystals with different orientations nucleate heterogeneously and little is known about how they interact to form coherent interfaces. The current work provides insights on how lithiated crystals interact, and the effects of coherency stresses at interface and grain boundaries.
A key limitation of the battery performance is the structural instability of the electrodes with continuous usage, resulting in its failure [1]. In a battery cycle the electrode microstructures typically operate in a two-phase window that comprises of a moving coherent interface. Theoretical studies [2-4] indicate that a moving coherent interface causes anisotropic expansion of electrode [2], and affects Li-ion kinetics [4]. These factors affect battery performance, and it is therefore important to explore how interfaces form.
In recent years, phase field models have been developed to describe microstructures in electrodes, and use Li-ion concentration as the order parameter [3-6]. While these models reveal insights on nucleation barrier [4] and stress generation at the coherent interface [6], the anisotropic properties arising from material crystallography is not captured. This limits understanding on the role of coherent interfaces on in-situ battery operation.
In the present study, we apply a phase field crystal (PFC) model to investigate the formation of coherent interfaces in LFP electrodes [7,8]. PFC models offer the advantage of describing crystal growth on atomic scales, while solving diffusive time equations [7]. This modeling method has been applied to describe grain boundary melting [9], graphene defect structures [8], and has been well established to solve for a range of crystal symmetries. Here, we define a two-dimensional spatial density field, ψ, of a group of atoms to represent the orthorhombic (rectangular in 2D) symmetry of LFP crystals. A challenge here is to model ψ (order parameter) as a function of Li-ion concentration to differentiate the lithiated and unlithiated phases of the electrode. The free energy of the system is then described as a function of ψ, and is equal to the sum of bulk and excess energy terms [8]. The excess energy terms incorporate the two-point and the three-point interactions, which set the lattice spacings (a = 6.008Å, c = 10.334Å), and the lattice symmetry [10]. Finally, we control the evolution of the density field by a generalized equation of motion for conserved fields [8].
The coherency stresses and diffusion barriers at interfaces are captured by the PFC simulation. The results provide insights and strategies for doping electrode materials so as to reduce coherency stresses in electrodes.
References
[1] X. H. Liu, Y. Liu, A. Kushima, S. Zhang, T. Zhu, J. Li, and J. Y. Huang. In Situ TEM Experiments of Electrochemical Lithiation and Delithiation of Individual Nanostructures. Adv. Energy Mater., 2, 722-741, 2012.
[2] A. Van der Ven, K. Garikipati, S. Kim, and M. Wagemaker. The Role of Coherency Strains on Phase Stability in LixFePO4: Needle Crystallites Minimize Coherency Strain and Overpotential. J. Electrochem. Soc., 156, A949, 2009.
[3] M. Tang, W. C. Carter, J. F. Belak, Y-M Chiang. Modeling the competing phase transition pathways in nanoscale olivine electrodes. Electrochimica Acta, 56(2) 969-976, 2010
[4] D. A. Cogswell and M. Z. Bazant. Theory of Nucleation in Phase-separating Nanoparticles. Nano Lett.13, 3036-3041, 2013.
[5] E. Bohn, T. Eckl, M. Kamlah, and R. McMeeking. A Model for Lithium Diffusion and Stress Generation in an Intercalation Storage Particle with Phase Change. J. Electrochem. Soc., 160, A1638-A1652, 2013.
[6] D.T O’Connor, M.J. Welland, W.K. Liu, and P.W. Voorhees. Phase transformation and fracture in single LixFePO4 cathode particles: a phase-field approach to Li-ion intercalation and fracture. Model. and Simul. in Mat. Sci. and Engg, 24(3), 035020, 2016.
[7] K. R. Elder, and M. Grant. Modeling elastic and plastic deformations in nonequilibrium processing using phase field crystals. Phys. Rev. E, 70(5), 051605, 2004.
[8] M. Seymour, and N. Provatas. Structural phase field crystal approach for modeling graphene and other two-dimensional structures. Phys. Rev. B, 93(3), 035447, 2016.
[9] J. Berry, K. R. Elder, and M. Grant. Melting at dislocations and grain boundaries: A phase field crystal study. Phys. Rev. B, 77, 224114, 2008.
[10] A. K. Padhi, K. S. Nanjundaswamy, and J. B. D. Goodenough. Phospho-olivines as positive-electrode materials for rechargeable lithium batteries. J. of the electrochemical soc., 144(4), 1188-1194, 1997