Hello Jia,
logK 9.3 is really a wrong value, which is the newest change we tried, and the right log K value is really the 16.57. In addition, we don't need any surface adsorption here, so I give a really small value to the porosity. I think I have control the condition well and they are same as I tried in Mingteq, but I get a wrong result. I can't get the solution from your suggestions. Did there any other problems I have made here? I'm looking forward to hearing from you.
Thank you again,
Cheng Zhang