where ψ chem and ψ mech describe the bulk chemical and mechanical energy density, respectively. The maximum molar density of Li lattice sites is C max = 49200 mols m -3 . c Li is the molar density of Li at a spatial position at a given time. The Li occupancy fraction θ can then be calculated as θ = c Li /C max . F is the total deformation gradient.