Физическая химия
УДК 541.123.7 А.И. Сечной
ИСПОЛЬЗОВАНИЕ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ ДЛЯ РАЗБИЕНИЯ ПОЛИЭДРОВ СОСТАВОВ ФИЗИКО-ХИМИЧЕСКИХ СИСТЕМ
В рамках создания общего алгоритма исследования равновесного соотношения смесей фаз решена задача разбиения полиэдров составов физико-химических систем с любым типом взаимодействия компонентов с использованием линейного программирования. Определены характерные особенности разбиения полиэдров физико-химических систем с твердыми растворами, с «выклинивающимися» соединениями и с образованием более двух соединений на бинарных элементах огранения.
Изучение стабильного комплекса взаимных физико-химических (ФХС) представляет большой практический интерес, особенно в приложении к многокомпонентным ФХС, так как полное изучение диаграмм фазового равновесия (ДФР) часто затруднено по различным причинам. Между тем, изучение стабильного комплекса или его отдельных элементов может дать много информации о смесях, получаемых в результате взаимодействия солей, образующих данную ФХС. Можно проводить предварительное определение элементов стабильного комплекса, используя условные тепловые эффекты реакций обменного разложения взаимных пар солей [1].
После задания полиэдра составов ФХС в многомерном пространстве и определения возможных соединений [2] следующим этапом моделирования взаимодействия компонентов является разбиение полиэдра данной ФХС. Полиэдры составов являются геометрическим отображением реальных соотношений фаз в ФХС. Отражением факта существования определенного химического соединения является особая точка на диаграмме состав - свойство. Полюс определенного химического соединения (вершина полиэдра составов) является основанием для разбиения ДФР на связные части - симплексы [3].
Комплекс К называется связным, если его невозможно разбить в сумму непустых подкомплексов Ь и М без общих симплексов. Комплекс К связен только тогда, когда для каждых его двух вершин а и е существует последовательность вершин: а = аь а2,..., ач = е, причем любые две соседние вершины этой последовательности служат вершинами одномерного комплекса из К [4].
Идея разбиения полиэдра составов на симплексы прямо следует из определения компоненты комплекса. Компонентой некоторого комплекса К называется такой связный его подкомплекс Ь, что К распадается в сумму непересекающихся подкомплексов Ь и М. Если Кь ..., Кр -совокупность всех компонент комплекса К, то они попарно не пересекаются и в сумме составляют весь комплекс К [4]. В [5] было введено понятие триангуляции для описания разбиения треугольного комплекса трехкомпонентных ФХС на вторичные треугольники - симплексы. Секущим элементом в этом случае является квазибинарная ФХС, на которой располагается сингулярная складка. Приведены примеры деления тройных ФХС на симплексы, образованные несколькими бинарными соединениями. Определено, что число и расположение секущих линий подчиняется следующим правилам:
a) так как двойные и тройные соединения располагаются в вершинах вторичных треугольников, то точки пересечения сингулярных секущих друг с другом и сторонами первичного треугольника должны соответствовать стабильным и лабильным формам реально существующих соединений;
b) каждому вторичному треугольнику принадлежит взаимная трехлучевая звезда с центром в эвтектической или эвтонической точке;
c) сингулярные секущие получают свое начало во всех полюсах двойных и тройных соединений системы.
Сочетание вторичных треугольников дает плоские древа или карты химических диаграмм, которые представляют геометрические изображения основных реакций соединения, вытеснения (замещения) и обменного разложения. Число вторичных треугольников на единицу больше числа секущих линий. Примеры диаграмм взаимного вытеснения приведены в [5].
Введено понятие изомерии древ, обусловленное различным расположением центра сингулярной звезды. При увеличении числа соединений, увеличивается число секущих элементов и число вторичных треугольников. Была установлена связь между главными элементами, обусловливающими строение химической диаграммы: числом двойных С2 и тройных С3 соединений, числом сингулярных секущих линий д, вторичных треугольников к, тройных эвтектических точек Е, седловидных точек (Ван-Рейна) т: д = т - С2 + 3С3, к = Е = 1 + С2 + 2С3.
На основании вышеизложенного в приложении к сингулярной триангуляции были сделаны выводы [5]: сингулярная триангуляция дает геометрическое изображение реакций растворения, соединения, замещения и обмена в трехкомпонентной ФХС; при одинаковом числе, но различном расположении секущих элементов, наблюдается топологическая изомерия химических реакций; ДФР «состав - свойство» имеют одинаковое топологическое строение.
Политоп составов определен как многомерный аналог многоугольников (полигонов) на плоскости и многогранников (полиэдров) в трехмерном пространстве. Указано, что «-мерный политоп ограничивает часть «-мерного пространства своими (« - 1)-мерными гранями (элементы огранения). Смежные (« - 1)-мерные грани стыкуются в плоских (« - 2)-мерных гранях и т.д. «-мерный политоп составов разбивается на группу простейших политопов той же мерности, получивших название симплексы. «-мерным симплексом £(«+1) называют выпуклую замкнутую фигуру, определяемую (« + 1) независимо расположенными точками (вершинами) [1].
Стабильный комплекс взаимных ФХС определен как комплекс геометрических элементов, определяемый стабильными диагоналями в полиэдрах отдельных простейших взаимных ФХС огранения. Стабильным комплексом полиэдр составов разделяется на ряд стабильных клеток или ячеек, совокупность которых может быть рассматриваема как «древо» фигур [1].
Была сформулирована задача изучения стабильного комплекса взаимных ФХС как совокупности трех сторон проблемы, связанных между собой: 1) изучение геометрической структуры комплекса; 2) определение типов реакций обменного разложения, связанных с комплексом; 3) установление общих термохимических соотношений в системе [1].
В [6-10] отмечено, что особенно важно уметь предсказывать стабильные элементы ДФР «состав - свойство», тем самым и направление протекающих в системе реакций, на основании свойств входящих в нее компонентов. Определено направление реакции обмена во взаимной ФХС АХ + ВУ —— АУ + ВХ в связи со свойствами входящих в уравнение реакции атомов и ионов: массой и радиусами атомов, температурой плавления и т.д. На основании имеющихся экспериментальных данных показано, что наиболее пригодной закономерностью для определения направления реакции обмена во взаимных ФХС служит удельная энтальпия реакции АГИ, так как изменение энтропии АД конденсированных фаз, как правило, незначительно и величина АгО , в основном, определяется величиной АГИ реакции [3].
Рассмотрена топология трехкомпонентных взаимных ФХС с одним и двумя соединениями на боковых сторонах, а также с одним и двумя гетеросоединениями, расположенными на стабильной или лабильной диагонали. Сделано заключение о том, что иногда можно ограничиться изучением лишь сингулярного комплекса или его отдельных элементов, чтобы составить суждение о взаимодействии компонентов в ФХС [11-13].
Рассмотрена топология четырехкомпонентных взаимных ФХС с образованием одного двойного соединения на одном из вертикальных ребер призмы составов и на бинарной стороне одного из пинакоидов. Показано, что при наличии одного соединения на ребре, грани или внутри призмы имеется 17 основных видов тетраэдрации, причем призма составов может разбиваться на 4-8 вторичных тетраэдров, что говорит о сложных отношениях компонентов при образовании комплексов. Приведенные схемы древ дают возможность наглядно иллюстрировать разнообразие химического взаимодействия компонентов [6,7,10,13-16].
В [3] применена булева алгебра для описания триангуляции комплекса К« с множеством вершин Е = еь ..., ем. Введя высказывания: 1) Л^) - «подмножество вершин Y порождает «мерный симплекс триангуляции комплекса К«»; 2) Аг■ - «е1 не принадлежит Y», получили логическое выражение: Л(У)^с (А v А2 v-А ), 1 < 12 <... < 4, где х - конъюнкция, которая распространяется на все совокупности попарно смежных вершин е1 ,...,е- , таких, что натяну-
тый на них (к - 1)-мерный симплекс не принадлежит комплексу, тогда как его граница принадлежит данному комплексу КРаскрыв скобки, Л^) приводится к нормальной дизъюнктивной
т
форме v0l, в которой каждый член 0, имеет вид: Д с А с ... С А
,.1 о 1 к
т
Так как истинность высказывания 0, обеспечивает истинность всей дизъюнкции v0,, то
I °1
множество вершин Y, для которого высказывание еіо ї Х с ец ї Х с ... С еік ї Х истинно, порождает в К (Ы- к - 2)-мерный симплекс, причем Х = (Е\ {еіо ,еІ1е^ .
Поставив в соответствие дизъюнкции - сложение, конъюнкции - умножение, равнозначности - равенство, упорядочив множество вершин комплекса Кп, используя закон поглощения,
Р
получили минимальную дизъюнктивную форму v0,, для которой никакие два ее члена не по-
,°1
Р
глощаются один другим. Поскольку вся дизъюнкция v01 равнозначна Л^), то ее членам вза-
I °1
имно однозначно соответствуют всевозможные симплексы в Кп максимальной размерности. Выписав после этого для каждого слагаемого полученной суммы номера вершин, не входящие в него, получают все искомые симплексы. Показано, что минимальная дизъюнктивная форма содержит только слагаемые с числом символов, равным N - п - 1, а число слагаемых равно числу п-мерных симплексов комплекса К [3].
Несмотря на наличие большого числа способов, был выявлен ряд ФХС (с твердыми растворами, с «выклинивающимися» соединениями, с образованием более двух соединений на бинарных элементах огранения), для которых оказалось невозможным провести разбиение полиэдров составов. Суть предлагаемого метода, учитывающего все возможные типы ФХС, заключается в анализе числа сопряженных вершин полиэдра составов, являющихся геометрическим отображением физико-химического взаимодействия компонентов.
Блок-схема алгоритма разбиения полиэдров Х многокомпонентных ФХС представлена на рис. 1. Составляется матрица смежности (МС) компонентов, которая представляет собой компактную форму записи схемы сопряжений вершин полиэдра составов в виде квадратной таблицы, в которой по строкам и столбцам единицей («1») отмечены сопряженные компоненты и фазы исследуемой ФХС, а нулем («О») - несопряженные.
1. В первой строке МС выбирается компонент С1 (являющийся геометрическим аналогом фазы постоянного или переменного состава), соответствующий вершине 1 полиэдра.
1.1. Определяются сопряжения вершины 1 с другими вершинами полиэдра, в результате чего первая строка МС заполняется некоторой последовательностью из «1» и «о».
1.2. В первой строке МС отмечаются первая і и вторая] вершины, с которыми связана выбранная вершина 1, и определяется наличие или отсутствие связи между ними.
1.2.1. Если вершины і и] связаны между собой («1»), то делается вывод об образовании 3вершинного связанного графа {1, і, ]}, являющегося геометрическим отображением фазового симплекса С1 - Сі - С].
1.2.2. Если вершины не связаны («О»), то осуществляется переход к вершине ] + 1.
1.3. В случае реализации п. 1.2.1 проводится анализ размерности графа {1, і, ]} сравнением числа вершин, входящих в граф {1, і, ]}, с числом независимых компонентов ФХС.
1.3.1. Если получен набор связанных графов набора А2 , то проводится операция блока ^5 по поиску внутренних секущих (рис. 1).
1.3.2. Если получен набор связанных графов набора А3 , то проводится операция блока ^6 по выбору выявленных ранее внутренних секущих.
1.3.3. Если получен набор связанных графов набора А1+, то данный граф {1, і, ]} является геометрическим отображением одного из искомых фазовых симплексов.
1.4. Для следующих вершин 1, і и] + 1 повторяются операции п.п. 1.2.
1.4.1. Если между парой вершин, не входящих в {1, і, ]} и {1, і,7+1}, т.е. между] и] + 1, отсутствует связь, то осуществляется переход от вершины 7 + 1 к вершине 7 + 2.
1.4.2. Если вершины ] и ] + 1 связаны между собой, то делается вывод об образовании 4вершинного связанного графа {1, і, ], 7+1}, являющегося геометрическим отображением фазового симплекса С1 - Сі - С] - С7+1.
1.5. Повторяется операция п. 1.3.
А
В
Б
О
Н
1.6. Операции п.п. 1.1-1.5 повторяются до тех пор, пока не будут рассмотрены все вершины полиэдра, сопряженные с компонентом, соответствующим первой строке МС.
2. Во второй строке МС выбирается компонент С2, соответствующий вершине 2 полиэдра.
В последующем, рассматриваются все 1
вершины полиэдра, для которых повторяются операции п.п. 1.1-1.6.
В результате проведенной процедуры выявляются все возможные комбинации вершин полиэдра, т.е. генерируются наборы А1+ и/или А2+ равновершинных связанных графов, из которых формируется набор фазовых симплексов и «древо фаз» исследуемой ФХС [1,5].
Рассмотрим разбиение полиэдра составов четырехкомпонентной ФХС Ыа, К, Са // С1,
Мо04 [17] с использованием разработанного алгоритма. МС (табл. 1) составлена на основе схемы связи вершин полиэдра составов, представленного на рис. 2 а.
МС составлена в рациональной форме, в которой вершины полиэдра расположены в порядке увеличения числа связей (последний столбец МС). Для простейших ФХС оказывается возможным определить симплексы непосредственно из таких МС, составленных в рациональной форме. Однако, если в ФХС присутствуют внутренние секущие (появление графов набора А2+), это невозможно сделать без проведения дополнительной операции поиска и выбора реальных секущих. Поиск внутренних секущих осуществляется следующим образом. К каждому графу набора А2+ применяется логическая операция
К
2 (Д2і ^ Д2] ) \ ( Д2і ° Д2] )
і,=1 ;і * j
с учетом закона
поглощения [3]. Из полученного в результате набора 2-вершинных связанных графов А4+ выбираются те графы, которые не принадлежат элементам огранения, иными словами, проходят в объеме полиэдра составов. Полученные внутренние секущие вносятся в исходную МС, а затем повторяются все операции п.п. 1.1-1.6 алгоритма и далее.
М
N
О
Р и с. 1. Блок - схема алгоритма разбиения полиэдров составов МК ФХС методом последовательного попарного перебора вершин, соответствующих сопряженным компонентам
Т а б л и ц а 1
МС четырехкомпонентной ФХС Ма, К, Са // С1, Мо04
Вещество № вершины 2 5 8 4 9 1 3 7 6 £1
Ыа2Мо04 2 1 О 1 О О О О 1 1 4
СаС12 5 1 О О 1 1 О О 1 4
ЫаКМо04 8 1 1 О О О 1 1 5
К2Мо04 4 1 О О 1 1 1 5
КСаС13 9 1 1 1 О 1 5
ЫаС1 1 1 1 1 1 6
КС1 3 1 1 1 6
Ыа3С1Мо04 7 1 1 7
СаМо04 6 1 9
3 (
2 8> 7 V4
7 \4 3
а
46
7 6 7 6
б
Р и с. 2. Моделирование четырехкомпонентной взаимной ФХС Иа, К, Са // С1, Мо04: а) полиэдр составов; б) схема совмещения фазовых симплексов
Если в результате поиска на основе новой МС был получен набор А3+ равновершинных связанных графов, то это свидетельствует о том, что одна или несколько из введенных внутренних секущих пересекаются в объеме полиэдра составов. Это приводит к нарушению правила фаз, так как число сосуществующих фаз превышает мерность ФХС. В этом случае применяется операция выбора тех внутренних секущих, введение которых не приводит к появлению набора А3 . Учитываются вершины, которые входят в каждый симплекс набора А3 . Вершины, которые входят в каждый из этих симплексов, исключаются из них. Полученный остаток -графы набора А2+, добавляются теперь попарно к каждому из симплексов набора А3+, в результате чего получается только искомый набор А/ равновершинных связанных симплексов.
Схема поиска наборов связанных графов в четырехкомпонентной взаимной ФХС Иа, К, Са // С1, Мо04 методом последовательного перебора вершин полиэдра приведена в табл. 2.
Т а б л и ц а 2
Определение симплексов в системе Иа, К, Са // С1, Мо04
2 25 24 28 29 21 23 27 <78> 26 <68> П4 2678
278 268 <67>
5 58 54 59 51 <19> 53 57 56 <69> пз 1569
159 569 <16>
8 84 89 81 83 87 <47> 86 <46> П6 4678
478 468 <67>
4 49 41 43 47 <37> 46 <36> П5 3467
347 346 <67>
9 91 93 <13> 97 96 <36> П2 1369
139 369 <16>
1 13 17 <37> 16 <36> п 1367
137 136 <67>
В результате применения описанной выше операции поиска набора связанных графов А/ было получено шесть 4-вершинных связанных графов:
П {1, 3, 6, 7}, га {1, 3, 6, 9}, га {1, 5, 6, 9}, га {2, 6, 7, 8}, га {3, 4, 6, 7}, ^6 {4, 6, 7, 8}. Применение данного алгоритма не привело к формированию набора связанных графов А3+, что говорит об отсутствии в теле полиэдра дополнительных связей (внутренних секущих), приводящих к нарушению правила фаз.
Эти наборы связанных графов (симплексов) являются геометрическим отображением набора фазовых симплексов:
П = ЫаС! - КС! - СаМо04 - ЫаъаМо04; п = ЫаС! - КС! - СаМо04 - КСаС!ъ;
^3 = ЫаС! - СаСІ2 - СаМо04 - КСаС!3; М = Ш2Мо04 - СаМо04 - ЫаъС1Мо04 - ЫаКМо04,
^5 = КС! - К2Мо04 - СаМо04 - ШзС1Мо04; ^6 = К2Мо04 - СаМо04 - ШзС1Мо04 - ЫаКМо04.
Учитывая, что симплексы правильно примыкают друг к другу и определяют схему разбиения, можно составить «древо фаз» [5] или совмещение фазовых симплексов (рис. 2 Ь), составляющих фазовый комплекс ФХС Ыа, К, Са // С!, Мо04.
Если в ФХС образуются твердые растворы, то задача поиска набора фазовых симплексов решается аналогично с учетом кратности твердых растворов.
Следует отметить, что представленный алгоритм поиска стабильных ассоциаций фаз может решить поставленную задачу поиска набора фазовых симплексов и в случае, если по ряду ФХС огранения отсутствуют надежные экспериментальные данные. В данном случае в качестве искомых внутренних секущих будут выступать различные элементы ФХС огранения, и, применяя соответствующую операцию их поиска, можно восстановить картину взаимодействия компонентов в ФХС меньшей мерности.
Данные по стабильным ассоциациям фаз в многокомпонентных ФХС являются ценной информацией для использования ее при планировании экспериментального исследования отношения фаз, включающем: а) определение области составов, которая соответствует материалу с заданным набором физико-химических свойств; Ь) изучение физических свойств образцов;
с) выбор оптимальных условий синтеза материалов.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Радищев В.П. Многокомпонентные системы. М.: ИОНХ АН СССР. 1963. 502 с. Деп. в ВИНИТИ 30.01.57. № Т - 15616 - 63.
2. Сечной А.И. Моделирование взаимодействия в системе ЬіЕ - ЫаЕ - КТ - А1Е3. // Современные неорганические фториды. Сборник трудов I Международного сибирского семинара 18Ш-2003 по современным неорганическим фторидам. Новосибирск: ИНХ СО РАН - ОАО НЗХК, 2003. С. 240-243.
3. Краева А.Г. Определение комплексов триангуляции и-мерных полиэдров // Прикладная многомерная геометрия. Труды МАИ. № 187. 1969. С. 76-82.
4. ПонтрягинЛ.С. Основы комбинаторной топологии. М.: Наука, 1986. 120 с.
5. КурнаковН.С. Введение в физико-химический анализ. М.-Л.: Изд-во АН СССР, 1940. 563 с
6. Палкин А.П. Взаимосвязь и развитие тройных и четверных взаимных систем в расплавленном состоянии. Харьков: Харьковский ГУ, 1960. 338 с.
7. Бережной А.С. Многокомпонентные системы окислов. Киев: Наук. думка, 1970. 544 с.
8. Зедгенидзе И.Г. Планирование эксперимента при исследовании многокомпонентных систем. М.: Наука, 1976. 390 с.
9. Палатник Л.С., Ландау А.И. Фазовые равновесия в многокомпонентных системах. Харьков: Харьковский ГУ, 1961. 406 с.
10. Лупейко Т.Г. Анализ солевых систем. Ростов: Ростовский ГУ, 1981. 144 с.
11. ПетровД.А. Тройные системы. М.: Изд-во АН СССР, 1953. 315 с.
12. Захаров А.М. Диаграммы состояния двойных и тройных систем. М.: Металлургия, 1978. 296 с.
13. Захаров А.М. Многокомпонентные металлические системы с промежуточными фазами. М.: Металлургия, 1985. 133 с.
14. ВоловикБ.Е., ЗахаровМ.В. Тройные и четверные системы. М.: Металлургиздат, 1948. 228 с.
15. Захаров А.М. Диаграммы состояния четверных систем. М.: Металлургия, 1966. 240 с.
16. ПосыпайкоВ.И. Методы исследования многокомпонентных солевых систем. М.: Наука, 1978. 255 с.
17. Сечной А.И., Гаркушин И.К., Трунин А.С. Дифференциация четырехкомпонентной взаимной системы из шести солей Ыа, К, Са // С!, Мо04 и схема описания химического взаимодействия // Журн. неорган. химии. 1988, Т. 33. № 3. С. 752-755.
Поступила 8.04.2003 г.