Комплексные алгоритмы анализа квантовых систем во внешних полях

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

Канонические преобразования гамильтониана к нормальной форме широко применяются для интегрирования и анализа уравнений движения нелинейных систем методами теории возмущений [1−7], включая современные модели, например, управления квантовомеханическими процессами [8, 9], регулярной и хаотической динамики заряженных частиц высоких энергий в прямых и изогнутых кристаллах [10], или квантовых вычислений и алгоритмов [11]. С другой стороны, ридберговские атомы с долгоживущими энергетическими уровнями, имеющими главные квантовые числа n ~ 50 и максимальные орбитальные и магнитные квантовые числа, в настоящее время интенсивно исследуются в микроволновых резонаторах как элементы будущего квантового компьютера [12]. Кроме того в экспериментах с ионами в радиочастотных ловушках изучаются водородоподобные ионы, например 9Ве+, с одним электроном на внешней оболочке и, соответственно, с двумя внутренними состояниями (сверхузкие или метастабильные уровни энергии), которые связываются с колебательными степенями свободы ионов в ловушках [12]. Поскольку гамильтониан водородоподобного атома связан с гамильтонианом гармонического осциллятора посредством преобразований типа Леви-Чивита [13], то метод квантования нормальных форм [14] применяется и для расчетов динамических характеристик атомных систем во внешних полях [15]. Применение различных систем компьютерной алгебры (MATHEMATICA, MAPLE и REDUCE) существенно расширяет возможности различных реализаций процедур нормализации [16 — 18] и квантования полиномиальных гамильтонианов [19 — 21], поэтому актуальной задачей является анализ и сравнение алгоритмов в рамках универсального псевдокода с учетом специфики языков символьных вычислений MATHEMATICA, MAPLE и REDUCE с целью повышения эффективности новых реализаций [22]. В разделе 1.1 мы сравниваем алгоритмы приведения класса полиномиальных гамильтонианов (ангармонического осциллятора) к нормальной форме методом Биркгофа-Густавсона и преобразованиями Ли (метод Депри-Хори), а также процедуры нахождения приближенных и точных интегралов движения и построения квантового аналога нормальной формы, представленные в виде универсального псевдокода. Рассмотрены примеры, позволяющие оценить эффективность работы вышеупомянутых алгоритмов и программ, реализованных на языках REDUCE, MAPLE и MATHEMATICA. Представлен расчет квазиклассического спектра двумерного атома водорода в электрическом поле удаленного точечного заряда и дано сравнение со спектром соответствующей квантовой задачи, полученным алгоритмом теории возмущений раздела 1.2.

Различные схемы теории возмущений для водородоподобного атома в поле удаленного точечного заряда и ее приложения в атомной и молекулярной физике изучается более 60 лет [23 — 32], поэтому такую интегрируемую систему естественно рассматривать как эталонную модель для проверки эффективности разрабатываемых символьных алгоритмов [22]. Так полиномиальные решения для невозмущенного водородоподобного атома в сфероидальной системе координат исследовали в работе [27], где было отмечено, что традиционные схемы теории возмущений [33] не приводят к правильным замкнутым аналитическим выражениям для решений, характеризуемых наборами сферических или параболических чисел, что якобы не позволило в работе [26] получить правильное выражение для спектра в пятом порядке теории возмущений. Последовательный, но более сложный в реализации и в дальнейшем использовании полученных решений алгоритм, обеспечивающий факто-ризованное асимптотическое представление решений по большому параметру, классифицированных сфероидальными квантовыми числами был предложен в исходя из возможности разделения переменных [28]. Там же было получено разложение спектра до шестого порядка по обратной степени расстояния до удаленного заряда. Однако достаточно эффективного алгоритма для вычисления волновой функции в виде разложений по малому параметру теории возмущений [34, 35] с коэффициентами, заданными в виде полиномов от параболических или сферических чисел в аналитическом виде, так и не было построено [32]. В настоящее время эта задача остается актуальной, поскольку наличие простых аналитических выражений для спектра и собственных функций водородоподобного атома в поле удаленного заряда необходимы для вычисления асимптотик эффективных потенциалов экзотических двухэлек-тронных атомных систем в адиабатическом представлении [32, 36], а также асимптотик сильно делокализованные состояний антипротонного атома гелия [37] и подобных им ридберговских состояний водородоподобных атомов, рассматриваемых как интегрируемые модели в специальных конфигурациях электрического и магнитных полей [38, 39].

