Главная » Просмотр файлов » Диссертация

Диссертация (1145383), страница 28

Файл №1145383 Диссертация (Единый аналитический и вычислительный подход к решению квантовой задачи трёх тел) 28 страницаДиссертация (1145383) страница 282019-06-29СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 28)

работы [102, 271] и приведённые в них ссылки) позволяют провеститакие тесты с хорошей точностью и достоверностью.Для удобства обсуждения, выпишем здесь гамильтониан именно для дан­ной системы. В атомных единицах он имеет вид:111 = − Δr1 − Δr2 −−+.22|r1 | |r2 | |r1 − r2 |(5.3)Здесь r1 , r2 – радиус-векторы электронов, – заряд ядра, масса ядра предпола­гается бесконечной. Для бесконечной массы ядра, координаты r1 , r2 совпадаютс обычными координатами Якоби x̃, ỹ. Далее будут использоваться обозначе­ния r1 , r2 , как общепринятые.

Уравнение для этого гамильтониана являетсяуравнением Шредингера (1.7), так что оно может быть исследовано развиты­ми в данной работе методами. Так как уравнения рассеяния для различныхполных угловых моментов независимы, они могут быть решены также незави­симо. В данной работе будет рассмотрено только рассеяние с нулевым полныммоментом.Проекция 0 гамильтониана (5.3) на подпространство нулевого полногоуглового момента определяется соотношением (1.29) в случае = 0:(5.4)0 = + (1 , 2 , ).Здесь кинетическая энергия задаётся как1 2=−−2 12(︂11+212 222)︂ (︂2+cot2)︂1 2−,2 22(5.5)где = |r |, = 1, 2, а – угол между векторами r1 и r2 . Для представ­ления кинетической энергии в таком виде сделана стандартная замена вол­новой функции Ψ(1 , 2 ) на Ψ(1 , 2 )/(1 2 ).

Такой выбор совпадает со стан­дартным в данной предметной области. Для волной функции, переопределён­ной таким образом, необходимо добавить дополнительное граничное условие183Ψ(1 , 0) = Ψ(0, 2 ) = 0, гарантирующее конечность исходной волновой функциипри нулевых значениях координат. Полный потенциал (1 , 2 , ) записываетсяв виде (1 , 2 , ) = −− + 12 (1 , 2 , ),1 2(5.6)где электрон-электронное взаимодействие12 (1 , 2 , ) =11= √︀ 2.|r1 − r2 |1 + 22 − 21 2 cos Поскольку электроны являются тождественными фермионами, волноваяфункция должна обладать определённой симметрией при их перестановке.

Сим­метризованная волновая функция Ψ определяется как Ψ = Ψ, где опера­тор симметризации задаётся стандартным образом1 = √ (1 + (−1) 12 ).2(5.7)Здесь значения =0 и =1 отвечают синглетному и триплетному состояниям,а оператор 12 переставляет электроны 1 и 2. Так как оператор симметризации12 коммутирует с гамильтонианом 0 , симметризованная волновая функцияудовлетворяет тому же самому уравнению Шредингера (0 − )Ψ = 0. Та­ким образом, все преобразования уравнения можно осуществлять в несиммет­ризованном виде, а затем применить оператор симметризации к полученнымуравнениям.Само по себе уравнение (5.4) не содержит никаких упрощений и приближе­ний: оно является первым из системы несвязанных уравнений для различныхполных угловых моментов.

Результаты его исследования приведены в разде­ле 5.3.3.5.3.2. Модель Темкина-ПоэтаВ данном разделе будет рассмотрена модель Темкина-Поэта [272, 273] (ТП­модель) рассеяния электрона на атоме водорода и одноэлектронном ионе. В184этой модели предполагается, что каждая электрон-ядерная подсистема обла­дает только нулевым моментом. Несмотря на такое упрощение, эта модель со­храняет многие существенные свойства и сложности исходной задачи, позво­ляя исследовать их с меньшими вычислительными затратами.

Для ТП-моде­ли рассеяния электрона на водороде, имеется обширная литература (см. на­пример [97, 274]), и доступны прецизионные референсные расчёты [102, 104].Вычисления процессов рассеяния в системе e-He+ для ТП-модели представле­ны меньше.

Тем не менее, можно назвать результаты, полученные CCC-мето­дом [275–277], методом R-матрицы [275, 276, 278] и с помощью решения времен­ного уравнения Шредингера [279].При проецировании уравнения (5.4) на -волновые состояния электронов,кинетическая энергия записывается в виде1 21 2=−−.2 12 2 22(5.8)Электрон-ядерные взаимодействия не изменяются, а электрон-электронное вза­имодействие приобретает вид12 (1 , 2 ) = 1/ max {1 , 2 } = √︀2(1 + 2 )2 +√︀(1 − 2 )2.(5.9)Второе представление необходимо для применения к потенциалу внешнего ком­плексного вращения. Полученный модельный оператор является двумерным,что объясняет относительную простоту его численного анализа.К уравнению Шредингера с этим модельным гамильтонианом применяет­ся метод расщепления потенциала, описанный в разделе 1.4.3.

