Аннотация: Целью работы является разработка компьютерной программы моделирования кластерной структуры расплавов простых металлов с кристаллической решеткой типа объемно-центрированной кубической (ОЦК) для сопоставления с экспериментальными данными рентгенодифракционного анализа. Для этого разработаны методы:
1. Расчета координационных чисел по кластеру для их последующего усреднения по образцу и сопоставления (аппроксимации) с экспериментальной радиальной функцией распределения атомов (РФРА). При этом выбор формы кластера осуществлялся в соответствии с принципом Кюри-Вульфа и Бравэ – о минимуме поверхностной энергии кристалла, находящегося в равновесии со своей жидкостью, и тем, что кристалл ограничивается атомными плоскостями с максимальной плотностью атомов. Для ОЦК решетки им соответствуют кластеры в форме: бипирамиды, призмы, трипирамиды, причем первой из названых соответствует наибольшее отношение объема к площади поверхности т.е. она есть наиболее вероятной.
2. Для указанных выше форм кластеров были определены дискретные функции формы и их непрерывные аналоги в различных направлениях трансляций.
3. По полученным функциям были определены Фурье – образы (интенсивности) для описания профилей дифракционных пиков для случая бипирамиды.
Использованием полученных зависимостей разработана компьютерная программа моделирования кластерной структуры, которая заключалась в аппроксимации экспериментальных данных рентгенодифракционных исследований расплавов щелочных металлов теоретической моделью. Результатами аппроксимации есть: определение среднего значения координационного числа, среднего межатомного расстояния, оптимального размера кластера, среднего расстояния между кластерами. Расчеты выполнены для многих щелочных металлов в широком интервале температур. Приведенные данные для рубидия при температуре 513К имеют следующие значения: среднее межатомное расстояние - 5,755 Å, оптимальный размер кластера - 12,756 Å, среднее координационное число - 1,325, среднее расстояние между кластерами - 0,518 Å.
Ключевые слова: координационные числа, оптимальный размер кластера.
Физико-математическиенауки
УДК 538.911
Дедоборец Александр Иосифович
Кандидат физико-математических наук, доцент,
доцент кафедры физики ДГАЭУ
Dedoborez O.
Candidate of physics and mathematics sciences,
associate professor,
associate professor the department of physics DSAEU
МОДЕЛИРОВАНИЕ ПРОЦЕССА КЛАСТЕРИЗАЦИИ В РАСПЛАВАХ МЕТАЛЛОВ С КРИСТАЛЛИЧЕСКОЙ РЕШЕТКОЙ ТИПА ОЦК
DESIGN OF PROCESSES OF CLUSTERIZATION IN FUSIONS OF METALS WITH THE CRYSTALLINE GRATE OF TYPE OF VCK
Аннотация: Целью работы является разработка компьютерной программы моделирования кластерной структуры расплавов простых металлов с кристаллической решеткой типа объемно-центрированной кубической (ОЦК) для сопоставления с экспериментальными данными рентгенодифракционного анализа. Для этого разработаны методы:
1. Расчета координационных чисел по кластеру для их последующего усреднения по образцу и сопоставления (аппроксимации) с экспериментальной радиальной функцией распределения атомов (РФРА). При этом выбор формы кластера осуществлялся в соответствии с принципом Кюри-Вульфа и Бравэ – о минимуме поверхностной энергии кристалла, находящегося в равновесии со своей жидкостью, и тем, что кристалл ограничивается атомными плоскостями с максимальной плотностью атомов. Для ОЦК решетки им соответствуют кластеры в форме: бипирамиды, призмы, трипирамиды, причем первой из названых соответствует наибольшее отношение объема к площади поверхности т.е. она есть наиболее вероятной.
2. Для указанных выше форм кластеров были определены дискретные функции формы и их непрерывные аналоги в различных направлениях трансляций.
3. По полученным функциям были определены Фурье – образы (интенсивности) для описания профилей дифракционных пиков для случая бипирамиды.
Использованием полученных зависимостей разработана компьютерная программа моделирования кластерной структуры, которая заключалась в аппроксимации экспериментальных данных рентгенодифракционных исследований расплавов щелочных металлов теоретической моделью. Результатами аппроксимации есть: определение среднего значения координационного числа, среднего межатомного расстояния, оптимального размера кластера, среднего расстояния между кластерами. Расчеты выполнены для многих щелочных металлов в широком интервале температур. Приведенные данные для рубидия при температуре 513К имеют следующие значения: среднее межатомное расстояние - 5,755 Å, оптимальный размер кластера - 12,756 Å, среднее координационное число - 1,325, среднее расстояние между кластерами - 0,518 Å.
Ключевые слова: координационные числа, оптимальный размер кластера.
Annotation: The purpose of work is development of the computer program of design of cluster structure of fusions of simple metals with the crystalline grate of type of body centered cubic (BCC) for comparison with experimental data of X – raying. For this purpose methods are worked out:
1. Calculation of coordination numbers on a cluster for their subsequent middling according to sample and comparison (approximations) with the experimental radial function of distribution of atoms (RFDA). Thus the choice of form of cluster was carried out in accordance with principle of Curie-Woulf and Brave - about a minimum of superficial energy of crystal, being in an equilibrium with the liquid, and that a crystal is limited to the atomic planes with the maximal closeness of atoms. For BCC of grate clusters correspond them in a form: bipyramids, prisms, threepyramids, moreover first from named corresponds most ratio of volume toward the area of surface, i.e. she is most credible.
2. For the forms of clusters indicated higher the discrete functions of form and their continuous analogues were certain in different directions of translations.
3. To the finding functions were certain Fourier- characters (intensities) for description of types of diffraction peaks for the case of bipyramids.
The use of finding dependences is work out the computer program of design of cluster structure, which consisted in approximation of experimental data of X - raying researches of fusions of alkaline metals a theoretical model.The results of approximation it is been: determination of mean value of coordination number, middle interatomic distance, optimal clustersize, middle distance between clusters. Calculations is executed for many alkaline metals in the wide interval of temperatures. Data for a rubidium at the temperature of 513К have next values: middle interatomic distance - 5,755 Å, optimal clustersize - 12,756 Å, middle coordination number - 1,325, middle c distance between clusters - 0,518 Å.
Key word: co-ordinating numbers, optimal clustersize.
Интерпретация данных дифракционных исследований на основе кластерной модели структуры жидкостей требует выполнения комплекса расчетных заданий - от выбора формы кластера к установлению связи между параметрами модели и данными дифракционного анализа. Соответствующие расчеты отображены в данной работе для кластеров со структурами объемноцентрированной кубической (ОЦК) решетки.
Выбор формы кластера определялся в соответствии с принципами Кюри-Вульфа о минимуме поверхностной энергии кристалла, который находится в равновесии со своей жидкостью и отвечающий принципу Браве, согласно которому кристалл ограничивается атомными плоскостями с максимальной плотностью атомов.
В ОЦК решетке такими являются атомные плоскости { 0, 0 , 0 , 101, 011, 10, 1 0, 10 , 01 , 0 1, 01} из которых могут быть построены кластеры в форме: бипирамиды, призмы, трипирамиды. Причем первая из названных отвечает наибольшим отношением объема к площади поверхности, то есть является наиболее вероятной.
При интерпретации данной радиальной функции распределения атомов (РФРА) в рамках кластерной модели необходимо определить методы расчета координационных чисел по кластеру для их следующего усреднения по образцу и сопоставления (аппроксимации) с экспериментальной РФРА. В отличие от бесконечного идеального кристалла координационное число атома, который отвечает к-ой координационной сфере, зависит от его положения, поэтому необходимо рассчитать среднее координационное число по кластеру для каждой координации. Очевидно, что это возможно следующим путем:
где -число атомов, которые имеют координационное число по каждой координационной сфере, - число атомов в кластере, в ребре которого υ- атомов.
В общем случае число атомов в кластере выражается кубической зависимостью от υ:
где α, β γ с определяется для каждой из форм кластеров. Числа:
также могут быть представлены аналогичной зависимостью
Но, если в первом случае, параметры определяются относительно просто, то для нахождения αк , βк , γк, δк необходимо определить значение Qk(υ) для разного числа атомов в ребре кластера (υ) по всем координационным сферам, которые использованы в расчетах. Для этого был создан алгоритм вычисления их с помощью ЭВМ, и после их определения поставленная задача могла быть решена для всех использованных форм кластеров.
Как известно, при рассмотрении дифракции на ограниченном объекте в структурном факторе целесообразно использовать функцию формы V(xp).
где xp - радиус-вектор, который соединяет некоторый произвольно выбран центральный атом из Р-ым атомом, S - фактор рассеяния. Функция формы
где - функция Эвальда;
где V - объем некоторой ограниченной области.
Были найдены дискретные формы для областей разной формы структуры ОЦК для кластеров в виде бипирамиды и их непрерывные аналоги. Также непрерывные функции формы найдены для всех форм кластеров структуры ОЦК.
Для бипирамиды:
( (четыре направления трансляции), где , а- большая диагональ бипирамиды.
(два направления трансляции), где .
Для обоих случаев 0 ≤ x ≤ L
Для них не составляет трудностей и нахождения Фурье- образов, для описания профилей дифракционных пиков (для случая бипирамиды)
В которых учитывалась и характерный размер кластера L по нормали к семейству отражающих плоскостей, а также, учтена кратность разных трансляций.
В первом случае трансляции
Во втором случае трансляции
Общая интенсивность –
где , d - межплоскостное расстояние.
Программа моделирования осуществлялась по такой схеме:
где τ- константа аподизации, ρ(х) - локальная атомная плотность. Причем
В правой части уравнения Gексп может быть найдена минимизацией квадратичной формы
Задавались исходные значения параметров zk, rk, σk, исходя из предположения о виде ближнего порядка соответственно положению дифракционного пика. σk - по скорости распространения звука и ближайшим межатомным расстоянием. Выходные zk определялись оценкой средних размеров упорядоченности в расплавах и по ее форме. Выбирая ту или другую форму области упорядоченности определяются такие значения ее средних размеров, чтобы отклонение теоретического профиля дифракционного пика от экспериментального было минимальным.
где модельная функция распределения кластеров по размерам
Коэффициент А определяется из условия нормировки
m=2-3 (наилучшей аппроксимации соответствует m=3, n=1).
Минимизацией (9) определяются rk, zk, σk – уточненные значения для 30 координационных сфер.
Разработанная компьютерная программа обеспечивала точность аппроксимации 0,05%. Для приведенного объекта [1] получено: среднее межатомное расстояние - 5,755 Å, оптимальный размер кластера - 12,756 Å, среднее координационное число - 1,325, среднее расстояние между кластерами - 0,518 Å.
Список литературы: