§ 3.3. Энергетическая зависимость фаз рассеяния при низких энергиях
При низких энергиях, когда импульс падающих частиц много меньше
обратного радиуса взаимодействия (kd << 1), зависимость фаз рассеяния от
энергии носит универсальный характер. Чтобы установить эту зависимость,
обратимся сначала к формулам (3.15) и (3.21). Исключив из них ƒℓ(k),
получаем
|
(3.29) |
Это строгое соотношение, связывающее фазу рассеяния ℓ-й парциальной
волны с соответствующей радиальной функцией Rℓ(r) во внутренней области
r
< d, аналогично соотношению (1.34) в общей формулировке задачи
потенциального рассеяния. Рассмотрим правую часть соотношения (3.29) при
kd << 1. Для этого нам понадобятся предельные выражения сферических
функций Бесселя и сферических функций Неймана при малых значениях
аргумента:
|
(3.30) |
|
(3.31) |
Рассматривая с учетом этих свойств функций jℓ(х) и nℓ(х) соотношение (3.11), видим, что функция Грина
(E,r,r') при
k → 0 теряет зависимость от энергии (импульса) частицы:
|
(3.32) |
Обращаясь теперь к интегральному уравнению (3.10), получаем, что при k → 0 решение этого уравнения Rℓ(r) ведет себя в зависимости от к так
же, как сферическая функция Бесселя jℓ(kr), т.е. пропорционально kℓ. Следовательно, из (3.29) вытекает
δ(k)|k→0 ~ k2ℓ+1 ~ Eℓ+1/2 . |
(3.33) |
Это и есть искомая зависимость фаз рассеяния от энергии (импульса)
частиц при малых энергиях. Она носит универсальный характер и присуща
рассеянию на любом потенциале конечного радиуса.
§ 3.4. Методы вычисления фаз рассеяния
Точное решение задачи рассеяния с целью вычисления фаз рассеяния
возможно лишь для отдельных искусственно придуманных потенциалов. На
практике, когда приходится иметь дело с реалистическими потенциалами,
фазы рассеяния всегда вычисляются приближенно, что связано либо с
использованием тех или иных физических аппроксимаций, либо с проведением
численного счета. Мы познакомимся с методами и того и другого родов.
а) Метод решения радиального уравнения Шредингера
Перепишем уравнение (3.3) в виде
|
(3.34) |
Здесь, как и при решении задачи о связанных состояниях частицы [1, с.114], функции
uℓ(r) связаны с Rℓ(r) соотношением
и при малых r ведут себя согласно степенному закону:
Поведение
uℓ(r) при асимптотически больших г получаем из (3.23):
|
(3.37) |
В общем случае уравнение (3.34) интегрируется численно, начиная от
r =
0; для «разгонки» численного интегрирования удобно использовать
свойство (3.36). Мы не будем вдаваться в чисто вычислительные аспекты
этой процедуры. Для нас важно, что задачей интегрирования является
построение всего «профиля» функции uℓ(r) с выходом при r → ∞ на
асимптотику (3.37); при этом должно быть соблюдено условие непрерывности
функции uℓ(r) и ее первой производной u'ℓ(r) во всей области 0 <
r
< ∞.
При рассеянии на потенциале конечного радиуса (пусть d − размеры
области взаимодействия) в задаче имеются два масштаба измерения
расстояний r: параметр d и длина волны
= 1/k. Соотношение (3.37), как и
соотношение (3.23), из которого оно получено, справедливы, если
выполняются два неравенства: kr >> 1, r > d. Если энергия частицы Е и
размеры области взаимодействия d таковы, что условие kr >> 1 не
выполняется, то асимптотическое поведение (3.37) начинается далеко за
границей области взаимодействия (r >> d). В этом случае не рационально
доводить численное интегрирование до асимптотически больших значений r >>
d, поскольку уже при r > d частица движется свободно и уравнение
(3.3) переходит в уравнение движения свободной частицы:
|
(3.38) |
При Е > 0 его линейно независимыми решениями являются сферическая
функция Бесселя jℓ(kr) и сферическая функция Неймана nℓ(kr), а общее
решение уравнения (3.34) при r > d имеет вид
Rℓ(r) = uℓ(r)/r = Aℓjℓ(kr) + Bℓnℓ(kr); r>d . |
(3.39) |
Вместо констант интегрирования Аℓ и Вℓ удобно выбрать две другие:
Aℓ = Cℓ cos δℓ; Вℓ = −Cℓ sin
δℓ. |
(3.40) |
Следовательно,
Rℓ(r) = uℓ(r)/r = Cℓ{cos
δℓ jℓ(kr) − sin δℓ
nℓ(kr)}; r > d . |
(3.41) |
С учетом (3.12) и (3.13) отсюда видно, что общее решение (3.39) имеет требуемую асимптотику (3.37).
Таким образом, практическая задача интегрирования уравнения (3.34)
заключается в том, чтобы, начиная с r = 0, выйти в область свободного
движения частицы r > d и «сшить» волновую функцию, найденную для
внутренней области, с функцией (3.41). Пусть
(r)
= /r
− решение
уравнения (3.3), (3.34), удовлетворяющее условию (3.36) и справедливое
во всей внутренней области 0 < r < d. Его «сшивание» с (3.41)
означает:
(d)
= Cℓ{cos δℓ jℓ(kd) − sin
δℓ nℓ(kd)},
d/dr|r=d
= kCℓ{cos δℓ j'ℓ(kd)
− sin δℓ n'ℓ(kd)}, |
(3.42) |
где j'ℓ(x) ≡ djℓ(x)/dx. Введем
логарифмическую производную волновой функции
(r) на границе области взаимодействия:
|
(3.43) |
Исключая из (3.42) константу Cℓ (а вместе с ней и произвольно
выбираемый при интегрировании уравнения (3.3) нормировочный множитель в
(r)), выражаем фазу рассеяния δℓ через логарифмическую производную ƒℓ:
|
(3.44) |
В качестве примера рассчитаем фазы рассеяния частицы на прямоугольной яме:
|
(3.45) |
В этом случае для получения волновой функции Rℓ(r) во внутренней области не требуется численного интегрирования:
(r) = const ·
jℓ(Kr) , |
(3.46) |
где
|
(3.47) |
Отсюда находим логарифмическую производную (3.43):
|
(3.48) |
а далее по формуле (3.44) и фазу рассеяния:
|
(3.49) |
В качестве другого примера рассмотрим рассеяние частицы абсолютно твердой сферой. На поверхности такой сферы (r = d) «внешняя» волновая функция обращается в нуль:
Подставляя сюда (3.41), находим
tg δℓ = jℓ(kd)/nℓ((kd) . |
(3.51) |
В частности, при малых энергиях частицы отсюда получаем
|
(3.52) |
б) Теория возмущений
Если V(r) мало, то радиальную функцию Rℓ(r), входящую в формулу
(3.29) для фазы и удовлетворяющую уравнению (3.10), можно заменить во
внутренней области (r < d) выражением
Rℓ(r) → iℓ(2ℓ + 1)jℓ(kr).
Тогда получаем приближенное выражение для фазы:
|
(3.53) |
а вместе с ним и для парциальной амплитуды рассеяния:
|
(3.54) |
Согласно (3.54), амплитуда рассеяния ƒ(0) вещественна. С точки зрения
оптической теоремы (3.28) такой результат, на первый взгляд,
бессмыслен: мнимая часть амплитуды рассеяния вперед Im ƒ(0) равна нулю, а
полное сечение σ, конечно, нет. Правда, если фазы рассеяния малы, то и
мнимые части парциальных амплитуд Im ƒℓ = (2ℓ + l)k-1 sin2
δℓ, интегральные парциальные сечения
σℓ = (4π/k2)(2ℓ + l) sin2 δℓ квадратичны относительно δℓ; поэтому можно сказать, что при малых δℓ, оптическая
теорема выполняется приближённо − с точностью до членов,
пропорциональных δℓ2.
Таким образом, условие применимости теории возмущений для фаз − это малость фаз:
Приближение (3.53) для фаз при соблюдении требования (3.55) будем
называть борновским приближением для фаз рассеяния. Согласно (3.53) и
(3.55), в этом приближении фазы вычисляются по формуле
|
(3.56) |
Она показывает, в частности, что если V(r) везде отрицательно (потенциал притяжения), то все борновские фазы
положительны; формула (3.23) показывает, что в этом случае
вся осцилляционная картина внешней волновой функции частицы как бы
«подтягивается» в области взаимодействия. При положительном V(r) все
борновские фазы отрицательны; волновая функция частицы в этом случае как
бы «выталкивается» из области взаимодействия.
Выясним условия, при которых соблюдается требование малости фаз
(3.55). Как и в §2.2, рассмотрим два предельных случая: длина волны
частицы =
k-1 много больше и много меньше, чем область взаимодействия.
Если
>> d, т.е. kd
<< 1, то, подставляя в (3.56) аппроксимацию (3.30), получаем
|
(3.57) |
где d − средний радиус области взаимодействия, a |V|
− средняя
амплитуда взаимодействия. Известно, что если сферически симметричный
потенциал притяжения аппроксимируется прямоугольной ямой, то условие
существования в нем связанного состояния − это
|
(3.58) |
Если же в потенциальной яме имеются два-три уровня, то параметр (2μ|V|d2)/ћ2, которому пропорциональна величина
|| составляет несколько единиц. Формула (3.57) показывает, что условие применимости борновского приближения (3.55)
может выполняться для фаз с большими ℓ и не выполняться для фаз с
малыми ℓ. Хуже всего борновское приближение годится для фазы s-волны. С
другой стороны, как видно из (3.25), относительный вклад s-волны в
амплитуду рассеяния меньше всего при самых малых углах рассеяния. Отсюда
следует важный качественный вывод: если энергия частицы и протяженность
потенциала конечного радиуса таковы, что kd << 1, то при рассеянии на
малые углы условия применимости борновского приближения всегда в среднем
более благоприятны, чем при рассеянии на большие углы.
Если
<< d, т. е. kd << 1, то для оценки фаз рассеяния следует подставить в (3.56) аппроксимацию (3.12):
|
(3.59) |
В этом случае условие применимости борновского приближения
|
(3.60) |
не содержит ℓ и полностью совпадает с тем, что для случая kd >> 1 было получено в § 2.2.
|