Variational reduced density matrix method in the doubly-occupied configuration interaction space using four-particle N-representability conditions: Application to the XXZ model of quantum magnetism | Publicación