Для удобства,приведём здесь уравнения, полученные для ТП-модели с его помощью. Реше­ние записывается в виде суммы четырёх слагаемых, аналогичных представле­нию (1.75), но в котором к каждому слагаемому применён оператор симметри­зации (5.7): Ψ = Ψ0 + Φ0 + Ψ1 + Φ1 .185(5.10)Функции Φ , = 0, 1, удовлетворяют уравнениям[︂(︂)︂]︂ + − − + 12 (1 , 2 ) − Φ (1 , 2 )2(︂ 1)︂= − − + 12 (1 , 2 ) ( − 1 )Ψ (1 , 2 ).1(5.11)Искажённая падающая волна Ψ0 (1 , 2 ) задана явным образом:⎧⎨ ^ ( ),1 < .0 0 1Ψ0 (1 , 2 ) = (2 )⎩ 0 0 ( , 1 ) + + ( , 1 ), 1 ≥ .0 0(5.12)Здесь (2 ) – нормализованная волновая функция связанного состояния мише­ни с энергией , = + 2 /2, параметр Зоммерфельда = −( − 1)/ , аостальные обозначения введены в разделе 1.4.3.Функция Ψ1 определяется как решение неоднородного уравнения(︂(︂)︂)︂ − + − + 12 (1 , 2 ) (1 − ) − Ψ1 (1 , 2 )21(︂)︂11=−Ψ0 (1 , 2 )ℎ(2 − 1 )(1 − ).1 2(5.13)Из-за множителя ℎ(2 − 1 )ℎ(1 − ) правая часть уравнения (5.13) отлична отнуля только в клинообразной области 1 > , 2 > 1 .

Благодаря комбинациифункций (2 )ℎ(2 − 1 ), правая часть быстро убывает по обеим координатам1 , 2 для любого . Эта комбинация множителей приводит также к тому, чтопри увеличении , правая часть убывает как (), т.е. экспоненциально. Врезультате, эффект от поправочные слагаемые Ψ1 , Φ1 также экспоненциальноубывает с ростом , см. рисунок 5.15.Как было указано в разделе 1.4.3, уравнения (5.11, 5.13) могут быть реше­ны с помощью метода комплексных вращений, поскольку они являются част­ным случаем рассмотренных там уравнений.

Наличие симметризации не раз­рушает возможность решения, поскольку функции в правой части уравненийостаются быстро убывающими по обеим координатам.После решения уравнений и вычисления волновой функции, необходимоопределить амплитуды рассеяния. Обозначим амплитуду рассеяния (упру­186гого или возбуждения), отвечающую начальному состоянию мишени () и ко­нечному (). Как было показано в разделе 1.4.5, амплитуды могут бытьвычислены с помощью проецирования волновой функции на двухчастичное свя­занное состояние (2 ).

Амплитуда находится как предел такой проекции набесконечности:Z√(︀)︀−1 ( ) = lim 2 0+ ( , 1 )2 (2 )Φ (1 , 2 ).∞1 →∞(5.14)0Симметризованное слагаемое здесь отсутствует, так как волновая функция свя­занного состояния экспоненциально убывает при увеличении 1 . При примене­нии метода внешнего комплексного вращения, волновая функция находится вобласти 1 , 2 ≤ .

Это означает, что интегрирование в (5.14) должно бытьограничено интервалом 2 ∈ [0, ], и значение правой части вычислено принекотором 1 ≤ . Вычисленные таким образом амплитуды стремятся к точ­ным значениям при → ∞. Сечения рассеяния определяются в терминахамплитуд как( ) = 4 | ( )|2 .(5.15)Сечения рассеяния могут быть также определены и другими способами.В частности, для вычисления сечений может использоваться оптическая теоре­ма [97]Z]︀⃒[︀]︀* [︀8 = − 2 ℑm Φ (, 2 ) Φ (1 , 2 ) ⃒1 = 2 .1(5.16)0Здесь индекс может означать полное сечение, сечения возбуждения или иони­зации, а – оператор проекции на соответствующие состояния [97].Ещё одна возможность вычислить амплитуды рассеяния – использоватьпредставление в терминах поверхностного интеграла [274]:[︂ ]︂Z√Φ(,)Φ(,)1212 ( ) = 2 2 12 (2 )0 ( 1 ) − 0 ( 1 ).11110(5.17)187Правая часть этого выражения должна быть вычислена при некотором 1 ≤ ,например, 1 = .ТП-модель: рассеяние электрона на водородеВ случае рассеяния электрона на водороде, = 1 и = 0.

Тогда 0 = 1,асимптотическое кулоновское взаимодействие отсутствует, и уравнение (5.11)для Φ0 совпадает с уравнениями, использованными в работах [97, 104, 274] дляисследования ТП-модели. Таким образом, полученные данные можно сравнитьс указанными результатами для проверки точности и стабильности разработан­ных методов.Опишем параметры использованной численной схемы.

Область в плоско­сти 1 , 2 разбивалась на прямоугольную сетку, образованную прямым произ­ведением одинаковых сеток по обеим координатам. Для каждой из координат,использовались пять КЭ для малых расстояний [0-4] а.е. и четыре КЭ суммар­ной длины 40 а.е.

для дискретизации области за точкой комплексного вращения. Промежуточная область делилась на элементы длиной 3 а.е. Таким образом,общее количество конечных элементов и размер глобальной МКЭ матрицы за­висят от радиуса комплексного вращения .

Для большинства вычислений,степень полиномов в МКЭ была выбрана равной семи. Максимальный размерматрицы, используемой в вычислениях, был равен 166464, а соответствующаястепень разрежённости – 2 · 10−4 . Значение угла комплексного вращения маловлияло на результаты и было выбрано равным Θ = 45∘ .На рисунке 5.5 показано поведение 1 и 5 синглетных сечений рассеянияв зависимости от радиуса расщепления , = .

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

Характеристики

Список файлов диссертации

Единый аналитический и вычислительный подход к решению квантовой задачи трёх тел
Свежие статьи
Популярно сейчас
Почему делать на заказ в разы дороже, чем купить готовую учебную работу на СтудИзбе? Наши учебные работы продаются каждый год, тогда как большинство заказов выполняются с нуля. Найдите подходящий учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6417
Авторов
на СтудИзбе
307
Средний доход
с одного платного файла
Обучение Подробнее