Уровни Ландау

Квантование циклотронных орбит

В квантовой механике энергии циклотронных орбит заряженных частиц в однородном магнитном поле квантуются до дискретных значений, поэтому известны как уровни Ландау . Эти уровни вырождены , с числом электронов на уровне, прямо пропорциональным силе приложенного магнитного поля. Он назван в честь советского физика Льва Ландау . [1]

Квантование Ландау способствует магнитной восприимчивости металлов, известной как диамагнетизм Ландау . В сильных магнитных полях квантование Ландау приводит к колебаниям электронных свойств материалов в зависимости от приложенного магнитного поля, известным как эффекты Де Гааза–Ван Альфена и Шубникова–де Гааза .

Квантование Ландау является ключевым элементом в объяснении целочисленного квантового эффекта Холла .

Вывод

Схема циклотронной орбиты частицы со скоростью v , которая является классической траекторией заряженной частицы (в данном случае положительного заряда) в однородном магнитном поле B. Квантование Ландау относится к квантовой заряженной частице в однородном магнитном поле.

Рассмотрим систему невзаимодействующих частиц с зарядом q и спином S, ограниченную областью A = L x L y в плоскости xy . Приложим однородное магнитное поле вдоль оси z . В единицах СИ гамильтониан этой системы (здесь эффекты спина пренебрегаются) равен Здесь — канонический оператор импульса , а — оператор для электромагнитного векторного потенциалапространстве положений ) . Б = ( 0 0 Б ) {\displaystyle \mathbf {B} ={\begin{pmatrix}0\\0\\B\end{pmatrix}}} ЧАС ^ = 1 2 м ( п ^ д А ^ ) 2 . {\displaystyle {\hat {H}}={\frac {1}{2m}}\left({\hat {\mathbf {p} }}-q{\hat {\mathbf {A} }}\right)^{2}.} п ^ {\textstyle {\hat {\mathbf {p} }}} А ^ {\textstyle {\hat {\mathbf {A} }}} А {\textstyle \mathbf {А} } А ^ = А {\textstyle {\hat {\mathbf {A} }}=\mathbf {A} }

Векторный потенциал связан с магнитным полем соотношением Б = × А . {\displaystyle \mathbf {B} =\mathbf {\nabla} \times \mathbf {A}.}

Существует некоторая свобода калибровки в выборе векторного потенциала для данного магнитного поля. Гамильтониан калибровочно инвариантен , что означает, что добавление градиента скалярного поля к A изменяет общую фазу волновой функции на величину, соответствующую скалярному полю. Но физические свойства не зависят от конкретного выбора калибровки.

В калибре Ландау

Из возможных решений для A часто используется калибровочная фиксация, введенная Львом Ландау для заряженных частиц в постоянном магнитном поле. [2]

Когда тогда возможно решение [3] в калибровке Ландау (не путать с калибровкой Ландау ). Б = ( 0 0 Б ) {\displaystyle \mathbf {B} ={\begin{pmatrix}0\\0\\B\end{pmatrix}}} А = ( 0 Б х 0 ) {\displaystyle \mathbf {A} ={\begin{pmatrix}0\\B\cdot x\\0\end{pmatrix}}} Р ξ {\displaystyle R_{\xi }}

В этой калибровке гамильтониан есть Оператор коммутирует с этим гамильтонианом, поскольку оператор отсутствует для этого выбора калибровки. Таким образом, оператор можно заменить его собственным значением . Поскольку не появляется в гамильтониане, а в кинетической энергии появляется только z-импульс, это движение вдоль z-направления является свободным движением. ЧАС ^ = п ^ х 2 2 м + 1 2 м ( п ^ у д Б х ^ ) 2 + п ^ з 2 2 м . {\displaystyle {\hat {H}}={\frac {{\hat {p}}_{x}^{2}}{2m}}+{\frac {1}{2m}}\left({\hat {p}}_{y}-qB{\hat {x}}\right)^{2}+{\frac {{\hat {p}}_{z}^{2}}{2m}}.} п ^ у {\displaystyle {\hat {p}}_{y}} у ^ {\displaystyle {\hat {y}}} п ^ у {\displaystyle {\hat {p}}_{y}} к у {\displaystyle \hbar k_{y}} з ^ {\displaystyle {\hat {z}}}

