Рівняння Дірака (графен)

Якщо врахувати тільки вклад найближчих сусідів у формування енергетичних зон, то гамільтоніан в наближенні сильного зв'язку для гексагональної ґратки приймає вигляд:

${\displaystyle H=-t\sum _{i\in \Lambda _{A}}\sum _{j=1}^{3}a^{\dagger }({\textbf {r}}_{i})b({\textbf {r}}_{i}+{\textbf {u}}_{j})-t\sum _{i\in \Lambda _{B}}\sum _{j=1}^{3}b^{\dagger }({\textbf {r}}_{i})a({\textbf {r}}_{i}+{\textbf {v}}_{j}),\qquad (1.1)}$

де ${\displaystyle t}$ — інтеграл перекриття між хвильовими функціями найближчих сусідів, який визначає також ймовірність переходу («стрибка») між сусідніми атомами (атомами з різних підрешіток), оператори ${\displaystyle a^{\dagger }({\textbf {r}}_{i})}$ та ${\displaystyle b^{\dagger }({\textbf {r}}_{i})}$ оператори народження, які діють на трикутних підрешітках кристалу ${\displaystyle \Lambda _{A}}$ та ${\displaystyle \Lambda _{B}}$ відповідно, ${\displaystyle a({\textbf {r}}_{i})}$ та ${\displaystyle b({\textbf {r}}_{i})}$ — оператори знищення. Вони задовольняють звичайній антикомутаційним співвідношенням для ферміонів:

