When studying the hydrogen atom, why do we seek simultaneous eigenfunctions of $\hat{L}^2$, $\hat{L}_z$, and $\hat{H}$? When solving the Schrödinger equation for the hydrogen atom, textbooks invariably work in a more constraint situation, whereby not only an eigenfunction for the Hamiltonian operator $\hat{H}$ is sought, but one which is simultaneously an eigenfunction for $\hat{L}^2$ and $\hat{L}_z$. My question is why we do this?
A similar question has been asked here, but the answers are unsatisfactory. Yes, I understand we can do it. Yes, I understand that we have lots of freedom in our choice of $\psi(\vec{x})$ if we merely solve for $\hat{H}$. But I want to know why this is the right way to proceed. As far as I understand, it is perfectly physically acceptable for a wave function to not be an eigenfunction of some operator, so why must the wave function for a hydrogen atom be an eigenfunction for $\hat{L}^2$ and $\hat{L}_z$?
 A: I'd like to elaborate some more on my comment. As you probably know, quantum mechanics take place in a Hilbert space $\mathcal H$, and every physical state is represented by a unit vector in $\mathcal H$. The time evolution of a given state $\vert\psi\rangle$ is given by the Schrödinger equation $i\hbar\vert\dot\psi\rangle=H\vert\psi\rangle$. Using some math we can show that if $\vert\psi(t=0)\rangle$ is an eigenstate of the Hamiltonian $H$, solving this equation becomes particularly easy (just multiply $\vert\psi(t=0)\rangle$ by an appropriate phase factor).  And if we can express any arbitrary state as a linear combination of eigenstates, it's still pretty easy (give each term a phase factor of its own). So solving the Schrödinger equation reduces to finding a basis of the Hilbert space made up entirely of eigenstates of the Hamiltonian (because then we can express any arbitrary state as a linear combination of eigenstates, and thus solve the equation as described above).
Now the main part of a lecture about the hydrogen atom will consist of finding this basis in different scenarios (with or without spin, with strong/weak/no electromagnetic fields, etc.). And it turns out that the Hamiltonian of the hydrogen atom is degenerate, so we have some free choices when looking for a basis. And it turns out that very conveniently, we have the choice to make the basis states into eigenstates of not only the Hamiltonian, but also additional physically relevant operators: $L^2,L_z,S^2$ and $S_z$.
This does not mean that the basis basis states are the only physically allowed states. Just that all physically allowed states can be expressed as a linear combination of these eigenstates. It's actually quite unlikely to find an atom in any of those basis eigenstates. For instance, the electron might have a defined energy quantum number $n=2$, defined absolute angular momentum number $l=1$, but the direction of the angular momentum might not be along our arbitrarily chosen $z$-axis, so the electron has no defined quantum number $m_l$. So maybe it's in the state $\vert\psi\rangle=\frac{1}{\sqrt2}(\vert n=2,l=1,m_l=1\rangle+\vert n=2,l=1,m_l=-1\rangle)$, which is not one of our basis states. But we can still use the time evolution of the basis states to calculate the time evolution of this other state, too.
There are also different choices for a basis, and every physically possible state can still be written as a linear combination of those. But the one that is found in the standard literature is the most convenient when trying to solve the Schrödinger equation.
A: It is very desirable for $\psi$ to be an eigenfunction of as many operators as possible
In fact, we would probably like it to be an eigenfunction of all angular momenta, but they do not commute, so we can't make it be an eigenfunction of all momenta, so we choose only one (usually $L_z$)
There are many reasons, but the main ones are:

*

*They are measurable

The Hamiltonian is the energy operator (roughly speaking). If a function is an eigenstate of the Hamiltonian, that means that $\phi_n$ has a defined energy. Energy is a quantity which is easily measurable, so we choose the Hamiltonian to be an important operator in our complete system of commuting operators (CSCO)
In the same way, angular momenta are easily measurable using magnetic fields.


*The fact that a function is an eigenstate of an operator allows us to "label" it with quantum numbers. We can set a state to be $|n\ l\ m_l\ s \ m_s\rangle$ because of that. You cannot have a quantum number if the function is not an eigenstate, because it wouldn't have an eigenvalue to label it with.

So having "something easily measurable to label with" is a good idea. Saying $n=1$ is good because it is easily measurable.
Plus, it is a really basic concept in physics, and we have a lot of intuition on it. If we say that the energy level is the first level, we quickly get an idea of how the electron is. This would not happen with weirder magnitudes.


