Bose-Einstein grand partition function