Preprints
https://doi.org/10.5194/egusphere-2026-4987
https://doi.org/10.5194/egusphere-2026-4987
01 Sep 2026
 | 01 Sep 2026
Status: this preprint is open for discussion and under review for Solid Earth (SE).

A unified explicit formula for temperature-, pressure-, and composition-driven reaction-front propagation

Liudmila Khakimova, Stefan Schmalholz, and Yury Podladchikov

Abstract. Metamorphic reactions, phase transitions, and fluid–rock interactions play key roles in many geodynamic processes. Such reactions can be driven by changes in temperature, pressure, or chemical composition, and hence by different thermodynamic variables, yet commonly propagate through rocks as moving reaction fronts. However, no unified theory predicts their propagation. Here, we show that reaction front propagation can be described within a single conservative moving-boundary framework. The theory couples diffusive transport with local equilibrium thermodynamics through the Rankine–Hugoniot condition and yields explicit analytical expressions for reaction-front position and velocity. Unlike classical Stefan-like solutions, which generally involve error functions and numerical solutions of a transcendental equation, our formulation directly relates front propagation to transport properties and equilibrium thermodynamics. The theory predicts the characteristic square-root-of-time, √t, scaling observed for temperature-, pressure-, and composition-driven fronts and shows that the diffusivity governing front propagation generally differs from the physical transport diffusivity. Instead, an effective front diffusivity combines diffusive transport with the thermodynamic transformation across the front. This combination is quantified by a dimensionless L number linking equilibrium thermodynamics to front dynamics. Consequently, slow front propagation does not necessarily imply slow transport but may reflect large changes in thermal energy or mass associated with the reaction. Furthermore, sharp reaction fronts can propagate under diffusion-dominated conditions because finite thermodynamic jumps are explicitly represented at the moving interface, unlike conventional advection–diffusion–reaction models, in which sharp fronts typically require advection-dominated transport. The framework provides a unified basis for predicting and interpreting reaction-front propagation in experiments and geological systems.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Share
Liudmila Khakimova, Stefan Schmalholz, and Yury Podladchikov

Status: open (until 13 Oct 2026)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
Liudmila Khakimova, Stefan Schmalholz, and Yury Podladchikov
Liudmila Khakimova, Stefan Schmalholz, and Yury Podladchikov
Metrics will be available soon.
Latest update: 01 Sep 2026
Download
Short summary
In “A unified explicit formula for temperature-, pressure-, and composition-driven reaction-front propagation”, we ask whether different changes in rocks share a common rule. Using mathematics, simulations and published data, we show that boundaries driven by heat, pressure or chemistry move with the square root of time. Slow-moving boundaries do not necessarily mean slow transport and sharp boundaries can form without strong fluid flow. This changes how we infer rates and timescales from rocks.
Share