OK, underwhelmed by responses, I figured out a crude means of getting the answer in React and treating the mineral as a kinetic reactant with the standard (1-Q/K) approach, using H(+) concentration instead of pH and then setting all relevant initial concentrations (H, K, Al, Si) to be on the order of 1e-9 mmolal in 1 kg of water. Is there a more elegant method?