Thank you Jia, this is incredibly helpful. I have been working through several issues using this simpler approach. I switched over to React some of the issues I am having. I reacted my initial solution with reactants from the source and I am having an issue I do am trying to diagnose. I am having charge balancing issues so I am using chloride for the make-up ion. When I run the reaction using 5 percent Ferrihydrite (swapped in for Fe+++ in the basis pane), the reaction seems to run just fine (see attached output file). When I run the reaction with greater than 5 percent ferrihydrite (see .rea file), I am getting a convergence error, with chloride being the residual. Not sure why adding more Ferrihydrite give me a convergence error, but I am thinking this is a real problem for me in the X1t model I am running. The aquifer has 10-15 percent ferrihydrite so I am trying to simulate this and cannot seem to. Do you have any suggestions? Overall, charge balance and residuals have been a real issue for me and so I am hoping you can help me troubleshoot?
React_output_sx w initial_0.05ferri.txt
10 percent ferrihydrite failure.rea