Create a SURFACE for sorption in phreeqpython

[ELarrea587]

Hello,
I am modeling the sorption of dissolved metals to ferrihydrite surfaces. Is there a way to use this code to access the SURFACE functionality of PHREEQC? Is there a method to create a surface or a work-around function? Thank you for time!

[martinwk]

You could always use pp.ip.run_str to define your surface in normal PHREEQC code. Only drawback is that you cannot access the properties of the surface directly. You can only assert the concentration of surface species in the solution if you access the species attribute of a Solution instance. But not the charge on the surface and all those other things phreeqc calculates about the surface for you.

[behzaadh]

You could always use pp.ip.run_str to define your surface in normal PHREEQC code. Only drawback is that you cannot access the properties of the surface directly. You can only assert the concentration of surface species in the solution if you access the species attribute of a Solution instance. But not the charge on the surface and all those other things phreeqc calculates about the surface for you.

Could you please share your code if you have an example?

[UKOCM001]

Hello!
Has anyone managed to do this, in terms of modelling sorption of metals to various surfaces with phreeqpython? Or has anyone got an alternative solution they are willing to share? Thank you!