Моделирование равновесных свойств квантовых систем методом Монте-Карло в расширенных ансамблях тема автореферата и диссертации по физике, 01.04.02 ВАК РФ

САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

ВОЗНЕСЕНСКИЙ Михаил Андреевич

МОДЕЛИРОВАНИЕ РАВНОВЕСНЫХ СВОЙСТВ КВАНТОВЫХ СИСТЕМ МЕТОДОМ МОНТЕ-КАРЛО В РАСШИРЕННЫХ АНСАМБЛЯХ

Специальность 01.04.02 —теоретическая физика

АВТОРЕФЕРАТ диссертации на соискание учёной степени кандидата физико-математических наук

003407652

На правах рукописи

1 О ДЕК 2009

Санкт-Петербург 2009

003487652

Работа выполнена на кафедре статистической физики физического факультета Санкт-Петербургского государственного университета.

Научный руководитель: доктор физико-математических наук,

профессор ВОРОНЦОВ-ВЕЛЬЯМИНОВ Павел Николаевич

Официальные оппоненты: доктор физико-математических наук,

профессор ПИСЬМАК Юрий Михайлович

доктор физико-математических наук, профессор КУЗЬМИН Владимир Леонидович

Ведущая организация: Объединенный институт высоких температур

Российской академии наук

О

часов на заседа-

Защита состоится « 2009 года в '

нии Совета Д 212.232.24 по защите докторских и кандидатских диссертаций при Санкт-Петербургском государственном университете по адресу: 199034, Санкт-Петербург, Средний пр. В. О., д. 41/43, ауд

С диссертацией можно ознакомиться в библиотеке Санкт-Петербургского государственного университета.

Л' » ]

Автореферат разослан ноября 2009 года.

Учёный секретарь диссертационного совета д. ф.-м. н., проф. Щёкин А.К.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы. Значение и возможности компьютерного моделирования в исследовании различных квантовых систем за последние время заметно возросли. Это связано с существенным ростом производительности вычислительных систем и развитием новых эффективных подходов.

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

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

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

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

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

Апробация работы. По материалам диссертации были сделаны доклады на следующих конференциях: International Conference on Strongly Coupled Coulomb Systems (20-25 июня 2005, Москва), «XIII Симпозиум по межмолекулярному взаимодействию и конформациям молекул» (19—23 июня 2006 года, Санкт-Петербург), International conference "PNP12 Physics with Non-Ideal Plasmas" (3-8 сентября 2006, Darmstadt, Германия), Conference on Computational Physics CCP2007 (5-8 сентября 2007, Brussels, Бельгия), XIII International Conference on Physics of Non-Ideal Plasmas (13 -18 сентября 2009, Черноголовка).

Структура и объём диссертации. Диссертация состоит из введения, четырёх глав, заключения и списка цитируемой литературы. Объём работы составляет 117 страниц, включая 48 рисунков и 6 таблиц. Список литературы содержит 59 наименований.

Публикации: По материалам диссертации опубликовано 6 печатных работ. Список публикаций приведён в конце автореферата.

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

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

В первой главе изложены методы, лежащие в основании подхода. Вначале выписывается общее выражение для статистической суммы системы тождественных частиц [1], далее изложен метод Монте-Карло (МК) интегралов по траекториям [2] и выписаны выражения для статистических сумм в вершинном приближении [3]. Описываются метод энтропического моделирования [4], связанный с ним алгоритм Ванга — Ландау (ВЛ) [5], и метод расширенного ансамбля [6].

Для системы N тождественных частиц, подчиняющихся статистике Ферми, статистическая сумма задаётся в виде суммы по классам (7 перестановок частиц в матрице плотности системы различимых частиц:

СОДЕРЖАНИЕ ДИССЕРТАЦИИ

*** {Р)=10м о?).

Здесь /} = (квТ) 1 — обратная температура в энергетических единицах, V с

= + " — спиновая часть статистической суммы, С„ - число циклов длинной V в данной перестановке, т.о. есть полное число

циклов в классе,

число элементов в С,

— координатная часть статистической суммы, соответствующая перестановке из класса б.

Статистические суммы представляются в виде интегралов по

траекториям в «-вершинном приближении (записано в безразмерных переменных):

, ( п и п \

ад-Ш к-кЦ-Р^)'-^""-'»]'

где гп+1 = — условие замыкания. Потен-

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

Метод расширенного ансамбля [6] позволяет в ходе одного расчёта получить отношение статистических сумм исследуемой и опорной систем. Пусть И = -/ЗН(д) есть приведённая конфигурационная часть классического гамильтониана Н(д), а /¡0 — гамильтониан системы, для которой статистическая сумма определена точно. Строится цепочка гамильтонианов \,\,..ки =И, в которой с изменением индекса т от 0 до М гамильтониан \ постепенно превращается в Ии, например, по линейному закону кт = (1 - Хт )й„ + ХтИм, где О = Х0 < Л, <... < Хм = 1. Каждому т соответствует конфигурационный интеграл 2т = |фехр(/гт) в каноническом ансамбле. Составляется расширенный и модифицированный ансамбль со статистической суммой:

м

т=0

