R_Li2S_s k_Li2S_s*max(eps_Li2S_s,eps)*(max(Li_1p,eps)^2*max(S_2m,eps)-Ksp_Li2S_s) "Reaction rate, Li2S(s) precipitation" R_S8_s k_S8_s*max(eps_S8_s,eps)*(S8-Ksp_S8_s) "Reaction rate, S8(s) precipitation"