Квантовая физика и нанотехнологии
.pdf
Глава 10. Движение частицы в поле потенциальной ступеньки
Тогда выражения для энергии частиц можно записать в виде:
E |
P12 |
2 2( 10)2 , |
|
||||||||
|
|
|
|||||||||
|
|
|
2m |
|
|
8m |
|
|
|
|
|
E |
|
P 2 |
|
2 2( |
30 |
) |
2 |
|
|||
|
3 |
|
|
|
, |
||||||
|
2m |
|
8m |
|
|
||||||
|
|
|
|
|
|
|
|
|
|||
E |
P 2 |
|
U 0 |
|
2 2( |
20 |
)2 |
||||
2 |
|
|
|
|
. |
||||||
|
|
8m |
|
||||||||
|
2m |
|
|
|
|
||||||
Для вычисления коэффициента прозрачности барьера (10.17) найдем выражения для неизвестных импульсов Р2, Р3, пользуясь выписанными значениями полной энергии и определением для
D 20 P2 / 10 P1 :
|
P [mE |
m2E 2 2 |
2(1 D)2 2 |
P 2 |
/ 4]1/ 2, |
(10.18) |
|||||||||
|
3 |
|
|
|
|
|
|
|
10 |
1 |
|
|
|
|
|
|
[m(E U |
|
) |
|
|
|
|
|
2 |
|
|
|
|||
P |
0 |
m2 |
(E U |
0 |
)2 |
2 2D2 |
P 2 |
/ 4]1/ 2. |
(10.19) |
||||||
2 |
|
|
|
|
|
|
|
|
|
10 |
1 |
|
|
||
Здесь выражения для импульсов получены в неявном виде через коэффициент пропускания D. Совместно формулы (10.17)— (10.19) позволяют вычислить коэффициент пропускания барьера в явном виде. Для падающих на потенциальную ступеньку частиц считаются известными E, P1 , 10 . При этом
EP12 2 2( 10)2 .
2m 8m
Полученные выражения для импульсов обеспечивают классический предел D = 1, когда % 0. Коэффициент прохождения D 0 при значениях полной энергии Е = U0 , в этом случае Р3 =
= Р1.
Как правило, квантовая составляющая частиц, падающих на барьер, мала, т.е.
0 2 2( 10)2 !! P12 Ek. 8m 2m
Тогда из (10.18) и (10.19) следует:
P3 P1(1 D)
0 / Ek , P2 
2m(Ek U 0), P1 
2mEk , (10.20)
т.е. величины отраженных импульсов составляют малую долю от величин падающих импульсов.
Полученные формулы для коэффициента прохождения D существенно отличаются от решения аналогичной задачи для квантовых частиц, падающих на потенциальную ступеньку, для которых
Глава 10. Движение частицы в поле потенциальной ступеньки
амплитуда плотности вероятности является постоянной во всем пространстве [1].
Покажем, чем отличается традиционное решение этой задачи от приведенного выше решения.
Для решения уравнения (10.1) предположим, что имеет место суперпозиция квантовых состояний падающих и отраженных частиц, т.е. их координаты не различимы, x1 x3 x
Пусть 1 (x) и 3 (x) являются решениями уравнения (10.1). Можно показать, что суперпозиция этих решений в виде:
13 x 1 3 − 2
1 3
тоже является решением этого уравнения, если E P12 /2m P22 /2m. Однако, чтобы координаты падающих и отраженных частиц были неразличимы вдоль направления движения, они должны быть делокализованы во всем полупространстве. Сделаем соответствующий переход в решениях для 1 , 3 , 2 ( xi % &), тогда получим граничные условия в виде:
10 30 2
10 30 20,
P1( 10 30) P2 20,
что и обеспечивает получение формул задачи 1 для R и D из [1], которые не зависят от постоянной Планка и, следовательно, не обеспечивают классического предела, когда % 0.
D 4P1P2 / (P1 P2)2.
Наличие классического предела для коэффициента прохождения D, конечно, не должно зависеть от крутизны границы потенциальной ступеньки. Дело в традиционном представлении о свободных квантовых частицах, когда свободной частице ставится в соответствие волна де Бройля.
Литература
1.Ландау Л.Д., Лифшиц Е.М. Квантовая механика. Нерелятивистская теория. — М.: Наука, 1974.
ГЛАВА 11
ТУННЕЛИРОВАНИЕ
Задача туннелирования решена в квазигидродинамическом представлении в гл. 7, когда рассматривался альфа-распад ядер. Туннелирование рассматривалось из потенциальной ямы с финитным движением частиц, когда поступательная составляющая частиц отсутствует. Рассмотрим туннелирование частиц в более общем виде, когда они налетают на прямоугольную стенку с некоторым импульсом.
Рассмотрим прохождение частиц через прямоугольный барьер конечной высоты U 0 (E ! U 0 ) и конечной ширины а (рис. 11.1).
Рис. 11.1.
Поскольку в области барьера макроскопический импульс частиц не существует, то необходимо решать систему уравнений (3.2), (3.6).
m div J 0, (3.2)
t
(J / ) |
|
J 2 |
|
2( )2 |
|
2 |
|
|
|
|
|
U |
|
|
|
. |
(3.6) |
t |
2m 2 |
8m |
|
|||||
|
|
|
4m |
|
||||
Для решения этой системы уравнений сделаем упрощающие
предложения. Положим, что при туннелировании в стационарных |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
J |
|
|
|
|
|
|||||
потенциальных полях |
|
|
|
|
|
|
0, т.е. остается постоянной энергия |
|||||||
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
t |
|
|
|
|
|
|
||||
туннелирующих частиц: |
|
|
|
|
|
|
|
|
|
|||||
|
J 2 |
|
|
|
|
|
|
|
2 |
|
2 |
|
||
E |
|
|
U 0 |
|
|
|
|
|
|
( )2 = const. |
(11.1) |
|||
|
2m 2 |
|
|
|
|
|
4m |
|
8m 2 |
|
||||
Глава 11. Туннелирование
Далее для упрощения записи формул будем рассматривать одномерные локализованные в пространстве частицы — «игольчатые» состояния. Если имеются поперечные составляющие энергии квантовых флуктуаций, то они «протаскиваются» через барьер без изменения, поскольку являются инвариантами движения. При этом в окончательных формулах учтем наличие поперечных составляющих.
Тогда из (11.1) следует, что J / зависит только от координаты x и
|
J |
|
2 2 |
|
|
2 |
|
2 |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
(U 0 |
E)2m. |
(11.2) |
||
|
|
2 x 2 |
|
|
|||||||||||
|
|
|
4 2 |
x |
|
|
|
||||||||
Для выполнения |
этого |
|
соотношения |
необходимо |
положить |
||||||||||
t (t) x (x).
Тогда уравнение (3.2) можно решить методом разделения переменных в квадратурах. Получаем
|
|
m |
x |
|
|
x |
|
J (x,t) exp(t / ) J 0 |
|
||
|
|
0 |
(x)dx . (11.3)
Здесь m / P12 , x = 0 — положение передней стенки барьера, на которой плотность потока вероятности частиц равна J0. Решение уравнения (3.2) можно записать и в другом виде:
|
. |
1 |
|
x |
|
1 |
|
m |
x |
|
J (x,t) J 0 |
exp 0 |
t m x (x)dx 3, J x (x) J 0 |
|
x (x)dx. (11.4) |
||||||
|
|
|||||||||
|
0 |
|
0 |
J x (x) |
3 |
|
|
0 |
||
|
/ |
|
|
|
2 |
|
|
|||
В том, что это одинаковые решения можно убедиться дифференцированием их по координате x. Из (11.4) следует, что поток плотности вероятности J(x,t) J 0 переносится по закону:
x |
|
|
t m x (x) dx. |
(11.5) |
|
0 |
J x (x) |
|
|
|
|
После подстановки решения (11.3) в уравнение (11.1) получаем систему уравнений:
|
J 2(x) |
|
2 |
d |
x |
|
2 |
2 |
|
d2 |
x |
|
|||
E U 0 |
x |
|
|
|
|
|
|
|
|
|
|
, |
|||
|
|
|
|
|
|
4m x |
|
|
|
||||||
|
2m x2(x) 8m x2 |
dx |
|
|
dx 2 |
||||||||||
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
(11.6) |
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
J x (x) J 0 |
|
x (x)dx. |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Глава 11. Туннелирование
Запишем эту систему в другом виде:
2J x(((J x( (J x(()2 42(J x( )2 52J x2 0, x (x) J x( , m
где 4 2 |
8m |
(U 0 E), 52 4P4 |
/ 4 |
(штрихом обозначены произ- |
|
|
|||||
|
2 |
1 |
|
|
|
водные по x). Решение J x |
найдено в виде: |
||||
J x (x) A exp(−x),
где [4 2 / 2 
44 / 4 52 ]1/2 , А — постоянная.
Далее будем вычислять туннельный коэффициент прохождения барьера D. Для этого можно использовать решения для плотности вероятности, полученные ранее, в областях перед барьером и после прохождения барьера. Присвоим падающим на барьер частицам, как и прежде индекс 1, отраженным — индекс 3, туннелирующим частицам в области барьера — индекс 2, частицам, прошедшим барьер, — индекс 4. Имеем:
1 10 cos2 (tP1 / m x1 x / x1),
3
2(x,t)
4
30 cos2( (tP3 / m x3 x / x3) 03),
|
t |
|
|
t |
|
|||
20 |
exp |
|
x , J 2 |
(x,t) J 20 |
exp |
|
x , |
|
|
|
|||||||
|
|
|
|
|
|
|||
40 cos2( (tP4 / m x4 x / x4) 04).
Законы сохранения энергии для свободных частиц можно представить в виде:
E |
P 2 |
|
2 |
2 2 |
, |
|||
1 |
10 |
|
||||||
2m |
8m |
|||||||
|
|
|
|
|
||||
E |
|
|
P 2 |
|
2 |
2 2 |
, |
|
|
|
3 |
30 |
|
||||
|
|
2m |
|
|||||
|
|
|
|
8m |
|
|||
E |
|
P 2 |
|
2 |
2 2 |
. |
||
|
4 |
40 |
|
|||||
|
2m |
|
||||||
|
|
|
|
8m |
|
|||
Граничные условия запишем, как и ранее, в виде:
P1 1(t0, x0) P3 3(t0, x0) J 2(t0, x0) J 2(t0, x0 a) P4 4(t0, x0 a),
1(t0, x0) 3(t0, x0) 2(t0, x0) 2(t0, x0 a) 4(t0, x0 a).
Полагаем далее x0 0, t0 0, поскольку в решении для туннелирующих частиц передняя стенка барьера привязана к координате х0 = 0.
Глава 11. Туннелирование
Решая выписанную систему уравнений совместно с соотношениями для законов сохранения энергии свободных частиц, можно вычислить D:
D |
P4 4 |
(a) |
. |
(11.7) |
|
P1 10 |
|||||
|
|
|
|||
Из граничных условий следует:
4(a) 2(a) 20e a ( 10 30)e a.
Тогда D равно:
D P4( 10 30) exp( a).
P1 10
При a % 0, D % 1, 30 % 0 должно быть P4 = P1. Тогда из закона сохранения энергии следует: 40 10 . Из решения для плотности
вероятности 4( a) 10 cos 42 .
Выпишем окончательные формулы для расчета коэффициента туннельного прохождения барьера D:
D |
10 30 exp( a ) cos 62. |
(11.8) |
|
10 |
|
|
4 |
|
Из формулы (11.8) можно видеть, что частицы, протуннелировавшие через барьер, сохраняют энергию, импульс, а волна плотности вероятности меняет фазу.
Из закона сохранения энергии отраженных частиц и определения R P3 30 / P1 10 находим выражение для 30 :
|
|
|
1 . |
|
|
|
11/ 2 |
|
||
30 |
|
|
2m(E ) |
4m2(E )2 27 (1 |
D)2 P 2 2 |
, (11.9) |
||||
|
|
|
0 |
3 |
||||||
|
|
|
||||||||
|
|
7 |
2 / |
|
|
1 10 |
2 |
|
||
|
. |
|
|
|
11/ 2 |
|
|
4 |
/ 2 |
1 / 4 52 |
/ 44 |
, |
|||
01 |
3 |
||||||
|
/ |
|
|
|
2 |
|
7 / 2.
2 0 E), 52 4P14 / 4.
Спомощью формул (11.8) и (11.9) вычисляется коэффициент прозрачности барьера, если заданы E, P1 , 10 . Прозрачность барье-
ра в классическом пределе должна равняться нулю, действительно, если устремить формально % 0, получаем % & и D % 0.8m
Глава 11. Туннелирование
Традиционная формула для рассматриваемого барьера (см. [1], с. 104) имеет вид:
D |
|
|
|
4k2 |
48 |
|
|
|
|
P / , 48 (2m(U |
|
E |
|
|
|
|
|
|
|
1 |
|
|
|
|
, k |
0 |
k |
) / 2)1/ 2 |
. (11.10) |
||||
|
|
|
|
|
|
|
|
|
||||||||
0 |
(k2 |
4 )2 sh2a4 |
0 |
4k2 |
4 |
1 |
1 |
|
|
|
||||||
|
1 |
|
8 |
|
|
|
1 |
8 |
|
|
|
|
|
|
|
|
|
|
|
|
P 2 |
|
|
|
|
|
|
|
|
|
|
|
|
Здесь Ek |
|
1 |
— кинетическая энергия частицы. Сравнивая фор- |
|||||||||||||
|
||||||||||||||||
|
|
|
|
2m |
|
|
|
|
|
|
|
|
|
|
|
|
мулы (11.10), (11.8) и (11.9), можно видеть, что туннелирование частиц с пространственно структурированной плотностью вероятности отличается от традиционного. Предэкспоненциальный множитель в формуле (11.8) меняется от 1 до 2. В показатель экспоненты входит параметр 4, который запишем виде, удобном для сравнения:
|
|
|
|
|
|
2 |
|
|
|
4 2 2m U |
E |
k |
|
10 |
|
2 2 . |
|||
|
|||||||||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
8m |
|
|
|
Можно видеть, что 4 ! 24 0 , однако в формулу для входит па- |
|||||
раметр 5, который равен |
|
|
|||
|
52 |
|
Ek2 |
|
|
44 |
4(U 0 E) |
||||
|
|||||
и который несколько увеличивает коэффициент . В связи с этим коэффициенты прохождения барьера необходимо сравнивать для конкретных случаев. Сравним эти формулы для случая, когда как обычно a4 0 1. Тогда из (11.10) имеем:
D |
16k2 |
48 |
|
||
|
1 |
|
e 2a40 . |
(11.11) |
|
|
|
|
|||
0 |
(k2 |
|
48)2 |
|
|
|
|
||||
|
1 |
|
|
|
|
Проведем численное сравнение коэффициентов прохождения. Рассмотрим «игольчатые» состояния, когда 0 и 52 / 44 !! 1, тогда 4 24 0 и показатели экспонент совпадают в формулах для D и D0. Тогда различие в коэффициентах прохождения будет определяться отношением предэкспоненциальных множителей.
Пусть k12 / 4 20 Ek / (U 0 Ek ) 1, тогда D / D0 9 0,5 ! 1. Ситуация радикально меняется, когда P1 k1 % 0 и туннелирование воз-
можно в основном за счет энергии квантовых флуктуаций частиц, которые не учитывает традиционный подход, тогда:
D / D0 9 420 / 8k12 1.
Таким образом, коэффициенты туннельного прохождения частиц с пространственно структурированной плотностью вероятно-
Глава 11. Туннелирование
сти практически всегда отличаются от значений, вычисляемых по традиционным формулам.
Время туннелирования частиц с пространственно структурированной плотностью вероятности имеет вполне определенную величину и в соответствии с формулой (11.5) и m / P12 равно:
t a 0 / P1, |
(11.12) |
где 0 — время пролета барьера шириной а с импульсом Р1:
0 am.
P1
Можно видеть, что туннелирование может быть медленным процессом, когда импульс падающей на барьер частицы мал.
Литература
1.Ландау Л.Д., Лифшиц Е.М. Квантовая механика. Нерелятивистская теория. — М.: Наука, 1974, с. 103.
ЗАКЛЮЧЕНИЕ
Квазигидродинамическое представление для инфинитного движения квантовых частиц с физическими переменными позволяет получить новые решения, меняющие представления о них. «Цена» физического вида исходной системы уравнений — нелинейность одного из них. Для теоретической физики значение аналитических решений квантовых задач имеет неоценимое значение. Однако в век компьютерных технологий решение квантовых задач в квазигидродинамическом представлении не должно представлять особого труда, тем более, что транспортные задачи в конкретных квантовых приборах в силу сложности решаются с применением компьютеров.
Квантовые частицы, совершающие инфинитное движение, помимо классической кинетической энергии всегда обладают вполне определенной квантовой энергией — энергией квантовых флуктуаций. При таком подходе устраняются противоречия традиционной теории для транспортных явлений — все новые квантовые формулы и уравнения имеют классический предел. В традиционной квантовой механике этим представлениям наиболее близка идея волновых пакетов для свободных частиц, поскольку в их основе лежит представление о флуктуациях импульса частиц. Однако энергия этих флуктуаций, которую переносят свободные частицы, помимо всего прочего не учитывалась в законе сохранения энергии.
Формулы для коэффициентов прохождения барьеров получены в неявном виде. Формально это связано с нелинейностью одного из квантовых уравнений движения. По существу процесс прохождения барьеров является самосогласованным для распределения плотности вероятности — меняется статистическое волновое поле частицы при сохранении инвариантов движения. В приведенных традиционных решениях квантовых задач это обстоятельство не учитывается.
Учет пространственного структурирования плотности вероятности свободных квантовых частиц дает, на наш взгляд, более правильные и нередко новые соотношения для транспортных явлений. В конечном счете, использование квазигидродинамического представления оправдано, если получены новые результаты, которые подтверждаются экспериментально или могут иметь экспериментальное подтверждение. В частности, экспериментальное доказательство резонансной накачки квантовой составляющей энергии инфинитного движения частиц открывает принципиально новый подход к решению некоторых прикладных проблем, явится еще одним доказательством волновой природы инфинитного движения квантовых частиц.
Приложение 1. Вывод квантовых уравнений движения в квазигидродинамическом представлении
ПРИЛОЖЕНИЕ 1
ВЫВОД КВАНТОВЫХ УРАВНЕНИЙ ДВИЖЕНИЯ В КВАЗИГИДРОДИНАМИЧЕСКОМ ПРЕДСТАВЛЕНИИ
Выпишем необходимые исходные уравнения Шредингера для дальнейших вычислений:
|
|
|
|
* |
|
|
i |
t |
H , |
i |
t |
H *. |
(П1.1) |
|
|
|
|
|
Выберем оператор Гамильтона для частицы массой m, которая двигается в произвольном потенциальном поле U = U(r,t):
|
|
2 |
U (r,t). |
(П1.2) |
|
H |
|
||||
2m |
|||||
|
|
|
|
Умножим первое уравнение Шредингера на * , а второе на и вычитая из первого уравнения второе, получим:
|
|
|
|
2 |
2 |
||
i |
t |
*: : * |
|
( * *) |
|
div ( * *), |
|
|
|
||||||
|
|
|
2m |
2m |
|||
где * . Введем плотность потока вероятности J/m
J i ( * * ).
m2m
Тогда получаем закон сохранения массы частиц в виде:
m div J 0.
t
Далее вычислим производную:
(J |
/ ) |
|
i |
|
|
|
(ln * ln ) |
|
|
|||||||||||
|
|
|
|
2 |
t |
|
|
|||||||||||||
|
t |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
i |
1 * |
|
|
|
1 |
|
1 |
||||||||||||
|
|
|
|
|
: * *: |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
t |
|
|
|
|
|
||||||||||||
|
2 |
* |
|
|
|
|
t |
|
2 |
|
|
|||||||||
(П1.3)
(П1.4)
Здесь использованы выписанные выше уравнения Шредингера. Вставив явный вид оператора Гамильтона, получаем:
(J / ) |
|
2 |
* * |
|
||
|
U |
|
|
. |
(П1.5) |
|
t |
4m |
|
||||
|
|
|
||||
