Выделю в отдельную тему.
Так как так и не нашёл никакой литературы, начал разбираться по статье, которую нашёл. Ссылка вот:
http://tddft.org/TDDFT2004/PracticalSessions/papers/kohn_pr123.1242_1.pdfНа первой же странице заголовок Cyclotron Resonance. Он пишет уравнение на операторы:
![$$\frac{d\hat{\mathbf P}}{dt}=\frac{i}{\hbar}\left[\hat{H}\ \hat{\mathbf P}\right].$$ $$\frac{d\hat{\mathbf P}}{dt}=\frac{i}{\hbar}\left[\hat{H}\ \hat{\mathbf P}\right].$$](https://dxdy-01.korotkov.co.uk/f/8/8/0/880b29a5f7a63523e6ebfc6c5f4b7b9282.png)
Начинаю покомпонентно коммутировать.
Значит гамильтониан

где последний член электрон-электронное взаимодействие. Автор выбрал следующий импульс (в калибровке Ландау):

Над гамильтонианом ставлю шляпу, над магнитным полем - нет. К непониманию привести не должно. Надо прокоммутировать гамильтониан и импульс. Очевидно, что коммутатор с потенциальной энергией будет равен нулю, поэтому пот. энергию не имеет смысла учитывать при коммутации. Т.к.

и

коммутируют, то имеет смысл рассматривать только коммутацию импульса с "удлиннёной" y-компонентой квадрата импульса в гамильтониане. Поэтому, как я понимаю надо лишь прокоммутировать квадрат импульса (только удлинённую у-компоненту), который стоит в гамильтониане, и, собственно, сам импульс. Однако, когда я это проделал, я получил вовсе не то, что получено в статье под номером

. Вот мои расчёты (знак суммы с индексом i временно опустил: вряд ли это на что-то повлияет). Так же, кстати, очевидно, что смысл имеет коммутировать только x-компоненту импульса, потому что коммутаторы с у и z компонентой будут равны нулю.
![$$[p^2\ p_x]\psi=\left(\left(\hat{p}_y +\frac{e}{c}Hx\right)^2\ \hat{p}_x -\hat{p}_x\ \left(\hat{p}_y +\frac{e}{c}Hx\right)^2\right)\psi=\text{раскрыл и подставил импульс}=$$ $$[p^2\ p_x]\psi=\left(\left(\hat{p}_y +\frac{e}{c}Hx\right)^2\ \hat{p}_x -\hat{p}_x\ \left(\hat{p}_y +\frac{e}{c}Hx\right)^2\right)\psi=\text{раскрыл и подставил импульс}=$$](https://dxdy-02.korotkov.co.uk/f/5/2/6/526400548d3edbf232219b36fd9f4a3482.png)

Что (после домножения на

, конечно) не приводится к виду уравнения

. Как быть? Не пойму где ошибка.