Безрадиационные переходы в малочастичных системах с кулоновским взаимодействием тема автореферата и диссертации по физике, 01.04.16 ВАК РФ
Монахов, Дмитрий Евгеньевич
АВТОР
|
||||
кандидата физико-математических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Дубна
МЕСТО ЗАЩИТЫ
|
||||
1999
ГОД ЗАЩИТЫ
|
|
01.04.16
КОД ВАК РФ
|
||
|
■ // хпз
у - ^ .
ОБЪЕДИНЕННЫЙ ИНСТИТУТ ЯДЕРНЫХ ИССЛЕДОВАНИЙ Лаборатория теоретической физики им Н. Н. Боголюбова
На правах рукописи УДК 539.189
МОНАХОВ ДМИТРИЙ ЕВГЕНИЕВИЧ
БЕЗРАДИАЦИОННЫЕ ПЕРЕХОДЫ В МАЛОЧАСТИЧНЫХ СИСТЕМАХ С КУЛОНОВСКИМ ВЗАИМОДЕЙСТВИЕМ
(01.04.16 - физика ядра и элементарных частиц)
ДИССЕРТАЦИЯ на соискание ученой степени кандидата физико-математических наук
Научный руководитель:
доктор физико-математических наук
В.Б. Беляев
ДУБНА, 1999 г.
СОДЕРЖАНИЕ
Введение 3
Глава 1 Тройные столкновения е- р 7Ве в солнечной
плазме 15
1.1 Реакции двойного и тройного соударения 15
1.2 Скорость реакции 17
1.3 Адиабатическое приближение 20
1.4 Ядерные волновые функции 23
1.4.А Состояние рассеяния........................................23
1.4.В Связанное состояние.......... ^ .............24
1.5 Результаты и их обсуждение 25 Глава 2 Оже распад антипротонного атома гелия 30
2.1 Основной формализм 31
2.2 Вариационный метод 32
2.3 Анализ матричного элемента перехода 33
2.4 Метод связанных каналов 34
2.5 Волновая функция конечного состояния 36
2.6 Результаты численного расчета 37
2.6.А Численная процедура........................................37
2.6.В Короткоживущие состояния................................43
2.6.С Скорости перехода..........................................45
2.7 Результаты и выводы 50
Глава 3 Основные свойства трех дейтронной мюонной молекулы 56
3.1 Энергия и размер Бз// 57
3.2 Ядерный переход 65
3.3 Скорость Оже распада 68
3.4 Выводы 71
А Вывод формулы для скорости тройной реакции 73
В Двухуровневая модель 76
С Вычисление матричных элементов Оже перехода в
Нере 78
О Вычисление Оже перехода в Бзд 82
Введение
В различных областях современной физики возникает большое количество задач, требующих для своего решения расчета квантово-механических систем из нескольких (п > 3) частиц. Для их описания статистические методы из-за малого количества тел неприемлемы, в то же время применение стандартных методов квантовой теории может привести к неинтегрируемым в общем виде уравнениям, которые даже могут не иметь единственного решения. Такой класс задач стимулировал появление квантовой теории для трех и более частиц, бурное развитие которой началось с созданием Л.Д. Фадде-евым математически корректной квантовой теории трехчастичного рассеяния [1]. В дальнейшем этот подход был обобщен на задачу четырех тел, а так же на случай кулоновских потенциалов [2].
Успешное практическое применение уравнений Фаддеева связано в основном с нуклонными системами, где взаимодействий частиц описывается короткодействующими потенциалами. Включение в уравнения Фаддеева кулоновского взаимодействия влечет за собой появление значительных вычислительных трудностей. Например, наибольший успех в расчете систем с кулоновским взаимодействием в рамках фад-деевского подхода, как для задачи на связанные состояния [3], так и для задачи рассеяния [4], был достигнут в работах, где авторам пришлось ограничиться случаем двух равных масс и полным моментом системы равным нулю. И это несмотря на то, что вычисления производились на одном из самых мощных компьютеров. Такая
ситуация связана с тем, что приходиться решать систему трехмерных уравнений. Кроме того в задаче рассеяния, учитывая дальнодей-ствующий характер кулоновского потенциала, приходится проводить интегрирование до достаточно больших расстояний для достижения правильного асимптотического поведения. Следует отметить, что как при увеличении количества частиц, находящихся в рассмотрении, так и при увеличении полного момента системы количество уравнений возрастает, что влечет за собой новые вычислительные трудности.
Однако для решения некоторых конкретных физических задач можно указать физически адекватные приближения, использование которых позволяет свести машинный счет к минимуму или по крайне мере сильно сократить. Не претендуя на указание универсального рецепта, в настоящей работе будут продемонстрированы задачи, содержащие кулоновское взаимодействие, где удается получить результат, не прибегая к решению точных уравнений типа Фаддеева. Одной из таких задач является безрадиационный синтез 8В в солнечной плазме. Исследование этого процесса крайне актуально в связи с проблемой солнечных нейтрино [5], которая выражается в противоречии экспериментальных данных стандартной модели солнца.
Так называемая стандартная модель солнца разрабатывалась в течение десятилетий коллективными усилиями многих авторов. К настоящему времени это вполне законченная теория, которая объясняет все, за небольшим исключением, наблюдаемые характеристики солнца. Одним из таких исключений является предсказанный ею поток солнечных нейтрино, который в несколько раз превосходит наблюдаемый экспериментально в земных условиях (детальное обсуждение этой проблемы можно найти в [5]). Обсуждается несколько возможных причин такого несоответствия, которые связаны либо с ядерными реакциями, производящими нейтрино, либо с их прохождением к земле.
Выход солнечных нейтрино происходит в ядерных реакциях, протекающих в плазме, состоящей из легких ядер (в основном протонов) и электронов. Последовательность этих реакций начинается с протон-протонных соударений и поэтому называется рр-циклом. В стандартной модели предполагается, что этот цикл состоит из следующих реакций:
р + р 2Н + е+ + V р + р + е- 2Н + ^ р + 2Н 3Не + 7 3Не + 3Не 4Не + р + р 3Не + 4Не 7Ве + 7, затем 7Ве может "выгореть" по одному из двух каналов:
е- + 7Ве 71л + V р + 7Ве 8В + у
р + 71л 8Ве + 7 8В ->• 8Ве* + е+ +
8Ве 4Не + 4Не 8Ве* 4Не + 4Не.
у
Распад 7Ве в такой последовательности представляет особый интерес, так как анализ экспериментальных данных приводит к парадоксальному результату: выход нейтрино при захвате электронов ядрами 7Ве должен быть сильно подавлен или даже отрицателен [6, 7].
Для разрешения упомянутого выше парадокса необходимо проанализировать все возможные каналы уничтожения ядер 7Ве, в частности сечение образования 8В в результате захвата протонов ядром 7Ве при низких 20 кэВ) энергиях, соответствующих гамов-скому окну. Экспериментальные данные из-за большого кулоновского барьера между ядрами имеются только для более высоких энергий
120 кэВ) [8], и в дальнейшем их приходится экстраполировать в область низких энергий, что может привести к значительной ошибке. Теоретическое исследование энергетической зависимости сечения про-
цесса 7Ве(р,7)8В было впервые проведено в работах [9, 10], где был произведен расчет Е1 перехода только с б-волновым входным каналом. Позднее были проведены более сложные микроскопические расчеты, в частности, основанный на методе генератора координаты трех-кластерный расчет [11] и работы [12, 13, 16], где решалась восьми нуклонная задача в рамках метода резонирующих групп.
В последние годы большое внимание было посвящено точному учету влияния эффектов электронного экранирования на протекание ядерных реакций в звездной плазме (смотри [14, 15, 16, 17]). Обычно такие эффекты рассматриваются в приближении электростатического экранирования. Это приближение как в классической, так и в квантовой реализации адекватно описывает плазму с большими скоростями сталкивающихся ядер, когда электронная плотность в процессе соударения остается постоянной. Было показано [17, 18], что общепринятая теория Дебая-Хьюкеля не может быть применена для вычисления экранировок в не очень плотных звездах, таких, например, как Солнце. Это стимулировало появление ряда работ [19, 20, 21, 22, 23], где проводится квантовомеханический расчет эффекта электронного экранирования и его влияния на скорость ядерных реакций в условиях солнечной плазмы. Например, в [20] показано, что наличие связанного на ядре 7Ве электрона увеличивает скорость реакции 7Ве(р, т)8В на 10%.
В главе 1 впервые рассмотрена реакция безрадиационного захвата протона ядром 7Ве, когда в начальном состоянии находится три частицы: электрон, протон и ядро 7Ве. В таком подходе, в отличие от стандартной модели, рассматриваются реакции, где начальное состояние может представлять собой не только две (бинарные реакции) , но и три частицы (трехчастичные реакции) в непрерывном спектре. Между этими двумя подходами существуют по крайне мере два су-
щественных отличия, которые условно можно классифицировать как кинематические и динамические. Первое представляет собой разницу в правилах отбора в двойных и тройных реакциях, второе вытекает из взаимозависимости различных ядерных процессов. Так некоторые бинарные ядерные реакции, запрещенные законами сохранения (спин, изоспин, четность ...), возможно смогут протекать при наличии третьей частицы, также участвующей в законах сохранения. Таким образом, трехчастичный механизм может снимать кинематические ограничения на отдельные бинарные реакции и играть существенную роль в ядерных процессах внутри звезд, где плотность вещества достаточно велика. Динамика трехчастичного движения тоже может привести к иной физический картине: процессы, рассматриваемые как независимые в двухчастичном подходе, могут стать зависимыми при включении трехчастичной динамики. Например, процессы е~ + 7Ве и р + 7Ве в этом случае станут зависеть друг от друга, поскольку они могут начинаться из одного начального состояния е~ + р + 7Ве.
В главе 1 рассмотрен ряд реакций с ядром 7Ве, которые не входят в стандартный рр-цикл. Проведен анализ влияния таких реакций на эволюцию ядра 7Ве в солнечной плазме и обсуждена их роль в интерпретации экспериментальных данных по солнечным нейтрино. В качестве примера в адиабатическом приближении рассмотрен процесс безрадиационного производства ядер 8В в тройной реакции р + е~ + 7Ве -4 8В+е . В центре Солнца скорость такой реакции составляет 10~4 от аналогичного бинарного процесса р + 7Ве —У 8В + 7.
В разделе 1.1 качественно исследована роль реакций тройного типа, в разделе 1.2 дано описание использованного формализма. Разделы 1.3 и 1.4 содержат вывод вспомогательных величин для вычисления скорости реакции. Результаты и их обсуждение отнесены к разделу 1.5. Некоторые детали относительно вывода скорости реакции даны в приложении А.
В главе 2 проведено вычисление Оже распада антипротонного атома гелия Нере, который представляет собой обычный атом гелия, где один из электронов заменен на антипротон. В последние годы к его изучению обращено большое внимание как теоретиков, так и экспериментаторов, поскольку он является первой и единственной малотельной системой атомного масштаба, содержащей антипротон [24, 25]. Существование системы такого типа (адронные атомы НеК~е, Не7г_е) было качественно предсказано в шестидесятых [26], а первые расчеты этих систем были сделаны в работе [27]. Антипротонный гелий, обладающий черезвычайно большим угловым моментом, живет микросекунды, что гораздо больше обычного времени жизни антипротона в среде (около Ю-"12 с).
Продолжительное время жизни антипротонного гелия объясняется следующими причинами. Большое значение полного углового момента ~ 30 сосредоточено в основном в угловом моменте между тяжелыми частицами, что является причиной центробежного барьера, который подавляет процессы аннигиляции. Экранирование антипротона электроном подавляет в свою очередь процессы девозбуждения за счет столкновений с окружающими атомами обычного гелия. Таким образом, основным каналом распада становится процесс испускания электрона, т.е. Оже переход.
Сравнивал спектры антипротонного гелия Нере и иона Hep, можно заметить, что для наиболее низких долгоживущих уровней Оже переход будет сильно подавлен в силу большой мулыпиполярности перехода До, т.е. углового момента вылетевшего электрона. Очевидно, что электрон должен унести по крайне мере разницу между полным угловым моментом L и наибольшем моментом энергетически разрешенного состояния иона остатка Hep.
Скорость Оже перехода существенным образом зависит от муль-типолярности перехода До, при увеличении До на единицу скорость
перехода уменьшается на три порядка. Например, скорость Оже перехода уменьшается, начиная со значения 108 с-1 при Ао = 3 до 105 с-1 для Ао = 4. Эта оценка была получена еще в работе [27] и подтверждена расчетами [28, 29]. Так как радиационное время жизни антипротонного гелия составляет примерно 10~6 с [29, 30, 31], то его можно считать метастабильным в состояниях с Ао > 4. Таким образом, метастабильный антипротонный гелий испытывает серию радиационных дипольных переходов до тех пор, пока не достигнет уровня с Ао — 3, где возможен сравнительно быстрый Оже распад, после которого немедленно следует аннигиляция из-за очень короткого (порядка Ю-12 с) времени жизни иона Hep.
Описанная выше схема эволюции антипротонного гелия является общепринятой и подтверждена на эксперименте [32], в котором используется метод индуцированной лазером резонансной спектроскопии. В этом методе лазерный импульс с резонансной длинной волны индуцирует переход между уровнями со значительной разницей во времени жизни, в результате чего наблюдается резкий аннигиляцион-ный пик. Благодаря этому методу, получена детальная информация об уровнях, их временах жизни и заселенности.
Ввиду вышеизложенных аргументов ясно, что для правильного количественного описания процессов в антипротонном гелии необходим точный расчет скоростей Оже распада. Скорости радиационного распада вычислены достаточно точно в работах [29, 30, 31], они хорошо согласуются друг с другом и с последними экспериментальными данными [33], в то время как вычисление скоростей Оже перехода представляет собой более сложную задачу. Наиболее интересны распады из состояний с мультиполярностью перехода Ао = 3,4. Это объясняется тем, что каскад радиационных переходов в антипротонном гелии прерывается на состоянии с Ао = 3, а в состоянии с Ао = 4 Оже и радиационные переходы являются конкурирующими, особенно
в случае 3Нере.
В главе 2 будет показано, что наибольший вклад в матричный элемент Оже перехода дают очень малые компоненты волновой функции. Именно поэтому в первых расчетах, где использовались сравнительно простые волновые функции, были получены только оценки для скорости Оже распада. Например, в самом первом расчете [27] вариационная волновая функция представляла собой произведение антипротонной волновой функции и я-волновой функции электрона. Уже в работе [28] была использована более сложная волновая функция, где были учтены слагаемые с разными угловыми моментами электрона и тяжелых частиц.
Исследование вкладов в матричный элемент перехода от малых компонент вариационной волновой функции было проведено в работах [29, 34, 35], где также были получены скорости Оже распада для некоторых состояний с мультиполярностью Ао = 3,4. Тот же самый метод совместно с теорией возмущений был использован в [35] при рассмотрении изтопической зависимости скоростей распада. Более совершенная вариационная волновая функция была использована в [36] для вычисления сходящихся скоростей распада ряда метастабильных состояний. Два типа вычислений представлены в [37]: в первом из них используется волновая функция из [36], а полученный результат сравнивается с расчетом проделанном по методу Борна-Оппенгеймера. Сходимость скорости Оже распада, расчитанной с волновой функцией из [36], сомнительна [37], так как подынтегральное выражение в матричном элементе перехода осциллирует. Результат полученный исходя из волновой функции вычисленно в приближении Борна-Оппенгеймера признан неудовлетворительным [37].
Следует отметить, что малые компоненты волновой функции, которые так необходимы для точного вычисления скорости Оже
распада, не могут быть найдены прямым вариационным расчетом с достаточной точностью. При вычислении матричного элемента перехода, главный вклад в интеграл дают большие р, где волновая функция экспоненциально мала. К тому же вклады от разных компонент взаимно компенсируют друг друга, так как имеют противоположный знак [29, 34].
В главе 2 были вычислены скорости Оже распада для долгоживу-щих состояний антипротонного атома гелия 3'4Нере. Для определения волновой функции был использован метод, являющийся комбинацией вариационного подхода и метода связанных каналов, что позволило получить более точное описание малых компонент волновой функции, играющих решающую роль в определении скорости Оже перехода. Показано, что некоторые метастабильные состояния обладают энергией близкой к энергии, так называемых короткоживущих состояний, обладающих особенным строением волновой функции. Вычисление скоростей распада для таких метастабильных состояний требует учета эффекта смешивания волновых функций.
Основной формализм дан в разделе 2.1, вариационный метод описан в 2.2, раздел 2.3 содержит анализ вкладов в матричный элемент перехода от различных компонент волновой функции, метод связанных каналов и волновая функция конечного состояния рассмотрены в разделах 2.4 и 2.5 соответственно