add finite mu_B into the calculations