Therefore, we can choose an arbitrary value for C 2 and get the corresponding value for By putting these results beck into Eq. ( ) and using the perturbation expansion theory and performing the integrals, the grand potential is obtained as follows. Note that, our perturbation parameter is l B /D which is an small dimensionless parameter. We can expand the grand potential in terms of the perturbation parameter and keep the leading order term as follows,