The choice of a thermodynamic formulation dramatically affects modelled chemical zoning in minerals

Quantifying natural processes that shape our planet is a key to understanding the geological observations. Many phenomena in the Earth are not in thermodynamic equilibrium. Cooling of the Earth, mantle convection, mountain building are examples of dynamic processes that evolve in time and space and are driven by gradients. During those irreversible processes, entropy is produced. In petrology, several thermodynamic approaches have been suggested to quantify systems under chemical and mechanical gradients. Yet, their thermodynamic admissibility has not been investigated in detail. Here, we focus on a fundamental, though not yet unequivocally answered, question: which thermodynamic formulation for petrological systems under gradients is appropriate – mass or molar?  We provide a comparison of both thermodynamic formulations for chemical diffusion flux, applying the positive entropy production principle as a necessary admissibility condition. Furthermore, we show that the inappropriate solution has dramatic consequences for understanding the key processes in petrology, such as chemical diffusion in the presence of stress gradients.