We show that the integral cross sections for state-to-state quantum scattering of cold molecules in a magnetic field can be efficiently computed using the total angular momentum representation despite the presence of unphysical Zeeman states in the eigenspectrum of the asymptotic Hamiltonian. We demonstrate that the unphysical states arise due to the incompleteness of the space-fixed total angular momentum basis caused by using a fixed cutoff value Jmax for the total angular momentum of the collision complex J. As a result, certain orbital angular momentum (l) basis states lack the full range of J values required by the angular momentum addition rules, resulting in the appearance of unphysical states. We find that by augmenting the basis with a full range of J-states for every l, it is possible to completely eliminate the unphysical states from quantum scattering calculations on molecular collisions in external magnetic fields. To illustrate the procedure, we use the augmented basis sets to calculate the state-to-state cross sections for rotational and spin relaxation in cold collisions of 40CaH(X2Σ+, v = 0, N = 1, MN = 1, MS = 1/2) molecules with 4He atoms in a magnetic field. We find excellent agreement with benchmark calculations, validating our proposed procedure. We find that N-conserving spin relaxation from the highest-energy to the lowest-energy Zeeman state of the N = 1 manifold, |1112〉→|1-1-12〉 is nearly completely suppressed due to the lack of spin-rotation coupling between the fully spin-stretched Zeeman states. Our results demonstrate the possibility of rigorous, computationally efficient, and unphysical state-free quantum calculations on cold molecular collisions and on near-threshold energy levels of strongly anisotropic atom-molecule collision complexes in an external magnetic field.