RUDN Journal of MlPh 2018 Vol. 26 No. 1 49-57
Вестник РУДН. Серия МИФ http://joumals.rudn.ru/miph
УДК 519.633, 538.9
DOI: 10.22363/2312-9735-2018-26-1-49-57
Об одном методе исследования самосогласованной нелинейной краевой задачи на собственные значения с растущими потенциалами И. В. Амирханов, Н. Р. Саркар
Лаборатория информационных технологий Объединённый институт ядерных исследований ул. Жолио-Кюри, д. 6, Дубна, Московская область, Россия, 141980
Один из распространённых методов исследования многочастичных задач в рамках вариационного подхода — переход к нелинейной одночастичной задаче путём введения самосогласованного поля, зависящего от состояний этих частиц. В работе рассматривается нелинейная краевая задача на собственные значения для уравнения Шрёдингера с растущим потенциалом, включающим зависимость от волновой функции и степенную зависимость от координаты г", п = 1, 2, 3 .. .. При п = 2 краевая задача для уравнения Шрёдингера (линейная задача) имеет точное решение. Для чётных степеней n показано, что решения такой задачи можно выразить через решения соответствующей линейной задачи, причём при п = 2 решение удаётся получить в явном виде. Получаемый для п = 2 набор решений характеризуется эквидистантностью расстояний между соседними собственными значениями. Показано, что решение нелинейной задачи отличается от решения линейной сдвигом собственных значений. В случае потенциала выше квадратичного, появляются новые растущие потенциалы меньшей степени. Для случая нечётных значений п обсуждается переход от интегро-дифференциальной формулировки задачи к системе дифференциальных уравнений, которая может быть решена численно на основе метода последовательных приближений, подтвердивший свою эффективность при исследовании модели полярона.
Ключевые слова: автолокализация, собственные значения, полярон, растущие потенциалы, нелинейная краевая задача
1. Введение
При исследовании решений многочастичных задач применяя вариационный принцип и путём введения самосогласованного поля, можно свести их к одночастичной задаче. Самосогласованные поля сами зависят от состояний этих частиц, и тем самым задача становиться нелинейной.
В качестве характерного примера можно привести концепцию полярона, являющуюся одной из классических моделей квантовой теории поля, которая имеет многочисленные приложения в физике конденсированных сред. Отметим, что проблема полярона была первоначально сформулирована как задача автолокализованного электрона в ионном кристалле. В настоящее время исследованию автолокализован-ных состояний в конденсированных средах посвящено много оригинальных работ, обзоров и монографий [1-5].
Эффект автолокализации в жидкостях приводит к образованию в них сольватиро-ванных электронов, играющих важную роль во многих химических процессах [6,7]. Под действием облучения вода (или иная среда) переходит в особое состояние, характеризуемое специальными физическими и химическими свойствами [8].
Подобные задачи возникают в рамках квантово-механической задачи двух тел, и её изучение является актуальной проблемой современной физики элементарных частиц. Например, в нерелятивистской потенциальной модели описание спектра тяжёлых кваркониев [9] сводится к решению уравнения Шрёдингера для двух тел.
Статья поступила в редакцию 9 января 2018 г.
Работа выполнена при финансовой поддержке гранта РФФИ, №17-01-00661^.
Авторы выражают благодарность Земляной Е.В., Дикусару Н.Д и Тухлиеву З.К. за полезное обсуждение и ценные замечания.
#
Релятивистское обобщение этой модели в рамках КХД, необходимое для единообразного описания спектров лёгких и тяжёлых кваркониев, приводит к релятивистским вариантам уравнения Шрёдингера. В качестве эффективного потенциала межквар-кового взаимодействия обычно используется сочетание растущего и кулоновского потенциалов. Численные исследования задач на собственные значения с кулоновским и линейно растущими потенциалами приведены в работах [10,11].
Другой пример использования самосогласованного поля — уравнение Хартри для исследования задач в ядерной или атомной физике, в том числе для системы электронов [12].
н
2
N
- 2т А + )+ £
t=ij=t
е2ф^)
^ = £j ^, 1 < j < N,
где потенциалы определяются из системы N уравнений Пуассона
АФ^ = —4^*^, 1 < i < N,
(1)
(2)
где Уяд(г^) — потенциал взаимодействия ] того электрона в кулоновском поле ядра.
Полученная таким способом потенциальная энергия (2) называется самосогласованным полем Хартри. Путём введения самосогласованного поля (2) многоэлектронная задача сводится к одноэлектронной, т.е. к решению уравнения Шрёдингера (1), содержащего координаты только одного электрона. Самосогласованное поле, правильно учитывающее симметрию перестановки частиц, было получено Фоком [12]. Система нелинейных интегро-дифференциальных уравнений Фока отличается от системы Хартри с присутствием обменного интеграла, обусловленным симметризацией координатных функций.
В данной работе предложен метод исследования решений нелинейной краевой задачи на собственные значения для уравнения Шрёдингера с растущим потенциалом, включающим зависимость от волновой функции.
где
2. Постановка задачи
Изучение стационарных автолокализованных состояний в среде, как правило, сводится к исследованию решений нелинейной краевой задачи в рамках уравнения Шрёдингера, которое в безразмерных переменных имеет вид
А + Е — сФ (М2) ф (г) = 0, (3)
Ф = JdflV (|г - п|) (fi)|2, (4)
при условии сохранения нормировки J df*("г)|2 = 1, а также граничные условия имеют вид:
ф(0) = const, ф (г) -> 0. (5)
Потенциал (4) является самосогласованным, т.е. так как зависит от самого решения уравнения (3).
Характер такой зависимости, а также форма функции V определяются конкретной моделью. Так, при исследовании состояний полярона V (|г — fi|) = 1/|г — fi|. Тогда
функция Ф удовлетворяет дифференциальному уравнению
ДФ = -4^|2. (6)
Отметим, что с вычислительной точки зрения обычно более удобно исследовать систему дифференциальных уравнений (3) и (6) с граничными условиями (5), а не систему интегро-дифференциальных уравнений (3)—(5). В работе проводится исследование нелинейной краевой задачи (3)—(5), где потенциал V имеет вид
V (| r-fi|) = |f-ri|n , п = 1, 2, 3,4,.... (7)
3. Методы исследования и анализ полученных результатов
Исследование свойств решений задачи (3)-(5), (7) будем проводить отдельно для чётных и нечётных значений.
1. Рассмотрим случай, когда п = 2, 4, 6,....
Сферически симметричное решение (3) ищем в виде ф (г) = (ф(г)/г)Уоо, с граничными условиями ф(0) = = 0, где Уоо = 1/ч/4л\ Тогда (4) примет следующий вид
сю к к
1 с- 32<
0 -к 0
При п = 2, после интегрирования по углам получим
Ф = 4^У dri^2(ri) f d(pi f d(91 Sin °i |г - ri Г . (8)
СС
Ф = A^2 + A2, где Ai = J dr i^2(r i), (9)
с учётом нормировки А\ = 1, А2 = ^ ёг ^ф2 (г 1).
о
Уравнение (3) для сферически симметричных решений в безразмерных переменных перепишем в виде
2
d
—2 + Ет - сФ d 2
'■Pm(r) = 0 m = 0,1, 2,..., (10)
с граничными условиями фто(0) = = 0.
Решая уравнение (10) с этими же граничными условиями, находим собственные значения Ет. Далее будем сравнивать аналитическое решения уравнения (10) с потенциалом Ф = г2 (линейная задача) с решениями нелинейной задачи с потенциалом (9). Решения ищем в виде
фг = Игфг, (11)
где ^¿(г = 0,1,2,...) — нормировочные постоянные, которые имеют вид
Щ = 1 , г = 0,1,2,...,
/
/ v?2dr
о
(ро = Г е-13г ; ( = г(г2 - а) е-13г ; ( = г(г2 - й2)(г 2 - аэ) е-^г ; (з = г(г2 - а_4)(г2 - а5)(г2 - аб)е-^г ; (р4 = т(г2 - аг)(г2 - ав)(г2 - ад)(г2 - а о) .
Подставляя выражения (11) для разных % = 0,1,2,... в уравнение (10) с потенциалом Ф = г2, получим систему алгебраических уравнений для нахождения собственных значений Ет(т = 0,1, 2, 3,4) и параметров а^ = 1, 2,..., 10. Решая эту систему уравнений, получим
Ео = 37, Е = 77, Е2 = 117, Ез = 157, Е4 = 197, 1 1 , ^ 1 , ^ 1
а = 27(3)' а2 = 2^(5 - аз = ^(5 + л/10), а4 = ^(14, 0657980),
а5 = —Ц(1, 3326518), аб = тг^(5, 6015501), а7 = —^(20, 3648752),
2ус 2ус 2ус
а8 = —^(10,2747750), ад = ^(4,3132752), а о = 7^(1 0470521), 3 = .
2у с 2у с 2 у с 2
Для нелинейной задачи с потенциалом Фто = г2 + А2т, где А2т = ^ ёгг2(^(г),
о
собственные значения Ет сдвигаются на сА2ТО. Тогда имеем Ео = Ео + СА20, Е = Е + СА21, ^2 = Е2 + СА22, Ез = Ез + СА23, Е4 = Е4 + СА24.
Таким образом, решение нелинейной краевой задачи с потенциалом (9) свели к решению линейной задачи с потенциалом Ф = г2. Действительно, если потенциал (9) поставить в уравнение (10), то получаем линейную задачу со сдвинутым собственным значением Ет = Ет + сА2ТО.
Так как уравнение (10) с потенциалом Ф = г2 имеет аналитическое решение, то и нелинейная задача также имеет аналитическое решение (см. табл. 1). Интегралы
а,™ = /*-
с функциями (11) вычисляются аналитически. Так как вычисления интегралов громоздкие, мы их опускаем, в таблице приведены только окончательные результаты. Из полученных результатов следует, что расстояния между соседними собственными значениями нелинейной задачи эквидистантны (см. последние колонки табл. 1). А также проверено, что решения (11) pn ортогональны, т.е.
сс
/ drptpm = 0, I = т.
При п = 4 после интегрирования по углам (8) получим
с
10 {
Ф = уЛ2Г2 + Лхг4 + A3, где А3 = drirfp2 (г 1). (12)
Таблица 1
Собственные значения линейной задачи Е, нелинейной задачи Е и значения интеграла А2т при п =2, т = 0,1, 2, 3,4
е т 0.25 1.0 4.0
Собственные значения линейной задачи Е 0 1.5 3 6
1 3.5 7 14
2 5.5 11 22
3 7.5 15 30
4 9,5 19 38
Значения интеграла 0 3 1,5 0,75
1 7 3,5 1,75
2 11 5,5 2,75
3 15 7,5 3,75
4 19 9,5 4,75
Собственные значения нелинейной задачи Е 0 2,25 4,50 9
1 5,25 10,5 21
2 8,25 16,5 33
3 11,25 22,5 45
4 14,25 28,5 57
Расстояние между соседними собственными значениями 3 6 12
3 6 12
3 6 12
3 6 12
При п =6 получим
сю
Ф = 7Лзг2 + 7Л2Г4 + Лхг6 + Л4, где Л4 = J ёгхг^^2 (г 1) (13)
0
и т.д.
В этом случае кроме потенциалов Ф = г4, Ф = г6 исходной задачи в нелинейной задаче появляются новые растущие потенциалы меньшей степени (см. (12) и (13)). Однако найти аналитические решения не удалось.
2. Теперь рассмотрим случаи, когда п принимает нечётные значения п = 1, 3, 5,... Пусть п = 1. Тогда
с
Фх(г) = У ёг1 |г — и | (п)|2 .
0
Эта функция обладает следующими свойствами:
с
Ф1(0)=/ёг1П (гх)|2 , Фх (г) -> г. (14)
0
Действуя на функцию Фх (г^ оператором А, получим
ДФх = 2Ф2, (15)
где
СЮ
Ф2 = / drïp^TT |Ф (ri)|2 . J |г - п|
о
Действуя на функцию Ф2 (г) оператором А, получим
ДФ2 = -4^|2, (16)
где
С
Ф2(0) = f dri - |ф (fi)|2 , Ф2 (г) -► -. (17)
J ri Г^Ю Г
0i
Таким образом, вместо системы интегро-дифференциальных уравнений (3)—(5) получили систему трёх дифференциальных уравнений (3), (15), (16).
Сферические симметричные решения системы (3), (15), (16) ищем в виде
р(г) V (г) V2 (г) ГГ
ф (г) = fSrLy,30, Ф1 (г) = vl^lYo0, Ф2 (г) = ^Y00, где Y0о = J-p.
г г г \ 4ж
Подставляя в систему уравнений (3), (15), (16), получим
p'' + Ер — —= Vip = 0, V/' = 2V2, V2' = — ^ (18)
V4n г г
с граничными условиями
p(0) = 0, р(то) = 0, Vi(0) = 0, Vi(0 -> г2, V2(0) = 0, V2'(to) = 0. (19)
При написании граничных условий (19) учтены свойства (14) и (17).
Для решения задачи (18), (19) можно использовать алгоритм на основе метода последовательных приближений, успешно апробированный для решения системы уравнений полярона в работе [13]:
1) решаем третье уравнение системы (18) при заданном (в начальном приближении) pk и найдём V2 ;
2) решаем второе уравнение для нахождения Vi ;
3) решая первое уравнение, найдём в следующем приближении pfc+i;
4) повторяем весь алгоритм до самосогласования, т.е. max \^pk+1(г) — pk(г)\ <
0<г<ю 1 1
£, где £ — заданная малая положительная величина, влияющая на точность получаемого численного решения.
Точность определяется порядком аппроксимации и шагом дискретной сетки.
Пусть п = 3. Тогда Ф3(г) = /С dri |г — Г[|3 |ф (ri)|2, где
СЮ
Фз(0) = I dr^ |Ф (ri)|2 , Фз (г)-► г3. (20)
I Г^-Ю
0
Действуя на функцию Фз (г) оператором А несколько раз, получим
ДФ3 = 12Ф1, АФ1 = 2Ф2, ДФ2 = —4^|ф|2. (21)
Таким образом, получили систему четырёх дифференциальных уравнений. Аналогично, при п = 5, получим систему пяти дифференциальных уравнений и т.д. Для численного решения этих систем необходимо обобщить вышеописанный алгоритм на случаю большего числа уравнений.
4. Выводы
Из полученных результатов можно сделать следующие выводы:
а) При исследовании обсуждаемой нелинейной краевой задачи (3)-(5), (7) в случае чётных значений п, эта задача сведена к решению линейной задачи. При п = 2 получено аналитическое решение нелинейной задачи, характеризующееся эквидистантностью расстояний между соседними собственными значениями спектра.
б) Для случая нечётных значений на основе исходной системы интегро-дифференциальных уравнений формулируется система дифференциальных уравнений, численное решение которой может быть проведено на основе метода последовательных приближений, что является предметом дальнейшей работы.
Литература
1. Пекар С. И. Исследования по электронной теории кристаллов. — М.: ГИТГЛ, 1951. — 256 с.
2. Каширина Н. И., Лахно В. Д. Математическое моделирование автолокализован-ных состояний в конденсированных средах. — М.: Физматлит, 2013. — 292 с.
3. Решение уравнений ЛЛП в теории биполярона / И. В. Амирханов, И. В. Пузы-нин, Т. А. Стриж, В. Д. Лахно // Известия АН, серия физическая. — 1995. — Т. 59, № 8. — С. 106-110.
4. Численное исследование квантово-полевой модели бинуклона сильной связи / И. В. Амирханов, Е. В. Земляная, В. Д. Лахно, И. В. Пузынин, Т. П. Пузынина, Т. А. Стриж // Препринт ОИЯИ, Дубна. — 1996. — № Р11-96-268.
5. Численное исследование нелинейной самосогласованной задачи на собственные значения в обобщенной модели полярона / И. В. Амирханов, В. Д. Лахно, И. В. Пузынин, Т. А. Стриж, В. К. Федянин // Препринт, Биологические исследования АН СССР, Пущино. — 1988. — 23 с.
6. Томпсон Д. Электроны в жидком аммиаке. — М.: Мир, 1979. — 138 с.
7. Численное исследование нелинейной самосогласованной задачи на собственные значения в обобщенной модели сольватированного электрона / И. В. Амирханов, И. В. Пузынин, Т. А. Стриж, О. В. Васильев, В. Д. Лахно // Препринт, Биологические исследования СССР, Пущино. — 1990. — 24 с.
8. Поляронная модель формирования состояний гидратированного электрона / В. Д. Лахно, А. В. Волохова, Е. В. Земляная, И. В. Амирханов, И. В. Пузынин, Т. П. Пузынина // Поверхность. Рентгеновские, синхротронные и нейтронные исследования. — 2015. — № 1. — С. 1-6.
9. Быков А. А., Дремин И. М., Леонидов А. В. Потенциальные модели кваркония // Успехи физических наук. — 1984. — Т. 143. — С. 3.
10. О некоторых проблемах численного исследования задачи на собственные значения в импульсном представлении / И. В. Амирханов, Е. В. Земляная, И. В. Пузынин, Т. П. Пузынина, Т. А. Стриж // Математическое моделирование. — 1997. — Т. 9, № 10. — С. 111-119.
11. Численное исследование релятивистских уравнений на связанные состояния с кулоновским и линейным потенциалами / И. В. Амирханов, Е. В. Земляная, И. В. Пузынин, Т. П. Пузынина, Т. А. Стриж // Математическое моделирование. — 2000. — Т. 12, № 12. — С. 79-96.
12. Поттер Д. Вычислительные методы в физике. — М.: Мир, 1975. — 387 с.
13. Амирханов И. В. и др. Численное исследование динамики поляронных состояний // Вестник тверского университета. Серия: Прикладная математика. — 2009. — № 17. — С. 5-14.
UDC 519.633, 538.9
DOI: 10.22363/2312-9735-2018-26-1-49-57
On a Method of Investigation of the Self-Consistent Nonlinear Boundary-Value Problem for Eigen-Values with Growing Potentials
I. V. Amirkhanov, N. R. Sarker
Laboratory of Information Technologies Joint Institute for Nuclear Research 6 Joliot-Curie Str., Dubna, Moscow region, Russia, 141980
One of the most common methods for investigating multiparticle problems in the framework of the variational approach is the transition to a nonlinear one-particle problem by introducing a self-consistent field that depends on the states of these particles. The paper considers a nonlinear boundary value eigenvalue problem for the Schrodinger equation with a growing potential including a dependence on the wave function and a power dependence on the coordinate V = rn where n =1, 2, 3.... For n = 2, the boundary value problem for the Schrodinger equation (linear problem) has an exact solution. For even powers of n, it is shown that solutions of such a problem can be expressed in terms of solutions corresponding to the linear problem, and for n = 2 the solution can be obtained in explicit form. The set of solutions obtained for n = 2 is characterized by equal distances between neighboring eigenvalues. It is shown that the solution of the nonlinear problem differs from the solution of the linear problem by the shift of the eigenvalues. In the case of a potential higher than the quadratic one, new growing potentials of a lesser degree appear. For the case of odd values of n, the transition is discussed, from the integro-differential formulation of the problem to a system of differential equations which can be solved numerically on the basis of the method of successive approximations, which has proved its effectiveness in the study of the polaron model.
Key words and phrases: self-localization, eigenvalues, polaron, growing potentials, nonlinear boundary value problem
References
1. S. I. Pekar, Studies on the Electronic Theory of Crystals, GITGL, Moscow, 1951, in Russian.
2. N. I. Kashirina, V. D. Lakhno, Mathematical Modeling of Autolocalized States in Condensed Media, Fizmatlit, Moscow, 2013, in Russian.
3. I. V. Amirkhanov, I. V. Puzynin, T. A. Strizh, V. D. Lakhno, Solution of LLP Equations in Bipolaron Theory, Bulletin of the Academy of Sciences, the series of physical 59 (8) (1995) 106-110, in Russian.
4. I. V. Amirkhanov, E. V. Zemlyanaya, V. D. Lakhno, I. V. Puzynin, T. P. Puzynina, T. A. Strizh, Numerical Investigation of the Quantum Field Model of the Strong-Binding Binucleon, Mathematical Modelling 8 (1997) 51-59, in Russian.
5. I. V. Amirkhanov, V. D. Lakhno, I. V. Puzynin, T. A. Strizh, V. K. Fedyanin, Numerical Study of a Nonlinear Self-Consistent Eigenvalue Problem in the Generalized Polaron Model, 1988, in Russian.
6. J. Thompson, Electrons in Liquid Ammonia, Mir, Moscow, 1979, in Russian.
7. I. V. Amirkhanov, I. V. Puzynin, T. A. Strizh, O. V. Vasilyev, V. D. Lakhno, Numerical Investigation of a Nonlinear Self-Consistent Eigenvalue Problem in the Generalized Model of a Solvated Electron, 1990, in Russian.
8. V. D. Lakhno, A. V. Volokhova, E. V. Zemlyanaya, I. V. Amirkhanov, I. V. Puzynin, T. P. Puzynina, Polaron Model of the Formation of Hydrated Electron States, Surface. X-ray, synchrotron and neutron studies (1) (2015) 1-6, in Russian.
9. A. A. Bykov, I. M. Dremin, A. V. Leonidov, Potential models of quarkonium, Successes of Physical Sciences 143 (1984) 3, in Russian.
10. I. V. Amirkhanov, E. V. Zemlyanaya, I. V. Puzynin, T. P. Puzynina, T. A. Strizh, On Some Problems of Numerical Investigation of the Eigenvalue Problem in the Momentum Representation, Mathematical modeling 9 (10) (1997) 111-119, in Russian.
11. I. V. Amirkhanov, E. V. Zemlyanaya, I. V. Puzynin, T. P. Puzynina, T. A. Strizh, Numerical Investigation of Relativistic Equations for Bound States with Coulomb and Linear Potentials, Math modeling 12 (12) (2000) 79-96, in Russian.
12. D. Potter, Computational Methods in Physics, Mir, Moscow, 1975, in Russian.
13. I. V. Amirkhanov, et al., Numerical Study of the Dynamics of Polaron States, Bulletin of Tver University. Series: Applied Mathematics (17) (2009) 5-14, in Russian.
© Амирханов И. В., Саркар Н.Р., 2018
Для цитирования:
Амирханов И. В., Саркар Н. Р. Об одном методе исследования самосогласованной нелинейной краевой задачи на собственные значения с растущими потенциалами // Вестник Российского университета дружбы народов. Серия: Математика. Информатика. Физика. — 2018. — Т. 26, № 1. — С. 49-57. — DOI: 10.22363/23129735-2018-26-1-49-57.
For citation:
Amirkhanov I. V., Sarker N. R. On a Method of Investigation of the Self-Consistent Nonlinear Boundary-Value Problem for Eigen-Values with Growing Potentials, RUDN Journal of Mathematics, Information Sciences and Physics 26 (1) (2018) 49-57. DOI: 10.22363/2312-9735-2018-26-1-49-57. In Russian.
Сведения об авторах:
Амирханов Илькизар Валиевич — старший научный сотрудник, кандидат физико-математических наук, начальник сектора Научного отдела вычислительной физики Лаборатории информационных технологий Объединённого института ядерных исследований, г. Дубна (e-mail: camir@jinr.ru, тел.: +7 (49621) 62547) Саркар Нил Ратан (Бангладеш) — кандидат физико-математических наук, старший научный сотрудник Научного отдела вычислительной физики Лаборатории информационных технологий Объединённого института ядерных исследований, г. Дубна (e-mail: sarker@jinr.ru, тел.: +7 (49621) 62547)
Information about the authors:
Amirkhanov I. V. — Senior Researcher, Candidate of Physical and Mathematical Sciences, Head of Sector "Scientific Division of Computational Physics". Laboratory of Information Technologies of the Joint Institute for Nuclear Research, Dubna (e-mail: camir@jinr.ru, phone: +7 (49621) 62547)
Sarker Nil Ratan — Candidate of Physical and Mathematical Sciences, Senior Researcher "Scientific Division of Computational Physics". Laboratory of Information Technologies of the Joint Institute for Nuclear Research, Dubna (e-mail: sarker@jinr.ru, phone: +7 (49621) 62547)