To treat correlations resulting from electrostatic interactions between charged species on the same footing as those resulting from hard sphere excluded volume interactions Mean Spherical Approximation was employed to solve the Ornstein-Zernike equation with respect to electrostatic direct correlation function. Taylor expansion of electrostatic free energy was cut after the second order. Then the electrostatic correlation component of the free energy (F ex corr ) is