Make solution class use SOLUTION.H.activity_coefficient to calculate the H+ activity.