1Назначение
ElectroChemLab — программный комплекс для количественного описания растворов электролитов: он рассчитывает равновесный состав сложной многокомпонентной системы, извлекает термодинамические и транспортные параметры из потенциометрического и кондуктометрического эксперимента и содержит полноценный раздел математической обработки данных, работающий и сам по себе.
Комплекс отвечает на два встречных вопроса. Прямая задача: дан состав раствора и набор констант — каким будет равновесный состав, ионная сила, коэффициенты активности, ЭДС ячейки, удельная электропроводность, вид кривой титрования. Обратная задача: дан измеренный сигнал — каковы константы устойчивости, произведения растворимости, предельные подвижности ионов, стандартный потенциал, параметры модели активности; и с какой достоверностью они определены по этим данным.
Расчётные методы комплекса не написаны заново «по учебнику»: в их основе — рабочие модули, годами применявшиеся в реальной практике обработки данных, перенесённые сюда с восстановлением исходных формул по первоисточникам. Найденные при переносе расхождения и опечатки — как в рабочих модулях, так и в самих публикациях — исправлены, а каждое исправление описано в справке соответствующего раздела и закреплено проверкой (§9).
1.1 Научная область и решаемые задачи
Предметная область — химия растворов и электрохимия на стыке с аналитической и физической химией: ионные равновесия (комплексообразование, кислотно-основные реакции, осаждение, ассоциация), термодинамика растворов электролитов и явления переноса в них. Величины, которые комплекс определяет из эксперимента:
- константы устойчивости комплексов $\lg K$ (ступенчатые и общие $\lg\beta$);
- константы кислотности $\mathrm{p}K_a$ и произведения растворимости $\mathrm{p}K_s$;
- константы ассоциации ионных пар — как обычные равновесия, а не как отдельная теория;
- стандартный потенциал ячейки $E^{\circ}$ и параметры активности — расстояние наибольшего сближения $\mathring{a}$ и коэффициент $b$;
- предельные молярные электропроводности ионов $\lambda^{0}$ и их радиусы;
- точка эквивалентности и концентрация определяемого вещества — из кривой титрования, без всякой модели равновесий.
Отдельно стоит раздел математической обработки: он не привязан к электрохимии и применим в физико-химическом эксперименте любого направления — кинетическом, спектрофотометрическом, хроматографическом, термохимическом (см. §5).
1.2 Кому адресован
Комплекс рассчитан на специалистов-химиков любого уровня подготовки, и уровень задаёт не сложность интерфейса, а глубина, на которую пользователь захочет спуститься.
| Уровень | Типичная задача | Что даёт комплекс |
|---|---|---|
| Студент практикум, курсовая, диплом |
Рассчитать состав раствора, построить кривую титрования, найти точку эквивалентности и концентрацию, увидеть, чем сглаживание отличается от аппроксимации | Пошаговый ввод «шаг N из 5», готовые примеры в комплекте, справка с выводом формул и разбором каждого метода |
| Аспирант, научный сотрудник исследование системы |
Определить константы по собственным измерениям, проверить гипотезу о составе комплексов, оценить достоверность параметров | Обратная задача по всей кривой, 11 методов поиска решения, доверительные интервалы, диагностика неидентифицируемости, сравнение моделей по AICc/BIC/PRESS |
| Аналитик заводской или испытательной лаборатории |
Рутинное определение концентрации по титрованию с оценкой погрешности | Два автономных окна анализа, которым не нужна ни модель равновесий, ни модель проводимости: только измеренная кривая — и результат с $\sigma$ |
| Преподаватель | Показать явление, а не результат | Режим моделирования: кривая строится по заданным константам, и видно, как её форма отзывается на каждый параметр |
2Как устроен комплекс
Комплекс поставляется как оболочка ElectroChemLab, в которой все задачи открываются из одного меню, и как четыре самостоятельные программы — на случай, когда нужна ровно одна из них. Это не разные сборки одного и того же: оболочка встраивает в свои вкладки те же самые окна, что запускаются отдельно, поэтому поведение и результат совпадают в точности.
2.1 Программы и модули
| Модуль | Что делает | Файл |
|---|---|---|
| Равновесия | Равновесный состав многокомпонентного раствора: концентрации всех частиц, ионная сила, коэффициенты активности, баланс заряда — для серии растворов сразу | *.eqp.xml |
| Потенциометрия | Обработка ЭДС по уравнению Нернста поверх решателя: поиск $E^{\circ}$, констант и параметров активности по серии растворов | *.ptp.xml |
| Потенциометрическое титрование ▸ Моделирование и обработка |
Кривая титрования целиком: смешение объёмов, равновесие в каждой точке, подгонка констант по всей кривой; осаждение — по той же схеме | *.ptp.xml |
| Потенциометрическое титрование ▸ Анализ |
Точка эквивалентности прямо из измеренной кривой: производные, пересечение ветвей $\Delta V/\Delta E$, четыре функции Грана — без модели равновесий | *.pta.xml |
| Кондуктометрия | Удельная и эквивалентная электропроводность: равновесный состав даёт свободные ионы, модель переноса — проводимость; поиск $\lambda^{0}$, радиусов и констант | *.cnd.xml |
| Кондуктометрическое титрование ▸ Моделирование и оптимизация |
Кондуктограмма поверх спесиации: показание прибора → κ, подгонка по всей кривой | *.cnt.xml |
| Кондуктометрическое титрование ▸ Анализ |
Концентрация по излому: пересечение линейных ветвей с переносом погрешности — нужна только измеренная кондуктограмма | *.cna.xml |
| Приготовление буферных растворов | Объёмы наличных растворов под заданный pH вместо справочной рецептуры, рассчитанной на круглые концентрации: пересчёт точен, потому что рецептура фиксирует количества вещества, а не объёмы. Справочник из девятнадцати систем; подбор ведётся по наименьшему объёму — им определяется погрешность рецепта | *.buf.xml |
| Математическая обработка восемь подпунктов |
Сортировка, сглаживание, интерполяция, сплайны, аппроксимация, 2D аппроксимация, дифференцирование, интегрирование | *.srt.xml … *.itg.xmlпо одному типу на подпункт |
2.2 Что у них общего
Модули не повторяют друг друга — они пользуются одними и теми же частями. Ввод исходных данных (растворитель, частицы, матрица стехиометрии, концентрации) — это один экранный модуль на все три химические программы; решатель равновесий — один; набор из одиннадцати методов поиска решения вместе с расчётом доверительных интервалов — один; окна математической обработки в оболочке и в отдельной программе — тоже одни и те же, а не копии.
Единообразен и сам интерфейс, и это заметно сокращает время освоения: ввод исходных данных везде идёт по шагам с рейлом слева, таблица во всех окнах одна и та же (вставка блока из Excel по Ctrl+V с ростом таблицы под размер вставленного, отмена последнего действия, пять форматов показа чисел), графики устроены одинаково, а у каждого окна есть свой пункт в меню «Справка». Освоив один модуль, пользователь уже умеет работать с остальными.
2.3 Комплекс в числах
Сводка по всему, что в комплексе посчитано и описано, — с указанием, где искать подробности.
| Группа | Сколько | Что именно | Подробно |
|---|---|---|---|
| Химическая часть | |||
| Формализм равновесий | 1 | Матричное описание Бринкли — на комплексообразование, кислотно-основные равновесия, ассоциацию и осаждение сразу | §3 |
| Модели коэффициентов активности | 4 | Предельный закон, уравнение Дэвиса, расширенное уравнение Дебая — Хюккеля и MSA; модель выбирается пользователем, а её параметры можно подбирать по эксперименту | §3.4 |
| Способы оценить достоверность решения | 3 | Доверительные интервалы, корреляционная матрица параметров и сравнение моделей активности по качеству подгонки | §4.5 |
| Модели электропроводности | 12 | 8 классических для симметричного бинарного электролита и 4 для смеси любого валентного типа | §4.2 |
| Методы поиска решения | 11 | 10 локальных из разных семейств и 1 глобальный; плюс доуточнение и расчёт доверительных интервалов | §4.1 |
| Оценки точки эквивалентности | 5 | Вершина параболы $dE/dV$, пересечение ветвей $\Delta V/\Delta E$, две функции Грана и вторая производная (как иллюстрация) | §4.4 |
| Функции Грана | 4 | Первый метод Грана и три функции второго метода (сильная кислота, слабая кислота, избыток титранта) | §4.4 |
| Математическая часть | |||
| Численные методы обработки | 62 | 13 + 9 + 8 + 8 + 5 + 9 + 10 по семи вычислительным подпунктам, плюс сортировка в двух режимах | §5.2 |
| Критерии выбора модели | 5 | Скорректированный $R^2$, AICc, BIC, GCV, PRESS — плюс число смен знака невязки, не зависящее от числа параметров | §5.4 |
| Документация | |||
| Страниц справки | 7 | По одной на каждый метод комплекса плюс эта, общая | §8 |
| Ссылок на литературу | 226 | Оригинальные источники каждого метода — с пометкой, где именно в программе использован каждый | §8 |
| Автоматических тестов | 1117 | Выполняются при каждой сборке | §9 |
3Матричный формализм: единая запись любых равновесий
Это основа всей химической части комплекса и главная причина, по которой разные на вид задачи решаются одной программой. Комплекс не различает «расчёт комплексообразования», «расчёт pH», «расчёт растворимости» как разные виды расчёта: он различает строки одной таблицы.
3.1 Базис и матрица стехиометрии
Из всех частиц раствора выбирается базис — минимальный набор, через который линейно выражаются остальные. Каждая производная частица $i$ задаётся строкой стехиометрических коэффициентов $\nu_{ij}$ по базисным частицам $j$ и своей константой образования $\lg K_i$. Закон действующих масс тогда записывается для всех частиц единообразно:
Условия материального баланса замыкают систему: аналитическая концентрация каждой базисной составляющей равна сумме её вкладов во все частицы,
Вся химия системы, таким образом, — это матрица $\nu$, столбец $\lg K$ и столбец зарядов. Вот как в одной такой таблице уживаются вещи, о которых учебник говорит в разных главах:
| Частица | Ag+ | Cl− | NH3 | H+ | z | lg K |
|---|---|---|---|---|---|---|
| AgCl2− — комплекс | 1 | 2 | 0 | 0 | −1 | 5.25 |
| Ag(NH3)2+ — комплекс | 1 | 0 | 2 | 0 | +1 | 7.23 |
| NH4+ — протонирование | 0 | 0 | 1 | 1 | +1 | 9.25 |
| OH− — автопротолиз воды | 0 | 0 | 0 | −1 | −1 | −14.00 |
| AgCl(тв) — осадок | 1 | 1 | 0 | 0 | 0 | −9.75 |
Одна таблица: комплексы, кислотно-основное равновесие, ионное произведение воды и осадок. Осадок отличается от прочих строк только пометкой «(тв)» в названии, а его константа — это $\lg K = -\mathrm{p}K_s$.
3.2 Как решается система
Система решается по формализму Бринкли методом Ньютона относительно логарифмов активностей базисных частиц. Сходимость двухуровневая, и это существенно: внутренний цикл Ньютона работает при замороженных коэффициентах активности (с дроблением шага по норме невязки), а внешний пересчитывает ионную силу и коэффициенты активности по уже найденному составу. Коэффициенты включаются только после того, как сошёлся первый, идеальный проход, — так система не разваливается на старте.
Коэффициенты активности рассчитываются по выбранной пользователем модели — их в комплексе четыре (§3.4). Значение по умолчанию, оно же историческое, — расширенное уравнение Дебая — Хюккеля:
причём $A$ и $B$ вычисляются из температуры, диэлектрической проницаемости и плотности растворителя, а $\mathring{a}$ и $b$ могут быть подобраны по эксперименту наравне с константами. Критерий сходимости внешнего цикла подстраивается под вид модели: там, где коэффициенты зависят от концентрации каждого сорта частиц по отдельности, а не только от ионной силы, к нему добавляется норма изменения самих коэффициентов — иначе ионная сила могла бы сойтись при ещё не сошедшихся $\gamma$. Относительная сходимость — $10^{-7}$ по составу и $10^{-5}$ по ионной силе; предусмотрены ограничение аргумента экспоненты и порог вырожденности, отнесённый к масштабу матрицы, — концентрации частиц в одном растворе различаются на много порядков, и абсолютный порог здесь неприменим.
Расчёт ведётся сразу по серии растворов, и аналитические концентрации задаются одним из двух способов на выбор: либо напрямую по базисным составляющим, либо через матрицу компонентов — тогда пользователь вводит концентрации реальных навесок (соль, кислота, буфер), а программа сама пересчитывает их в базисные. Второй способ ближе к тому, как готовится раствор на практике. Для каждого раствора возвращаются концентрации всех частиц, ионная сила, коэффициенты активности, логарифмы активностей базисных частиц, число итераций и невязка баланса заряда — последняя служит независимой проверкой: она обязана быть нулевой с точностью расчёта, и её отличие от нуля сразу указывает на ошибку в зарядах или в стехиометрии.
3.3 Чем ценно единообразие
Единая запись — не вопрос вкуса, у неё есть измеримые следствия.
- Новый вид химии не требует нового кода. Ассоциация ионных пар в кондуктометрии — это обычная строка матрицы; заряженный ассоциат остаётся в списке ионов и продолжает проводить ток, незаряженный выбывает из него сам собой.
- Осадок вводится как обычная частица. Гетерогенное равновесие решается надстройкой над тем же решателем, а не отдельной программой; при этом учитывается собственная остаточная растворимость осадка — после полного осаждения именно она определяет ЭДС, и расчёт без неё промахивается более чем на 100 мВ.
- Один и тот же файл питает разные методы. Описание системы, введённое для расчёта равновесий, целиком входит в потенциометрический и кондуктометрический проект — переописывать химию не нужно.
- Проверяемость. Матричная запись позволяет сравнивать независимые пути расчёта между собой: этим и пользуется набор автоматических тестов (§9).
3.4 Коэффициенты активности: четыре модели на выбор
Ион-ионное взаимодействие входит в расчёт равновесия косвенно, но входит всегда: коэффициенты активности стоят в законе действующих масс, а ионная сила зависит от найденного состава. Точность всего расчёта определяется в конечном счёте адекватностью того уравнения, которым описаны коэффициенты активности, — и потому модель в комплексе не вшита, а выбирается на первом шаге ввода данных. Выбор относится ко всей задаче и доходит до всех приложений автоматически: он хранится в равновесной части проекта.
| Модель | Уравнение или основа | Параметры | Область |
|---|---|---|---|
| Предельный закон Дебая — Хюккеля | $\lg\gamma_i = -A z_i^2\sqrt{I}$ — линеаризованное уравнение Пуассона — Больцмана для точечных ионов | нет | $I \lesssim 0.001$; эталон предельного перехода |
| Уравнение Дэвиса | $\lg\gamma_i = -A z_i^2\!\left(\dfrac{\sqrt{I}}{1+\sqrt{I}} - 0{,}3\,I\right)$ — размер иона принят одинаковым для всех | нет | $I \lesssim 0.1$–0.5 |
| Расширенное уравнение Дебая — Хюккеля | $\lg\gamma_i = -\dfrac{A z_i^2\sqrt{I}}{1 + B\mathring{a}\sqrt{I}} + b\,I$ — конечный размер иона плюс эмпирическая линейная поправка | $\mathring{a}$, $b$ — общие или по каждому иону | $I \lesssim 0.5$–1 |
| MSA (среднесферическое приближение) | решение уравнения Орнштейна — Цернике для смеси заряженных твёрдых сфер: электростатика плюс уравнение состояния твёрдых сфер | радиусы ионов | $I \lesssim 1$ и выше |
Разница между ступенями — не академическая. Для $\mathrm{MgCl_2}$ в воде при 25 °C коэффициент активности двухзарядного катиона выходит таким:
| $I$, моль/л | предельный закон | Дэвис | расширенное ур. | MSA | вывод |
|---|---|---|---|---|---|
| 0.01 | 0.625 | 0.661 | 0.660 | 0.678 | три «разумные» согласны до 3 % — выбор безразличен |
| 0.30 | 0.076 | 0.289 | 0.224 | 0.310 | расхождение 38 %, предельный закон — вчетверо |
Отсюда и постановка задачи: поднять порог применимости расчёта хотя бы до $I = 0{,}3$ и дать возможность узнать, насколько результат обязан модели. Второе рассуждением не решается — только счётом одного и того же разными моделями, и такой инструмент в комплексе есть (§4.5).
Отдельного упоминания стоит термодинамическая согласованность. Набор индивидуальных коэффициентов активности порождается избыточной энергией Гиббса тогда и только тогда, когда $\ln\gamma_i = z_i^2 g(I)$. У расширенного уравнения это нарушено дважды — линейный член не масштабирован зарядом, и разные размеры разных ионов делают зависимость от сорта частицы шире множителя $z_i^2$; такой энергии тогда не существует вовсе, а с нею не существует ни осмотического коэффициента, ни активности растворителя. На сам равновесный состав это не влияет (он определяется произведением коэффициентов по зарядо-сбалансированной реакции и к такому произволу нечувствителен), но влияет на всё, что выводится из растворителя. Это и есть главный довод в пользу MSA: ион-специфичность в рамках Дебая — Хюккеля покупается ценой согласованности, а MSA даёт индивидуальные размеры, её не теряя.
Как выбор модели проявляется в каждом приложении
| Приложение | Через что действует модель | Что из этого следует |
|---|---|---|
| Равновесия | Прямо на состав: сдвиг равновесий внутри заданного баланса | Пункт рельса «Сравнение моделей» считает проект несколькими моделями сразу. Мера расхождения — сдвиг $\lg c$, то есть ровно то, на сколько пришлось бы подвинуть $\lg K$, чтобы одна модель дала ответ другой |
| Потенциометрия | ЭДС измеряет активность, а не концентрацию: $a_i = c_i\gamma_i$ | Сдвиг зависит от ионной силы, а она по серии меняется, — поэтому в $E^{\circ}$ он не поглощается и проявляется ходом невязок по $I$. Ось «ионная сила» на графике невязок — самый чувствительный признак неверной модели |
| Осаждение и титрование | Условие насыщения $\prod_j a_j^{s_j} = K_{sp}$ записано через активности | Осадок считается той же моделью, что и гомогенная часть, — иначе точки до скачка и после отвечали бы разным моделям. Найденное $\mathrm{p}K_{sp}$ неотделимо от модели, при которой получено |
| Кондуктометрия | Через равновесный состав — на концентрации свободных ионов (§4.3), а под MSA ещё и напрямую | Один радиус служит и $\gamma$, и $\Lambda$. В совместном анализе параметров становится не больше, а меньше, и появляется проверяемое утверждение: размер иона, найденный по электропроводности, обязан воспроизводить коэффициенты активности |
Подробное изложение — лестница моделей с выводами и границами, разбор термодинамической согласованности, что отвергнуто и почему, устройство MSA, неводные среды и порядок выбора модели шестью шагами — в справке «Расчёт равновесного состава», раздел 5 (§8).
4Электрохимические методы: строгость равновесия и полнота переноса
Комплекс опирается на два метода, которые дополняют друг друга по самой своей природе. Потенциометрия измеряет равновесную величину и потому термодинамически строга. Электропроводность — величина неравновесная, зато она чувствительна к тому, чего потенциометрия не видит вовсе: к подвижности частиц и к присутствию незаряженных ассоциатов.
4.1 Равновесная потенциометрия: термодинамическая строгость
ЭДС ячейки связана с активностями потенциалопределяющих ионов уравнением Нернста
где $\sigma = \pm 1$ задаёт полярность электрода. Строгость метода в том, что здесь нет модельных допущений о механизме: измеряется именно активность, и связь её с потенциалом — точное термодинамическое соотношение. Всё модельное сосредоточено в одном месте — в коэффициентах активности, и их параметры при желании определяются из того же эксперимента.
Поддержаны ячейки с переносом (один потенциалопределяющий ион) и без переноса (два); диффузионный потенциал считается полностью устранённым. Повторные измерения в точке дают вес $1/\text{дисперсия}$, то есть точки с худшей воспроизводимостью автоматически получают меньший вклад.
| Метод | Семейство | Когда брать |
|---|---|---|
| Нелдера — Мида | прямой поиск производные не нужны |
Первый запуск: надёжен, нетребователен к начальному приближению |
| Хука — Дживса | Параметры очень разного масштаба (шаг подстраивается под каждый) | |
| Быстрый поиск (вращающиеся направления) | Изогнутый овраг: направления поворачиваются вслед за долиной | |
| Пауэлла | сопряжённые направления | Гладкая задача без производных, требуется быстрая сходимость у минимума |
| Брента (метод главных осей) | Трудная задача: базис перестраивается по главным осям самой задачи | |
| BFGS | градиентные | Много параметров, нужна оценка кривизны без второй производной |
| Сопряжённых градиентов (Полака — Рибьера) | То же при большом числе параметров: память растёт линейно | |
| Гаусса — Ньютона | по вторым производным и по невязкам |
Самый быстрый на гладкой задаче наименьших квадратов у самого минимума |
| Левенберга — Марквардта | Далёкое начальное приближение: шаг демпфируется автоматически | |
| Демпфированного Ньютона | Окончательное уточнение и проверка, что найден минимум, а не седло | |
| Дифференциальная эволюция | глобальный | Подозрение на несколько минимумов: ищет по всей области сразу. Дороже локальных на один-два порядка, поэтому включается осознанно |
4.2 Электропроводность: вся лестница теорий от классики до современного уровня
Электропроводность неравновесна, и её описание — это описание торможения иона двумя механизмами: релаксационным (ионная атмосфера не успевает перестроиться и тянет ион назад) и электрофоретическим (атмосфера движется навстречу, увлекая растворитель). Теория этих эффектов развивалась почти столетие, и комплекс содержит не одну «правильную» модель, а всю лестницу — от предельного закона до современных подходов.
Нижняя ступень лестницы — предельный закон Онзагера, где оба эффекта дают вклад, пропорциональный корню из концентрации:
причём слагаемое $\alpha\Lambda^{0}$ отвечает релаксации, а $\beta$ — электрофорезу. Дальнейшие приближения добавляют члены следующих порядков,
и различаются как раз тем, как получены коэффициенты $E$, $J_1$, $J_2$ и какое расстояние наибольшего сближения в них входит. Именно здесь классические модели расходятся между собой, и именно здесь у них общее ограничение: все они выведены для симметричного бинарного электролита, то есть для двух сортов ионов с равными по модулю зарядами. Реальный раствор с комплексами, ассоциатами или несколькими солями под это условие не подходит — и поэтому в интерфейсе предлагаются только четыре современные модели, работающие со смесью любого состава; классическая лестница остаётся в расчётном ядре как эталон для их проверки.
| Модель | Область | Чем важна |
|---|---|---|
| Предельный закон Онзагера | симметричный бинарный | точка отсчёта: к нему обязана сходиться любая модель при бесконечном разбавлении |
| Ли — Уитон (симметричная форма) | симметричный бинарный | замкнутая форма, эталон сравнения для общих движков |
| Фуосс — Онзагер — Скиннер | симметричный бинарный | классика с контактными парами (при явном учёте ассоциации отключается, чтобы пары не учитывались дважды) |
| Питтс; Фуосс — Хсиа | симметричный бинарный | два конкурирующих семейства высших приближений |
| Чен; Фуосс — Джастис — Чен | симметричный бинарный | формы, принятые в обзорной литературе по обработке данных |
| Квинт — Вьяллар (чистый электролит) | бинарный любого типа | замкнутая форма для несимметричных солей |
| Квинт — Вьяллар (смесь) | смесь любого типа | спектральный подход; работает с заряженным ассоциатом как с обычным ионом |
| Ли — Уитон (полная смесевая форма) | смесь любого типа | разложение по степеням с перекалиброванным спектральным каркасом |
| Ли — Уитон (прямая форма) | смесь любого типа | полное уравнение без разложения; индивидуальные расстояния для каждого иона |
| MSA-перенос | смесь любого типа | современный подход: индивидуальные диаметры ионов вводятся изначально, а не подгоняются одним общим параметром |
Насколько это работает — вопрос не мнения, а сверки с опубликованными таблицами. Модель MSA воспроизводит данные по MgCl2 с отклонением 0.18 %, по LaCl3 — менее 0.6 %, а смесь NaCl + MgCl2 — вообще без подгоночных параметров. Границы применимости названы честно и проверены: ионная сила примерно до 1–2 M для однозарядных, до 0.5 M для 2:1 и до 0.06 M для 3:1.
4.3 Шов между разделами: модель работает по свободным ионам
Здесь сходятся §3 и §4, и это архитектурное решение всего комплекса. Ток переносят свободные ионы, а не аналитические концентрации; доля свободных ионов — результат решения равновесной задачи. Поэтому расчёт идёт цепочкой: матрица стехиометрии → решатель → свободные концентрации и ионная сила → модель переноса → κ. Ассоциация при этом не является особой теорией внутри модели проводимости: она вводится обычной строкой матрицы со своей $\lg K$.
4.4 Титрование: два разных вопроса — четыре окна
Титрование присутствует в обоих методах, и в каждом оно разделено на две самостоятельные задачи, потому что задачи эти требуют разного.
- Моделирование и обработка кривой целиком. Смешение объёмов считается физически (аликвота + разбавитель + добавленный титрант), в каждой точке решается равновесие, и по всей кривой ищутся константы. Здесь кривая — источник информации о системе: в одной кривой сотни точек при разном соотношении реагентов, и этого часто достаточно, чтобы определить то, что по отдельным растворам не определяется.
- Анализ: точка эквивалентности и концентрация. Нужна только измеренная кривая. В потенциометрии эквивалентность ищется по скачку — производной, пересечению ветвей $\Delta V/\Delta E$ и четырьмя функциями Грана; программа сама судит, какой оценке верить, и показывает две, когда данные не позволяют выбрать. В кондуктометрии ищется излом — пересечение прямых ветвей, а погрешность получается переносом ковариаций обеих подгонок, то есть это результат с $\sigma$, а не оценка на глаз.
На поставляемом примере титрования соляной кислоты три независимые оценки дают 10.0437 ± 0.021, 10.0144 и 10.0001 мл при истинном значении 10.00 — они сходятся в пределах одного деления бюретки, и именно это согласие независимых способов является настоящей проверкой результата.
4.5 Достоверность решения: интервалы, корреляции, выбор модели
Подгонка всегда что-нибудь находит. Вопрос в том, что именно найдено: величина, определённая данными, или число, которое лишь удачно легло на кривую. Комплекс отвечает на этот вопрос тремя инструментами, и они отвечают на три разных вопроса — ни один не заменяет остальных.
| Инструмент | На какой вопрос отвечает | Где смотреть |
|---|---|---|
| Доверительный интервал | Насколько плохо определён параметр сам по себе | вкладка «Решение» |
| Корреляционная матрица | Не определена ли вместо двух параметров лишь их комбинация | вкладка «Корреляции» |
| Сравнение моделей активности | Верна ли сама модель, в рамках которой всё это найдено | вкладка «Сравнение моделей» |
Корреляционная матрица
Второй вопрос по интервалам не виден вовсе. Оба интервала могут выглядеть приемлемо, а поверхность целевой функции быть длинным узким оврагом, вдоль которого решение гуляет свободно: сдвиньте один параметр — и второй сдвинется вслед, вернув сумму квадратов почти к прежнему значению. Именно так выглядит задача, в которой определена комбинация, а не составляющие.
Матрица получается из того же обратного гессиана, которым считаются интервалы,
причём множитель $\sigma^2$ здесь сокращается — поэтому вывод корреляций не стоит ни одного дополнительного вычисления целевой функции и доступен после любого расчёта.
| $|r|$ | Что это значит | Что делать |
|---|---|---|
| < 0.9 | параметры разделяются приемлемо | — |
| 0.9 – 0.99 | разделяются плохо: настоящая неопределённость каждого шире, чем следует из его интервала | расширить диапазон условий опыта; приводить корреляцию вместе с результатом |
| ≥ 0.99 | практически неразделимы: данные определяют только комбинацию | закрепить один параметр независимыми данными либо не истолковывать значения по отдельности |
В задачах этого комплекса высокая корреляция — норма, а не исключение, и знать её типичные источники полезнее, чем удивляться им:
- $\lg K \leftrightarrow E^{\circ}$. Константа сдвигает расчётные активности почти равномерно по серии, а стандартная ЭДС сдвигает всю кривую; на данных, снятых в узком диапазоне ионной силы, различить их нечем.
- $\lambda^{0} \leftrightarrow \lg K$ в кондуктометрии: и рост подвижности, и ослабление ассоциации повышают электропроводность. Это фундаментальное свойство задачи, а не дефект подгонки, — и именно оно делает необходимыми одиннадцать методов поиска решения (§4.1).
- Параметры модели активности $\leftrightarrow$ и то, и другое. Прямая причина, по которой задание $\mathring{a}$ и $b$ по отдельным ионам выключено по умолчанию: каждый добавленный параметр отбирает определённость у констант.
Ячейки с $|r| \ge 0.9$ подсвечиваются, с $|r| \ge 0.99$ выделяются; прочерк означает, что величины не существует — дисперсия по этому параметру неположительна, то есть найденная точка не является минимумом по нему. Это само по себе диагноз. Матрица выгружается в Excel отдельным листом и копируется в буфер наравне с остальными таблицами.
Сравнение моделей по качеству подгонки
Третий вопрос — о самой модели — данные решают сами, если есть измеренная величина. Вкладка «Сравнение моделей» подгоняет задачу заново каждой отмеченной моделью коэффициентов активности и раскладывает результаты рядом.
Три исхода сравнения требуют разных действий, и худший из них — не тот, о котором думают:
- Одна модель выигрывает решающе — берите её и указывайте вместе с константами.
- Разрыв невелик — берите лучшую, но проверьте, сохраняется ли преимущество при других весах и на другой части серии.
- Данные модели не различают, а состав от модели зависит. Самый частый и самый неприятный случай: эксперимент не выбирает модель, но ответ от неё зависит. Честный выход — привести результат для нескольких моделей и указать разброс. Это не слабость работы, а измеренная величина неопределённости, которую иначе пришлось бы замалчивать.
Инструменты доступны в окнах оптимизации потенциометрии и кондуктометрии — там, где есть эксперимент. В «Равновесиях» сравнивается только состав: опыта там нет, и сумме квадратов взяться неоткуда. Подробный разбор обеих вкладок — в справке «Расчёт равновесного состава», раздел 5.9, и в справке «Потенциометрия и поиск решения», раздел 7.
5Математическая обработка: полноценный раздел, а не приложение
Эксперимент редко даёт числа, готовые к расчёту. Их нужно упорядочить, очистить от случайной погрешности, снять с них производную или интеграл, подобрать закон, оценить достоверность и представить результат. В комплексе это отдельный раздел из восьми подпунктов и 62 численных методов, а не набор кнопок внутри химических окон.
5.1 Три места в работе
| Этап | Что делается | Подпункты |
|---|---|---|
| До расчёта подготовка первичного материала |
Упорядочить измерения, отбраковать и сгладить шум, восстановить значение между точками. Сюда же — получение физико-химических параметров растворителя: плотность, вязкость и диэлектрическая проницаемость обычно заданы справочной таблицей от состава и температуры сразу, и значение при нужных условиях снимается двумерной аппроксимацией | сортировка, сглаживание, интерполяция, сплайны, 2D аппроксимация |
| Во время расчёта | Производная кривой титрования, площадь под пиком, восстановление плавной кривой по разреженным точкам | дифференцирование, интегрирование, сплайны |
| После расчёта анализ результата |
Проверить невязки на систематичность, подобрать зависимость найденного параметра от температуры или состава, сравнить модели между собой, подготовить итоговые таблицы и графики | аппроксимация, типовые уравнения термодинамики, 2D аппроксимация |
5.2 Восемь подпунктов, 62 метода
| Подпункт | Методов | Что внутри |
|---|---|---|
| Сортировка | 2 режима | Вектор или таблица целиком: сортировка по выбранной колонке с переносом всей строки, устойчивая (равные ключи сохраняют исходный порядок) |
| Сглаживание | 13 | От скользящего среднего и фильтра Савицкого — Голея до устойчивых к выбросам (медиана, Тьюки, LOWESS) и сохраняющего скачок TVD; выбор — по численным характеристикам, а не на глаз |
| Интерполяция | 9 | Все проходят точно через узлы: три формы одного многочлена, барицентрическая форма, локальные кубики Акимы и PCHIP, рациональные Тиле и Флотера — Хормана; производные и обратная задача |
| Сплайны | 8 | Кубический с пятью краевыми условиями, эрмитов, Катмулла — Рома, кардинальный, B-сплайн, с натяжением, Акимы, монотонный; выбор по классу гладкости и локальности |
| Аппроксимация | 8 | Кривая проводится не через точки: полином, ортогональный полином Форсайта, МНК по набору функций, типовые уравнения термодинамики, кривая Безье и три сплайновых метода; веса, устойчивость к выбросам, нормы $L^1$ и $L^\infty$ |
| 2D аппроксимация | 5 | Функция двух переменных, заданная таблицей: билинейная, бикубическая, тензорный сплайн (в том числе сглаживающий), множественная регрессия со статистикой, двумерный полином |
| Дифференцирование | 9 | Задача некорректна, поэтому окно устроено как инструмент сравнения: разностные схемы, локальный многочлен, регуляризация Тихонова и TVD, производные построенных кривых; L-кривая и самопроверка методом Монте-Карло |
| Интегрирование | 10 | По таблице и по кривой: от прямоугольников и трапеций до Ромберга, Гаусса — Лежандра, Гаусса — Кронрода и адаптивных схем; с оценкой погрешности метода и отдельно — неопределённости от шума |
5.3 Тот же раздел — как самостоятельный инструмент
Раздел поставляется и отдельной программой DataAnalysisLab, в которой нет ни одной химической строки: ни решателя равновесий, ни уравнения Нернста, ни моделей проводимости. Это не декларация — отсутствие таких связей проверяется автоматическим тестом при каждой сборке, а сам файл программы занимает 5.8 МБ против 26.5 МБ у полной оболочки.
Практический смысл прост: те же методы применимы в физико-химическом эксперименте любого направления — кинетическом (порядок реакции по производной), спектрофотометрическом (площадь и положение полосы), хроматографическом (площадь пика), термохимическом (тепловой эффект как интеграл), — и коллеге, которому нужна только обработка, не нужно объяснять, что такое базисные частицы.
5.4 Статистика: у каждого числа есть его погрешность
Сквозное правило комплекса: результат без оценки достоверности результатом не считается. Поэтому рядом с каждым найденным значением стоит его погрешность, а рядом с каждой подобранной моделью — признаки того, годится ли она.
- Параметры — доверительные интервалы из обратного гессиана, с проверкой его положительной определённости; непройденная проверка означает, что параметры по этим данным неразличимы, и об этом говорится прямо.
- Кривая — доверительный интервал (о положении кривой) и коридор ошибок (о разбросе отдельных измерений) — это разные вещи, и они не смешиваются.
- Выбор модели — сумма квадратов монотонно падает с числом параметров, поэтому по ней выбирать нельзя. Считаются AICc (с обязательной поправкой на малую выборку — в химическом опыте точек мало), BIC, GCV, PRESS, а также число смен знака невязки — единственная мера, ничего не знающая о числе параметров: мало смен означает, что модель неверна по форме, а не что шум велик.
- Коллинеарность — VIF и число обусловленности: они показывают, когда красивые коэффициенты нельзя истолковывать порознь.
- Регрессионный анализ — значимость каждого коэффициента ($t$, $p$) и дисперсионный анализ модели в целом.
5.5 Наглядный графический анализ
График в комплексе — рабочий инструмент, а не иллюстрация к отчёту. Всюду, где есть расчёт, есть и графики с выбором того, что отложено по осям: экспериментальные и расчётные кривые, невязки (по номеру раствора, по концентрации, по ионной силе или по объёму титранта — систематический ход виден сразу), вклады отдельных ионов в электропроводность, равновесные концентрации выбранных частиц в логарифмической шкале, ход сходимости оптимизации, доверительные полосы, а для функции двух переменных — сечения, изолинии и карта невязок. Экспериментальные точки рисуются полыми кружками, расчётные — залитыми и мельче, поэтому совпадение видно как точка внутри кольца, а не как слипшееся пятно.
5.6 Форматы, сохранение и перенос данных
Данные не запираются внутри программы.
- Проект — открытый XML, четырнадцать типов по виду задачи (
*.eqp.xml,*.ptp.xml,*.cnd.xml,*.spl.xml, …). Файл читается человеком и текстовыми средствами. - Файл хранит рецепт, а не готовые числа. Сглаженные и подогнанные кривые пересчитываются при открытии, поэтому файл в принципе не может разойтись с алгоритмом — важное свойство для данных, которые хранятся годами.
- Ввод — вставкой из буфера (Ctrl+V) прямо из Excel или
лабораторного журнала: таблица сама вырастает под размер вставленного блока, а
числа разбираются в свободном формате (
0,5,1.2E-3,1.2*10^-3). - Вывод — копирование любой таблицы в буфер как TSV (одним пунктом меню, без
диалогов) и экспорт результатов в Excel (
.xlsx) в химических модулях; графики копируются как изображение средствами окна.
6Требования и особенности применения
| Что | Требуется |
|---|---|
| Операционная система | Windows 10 или 11, 64 бита |
| Среда выполнения | .NET 10 Desktop Runtime — распространяется бесплатно корпорацией Microsoft |
| Для показа справки | Microsoft Edge WebView2 Runtime (входит в состав актуальной Windows). Без него работают все расчёты, недоступны только страницы справки |
| Процессор и память | Любой современный; расчёт однопоточный, память — десятки мегабайт. Отдельная видеокарта не нужна |
| Экран | От 1280×800. Масштабирование 100–200 % поддержано полностью: окна объявлены DPI-осведомлёнными, размеры и графики пересчитываются под монитор |
| Установка | Не требуется. Каждая программа — один файл .exe (5.8–26.5 МБ): скопировать и запустить, права администратора не нужны |
| Сеть | Не используется вовсе. Расчёты локальны, справка со всеми формулами работает офлайн — комплекс пригоден для лабораторной машины без интернета и для закрытого контура |
| Где хранятся настройки | В профиле пользователя
(%LOCALAPPDATA%); в системные каталоги и реестр программа не пишет |
Что считается быстро, а что долго
Знать это полезно при планировании работы.
- Мгновенно — расчёт равновесия в одной точке (миллисекунды), любой метод обработки данных на таблице в сотни точек, весь анализ кривой титрования.
- Секунды — подгонка локальным методом по серии растворов или по кривой: равновесие кэшируется, и при изменении только $\lambda^{0}$ или $E^{\circ}$ решатель повторно не запускается.
- Заметно дольше — глобальная дифференциальная эволюция (каждая проба меняет константы, кэш не срабатывает, и полный расчёт идёт на каждую пробу) и задачи с осадком (в каждой точке дополнительно решается уравнение насыщения). Это плата за надёжность, и включаются такие режимы осознанно.
Границы объёма данных
Ограничения заданы не произвольно, а вытекают из численных схем, и в этих случаях программа отказывается считать с объяснением, а не выдаёт правдоподобный неверный результат: регуляризованное дифференцирование решает плотную задачу и потому ограничено 600 точками; погрешность интеграла от шума считается точным прогоном до 400 точек, выше — статистической оценкой. Границы применимости моделей электропроводности по ионной силе названы в §4.2.
7С чего начать
В комплекте есть готовые примеры (папка docs/examples): они открываются
как обычные проекты и заодно служат образцом формата файла. Ниже — три коротких маршрута.
- Откройте модуль Равновесия и пройдите ввод данных по шагам: растворитель → частицы → матрица стехиометрии и $\lg K$ → концентрации.
- Нажмите «Готово» — таблица состава и графики появятся сразу.
- Измените одну константу и повторите: видно, какая частица за что отвечает.
- Откройте Титрование ▸ Анализ (потенциометрическое или кондуктометрическое — по вашему методу).
- Вставьте два столбца из журнала или из Excel через Ctrl+V; таблица вырастет сама.
- Пройдите этапы по рейлу слева и прочтите итог: точка эквивалентности и концентрация — каждая со своей погрешностью, плюс сравнение независимых оценок.
- Опишите систему в модуле Равновесия и убедитесь, что прямой расчёт даёт разумный состав.
- Перейдите в Потенциометрию или Кондуктометрию, введите измерения и веса.
- На вкладке оптимизации отметьте искомые параметры, задайте границы и запустите поиск. Начните с метода Нелдера — Мида, затем повторите другим методом: совпадение результатов разных методов — первая проверка, а расхождение говорит о том, что минимум неглубок или параметры коррелируют.
- Посмотрите доверительные интервалы и график невязок: систематический ход невязки означает, что дело не в точности данных, а в модели.
8Карта справок комплекса
Эта страница — общая. Подробности каждого метода, с выводом формул и разбором численных схем, изложены в отдельных страницах; все они открываются из меню «Справка» и содержат в сумме 226 ссылок на литературу — от оригинальных статей 1920-х годов до работ последних лет, с указанием, где именно в программе использован каждый источник.
| Страница | О чём | Ссылок |
|---|---|---|
| Расчёт равновесного состава | Формализм Бринкли, вывод системы уравнений, коэффициенты активности, порядок работы | 16 |
| Потенциометрия и поиск решения | Уравнение Нернста, типы ячеек, веса, все одиннадцать методов оптимизации по отдельности, доверительные интервалы | 10 |
| Потенциометрическое титрование | Теория скачка, осаждение, обработка кривой целиком и анализ по эксперименту: производные, ветви $\Delta V/\Delta E$, функции Грана | 9 |
| Кондуктометрия — модели и решение | Два эффекта торможения, все модели электропроводности с областями применимости, идеология планирования опыта | 21 |
| Кондуктометрическое титрование | Почему проводимость линейна по объёму, ячейка и фон растворителя, разбавление, вывод формул анализа излома | 8 |
| Математическая обработка | Все восемь подпунктов: каждый метод с формулой, областью применения и оригинальной ссылкой | 162 |
9Чем проверен комплекс
Программа, считающая константы, ценна ровно настолько, насколько ей можно верить. Проверка здесь ведётся четырьмя независимыми способами, и все они автоматические — они выполняются заново при каждой сборке.
- Более тысячи автоматических тестов (1117 на момент выпуска этой версии). Они проверяют не «программа не упала», а свойства результата: что метод воспроизводит точное решение там, где оно известно; что предельный переход даёт классическую формулу; что округление формата не портит данные.
- Сверка с опубликованными таблицами. Модели электропроводности сопоставлены с экспериментальными данными оригинальных статей — по хлориду магния, хлориду лантана, смеси хлоридов натрия и магния, сульфату магния, хлориду кадмия, сульфатам натрия и калия, — с теми же отклонениями, что названы в §4.2.
- Тождества между независимыми путями. Там, где одну величину можно получить двумя способами, совпадение проверяется как обязательное: смесевая модель обязана сводиться к бинарной, любая модель — к предельному закону Онзагера при бесконечном разбавлении, все четыре формы интерполяционного многочлена — совпадать между собой, сечение двумерной поверхности вдоль узловой линии — в точности равняться одномерному сплайну по тому же столбцу.
- Воспроизводимость. Методы, использующие случайные числа (дифференциальная эволюция, оценка следа, самопроверка Монте-Карло), работают от фиксированного зерна: повторный запуск на тех же данных даёт тот же ответ, и результат в статье можно воспроизвести спустя годы.
10Чего комплекс пока не делает
Границы возможностей названы здесь намеренно: специалисту важнее знать их заранее, чем обнаружить в середине работы.
- Кинетика. Комплекс описывает равновесие и стационарный перенос; скорости реакций, диффузионные задачи и нестационарные методы (вольтамперометрия, хроноамперометрия, импеданс) не рассматриваются.
- Только водные и однородные по составу растворители. Свойства воды заданы формулами от температуры; для другого растворителя его диэлектрическую проницаемость и вязкость нужно задать самому (что предусмотрено полями ввода).
- Двумерная аппроксимация — только по регулярной сетке. Рассеянные точки произвольного расположения не поддерживаются.
- Расчёт по свободным ионам в моделях проводимости ограничен областью применимости самих моделей по ионной силе (§4.2); за её пределами программа посчитает, но результат не будет физически осмысленным.
- Экспорт в Excel реализован в химических модулях; в разделе математической обработки пока доступно копирование таблиц через буфер обмена.
- Уравнений Питцера и SIT нет, и это решение, а не пробел. Первые вобрали ионную ассоциацию в свои коэффициенты и при работе с явными комплексами дали бы двойной учёт; вторая опирается на водную базу параметров, тогда как комплекс рассчитан на любые среды (§3.4). Для концентрированных водных растворов, где Питцер — стандарт, нужен именно он, и здесь его нет.
- Комплекс не заменяет постановку опыта. Ни один численный метод не восстановит информацию, которой в данных нет: если параметры коррелируют, а точек мало, честный ответ — «по этим данным неразличимо», и программа его даёт.