We revisit the problem of the kink in the charge radius shift of neutron-rich even isotopes near the N = 126 shell closure. We show that the ability of a Skyrme force to reproduce the isotope shift is determined by the occupation of the neutron 1i 11/2 orbital beyond N = 126 and the corresponding change it causes to deeply-bound protons orbitals with a principal quantum number of 1. Given the observed position of the single-particle energies, one must either ensure occupation is allowed through correlations, or not demand that the single-particle energies agree with experimental values at the mean-field level.PACS numbers: 21.10. Ft, 21.60.Jz, 27.80.+w The evolution of charge radii across the isotope chart is one of the most basic nuclear structure observables, as it provides a particularly useful characterisation of the proton distribution that can be accessed by a variety of experiments [1]. We define the charge radius isotope shift as the differences between the mean squared charge radius, r 2 ch , of a series of isotopes and that of a given reference isotope (208 Pb for lead and 210 Po for polonium). Theoretically, charge radii have been traditionally studied within the droplet model [2], which captures qualitatively their evolution across the nuclear chart. Yet, in some specific cases, quantum shell effects dominate the density distribution and provide a departure from smooth systematic trends. Perhaps the most well-known example of these abrupt changes is the kink in the isotope shift of even lead isotopes as one passes through the N = 126 magic number.A summary of the experimental results [1, 3] around this neutron shell closure is given in Figures 1 and 2 for lead and polonium, respectively. For lead, the fitted solid line is of the form:Linear regression givesi.e. the slope of the shift is observed to double at A = 208 to a very good approximation. This abrupt change in change radius cannot be explained within the droplet model [2]. Recent experimental work using laser spectroscopy techniques have identified a similar kink in neutron-rich polonium isotopes above the Z = 82 shell closure [3] (see also Fig. 2) in agreement with older radon and radium data [4]. Moreover, the details of the changes of proton and neutron radii in isotopes around 208 Pb have also gained renewed interested in the form of the neutron skin, which correlates strongly with nuclear matter properties [5,6]. On the neutron deficient side, the onset of deformation is also probed by measurements of isotope shifts [7]. In this letter, we propose a new mechanism to explain the existence of the kink in lead isotope shifts using density functional calculations supplemented by pairing effects.On the theoretical side, mean-field models, or equivalently density functional theories have been widely applied to the systematic study of all observed and hypothesized nuclear isotopes [8]. Both the Skyrme-HartreeFock (SHF) and the Relativistic Mean-Field (RMF) approaches are able to give broadly good descriptions of many nuclear properties acros...