In upcoming work, we will investigate the performance of the methods introduced herein in more detail for both transition metal and lanthanoid complex in order to establish trends and develop an intuitive understanding of the field dependence of NMR shifts in paramagnetic molecules. Furthermore, we plan the implementation of the approaches developed herein beyond model Hamiltonians, e.g. at the level of the ab initio QDPT method.