Исследуйте механизм квантового туннелирования в сканирующем туннельном микроскопе (СТМ): от чего зависят поток туннелирующих электронов, пространственное разрешение и влияние температуры, и как моделирование потенциала контакта улучшит объяснение результатов
Кратко и по делу. Механизм (основная идея) - Туннелирование — квантово-механическое прохождение электронов через потенциальный барьер вакуума между остриём зонда и образцом; текущее значение определяется скоростью передачи состояний и коэффициентом передачи барьера. - Общая формула (в упрощённом виде, схема Бордина/ЗГФ): I∝∫−∞+∞[f(E−eV)−f(E)] ns(E) nt(E−eV) T(E) dE,
I \propto \int_{-\infty}^{+\infty}\bigl[f(E-eV)-f(E)\bigr]\; n_s(E)\,n_t(E-eV)\;T(E)\;dE, I∝∫−∞+∞[f(E−eV)−f(E)]ns(E)nt(E−eV)T(E)dE,
где f(E)=(e(E−μ)/kBT+1)−1f(E)=(e^{(E-\mu)/k_BT}+1)^{-1}f(E)=(e(E−μ)/kBT+1)−1 — функция Ферми, ns,tn_{s,t}ns,t — плотности состояний образца и зонда, T(E)T(E)T(E) — коэффициент пропускания барьера. Коэффициент пропускания (WKB-приближение) - Для барьера толщиной ddd и локальной высотой U(z)U(z)U(z) в WKB: T(E)≈exp(−2∫z1z2κ(z) dz),κ(z)=2mℏ2(U(z)−E).
T(E)\approx \exp\Bigl(-2\int_{z_1}^{z_2}\kappa(z)\,dz\Bigr),\qquad \kappa(z)=\sqrt{\frac{2m}{\hbar^2}\bigl(U(z)-E\bigr)}. T(E)≈exp(−2∫z1z2κ(z)dz),κ(z)=ℏ22m(U(z)−E).
- Для приблизительно прямоугольного вакуумного барьера высотой ϕ\phiϕ: T(E)≈exp(−2κd),κ=2m(ϕ−E)ℏ2.
T(E)\approx\exp(-2\kappa d),\qquad \kappa=\sqrt{\frac{2m(\phi-E)}{\hbar^2}}. T(E)≈exp(−2κd),κ=ℏ22m(ϕ−E).
Из этого вытекает экспоненциальная зависимость тока от расстояния ddd: при прочих равных I∝exp(−2κd)I\propto \exp(-2\kappa d)I∝exp(−2κd). Зависимость тока от параметров - Расстояние ddd: самое сильное влияние (экспоненциально). Небольшое изменение ddd на ∼1 A˚\sim 1\,\text{Å}∼1A˚ даёт порядок величины изменения тока при обычных ϕ\phiϕ. - Рабочая функция / барьер ϕ\phiϕ: влияет через κ\kappaκ. При ϕ≈4 eV\phi\approx 4\ \text{eV}ϕ≈4eV получается примерно κ∼1 A˚−1\kappa\sim 1\ \text{Å}^{-1}κ∼1A˚−1. - Напряжение VVV: смещает энергетические окна интеграла и позволяет измерять спектроскопию (dI/dV ~ LDOS). - Плотности состояний ns,ntn_s,n_tns,nt и матричный элемент связи MnmM_{nm}Mnm (зависит от формы орбиталей острия и поверхности) — определяют энергетическую зависимость и симметрию тока. - Неупругие процессы (фононы, спины) добавляют вклад в ток и меняют спектральную форму. Пространственное разрешение - Латеральное разрешение определяется экспоненциальным спадом волновой функции в вакууме и симметрией/локализацией орбиталей острия: - Чем больше κ\kappaκ, тем сильнее локализуется вклад от точечных выступов и тем лучше разрешение. - Острие с атомарно острым апексом и орбиталями высокой локализации (например, sss- или pzp_zpz-характер) даёт лучший контраст, апекс радиус важен. - В приближении Тёрсоффа–Хэмана для простого сферического острия: I(V)∝∫EFEF+eVρs(r0,E) dE,
I(V)\propto\int_{E_F}^{E_F+eV}\rho_s(\mathbf{r}_0,E)\,dE, I(V)∝∫EFEF+eVρs(r0,E)dE,
где ρs(r0,E)\rho_s(\mathbf{r}_0,E)ρs(r0,E) — LDOS в точке r0 \mathbf{r}_0r0 на поверхности прямо под апексом; это объясняет атомарное разрешение при локальной чувствительности LDOS. - Практически: латеральное разрешение порядка единиц ангстрем при стандартных ϕ\phiϕ и расстояниях. Влияние температуры - Энергетическое «смешивание» уровней определяется kBTk_B TkBT. Энергетическая разрешающая способность примерно ограничена масштабом kBTk_B TkBT: ΔEthermal∼kBT.
\Delta E_{\text{thermal}}\sim k_B T. ΔEthermal∼kBT.
Примеры: kBT≈25 meVk_B T\approx 25\ \text{meV}kBT≈25meV при 300 K300\ \text{K}300K, kBT≈0.34 meVk_B T\approx 0.34\ \text{meV}kBT≈0.34meV при 4 K4\ \text{K}4K. - При повышенной температуре: размываются спектральные особенности (dI/dV), усиливаются тепловые флуктуации расстояния (дрожание зонда), увеличивается вклад неупругих каналов. - Температура также влияет на фононную активность и внетуннельные процессы (IETS). Как моделирование потенциала контакта улучшает объяснение результатов - Простые модели (ровный прямоугольный барьер, Tersoff–Hamann) объясняют базовую экспоненциальную зависимость и LDOS-зондирование, но не учитывают: - форму реального барьера — электростатическое смещение под напряжением, билинейную кривизну, изображенийный потенциал (Uim(z)∝−1/zU_{\rm im}(z)\propto -1/zUim(z)∝−1/z) ; - атомарную структуру апекса и многорежимные орбитали зонда/поверхности; - пространственную вариацию рабочего потенциала (локальная зарядовая перераспределение, бэндбендинг); - неупругие процессы и многоканальную передачу. - Улучшения при учёте реального потенциала даёт: - Правильную величину и зависимость тока от смещения и напряжения (асимметрия I–V, ступени из-за локальных состояний). - Корректную форму барьера (включая image-potential) — меняется κ(z)\kappa(z)κ(z) и, следовательно, дистанционная и энергетическая зависимость тока. - Описание локальных эффектов (например, резонансы адатомов, молекулярные орбитали) и электростатического поля между зондом и образцом. - Методы моделирования: - Tersoff–Hamann: простая интерпретация STM как зонд LDOS (быстро, полезно для STM-изображений). - Bardeen: выражение через матричные элементы наплавленного волнового поля; полезно для учёта орбитальной симметрии. Mmn=ℏ22m∫S(ψm∗∇ψn−ψn∇ψm∗)⋅dS.
M_{mn}=\frac{\hbar^2}{2m}\int_S\bigl(\psi_m^*\nabla\psi_n-\psi_n\nabla\psi_m^*\bigr)\cdot d\mathbf{S}. Mmn=2mℏ2∫S(ψm∗∇ψn−ψn∇ψm∗)⋅dS.
- WKB + корректировки (image potential) — для оценки T(E)T(E)T(E) с координатной зависимостью барьера. - DFT + NEGF (или DFT+transport): полное квантово-химическое моделирование структуры контакта, локальных потенциалов и передачи при конечном VVV и TTT — рекомендуется для количественного сравнения с экспериментом. - Tight-binding / многорежимные модели и включение фононов для описания неупругих вкладов (IETS). - Практический эффект: более точное моделирование позволяет разделить вклад геометрии (форма острия), локальной электронной структуры (LDOS, резонансы) и элекростатических эффектов (барьер/смещение), что улучшает интерпретацию карт и спектров STM, особенно при больших напряжениях и сложных системах (молекулы, полупроводники). Короткие выводы - Ток STM экспоненциально зависит от расстояния и барьера; спектрально — от DOS и матричных элементов. - Разрешение определяется экспоненциальным спадом волновой функции и формой орбиталей острия. - Температура ограничивает энергетическое разрешение через kBTk_B TkBT и усиливает флуктуации и неупругие процессы. - Подробное моделирование потенциального профиля (image-potential, поле при смещении, атомная структура апекса) и использование DFT+NEGF/включение неупругих каналов существенно повышает согласование теории и эксперимента и проясняет физику наблюдаемых эффектов.
Механизм (основная идея)
- Туннелирование — квантово-механическое прохождение электронов через потенциальный барьер вакуума между остриём зонда и образцом; текущее значение определяется скоростью передачи состояний и коэффициентом передачи барьера.
- Общая формула (в упрощённом виде, схема Бордина/ЗГФ):
I∝∫−∞+∞[f(E−eV)−f(E)] ns(E) nt(E−eV) T(E) dE, I \propto \int_{-\infty}^{+\infty}\bigl[f(E-eV)-f(E)\bigr]\; n_s(E)\,n_t(E-eV)\;T(E)\;dE,
I∝∫−∞+∞ [f(E−eV)−f(E)]ns (E)nt (E−eV)T(E)dE, где f(E)=(e(E−μ)/kBT+1)−1f(E)=(e^{(E-\mu)/k_BT}+1)^{-1}f(E)=(e(E−μ)/kB T+1)−1 — функция Ферми, ns,tn_{s,t}ns,t — плотности состояний образца и зонда, T(E)T(E)T(E) — коэффициент пропускания барьера.
Коэффициент пропускания (WKB-приближение)
- Для барьера толщиной ddd и локальной высотой U(z)U(z)U(z) в WKB:
T(E)≈exp(−2∫z1z2κ(z) dz),κ(z)=2mℏ2(U(z)−E). T(E)\approx \exp\Bigl(-2\int_{z_1}^{z_2}\kappa(z)\,dz\Bigr),\qquad
\kappa(z)=\sqrt{\frac{2m}{\hbar^2}\bigl(U(z)-E\bigr)}.
T(E)≈exp(−2∫z1 z2 κ(z)dz),κ(z)=ℏ22m (U(z)−E) . - Для приблизительно прямоугольного вакуумного барьера высотой ϕ\phiϕ:
T(E)≈exp(−2κd),κ=2m(ϕ−E)ℏ2. T(E)\approx\exp(-2\kappa d),\qquad
\kappa=\sqrt{\frac{2m(\phi-E)}{\hbar^2}}.
T(E)≈exp(−2κd),κ=ℏ22m(ϕ−E) . Из этого вытекает экспоненциальная зависимость тока от расстояния ddd: при прочих равных I∝exp(−2κd)I\propto \exp(-2\kappa d)I∝exp(−2κd).
Зависимость тока от параметров
- Расстояние ddd: самое сильное влияние (экспоненциально). Небольшое изменение ddd на ∼1 A˚\sim 1\,\text{Å}∼1A˚ даёт порядок величины изменения тока при обычных ϕ\phiϕ.
- Рабочая функция / барьер ϕ\phiϕ: влияет через κ\kappaκ. При ϕ≈4 eV\phi\approx 4\ \text{eV}ϕ≈4 eV получается примерно κ∼1 A˚−1\kappa\sim 1\ \text{Å}^{-1}κ∼1 A˚−1.
- Напряжение VVV: смещает энергетические окна интеграла и позволяет измерять спектроскопию (dI/dV ~ LDOS).
- Плотности состояний ns,ntn_s,n_tns ,nt и матричный элемент связи MnmM_{nm}Mnm (зависит от формы орбиталей острия и поверхности) — определяют энергетическую зависимость и симметрию тока.
- Неупругие процессы (фононы, спины) добавляют вклад в ток и меняют спектральную форму.
Пространственное разрешение
- Латеральное разрешение определяется экспоненциальным спадом волновой функции в вакууме и симметрией/локализацией орбиталей острия:
- Чем больше κ\kappaκ, тем сильнее локализуется вклад от точечных выступов и тем лучше разрешение.
- Острие с атомарно острым апексом и орбиталями высокой локализации (например, sss- или pzp_zpz -характер) даёт лучший контраст, апекс радиус важен.
- В приближении Тёрсоффа–Хэмана для простого сферического острия:
I(V)∝∫EFEF+eVρs(r0,E) dE, I(V)\propto\int_{E_F}^{E_F+eV}\rho_s(\mathbf{r}_0,E)\,dE,
I(V)∝∫EF EF +eV ρs (r0 ,E)dE, где ρs(r0,E)\rho_s(\mathbf{r}_0,E)ρs (r0 ,E) — LDOS в точке r0 \mathbf{r}_0r0 на поверхности прямо под апексом; это объясняет атомарное разрешение при локальной чувствительности LDOS.
- Практически: латеральное разрешение порядка единиц ангстрем при стандартных ϕ\phiϕ и расстояниях.
Влияние температуры
- Энергетическое «смешивание» уровней определяется kBTk_B TkB T. Энергетическая разрешающая способность примерно ограничена масштабом kBTk_B TkB T:
ΔEthermal∼kBT. \Delta E_{\text{thermal}}\sim k_B T.
ΔEthermal ∼kB T. Примеры: kBT≈25 meVk_B T\approx 25\ \text{meV}kB T≈25 meV при 300 K300\ \text{K}300 K, kBT≈0.34 meVk_B T\approx 0.34\ \text{meV}kB T≈0.34 meV при 4 K4\ \text{K}4 K.
- При повышенной температуре: размываются спектральные особенности (dI/dV), усиливаются тепловые флуктуации расстояния (дрожание зонда), увеличивается вклад неупругих каналов.
- Температура также влияет на фононную активность и внетуннельные процессы (IETS).
Как моделирование потенциала контакта улучшает объяснение результатов
- Простые модели (ровный прямоугольный барьер, Tersoff–Hamann) объясняют базовую экспоненциальную зависимость и LDOS-зондирование, но не учитывают:
- форму реального барьера — электростатическое смещение под напряжением, билинейную кривизну, изображенийный потенциал (Uim(z)∝−1/zU_{\rm im}(z)\propto -1/zUim (z)∝−1/z) ;
- атомарную структуру апекса и многорежимные орбитали зонда/поверхности;
- пространственную вариацию рабочего потенциала (локальная зарядовая перераспределение, бэндбендинг);
- неупругие процессы и многоканальную передачу.
- Улучшения при учёте реального потенциала даёт:
- Правильную величину и зависимость тока от смещения и напряжения (асимметрия I–V, ступени из-за локальных состояний).
- Корректную форму барьера (включая image-potential) — меняется κ(z)\kappa(z)κ(z) и, следовательно, дистанционная и энергетическая зависимость тока.
- Описание локальных эффектов (например, резонансы адатомов, молекулярные орбитали) и электростатического поля между зондом и образцом.
- Методы моделирования:
- Tersoff–Hamann: простая интерпретация STM как зонд LDOS (быстро, полезно для STM-изображений).
- Bardeen: выражение через матричные элементы наплавленного волнового поля; полезно для учёта орбитальной симметрии.
Mmn=ℏ22m∫S(ψm∗∇ψn−ψn∇ψm∗)⋅dS. M_{mn}=\frac{\hbar^2}{2m}\int_S\bigl(\psi_m^*\nabla\psi_n-\psi_n\nabla\psi_m^*\bigr)\cdot d\mathbf{S}.
Mmn =2mℏ2 ∫S (ψm∗ ∇ψn −ψn ∇ψm∗ )⋅dS. - WKB + корректировки (image potential) — для оценки T(E)T(E)T(E) с координатной зависимостью барьера.
- DFT + NEGF (или DFT+transport): полное квантово-химическое моделирование структуры контакта, локальных потенциалов и передачи при конечном VVV и TTT — рекомендуется для количественного сравнения с экспериментом.
- Tight-binding / многорежимные модели и включение фононов для описания неупругих вкладов (IETS).
- Практический эффект: более точное моделирование позволяет разделить вклад геометрии (форма острия), локальной электронной структуры (LDOS, резонансы) и элекростатических эффектов (барьер/смещение), что улучшает интерпретацию карт и спектров STM, особенно при больших напряжениях и сложных системах (молекулы, полупроводники).
Короткие выводы
- Ток STM экспоненциально зависит от расстояния и барьера; спектрально — от DOS и матричных элементов.
- Разрешение определяется экспоненциальным спадом волновой функции и формой орбиталей острия.
- Температура ограничивает энергетическое разрешение через kBTk_B TkB T и усиливает флуктуации и неупругие процессы.
- Подробное моделирование потенциального профиля (image-potential, поле при смещении, атомная структура апекса) и использование DFT+NEGF/включение неупругих каналов существенно повышает согласование теории и эксперимента и проясняет физику наблюдаемых эффектов.