For an operator which is a finite sum of products of finitely many Toeplitz operators on the harmonic Bergman space over the half-space, we study the problem: Does the boundary vanishing property of the Berezin transform imply compactness? This is motivated by the Axler-Zheng theorem for analytic Bergman spaces, but the answer would not be yes in general, because the Berezin transform annihilates the commutator of any pair of Toeplitz operators. Nevertheless we show that the answer is yes for certain subclasses of operators. In order to do so, we first find a sufficient condition on general operators and use it to reduce the problem to whether the Berezin transform is one-to-one on related subclasses.