My first understanding of 'van Vleck paramagnetism' refers to open shell configurations short of half filling by one electron. In LS coupling we then have L=S and by the third Hund's rule we obtain a singlet ground state. There is nevertheless a finite paramagnetic susceptibility because of the Zeeman Hamiltonian couples singlet and triplet states (by nonzero off-diagonal elements).
Having started to deal with crystal field states I have learned that more generally, non-diagonal Terms of Zeeman-Hamiltonians are being called 'van Vleck terms'. In case they couple non-degenerate states this produces a second order contribution to the paramagnetic moment and susceptibility. As I see it, this is what constitutes the so-called 'single ion anisotropy'.
However, I have come across abstracts & papers referring to e.g. Tm as a 'van Vleck ion'. Does this imply a singlet ground state (how would Tm do that ?) or more generally express the notion of an 'ion with single ion anisotropy' (due to van Vleck terms)?