function alpha = calc_alpha_z1(params,V,z)
% params: [kf, kr];
R = 8.314;
T = 310;
F = 96485;

% alpha = params(:,1).*exp(1*F*V/R/T);
alpha = bsxfun(@times, params(:,1)',exp(z*F*V/R/T));

end