晶体电荷密度计算