Журнал физической химии, 2020, T. 94, № 6, стр. 807-811

Оптимизация экспериментального определения тетраэдрации изобарно-изотермических фазовых диаграмм четырехкомпонентных систем

В. А. Шестаков a*, Е. В. Грачев a, В. И. Косяков a

a Российская академия наук, Сибирское отделение, Институт неорганической химии им. А.В. Николаева
Новосибирск, Россия

* E-mail: vsh@niic.nsc.ru

Поступила в редакцию 30.04.2019
После доработки 14.05.2019
Принята к публикации 15.10.2019

Полный текст (PDF)

Аннотация

Разработаны алгоритм и программы для оптимизации экспериментального исследования изобарно-изотермических фазовых диаграмм четырехкомпонентных систем с фазами постоянного состава или с небольшой областью гомогенности. Для этого разработаны: (а) программа перечисления фазовых диаграмм с фазами заданного состава, (б) алгоритм и программное обеспечение итерационного исследования фазовых диаграмм; на каждой итерации экспериментатор вводит в программу данные проведенного эксперимента и получает перечень составов, предпочтительных для дальнейшего исследования. Приведены результаты имитации построения фазовых диаграмм с одним бинарным, тремя тройными и одним четверным соединениям. Установлено, что из 10 приведенных примеров число итераций, достаточных для построения диаграммы, равно 1 (один пример), 2 (четыре примера), 3 (два примера), 4 (один пример), 5 (два примера).

Ключевые слова: четырехкомпонентные системы, фазовые диаграммы, топология фазовых диаграмм, экспериментальное исследование фазовых диаграмм

Основной объем информации о фазовых диаграммах в литературных источниках и специализированных базах данных относится к однокомпонентным, бинарным и тройным системам, достаточно хорошо разработана теория, описывающая геометрию и топологию таких диаграмм, а четырехкомпонентные системы изучены гораздо хуже (например, [18]). Существует несколько причин большого дефицита информации о четверных фазовых диаграммах. Основная из них связана с увеличением размерности. Симплекс составов четверных систем является тетраэдром, и геометрические образы фазовых диаграмм представляют собой фигуры в 3D-пространстве при P, Т = const, в 4D-пространстве при P = const или T = const и в 5D-пространстве при переменных Р и Т. Это приводит к усложнению традиционного описания диаграмм в виде двумерных изображений их проекций и сечений в литературных источниках, либо в виде трехмерных компьютерных образов. Кроме того, в реальных четырехкомпонентных системах присутствуют все фазы, которые содержатся в четырех ограняющих тройных системах, а также четырехкомпонентные фазы, внутри тетраэдра. Увеличение сложности химического строения четверных систем также приводит к дополнительным трудностям их описания. Другая причина дефицита информации о четверных фазовых диаграммах связана с тем, что для экспериментального определения их тетрангуляции требуется, как правило, проведение большого количества экспериментальной работы. Цель данной работы – построение алгоритма, позволяющего существенно снизить необходимое для этого количество экспериментов.

Пространство фазовой диаграммы является декартовым произведением подпространства составов на подпространство термодинамических параметров (одномерного Т) или (двумерного Р, Т) [9]. Подпространство параметров разделено на участки, в каждом из которых расположена своя изобарно-изотермическая фазовая диаграмма. Поэтому именно такие диаграммы – первостепенные объекты исследования. В настоящей работе рассматриваются системы, в которых присутствуют только твердые фазы постоянного состава. Изобарно-изотермическая диаграмма такой системы представляет собой концентрационный тетраэдр, разбитый на внутренние элементарные тетраэдры. Особенности строения таких конструкций (тетраэдраций) исследованы в работе [10]. Предложены принципы классификации тетраэдраций, описан алгоритм их построения. Разработанная программа позволяет построить исчерпывающее множество тетраэдраций с заданным количеством бинарных, тройных и четверных соединений и определить их топологические характеристики. Отметим, что в рассмотренной в [10] задаче перечисления тетраэдраций не принимается во внимание состав соединений.

ПОСТРОЕНИЕ ДИАГРАММ

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

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

1. Исходная информация

• Закодированный список химических соединений и их составы.

