УДК 539.192
Ю. Б. Малыханов, С. В. Евсеев, И. Н. Еремкин
РАСЧЕТ ЭНЕРГИИ АТОМОВ В КОНФИГУРАЦИЯХ С ТРЕМЯ ОТКРЫТЫМИ ОБОЛОЧКАМИ В АЛГЕБРАИЧЕСКОМ ВАРИАНТЕ МЕТОДА ХАРТРИ - ФОКА
Аннотация. В рамках атомной теории Хартри - Фока - Рутана (метод Рутана -Багуса) выполнены высокоточные расчеты аналитических хартри-фоковских орбиталей и энергий атомов в возбужденных конфигурациях с тремя открытыми оболочками разной симметрии. Расчеты проводились в алгебраическом приближении с использованием в качестве базисных функций атомных орбиталей слэтеровского типа. Нелинейные параметры атомных орбиталей (орбитальные экспоненты) оптимизировались с помощью методов первого и второго порядка с высокой точностью. Благодаря этому удалось достичь точности выполнения вириального отношения 10-14-10-17, а вычисленные энергии термов атомов близки к хартри-фоковскому пределу.
Ключевые слова: атомная орбиталь, метод Хартри - Фока, энергия, атом, открытая оболочка.
Abstract. According to the nuclear theory of Roothaan - Hartree - Fock theory (a method of Roothaan - Bagus) the authors have conducted high-precision calculations of analytical Hartree - Focks orbitals and energies of atoms in excited configurations with three open shells of different symmetry. Calculations have been performed in algebraic approach applying nuclear orbitals of Slater type as basic functions. Non-linear parameters (orbital exhibitors) have been optimised by means of methods of the first and second order with split-hair accuracy. Due to this method it is possible to reach precision in performance of the virial relation 10-14-10-17, and the calculated energies of terms of atoms are close to Hartree - Fock limit.
Key words: atomic orbital, Hartree - Fock method, energy, atom, open shell.
Введение
Расчеты энергии и волновых функций атомов, необходимые для конкретных приложений, по-прежнему остаются актуальными, и интерес к таким расчетам со временем не ослабевает [1-7]. В случае тяжелых атомов подобные расчеты являются вполне выполнимыми на основе метода Хартри - Фока (ХФ). Уравнения ХФ удается решить численно практически для любых атомов [7, 8]. Ценность подобных расчетов заключается в нахождении с высокой точностью значений энергии термов, и эти результаты являются эталонными для аналитических (алгебраических) методов решения уравнений ХФ. Однако одноэлектронные функции, представляющие собой обширные таблицы, мало пригодны для хранения и последующего использования для конкретных приложений, если последние не были заложены в программу.
Существенно более удобным с точки зрения дальнейших приложений является алгебраический подход к решению уравнений Хартри - Фока, когда хартри-фоковские орбитали ищутся в аналитическом виде в форме разложения по функциям заданного вида - атомным орбиталям (АО). Выбор подходящего базиса АО как по размеру, так и по составу содержит большой произвол и зачастую требует многочисленных пробных расчетов. Именно благодаря этому произволу можно выбрать оптимальный базисный набор, с помо-
щью которого удается достичь хартри-фоковский предел рассчитываемой энергии и других свойств атомов.
1. Уравнения самосогласованного поля в методе Хартри - Фока - Рутана
В расчетах электронных оболочек атомов в алгебраическом приближении наибольшее распространение получила атомная теория Хартри - Фока -Рутана (ХФР) (метод Рутана - Багуса) [9], который применим для расчета атомов с любым числом открытых оболочек разной симметрии. Уравнения самосогласованного поля (ССП) в этом методе выводятся на основе выражения общего вида для энергии атома, в котором выполнено интегрирование по спин-угловым переменным. В связи с этим метод применим только для расчета энергии свободных атомов, которая не зависит от магнитного квантового числа электронов. Полученные в рамках этого метода аналитические харт-ри-фоковские орбитали можно использовать в расчетах электрической и динамической поляризуемостей, частот и сил осцилляторов электронных переходов [10-12].
Выражение для энергии (формула (12) из [9]) в алгебраическом приближении удобно записать с использованием полной матрицы плотности
БХ = Ос + °0 [13]
*=Г
-ЯР
(ы£бТ - бТ8абТ)Я(Ац) (Бц - БцБц (- N0)(нЦ- N0)
(1)
которое при варьировании автоматически приводит к единой матрице Фока. Здесь Ыс , Nо - число электронов заполненных и открытых оболочек типа X . В случае атомов в этой формуле индексы X и ц эквивалентны орбитальному квантовому числу I. Матрицы |р^Х,ц^(Б)| и |^Х,^(Б)| определе-
ны соотношениями
(
Р (М(б )} = Г Б,
1
А + Ц
Л
Г 4^
у=|А-|і|
рд.цг*
Г 2(а.і)< 1 А+ц
Г с^Ар^.цг* - — Г
у=0 у=|А-|і|
рд.цг*
(-)
в которых коэффициенты ./уХц и КУХц задают рассматриваемый терм, Ахц -численные коэффициенты заданного вида [9]. Кулоновский /Хрд цга и об-
менный интегралы вычислены с радиальными частями базисных
функций (АО).
Задача состоит в нахождении минимума энергии (1) как функции мно-
X X
гих переменных - элементов матриц плотности Бс и Бо и совокупности нелинейных параметров, от которых зависят АО (орбитальные экспоненты). Уравнения ССП можно получить из условия экстремума энергии 8Е = 0, где 8Е следует вычислять по формуле
8* = Г Яр Г §* С,
А дБТ і д^'
(3)
где С,і - нелинейные параметры АО (орбитальные экспоненты). В расчетах атомов наиболее эффективным является базис АО слэтеровского (экспоненциального) типа. а в расчетах молекул - базис гауссовского типа. Производные в формуле (3) легко найти исходя из (1) [13]:
дЕ
А
дБ
= ЕА= НА+ г|р (А’ц)БЦ)-
(нС- N0)
(нС- N0)
О (А.ц)( ) +
БаО(А,ц) ( ) + о(А,ц) ( )бА$
| - г£ - Г І«1'* —Г
-¿рБо (ууО(^ц) ( )) + 25рБхО(,ц) ( )бХ (§Х )/(Ых -Ыо) где V;А = ЭЛ/Э^I .
XXX
Принимая во внимание, что 8Б^ = 8Бс + 8Бо , учтем дополнительные
X X
условия, налагаемые на вариации 8Бс и 8Бо, которые следуют из свойств ортонормировки орбиталей, записанных в терминах матриц плотности:
бСSхбС = ыСбС, бОSхбО = N0Бо, бСSхбО = бОSхбС = 0, X = 0,1,2,...
Техника учета этих условий основана на разбиении матриц 8Бс и
5Бо на положительные и отрицательные проекции [13]:
8БС = (8БХ )- + (б£)+, 5Б° = (5Б0)(-) + (5Б&)(
где
(5Ба )- =
(5бА )8А БС + БС 8А (5бА ) —2- БС 8А (5бС )8А бС
N1
/ NC, (4.а)
(5БС)+ =-БС 58а БС / ^ =-Г БС ('і8А)бС5С
(4.б)
Проекции матрицы 8Бо вычисляются аналогично, если произвести за-
X XX X
мену Вс ^ £>о, Ыс ^ Ыо. С учетом этого формула (3) примет вид
5Е = Г|хр5БС (Ч + Яр5БА{¥>
+Г
і
Г™ А^А^А
Бс Г Бс Бо Г Бі
’С
А гДт^А \
V
NC
+ -
о1
'0
N0
(5)
где минус-проекция | вычисляется по формуле (4,а), в которой следует
X X
переставить местами матрицы Бс и 8 . Аналогично определяется минус-
*А
(X \ X
Е I . В выражении (5) вариации 8Бс, 8Бо и 8Сг- являются независимыми и из условия 8Е = 0 находим систему уравнений:
(А) = 0, (А)( ) = 0;
(6,а)
дСі N
Хр('і8 А)БС% Б —1ІХр('і8 АЫГ Ъ
¡с N0 ' '
Л = о
о = °.
(6,б)
Уравнения ССП (6,а) можно переписать в форме коммутационных соотношений:
ГА БС 8а = 8АБСГА, ГАБО8А = 8АбОГА.
Здесь Е - единая матрица Фока в методе Рутана - Багуса, с использованием которой уравнения ССП (6,а) для орбитальных коэффициентов заполненных и открытых оболочек в форме задачи на собственные значения имеют вид [13]
і7іА/~іА АоА/~іА
Г сі = гі 8 сі.
(7)
индекс 1 нумерует орбитали заполненных и открытых оболочек.
Уравнения (6,б) представляют собой необходимые условия минимума энергии как функции орбитальных экспонент, левая часть которого означает
полную производную (градиент).
й Су
Достаточным условием минимума энергии является положительная определенность матрицы вторых производных й Е / йСуйСj, выражение которой также можно найти. В конкретных расчетах вычисление точной второй
производной энергии по орбитальным экспонентам представляет собой очень сложную задачу, но ее можно решить более простым способом. С достаточной точностью вторую производную можно вычислить численным дифференцированием точной первой производной [13].
2. О решении уравнений Хартри - Фока - Рутана
Решение уравнений Хартри - Фока в алгебраическом приближении сводится к решению двух взаимосвязанных задач - непосредственно решению уравнений ССП для орбитальных коэффициентов, представляющих собой обобщенную задачу на собственные значения, и к нахождению оптимальных значений нелинейных параметров АО, обеспечивающих минимум рассчитываемой энергии атома. Уравнения ССП (7) традиционно решаются циклическим методом Рутана путем диагонализации матриц Фока методом Якоби. Условия сходимости циклического метода Рутана при решении уравнений ССП теоретически не доказаны и на практике сходимость существенно зависит от начальных значений орбитальных экспонент, которые выбираются достаточно произвольно. Кроме того, сходимость метода Рутана зависит также от выбора начальных значений для орбитальных коэффициентов занятых орбиталей, которые в нашей программе для удобства берутся равными нулю. Однако для некоторых атомов с таким начальным приближением циклический метод с очень высокой степенью согласования может не сходиться, и тогда используются матрицы плотности, полученные для предыдущего атома или другого терма этого атома. В программе в качестве критерия окончания циклического процесса используется равенство нулю (с точностью 10-11-10-13)
производной , которая представляет собой необходимое условие минимума энергии как функции матрицы плотности и означает выполнение уравнений (6,а) с заданной точностью.
Существенно более сложной задачей является высокоточная оптимизация орбитальных экспонент АО слэтеровского типа, которые использовались в наших расчетах. Традиционно для оптимизации орбитальных экспонент используются методы прямой минимизации (методы нулевого порядка), например метод сопряженных направлений Пауэлла [2, 14]. Как правило, эти методы эвристические по своей природе и имеют слабое математическое обоснование [15]. Важным качеством таких методов является слабая чувствительность к начальному приближению оптимизируемых параметров. Однако эти методы имеют достаточно медленную сходимость, дают низкую точность и требуют больших затрат машинного времени. Большой проблемой при реализации прямых методов является выбор критерия остановки счета, данная проблема в таких методах строго не может быть решена в принципе. Конкретные расчеты атомов, выполненные с использованием прямых методов, позволяют получить вириальное отношение с точностью всего 10-6-10-8 [1, 2, 14].
Другое дело - применение для оптимизации орбитальных экспонент методов минимизации первого и второго порядков, основанных на использо-
й dE d2 E В
вании выражений для ------ и ------- [13]. В таких алгоритмах имеется стро-
dCi d Cid С j
гий критерий окончания счета - необходимое условие минимума энергии
йЕ
----= 0 . Методы первого порядка, основанные на знании явного вида только
і
с1Е
, в случае нашей задачи имеют слабую сходимость (например, метод а Су
наискорейшего спуска), поэтому потребовался выход за рамки этих методов. На базе точных первых производных энергии (6,б) можно построить алгоритмы второго порядка, в которых матрица вторых производных энергии (матрица Гессе) заменяется на ее аппроксимацию. Такие методы называются квазиньютоновскими и относятся к методам переменной метрики. В наших расчетах [13] наиболее эффективным из них оказался метод Мурта-га - Саджента [16]. Скорость сходимости данного метода существенно выше, чем градиентных методов, хотя оба они используют только . Нако-
а С
нец, было показано, что особенно высокую скорость сходимости имеет метод Ньютона, даже если матрица Гессе вычисляется не по точной формуле, а путем численного дифференцирования . На каждом шаге итерационного
а С
процесса проверяется положительная определенность матрицы Гессе и в случае отсутствия таковой принимаются меры по восстановлению положительной определенности. В итоге была использована многоуровневая схема оптимизации орбитальных экспонент, позволяющая выполнить оптимизацию с любой наперед заданной точностью (в программе окончание счета происходит
при-----=10-12-10-14 а. е.):
а С1
- с помощью прямого метода (Нелдера и Мида [15]) достаточно произвольно выбранные орбитальные экспоненты оптимизируются с точностью
10-2-10-4;
- на полученных экспонентах запускается квазиньютоновский метод
йЕ
(Муртага - Саджента) и градиент------- доводится до значения 10-8-10-10;
а Су
- далее расчет ведется по методу Ньютона вплоть до достижения желаемой точности (в наших расчетах 10-12-10-14). Для проведения столь глубокой оптимизации Су требуется решение уравнений ССП (7) с аналогичной точностью. При этом вириальное отношение выполняется с высокой точностью (10-14-10-17). В ряде случаев первый пункт расчетной схемы можно пропустить, сразу запустив метод Муртага - Саджента. При использовании методов оптимизации нулевого порядка, например [1-6, 14, 17], вириальное отношение выполняется с точностью 10-6-10-8 или его значение авторами вообще умалчивается, что затрудняет сравнение результатов, полученных разными авторами.
3. Результаты и обсуждение
Применяя последовательно вышеизложенный алгоритм, мы получали такой базис, добавление в который АО любого типа уже не давало заметного улучшения энергии. Данный способ требует большого числа пробных расче-
тов, но получаемый базис является достаточно полным. Следует признать, что решение уравнений Хартри - Фока в алгебраическом приближении с высокой точностью для тяжелых атомов представляет собой достаточно сложную задачу. Полученные в результате оптимизационных расчетов энергии атомов с тремя открытыми оболочками приведены в табл. 1 и 2. Для сравнения в таблицы включены данные для энергии атомов, вычисленные другими авторами в результате численного решения уравнений Хартри - Фока [5]. Данные для атомов У-ИЬ с тремя открытыми оболочками 4d'ld5sl5pl -типа получены впервые нами.
Таблица 1
Сравнение рассчитанных (а. е.) в алгебраическом приближении значений энергии атомов 8с-№ с тремя открытыми оболочками 3dNd4s14p1 -типа с данными численного решения уравнений Хартри - Фока
Конфигурация атома Параметр Численное значение
21 8с [Лг^йЧіЧр1 (^ Т.Е К.Е Р.Е и/Т + 2 Т.Е“ 759,700622994082911 759,700622994082958 1519,40124598816587 2,1(-18) 759,700623
22 Ті [Лг^йЧіЧр1 (5в) Т.Е К.Е Р.Е и/Т + 2 Т.Еа 848,372272345411041 848,372272345411011 1696,74454469082205 -1,4(-17) 848,372273
23 V [Лфd34s14p1 (6в) Т.Е К.Е Р.Е и/Т + 2 Т.Еа 942,851303617011978 942,851303617012010 1885,70260723402399 3,3(-17) 942,851304
25 Мп [Лx]3d54s14p1 (8Р) Т.Е К.Е Р.Е и/Т + 2 Т.Еа 1149,82705689962773 1149,82705689962787 2299,65411379925560 -6,5(-17) 1149,827058
26 Бе [Лг^йЧіЧр1 (7Б) Т.Е К.Е Р.Е и/Т + 2 Т.Еа 1262,40361656886867 1262,40361656886850 2524,80723313773717 -2,2(-16) 1262,403618
27 Со [Лг^йЧіЧр1 (6Б) Т.Е К.Е Р.Е и/Т + 2 Т.Еа 1381,35622339477343 1381,35622339477328 2762,71244678954672 -1,1(-16) 1381,356226
28 Мі [Лг]3 йЧ^Чр1 (5в) Т.Е К.Е Р.Е и/Т + 2 Т.Еа 1506,79787195283829 1506,79787195286599 3013,59574390570428 1,8(-14) 1506,797875
Таблица 2
Рассчитанные (а. е.) в алгебраическом приближении энергии атомов У-ИЬ с тремя открытыми оболочками 4d'ld5sl5pl-типa
Конфигурация атома Параметр Численное значение
39 У [Кг]4й15і15р1 (4Б) Т.Е К.Е Р.Е и/Т + 2 3331,65920248460181 3331,65920248460217 6663,31840496920399 1,0(-16)
40 гг [Кг]4й25і15р1 (5в) Т.Е К.Е Р.Е и/Т + 2 3538,97424042552964 3538,97424042552859 7077,94848085105822 3,0(-16)
4ШЪ [Кг]4й35і15р1 (6в) Т.Е К.Е Р.Е и/Т + 2 3753,53370367425193 3753,53370367425040 7507,06740734850232 -4,1(-16)
43 Тс [Кг]4й55і15р1 (8Р) Т.Е К.Е Р.Е и/Т + 2 4204,76637000578629 4204,76637000578842 8409,53274001157471 -4,1(-16)
44 Яи [Кг]4й65і15р1 (7Б) Т.Е К.Е Р.Е и/Т + 2 4441,46540234801667 4441,46540234802064 8882,93080469603731 9,0(-16)
45 ЯИ [Кг]4й75і15р1 (6Б) Т.Е К.Е Р.Е и/Т + 2 4685,75878077752308 4685,75878077752114 9371,51756155504422 -4,2(-16)
46 ра [Кг]4й85і15р1 (5в) Т.Е К.Е Р.Е и/Т + 2 4937,72536682913386 4937,72536682913309 9875,45073365826695 -1,6(-16)
Примечание. Т.Е - полная энергия; К.Е - кинетическая энергия; - Р.Е - потенциальная энергия; и/Т + 2 - точность выполнения теоремы вириала; - Т.Еа - полная энергия, найденная путем численного решения уравнений Хартри - Фока [5]; [Аг] - 1^22522р63523р6; [Кг] - [Аг]4/4р6.
В плане сравнения результатов наших расчетов с результатами численного решения уравнений Хартри - Фока [5] можно сказать следующее: приведенные в [5] значения энергии с точностью 9-10 значащих цифр будут обоснованными, если теорема вириала выполняется с точностью 10-11-10-12, которую численные методы не дают [8]. Поэтому в приведенных в работе [5] значениях энергии последние 1-2 цифры являются результатом накопления ошибок округления. Выполненные в работе расчеты демонстрируют огромные возможности ньютоновских методов оптимизации в решении уравнений Хартри - Фока в алгебраическом приближении.
При проведении нами оптимизации базисных функций АО атомов с тремя открытыми оболочками было обращено внимание на то, что можно добиться практически линейной зависимости оптимизированных орбитальных экспонент однотипных АО от заряда ядра атомов в методе ХФР. Нали-
чие линейных зависимостей позволяет непосредственно получать волновые функции для любого атома с достаточной точностью без оптимизации орбитальных экспонент АО. В качестве примера на рис. 1 представлены линейные зависимости орбитальных экспонент (АО) 5-типа от заряда ядра Ъ атома. Такую же линейную зависимость имеют орбитальные экспоненты АО р- и й-типа. Полученные зависимости имеют не только теоретический интерес, но и очевидную практическую ценность.
Экспоненты АО s-функций
Рис. 1. Зависимость экспонент АО t>i 5-функций последовательности атомов Y-Rh (Z 39-46) с тремя открытыми оболочками 4^¥й5515/>1-типа от заряда ядра
Изложенный метод решения уравнений Хартри - Фока в алгебраическом приближении позволяет выполнить расчет любого атома периодической системы с высокой точностью с использованием персонального компьютера средней мощности. При этом никаких ограничений по числу открытых оболочек (разной симметрии) и размеру используемого базисного набора в программе не заложено.
Список литературы
1. Bunge, C. F. Roothaan - Hartree - Fock ground-state atomic wave functions / C. F. Bunge, J. A. Barrientos, A. V. Bunge // Phys. Rev. A. - 1992. - V. 46. - P. 36913696.
2. Koga, T. Roothaan - Hartree - Fock wave functions for atoms with Z<54 / T. Koga, H. Tatewaki, A. J. Thakkar // Phys. Rev. A. - 1993. - V. 47, № 5. - P. 4510-4512.
3. Tatewaki, H. Numerical Hartree - Fock energies of low-lying excited states of neutral atoms with Z < 18 / H. Tatewaki, T. Koga, Y. Sakai, A. J. Thakkar // J. Chem. Phys. -
1994. - V. 101, № 6. - P. 4945-4948.
4. Tatewaki, H. Numerical Hartree-Fock energies of low-lying excited states of neutral atoms with 19 < Z< 36 / H. Tatewaki, T. Koga // Chem. Phys. Lett. - 1994. - V. 228, № 6. - P. 562-567.
5. Tatewaki, H. Valence 4p functions for the first-row transition metal atoms / H. Tatewaki, T. Koga // J. Chem. Phys. - 2000. - V. 105, № 6. - P. 55-61.
6. Bunge, C. F. Slater-type orbital expansions and expectation values for Z = 2-54 I C. F. Bunge, J. A. Barrientos, A. V. Bunge II At. Data and Nucl. Data Tables. - 1993. -V. 53. - P. 113-162.
7. Koga, T. Chemically reliable uncontracted Gaussian-type basis sets for atoms H to Lr I T. Koga, H. Tatewaki, T. Shimazaki II Chem. Phys. Lett. - 2000. - V. 328. -P. 473-482.
8. Froese-Fisher, C. The Hartree-Fock Method for Atoms I C. Froese-Fisher. - New York, Wiley, 1977.
9. Roothaan, C. C. J. Atomic Self-Consistent Field Calculation by the Expansion Method I C. C. J. Roothaan, P. S. Bagus II Method in computational physics. - New-York : Academic Press. - 1963. - V. 2. - P. 47-94.
10. Малыханов, Ю. Б. Pасчет спектроскопических характеристик атомов группы железа в приближении Хартри - Фока - Pутана I Ю. Б. Малыханов, И. Н. Еремкин II Труды XXIV Съезда по спектроскопии. - Троицк : Тровант, 2010. - Т. 2. -С. 515-516.
11. Малыханов, Ю. Б. Pасчет оптической поляризуемости атомов с двумя открытыми оболочками в приближении Хартри - Фока I Ю. Б. Малыханов, И. Н. Еремкин II Журнал прикладной спектроскопии. - 2008. - Т. 75, № 4. -С. 458-462.
12. Малыханов, Ю. Б. Pасчет коэффициентов векторной связи для атомов с одной и двумя открытыми оболочками в приближении Хартри - Фока I Ю. Б. Малыханов, И. Н. Еремкин II Журнал прикладной спектроскопии. - 2011. - Т. 78, № 3. - С. 325-332.
13. Малыханов, Ю. Б. Pасчеты энергий возбужденных состояний атомов с открытой оболочкой в приближении Хартри - Фока I Ю. Б. Малыханов, С. А. Pома-нов II Журнал прикладной спектроскопии. - 2005. - Т. 72, № 3. - С. 292-299.
14. Koga, T. Improved Roothaan - Hatree - Fock wave functions for atoms and ions with N < 54 I T. Koga, S. Watanabe, K. Kanayama, R. Yasuda II J. Chem. Phys. -
1995. - V. 103, № 8. - P. 3000-3005.
15. Гуснин, С. Ю. Минимизация в инженерных расчетах на ЭВМ. Библиотека программ I С. Ю. Гуснин, Г. А. Омельянов, Г. В. Pезников и др. - М. : Машиностроение, 1981. - 120 с.
16. Murtagh, B. A. Variable metric gradient method for minimization procedures I
B. A. Murtagh, R. W. H. Sargent II Comput. J. - 1970. - V. 13, № 1. - P. 185-190.
17. Koga, T. Roothaan - Hartree - Fock wave functions for ions with Z < 54 I T. Koga, H. Seki, A. J. Thakkar and H. Tatewaki II J. Phys. B: At. Mol. Opt. Phys. - 1993. -V. 26. - P. 2529-2532.
Малыханов Юрий Борисович
доктор физико-математических наук, профессор, кафедра физики и методики обучения физике, Мордовский государственный педагогический институт им. М. Е. Евсевьева (г. Саранск)
E-mail: [email protected]
Евсеев Сергей Васильевич начальник редакционно-издательского центра, Мордовский государственный педагогический институт им. М. Е. Евсевьева (г. Саранск)
E-mail: [email protected]
Malykhanov Yury Borisovich Doctor of physical and mathematical sciences, professor, sub-department of physics and physics teaching methods, Mordovia State Pedagogical Institute named after M. E. Evsevyev (Saransk)
Evseev Sergey Vasylyevich Director of editorial and publishing center, Mordovia State Pedagogical Institute named after Evsevyev (Saransk)
Еремкин Игорь Николаевич Eryomkin Igor Nikolaevich
кандидат физико-математических наук, Candidate of physical and mathematical
инженер, ОАО «Орбита» (г. Саранск) sciences, engineer, “Orbita” plc. (Saransk)
E-mail: [email protected]
УДК 539.192 Малыханов, Ю. Б.
Расчет энергии атомов в конфигурациях с тремя открытыми оболочками в алгебраическом варианте метода Хартри - Фока I Ю. Б. Малыханов, С. В. Евсеев, И. H. Еремкин II Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2011. - № 3 (19). -
С. 120-130.