Основное состояние задается волновой функцией
\[\Psi = e^{i\varphi} \Psi_g.\]Значит вектор состояния можно написать из определения
Гамильтониан должен действовать на вектор состояния также, как он действует на волновую функцию:
\[ \hat{\mathcal{H}} \Psi = \hat{\mathcal{H}} (\alpha \Psi_e + \beta \Psi_g) = \alpha E_e \Psi_e + \beta E_g \Psi_g.\]Отсюда несложно написать гамильтониан в виде матрицы:
\[ \hat{\mathcal{H}}\begin{pmatrix} \alpha \\ \beta \end{pmatrix} = \begin{pmatrix} \alpha E_e \\ \beta E_g \end{pmatrix}.\]
Перепишем уравнение Шредингера покомпонентно:
\[i\hbar \begin{pmatrix} \dot{\alpha} \\ \dot{\beta} \end{pmatrix} = \begin{pmatrix}
E_e \alpha \\ E_g \beta \end{pmatrix}.\]
Получаем систему независимых уравнений на $\alpha$ и $\beta$, которые имеют решения в виде комплексных экспонент:
\[ \alpha(t) = \alpha(0) e^{-\frac{iE_et}{\hbar}}, \quad \beta(t) = \beta(0) e^{-\frac{iE_gt}{\hbar}}.\]
Выполним непосредственное вычисление
\[ \langle x(t) \rangle = \frac{1}{2} \begin{pmatrix} e^{\frac{iE_et}{\hbar}} & e^{\frac{iE_gt}{\hbar}} \end{pmatrix} \begin{pmatrix} 0 & a \\ a & 0 \end{pmatrix} \begin{pmatrix} e^{-\frac{iE_et}{\hbar}} \\ e^{-\frac{iE_gt}{\hbar}} \end{pmatrix} = \frac{a}{2} \left( e^{\frac{it}{\hbar} (E_e - E_g) } + e^{-\frac{it}{\hbar} (E_e - E_g) } \right) = a \cos \frac{(E_e-E_g)t}{\hbar} \]