add a line for population analysis to get the charge distribution