• Коды ограняющих систем (КОС). Обозначим компоненты четырехкомпонентной системы через k1, k2, k3, k4. Перечислим ограняющие трехкомпонентные системы (грани тетраэдра) k1, k2, k3; k1, k2, k4; k1, k3, k4; k2, k3, k4 и двухкомпонентные подсистемы (ребра тетраэдра) k1, k2; k1, k3; k1, k4; k2, k3; k2, k4; k3, k4. Количества трехкомпонентных соединений в соответствующих трехкомпонентных подсистемах и количества двухкомпонентных соединений в двойных подсистемах назовем кодом ограняющих систем (КОС). Таким образом он будет состоять из 10 натуральных чисел.

Принципы кодирования соединений и ограняющих систем описаны в работе [10].

2. Построение всех возможных тетраэдраций. По списку соединений и их составов программа вычисляет все тетраэдрации. Алгоритм решения этой задачи без учета состава соединений описан в [10]. Программа, использующая этот алгоритм, была адаптирована для перечисления тетраэдраций в системе с соединениями заданного состава. Модифицированная программа выдает список Р фигуративных точек, описывающих предпочтительные составы экспериментальных образцов (при этом наиболее информативные точки находятся в начале списка, в этих точках может реализоваться наибольшее число фазовых комплексов). В указанный список входят точки трех видов: (1) точки, лежащие на пересечении двух конод, (2) точки, лежащие на пересечении конод с треугольниками, образованными точками, отвечающими составам компонентов и соединений, (3) точки, расположенные в центрах всех возможных элементарных тетраэдров, образованных всеми конодами. Программа позволяет в процессе итерационной процедуры выбрать из списка фигуративных точек минимальное множество фигуративных точек, достаточное для построения любой тетрангуляции.

3. Эксперимент 1. Пользователь выбирает из множества Р произвольную точку p1. Программа предлагает список всех возможных двух-, трех- и четырехфазных равновесий в этой точке. Пользователь синтезирует образец и определяет его фазовый состав. Затем выбирает из списка равновесий то из них, которое удовлетворяет полученному экспериментальному фазовому составу. Выбранное равновесие вводится в программу.

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

5. Эксперименты 2, 3, … Пункты 3 и 4 данного алгоритма повторяются до тех пор, пока не будет построена единственная тетраэдрация, удовлетворяющая экспериментальным данным.

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

ПРИМЕРЫ

Как и в [10], будем обозначать бинарные соединения, расположенные на ребрах тетраэдра, символами сi, тройные, расположенные на гранях тетраэдра, – символами ti, четверные соединения – символами si. Строение диаграммы будем описывать с помощью двух графов [12]. Граф смежности вершин показывает проекцию трехмерного компьютерного изображения диаграммы на плоскость рисунка. Граф смежности тетраэдров показывает относительное расположение внутренних тетраэдров внутри концентрационного симплекса. Каждый внутренний тетраэдр помечен строкой из четырех символов компонентов и соединений, находящихся в его вершинах. Два тетраэдра являются смежными, если они имеют общую грань. Вершины графа смежности тетраэдров помечены обозначениями тетраэдров.

Для тестирования предлагаемого метода проведена имитация процесса построения фазовых диаграмм систем, содержащих пять соединений – с1, t1, t2, t3, s1. Координаты точек в пространстве составов будем выражать через мольные доли компонентов k2, k3, k4. Проведено 10 виртуальных экспериментов по построению фазовых диаграмм. В табл. 1 приведены коды ограняющих систем и составы соединений, выраженные через формулы соединений, где А = k1, B = k2, C = k3, D = k4, в рассмотренных фазовых диаграммах. Показано количество возможных тетраэдраций после каждой итерации, а также количество (n) экспериментально исследованных образцов в каждом примере.

Таблица 1.  

Тестовые эксперименты по построению тетраэдраций (n количество образцов)

КОС c1 t1 t2 t3 s1 Итерации n
1 3000000001 CD2 ABC2 A3BC3 A2B2C3 AB4CD3 38, 7, 4, 1 3
2 3000000100 BC3 ABC3 A3BC4 A2B2C3 AB4CD2 35, 6, 1 2
3 2100000001 СD A2B2C7 A2B2C AB2D6 ABCD4 300, 3, 1 2
4 2100000010 B2D2 ABC3 A3BC3 A2B2D7 AB4CD3 42, 24, 17, 5, 2, 1 5
5 2100000100 B4C AB3C9 ABC AB2D8 AB4CD2 272, 49, 5, 1 3
6 2100100000 AB2 A3B2C5 A2B2C3 A3BD3 AB4CD3 134, 103, 1 2
7 1110000001 CD4 A2B2C7 A2B2D AC2D6 A2BCD4 405, 1 1
8 1110001000 A3D4 A3BC3 A2BD4 A2C2D7 A2BCD4 299, 215, 11, 4, 2, 1 5
9 3000000001 C2D5 A3B2C12 AB3C9 A2B2C3 ABCD4 18, 2, 1 2
10 3000000100 B7C2 A2B2C3 A2B2C7 A3B2C5 AB4CD3 28, 7, 4, 2, 1 4

