In this paper, we propose a theoretical method to deduce the charge polarization (CP) and most probable charge for fission fragments for the selected range of mass numbers based on a quantum many-body framework, namely, a constrained Skyrme Hartree–Fock+BCS model. We investigate the CP on the low energy neutron-induced fission of [Formula: see text]U as a representative reaction. We have found that the calculated CP has a noticeable dip that is clearly affected by the spherical shell structure around [Formula: see text]Sn and deformed one around [Formula: see text]144 and a peak for their counterparts, which are common features to those of widely-used Wahl’s systematics, but otherwise, our CP values converge to zero. Especially, we could not see a deep valley and a big hump in the region of symmetric fission at [Formula: see text]112 and 124 present in Wahl’s systematics. The gradient of our CP near the symmetric fission is negative concerning the mass number of fission fragments, which is consistent with other previous theoretical studies; however, that is the opposite of the gradient in Wahl’s systematics. This paper shows the difference between the charge distribution of fission fragments, which has been used practically so far, in comparison with the calculation by microscopic nuclear theory.