Hyperfine and quadrupole interactions for Dy isotopes in DyPc$_2$ molecules

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.