Ploting pe (or Eh) vs reaction temperature

[nesa]

Hi there. I am trying to model Pb, Zn, Cu, Fe sulfide minerals dissolution in various acids (e.g. H2SO4, HCl, HNO3) in temperature range. My question is how to plot pe on y-axis versus reaction temperature on x-axis. Actually, I need a statement in USER_GRAPH block for graph_y ... pe ...

Maybe it is posted but I couldn't find any topic about that.

Regards

[dlparkhurst]

pe can be plotted as -LA("e-"). Alternatively, you can plot SI("O2(g)") [fugacity of O2(g)] as a measure of redox condition.