Saturation index

[tianshiyu]

First, I want to thank dlparkhurst for helping me. But I couldn't reply to you in my last post, I had to created a new post.Now I have another problem need to consult you. I use EQUILIBRIUM_PHASES to make the SI=0(hydroxyapatite), but hydroxyapatite only settle at a High supersaturation and I didn't know the value. So in the simulation process, the hydroxyapatite will settled at a low pH(pH6-7) than reality(pH7-9). Is there any way to make the simulation more realistic?
thank you again for your reply.

[Pak]

Maybe changing the themodynamic equilibrium constant? not sure.

[tianshiyu]

But the themodynamic equilibrium constant is the database default. Is there any method to change it？

[dlparkhurst]

You can use a PHASES data block in your input file that redefines the phase with a different equilibrium constant. Just copy the definition from the database and put it in your input file.

The equilibrium constant for hydroxyapatite does vary among the databases, so varying within that range may be reasonable. If you need to change the log K more than that, you may be missing something.

[tianshiyu]

Thank you a lot. I think I might have a clue!