*The Hamiltonian is more than that

The Hamiltonian, besides the energy, is also responsible of the time evoluition. As it is involved in the Schrödinger's equation, the Hamiltonian rules the time evolution. If something commutes with the Hamiltonian, then that quantity is conserved over time. So if $[H, L_z]=O$, then the value you measure for $L_z$ is conserved over time. So a state with $m_s=+1$, for example, will keep that value over time. That's why the Hamiltonian is important.
The angular momenta are important as well, because things that commute with the angular momentum are invariant under rotations, which is useful, because it tells us if things are symmetric or not. You know that symmetry plays an important role when simplifying problems.


*Commutation itself

The very fact that two operators commute is probably always a good thing. Commutation means that you can measure $H$ and then $L_z$ and vice-versa, and the result is the same. Translation: if you measure energy, it doesn't disrupt the system and you can measure $L_z$ afterwards. If they do not commute, measuring one changes the state for the next measurement.
So, you want them to commute, because you want to measure both things of the same state. It happens to be that "knowing" $L_z$ is not enough to determine the state, there are many states with the same value of $m_l$. We need more information.
However, given a certain energy, with a certain momentum, and certain spin, and so on, the state is unambiguously determined.
That's why we need a complete set of operators. And they must commute so that we can measure all of them without disrupting the system under measurement. That's why we seek for a CSCO.
And, of course, if you need 5 operators to determine your system, you'd better choose operators whose magitudes are easily measurable and have a good physical meaning.
And there are probably more reasons, but I can't remember them all now. Feel free to complete my answer in comments.
A: The question that you have shared does deal with this point, but I thought I'd elaborate on it a little more. When we solve a quantum mechanical system, we would like to work with an eigenbasis that is labelled by quantum numbers that completely specify the state of the system. Essentially, this comes down to finding something called a Complete Set of Commuting Observables.
The wikipedia article linked gives a good introduction, but basically a "CSCO" is the largest set of operators that can be formed such that we can speak of the eigenvalues of all of these operators simultaneously. (For example, $\hat{x}$ and $\hat{p}$ would never be in a CSCO since there exists no eigenbasis that is simultaneously an eigenbasis of $\hat{x}$ and $\hat{p}$.)
In the case of the hydrogen atom, there is a large amount of degeneracy. Solving the Schrodinger equation, it can be shown that for a state of definite energy $n$, the degeneracy is $n^2$. (I'll come back to this further down.) But all of these $n^2$ states are not equivalent. While they all do have the same energy, they differ in their values of total and azimuthal angular momentum. Finding an eigenbasis that is a simultaneous eigenbasis of all of these operators would allow us to label each of these degenerate states with two more numbers, which makes them "unique".
But, you may ask, how do we know that these three operators are sufficient to form the  CSCO? To my understanding, there's no way for us to know, apart from experiment. It turns out that $\{\hat{H}, \hat{L}^2, \hat{L}_z\}$ don't form a CSCO on their own! To completely specify a state of an electron in the hydrogen atom one is also required to specify the Spin of the electron, and so the true CSCO is $\{\hat{H}, \hat{L}^2, \hat{L}_z, \hat{S}_z\}$. (This is why the true degeneracy of the hydrogen atom is $2n^2$, as there are two values of spin possible per state.)
So as I see it, the states of the Hydrogen atom don't have to be eigenstates of all three operators. But if we would like each physically distinguishable state to be represented by a unique vector in the Hilbert Space, then we need to identify each eigenvector by the largest set of eigenvalues that uniquely specifies it. Ignoring spin, these eigenvalues are $|n l m\rangle$.
A: " so why must the wave function for a hydrogen atom be an eigenfunction for $\hat{L}^2$ and $\hat{L}_z$?" They must not, but since the $\hat{H}$, $\hat{L}^2$, $\hat{L}_z$ do commute with each other, they can have some mutual eigenfunctions.You do not have to choose these specific mutual eigenfunctions, you can choose some others, but I think it is obviously very convenient to work with such functions that are simultaneously eigenfunctions of all these three operators(and this is why people usually choose these ones).
A: The nuclear potential is rotationally symmetric so all eigenfunctions must be simultaneously eigenfunctions of angular momentum.