где 5т — так называемые балансировочные параметры. В расширенном ансамбле организуется МК-блуждание согласно стандартному алгоритму Метрополиса. В ходе описанной МК-процедуры подсчитывается число раз

6

пт, которое система побывала в подансамбле с индексом т. В результате получается оценка вероятности состояния характеризуемого значением %т: рт-пт1 пмс, где пмс — полное число МК-шагов. С другой стороны,

Рт ~12, и, следовательно

7-7 лм о е

"о

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

Во второй главе в начале описывается использование ВЛ-алгоритма в методе расширенных ансамблей. Следуя алгоритму, вводятся два массива размером М+1, один для весов , второй для чисел посещения пт. В начале расчёта все 5т и пт берутся равными нулю. В расширенном ансамбле запускается блуждание. Вероятность перехода определяется выражением

А—и = тт(1,ехр(/гт±1 -Ит- 5т±1 + 5И)). Вне зависимости от того, состоялся переход или нет, в счётчик пт, соответствующий текущему подансамблю, добавляется единица, а значение 5„ увеличивается на величину АБ. Таким образом, последним действием мы уменьшаем вес текущего подансамбля, а, следовательно, вероятность попасть в него в следующий раз. Описанным образом выполняется серия, обычно порядка 106-И07 шагов. Такого количества вполне достаточно, по крайней мере, для рассматриваемых в этой работе систем, для формирования равномерного распределения чисел посещения подансамблей пт. Далее, величина ЛБ домножается на константу а < 1 и запускается новая серия, тем самым начинается более тонкая настройка балансировочных параметров.

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

6 в ь

Рис. 1. Зависимость от обратной темпера- Рис. 2. Зависимость от обратной температуры фактора /2 системы двух невзаимо- туры энергии двух невзаимодействующих действующих бесспиновых частиц в гар- бесспиновых частиц в гармоническом по-моническом поле. Сплошная линия — ле- Сплошная линия — точные формулы, точные формулы, штрихованная - ко- штрихованная - конечномерное прибли-

нечномерное приближение.

жение. Е0 =2 для с!=\, Е0 =4 для <) = Ъ.

для которой существуют точные выражения для равновесных средних [7]. Для N = 2 имеется два класса перестановок, которые обозначены как 20 (две отдельные частицы) и 01 (цикл из двух частиц), Л/(20) = Л/(01) = 1.