Гамильтониан можно также записать проще, заметив, что циклотронная частота равна , что дает Это в точности гамильтониан для квантового гармонического осциллятора , за исключением того, что минимум потенциала смещен в координатном пространстве на . ω с = д Б / м {\displaystyle \omega _{c}=qB/m} ЧАС ^ = п ^ х 2 2 м + 1 2 м ω с 2 ( х ^ к у м ω с ) 2 + п ^ з 2 2 м . {\displaystyle {\hat {H}}={\frac {{\hat {p}}_{x}^{2}}{2m}}+{\frac {1}{2}}m\omega _{\rm {c}}^{2}\left({\hat {x}}-{\frac {\hbar k_{y}}{m\omega _{\rm {c}}}}\right)^{2}+{\frac {{\hat {p}}_{z}^{2}}{2m}}.} х 0 = к у / м ω с {\displaystyle x_{0}=\hbar k_{y}/m\omega _{c}}

Чтобы найти энергии, обратите внимание, что трансляция потенциала гармонического осциллятора не влияет на энергии. Таким образом, энергии этой системы идентичны энергиям стандартного квантового гармонического осциллятора , [4] Энергия не зависит от квантового числа , поэтому будет конечное число вырождений (если частица помещена в неограниченное пространство, это вырождение будет соответствовать непрерывной последовательности ). Значение непрерывно, если частица неограничена в направлении z, и дискретно, если частица также ограничена в направлении z. Каждый набор волновых функций с одинаковым значением называется уровнем Ландау . Э н = ω с ( н + 1 2 ) + п з 2 2 м , н 0. {\displaystyle E_{n}=\hbar \omega _{\rm {c}}\left(n+{\frac {1}{2}}\right)+{\frac {p_{z}^{2}}{2m}},\quad n\geq 0.} к у {\displaystyle k_{y}} п у {\displaystyle p_{y}} п з {\displaystyle p_{z}} н {\displaystyle n}

Для волновых функций напомним, что коммутирует с гамильтонианом. Тогда волновая функция разлагается на произведение собственных состояний импульса в направлении и собственных состояний гармонического осциллятора, смещенных на величину в направлении: где . В целом состояние электрона характеризуется квантовыми числами, , и . п ^ у {\displaystyle {\hat {p}}_{y}} у {\displaystyle у} | ϕ н {\displaystyle |\phi _{n}\rangle } х 0 {\displaystyle x_{0}} х {\displaystyle x} Ψ ( х , у , з ) = е я ( к у у + к з з ) ϕ н ( х х 0 ) {\displaystyle \Psi (x,y,z)=e^{i(k_{y}y+k_{z}z)}\phi _{n}(x-x_{0})} k z = p z / {\displaystyle k_{z}=p_{z}/\hbar } n {\displaystyle n} k y {\displaystyle k_{y}} k z {\displaystyle k_{z}}

В симметричной калибровке

Вывод трактовал x и y как асимметричные. Однако, в силу симметрии системы, нет физической величины, которая различала бы эти координаты. Тот же результат можно было бы получить, если бы подходящим образом поменяли местами x и y .

Более адекватным выбором калибра является симметричный калибр, который относится к выбору A ^ = 1 2 B × r ^ = 1 2 ( B y B x 0 ) . {\displaystyle {\hat {\mathbf {A} }}={\frac {1}{2}}\mathbf {B} \times {\hat {\mathbf {r} }}={\frac {1}{2}}{\begin{pmatrix}-By\\Bx\\0\end{pmatrix}}.}

В терминах безразмерных длин и энергий гамильтониан можно выразить как H ^ = 1 2 [ ( i x + y 2 ) 2 + ( i y x 2 ) 2 ] {\displaystyle {\hat {H}}={\frac {1}{2}}\left[\left(-i{\frac {\partial }{\partial x}}+{\frac {y}{2}}\right)^{2}+\left(-i{\frac {\partial }{\partial y}}-{\frac {x}{2}}\right)^{2}\right]}

Правильные единицы измерения можно восстановить, введя множители и . q , , B {\displaystyle q,\hbar ,\mathbf {B} } m {\displaystyle m}

