Dissolution kinetics as a function of the Gibbs free energy of reaction: An experimental study based on albite feldspar

Here we report on an experimental investigation of the relation between the dissolution rate of albite feldspar and the Gibbs free energy of reaction, Delta G(r). The experiments were carried out in a continuously stirred flow-through reactor at 150 degrees C and pH((150) (degrees C)) 9.2. The dissolution rates R are based on steady-state Si and Al concentrations and sample mass loss. The overall relation between Delta G(r) and R was determined over a free energy range of -150 < Delta G(r) < -15.6 kJ mol(-1). The data define a continuous and highly non-linear, sigmoidal relation between R and Delta G(r) that is characterized by three distinct free energy regions. The region furthest from equilibrium, delimited by -150 < Delta G(r) < -70 kJ mol(-1), represents an extensive dissolution rate plateau with an average rate (R) over bar = 1.0 x 10(-8) mol m(-2) s(-1). In this free energy range the rates of dissolution are constant and independent of Delta G(r), as well as [Si] and [Al]. The free energy range delimited by -70 <= Delta G(r) <= -25 kJ mol(-1), referred to as the 'transition equilibrium' region, is characterized by a sharp decrease in dissolution rates with increasing Delta G(r,) indicating a very strong inverse dependence of the rates on free energy. Dissolution nearest equilibrium, defined by Delta G(r) > -25 kJ mol(-1), represents the 'near equilibrium' region where the rates decrease as chemical equilibrium is approached, but with a much weaker dependence on Delta G(r). The lowest rate measured in this study, R = 6.2 x 10(-11) mol m(-2) s(-1) at Delta G(r) = -16.3 kJ mol(-1), is more than two orders of magnitude slower than the plateau rate. The data have been fitted to a rate equation (adapted from Burch et al. [Burch, T. E., Nagy, K. L., Lasaga, A. C., 1993. Free energy dependence of albite dissolution kinetics at 80 degrees C and pH 8.8. Chem. Geol. 105, 137-162)) that represents the sum of two parallel reactions R = k(1) [1 - exp(-ng(m1))] + k(2)[1 - exp(-g)](m2), where k(1) and k(2) are rate constants that have been determined by regression, with values 1.02 x 10(-8) and 1.80 x 10(-10) mol m(-2) s(-1), g vertical bar Delta G(r)vertical bar/RT is a dimensionless number, and n, m(1), and m(2) are adjustable fitted parameters (n = 7.98 x 10(-5), m(1) = 3.81 and m(2) = 1. 17). Based on measurements of the temporal evolution of R-si and R-Al for each experiment, steady-state dissolution rates appear to be congruent at all Delta G(r). In contrast, non-steady-state dissolution is incongruent, and is related to Delta G(r.) Scanning electron microscopy (SEM) images of post-reaction grain surfaces indicate that dissolution close to equilibrium (Delta G(r) > -25 kJ mol(-1)) resulted in the precipitation of a secondary crystalline phase, but there are no indications that this altered the measured R-Delta G(r) relation. (c) 2005 Elsevier Inc. All rights reserved.