Как известно, группа 0(4) есть группа симметрии трехмерного атома водорода, а группа вращений 0(4,2) в шестимерном псевдоевклидовом пространстве является группой динамической симметрии [40, 41]. В самом деле, одно бесконечномерное неприводимое представление группы вращений связано масштабным преобразованием с волновыми функциями дискретного спектра атома водорода. Это обстоятельство позволяет формулировать алгебраическую схему теории возмущений в рамках алгебры группы 0(4,2) для атома водорода в трехмерном пространстве, заменяя полиномиальное возмущение подходящей комбинацией операторов1. Для построения такой схемы теории.

1 Группа 0(3) есть группа симметрии двумерного атома водорода, группой динамической симметрии которого является группа вращений 0(3,2) в пятимерном псевдоевклидовом пространстве. возмущений принято использовать представление алгебры 0(2,1), исходя из возможности разделения переменных для эффекта Штарка [42] или задачи двух кулоновских центров [28]. Однако, в случае полиномиального возмущения общего типа, например для атома водорода в электрическом и магнитном полях или в неоднородном электрическом поле, когда переменные не разделяются, такое представление неудобно и требуется разработка общей схемы теории возмущений [31]. Для решения спектральных задач удобнее выбрать нормировку собственного вектора таким образом, чтобы она не содержала дробных степеней квантовых чисел, участвующих в процессе алгебраических вычислений в символьном виде. В разделе 1.2 дана формулировка наиболее эффективной схемы теории возмущений, которая не использует идею явного разделения переменных и допускает дальнейшее обобщение при рассмотрении вырожденной теории возмущений. Для того, чтобы изложение материала было максимально наглядным, мы демонстрируем работу предлагаемой схемы на примере типа эффекта Штарка, который позволяет получить правильные волновые функции нулевого приближения в замкнутом виде, а также в некоторых специальных конфигурациях электрического и магнитного полей.

В течение последних десятилетий большое внимание уделяется исследованию различных квантово-оптических моделей с нелинейными гамильтонианами, выраженных в терминах генераторов алгебры Ли, позволяющих изучать новые физические эффекты и явления. Для анализа таких моделей используются в основном численные расчеты [43 — 45], потому что стандартные Ли-алгебраические методы, хорошо приспособленные к решению проблем с линейными (по генераторам алгебры Ли) гамильтонианами неэффективными [46 — 48], а большинство других аналитических методов (алгебраический анзац Бете [49] и др.) требуют сложных вычислений и не всегда дают простые аналитические выражения для физических величин. С другой стороны, стандартные численные схемы ограничены компьютерными ресурсами и не приспособлены для исследования особенностей динамики модели [50 — 52].

Новый универсальный Ли алгебраический подход, улучшающий и аналитические и численные решения физических проблем, был предложен в [53] и разработан в [50, 51, 54, 55] для некоторых нелинейных квантовых моделей, с инвариантными гамильтонианами Н относительно группы Gi: [G^ Н] = 0. Этот подход основан на переформулировке изучаемых моделей в терминах (введенный в [51, 53, 55]) полиномиальных алгебр Ли (PLA) gpd как алгебр динамической симметрии gD: gD — полностью описывающих динамику модели. В разделе 1.3 представлены соответствующие алгоритмы, реализующие данный подход для класса интегрируемых полиномиальных моделей квантовой оптики в подпространстве конечной размерности, а также проведено сравнение точных квантовых расчетов с различными квазиклассическими приближениями.

