AbstractThe general tensorial form of the hyperfine interaction operator in the formalism of second quantization is presented. Both diagonal and off-diagnonal matrix elements of the above-mentioned operator are found using an approach based on a combination of second quantization in the coupled tensorial form, angular momentum theory in three spaces (orbital, spin and quasispin) and a generalised graphical technique. This methodology allows us to account for correlation effects efficiently and, therefore, to study the hyperfine interactions in complex many-electron atoms, those with openf-shells included, in a practical manner. All this will lead us to design an efficient program for large scale calculations of hyperfine structure and isotope shift.