Компьютерное моделирование комплексометрического титрования смесей металлов Д.В. Калашников, В.И. Вершинин, Омский государственный университет, кафедра аналитической химии и химии нефти Постановка задачи и выбор модели.
    Моделирование комплексометрического титрования металлов позволяет в ходе вычислительного эксперимента провести априорную оптимизацию, что особенно важно при создании новых методик анализа. Известен ряд моделей и компьютерных программ на их основе [1-3], однако возможность титрования смесей металлов в подобных моделях обычно не предусматривалась(1) Именно этот случай имеет первостепенное значение для комплексометрии: методики, не требующие разделения определяемых металлов, как правило, менее трудоемки и дают более точные результаты [5]. Ранее мы использовали системы уравнений материального баланса для моделирования титрования одного металла произвольным лигандом в присутствии маскирующего вещества [6]. Данный подход позволяет обходиться без грубых упрощений (условных констант) и правильно моделировать кривые титрования, имеющие несколько скачков [7]. Целью настоящей работы было создание на этой же основе модели для комплексометрического анализа смесей металлов, в том числе в условиях маскирования.
    В ходе данной работы необходимо было: построить модель, приближенную к реальным условиям комплексометрического анализа; разработать алгоритм расчетов; составить реализующую его программу для IBM-совместимых компьютеров и подтвердить правильность ее работы сопоставлением результатов расчета с литературными данными, относящимися к комплексонометрии.
    По нашему мнению, необходимым и достаточным приближением можно считать следующую модель. Пусть при температуре 250С в водном растворе с нулевой ионной силой и произвольным значением рН присутствуют ионы нескольких (не более трех) металлов М1, М2, М3 и может присутствовать маскирующее вещество R, дающее с имеющимися металлами равновесную смесь комплексов. В ходе титрования раствор лиганда Y вводится отдельными небольшими порциями, причем после каждого прибавления достигается химическое равновесие. М1, М2 и М3 подвержены ступенчатому гидролизу, а R и Y ступенчато протонируются. Максимально возможное координационное число металлов - 6; дентатность лигандов может быть любой, основность их - не более 4. Наличие прочих компонентов раствора в модели не учитывается, т.е. требуется химическая инертность примесей. В системе не образуются смешанные, протонированные и полиядерные комплексы, не протекают окислительно-восстановительные реакции, не выпадают осадки. Начальные общие концентрации компонентов, рН раствора и константы равновесий известны. рН раствора и его температура в ходе титрования остаются постоянными.
    В наиболее сложном случае (присутствуют 3 металла и 1 маскирующее вещество) условие материального баланса приводит к системе из пяти уравнений, а именно:
    где СМ1, СМ2, СМ3, СY, СR - текущие общие концентрации M1, M2, M3, Y, R соответственно.
    После подстановки констант и алгебраических преобразований получается пригодная для численного решения система из пяти нелинейных уравнений с пятью неизвестными, которыми являются равновесные концентрации [M1], [M2], [M3], [Y] и [R]. Аналогичная система для более простого случая приведена нами в [6]. ............