function [pH] = pHfunc_Chain(H_plus)

pH = -log10(abs(H_plus));

end