Grand partition function