Crystal Melting, BPS Quivers and Plethystics

We study the refined and unrefined crystal/BPS partition functions of D6-D2-D0 brane bound states for all toric Calabi-Yau threefolds without compact 4-cycles and some non-toric examples. They can be written as products of (generalized) MacMahon functions. We check our expressions and use them as vacuum characters to study the gluings. We then consider the wall crossings and discuss possible crystal descriptions for different chambers. We also express the partition functions in terms of plethystic exponentials. For $\mathbb{C}^3$ and tripled affine quivers, we find their connections to nilpotent Kac polynomials. Similarly, the partition functions of D4-D2-D0 brane bound states can be obtained by replacing the (generalized) MacMahon functions with the inverse of (generalized) Euler functions.