We provide lower bounds for the sum of the negative eigenvalues of the operator |σ · pA| 2s − Cs/|x| 2s + V in three dimensions, where s ∈ (0, 1], covering the interesting physical cases s = 1 and s = 1/2. Here σ is the vector of Pauli matrices, pA = p − A, with p = −i∇ the three-dimensional momentum operator and A a given magnetic vector potential, and Cs is the critical Hardy constant, that is, the optimal constant in the Hardy inequality |p| 2s ≥ Cs/|x| 2s . If spin is neglected, results of this type are known in the literature as Hardy-Lieb-Thirring inequalities, which bound the sum of negative eigenvalues from below by −Ms V 1+3/(2s) − , for a positive constant Ms. The inclusion of magnetic fields in this case follows from the non-magnetic case by diamagnetism. The addition of spin, however, offers extra challenges that make the result more elusive. It is the purpose of this article to resolve this problem by providing simple bounds for the sum of the negative eigenvalues of the operator in question. In particular, for 1/2 ≤ s ≤ 1 we are able to express the bound purely in terms of the magnetic field energy B 2 2 and integrals of powers of the negative part of V .