Compute the real-space charge (and spin) density of the crystal. Returns Gaussian cube file(s) for visualization plus a compact JSON summary (spin channels, cube names). Useful for bonding, charge transfer, and magnetic density distributions.