Ионизация атомной мишени электронным ударом с одновременным измерением на совпадение продуктов реакции широко и интенсивно используется для исследования электронной структуры и механизмов ионизации атомов [56]. Режим большой передачи импульса от налетающего электрона одному из электронов атомной мишени, когда измеряется на совпадение быстрая электронная пара, особенно выгоден для изучения структуры волновой функции атома. В этом случае ионизационный механизм хорошо описывается моделью квазиупругого удара, и соответствующее дифференциальное сечение пропорционально спектральной функции ионизационной дырки в пространстве импульсов [57]. Однако такая простая картина процесса справедлива только при достаточно высоких энергиях налетающего электрона (3−5 кэВ), когда дифференциальные сечения малы и, следовательно, трудно измеримы в эксперименте. Поэтому довольно часто эксперименты проводятся при меньших энергиях (0.5−1 кэВ), что с точки зрения теории требует учета поправок к модели квазиупругого удара. На примере реакций однократной ионизации, или так называемых (е, 2е) — реакций, эти поправки можно условно разбить на две категории: эйкональные поправки за счет эффекта среднего поля в атоме [58] и поправки за счет кулоновского взаимодействия электронов в конечном состоянии [59].

Экспериментально (е, 2е)~ реакции на различных мишенях в кинематике квазиупругого удара исследуются в течение длительного времени, с конца 60-х годов, после появления первых теоретических работ на эту тему [60]. Позже теоретически было предсказано, что (е, 3е)~ реакции, когда в результате встряски ион-остаток испускает третий медленный электрон, дают эксклюзивную информацию об электронных корреляциях в мишени. И даже если измеряется не угол медленного испущенного электрона, а лишь его энергия (т.н. (е, 3−1е) — реакция), то и в этом случае мы можем изучать вклад высших парциальных волн полной волновой функции атома-мишени в дифференциальное сечение, а также радиальные корреляции электронов [61].

Недавно, основываясь на теоретических предсказаниях [61], осуществлён первый (е, 3−1е) — эксперимент на атоме Не в симметричной компланарной геометрии [62]. Измерения проводились при значениях энергии падающего электрона — 580 эВ, энергий обоих быстрых электронов Е — i?2 = Е = 250 эВ и энергии медленного электрона Е^ = 1 эВ. Отметим, что медленный испущенный электрон не детектировался в эксперименте, так как в случае атома гелия данное обстоятельство не принципиально-энергия этого электрона фиксирована законом сохранения. Четырехкратное дифференциальное сечение (в относительных единицах) d4a/dEidE2dQidD, 2 измерялось как функция угла разлета быстрых конечных электронов 6 — 92 = 9 относительно вектор-импульса начального электрона pq. Было установлено, что экспериментальное сечение имеет некоторую характерную угловую структуру, которая позволяет говорить о наблюдении радиальных корреляций, однако оно оказалось примерно на 10° смещено в сторону больших углов разлета по сравнению с теоретическими расчетами в рамках приближения плоских волн [61], что указывает на необходимость учета поправок.

В разделах 3.1 и 3.2 приводится последовательная схема учета поправок к механизму квазиупругого удара в реакциях двухкратной ионизации атомов электронным ударом. На основе приведенного формализма делаются приближения, которые отвечают случаю эйкональных [58] и квазиклассических.

59] поправок, хорошо известных в теории квазиупругих (е, 2е) — реакций. В рамках сделанных приближений выполняются расчеты дифференциальных сечений (е, 3−1е) — реакций в гелии и проводится анализ и сравнение результатов с экспериментальными данными [62].