Рассмотрим операторов a ^ = 1 2 [ ( x 2 + x ) i ( y 2 + y ) ] a ^ = 1 2 [ ( x 2 x ) + i ( y 2 y ) ] b ^ = 1 2 [ ( x 2 + x ) + i ( y 2 + y ) ] b ^ = 1 2 [ ( x 2 x ) i ( y 2 y ) ] {\displaystyle {\begin{aligned}{\hat {a}}&={\frac {1}{\sqrt {2}}}\left[\left({\frac {x}{2}}+{\frac {\partial }{\partial x}}\right)-i\left({\frac {y}{2}}+{\frac {\partial }{\partial y}}\right)\right]\\{\hat {a}}^{\dagger }&={\frac {1}{\sqrt {2}}}\left[\left({\frac {x}{2}}-{\frac {\partial }{\partial x}}\right)+i\left({\frac {y}{2}}-{\frac {\partial }{\partial y}}\right)\right]\\{\hat {b}}&={\frac {1}{\sqrt {2}}}\left[\left({\frac {x}{2}}+{\frac {\partial }{\partial x}}\right)+i\left({\frac {y}{2}}+{\frac {\partial }{\partial y}}\right)\right]\\{\hat {b}}^{\dagger }&={\frac {1}{\sqrt {2}}}\left[\left({\frac {x}{2}}-{\frac {\partial }{\partial x}}\right)-i\left({\frac {y}{2}}-{\frac {\partial }{\partial y}}\right)\right]\end{aligned}}}

Эти операторы следуют определенным коммутационным соотношениям [ a ^ , a ^ ] = [ b ^ , b ^ ] = 1. {\displaystyle [{\hat {a}},{\hat {a}}^{\dagger }]=[{\hat {b}},{\hat {b}}^{\dagger }]=1.}

В терминах приведенных выше операторов гамильтониан можно записать как, где мы снова ввели единицы измерения. H ^ = ω c ( a ^ a ^ + 1 2 ) , {\displaystyle {\hat {H}}=\hbar \omega _{\rm {c}}\left({\hat {a}}^{\dagger }{\hat {a}}+{\frac {1}{2}}\right),}

Индекс уровня Ландау является собственным значением оператора . n {\displaystyle n} N ^ = a ^ a ^ {\displaystyle {\hat {N}}={\hat {a}}^{\dagger }{\hat {a}}}

Применение увеличивается на одну единицу при сохранении , тогда как применение одновременно увеличивает и уменьшает на одну единицу. Аналогия с квантовым гармоническим осциллятором дает решения , где и b ^ {\displaystyle {\hat {b}}^{\dagger }} m z {\displaystyle m_{z}} n {\displaystyle n} a ^ {\displaystyle {\hat {a}}^{\dagger }} n {\displaystyle n} m z {\displaystyle m_{z}} H ^ | n , m z = E n | n , m z , {\displaystyle {\hat {H}}|n,m_{z}\rangle =E_{n}|n,m_{z}\rangle ,} E n = ω c ( n + 1 2 ) {\displaystyle E_{n}=\hbar \omega _{\rm {c}}\left(n+{\frac {1}{2}}\right)} | n , m z = ( b ^ ) m z + n ( m z + n ) ! ( a ^ ) n n ! | 0 , 0 . {\displaystyle |n,m_{z}\rangle ={\frac {({\hat {b}}^{\dagger })^{m_{z}+n}}{\sqrt {(m_{z}+n)!}}}{\frac {({\hat {a}}^{\dagger })^{n}}{\sqrt {n!}}}|0,0\rangle .}

Можно убедиться, что указанные выше состояния соответствуют выбору волновых функций, пропорциональных , где . ψ n , m z ( x , y ) = ( w w ¯ 4 ) n w n + m z e | w | 2 / 4 {\displaystyle \psi _{n,m_{z}}(x,y)=\left({\frac {\partial }{\partial w}}-{\frac {\bar {w}}{4}}\right)^{n}w^{n+m_{z}}e^{-|w|^{2}/4}} w = x i y {\displaystyle w=x-iy}

В частности, низший уровень Ландау состоит из произвольных аналитических функций, умноженных на гауссиану, . n = 0 {\displaystyle n=0} ψ ( x , y ) = f ( w ) e | w | 2 / 4 {\displaystyle \psi (x,y)=f(w)e^{-|w|^{2}/4}}

Вырождение уровней Ландау

В калибре Ландау