${\displaystyle [a({\textbf {r}}_{i}),a^{\dagger }({\textbf {r}}_{i^{'}})]_{+}=[b({\textbf {r}}_{i}),b^{\dagger }({\textbf {r}}_{i^{'}})]_{+}=\delta _{ii^{'}}.\qquad (1.2)}$

Шість векторів ${\displaystyle {\textbf {u}}_{i}}$ та ${\displaystyle {\textbf {v}}_{i}}$ вказують на найближчі вузли від вибранного центрального атому і задаються відношеннями

${\displaystyle {\textbf {u}}_{1}=(-d,0),\,{\textbf {u}}_{2}=\left({\frac {1}{2}}d,{\frac {\sqrt {3}}{2}}d\right),\,{\textbf {u}}_{3}=\left({\frac {1}{2}}d,-{\frac {\sqrt {3}}{2}}d\right),\qquad (1.3)}$

${\displaystyle {\textbf {v}}_{1}=(d,0),\,{\textbf {v}}_{2}=\left(-{\frac {1}{2}}d,-{\frac {\sqrt {3}}{2}}d\right),\,{\textbf {v}}_{3}=\left(-{\frac {1}{2}}d,{\frac {\sqrt {3}}{2}}d\right).\qquad (1.4)}$

Перетворення Фур'є операторів народження та знищення

${\displaystyle a({\textbf {r}}_{i})=\int \limits _{BZ}{\frac {d^{2}k}{(2\pi )^{2}}}e^{i{\textbf {k}}{\textbf {r}}_{i}}{\tilde {a}}({\textbf {k}}),b({\textbf {r}}_{i})=\int \limits _{BZ}{\frac {d^{2}k}{(2\pi )^{2}}}e^{i{\textbf {k}}{\textbf {r}}_{i}}{\tilde {b}}({\textbf {k}}),\qquad (1.5)}$

де інтегрирування по хвильовим векторам проводиться з першої зони Бриллюена, дозволяє записати гамильтониан у вигляді

${\displaystyle H=\int \limits _{BZ}{\frac {d^{2}k}{(2\pi )^{2}}}{\tilde {\psi }}^{\dagger }({\textbf {k}}){\tilde {H}}{\tilde {\psi }}({\textbf {k}}),\qquad (1.6)}$

де прийняті такі позначення:

${\displaystyle {\tilde {\psi }}({\textbf {k}})=\left({\tilde {a}}({\textbf {k}}),{\tilde {b}}({\textbf {k}})\right)^{T},\,{\tilde {\psi }}^{\dagger }({\textbf {k}})=\left({\tilde {a}}^{\dagger }({\textbf {k}}),{\tilde {b}}^{\dagger }({\textbf {k}})\right),\qquad (1.7)}$

та

${\displaystyle {\tilde {H}}=\left({\begin{array}{cc}0&-t\sum _{j=1}^{3}e^{i{\textbf {k}}{\textbf {u}}_{j}}\\-t\sum _{j=1}^{3}e^{i{\textbf {k}}{\textbf {v}}_{j}}&0\\\end{array}}\right).\qquad (1.8)}$

Вираз (1.6) можна отримати якщо підставити (1.5) в (1.1). Розглянемо суму

${\displaystyle \sum _{i\in \Lambda _{A}}\sum _{j=1}^{3}a^{\dagger }({\textbf {r}}_{i})b({\textbf {r}}_{i}+{\textbf {u}}_{j}),\qquad (1.9)}$

яку, використавши співвідношення (1.5) можна записати у вигляді

${\displaystyle \sum _{i\in \Lambda _{A}}\sum _{j=1}^{3}\int \limits _{BZ}{\frac {d^{2}k}{(2\pi )^{2}}}e^{-i{\textbf {k}}{\textbf {r}}_{i}}{\tilde {a\dagger }}({\textbf {k}})\int \limits _{BZ}{\frac {d^{2}k^{'}}{(2\pi )^{2}}}e^{i{\textbf {k}}^{'}({\textbf {r}}_{i}+{\textbf {u}}_{j})}{\tilde {b}}({\textbf {k}}^{'}),\qquad (1.10)}$

або

${\displaystyle \int \limits _{BZ}{\frac {d^{2}k}{(2\pi )^{2}}}{\tilde {a\dagger }}({\textbf {k}})\int \limits _{BZ}{\frac {d^{2}k^{'}}{(2\pi )^{2}}}\sum _{i\in \Lambda _{A}}e^{-i{\textbf {k}}{\textbf {r}}_{i}+i{\textbf {k}}^{'}{\textbf {r}}_{i}}\sum _{j=1}^{3}e^{i{\textbf {k}}^{'}{\textbf {u}}_{j}}{\tilde {b}}({\textbf {k}}^{'}).\qquad (1.11)}$

Використавши співвідношення

${\displaystyle \sum _{i\in \Lambda _{A}}e^{-i{\textbf {k}}{\textbf {r}}_{i}+i{\textbf {k}}^{'}{\textbf {r}}_{i}}=(2\pi )^{2}\delta \left({\textbf {k}}^{'}-{\textbf {k}}\right),\qquad (1.12)}$

знаходимо після інтегрування по ${\displaystyle {\textbf {k}}^{'}}$ вираз

${\displaystyle \int \limits _{BZ}{\frac {d^{2}k}{(2\pi )^{2}}}{\tilde {a\dagger }}({\textbf {k}})\sum _{j=1}^{3}e^{i{\textbf {k}}{\textbf {u}}_{j}}{\tilde {b}}({\textbf {k}}).\qquad (1.13)}$

Аналогічне перетворення другої суми в гамільтоніані (1.1) приводить до бажаного результату (1.6).

Власне значення гамільтоніану (1.8) приймає значення

${\displaystyle E=\pm t{\sqrt {\sum _{j=1}^{3}e^{i{\textbf {k}}{\textbf {u}}_{j}}\sum _{j^{'}=1}^{3}e^{i{\textbf {k}}{\textbf {v}}_{j^{'}}}}}=\pm t{\sqrt {\left(e^{-ik_{x}d}+2e^{ik_{x}d/2}\cos {{\frac {\sqrt {3}}{2}}dk_{y}}\right)\left(e^{ik_{x}d}+2e^{-ik_{x}d/2}\cos {{\frac {\sqrt {3}}{2}}dk_{y}}\right)}}=}$

${\displaystyle \pm t{\sqrt {\left(1+2e^{i3k_{x}d/2}\cos {{\frac {\sqrt {3}}{2}}dk_{y}}\right)\left(1+2e^{-i3k_{x}d/2}\cos {{\frac {\sqrt {3}}{2}}dk_{y}}\right)}}=\pm t{\sqrt {1+4\cos \left({\frac {\sqrt {3}}{2}}k_{y}d\right)\left[\cos \left({\frac {3}{2}}k_{x}d\right)+\cos \left({\frac {\sqrt {3}}{2}}k_{y}d\right)\right]}},\qquad (1.14)}$

яке визначає зонну структуру графена.[2]

Зони (1.14) з позитивною енергією (квазічастки - електрони) та з негативною енергією (квазічастки - дірки) перетинаються в шести точках, які називаються діраковськими точками, оскільки поблизу них енергетичний спектр приймає "лінійну" залежність від хвильового вектора. Координати цих точок рівні

${\displaystyle \left(0,{\frac {4\pi }{3{\sqrt {3}}d}}\right),\,\left(0,-{\frac {4\pi }{3{\sqrt {3}}d}}\right),\,\left({\frac {2\pi }{3d}},{\frac {2\pi }{3{\sqrt {3}}d}}\right),\,\left({\frac {2\pi }{3d}},-{\frac {2\pi }{3{\sqrt {3}}d}}\right),\,\left(-{\frac {2\pi }{3d}},{\frac {2\pi }{3{\sqrt {3}}d}}\right),\,\left(-{\frac {2\pi }{3d}},{\frac {-2\pi }{3{\sqrt {3}}d}}\right).\qquad (2.1)}$

Дві незалежні долини можна вибрати так, що вершини валентних зон будуть знаходитися в діраковських точах з координатами

${\displaystyle {\textbf {K}}^{\pm }=\left(0,\pm {\frac {4\pi }{3{\sqrt {3}}d}}\right).\qquad (2.2)}$

Розглянемо недіагональний елемент ^{${\displaystyle {\tilde {H}}_{12}}$} оператора Гамільтона (1.8). Разкладемо його поблизу точок Дірака (2.2) по малому параметру d

${\displaystyle \lim _{d\rightarrow 0}d^{-1}{\tilde {H}}_{12}|_{{\textbf {k}}={\textbf {K}}^{\pm }+{\boldsymbol {\kappa }}}=-t\lim _{d\rightarrow 0}d^{-1}\left(e^{-i\kappa _{x}d}+2e^{i\kappa _{x}d/2}\cos {\frac {{\sqrt {3}}d}{2}}\left(\pm {\frac {4\pi }{3{\sqrt {3}}d}}+\kappa _{y}\right)\right)={\frac {3t}{2}}(i\kappa _{x}\pm \kappa _{y}).\qquad (2.3)}$

Для ^{${\displaystyle {\tilde {H}}_{21}}$} розклад обчислюється аналогічним чином, тому в результаті можна записати гамільтоніан для квазічасток поблизу точок Дірака у вигляді

${\displaystyle \left({\begin{array}{cc}H^{+}&0\\0&H^{-}\\\end{array}}\right)=\hbar v_{F}(\alpha ^{1}\kappa _{x}+\alpha ^{2}\kappa _{y}),\qquad (2.4)}$

де швидкість Фермі ${\displaystyle v_{F}=3td\hbar ^{-1}/2}$ та

${\displaystyle \alpha ^{1}=-\left({\begin{array}{cc}\sigma ^{2}&0\\0&\sigma ^{2}\\\end{array}}\right),\,\alpha ^{2}=\left({\begin{array}{cc}\sigma ^{1}&0\\0&-\sigma ^{1}\\\end{array}}\right).\qquad (2.5)}$

Тут ^{${\displaystyle \sigma ^{1}}$} та ^{${\displaystyle \sigma ^{2}}$} — матриці Паулі.

Якщо тепер перейти до координатного представлення, зробивши фур'є перетворення гамільтоніану (2.4), то приходимо до гамільтоніану в рівнянні Дірака для квазічасток в графені

${\displaystyle H=-i\hbar v_{F}(\alpha ^{1}\partial _{x}+\alpha ^{2}\partial _{y}).\qquad (2.6)}$

Розв'язком рівняння Дірака для графену ${\displaystyle H\psi =E\psi }$ буде чотирьохкомпонентний стовбчик вигляду

${\displaystyle \psi =(\psi _{A}^{+},\psi _{B}^{+},\psi _{A}^{-},\psi _{B}^{-})^{T},\qquad (2.7)}$

де індекси ^{${\displaystyle A}$} та ^{${\displaystyle B}$} відповідають двом підрешіткам кристалу, а знаки «+» та «-» позначають нееквівалентні точки Дірака в k-просторі.[2]