# React script, saved Tue Jul 09 2013 by bfarrell data = thermo_ladder.tdat time start = 0 day, end = 10 day temperature = 25 C decouple Fe+++ decouple Mn+++ decouple MnO2(aq) decouple Birn(aq) decouple MnO4- decouple MnO4-- H2O = 1 free kg pH = 8.7 swap O2(g) for O2(aq) O2(g) = .2 fugacity Na+ = .7 molal Cl- = .7 molal balance on Cl- Mn++ = 20 umolal swap Manganite for Mn+++ Manganite = 1e-9 free g swap Fe(OH)3(ppd) for Fe+++ Fe(OH)3(ppd) = 1 free mmol fix pH fix fugacity of O2(g) kinetic redox-1 rxn = "Mn++ + .25*O2(aq) + 1.5*H2O -> Manganite + 2*H+" rate_con = 2e4 surface = 2e6 catalyst = Fe(OH)3(ppd) mpower(Mn++) = 1 mpower(O2(aq)) = 1 apower(OH-) = 2 suppress Bixbyite Hausmannite Hematite Goethite suppress Fe(OH)3(soil) suffix _catalysis