Эффективное адиабатическое приближение на основе на канонического преобразования недавно рассматривалось на примере трехчастичных расчетов упругого рассеяния и энергий слабосвязанных состояний [63]. Однако исчерпывающий анализ асимптотических результатов этой процедуры, который исключительно важен в вычислениях резонансных состояний как слабосвязанных, так и нулевых энергий (отсчитанных от парного порога) не был выполнен. Для исследования и демонстрации эффективности такого канонического адиабатического приближения мы выбираем эталонную модель трех частиц на прямой с притягивающими потенциалами нулевого радиуса, так как в этом случае энергия и фазы рассеяния известны в элементарных функциях [64, 101]. Используя асимптотическое поведение результирующего дальнодействующего потенциала, мы демонстрируем в разделе 2.3 сходимость адиабатического разложения при конечном значении радиальной переменной. Мы показываем, что каноническое адиабатическое приближение обеспечивает адиабатическое разделение переменных в задаче упругого рассеяния 2 + 1, так же как и правильную асимптотику двумерной волновой функции. Мы получаем правильные верхние и нижние оценки при вычислении энергии состояния с нулевой энергией и подходящее поведение фазовых сдвигов упругого рассеяния.

В математическом моделировании физических процессов важную роль играют точно решаемые модели, когда можно явно определить спектр, собственные функции, резонансы, параметры рассеяния и т. п. К числу таких моделей относится описание движения частиц в потенциале, сосредоточенном на некотором дискретном множестве точек (5 — потенциал). Большинство работ на эту тему сделано в квантовой механике.

Одноцентровый гамильтониан с точечным взаимодействием в трёхмерном случае впервые был исследован Бете и Пайерлсом в 1935 г. при рассмотрении дейтрона [65]. Ферми в 1936 г. при исследовании движения нейтронов в водородной среде ввёл 5 — потенциал, получивший позже название псевдопотенциала Ферми [66]. Обширный обзор приложений в атомной физике дан в монографии Демкова и Островского [67].

Несмотря на огромный объём физических и математических работ по <5 -потенциалам (см., например, монографию [68]), практически в стороне осталась одна достаточно простая тема — можно ли при решении задачи рассеяния приблизить произвольный ограниченный потенциал конечной решёткой, состоящей из 5 — функций. Поскольку задача рассеяния на узлах решетки решается точно, то такой подход имеет преимущества по сравнению с традиционным разложением неизвестной волновой функции в ряд по некоторому ортонормированному базису. В самом деле, при разложении волновой функции в квантовой задаче рассеяния для системы нескольких частиц возникают известные проблемы [69, 70, 101]. Во-первых, практически очень тяжело оценивать вклад отброшенной части базиса. Во-вторых, асимптотика полученного решения, как правило, отличается от асимптотики точного решения. В-третьих, уточнение любых спектральных характеристик (энергии связанных состояний, ширины и положения резонансов, положения угловых пиков и т. п.) требует значительного увеличения размеров базиса.

Ответ на вопрос: можно ли приблизить кусочно-непрерывный ограниченный потенциал конечной суммой сингулярных S — функций и использовать такое представление в качестве алгоритма решения эллиптического дифференциального уравнения? — далеко не очевиден, поскольку в самом уравнении появляется произведение функций, сингулярных в одной точке (сингулярность самой функции накладывается на сильную сингулярность дельта-функции). В работе [104] на основе численного эксперимента показано, что в одном измерении такое приближение хорошо работает и, при увеличении числа дельта-функций, сходится к точному решению. Попутно обнаружено, что сетка равных площадей, предложенная в [67], менее эффективна, чем простая сетка равных интервалов.

В трехмерном пространстве сингулярность волновой функции значительно сильнее, чем простой излом в одномерном случае. Хотя амплитуда рассеяния и здесь определяется специальными граничными условиями, однако приходится отказаться от требования ограниченности волновой функции [68]. В разделе 3.1 предложен способ преодолеть эту трудность и получена математическая схема, которая позволяет в принципе приближать потенциал 5 — функциями в многомерном пространстве и свести исходную задачу нахождения амплитуды рассеяния к решению неоднородной системы линейных алгебраических уравнений.

