КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ РЕЗОНАНСНОГО КОГЕРЕНТНОГО ВОЗБУЖДЕНИЯ ИОНОВ 17+Ar С УЧЕТОМ ТОНКОЙ СТРУКТУРЫ ЭНЕРГЕТИЧЕСКИХ УРОВНЕЙ

9. Leawitt R.P. Pressure broadening and shifting in microwave and infrared spectra of molecules of arbitrary symmetry: An irreducible tensor approache // J. Chem. Phys. - 1980. - V. 73. - № 12. -P. 5432-5450.

10. Черкасов М.Р. Формализм квантовомеханического оператора Лиувилля в расчетах релаксационных параметров. - Томск: Изд-во ИОА СО АН СССР, 1975. - 47 c. (препринт № 26).

11. Гроссман В.Э., Броуэлл Э.В., Быков А.Д. и др. Экспериментальное и теоретическое исследование сдвигов линий Н2О давлением N2, О2 и воздуха // Оптика атмосферы. - 1990. - Т. 3. - № 7. - С. 675-690.

12. Быков А.Д., Стройнова В.Н. Сдвиг центров линий поглощения Н2О в диапазоне 13550-13950 см-1 // Оптика атмосферы и океана. - 2006. - Т. 19. - № 1. - С. 31-38.

13. Стройнова В.Н. Расчеты коэффициентов уширения линий поглощения Н2О в диапазоне 13550-13950 см-1 // Оптика атмосферы и океана. - 2005. - Т. 18. - № 9. - С. 820-824.

14. Быков А.Д., Лазарев В.В., Пономарев Ю.Н., Стройнова В.Н. и др. Сдвиги линий поглощения Н2О в полосе v+3 v3, индуциро-

ванные давлением благородных газов // Оптика атмосферы и океана. - 1994. - Т. 4. - № 9. - С. 1207-1219.

15. Варшалович Д.А., Москалев А.Н., Херсонский В.К. Квантовая теория углового момента. - Л.: Наука, 1975. - 439 c.

16. Gray C.G., Gubbins K.E. Theory of molecular fluids. - New York: Oxford University Press, 1984. - 626 p.

17. Mengel M., Jensen Per. A theoretical study of the Stark effect in tri-atomic molecules: Application to Н2О // J. Mol. Spectrosc. - 1995. - V. 169. - P. 73-91.

18. Физические величины / Справочник под ред. И.С. Григорьева,

Е.З. Мельникова. - М.: Энергоатомиздат, 1991. - 891 c.

19. Emelyanov D.S., Stroinova V.N. H2O halfwidths and lineshifts in 3v+v3 band induced by Ar pressure // Proc. SPIE. - 2006. -V. 6580. - 65800H (7 pages).

Поступила 14.12.2007г.

УДК 538.97:539.186:539.184

КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ РЕЗОНАНСНОГО КОГЕРЕНТНОГО ВОЗБУЖДЕНИЯ ИОНОВ 17+Ar С УЧЕТОМ ТОНКОЙ СТРУКТУРЫ ЭНЕРГЕТИЧЕСКИХ УРОВНЕЙ

А.А. Бабаев, Ю.Л. Пивоваров

Томский политехнический университет E-mail: krass58ad@mail.ru

Построена компьютерная модель, позволяющая изучать влияние тонкой структуры на характеристики резонансного когерентного возбуждения релятивистских водородоподобных ионов при плоскостном каналировании. Модель использована для описания резонансного когерентного возбуждения при плоскостном каналировании 390 МэВ/нуклон ионов 17+Ar. Показано, что резонансные кривые характеризуются двумя близко расположенными пиками, обусловленными переходами электрона иона из основного состояния в состояния, соответствующие компонентам тонкой структуры возбужденного состояния. Получено достаточно хорошее согласие с экспериментом.

При плоскостном каналировании водородопо-добных ионов ион движется в периодическом поле плоскостей кристалла. В системе покоя иона на орбитальный электрон действует периодическое во времени электрическое поле. Если частота такого воздействия совпадает с частотой перехода между уровнями энергии орбитального электрона, возможен резонансный переход электрона в соответствующее возбужденное состояние [1] - резонансное когерентное возбуждение (resonant coherent excitation, RCE).

Ранее исследовалось, главным образом, RCE при каналировании нерелятивистских ионов (см., например, [2-5]). Проведены первые эксперименты по наблюдению RCE тяжелых релятивистских ионов. Существуют экспериментальные данные [6] и теоретические оценки [7], указывающие на значительное влияние тонкой структуры уровней энергии орбитального электрода на форму резонансной кривой. Таким образом, представляет интерес полное теоретическое исследование и компьютерное моделирование RCE тяжелых реляти-

