here use either solubility_vec[i] or teble_limits[i] instead of c_aq_max[i]