работа содержит 13 страниц, 5 рисунков, 7 источников. Актуальность (стр. 2 из 3)
что означает невозможность разрыва связи пары при прохождении барьера, и
где функции
– собственные функции уравнения Шредингера с потенциалом 
, собственные значения этого уравнения 
. Коэффициенты S и R определяют вероятности прохождения и отражения барьера соответственно: 
где импульсы
.Классические модели диффузии в твердом теле предполагают абсолютную непрозрачность внутрикристаллических барьеров при энергии частицы E, меньшей высоты барьера Emax. То есть вероятность прохождения барьера P для одной частицы пропорциональна единичной функции:
При этом диффузию как функцию энергии можно записать в виде: D=D0P(E), где D0 имеет масштаб газовой диффузии [7]. После термодинамического усреднения получаем известную формулу эмпирической зависимости диффузии от температуры:
.В квантовом случае вероятность P(E) нужно заменить на полную вероятность прохождения через барьер:
. Соответственно: 
.Температура
берется в единицах энергии. Пример расчета функции 
приведен в [6]. Дальнейшие расчеты при данных параметрах показали, что превышения молекулярной резонансной диффузии над одночастичной не наблюдается.Расчеты, пример которых упомянут выше, чрезвычайно сложны и продолжительны. В [6] использовались параметры потенциалов твердого тела «по возможности», то есть такие, с которыми еще могла работать программа, написанная для машины с восьмиразрядным процессором. При этом некоторая часть резонансов уже и здесь не могла быть рассчитана точно, и для учета их вкладов использовалась полуэмпирическая модель экстраполяции от точно просчитанных резонансов. Параметры потенциалов реальных твердых тел задают еще более сложную перспективу для расчетов. Имеет смысл для предварительной оценки целесообразности таких расчетов в использовании неких приближений, которые, конечно, огрубляют искомый результат, но, вместе с тем, могут дать хотя бы предварительную картину.
Уравнение (2) соответствует двумерной задаче рассеяния на потенциальном барьере, примерный вид которого приведен на рис. 1.
Рисунок 1 – Суммарный потенциал в задаче на резонансную диффузию. Случай осцилляторного потенциала взаимодействия частиц и гауссового – барьера. Стрелками указаны метастабильные состояния
Для получения предварительных результатов заменим эту задачу на более простую: одномерное рассеяние на двух прямоугольных барьерах. Основными параметрами будут ширины и высоты барьеров, для масштабирования будут применяться преобразования, аналогичные (1). Коэффициент C выбран, по аналогии с [6], как для двухатомной молекулы бериллия.
В качестве одного из учитываемых параметров используется также глубина ямы между двумя барьерами. Этот параметр необходимо учитывать, так как в некоторых случаях, например, когда рассматривается рассеяние на двумерном барьере из двух гауссовых потенциалов, резонансное прохождение при низких энергиях отсутствует. В работе исследуется зависимость для задачи, когда потенциал связи частиц выбран осцилляторный, а барьерный потенциал – гауссовой формы. Для этого случая минимальная энергия, при которой возможно резонансное прохождение, определяется как сумма:
, (3)где
– суммарный потенциал. Точка 
является решением уравнения на поиск минимум потенциальной энергии 
.Формула (3) формируется следующим образом: первое слагаемое – значение потенциала в точке минимума, в двух других выражения под корнем – «частота» двух перпендикулярных осцилляторных потенциалов, которыми можно аппроксимировать поведение суммарного потенциала в точке минимума. В целом два последних слагаемых – выражение для минимальной энергии осциллятора
. Этот параметр является достаточно важным – он определяет, кроме вышеуказанного, в спорных случаях возможность резонансной диффузии в целом. Резонансная диффузия возможна, если только (3) меньше высоты барьера.Расчет в выбранном приближении
При преобразовании из (1) в (2) потенциалы преобразуются следующим образом. Коэффициент преобразования
выбирается, чтобы преобразовать потенциал 
в наиболее простой вид, например, 
. Предполагая, данный случай за образец, получаем коэффициент преобразования, как было сказано выше, в виде 
. Можно заметить, что этот параметр имеет размерность длины, то есть, он может быть переписан в виде 
, где 
.В случае, если, как и в [6], рассматривается молекула бериллия, этот параметр равен 0.16 Å. Потенциал
в (2) будет иметь вид 
, 
,параметр
в данном случае равен 400 К [6]. В (1) множитель 
также появиться перед энергией, и в (2) энергия уже будет выражаться в этих единицах.Задача прохождения атома через тот же барьер, то есть одной частицы с массой
, описывается с помощью уравнения Шредингера 
.Для корректного сравнения это уравнение преобразуется таким образом, чтобы перед энергией снова появился тот же множитель
. Очевидно, это возможно, если использовать преобразование 
.Потенциал, таким образом, будет иметь преобразованный вид:
.Потенциал для одной частицы, таким образом «тоньше» на
, чем потенциал в (2) для двухчастичной задачи.Численная схема и результаты расчетов
Для расчетов использовалась аналитическая формула для вероятности прохождения через два одинаковых прямоугольных барьера
, и аналитическая формула вероятности прохождения одного барьера 
, как функции от энергии. Характерный вид кривой вероятности показан на рис. 2. 
Рисунок 2 – Пример расчета вероятности прохождения двух барьеров. Параметры
, 
и 
Параметрами были ширина барьера
, высота 
и расстояние между барьерами 
. Обе функции интегрировались вместе с весовой функцией 
. Интегрирование производилось численно с помощью программы, составленной в программной среде Fortran.