I have been trying to show that $$\|M+N\|_1 \le \|M\|_1 + \|N\|_1$$ However, I seem to be missing some fundamental property of either how the trace or the square root acts on these sums of matrices, or how the Hilbert-Schmidt inner product can be used. I can expand it easily, getting \begin{align*} \|M+N\|_1 &= {\rm tr}|M+N|\\ &={\rm tr}\Big[\sqrt{(M+N)(M+N)^\dagger}\Big]\\ &={\rm tr}\Big[\sqrt{MM^{\dagger}+MN^{\dagger}+NM^{\dagger}+NN^{\dagger}}\Big] \end{align*}
Now obviously you can't just square both sides, as the square operations doesn't distribute over the trace function, nor does the square root over the sum, or this would be trivial using Hilbert-Schmidt inner product. So what am I missing here?