Динамика заряженной частицы в скрещенных постоянном магнитном и время-зависящем (переменном) электрическом полях [8] или в специальной метрики искривленного пространства-времени [71] проявляет много интересных эффектов. Например, вычисления классической траектории электрона показывают, что электромагнитное излучение релятивистского электрона в постоянном магнитном поле и сверхкоротком интенсивном лазерном импульсе имеет богатую структуру в направлениях как магнитного, так и электрического полей [72]. Если резонансное условие между циклотронной частотой и лазерной частотой имеет место, то радиационные потери минимизируются и экстремальное ускорение электрона становится возможным на очень короткой длине пробега электрона [73]. Решение уравнения Дирака для атомного электрона в комбинированном сильном магнитном и лазерном поле приводит к кольцевому распределению в резонансном режиме [74], подобно движению волнового пакета свободного электрона в постоянном магнитном и переменном электрическом полях[8]. Наложение статического магнитного поля может усилить скорость ионизации или стабилизировать квантовую систему при воздействии интенсивного лазерного импульса [75]. Тот же эффект наблюдается для модели системы с взаимодействием нулевого радиуса [76]. Более того, не только одноэлектронной, но и многоэлектронной динамикой можно управлять используя конфигурации скрещенных магнитных и электрических полей [77]. Эти наблюдения показывают большие потенциальные возможности контроля динамики классических и квантовых систем, используя оптимизацию различных параметров в комбинации магнитного и электрических полей.

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

Для аппроксимации решения по временной переменной обычно применяют схемы типа Кранка-Николсона [78 — 80]. Для аппроксимации решения по пространственным переменным обычно используют разложение решения по ортогональному базису с выделением радиальной переменной [81], при этом радиальное решение в предположении его достаточно гладкого поведения аппроксимируют методом конечных элементов [82, 83]. Для быстроосцилирую-щего решения можно выделить аналитически быстроменяющийся фазовый множитель, а для нахождения достаточно гладкой огибающей, применить схему Кранка-Николсона и затем вернуться к исходному решению. В разделе 3.2 такая адаптивная схема реализована, а ее эффективность продемонстрирована на примере эталонной точно решаемой модели двумерного осциллятора в электрическом поле [8]. Ослабление возмущений искомого решения на границе можно обеспечить искусственным введением мягкой поглощающей диафрагмы [84] или комплексного скейлинга в асимптотической области [85]. Ранее также был предложен так называемый метод расщепления волновой функции [86], который использует представление ее в виде суммы асимптотической и взаимодействующей частей. При помощи граничных условий излучения [87] можно физически корректно разрешить этот вопрос. В случае, если наибольший интерес представляет именно поведение вытекшего излучения, перспективным является использование адаптивной схемы в расширяющейся системе координат, представленной в разделе 3.2. Здесь показано, что для интенсивного сверхкороткого лазерного импульса, огибающая волнового пакета определяет распределение по импульсам плотности вероятности ионизации электрона.

Современные лазерные эксперименты стимулируют компьютерное моделирование динамики экзотических малочастичных кулоновских систем. В этом случае конструирование алгоритмов расщепления по физическим параметрам позволяет экономить ресурс ЭВМ [88]. Имеется два важных требования к разрабатываемым схемам и алгоритмам. Они должны быть стабильными и обеспечивать высокую точность по временной и пространственным переменным. С ростом интенсивности лазерных импульсов схема Кранка-Николсона становится неэффективной для аппроксимации оператора эволюции [89], также как и традиционные алгоритмы ее расщепления по физическим параметрам [90]. Поэтому в разделе 3.3. для решения нестационарного уравнения Шредингера разработаны алгоритмы расщепления, основанные на унитарных разложениях оператора эволюции, аппроксимирующих импульсы электрического поля более короткие, чем период классической орбиты электрона в атоме. Этот режим представляет значительный концептуальный интерес, так как перекидывает мост между фотоионизацией и ударной ионизацией атома частицами высоких энергий (см. раздел 2.2). Даны примеры расчета динамики населенностей атома водорода в импульсном резонансном лазерном поле и постоянном магнитном поле, для которых возможны эффекты стабилизации.

Общая характеристика работы.

Актуальность проблемы.

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

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

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

