If two operators commute, do they have "a mutual set of eigenfunctions", or "the same set of eigenfunctions"? My quantum chemistry book uses these as if they are interchangeable, but they do not seem to be the same in a very significant way.
One direct consequence of my confusion with respect to this arises when considering angular momentum operators and the fact that: $$[\hat{L}^2, \hat{L}_x] = [\hat{L}^2, \hat{L}_y] = [\hat{L}^2, \hat{L}_z] = 0$$ which implies that $\hat{L}^2$ shares a mutual set of eigenfunctions with $\hat{L}_x, \hat{L}_y,$ and $\hat{L}_z$. However, the spherical harmonics (which I thought were the only eigenfunctions of $\hat{L}^2$) are only eigenfunctions of $\hat{L}_z$ and $\hat{L}^2$ (when considering these 4 operators)! Thus what are the eigenfunctions that $\hat{L}^2$ shares with $\hat{L}_x$ and $\hat{L}_y$ since we know there must be some from the commutation relation!? (I understand that we can redefine which axis is x, y, and z, but my point is that only one of the three axes can have its angular momentum component operator have the spherical harmonics as eigenfunctions regardless of how you define your axes).
$\hat{L}_x$, $\hat{L}_y$, and $\hat{L}_z$ do not commute with each other, yet they all three commute with a fourth common operator as already mentioned, $\hat{L}^2$! It doesn't make sense to me how it's possible for $A$ to commute with $B$ and $B$ to commute with $C$ yet $A$ not to commute with $C$.
If two operators don't commute, can they still share, for example, 1 eigenfunction, or must they not share any eigenfunctions at all?
Any answer that doesn't presuppose any extensive math background beyond differential equations and basic linear algebra would be greatly appreciated! I'm taking first semester quantum chemistry. Thanks!