Эффекты уровней Ландау можно наблюдать только тогда, когда средняя тепловая энергия kT меньше расстояния между энергетическими уровнями , что означает низкие температуры и сильные магнитные поля. k T ω c {\displaystyle kT\ll \hbar \omega _{c}}

Каждый уровень Ландау вырожден из-за второго квантового числа , которое может принимать значения, где — целое число. Допустимые значения дополнительно ограничены условием, что центр силы осциллятора, , должен физически находиться внутри системы, . Это дает следующий диапазон для , k y {\displaystyle k_{y}} k y = 2 π N L y , {\displaystyle k_{y}={\frac {2\pi N}{L_{y}}},} N {\displaystyle N} N {\displaystyle N} x 0 {\displaystyle x_{0}} 0 x 0 < L x {\displaystyle 0\leq x_{0}<L_{x}} N {\displaystyle N} 0 N < m ω c L x L y 2 π . {\displaystyle 0\leq N<{\frac {m\omega _{\rm {c}}L_{x}L_{y}}{2\pi \hbar }}.}

Для частиц с зарядом верхнюю границу можно просто записать как отношение потоков , где — фундаментальный квант магнитного потока , а — поток через систему (с площадью ). q = Z e {\displaystyle q=Ze} N {\displaystyle N} Z B L x L y ( h / e ) = Z Φ Φ 0 , {\displaystyle {\frac {ZBL_{x}L_{y}}{(h/e)}}=Z{\frac {\Phi }{\Phi _{0}}},} Φ 0 = h / e {\displaystyle \Phi _{0}=h/e} Φ = B A {\displaystyle \Phi =BA} A = L x L y {\displaystyle A=L_{x}L_{y}}

Таким образом, для частиц со спином максимальное число частиц на уровень Ландау равно что для электронов (где и ) дает два доступных состояния для каждого кванта потока, который проникает в систему. S {\displaystyle S} D {\displaystyle D} D = Z ( 2 S + 1 ) Φ Φ 0 , {\displaystyle D=Z(2S+1){\frac {\Phi }{\Phi _{0}}},} Z = 1 {\displaystyle Z=1} S = 1 / 2 {\displaystyle S=1/2} D = 2 Φ / Φ 0 {\displaystyle D=2\Phi /\Phi _{0}}

Вышеизложенное дает лишь приблизительное представление об эффектах геометрии конечных размеров. Строго говоря, использование стандартного решения гармонического осциллятора справедливо только для систем, неограниченных в -направлении (бесконечные полосы). Если размер конечен, граничные условия в этом направлении приводят к нестандартным условиям квантования магнитного поля, включающим (в принципе) оба решения уравнения Эрмита. Заполнение этих уровней многими электронами все еще [5] является активной областью исследований. x {\displaystyle x} L x {\displaystyle L_{x}}

В общем, уровни Ландау наблюдаются в электронных системах. По мере увеличения магнитного поля все больше и больше электронов могут поместиться на данном уровне Ландау. Заполнение самого высокого уровня Ландау варьируется от полностью заполненного до полностью пустого, что приводит к колебаниям в различных электронных свойствах (см. эффект Де Гааза–Ван Альфена и эффект Шубникова–де Гааза ).

Если включить расщепление Зеемана , каждый уровень Ландау расщепляется на пару, одну для электронов со спином вверх, а другую для электронов со спином вниз. Тогда заполнение каждого уровня Ландау со спином равно просто отношению потоков . Расщепление Зеемана оказывает значительное влияние на уровни Ландау, поскольку их энергетические масштабы одинаковы, . Однако энергия Ферми и энергия основного состояния остаются примерно одинаковыми в системе со многими заполненными уровнями, поскольку пары расщепленных энергетических уровней компенсируют друг друга при суммировании. D = Φ / Φ 0 {\displaystyle D=\Phi /\Phi _{0}} 2 μ B B = ω c {\displaystyle 2\mu _{B}B=\hbar \omega _{c}}

