Spinodal decomposition and domain coarsening in a multilayer Cahn-Hilliard model for lithium intercalation in graphite

During the intercalation of lithium in layered host materials such as graphite, lithium atoms can move within the plane between two neighboring graphene sheets, but cannot cross the sheets. Repulsive interactions between atoms in different layers lead to the existence of ordered phases called "stages,"with stage n consisting of one filled layer out of n, the others being empty. Such systems can be conveniently described by a multilayer Cahn-Hilliard model, which can be seen as a mean-field approximation of a lattice-gas model with intra- and interlayer interactions between lithium atoms. In this paper, the dynamics of stage formation after a rapid quench to lower temperature is analyzed, both by a linear stability analysis and by numerical simulation of the full equations. In particular, the competition between stages 2 and 3 is studied in detail. The linear stability analysis predicts that stage 2 always grows the fastest, even in the composition range where stage 3 is the stable equilibrium state. This is borne out by the numerical simulations, which show that stage 3 emerges only during the nonlinear coarsening of stage 2. Some consequences of this finding for the charge-discharge dynamics of electrodes in batteries are briefly discussed.