Во введении было показано, что для построения дискретных моделей и комплексных алгоритмов применяются алгебры динамической симметрии, канонические преобразования, символьные и численные методы. Использование символьных методов позволяет избежать ошибок округления, получить решение поставленной задачи для необходимой области определения исходных параметров, а также рассматривать сингулярные функции и изучать их асимптотическое поведение. Однако промежуточные результаты символьных вычислений, как правило, имеют большой объём, а их преобразование к искомой форме требует значительных компьютерных ресурсов. Поэтому, необходима разработка эффективных комплексных алгоритмов, сочетающих достоинства символьных и численных методов и позволяющих визуализировать исследуемые решения, зависящие от набора физических параметров. Для тестирования разрабатываемых алгоритмов и контроля эффективности их работы необходимо построение эталонных моделей [91].

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

В связи с этим были поставлены и решены следующие задачи по разработке, анализу и апробации комплексных дискретных моделей и алгоритмов:

• нормализации и квантования полиномиальных гамильтонианов.

• анализ моделей ионизации в трехчастичных квантовых системах.

• решения задачи рассеяния плоских волн на системе квантовых точек.

• расчета населённостей водородоподобных атомов во внешних импульсных полях.

• построения системы эталонных моделей.

• визуализации физических характеристик изучаемых процессов.

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

Практическая значимость. Разработаны комбинированные программы символьного решения задач нормализации и квантования полиномиальных гамильтонианов, реализованные на языках MAPLE, MATHEMATICA, REDUCEпроизведено сравнение эффективности алгоритмов для нормализации полиномиальных гамильтонианов классической механики методами Биркгофа-Густавсона и Депри-Хори, а также их различных реализаций на указанных языках. Вычислены мультипольные поправки к энергетическому спектру и волновой функции для двухи трехмерного водородоподобного атома в однородном и неоднородном электрических полях, реализована модель среднего поля для вычисления сечений процесса ионизации атома гелия. Для нестационарного уравнения Шредингера разработаны программы, реализованные на языках FORTRAN и MAPLE, для расчета динамики населенностей атома водорода в импульсном резонансном лазерном поле и постоянном магнитном поле, для которых возможны эффекты стабилизации. Применение технологии комплексных символьных и численных алгоритмов включает новые возможности визуализации искомых решений, в частности, сложных движений и эволюции волновых пакетов в классических и квантовых системах, а также построения эталонных моделей для тестирования алгоритмов и контроля эффективности их работы. Разработанную комбинированную методику и программы можно применять для расчета характеристик квантовых систем во внешних полях.

Работы выполнялись в рамках совместных проектов по теоретическим и экспериментальным исследованиям квантовых систем во внешних полях, поддержанных грантами РФФИ: 00−02−16 337-а «Исследование эффектов, связанных с явлением динамического хаоса при прохождении частиц большой энергии через кристаллы», 00−02−81 023-Бел2000а «Интегрируемые модели в квантовой механике и квантовой оптике», 03−02−16 263-а «Явление динамического хаоса во взаимодействии частиц большой энергии с кристаллическими структурами» и ФНТП РФ 40.018.1.1.1314 (2003;2004 гг.) «Разработка новых методов диагностики биологических структур на основе прецизионных оптических и спектроскопических измерений» .

Апробация результатов. Результаты диссертационной работы докладывались на 19 международных конференциях:

