Поверхность. Рентгеновские, синхротронные и нейтронные исследования, 2022, № 10, стр. 107-112
Классическое и квантовое описания эффекта каналирования как взаимно дополняющие приближения
Н. П. Калашников a, *, А. С. Ольчак a, **
a Национальный исследовательский ядерный университет “МИФИ”
115409 Москва, Россия
* E-mail: kalash@mephi.ru
** E-mail: asolchak@mephi.ru
Поступила в редакцию 28.12.2021
После доработки 18.02.2022
Принята к публикации 21.02.2022
- EDN: HKTZKU
- DOI: 10.31857/S1028096022080088
Аннотация
Работа продолжает серию исследований, посвященную разным аспектам каналирования релятивистских электронов в монокристаллах. Движение заряженной частицы в режиме каналирования удобно рассматривать в так называемой сопутствующей системе отсчета, движущейся вдоль направления каналирования со скоростью, равной продольной компоненте скорости каналированной частицы. В такой системе движение частицы финитно и подобно колебательному движению в случае одномерного потенциала (при плоскостном каналировании) или двумерному финитному движению в центральном поле (при аксиальном каналировании). Движение электронов достаточно больших (релятивистских) энергий можно рассматривать как в квантовом, так и в классическом приближении. При классическом рассмотрении удается достаточно просто, аналитически рассчитать интенсивность возникающего электромагнитного излучения, его спектральные характеристики и даже характерные времена жизни квантовых каналированных состояний и вероятности переходов между ними, что непосредственно в квантовом подходе удается сделать только численно. В настоящей работе метод упрощенного аналитического рассмотрения применен к расчету спектральных характеристик и интенсивности излучения, возникающего как при плоскостном, так и при аксиальном каналировании электронов ультрарелятивистских энергий (до нескольких ГэВ). Показано, что при прохождении ориентированной монокристаллической мишени толщиной несколько миллиметров – это излучение способно привести к конверсии значительной части энергии электронного пучка в гамма-кванты высокой энергии.
ВВЕДЕНИЕ
Движение заряженной частицы в поле кристаллической плоскости (плоскостное каналирование) или атомной цепочки (аксиальное каналирование), а также электромагнитное излучение, возникающее при таком движении, можно теоретически рассматривать как в классическом приближении, так и в квантовом. В литературе можно найти примеры обоих подходов [1–9]. Традиционно считается, что при относительно невысоких энергиях (E < 100 МэВ), когда число квантовых связанных каналированных состояний невелико, следует применять квантовый подход – решать релятивистское уравнение Шредингера для определения спектра состояний, применять квантовую электродинамическую теорию возмущений для расчета матричных элементов радиационных переходов. Для приближенных к реальности моделей усредненных потенциалов выполнить такие расчеты можно только численно, что не способствует ясному пониманию физики таких процессов.
При энергии от сотен МэВ и выше, когда число квантовых состояний становится велико, можно воспользоваться классической (не квантовой) релятивистской механикой и классической электродинамикой. В этом приближении многие вычисления можно проделать аналитически, хотя для “формально усредненных” потенциалов даже классические расчеты иначе как численно не выполнишь. “Формальное усреднение” – это усреднение атомных потенциалов по соответствующей плоскости или оси. Если отталкиваться от известных в литературе приближений атомных потенциалов (потенциал Мольер, приближение Хартри–Фока и других), то после формального усреднения для осевого потенциала можно получить функции, в лучшем случае выражаемые через сложные логарифмические (так называемый стандартный потенциал Линдхарда [1]) или специальные (Макдональда [2]) функции, работать с которыми можно только численно.
С другой стороны, точность известных экспериментов по измерению спектров излучения каналированных частиц высоких энергий, как правило, не настолько велика, чтобы можно было заметить отличия близких по форме функций потенциальной энергии от координат. Расчеты в разных приближениях и с разными модельными потенциалами дают сходные результаты и качественно сходное совпадение с экспериментами [4, 9, 10]. В этой связи возникает вопрос: если численный расчет не позволяет качественно улучшить совпадение результата с экспериментом, не лучше ли провести расчет хотя и приближенно, но физически прозрачно и понятно, т.е. аналитически? Численный расчет – это, безусловно, мощнейший математический инструмент, но понять и почувствовать физический смысл рассчитываемого эффекта он не помогает.
В настоящей работе авторы предлагают примеры расчета спектра состояний и интенсивности излучения каналированной частицы, сознательно выполненные для максимально простых (математически) модельных потенциалов с использованием классического и упрощенного квантового подходов. Целью является демонстрация возможности и желательности использования приближенных методов при изучении как самого явления каналирования, так и излучения при каналировании.
Для удобства движение частиц рассматривается как в лабораторной системе отсчета, так и в так называемой сопутствующей системе отсчета [11–13], движущейся вдоль направления каналирования со скоростью, равной продольной компоненте скорости каналированной частицы. В сопутствующей системе отсчета движение частицы финитно и подобно колебательному движению в случае одномерного потенциала (при плоскостном каналировании) или двумерному финитному движению по орбитам в центральном поле (при аксиальном каналировании). Используя удобные приближенные модели усредненного потенциала, можно достаточно просто, аналитически рассчитать спектральные характеристики возникающего электромагнитного излучения, его интенсивность и характерные времена жизни квантовых каналированных состояний. Покажем это на нескольких примерах.
ПЛОСКОСТНОЕ КАНАЛИРОВАНИЕ – ОДНОМЕРНАЯ МОДЕЛЬ АТОМА
В сопутствующей системе отсчета каналированные в плоскостном канале частицы совершают финитные колебания между соседними ионными плоскостями (если они заряжены положительно) или вблизи одной из таких плоскостей (если они заряжены отрицательно). По сути, электрон в режиме плоскостного каналирования представляет собой одномерную модель атома. Аналитически несложно решить задачу о движении частиц в случае параболического потенциала, имеющего минимум посередине между соседними ионными плоскостями и достигающего максимумов точно на плоскостях (вполне реалистичное приближение для каналирования положительно заряженных частиц – позитронов или протонов):
(1)
$U = {{k{{x}^{2}}} \mathord{\left/ {\vphantom {{k{{x}^{2}}} 2}} \right. \kern-0em} 2} = {{4{{U}_{0}}{{x}^{2}}} \mathord{\left/ {\vphantom {{4{{U}_{0}}{{x}^{2}}} {{{d}^{2}},}}} \right. \kern-0em} {{{d}^{2}},}}$(2)
${{\omega }_{{{\text{кл}}}}} = c\sqrt {{k \mathord{\left/ {\vphantom {k E}} \right. \kern-0em} E}} = ({{2с} \mathord{\left/ {\vphantom {{2с} d}} \right. \kern-0em} d}){{({{2{{U}_{0}}} \mathord{\left/ {\vphantom {{2{{U}_{0}}} Е}} \right. \kern-0em} Е})}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}},$В квантовом приближении уровни поперечной энергии связанного одномерного движения можно определить из одномерного (в сопутствующей системе отсчета) релятивистского уравнения Шредингера с потенциалом атомной плоскости U(x), умноженным на лоренц-фактор:
(3)
${{\hbar {{с}^{2}}{{d}^{2}}\psi (x)} \mathord{\left/ {\vphantom {{\hbar {{с}^{2}}{{d}^{2}}\psi (x)} {d{{x}^{2}}}}} \right. \kern-0em} {d{{x}^{2}}}} + 2E({{E}_{{x,n}}} - U(x))\psi (x) = 0,$(4)
$\oint {pdx} = 2\pi \hbar (n + {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}),\,\,\,\,n = 0,1,2,3, \ldots ,$(5)
$N \sim {{{{U}_{0}}} \mathord{\left/ {\vphantom {{{{U}_{0}}} {\hbar {{\omega }_{{{\text{кл}}}}}}}} \right. \kern-0em} {\hbar {{\omega }_{{{\text{кл}}}}}}} \sim ({{{{E{{U}_{0}}} \mathord{\left/ {\vphantom {{E{{U}_{0}}} {2)}}} \right. \kern-0em} {2)}}}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}({d \mathord{\left/ {\vphantom {d {2с\hbar }}} \right. \kern-0em} {2с\hbar }}),$(6)
$U(x) \approx - {\kern 1pt} {{U}_{0}}(1 - {x \mathord{\left/ {\vphantom {x b}} \right. \kern-0em} b})\,\,{\text{при}}\,\,х < b \sim {d \mathord{\left/ {\vphantom {d 4}} \right. \kern-0em} 4},$УПРОЩЕННЫЙ ЛИНЕЙНЫЙ ПОТЕНЦИАЛ И КЛАССИЧЕСКОЕ ПРИБЛИЖЕНИЕ В ТЕОРИИ ИЗЛУЧЕНИЯ
Учитывая слабую зависимость спектральных характеристик связанных состояний частиц в плоскостных каналах от геометрии усредненного потенциала, для оценки интенсивности электромагнитного излучения, испускаемого каналированным электроном, воспользуемся потенциалом (6), достаточно реалистичным и позволяющим выполнить такой расчет аналитически. Напомним, что “точное” квантовое решение этой задачи требует расчета волновых функций c помощью релятивистского уравнения Шредингера (3) и вычисления матричных элементов для переходов между уровнями с найденными функциями, что можно сделать только численно даже с самыми простыми потенциалами. Однако, если воспользоваться несколько упрощенным квантово-классическим подходом, то эта задача решается аналитически и очень просто. Заметим, что в случае потенциала (6) электрон (в классическом приближении) движется под действием практически постоянной по величине силы, испытывая (в сопутствующей системе отсчета) ускорение:
(7)
$\left| w \right| = {{\left| {{{дU} \mathord{\left/ {\vphantom {{дU} {дх}}} \right. \kern-0em} {дх}}} \right|} \mathord{\left/ {\vphantom {{\left| {{{дU} \mathord{\left/ {\vphantom {{дU} {дх}}} \right. \kern-0em} {дх}}} \right|} Е}} \right. \kern-0em} Е} \approx {{{{U}_{0}}} \mathord{\left/ {\vphantom {{{{U}_{0}}} {bm.}}} \right. \kern-0em} {bm.}}$(8)
$I = ({{2k{{e}^{2}}} \mathord{\left/ {\vphantom {{2k{{e}^{2}}} {3{{c}^{3}}}}} \right. \kern-0em} {3{{c}^{3}}}}){{w}^{2}} = {{2k{{e}^{2}}U_{0}^{2}} \mathord{\left/ {\vphantom {{2k{{e}^{2}}U_{0}^{2}} {3{{b}^{2}}{{m}^{2}}{{c}^{3}},}}} \right. \kern-0em} {3{{b}^{2}}{{m}^{2}}{{c}^{3}},}}$(9)
$t \sim {{{{U}_{0}}} \mathord{\left/ {\vphantom {{{{U}_{0}}} I}} \right. \kern-0em} I} \sim {{{{b}^{2}}{{m}^{2}}{{c}^{3}}} \mathord{\left/ {\vphantom {{{{b}^{2}}{{m}^{2}}{{c}^{3}}} {k{{e}^{2}}{{U}_{0}}.}}} \right. \kern-0em} {k{{e}^{2}}{{U}_{0}}.}}$В процессе потери поперечной энергии электрон будет терять и энергию полную, причем при каждом переходе с изменением поперечной энергии на ∆ε будет испускаться фотон в направлении вперед (эффективные углы излучения θ ~ mc2/E) c энергией в лабораторной системе отсчета, сильно увеличенной за счет эффекта Доплера [2, 3, 8, 9, 13, 17, 18]:
(10)
$\begin{gathered} h\nu \sim {{({E \mathord{\left/ {\vphantom {E {m{{c}^{2}}}}} \right. \kern-0em} {m{{c}^{2}}}})}^{2}}\Delta \varepsilon ~\sim {{({E \mathord{\left/ {\vphantom {E {m{{c}^{2}}}}} \right. \kern-0em} {m{{c}^{2}}}})}^{2}}({{{{U}_{0}}} \mathord{\left/ {\vphantom {{{{U}_{0}}} N}} \right. \kern-0em} N}) \sim \\ \sim 4hс{{({E \mathord{\left/ {\vphantom {E {m{{c}^{2}}}}} \right. \kern-0em} {m{{c}^{2}}}})}^{2}}{{({{{{U}_{0}}} \mathord{\left/ {\vphantom {{{{U}_{0}}} Е}} \right. \kern-0em} Е})}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}. \\ \end{gathered} $(11)
$\Delta Е \sim {{({E \mathord{\left/ {\vphantom {E {m{{c}^{2}}}}} \right. \kern-0em} {m{{c}^{2}}}})}^{2}}{{U}_{0}}.$АКСИАЛЬНОЕ КАНАЛИРОВАНИЕ ЭЛЕКТРОНОВ – ДВУМЕРНАЯ МОДЕЛЬ АТОМА
В сопутствующей системе отсчета движение электрона вокруг притягивающей его ионной цепочки подобно двумерному финитному движению по орбитам в центральном поле. В литературе известны весьма разные модельные потенциалы, использовавшиеся для описания такого движения. Исторически первой моделью усредненного потенциала атомной цепочки был так называемый стандартный потенциал Линдхарда [1]:
(12)
$U(\rho ){\text{ }} = - \frac{{Z{{e}^{2}}}}{d}{\text{ln}}\left[ {\frac{{3R_{{{\text{T}} - {\text{F}}}}^{2}}}{{{{\rho }^{2}}}}} \right] + 1,$В [2] выполнено точное усреднение потенциала цепочки атомов, описываемых экспоненциально экранированным кулоновским потенциалом, приводящее к специальным функциям Макдональда, дальнейшая работа с которыми возможна только численно. В [4, 8, 19, 20] предлагали аппроксимировать непрерывный потенциал атомной цепочки двумерным кулоновским потенциалом:
(13)
$U\left( \rho \right) = - {\kern 1pt} c\frac{{Z{{e}^{2}}{{R}_{{{\text{T}} - {\text{F}}}}}{\text{\;\;\;}}}}{{\rho d}}~,$(14)
$\begin{gathered} {{\psi }_{{nm}}}(\rho ,\varphi ) = {\text{exp}}\left( {im\varphi } \right)\frac{1}{{\sqrt \rho }} \times \\ \times \,\,{\text{exp}}\left( { - {\kern 1pt} \gamma \frac{\rho }{a}} \right){{\left( {\frac{\rho }{a}} \right)}^{{\left| m \right|\,{\kern 1pt} + \,{\kern 1pt} \frac{1}{2}}}}\sum\limits_{k{\kern 1pt} \, = \,{\kern 1pt} 0}^n {{{a}_{k}}{{{\left( {\frac{\rho }{a}} \right)}}^{k}}} . \\ \end{gathered} $(15)
$U(\rho ) = \left\{ {\begin{array}{*{20}{c}} { - {\kern 1pt} {{U}_{0}}~}&{0 \leqslant \rho \leqslant {{R}_{{{\text{T}} - {\text{F}}}}}} \\ {~0}&{~\rho > {{R}_{{{\text{T}} - {\text{F}}}}},} \end{array}} \right.$(17)
$ - \frac{{{{\hbar }^{2}}}}{{2\mu }}\left( {\frac{{{{d}^{2}}R}}{{d{{\rho }^{2}}}} + \frac{1}{\rho }\frac{{dR}}{{d\rho }} - \frac{{{{m}^{2}}}}{{{{\rho }^{2}}}}R} \right) + U(\rho )R = \varepsilon R.$(18)
${{\varepsilon }_{{n,m}}} = {{{{\pi }^{2}}{{n}^{2}}} \mathord{\left/ {\vphantom {{{{\pi }^{2}}{{n}^{2}}} {2ER_{{{\text{T}} - {\text{F}}}}^{2}}}} \right. \kern-0em} {2ER_{{{\text{T}} - {\text{F}}}}^{2}}} + {{{{m}^{2}}} \mathord{\left/ {\vphantom {{{{m}^{2}}} {2ER_{{{\text{T}} - {\text{F}}}}^{2} - {{U}_{0}}}}} \right. \kern-0em} {2ER_{{{\text{T}} - {\text{F}}}}^{2} - {{U}_{0}}}}.$ИЗЛУЧЕНИЕ ПРИ АКСИАЛЬНОМ КАНАЛИРОВАНИИ. КЛАССИЧЕСКАЯ ОЦЕНКА ИНТЕНСИВНОСТИ
Для оценки интенсивности излучения электрона при аксиальном каналировании удобно воспользоваться упрощенным, но сохраняющим физическую суть рассмотрением, аналогичным случаю плоскостного каналирования. Для выбора самой простой, но реалистичной модели усредненного потенциала атомной цепочки приведем несколько очевидных физических соображений.
Потенциал каждой ионной цепочки для каналированного электрона является притягивающим и убывает при приближении к цепочке. Вблизи каждой цепочки формально усредненный потенциал покоящихся строго на оси цепочки ионов с точечными ядрами в центре “расходится” по логарифмическому закону: U(ρ) → –∞ при ρ → 0 [1, 2]. Понятно, однако, что это не физическая расходимость. Реальный потенциал всегда ограничен (конечен), если учесть хотя бы неустранимое тепловое движение ионов и дополнительно усреднить по нему. Потенциалы соседних цепочек перекрываются, благодаря чему устраняется нефизическая логарифмическая расходимость [1, 2] формально усредненных потенциалов ионных цепочек на большом удалении от осей.
Форма реального потенциала “пучка” параллельных ионных цепочек – сложная искривленная поверхность с коническими углублениями, имеющими вершины на осях цепочек. Недалеко от оси цепочки потенциал будет расти практически линейно, как и потенциал (6) в плоскостном случае, но с заменой переменной х на ρ: U0(1 – ρ/b), где U0 – эффективная глубина потенциальной ямы, а b ≤ d/2 – параметр, который при расчете можно подобрать так, чтобы обеспечить лучшее совпадение с реальным.
В случае аксиального, линейно растущего с удалением от оси потенциала частица (в сопутствующей системе отсчета) может, в частности, двигаться по круговым финитным траекториям разного радиуса, но всегда с одинаковым центростремительным ускорением, которое обеспечивает этот потенциал: |w| = U0/bm. Соответственно, это движение будет сопровождаться электромагнитным излучением постоянной интенсивности, как и в плоскостном случае, определяемой той же формулой (8). Спектральные характеристики (10), характерное время потери поперечной энергии (9), характерная полная потеря энергии (11) будут такими же, как и в плоскостном случае, с заменой глубины плоскостного канала на глубину осевого и заменой параметра b. Численная оценка длины, на которой электрон опустится на дно потенциальной ямы (несколько миллиметров), тоже остается прежней. Интересно, что расчет с потенциалом кулоновского типа (13), естественно, дает оценку времени падения электрона на ось того же порядка (~3 × 10–11 c), что и оценка времени падения электрона на ядро в модели атома Резерфорда [22], и обе они близки к оценке (9).
ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ И ЗАКЛЮЧЕНИЕ
Применение упрощенного, квантово-классического подхода для теоретического рассмотрения задачи об интенсивности и спектральных характеристиках электромагнитного излучения, возникающего при каналировании электронов в монокристаллах, позволяет аналитически рассчитать основные характеристики излучения для электронов с энергией Е ~108–109 эВ, когда количество уровней поперечной энергии в случае усредненного потенциала велико и можно уже использовать классическое (не квантовое) приближение, одновременно применяя нерелятивистское приближение для описания поперечного движения электрона в канале.
Проведенный анализ позволяет установить, что спектральные и энергетические характеристики излучения, возникающего при плоскостном и при аксиальном каналировании в рассматриваемом диапазоне значений энергии, качественно близки. Интегральная интенсивность излучения растет пропорционально квадрату полной начальной энергии электрона. При энергии электронов несколько ГэВ серия радиационных переходов с верхних связанных уровней каналированного движения на нижние способна привести к потере энергии, сравнимой с начальной энергией электрона, на монокристалле толщиной несколько миллиметров (l ≤ 1 см). Оценка толщины кристалла, способного привести к почти полной конверсии энергии электрона в энергию жесткого гамма-излучения (l ~ 1 см), слабо зависит от типа кристалла и ориентации кристаллической решетки. Излучение, сопровождающее практически любой режим каналирования, должно иметь сходные количественные характеристики.
Список литературы
Линдхард Й. // УФН. 1969. Т. 99. № 2. С. 249.
Калашников Н.П. Когерентные взаимодействия заряженных частиц в монокристаллах. М.: Атомиздат, 1981. 224 с.
Барышевский В.Г. Каналирование, излучение и реакции в кристаллах при высоких энергиях. Минск: Изд. БГУ им. В.И. Ленина, 1982. 256 с.
Воробьев С.А. Каналирование электронных пучков. М.: Энергоатомиздат, 1984. 96 с.
Базылев В.А., Живаго Н.К. Излучение быстрых частиц в веществе и во внешних полях. М.: Наука, 1987. 272 с.
Байер В.Н., Катков В.М., Страховенко В.М. Электромагнитные процессы при высокой энергии в ориентированных кристаллах. Новосибирск: Наука, 1989. 400 с.
Ахиезер А.И., Шульга Н.Ф. Электродинамика высоких энергий в веществе. М.: Наука, 1993. 344 с.
Рябов В.А. Эффект каналирования. М.: Энергоатомиздат, 1994. 240 с.
Оцуки Е.-Х. Взаимодействие заряженных частиц с твердыми телами. М.: Мир, 1985. 280 с.
Andersen J.U., Bonderup E., Loegsgaard E. et al. // Nucl. Instrum. Methods. 1982. V. 194. P. 209.
Kalashnikov N.P., Olchak A.S. // Nucl. Instrum. Methods Phys. Res. B. 2015. V. 355. P. 121.
Калашников Н.П., Ольчак А.С. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2018. № 4. С. 25.
Загайнов В.А., Калашников Н.П., Ольчак А.С. // Поверхность. Рентген., синхротр. и нейтрон. исслед. 2020. № 3. С. 109.
Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. III Квантовая механика. Нерелятивистская теория. М.: Наука, 1977. 768 с.
Калашников Н.П., Ольчак А.С. Взаимодействие ядерных излучений с монокристаллами. / Ред. Рязанов М.И. М.: МИФИ, 1979. 58 с.
Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. II. Теория поля. М.: Наука, 1973. 504 с.
Калашников Н.П., Мамонов М.Н., Ольчак А.С., Стриханов М.Н. // Физика твердого тела. 1983. Т. 25. № 1. С. 190.
Kalashnikov N.P., Olchak A.S., Khangulian E.V. // Nucl. Instrum. Methods Phys. Res. B. 2013. V. 309. P. 67.
Kaplin V.V., Plotnikov S.V., Tsekhanovsky I.A. // Phys. Lett. A. 1975. V. 54. № 6. P. 447.
Рябов В.А. // ЖЭТФ. 1982. Т. 8. С. 1176.
Каплин В.В., Воробьев С.А. // ЖЭТФ. 1977. Т. 73. № 2. С. 583.
Rutherford E. // Philosoph. Magazine. Ser. 6. V. 21. Iss. 125. P. 669.
Дополнительные материалы отсутствуют.
Инструменты
Поверхность. Рентгеновские, синхротронные и нейтронные исследования