# Something special about energy eigenstates when it comes to time evolution?

The reason as to why the energy eigenstates are special is because of the Schrodinger's equation $$i \hbar \frac{\partial}{\partial t} \Psi(x,t) = \hat{H} \Psi(x,t)$$ If $\hat{H}$ is time independent, the solution of the equation is given by $$\Psi(x,t)= e^{-\frac{i}{\hbar}\hat{H} t} \, \Psi(x,0)$$

The fact that energy eigenstates are special with respect to the time evolution is precisely due to the fact that the Hamiltonian controls the time evolution of the system as given by the above equation.

If you expand your initial state $\Psi(x,0)$ in terms of the eigenstates of the Hamiltonian, then the action of the exponential operator on those states is easily given. Suppose $\psi_n(x)$ represents a state of energy $E_n$ and suppose you can expand your initial state as $$\Psi(x,0) = \sum_n c_n \psi_n(x)$$ where $c_n$ are constants (independent of time). Then the full solution to the Schrodinger's equation is $$\Psi(x,t) = \sum_n c_n e^{-\frac{i}{\hbar} E_n t} \psi_n(x)$$ Now because $c_n$ are independent of time, the probability of the system being in a particular energy eigenstate is same for all time. You can see this by calculating $|\langle \psi_m | \Psi(x,t) \rangle|^2=|c_m|^2$

The reason why for other operators this may not be true for a random operator, say $\hat{B}$ is because the eigenstates of this operator may not evolve independently. Upon time evolution, they may mix into each other. In other words, if you expand a random initial state in terms of the eigenstates of the operator $\hat{B}$ and then look at the projection onto an eigenstate after some time evolution, the coefficient will change. With a little work you can figure out that the condition that the coefficients remain time independent is that the operator $\hat{B}$ commutes with the hamiltonian, i.e., $[\hat{H},\hat{B}]=0$.

P.S - As pointed out in the comments by @ACuriousMind, although $[\hat{H},\hat{B}]=0$ is necessary, it may not be sufficient to ensure that the coefficients are time independent. Ultimately what you need is a common eigenbasis for the operators $\hat{B}$ and $\hat{H}$ to ensure that the time evolution does not mix the states.

An energy eigenstate is an eigenstate of the Hamiltonian $H$. If the Hamiltonian is not time-dependent, then the time evolution operator is $U(t) = \mathrm{e}^{\mathrm{i}Ht}$, so every energy eigenstate is also an eigenstate of the time evolution operator, meaning it remains the same state under time evolution.

No other operator exhibits this peculiar feature for its eigenstates because the Hamiltonian is special in being connected to time evolution through the Schrödinger equation.