Опишем пример 1, сводные данные для которого приведены в табл. 2. Здесь N – номер итерации, Мт – число возможных тетраэдраций до данной итерации, Х – состав выбранной фигуративной точки, Е – равновесие, выбранное из всех возможных в данной точке, в соответствии с результатом виртуального эксперимента.

Таблица 2.  

Описание примера 1

N Мт Х Е
1 38 0.0200, 0.4000, 0.2000 k4, t1
2 7 0.2119, 0.3560, 0.0762 k1, k4, t2, t3
3 4 0.3999, 0.2001, 0.2999 k2, c1, t1, s1

Единственная оставшаяся тетраэдрация, представлена на рис. 1. Фазовая диаграмма является трехмерной конструкцией, имеющей 9 вершин, 29 ребер, 36 граней (из них 24 внутренних и 12 внешних), 15 тетраэдров. Строение диаграммы определено в результате трех итераций, т.е. для решения задачи при имитированном фазовом составе образцов оказалось достаточным приготовить три равновесных образца, с указанными в таблице составами.

Рис. 1.

Граф смежности вершин (а) и граф смежности тетраэдров (б) фазовой диаграммы, построенной в примере 1.

ЗАКЛЮЧЕНИЕ

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

Результаты тестирования продемонстрировали эффективность предложенного алгоритма исследований. Действительно, в примере 7 из 405 возможных тетраэдраций оказалось возможным выбрать единственную правильную тетраэдрацию на основе исследования одного образца. В остальных примерах потребовалось исследовать от 2 до 5 образцов. Отметим, что переход от одного из этапов к следующему происходит в результате произвольного выбора состава образца из предложенной таблицы, поэтому при повторной имитации исследования одной и той же диаграммы (с другим набором выбранных для расчета фигуративных точек) количество образцов может быть разным, но, скорее всего, небольшим.

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

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

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

  1. Курнаков Н.С. Введение в физико-химический анализ. Л.: ОНТИ; Химтеорет, 1936. 193 с.

  2. Аносов В.Я., Озерова Н.И., Фиалков Ю.Я. Основы физико-химического анализа. М.: Наука, 1976. 504 с.

  3. Халдояниди К.А. Фазовые диаграммы гетерогенных систем с трансформациями. Новосибирск: ИНХ СО РАН, 2004. 382 с.

  4. Райнз Ф. Диаграммы фазового равновесия в металлургии. М.: Металлургия, 1960. 376 с.

  5. Михеева В.И. Метод физико-химического анализа в неорганическом синтезе. М.: Наука, 1975. 272 с.

  6. MSI Eureka http://www.msiport.com/

  7. NIST Standard Reference Database 31 http://www.nist.gov/srd/nist31.cfm

  8. Киселева Н.Н. Компьютерное конструирование неорганических соединений: Использование баз данных и методов искусственного интеллекта. М.: Наука, 2005. 289 с.

  9. Краева А.Г. Методы изучения комбинаторных свойств диаграмм состояния // Сб. Вопросы геохимии и алгоритмы качественной теории фазовых превращений в многокомпонентных системах (диаграммы состояния). Институт геологии и геофизики СО АН СССР, Новосибирск, 1975. С. 10–47.

  10. Косяков В.И., Шестаков В.А., Грачев Е.В., Комаров В.Ю. // Журн. неорган. химии. 2016. Т. 61. № 10. С. 1325.

  11. Kosyakov V.I., Shestakov V.A., Grachev E.V. // MATCH Commun. Math. Comput. Chem. 2013. V. 69. № 3. P. 795.

  12. Lutsyk V., Vorobjeva V., Parfenova M. // Adv. Mater. Research. 2013. V. 704. P. 55.

Дополнительные материалы отсутствуют.