  Nuclear spin levels play an important role in understanding magnetization
dynamics and implementation and control of quantum bits in lanthanide-based
single-molecule magnets. We investigate the hyperfine and nuclear quadrupole
interactions for $^{161}$Dy and $^{163}$Dy nucleus in anionic DyPc$_2$
(Pc=phthalocyanine) single-molecule magnets, using multiconfigurational
ab-initio methods (beyond density-functional theory) including spin-orbit
interaction. The two isotopes of Dy are chosen because the others have zero
nuclear spin. Both isotopes have the nuclear spin $I=5/2$, although the
magnitude and sign of the nuclear magnetic moment differ from each other. The
large energy gap between the electronic ground and first-excited Kramers
doublets, allows us to map the microscopic hyperfine and quadrupole interaction
Hamiltonian onto an effective Hamiltonian with an electronic pseudo-spin
$S_{\rm eff}=1/2$ that corresponds to the ground Kramers doublet. Our ab-initio
calculations show that the coupling between the nuclear spin and electronic
orbital angular momentum contributes the most to the hyperfine interaction and
that both the hyperfine and nuclear quadrupole interactions for $^{161}$Dy and
$^{163}$Dy nucleus are much smaller than those for $^{159}$Tb nucleus in
TbPc$_2$ single-molecule magnets. The calculated separations of the
electronic-nuclear levels are comparable to experimental data reported for
$^{163}$DyPc$_2$. We demonstrate that hyperfine interaction for Dy Kramers ion
leads to tunnel splitting (or quantum tunneling of magnetization) at zero
field. This effect does not occur for TbPc$_2$ single-molecule magnets. The
magnetic field values of the avoided level crossings for $^{161}$DyPc$_2$ and
$^{163}$DyPc$_2$ are found to be noticeably different, which can be observed
from experiment.
