Vibrational contributions to indirect spin-spin coupling constants calculated via variational anharmonic approaches