вистских ионов с учетом тонкой структуры уровней энергии орбитального электрона.

Рассмотрим движение релятивистского водоро-доподобного иона вдоль кристаллографической плоскости в режиме плоскостного каналирования. Пусть направление оси ОХ прямоугольной системы координат перпендикулярно плоскостям кана-лирования, а направления осей OY и OZ совпадают с направлениями кристаллографических осей в плоскости каналирования. Электрическое поле кристалла в системе покоя иона представляется в виде суммы непрерывного потенциала и гармоник периодического электрического поля с частотами:

(

=lnJv

, cosœ ,sinф k-— +1——

A

где у- релятивистский фактор; v - скорость иона; ay и az - расстояния между атомами в плоскости; ф - угол между скоростью иона и направлением оси OY в плоскости; k и ! - целые числа. Резонансный

переход между г-м и /-м уровнями энергии орбитального электрона наблюдается, если для некоторой гармоники периодического электрического поля выполняется условие

(

Еи = = 2пЛуу

С08ф 81И ф К--г I-

Л

где Е/ - энергия перехода между г-м и /-м уровнями энергии; Й - постоянная Планка.

При вычислениях потенциал атома вещества решетки аппроксимирован потенциалом Мольер. Учитывается потенциал только двух ближайших к траектории иона плоскостей кристалла.

Система уровней энергии орбитального электрона тяжелого иона в отсутствие внешнего поля определяется главным квантовым числом уровня, тонкой структурой уровня и лэмбовским сдвигом. Внешнее электрическое поле плоскостей кристалла приводит к возникновению эффекта Штарка, проявляющегося в сдвиге и расщеплении уровней энергии, причем степень влияния эффекта Штарка на положение уровней энергии зависит от расстояния иона от центра канала. Пусть Иь - гамильтониан системы электрон-ион, не учитывающий спин-орбитальное взаимодействие, эффект Штарка и лэмбовский сдвиг, щь - собственные функции и еь - собственные значения энергии этого гамильтониана. Пусть V - оператор, отвечающий за возникновение лэмбовского сдвига, тонкой структуры уровней энергии электрона иона, действие непрерывного потенциала кристаллических плоскостей на электрон иона.

Волновые функции и собственные значения энергии системы с учетом оператора определяются методами теории возмущений из уравнения Шре-дингера

как мало отличающиеся от волновых функций щ и энергий еь. Рассматривая резонансный переход электрона из основного в первое возбужденное состояние, необходимо отметить, что под влиянием возмущения У1 основной уровень энергии не расщепляется, и ему соответствуют две волновые функции, различающиеся значением проекции полного момента орбитального электрона (±1/2). Первый возбужденный уровень расщепляется на четыре подуровня, которым соответствуют восемь различных волновых функций.

В системе покоя иона на электрон действует внешнее периодическое во времени электрическое поле кристалла (действием магнитного поля, появляющегося в системе покоя иона, пренебрегаем). Квантовая динамика электрона иона описывается уравнением типа Шредингера:

дТ

т— = ( иь + V + и) т,

д'

где и — оператор, описывающий взаимодействие электрона с внешним периодическим полем, со-

держит только ту гармонику, для которой выполняется условие резонанса (1). Под временем подразумевается собственное время в системе покоя иона. Уравнение (3) решается методами теории возмущений, взаимодействие и предполагается малым по сравнению с гамильтонианом Нь+У{. Волновая функция Т ищется как суперпозиция волновых функций щ, являющихся решением уравнения (2) с зависящими от времени коэффициентами:

Т = ^ aj (')щj exp(-iEj' / Й).

]=1

Подстановка волновой функции (4) в уравнение (3) дает систему уравнений на коэффициенты яД/):

т %=Е(щ1 \и

д'

т ^ д'

j=3 10

j=3

да_

т^. = щ|и|щт) + >2 |и|щт)' а2г'Дт'; т = 3-10.

Здесь < щ| и|щ> - матричный элемент волновых функций щг и щ с оператором периодического взаимодействия и, вычисленный в системе покоя иона в резонансном приближении; А-=(Е-Е1)/Н-уу(кусо5ф+к£тд)); ку и к1 - вектора обратной решетки направлений ОУ и OZ в плоскости. Индексы /=1,2 соответствуют волновым функциям основного состояния, /=3-10 - волновым функциям первого возбужденного состояния.

