High precision calculation of the isotope shift of the 3 2 S 1=2 -2 2 S 1=2 transition in lithium is presented. The wave function and matrix elements of relativistic operators are obtained by using recursion relations. Apart from the relativistic contribution, we obtain the nuclear polarizability correction for 11 Li. The resulting difference of the squared charge radii 11 The recent advances in precise spectroscopy of atomic systems make possible the determination of nuclear charge radii from isotope shift measurements [1][2][3][4]. In spite of the fact that the nuclear size is 5 orders of magnitude smaller than the atomic size, the achieved experimental precision for transition frequencies allows one to determine nuclear charge radii much more accurately than from electron scattering measurements. This, however, requires the theoretical calculations to be performed with high relative precision, for example, at least 10 ÿ6 for lithium isotopes. The accuracy of theoretical predictions achieved for hydrogen [5] leads to the determination of the proton charge radius [5,6], which is far more accurate than from the electron scattering measurements. A similar experimental accuracy achieved for deuterium allowed one to determine very accurately the deuteron radius [1]. Surprisingly, the atomic measurements lead to a slightly different value from the electron scattering determination, which stimulated a reanalysis of electron scattering data. What makes the deuteron different from other typical nuclei is its low binding energy of about 2.226 MeV. Such a soft nucleus is distorted by a surrounding electron, which leads to a shift in atomic transition frequencies. Even smaller is the twoneutron separation energy in 11 Li [7], which indicates the possible significance of nuclear structure effects on the isotope shift.In this Letter, we present significantly improved calculations of finite nuclear mass contributions to the isotope shift in the lithium 3 2 S 1=2 -2 2 S 1=2 transition. Such calculations have already been performed by Yan and Drake in Refs. [8][9][10] and were used to determine nuclear charge radii of various lithium isotopes on the basis of recent isotope shift measurements [3,4]. Our result for the relativistic recoil correction is about 10 times smaller than that reported in Refs. [8,9]. Apart from the known nonrelativistic, leading relativistic, and QED corrections, we include higher order recoil corrections and the nuclear polarizability effect E pol , the last being significant for the 11 Li nucleus. Finally, we use our combined results to calculate new nuclear charge radii for the lithium isotopes, taking the experimental isotope shift results from Refs. [3,4].Let us denote ÿ =M. The expansion of an energy level in the fine structure constant iswhere E fs is the nuclear finite size correction, and m, M, are the electron, nucleus, and the reduced mass, respectively. In the following, we obtain these expansion coefficients. The nonrelativistic Hamiltonian of the lithium atom in atomic units is