Более того, приведенный выше вывод в калибровке Ландау предполагал, что электрон ограничен в -направлении, что является соответствующей экспериментальной ситуацией — обнаруженной, например, в двумерных электронных газах. Тем не менее, это предположение не является существенным для результатов. Если электроны могут свободно двигаться вдоль -направления , волновая функция приобретает дополнительный мультипликативный член ; энергия, соответствующая этому свободному движению, , добавляется к обсуждаемому. Этот член затем заполняет разделение по энергии различных уровней Ландау, размывая эффект квантования. Тем не менее, движение в - -плоскости, перпендикулярной магнитному полю, по-прежнему квантуется. z {\displaystyle z} z {\displaystyle z} exp ( i k z z ) {\displaystyle \exp(ik_{z}z)} ( k z ) 2 / ( 2 m ) {\displaystyle (\hbar k_{z})^{2}/(2m)} E {\displaystyle E} x {\displaystyle x} y {\displaystyle y}

В симметричной калибровке

Каждый уровень Ландау имеет вырожденные орбитали, помеченные квантовыми числами в симметричной калибровке. Вырождение на единицу площади одинаково на каждом уровне Ландау. m z {\displaystyle m_{z}}

Z -компонента момента импульса равна L ^ z = i θ = ( b ^ b ^ a ^ a ^ ) {\displaystyle {\hat {L}}_{z}=-i\hbar {\frac {\partial }{\partial \theta }}=-\hbar ({\hat {b}}^{\dagger }{\hat {b}}-{\hat {a}}^{\dagger }{\hat {a}})}

Используя это свойство, мы выбрали собственные функции, которые диагонализируют и , Собственное значение обозначается как , где ясно, что на уровне Ландау. Однако оно может быть сколь угодно большим, что необходимо для получения бесконечного вырождения (или конечного вырождения на единицу площади), демонстрируемого системой. [ H ^ , L ^ z ] = 0 {\displaystyle [{\hat {H}},{\hat {L}}_{z}]=0} H ^ {\displaystyle {\hat {H}}} L ^ z {\displaystyle {\hat {L}}_{z}} L ^ z {\displaystyle {\hat {L}}_{z}} m z {\displaystyle -m_{z}\hbar } m z n {\displaystyle m_{z}\geq -n} n {\displaystyle n}

Релятивистский случай

Уровни Ландау в графене . Носители заряда в графене ведут себя как релятивистские безмассовые дираковские частицы.

Электрон, следующий уравнению Дирака в постоянном магнитном поле, может быть решен аналитически. [6] [7] Энергии определяются как E r e l = ± ( m c 2 ) 2 + ( c k z ) 2 + 2 ν ω c m c 2 {\displaystyle E_{\rm {rel}}=\pm {\sqrt {(mc^{2})^{2}+(c\hbar k_{z})^{2}+2\nu \hbar \omega _{\rm {c}}mc^{2}}}}

где c — скорость света, знак зависит от компонента частица-античастица, а ν — неотрицательное целое число. Из-за спина все уровни вырождены, за исключением основного состояния при ν = 0 .

Безмассовый двумерный случай можно смоделировать в однослойных материалах , таких как графен вблизи конусов Дирака , где собственные энергии задаются выражением [8], где скорость света необходимо заменить на скорость Ферми v F материала, а знак минус соответствует электронным дыркам . E g r a p h e n e = ± 2 ν e B v F 2 {\displaystyle E_{\rm {graphene}}=\pm {\sqrt {2\nu \hbar eBv_{\rm {F}}^{2}}}}

Магнитная восприимчивость ферми-газа

Ферми -газ (ансамбль невзаимодействующих фермионов ) является частью основы для понимания термодинамических свойств металлов. В 1930 году Ландау вывел оценку магнитной восприимчивости ферми-газа, известной как восприимчивость Ландау , которая постоянна для малых магнитных полей. Ландау также заметил, что восприимчивость осциллирует с высокой частотой для больших магнитных полей, [9] это физическое явление известно как эффект Де Гааза–Ван Альфена .

Двумерная решетка

Известно, что энергетический спектр сильной связи заряженных частиц в двумерной бесконечной решетке является самоподобным и фрактальны , как показано в бабочке Хофштадтера . Для целого отношения кванта магнитного потока и магнитного потока через ячейку решетки, можно восстановить уровни Ландау для больших целых чисел. [10]

Целочисленный квантовый эффект Холла

Энергетический спектр полупроводника в сильном магнитном поле образует уровни Ландау, которые можно обозначить целыми индексами. Кроме того, сопротивление Холла также демонстрирует дискретные уровни, обозначенные целым числом ν . Тот факт, что эти две величины связаны, можно показать разными способами, но проще всего увидеть из модели Друде : проводимость Холла зависит от электронной плотности n как

