We study the attractive interactions between rod-like charged polymers in solution that appear in the presence of multi-valence counterions. The counterions condensed to the rods exhibit both a strong transversal polarization and a longitudinal crystalline arrangement. At short distances between the rods, the fraction of condensed counterions increases, and the majority of these occupy the region between the rods, where they minimize their repulsive interactions by arranging themselves into packing structures. The attractive interaction is strongest for multivalent counterions. Our model takes into account the hard-core volume of the condensed counterions and their angular distribution around the rods. The hard core constraint strongly suppresses longitudinal charge fluctuations.PACS numbers: 61.20. Qg, 61.25.Hq, 87.15.Aa Strongly charged polymers precipitate from a dilute solution into compact structures when high-valence counterions (oppositely charged particles) are added to the solution [1][2][3][4][5][6]. The counterions experience strong electrostatic attractions to the backbone of the chains, and a finite fraction of them "condense", i.e., are found within short distance from the chains [7]. Counterions are more attracted to compact chains or aggregates of rod-like chains. This creates the possibility of a transition from single chains with a small number of condensed counterions to almost neutral aggregates of chains, or even mono-molecular collapse in the case of flexible polymers. These aggregates are stable only when the internal arrangement of the counterions within them, provides a strong enough cohesive energy.In this letter we study the attraction between two rodlike polyelectrolytes. We show that it is essential to include the size and angular degrees of freedom (around the rods) of the counterions as well as the discrete nature of the charge along the polyelectrolytes to find the origin and strength of the counterion mediated attraction. Our work suggests that these factors are also crucial in determining the collapse of flexible and semi-flexible polyelectrolytes recently studied in Refs. [8][9][10][11][12]. Experimental observations show that the size of the precipitating particles is indeed a relevant parameter in the problem [1,6].It has been argued that longitudinal charge fluctuations resulting from the thermal motion of point counterions induce attractions between rod-like polyelelctrolytes [13,14] and induces "buckling" of semiflexible polyelelctrolytes [11,12]. Here, we show that such charge fluctuation are suppressed when the hard core volume of monomers and counterions are taken into account. Instead, we find that the counterions arrangement around the rods create a non-zero transversal polarization as the distance between the chains decreases. At very short distances between the rods we find strong longitudinal correlations but only at very short wavelengths, implying a crystalline state along the rod, reinforcing, for the case of multivalent counterions, the attractive interactio...