Simulation of intercalation and phase transitions in nano-porous, polycrystalline agglomerates

Optimal microstructure design of battery materials is critical to enhance the performance of batteries for tailored applications such as high power cells. Accurate simulation of the thermodynamics, transport, and electrochemical reaction kinetics in commonly used polycrystalline battery materials remains a challenge. Here, we combine state-of-the-art multiphase field modelling with the smoothed boundary method to accurately simulate complex battery microstructures and multi-phase physics. The phase-field method is employed to parameterize complex open pore cathode microstructures and we present a formulation to impose galvanostatic charging conditions on the diffuse boundary representation. By extending the smoothed boundary method to the multiphase-field method, we build a simulation framework which is capable of simulating the coupled effects of intercalation, anisotropic diffusion, and phase transitions in arbitrary complex polycrystalline agglomerates. This method is directly compatible with voxel-based data, e.g. from X-ray tomog-raphy. The simulation framework is used to study the reversible phase transitions in LiX NiO2 in dense and nanoporous agglomerates. Based on the thermodynamic consistency of phase-field approaches with ab-initio simulations and the open circuit potential, we reconstruct the Gibbs free energies of four individual phases (H1, M, H2 and H3) from experimental cycling data. The results show remarkable agreement with previously published DFT results. From charge simulations, we discover a strong influence of particle morphology on the phase transition behaviour, in particular a shrinking core-like behaviour in dense polycrystalline structures and a particle-by-particle mosaic behavior in nanoporous samples. Overall, the proposed simulation framework enables the detailed study of phase transitions in intercalation materials to enhance microstructure design and fast charging protocols.