III — VI International Workshops on Computer Algebra in Scientific Computing, (Самарканд и Бухара, 5−11 октября 2000 г.- Констанц, 22−26 сентября, 2001 г.- Ялта, 22−27 сентября, 2002 г.- Пассау, 20−26 сентября, 2003 г.) — 6th IMACS Conference on Applications of Computer Algebra (Санкт-Петербург, 25−29 июня 2000 г.) — 2 международная конференция «Современные направления вычислительной физики», (Дубна, 24−29 июля 2000 г.) — XXIII International Colloquium Group Theoretical Methods In Physics (Дубна, 31 июля-5 августа 2000 г.) — международная конференция «Дифференциальные Уравнения и Системы Компьютерной Алгебры», (Брест, 18−23 сентября 2000 г.) — International Workshop on Computer Algebra and its Application to Physics, (Дубна, 28−30 июня 2001 г.) — VII Всероссийская научно-техническая конференция «Новые информационные технологии» (24 апреля 2002 г.- 24−25 марта 2004 г.). Sixth Workshop on Computer Algebra (Дубна, 27−28 мая 2002 г.) — VIII International Workshop on Advanced Computing and Analysis Techniques in Physics Research (Москва, 2428 июня 2002) — Fifth International Congress on Mathematical Modelling (Dubna, 30 сентября — 2 октября 2002 г.) — 12-th International Colloquium Quantum Groups and Integrable Systems (Прага, 12 — 14 июня 2003 г.) — XXIII International Conference on Photonic, Electronic and Atomic Collisions, (Стокгольм, 23−29 июля 2003 г.) — Focus Symposium on «Quantum Physics and Communication» (Дубна, 30 июля — 2 августа 2003 г.) — X International Conference on Symmetry Methods in Physics (Ереван, 13−19 августа 2003 г.) — International Workshop on Laser Physics and Photonics (Саратов, 1−4 октября, 2002 г.- 2−7 октября, 2003 г.) — а также на научных семинарах Научного центра прикладных исследований (НЦеПИ), Лаборатории информационных технологий (ЛИТ), Лаборатории теоретической физики (ЛТФ) Н. Н. Боголюбова Объединённого института ядерных исследований, физико-математического факультета Белгородского государственного университета (БелГУ), НИИЯФ МГУ, РУДН.

Публикации. Основное содержание диссертации отражено в 16 публикациях [92 — 107] в виде статей в журналах Программирование, Вычислительные методы и программирование, Журнал физической химии, Journal of physics A, Physics Letters A., SPIE, Journal of Computational Methods in Sciences and Engineering, в сборнике научных трудов «Новые информационные технологии», докладов в трудах международных конференций, препринтов и сообщений ОИЯИ, программы для расчета динамических характеристик системы ангармонических осцилляторов и водородоподобных атомов во внешних полях официально зарегистрированы в Роспатенте [108 -110]. Структура и объём диссертации. Диссертация состоит из введения, трёх глав, заключения и семи приложений. Объём диссертации — 160 страниц, 24 рисунка, 8 таблиц.

Список литературы

включает 141 наименование.

Заключение

.

В диссертационной работе разработаны дискретные модели, комплексные алгоритмы и программы для анализа и расчета динамических характеристик квантовых систем во внешних полях.

Актуальность темы

исследований определяется известными проблемами редукции, понижения размерности и построения устойчивых и экономичных схем численного решения параметрических спектральных и эволюционных краевых задач для многомерного уравнения Шредингера, возникающих при моделировании процессов ударной ионизации атомов, импульсного воздействия лазерных полей на атомы в магнитных ловушках, рассеяния света на системе квантовых точек, генерации второй гармоники оптического атома и трехфотонного рассеяния, как типичных представителей квантовых систем во внешних полях, имеющих важное прикладное значение. Например, их можно рассматривать как прототипы возможных физических процессов в интерференционной микроскопии [139] и в элементах квантового компьютера [12], или магнитных ловушках в эксперименте ATHENA CERN по созданию низкотемпературной плазмы из позитронов и антипротонов и сверхточной спектроскопии антипротонных атомов [141].

В диссертации для формулировки дискретных моделей и адаптивных алгоритмов использовались канонические преобразования, неприводимые представления алгебр динамической симметрии квантовой системы в сочетании со схемами квазиклассического квантования и теории возмущений, методами компьютерной алгебры и численными методами, такими как расщепления оператора эволюции по физическим параметрам, приведения к обыкновенным дифференциальным уравнениям Канторовича, алгебраической редукции Галёркина и конечных элементов.

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

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

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

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