This work presents a first time accurate calculation of the magnetic dipole hyperfine structure constants for the ground state and some low-lying excited states of Pb + . By comparing different levels of approximation with experimental results, we demonstrate the importance of correlation effects which reach beyond lower order relativistic many body perturbation theory. Employing relativistic coupled-cluster theory we obtain a quantitative understanding of the core-polarization and correlation effects inherent in this system and observe completely different trends compared to Ba + .Coupled-cluster theory has been used to study a wide range of many-body systems [1]. Although the non-relativistic version of this theory has been very successfully applied to a variety of light atoms and molecules [2], its extension to the relativistic regime is rather recent [3,4]. There have been relatively few theoretical studies of properties of heavy atomic systems based on the relativistic coupled-cluster (RCC) theory. Pb + (Z=82) is the heaviest atomic ion that has been trapped and cooled so far [5,6]. The magnetic dipole hyperfine constants have been measured for the 6p 2 P 1/2 and 6p 2 P 3/2 states of this ion [7] and these data can be compared with calculations of the corresponding quantities using RCC theory. Such comparisons would indeed constitute an important test of this theory. The non-linear RCC in the singles and doubles approximation with partial triples added in some cases has yielded results to an accuracy of about one percent for atoms and ions with a single s valence electron [8,9,10]. However, the correlation effects in Pb + are expected to be much stronger as it has a 6p valence electron and two 6s electrons in its outermost core orbital.The hyperfine structure constant (A)for the atomic state |JM can be expressed in terms of a reduced expectation value where r j is the radial position of the j th electron, α j is the Dirac matrix and Y10 is a vector spherical harmonic.We have used the RCC theory in to obtain the atomic wavefunctions. As pointed out in our earlier work [12] coupled-cluster theory is equivalent to all order manybody perturbation theory (MBPT). In the open-shell coupled-cluster theory [13,14] the many-body wavefunction for a system with single valence electron can be written aswhere a † v is the creation operator corresponding to a valence orbital 'v' and |Φ 0 is a closed-shell determinantal state built from occupied Dirac-Fock (DF) orbitals. T-and S v -are the closed and open shell excitation operators respectively. In this work both T-and S voperators are truncated beyond double excitations and triple excitations are added on the leading order MBPT level.Explicitly, the T-operator is defined as