We investigate the origin of the kink effect (KE) in the relativistic mean field theory by transforming the single-particle Dirac equation into a Schrödinger-like equation. It is found that relativistic self-consistent effects as well as contributions coming from the ρ meson determine the actual structure of the KE. However, the spin-orbit force generated by the ρ meson has no significant influence on the KE. On the other hand, in Refs. [3,4] the anomalous behavior of the charge radii of these isotopes has been studied in the relativistic mean field approach (RMFA). The Pb data are well described by this type of theory. In Ref.[4] it is supposed that the success of the RMFA in reproducing the KE is achieved due to the weak isovector dependence of the spin-orbit force generated by the RMFA.In Ref.[5] the problem of isotopic shifts has been investigated with exhaustive detail, both in the SHF and RMFA. The authors of Ref.[5] have proposed a new Skyrme functional SkI4 with a more general structure of the spin-orbit force than that of the standard Skyrme functional. They have also succeeded to reproduce the KE, but at the price of introducing an extra parameter for the two body spin-orbit potential.As the authors of Refs. [3,4] state, the standard RMFA provides an excellent description of the anomalous kink in the isotopic shifts about 208 P b without using extra fitting parameters. In the present paper we shall focus our interest into the problem of the origin of the KE in the RMFA paying special attention to the role of the isovector ρ meson in the kink structure.In the RMFA [6], the relativistic single-nucleon wave function ψ( r), for a nucleon with the rest mass M, satisfies the Dirac equationwhere E = M + ε and ε is the binding energy eigenvalue while S( r) and V ( r) are the scalar and vector potentials, respectively, defined as follows: where the central potential V cent and the spin-orbit potential V SO are given by: For simplicity, we shall mainly concentrate our discussion in the three Pb isotopes 206 P b, 208 P b and 210 P b. The most stable configuration of the 210 P b isotope in the two nonlinear models considered contains two neutrons in the 1i 11/2 states. In the L model, the two configurations in which two neutrons fill up the 1i 11/2 or 2g 9/2 single-particle states are almost degenerate (we shall refer to these two configurations as IC and GC, respectively).We approach the 210 P b ground state by the IC in both linear and nonlinear models, although the BCS formalism would be more appropriate in the first case. Nevertheless, we expect that our conclusions are also valid for the L model, at least, at a qualitative level.The self-consistent results for the charge radii corresponding to the exact models are summarized in Table I, case a, and in Fig. 1 (the r c are normalized in this figure to the 3