Решение (5) ищется в виде: а1(') = Л1в<а',

а2(') = Л2в'°',

I (О+Дт)'.

ат(') = Атв—т>-; т = 3-10:

где О - пока неизвестная величина.

Подставляя (6) в (5), получим систему уравне

ний:

10

ЙОА +Е Aj (щ |и|щ) = 0,

j=3 10

ЙОА, +ХAj Щ \иЩ) = 0,

j=3

Й(О + Дт) Ат + А Щ рЩ )' + +Л2 Щ2 |и|* = 0; т = 3 -10.

Система (7) дает десять различных значений О, которые будем обозначать О-, „/'=1-10. В решении (6) необходимо учесть все значения О , каждому из которых соответствует свой набор величин А, (первый индекс нумерует волновую функцию, к которой относится данный коэффициент, второй - частоту О ). Поэтому решение (6) необходимо переписать в виде:

а=Х А /

1=1 10

а2 = Х А /

1=1

ак =Х А

1=1

г<(а 1 +Ак),. к = з-ю.

Амплитуды А^ определяются из системы (7) с соответствующим Ц и начальным условием ах(0)=1, а;(0)=0, /=2-10, которое означает, что первоначально электрон иона находился в состоянии с определенным значением спина. Коэффициенты а/ нормируются условием:

XI а

= 1.

Результат (8), полученный при решении системы (5), не учитывает возможную ионизацию иона при взаимодействии с электронами кристалла. Вероятность ионизации определяется как [8]:

10 2 Рт = 1 - ехр(-£ р( х)с, (О |2)),

где р(х) - плотность электронов кристалла на траектории (траектория в данной модели предполагается прямолинейной, такая аппроксимация применима для больших энергий иона и тонких кристаллов); х - расстояние иона от центра канала; й -толщина кристалла; ст; - сечение ионизации иона из /-го состояния при взаимодействии с электроном кристалла (определяется по формуле Лотца); (|а;(/)|2) - среднее значение величины |а,(/)|2на данной траектории. Вероятность перехода иона в первое возбужденное состояние определяется как:

рех=(1 - р,оп) X а (4 )|:

4 - время, проведенное ионом в кристалле. Вероятность сохранения зарядового состояния иона:

Р = 1 - Р .

тп юп V ^ ^ /

Полученные значения вероятностей (9-11) необходимо усреднить по всем траекториям ионов в кристалле.

Рис. 1, 2 представляют результаты расчетов ЯСЕ при плоскостном каналировании 390 МэВ/нуклон ионов 17+Аг вдоль (220) плоскостей кристалла кремния при толщине 21 мкм. Исследовано ЯСЕ для первой гармоники периодического потенциала кристалла (к=1=1 в (1)).

Полученные резонансные кривые характеризуются двумя близко расположенными пиками, обусловленными переходами электрона из основного состояния в состояния, соответствующие различным компонентам тонкой структуры первого возбужденного состояния, различающихся значением полного момента электрона (1/2 или 3/2). Совпадение в расположении пиков, рис. 1 (вероятность

возбуждения), и провалов, рис. 2 (вероятность сохранения зарядового состояния), обусловлено тем, что электрон в возбужденном состоянии обладает меньшей энергией связи, и, следовательно, легче, чем электрон в основном состоянии, покидает ион при взаимодействии с электронами мишени.

Р А

ех 0,06

Р

Рис. 1.

Вероятность: a) возбуждения Peл; Ь) сохранения зарядового состояния иона Psurv при пролете его через кристалл Si в зависимости от угла (в градусах) между продольной скоростью иона и направлением оси OY в (220) плоскости

Р

вигу 0,40

0,35

0,30

0,25

0,20

0,15

0,10

0,05

3,310

3,315

3,320

3,325 ке

V

Рис. 2.

