Negative Binomial LDS via Polya-Gamma Augmentation for Neural Spike Count Modeling