Статистическая сумма 2^ = \210-2д{\12 представляется в виде 2^=220/2, где /2=[\-2011210\12, и производится расчёт отношения 2М / 210 методом расширенного ансамбля с настройкой балансировочных параметров по ВЛ-алгоритму. Использовано 5-и вершинное приближение, расширенный ансамбль составлен из 2-х подансамблей. Выполнены расчёты для системы в одномерном и трёхмерном случаях и проведено сравнение результатов с точными [7] (сплошные линии) и конечномерными [8] (штрихованные линии) зависимостями (рис. 1). Хорошее согласие результатов наблюдается до ¿ = 10 (Ь = рЬа>). При более низких температурах отличие отношений статистических сумм от единицы становится меньше точности расчётов. Возможность вычисления канонических средних показана на примере средней энергии. Энергия вычислена двумя способами — усреднением эстиматоров энергии (треугольники) и путём вычисления ко-

Рис. 3. Зависимость от обратной темпера- Рис. 4. Зависимости от обратной температуры статистических сумм системы ЭНСРГИИ »У* тождественных бесспиновых частиц с кулоновским отгалки-

двух тождественных бесспиновых частиц

вашем в гармоническом поле. Энергии

с кулоновским отталкиванием в гармони-

основных состоянии см. текст.

ческом поле.

нечно-разностной производной от логарифма статистической суммы (кружки) (рис. 2).

Удаётся проследить выход кривых на горизонтальный уровень основного состояния. Аналогичные расчёты проведены для системы трёх частиц. Получены следующие оценки для энергии основного состояния: для ¿ = 1, 5 = 1/2 — Е0 =2,47 + 0,04 (точное значение Е0 =2,5), для с! = Ъ, 5=1/2 — Е0 = 5,54 ± 0,03 (точное значение Е0 = 5,5).

В разделе 6 выполнены расчёты для двух тождественных частиц с кулоновским отталкиванием в гармоническом поле. Статистическая сумма

представляется в виде 2^ ^¿"о^о/ 2^-2д^20112й0^12 и вычисляются отношения 220 / 2°0 и 20112\х. В качестве опорной системы использована система без взаимодействия. Расширенный ансамбль составлен из 10-и подансамблей. Проведены расчёты в 5-и вершинном приближение в одномерном и трёхмерном случае (рис. 3). Получена оценка энергии основного состояния для с} = 1, 5 = 0 Е0 =2,82±0,02; для <1 = 1, 5 = 1/2 — Ед= 2,53±0,01; для ¿=3, 5 = 0 — Е0 =4,52±0,03; для ¿ = 5 = 1/2 —

N

/г = 1/2

0,4 0,6 0,В 1,0 1,2 1,4 1,6 1,8 2,0 2,2

ь

Рис. 5. Статистические суммы /?01 и Рис. 6. Отношение системы

?(А) „„„ „„„„.«,.„, тл^п^тт,» N - 5 тождественных невзаимодействующих частиц в трёхмерном гармоническом поле в сравнении с точными зависимостями.

Для системы двух тождественных частиц с кулоновским отталкиванием в кулоновском поле притяжения заряда 2.

Е0 =4,13 ±0,01 (рис. 4). Также выполнены аналогичные расчёты для системы трёх частиц с кулоновским отталкиванием в трёхмерном гармоническом поле, получены следующие оценки для 5 = 0 — Еа =12,28 + 0,04, для л- = 1 / 2 — £0 = 1 1,94 ±0,02.

В разделе 7 рассмотрена система двух фермионов с кулоновским отталкиванием в кулоновском поле притяжения центрального заряда 2 (атом гелия) (рис. 5). В качестве опорной системы использована система без взаимодействия в гармоническом поле. Использовано 100 вершинное приближение, расширенный ансамбль составлен из 22-х подансамблей. Полученные значения энергии основного Е0 =-2,83±0,05 и триплетного Е, =-2,03 ±0,07 состояний удовлетворительно согласуются с литературными данными Ей =-2,904 и Е0 =-2,175 [9]. Для решения проблемы сингулярности кулоновского потенциала использовано приближение Кола — Де Рэдта [10]. В рамках этого приближения повторены расчёты для системы двух частиц с кулоновским отталкиванием в одномерном гармониче-

2 20 N

£ 15

10

0,1 1,0 1,2 1,4 1.6 1,8 2,0 и

ь

и.

0,0 О X 0,4 0,6 О,В 1,0

Рис. 7. Построение аппроксимации Рис. 8. Построение аппроксимации -1п г[А) = Е0Ь + С/Ь для системы N = 5 Р-=Ей + А12 для системы N = 5 тождест-тождественных частиц с кулоновским от- венных частиц с кулоновским отталкива-талкиванием в трёхмерном гармоническом нием в трёхмерном гармоническом поле, поле.

ском поле. Уточнённое значение Е0 =2,70 ±0,02 оказывается в лучшем согласии с данными из [11] (£0 = 2,72).

В третьей главе предлагается обобщение метода на случай большего числа частиц. Вклады в статистическую сумму разделяются на две группы — положительную и отрицательную и применяется подход, аналогичный изложенному во второй главе. Вначале выполнено тестирование метода на системе N = 5 невзаимодействующих частиц в трёхмерном гармоническом поле. Использовано 5-ти вершинное приближение, расширенный ансамбль составлен из 10-и подансамблей. Проведён расчёт отношения 2^12^ (рис. 6).

Далее, проведены расчёты для системы N = 5 частиц с кулоновским отталкиванием в трёхмерном гармоническом поле. Статистическая сумма представлена в виде = и проведены вы-

числения отношений 2^ / и 2^ / 2\_. В качестве опорной взята система без взаимодействия. Для зависимости статистической суммы от обратной температуры использована аппроксимация -1п 2^=ЕйЬ + С!Ь

Рис. 9. Низкотемпературная термодина- Рис. 10. Свободные энергии систем = 6

мика для системы N-5 фермионов и N = 1 фермионов (5 = 1/2) с кулонов-

(^ = 1/2) с кулоновским отталкиванием в ским отталкиванием в трёхмерном гармо-

гармоническом поле. ническом поле.

(рис. 7), а для зависимости свободной энергии от температуры —

Р = Ч\пг[л) = Е0 + А?, где / = 1 /Ь (рис. 8).

Предложенные аппроксимации дают значения энергии основного состояния: Ед = 17,6 ± 0,1 для 5 = 0 и Ед =15,94±0,04 для $ = 1/2. На основе зависимости свободной энергии от температуры построена низкотемпературная термодинамика системы (рис. 9). При использовании параболической аппроксимации для свободной энергии энтропия и теплоёмкость совпадают. Аналогичные расчёты произведены для системы N = 6 и N = 1 частиц. Получены следующие значения энергий основных состояний: для N = 6 £0 = 24,6 ±0,3 для 5 = 0 и Е0 =21,25 + 0,05 для 5 = 1/2; для N = 7 £0 = 29,9 ±0,5 для 5 = 0 и £0=27,6±0,4 для 5 = 1/2 (рис. 10).

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

Выводы. В работе предлагается подход к изучению квантовых систем, основанный на сочетании метода интегралов по траекториям и метода Монте-Карло в расширенном ансамбле с настройкой весов по алгоритму Ванга — Ландау.

Предложенный подход даёт достаточно точные результаты при сравнительно небольших затратах машинного времени. Типичное время расчёта для одного значения температуры составляет от получаса в случае системы двух невзаимодействующих частиц до суток в случае атома гелия или пяти частиц с кулоновским отталкиванием в гармоническом поле на компьютере уровня Репйит IV (ЗГГц).

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

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

Выносимые на защиту основные результаты диссертации:

1. Разработан подход к изучению температурных свойств квантовых систем, сочетающий метод Монге-Карло интегралов по траекториям и метод расширенных ансамблей и дающий возможность ослабить проблему знака.

2. Проведено тестирование алгоритма на системах 2, 3 и 5 невзаимодействующих ферми-частиц в гармоническом внешнем поле.

3. Выполнены расчёты для 2, 3 и 5 - 7 ферми-частиц с кулоновским отталкиванием в гармоническом внешнем поле, а также для атома гелия.

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

Содержание диссертации изложено в следующих публикациях:

1. Р. N. Vorontsov-Velyaminov, M. A. Voznesenski, D. V. Malakhov, A. P. Lyubartsev and A V Broukhno. Path integral method in quantum statistics problems: generalized ensemble Monte Carlo and density functional approach. J. Phys. A: Math. Gen., 2006, v.39, Pp. 4711-4716.

2. M. А. Вознесенский, П. H. Воронцов-Вельяминов, А. П. Любарцев. Расчеты равновесных свойств квантовых систем методом Монте-Карло в расширенном ансамбле. // Вычислительные методы и программирование, 2008, Раздел 1, 170-183 (http://num-meth.srcc.msu.ruy).

3. Вознесенский М. А., Воронцов-Вельяминов П. Н., Любарцев А. П. Расчёты равновесных свойств квантовых систем с кулоновским взаимодействием методом Монте-Карло в расширенном ансамбле. // Вестник Санкт-Петербургского университета. Серия 4 физика и химия, выпуск 2, июнь 2009, 3-16.

4. Vorontsov-Velyamitiov P. N., Voznesenski M. A., Malakhov D. V., Lyubartsev A. P., Broukhno A. V. Path integral method in quantum statistics problems: Generalized ensemble Monte Carlo and density functional approach. // International Conference on Strongly Coupled Coulomb Systems. June 20-25 2005, Moscow, Russia. Book of abstracts, p. 95.

5. E. А. Поляков, M. А. Вознесенский, А. В. Брухно, А. П. Любарцев, П. H. Воронцов-Вельяминов. Квантово-статистический метод интегралов по траекториям. Расчёты равновесных свойств и возбужденных состояний квантовых систем методами Монте-Карло и функционала плотности. // XIII Симпозиум по межмолекулярному взаимодействию и конформациям молекул 19-23 нюня 2006, Санкт-Петербург. Тезисы докладов, с. 131.

6. Vorontsov-Velyaminov P. N., Voznesenskiy М. A., Polyakov Е. А., Lyubartsev А. P. Calculation of canonical properties and excited states by path integral numerical methods // XIII International Conference on Physics of Non-Ideal Plasmas. September 13-18 2009, Chernogolovka, Russia. Book of abstracts, p. 18.

Цитируемая литература:

1. Фейнман P. Статистическая механика. Курс лекций. М.: Мир. 1975.

2. Фейнман Р., Хиббс А. Квантовая механика и интегралы по траекториям. М.: Мир, 1968.

3. Замалин В. М., Норман Г. Э., Филинов В. С. Метод Монте-Карло в статистической термодинамике. М.: Наука, 1977.

4. Lee J. New Monte-Carlo algorithm: Entropic sampling // Phys. Rev. Lett., 1993, v. 71(2), p. 211-214.

5. Wang F., Landau D. P. Efficient, multiple-range random walk algorithm to calculate the density of states. Phys. Rev. Lett., 2001, v. 86, p. 2050-2053.

6. Lyubartsev A. P., Martsinovski A. A., Shevkunov S. V., Vorontsov-Velyaminov P. N. New approach to Monte Carlo calculation of the free energy: Method of expanded ensembles // J. Chem. Phys. 1992. v. 96(3), p. 1776-1783.

7. Vorontsov-Velyaminov P. N., Nesvit M. 0., Gorbunov R. I. Bead-Fourier path-integral Monte Carlo method applied to systems of indentical particles // Phys. Rev. E. 1997. v. 55(2), pp. 1979-1997.

8. Vorontsov-Velyaminov P. N., Lyubartsev A. P. Entropic sampling in the path integral Monte-Carlo method // J. Phys. A: Math, and Gen., 2003, v. 36(3), pp. 685-693.

9. C.L. Pekeris, l'S, 2'S and 2>S States of H~ and of He // Phys. Rev., 1962, v. 126(4), pp. 1470-1476.

10. Kole J. S., De Raedt H. Quantum Monte Carlo Method for Attractive Coulomb Potentials // Phys. Rev. E., 2001, v. 64, pp. 016704-1-016704-6.

11. ПоляковЕ. А., Воронцов-Вельяминов П. H. Квантовый газ во внешнем поле при конечных температурах. Точное выражение для плотности и возбужденный состояния // Вычислительные методы и программирование. 2007.Раздел I.e. 334-351. (http://num-meth.srcc.msu.ru/).

Подписано к печати 18.11.09. Формат 60 х 84 1Лб. Бумага офсетная. Гарнитура Тайме. Печать цифровая. Печ. л. 1,0. Тираж 100 экз. Заказ 4553. Отпечатано в Отделе оперативной полиграфии химического факультета СПбГУ. 198504, Санкт-Петербург, Старый Петергоф, Университетский пр., 26. Тел.: (812) 428-4043,428-6919.

Введение

1 Методы и алгоритмы

1.1 Выражение для статистической суммы системы тождественных частиц.

1.2 Метод интегралов по траекториям.

1.3 Метод энтропического моделирования.

1.4 Алгоритм Ванга — Ландау.

1.5 Метод расширенного ансамбля.

2 Метод расширенного ансамбля для малых квантовых систем

2.1 Алгоритм Ванга — Ландау в методе расширенных ансамблей

2.2 Эстиматоры.

2.3 Точные и конечномерные выражения для систем невзаимодействующих частиц в гармоническом поле.

2.4 Система невзаимодействующих частиц

2.4.1 Две частицы

2.4.2 Три частицы

2.5 Вычисление энергии

2.6 Система тождественных частиц с кулоновским отталкиванием в гармоническом поле.

2.6.1 Две частицы

2.6.2 Три частицы

2.7 Системы в кулоновском поле притяжения.

2.7.1 Аппроксимация для кулоновского потенциала

2.7.2 Две частицы с кулоновским отталкиванием в гармоническом поле.

2.7.3 Атом гелия.

3 Обобщение метода

3.1 Системы многих частиц.

3.2 Система невзаимодействующих частиц в гармоническом поле.

3.3 Система частиц с кулоновским отталкиванием в гармоническом внешнем поле.

3.3.1 Пять частиц.

3.3.2 Шесть и семь частиц.

3.4 Параллелизация вычислений.

4 Метод плотности состояний

4.1 Квантовая частица в гармоническом поле.

4.2 Система двух различимых квантовых частиц с кулоновским отталкиванием в гармоническом поле .ЮЗ

В далёком 1965 году, на заре компьютерной эры, Гордон Мур, впоследствии один из основателей корпорации Intel, обратил внимание на интересную закономерность в развитии компьютеров. Он заметил, что новые модели микросхем разрабатывались спустя более или менее одинаковые периоды после появления их предшественников, а ёмкость их при этом возрастала каждый раз примерно вдвое. Если такая тенденция продолжится, заключил Мур, то мощность вычислительных устройств со временем будет экспоненциально возрастать. Данная закономерность известна как закон Мура и оказывается справедливой и по сей день. Столь стремительный рост мощностей вычислительных систем открывает поистине фантастические возможности для компьютерного моделирования и численных эскпериментов, зарекомендовавших себя как один из верных способов решения задач из всех областей знаний.

Одним из наиболее мощных, но зачастую требовательных к ресурсам, стохастических методов является метод Монте-Карло. Идея использовать случайные явления в области приближённых вычислений возникла ещё в XIX веке, когда появилась работа Холла об определении числа 7г путём бросания иглы на разграфлённую параллельными линиями бумагу. Однако, годом рождения метода Монте-Карло принято считать 1949-й, когда вышла статья Метрополиса и Улама [1], содержащая основания и общее описание статистического подхода к решению интегральных и диффренциальных уравнений, возникающих в различных областях естественных наук.

Ранние методы квантового моделирования были направлены на решение уравнения Шредингера или уравнения Блоха. Они включают в себя такие методы, как метод Хартри-Фока, метод функционала плотности, диффузионный метод Монте-Карло и другие. Однако все эти методы позволяют найти лишь волновые функции чистых состояний, обычно основного, и не дают возможности получить теромодинамику квантовой системы при конечной температуре. Основанный на Фейнмановской трактовке квантовой механики [2, 3] метод Монте-Карло интегралов по траекториям позволяет решить эту проблему путем сведения статистической суммы квантовой системы к статистической сумме чисто классической системы. Используя групповое свойство, матрицу плотности можно представить в виде произведения высокотемпературных матриц плотности, для которых уже возможно использовать различные высокотемпературные приближения. В результате, в простейшем случае метода вершин квантовая частица представляется в виде траектории — замкнутого «полимера» с гармоническими связями между вершинами. Существенная положительная особенность данного подхода заключается в том, что здесь без каких-либо приближений полностью учитываются все межчастичные взаимодействия, и как следствие все корреляции.

Первым, кто использовал метод Монте-Карло для численного вычисления статистической суммы квантовой системы был Фосдик [4]. В его работе статистическая сумма записана в терминах интеграла Винера, который выражается через n-кратный интеграл Римана, который, в свою очередь, вычисляется методом Монте-Карло. Фосдик демонстрирует свой подход на двумерной системе двух тождественных взаимодействующих частиц в ящике и даже выписывает выражения для системы трёх частиц. Но в следующих работах [5, 6] он делает основной упор на использовании интегралов по траекториям для расчёта двухчастичной слэтеровской суммы в Не4 при низких температурах. Следующей в середине 70-х годов была серия работ группы Нормана и Филинова по плазме и ферми-системам [7, 8, 9, 10]. Фактически они первые ввели понятие метода вершин. Чуть позже, в 1979 году, Баркер [11] в рамках метода вершин провёл расчёт для энергии и функции распределения систем одной и двух различимых частиц в ящике. Эти авторы считаются основоположниками метода Монте-Карло интегралов по траекториям.

Существует большое число результативных исследований с использованием метода Монте-Карло интегралов по траекториям [12, 13], но, к сожалению на данный момент указанный метод успешно применялся лишь для ферми-систем с малым числом степеней свободы или для систем, удовлетворяющих статистике Бозе. Это связано с тем, что для ферми-частиц волновая функция, а равно и матрица плотности меняет знак (антисимметрична) при перестановке любых двух частиц. Это приводит к тому, что веса траекторий оказываются знакопеременными и не могут быть непосредственно интерпретированы как вероятности. Более того, при понижении температуры разность положительных и отрицательных вкладов в статистической сумме и в выражениях для средних экспоненциально уменьшается и становится трудно различимой на фоне статистического шума. Указанная проблема известна как проблема знака [14].

Проблема знакопеременности веса решалась путем перехода к усреднению по положительной весовой функции, в простейшем случае — по абсолютному значению [15, 16]: но проблема знака всё равно оставалась, выражаясь в стремлении (sgn(W)) + к drOW =

JdrOsgn (W) \W\ fdr\W\ (Osgn(MQ) /rfrsgnCW)!^! ' fdT\W\~ (sgn(H0)4

1) нулю.

Одна из первых попыток использовать метод Монте-Карло интегралов но траекториям для исследования термодинамических свойств систем тождественных частиц принадлежит Такахаши и Имаде [17, 18, 19]. Они предложили в методе вершин симметризовать (антисимметризовать) каждую высокотемпер-татурную матрицу плотности в разложении, таким образом на каждом шаге мнимого времени происходит учёт всех перестановок частиц. В качестве вероятности в методе существенной выборки Такахаши и Имада использовали абсолютное значние веса траектории, а учёт знаков производили при накоплении канонических средних. Данный подход, в силу геометрических особенностей конфигурационного пространства, позволил даже решить проблему знака в одномерном случае, но для систем с большей размерностью он оказался малоэффективен.

Через некоторое время, Холл [20, 21, 22, 23] предлагает ввести проекционный оператор, который исключает из рассмотрения траектории с весом меньше ошибки конечномерного приближения. Подход позволяет заметно уменьшить дисперсию ошибки вычислений методом Монте-Карло, но поскольку поверхности в фазовом пространстве, на которых веса траекторий равны пулю, неизвестны, приходится использовать приближенную проекцию.

Ныоман и Куки в 1992 году [24] разработали метод для ослаблания проблемы знака в системе двух фермионов, в котором на каждом МК-шаге усредняется несколько эквивалентных траекторий, полученных путём вращения исходной. Данный подход является численно-точным, но сформулировать его для большего числа частиц пока не удалось.

В это же время Цеперли предпринял попытку решить проблему знака, введя так называемый «метод интегралов по траекториям с ограничениями» [25, 26], идея которого отчасти перекликается с идеями Холла. Суть подхода заключается в исключении из усреднения отрицательных слагаемых статистической суммы путем ограничения области интегрирования. Этот метод действительно снимает проблему знака, но поверхности, ограничивающие область интегрирования, известны лишь приблизительно, они становятся крайне сложными в случае многочастичных систем с сильным взаимодействием, и правомочность такого ограничения оказывается спорной. Кроме того, как было показано Фи-линовым [27], данный метод приводит к неправильному поведению при низких температурах простой системы, имеющей аналитическое решение.

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

Следует также упомянуть работы других авторов: Карлсона и Калоса [30], предложивших концепцию зеркальных потенциалов, Коржениовского и др. [31], указавших на возможность использования метода Фейнмана-Каца для получения основного состояния многочастичных систем, Жанга [32], исследовавшего ферми-системы при конечных температурах в большом каноническом ансамбле, Корнея и Драммонда [33, 34], записавших гауссовы представления для ферми систем в фазовом пространстве и ряд других [35, 36, 37, 38].

Уже было отмечено, что метод Монте-Карло интегралов по траекториям позволяет перейти от квантовой к формально классической системе. Это открывает возможность использовать всю совокупность классических методов моделирования для решения квантовых задач.

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

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

Заключение

В работе предлагается подход к изучению квантовых систем, основанный на сочитании метода интегралов по траекториям и метода Монте-Карло в расширенном ансамбле с настройкой весов по алгоритму Ванга — Ландау.

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

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

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

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

Предложенный подход даёт достаточно точные результаты при сравнительно небольших затратах машинного времени. Типичное время расчёта для одного значения температуры составляет от получаса в случае системы двух невзаимодействующих частиц до суток в случае атома гелия или пяти частиц с кулоновским отталкиванием в гармоническом поле на компьютере Pentuim IV (3 ГГц).

В работе использовалось простейшее представление для высокотемпературной матрицы плотности. Подход не накладывает никаких ограничений на использование более сложных представлений для матрицы плотности, например [52, 53, 54, 55], на число вершин и типы МК шагов. Также возможно использование различных техник, ускоряющих сходимость интегралов по траекториям [56, 57, 58, 59]. Однако возможности комбинирования указанных подходов требуют отдельного исследования. Кроме этого нет никаких ограничений на размерность и на потенциалы взаимодействия в системе. Расчёты могут быть параллелизованы, что очень важно при моделировании многочастичных систем.

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

1. N. Metropolis and S. Ulam. Journal of the American Statistical Association, 44 (247):335-341, 1949.

2. P. Фейнман, А. Хибс. Квантовая механика и интегралы по траекториям. «Мир», Москва. 1968.

3. Р. Фейнман. Статистическая механика. «Мир», Москва. 1975.

4. L. D. Fosdick. Numerical estimation of the partition function in quantum statistics. Journal of Mathematical Physics, 3(6):1251-1264, 1962.

5. L. D. Fosdick and H. F. Jordan. Path-integral calculation of the two-particle slater sum for he4. Phys. Rev., 143(1):58-66, Mar 1966.

6. H. F. Jordan and L. D. Fosdick. Three-particle effects in the pair distribution function for he4 gas. Phys. Rev., 171(1):128-149, Jul 1968.

7. В. M. Замалин, Г. Э. Норман. Исследование термодинамики квантовых жидкостей методом Монте-Карло. XV Всесоюз. совещ. по физике низких температур. Тезисы докл. Тбилиси, ИФАН ГССР, 1968.

8. В. М. Замалин, Г. Э. Норман. О методе Монте-Карло в фейнмановской формулировке квантовой статистики. ЖВМиМФ, т. 13(2), 1973.

9. В. С. Филинов. Учёт вырождения электронов в псевдопотенциальной модели неидеальной плазмы. ТВТ, т.11(4), 1973.

10. В. М. Замалин, Г. Э. Норман, В. С. Филинов. Метод Монте-Карло в статистической термодинамике. «Наука», Москва. 1977.

11. J. A. Barker. A quantum-statistical monte carlo method; path integrals with boundary conditions. The Journal of Chemical Physics, 70(6):2914—2918, 1979.

12. D. M. Ceperley. Path integral Monte Carlo methods for fermions, in: Simulation in condensed matter, physics and chemistry (K. Binder and G. Cicotti eds.). Bologna, Itally, 1996.

13. H. Kleinert. Path integrals in quantum mechanics, statistics, polymer physics, and financial markets. World Scientific, Singapore, 2002.

14. J. E. Hirsch, R. L. Sugar, D. J. Scalapino, and R. Blankenbecler. Monte carlo simulations of one-dimensional fermion systems. Phys. Rev. В, 26(9):5033—5055, Nov 1982.

15. H. De Raedt and A. Lagendijk. Monte carlo calculation of the thermodynamic properties of a quantum model: A one-dimensional fermion lattice model. Phys. Rev. Lett., 46(2):77-80, Jan 1981.

16. P. N. Vorontsov-Velyaminov, M. 0. Nesvit, and R. I. Gorbunov. Bead-fourier path-integral monte carlo method applied to systems of identical particles. Phys. Rev. E, 55(2):1979-1997, Feb 1997.

17. M. Takahashi and M. Imada. Monte carlo calculation of quantum systems. Journal of the Physical Society of Japan, 53(3):963-974, 1984.

18. M. Imada. On the monte carlo method for fermions in multi-dimensional systems. Journal of the Physical Society of Japan, 53(9):2861—2864, 1984.

19. M. Imada and M. Takahashi. Quantum monte carlo simulation of a two-dimensional electron system melting of wigner crystal. Journal of the Physical Society of Japan, 53(11):3770-3781, 1984.

20. R. W. Hall. The treatment of exchange in path integral simulations via an approximate pseudopotential. The Journal of Chemical Physics, 89(7):4212-4215, 1988.

21. R. W. Hall. The exchange potential in path integral studies: Analytical justification. The Journal of Chemical Physics, 91(3):1926—1931, 1989.

22. R. W. Hall and M. R. Prince. Development, justification, and use of a projection operator in path integral calculations in continuous space. The Journal of Chemical Physics, 95(8):5999-6004, 1991.

23. R. W. Hall. Formally exact path integral monte carlo calculations using approximate projection operators. The Journal of Chemical Physics, 97(9):6481-6484, 1992.

24. W. H. Newman and A. Kuki. Improved methods for path integral monte carlo integration in fermionic systems. The Journal of Chemical Physics, 96(2): 14091417, 1992.

25. D. M. Ceperley. Path-integral calculations of normal liquid he3. Phys. Rev. Lett., 69(2):331-334, Jul 1992.

26. D. M. Ceperley. Path integrals in the theory of condensed helium. Rev. Mod. Phys., 67(2):279-355, Apr 1995.

27. V. S. Filinov. Cluster expansion for ideal fermi systems in the 'fixed-node approximation'. Journal of Physics A: Mathematical and General, 34(8):1665-1677, 2001.

28. С. Н. Мак and R. Egger. A multilevel blocking approach to the sign problem in real-time quantum monte carlo simulations. The Journal of Chemical Physics, 110(1):12-14, 1999.

29. R. Egger, L. Miihlbacher, and С. H. Мак. Path-integral monte carlo simulations without the sign problem: Multilevel blocking approach for effective actions. Phys. Rev. E, 61(5):5961-5966, May 2000.

30. J. Carlson and M. H. Kalos. Mirror potentials and the fermion problem. Phys. Rev. C, 32(5):1735—1741, Nov 1985.

31. A. Korzeniowski, J. L. Fry, D. E. Orr, and N. G. Fazleev. Feynman-kac path-integral calculation of the ground-state energies of atoms. Phys. Rev. Lett., 69(6):893-896, Aug 1992.

32. S. Zhang. Finite-temperature monte carlo calculations for systems with fermions. Phys. Rev. Lett., 83(14):2777-2780, Oct 1999.

33. J. F. Corney and P. D. Drummond. Gaussian quantum monte carlo methods for fermions and bosons. Phys. Rev. Lett., 93(26):260401, Dec 2004.

34. J. F. Corney and P. D. Drummond. Gaussian phase-space representations for fermions. Physical Review В (Condensed Matter and Materials Physics), 73(12):125112, 2006.

35. V. Elser. Eliminating the minus sign in monte carlo simulations of fermions. Phys. Rev. A, 34(3):2293-2301, Sep 1986.

36. S. B. Fahy and D. R. Hamann. Positive-projection monte carlo simulation: A new variational approach to strongly interacting fermion systems. Phys. Rev. Lett., 65(27):3437-3440, Dec 1990.

37. G. G. Batrouni and P. de Forcrand. Fermion sign problem: Decoupling transformation and simulation algorithm. Phys. Rev. B, 48(l):589-592, Jul 1993.

38. J. Carlson, J. E. Gubernatis, G. Ortiz, and S. Zhang. Issues and observations on applications of the constrained-path monte carlo method to many-fermion systems. Phys. Rev. B, 59(20):12788-12798, May 1999.

39. E. L. Pollock and D. M. Ceperley. Simulation of quantum many-body systems by path-integral methods. Phys. Rev. B, 30(5):2555-2568, Sep 1984.

40. A. P. Lyubartsev and P. N. Vorontsov-Velyaminov. Path-integral monte carlo method in quantum statistics for a system of n identical fermions. Phys. Rev. A, 48(6):4075-4083, Dec 1993.

41. I. R. McDonald and K. Singer. Machine calculation of thermodynamic properties of a simple fluid at supercritical temperatures. The Journal of Chemical Physics, 47(11) :4766-4772, 1967.

42. J. Lee. New monte carlo algorithm: Entropic sampling. Phys. Rev. Lett., 71(2):211-214, Jul 1993.

43. F. Wang and D. P. Landau. Efficient, multiple-range random walk algorithm to calculate the density of states. Phys. Rev. Lett., 86(10):2050-2053, Mar 2001.

44. A. P. Lyubartsev, A. A. Martsinovski, S. V. Shevkunov, and P. N. Vorontsov-Velyaminov. New approach to monte carlo calculation of the free energy: Method of expanded ensembles. The Journal of Chemical Physics, 96(3): 1776-1783, 1992.

45. А. П. Любарце в, А. А. Марциновский, П. H. Воронцов-Вельяминов, Т. В. Кузнецова. Журнал физической химии, том 67:254, 1993.

46. M. F. Herman, E. J. Bruskin, and B. J. Berne. On path integral monte carlo simulations. The Journal of Chemical Physics, 76(10):5150—5155, 1982.

47. P. N. Vorontsov-Velyaminov, S. D. Ivanov, and R. I. Gorbunov. Quantum gas in an external field: Exact grand canonical expressions and numerical treatment. Phys. Rev. E, 59(1): 168-176, Jan 1999.

48. P. N. Vorontsov-Velyaminov and A. P. Lyubartsev. Entropic sampling in the path integral monte carlo method. Journal of Physics A: Mathematical and General, 36(3):685-693, 2003.

49. J. S. Kole and H. DeRaedt. Phys. Rev. E., 64:016704, 2001.

50. E. А. Поляков, П. Н.Воронцов-Вельяминов. Квантовый газ во внешнем поле при конечных температурах. Точное выражение для плотности и возбужденный состояния. Вычислительные методы и программирование (http://num-meth.srcc.msu.ru/), 8:334-351, 2007.

51. С. L. Pekeris. lsl, 2sl, and 2s3 states of h- and of he. Phys. Rev., 126(4):1470-1476, May 1962.

52. M. Sprik, M. L. Klein, and D. Chandler. Staging: A sampling technique for the monte carlo evaluation of path integrals. Phys. Rev. B, 31 (7)-.4234-4244, Apr 1985.

53. L. Brualla, K. Sakkos, J. Boronat, and J. Casulleras. Higher order and infinite trotter-number extrapolations in path integral monte carlo. The Journal of Chemical Physics, 121(2):636-643, 2004.

54. C. Predescu. Existence of short-time approximations of any polynomial order for the computation of density matrices by path integral methods. Phys. Rev. E, 69(5):056701, May 2004.

55. С. Predescu. Fast sampling algorithm for lie-trotter products. Phys. Rev. E, 71(4):045701, Apr 2005.

56. A. Bogojevic, A. Balaz, and A. Belie. Systematically accelerated convergence of path integrals. Phys. Rev. Lett, 94(18):180403, May 2005.

57. A. Bogojevic, A. Balaz, and A. Belie. Systematic speedup of path integrals of a generic n -fold discretized theory. Phys. Rev. В, 72(6)-.064302, Aug 2005.

58. S. A. Chin. Higher-order splitting algorithms for solving the nonlinear schro-umlaut]dinger equation and their instabilities. Physical Review E (Statistical, Nonlinear, and Soft Matter Physics), 76(5):056708, 2007.

59. A. Balaz, A. Bogojevic, I. Vidanovic, and A. Pelster. Recursive schro-umlaut]dinger equation approach to faster converging path integrals. Physical Review E (Statistical, Nonlinear, and Soft Matter Physics), 79(3):036701, 2009.