Вероятность сохранения зарядового состояния иона Ps[„ при пролете через кристалл кремния в зависимости от энергии возмущающей гармоники

2

I=1

1=3

Результаты численного расчета демонстрируют хорошее согласие с экспериментами [6, 9], проведенными при тех же условиях, хотя имеется некоторое смещение теоретической кривой, рис. 2, относительно экспериментальных данных [6]. Расчетные резонансные кривые (рис. 1, 2) имеют четко выраженную дублетную структуру, обусловленную, как указано выше,

тонкой структурой первого возбужденного состояния. Наличие третьего пика на экспериментальных кривых [6, 9] требует дополнительного изучения. Возможно, его появление связано с влиянием кильватерного потенциала, возникающего в результате действия электрического поля иона на электроны кристалла. Этот вопрос будет исследован в дальнейшем.

СПИСОК ЛИТЕРАТУРЫ

1. Okorokov V.V., Tolchenkov D.L., Tolchenkov I.S. et al. The Coherent Excitation of Atoms Moving Through a Crystal Experimental Results // Physical Letters A. - 1973. - V. 43. - № 6. - P. 485-490.

2. Shindo S., Ohtsuki Y.H. Theory of the Okorokov Effect // Physical Review B. - 1976. - V. 14. - № 9. - P. 3929-3932.

3. Оцуки Ё.-Х. Взаимодействие заряженных частиц с твердыми телами. - М.: Мир, 1985. - 280 c.

4. Okorokov V.V. Comments on the Coherent Excitation: Some Problems // 8-th Japan-Russia Intern. Symp. on Interaction of Fast Charged Particles with Solids, Kyoto, Japan, 24-30 November 2002. - Kyoto: Quantum Science and Engineering Center, Kyoto University, 2003. - P. 207-211.

5. Кривошеев О.Е., Пивоваров Ю.Л. Резонансное когерентное возбуждение быстрых ионов Ne+9 в кристалле - компьютерный эксперимент // Письма в ЖЭТФ. - 1992. - Т. 56. - Вып. 5. - С. 246-250.

6. Azuma T., Ito T., Takabayashi Y. et al. Resonant Coherent Excitation of Hydrogen-Like Ar Ions to the n=3 States // Physica Scripta. -2001. - V. 92. - P. 61-65.

7. Pivovarov Yu.L. Coherent Excitation ofHydrogen-Like Relativistic Heavy Ions in a Crystal: Structure of Electronic Levels and Resonance Width // Nuclear Instruments Methods B. - 1998. - № 145. - P. 96-101.

8. Othsuki Y. H., Yamashita Y. Peak Profiles of the Okorokov Effects for Heavy Ions in a Crystal // Physical Review B. - 1980. - V. 22. -№3. - P. 1183-1189.

9. Azuma T., Ito T., Komaki K. et al. Impact Parameter Dependent Resonant Coherent Excitation of Relativistic Heavy Ions Planar Channeled in Crystals // Physical Review Letters. - 1999. - V. 83. - № 3. - P. 528-531.

Поступила 15.12.2006 г.

УДК 535.218

ЭФФЕКТИВНОСТЬ ПЕРЕДАЧИ ЭНЕРГИИ АДСОРБЦИИ И РЕКОМБИНАЦИИ АТОМОВ ТВЕРДОМУ ТЕЛУ ПРИ РАЗЛИЧНЫХ МЕХАНИЗМАХ ВОЗБУЖДЕНИЯ

Ю.И. Тюрин, В.Д. Хоружий, С.Х. Шигалугов*, Ю.А. Сивов, Т.В. Смекалина

Томский политехнический университет *Норильский индустриальный институт E-mail: tyurin@fnsm.tpu.edu.ru

Рассмотрены элементарные механизмы возбуждения электронной (ионной) подсистемы твердых тел атомными частицами тепловой энергии: ионизационный механизм генерации электронно-дырочных пар, адиабатический и неадиабатический переходы, многоквантовый колебательно-электронный переход. Установлены особенности каждого из механизмов в зависимости от условий возбуждения. Получены аналитические закономерности, позволяющие сравнить эффективность генерации электронно-возбужденных состояний при различных механизмах возбуждения.

Введение

Наряду с одноквантовыми механизмами коле-бательно-фононной, колебательно-электронной и колебательно-колебательной релаксации возбужденных вновь образованных локализованных состояний на поверхности твердого тела, имеют место конкурирующие каналы высокоэнергетического захвата твердым телом в одном элементарном акте значительной части энергии, выделившейся в актах адсорбции, рекомбинации и иных химических превращений на поверхности.

Высокоэнергетическая аккомодация особенно эффективна при передаче энергии на электронные степени свободы кристалла. Генерация электронно-возбужденных состояний в твердых телах атомами и

молекулами наглядно проявляется в явлениях гетерогенной хемилюминесценции (ГХЛ) и хемоэмис-сии, неравновесной хемопроводимости, вентильной хемоЭДС, хемомагнитоэлектрическом эффекте [1].

Рассмотрим элементарные механизмы возбуждения электронной подсистемы твердых тел атомными частицами тепловой энергии.

1. Ионизационный механизм генерации

электронно-дырочных пар

Модель ионизационного механизма генерации неравновесных электронно-дырочных пар неоднократно использовалась для объяснения наблюдавшихся экспериментально особенностей ГХЛ и родственных эффектов [1, 2].