ρ x y = B n e . {\displaystyle \rho _{xy}={\frac {B}{ne}}.}

Поскольку плато сопротивления определяется выражением

ρ x y = 2 π e 2 1 ν , {\displaystyle \rho _{xy}={\frac {2\pi \hbar }{e^{2}}}{\frac {1}{\nu }},}

требуемая плотность

n = B Φ 0 ν , {\displaystyle n={\frac {B}{\Phi _{0}}}\nu ,}

что является как раз той плотностью, которая требуется для заполнения уровня Ландау. Разрыв между различными уровнями Ландау вместе с большим вырождением каждого уровня делает сопротивление квантованным.

Смотрите также

Ссылки

  1. ^ Ландау, Л. (1930). «Diamagnetismus der Metalle» [Диамагнетизм металлов]. Zeitschrift für Physik (на немецком языке). 64 ( 9–10 ). Springer Science and Business Media LLC: 629–637 . Бибкод : 1930ZPhy...64..629L. дои : 10.1007/bf01397213. ISSN  1434-6001. S2CID  123206025.
  2. ^ "Заряд в магнитном поле" (PDF) . courses.physics.illinois.edu . Получено 11 марта 2023 г. .
  3. ^ Столь же правильным решением в калибровке Ландау было бы: . A = ( B y 0 0 ) T {\displaystyle \mathbf {A} ={\begin{pmatrix}-By&0&0\end{pmatrix}}^{T}}
  4. ^ Ландау, Л. Д .; Лифшиц, Э. М. (1977). Квантовая механика: нерелятивистская теория (3-е изд.). Амстердам: Butterworth Heinemann. С.  424–426 . ISBN 978-0-7506-3539-4. OCLC  846962062.
  5. ^ Михайлов, СА (2001). "Новый подход к основному состоянию квантовых холловских систем. Основные принципы". Physica B: Condensed Matter . 299 ( 1– 2): 6– 31. arXiv : cond-mat/0008227 . Bibcode :2001PhyB..299....6M. doi :10.1016/S0921-4526(00)00769-9. S2CID  118500817.
  6. ^ Раби, II (1928). «Свободный электрон в однородном Магнитфельде в теории Диракшена». Zeitschrift für Physik (на немецком языке). 49 ( 7–8 ): 507–511 . Бибкод : 1928ZPhy...49..507R. дои : 10.1007/BF01333634. ISSN  1434-6001. S2CID  121121095.
  7. ^ Берестецкий, В.Б.; Питаевский, Л.П.; Лифшиц, Э.М. (2012-12-02). Квантовая электродинамика: Том 4. Elsevier. ISBN 978-0-08-050346-2.
  8. ^ Инь, Лонг-Цзин; Бай, Кэ-Ке; Ван, Вэнь-Сяо; Ли, Си-Ю; Чжан, Юй; Хэ, Линь (2017). «Квантование Ландау фермионов Дирака в графене и его многослойных структурах». Frontiers of Physics . 12 (4): 127208. arXiv : 1703.04241 . Bibcode : 2017FrPhy..12l7208Y. doi : 10.1007/s11467-017-0655-0 . ISSN  2095-0462.
  9. ^ Ландау, Л. Д.; Лифшиц, Э. М. (22 октября 2013 г.). Статистическая физика: Том 5. Elsevier. С. 177. ISBN 978-0-08-057046-4.
  10. ^ Analytis, James G.; Blundell, Stephen J.; Ardavan, Arzhang (май 2004 г.). «Уровни Ландау, молекулярные орбитали и бабочка Хофштадтера в конечных системах». American Journal of Physics . 72 (5): 613– 618. Bibcode :2004AmJPh..72..613A. doi :10.1119/1.1615568. ISSN  0002-9505.
  • Лев Ландау (1930). «Diamagnetismus der Metalle» (PDF) (на немецком языке).{{cite web}}: CS1 maint: numeric names: authors list (link)

Дальнейшее чтение

  • Ландау, Л. Д.; и Лифшиц, Э. М.; (1977). Квантовая механика: нерелятивистская теория. Курс теоретической физики . Том 3 (3-е изд. Лондон: Pergamon Press). ISBN 0750635398 . 
Retrieved from "https://en.wikipedia.org/w/index.php?title=Landau_levels&oldid=1263711269"