Population-Level Inference for Galaxy Properties from Broadband Photometry with Neural Density Estimation