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