Читать статья по информационным технологиям: "Компьютерное моделирование комплексометрического титрования с учетом побочных реакций" Страница 2

(Назад) (Cкачать работу)

В программе "MODELCOM" (ОмГУ, 1996-1998) решение данной системы уравнений осуществляется с использованием метода покоординатного спуска на первом этапе и метода Ньютона-Канторовича для уточнения решения. Это позволяет обеспечить гарантированную сходимость на множестве допустимых входных параметров при достаточно простом алгоритмическом оформлении и высокой скорости вычислений: машины на базе процессора i80386DX-40 рассчитывают кривую титрования в течение нескольких секунд. Пользовательский интерфейс программы организован в соответствии с принципами "COMMON USER ACCESS", предусмотрены режимы тематической и контекстной помощи.

Входными параметрами являются: начальная концентрация металла (С0), концентрация титранта (СT) и его максимальный объем, начальный объем титруемого раствора (V0) и шаг ввода титранта. Учет побочных реакций требует ввода рН раствора и общей концентрации маскирующего вещества. Логарифмы констант учитываемых равновесий вводятся пользователем вручную либо запрашиваются в банке данных. Основной расчетной процедурой программы является построение интегральной и дифференциальных кривых титрования, то есть функций вида рМ = f(V), dpM/dV = f(V) и d2pM/dV2 = f(V). При этом учитывается разбавление раствора в ходе титрования. Предусмотрены расчет величины pM для любого момента титрования, исследование дифференциальной кривой на наличие экстремумов, получение данных по точкам эквивалентности, скачкам титрования и составу титруемых комплексов. Возможны одновременное построение на экране трех кривых титрования, отличающихся друг от друга по набору входных параметров, и выдача на печать результатов расчета и/или графиков.

На рис.1 показаны полученные компьютером кривые титрования ионов серебра с применением лигандов различной дентатности. Как и следует из теоретических соображений, на кривых проявляются как один (комплексон III), так и несколько (тиосульфат) скачков. ЭВМ прогнозирует, что при введении маскирующих веществ (например, аммиака) должен происходить сдвиг скачка, уменьшение его высоты, исчезновение одного или даже всех скачков (это зависит от относительной устойчивости разных MYi и MRi, а также величины избытка R). Естественно, исчезновение или сдвиг скачка указывают на неприемлемость соответствующих условий анализа.

Так как количественный расчет кривых титрования в подобных случаях упрошенными способами [1-3] провести нельзя, то для проверки правильности работы программы и адекватности выбранной нами модели сопоставляли некоторые кривые, рассчитанные с помощью программы Modelcom, с данными потенциометрического титрования тех же металлов в тех же условиях. Титрование вели при (202)0C, контролируя потенциал серебряного электрода на приборе рН-673 с погрешностью, меньшей 5 мВ, и усредняя данные по 3-5 параллельным опытам. Значения потенциалов переводили в значения pM по предварительно построенному градуировочному графику. Для серебряного электрода электродная функция была строго линейной в

Похожие работы