Журнал физической химии, 2021, T. 95, № 7, стр. 1118-1123

Изменение энергии гиббса распределительной системы фракционирования при разделении многокомпонентной смеси

В. П. Чижков a, В. Н. Бойцов a*

a Российская академия наук. Институт органической химии им. Н.Д. Зелинского
Москва, Россия

* E-mail: boysoft@list.ru

Поступила в редакцию 05.07.2020
После доработки 05.07.2020
Принята к публикации 14.09.2020

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

Аннотация

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

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

При термодинамической интерпретации эффективности процессов разделения важное значение приобретает термодинамическое обоснование и расчет функций, определяющих работу разделения [17]. Максимальная полезная работа при заданных температуре и давлении равна убыли энергии Гиббса термодинамической системы [8, 9]. Процесс распределения разделяемых компонентов между двумя фазами (жидкость–жидкость, жидкость–сорбент) рассматривается обычно как изобарно-изотермический [10, 11]. В [12] представлен расчет изменения энергии Гиббса в системах фракционной экстракции при разделении пары компонентов (А и В) исходной (перерабатываемой в каскаде [1315]) смеси. Цель настоящей работы состоит в распространении предложенного подхода на смеси, содержащей более чем одну пару разделяемых компонентов [7].

Рассмотрим в качестве примера процесс разделения компонентов A, B, C и D многокомпонентной смеси, содержащей две выбранные пары AB и CD, в распределительной системе (см. рис. 1) в условиях независимого распределения разделяемых компонентов между двумя несмешивающимися растворителями – экстрагентами E1 и E2. При использовании взаимонерастворимых экстрагентов в количествах $E_{1}^{0}$ и $E_{2}^{0}$ молей в систему с одной стороны подаются два бесконечно разбавленных раствора, содержащих разделяемые компоненты в количествах A01, B01, C01, D01 и A02, B02, C02, D02 молей, а с другой стороны выводятся два аналогичных раствора, содержащих A1, B1, C1, D1 и A2, B2, C2, D2 молей разделяемых компонентов. Данная система представляет собой ступень каскада, подобного описанному в [7, 1315]. Эффект однократного разделения выбранных пар AB и CD охарактеризуем, как это принято в общей теории процессов разделения [1, 16], коэффициентами разделения α и эффективными коэффициентами разделения β:

${{\alpha }_{{{\text{AB}}}}} = K_{{\text{B}}}^{*}{\text{/}}K_{{\text{A}}}^{*} = \frac{{{{B}_{1}}{\text{/}}{{B}_{0}}}}{{{{A}_{1}}{\text{/}}{{A}_{0}}}}{{\left[ {\frac{{1 - {{B}_{1}}{\text{/}}{{B}_{0}}}}{{1 - {{A}_{1}}{\text{/}}{{A}_{0}}}}} \right]}^{{ - 1}}} = \frac{{{{B}_{1}}{\text{/}}{{A}_{1}}}}{{{{B}_{2}}{\text{/}}{{A}_{2}}}},$
${{\alpha }_{{{\text{CD}}}}} = K_{{\text{D}}}^{*}{\text{/}}K_{{\text{C}}}^{*} = \frac{{{{D}_{1}}{\text{/}}{{D}_{0}}}}{{{{C}_{1}}{\text{/}}{{C}_{0}}}}{{\left[ {\frac{{1 - {{D}_{1}}{\text{/}}{{D}_{0}}}}{{1 - {{C}_{1}}{\text{/}}{{C}_{0}}}}} \right]}^{{ - 1}}} = \frac{{{{D}_{1}}{\text{/}}{{C}_{1}}}}{{{{D}_{2}}{\text{/}}{{C}_{2}}}},$
${{\beta }_{{{\text{AB}}}}} = \frac{{{{B}_{1}}{\text{/}}{{B}_{0}}}}{{{{A}_{1}}{\text{/}}{{A}_{0}}}},$
${{\beta }_{{{\text{CD}}}}} = \frac{{{{D}_{1}}{\text{/}}{{D}_{0}}}}{{{{C}_{1}}{\text{/}}{{C}_{0}}}}.$
Отсюда
${{B}_{1}}{\text{/}}{{B}_{0}} = ({{\alpha }_{{{\text{AB}}}}}--{{\beta }_{{{\text{AB}}}}}){\text{/}}({{\alpha }_{{{\text{AB}}}}}--1) = {{q}_{{\text{B}}}},$
${{B}_{2}}{\text{/}}{{B}_{0}} = ({{\beta }_{{{\text{AB}}}}}--1){\text{/}}({{\alpha }_{{{\text{AB}}}}}--1) = {{p}_{{\text{B}}}},$
${{A}_{1}}{\text{/}}{{A}_{0}} = ({{\alpha }_{{{\text{AB}}}}}--{{\beta }_{{{\text{AB}}}}}){\text{/}}{{\beta }_{{{\text{AB}}}}}({{\alpha }_{{{\text{AB}}}}}--1) = {{q}_{{\text{A}}}},$
${{A}_{2}}{\text{/}}{{A}_{0}} = {{\alpha }_{{{\text{AB}}}}}({{\beta }_{{{\text{AB}}}}}--1){\text{/}}{{\beta }_{{{\text{AB}}}}}({{\alpha }_{{{\text{AB}}}}}--1) = {{p}_{{\text{A}}}},$
${{D}_{1}}{\text{/}}{{D}_{0}} = ({{\alpha }_{{{\text{CD}}}}}--{{\beta }_{{{\text{CD}}}}}){\text{/}}({{\alpha }_{{{\text{CD}}}}}--1) = {{q}_{{\text{D}}}},$
${{D}_{2}}{\text{/}}{{D}_{0}} = ({{\beta }_{{{\text{CD}}}}}--1){\text{/}}({{\alpha }_{{{\text{CD}}}}}--1) = {{p}_{{\text{D}}}},$
${{C}_{1}}{\text{/}}{{C}_{0}} = ({{\alpha }_{{{\text{CD}}}}}--{{\beta }_{{{\text{CD}}}}}){\text{/}}{{\beta }_{{{\text{CD}}}}}({{\alpha }_{{{\text{CD}}}}}--1) = {{q}_{{\text{C}}}},$
${{C}_{2}}{\text{/}}{{C}_{0}} = {{a}_{{{\text{CD}}}}}({{\beta }_{{{\text{CD}}}}}--1){\text{/}}{{\beta }_{{{\text{CD}}}}}({{\alpha }_{{{\text{CD}}}}}--1) = {{p}_{{\text{C}}}}.$
Здесь αCD > αAB [7], 1 < β < α [16] и
${{A}_{0}} = {{A}_{1}} + {{A}_{2}} = {{A}_{{01}}} + {{A}_{{02}}},$
${{B}_{0}} = {{B}_{1}} + {{B}_{2}} = {{B}_{{01}}} + {{B}_{{02}}},$
${{С}_{0}} = {{С}_{1}} + {{С}_{2}} = {{С}_{{01}}} + {{С}_{{02}}},$
${{D}_{0}} = {{D}_{1}} + {{D}_{2}} = {{D}_{{01}}} + {{D}_{{02}}}.$
Обозначения p и q широко используются в жидкостной экстракции [13]. Смысл этих обозначений следует из соотношений:
$p = 1--q = \varepsilon {\text{/}}(1 + \varepsilon ) = 1{\text{/}}(1 + k),$
где ε – факторы экстракции [13] и k – факторы хроматографического удерживания [11, 17] компонентов A, B, C и D.

Рис. 1.

Распределительная система фракционирования.

Термодинамические коэффициенты распределения [10, 11]

$K_{{\text{B}}}^{*} \cong \frac{{{{B}_{1}}{\text{/}}E_{1}^{0}}}{{{{B}_{2}}{\text{/}}E_{2}^{0}}},\quad K_{{\text{A}}}^{*} \cong \frac{{{{A}_{1}}{\text{/}}E_{1}^{0}}}{{{{A}_{2}}{\text{/}}E_{2}^{0}}},$
$K_{{\text{C}}}^{*} \cong \frac{{{{C}_{1}}{\text{/}}E_{1}^{0}}}{{{{C}_{2}}{\text{/}}E_{2}^{0}}},\quad K_{{\text{D}}}^{*} \cong \frac{{{{D}_{1}}{\text{/}}E_{1}^{0}}}{{{{D}_{2}}{\text{/}}E_{2}^{0}}}.$
Далее будем постулировать
${{(K_{{\text{B}}}^{*}K_{{\text{A}}}^{*})}^{{1/2}}} = {{(K_{{\text{D}}}^{*}K_{{\text{C}}}^{*})}^{{1/2}}} = K_{{\text{F}}}^{*}$
и
${{({{k}_{{\text{B}}}}{{k}_{{\text{A}}}})}^{{1/2}}} = {{({{k}_{{\text{D}}}}{{k}_{{\text{C}}}})}^{{1/2}}} = {{k}_{{\text{F}}}},$
где F – воображаемый компонент [17] (kD > kB > > kF > kA > kC) и
$k = \frac{{E_{1}^{0}}}{{E_{2}^{0}}}K{\kern 1pt} {\text{*}}.$
Постулируемые равенства наблюдаются, например [17, 18], при хроматографическом разделении веществ гомологического ряда. Для A, B, C и D, подбирая экстрагенты и условия разделения, возможно, как и в случае отдельного разделения пар AB и CD [13], установить, по крайней мере, в принципе
$E_{1}^{0}{\text{/}}E_{2}^{0} = 1{\text{/}}{{(K_{{\text{B}}}^{*}K_{{\text{A}}}^{*})}^{{1/2}}} = 1{\text{/}}{{(K_{{\text{D}}}^{*}K_{{\text{C}}}^{*})}^{{1/2}}}.$
Эти равенства отвечают симметричному разделению [13]. При симметричном разделении получают наилучшие результаты [13].

Для рассматриваемой системы (см. рис. 1) принимаем [10, 12]

$\begin{gathered} \frac{{E_{1}^{0}}}{{E_{2}^{0}}} = {{k}_{0}} = \frac{{{{A}_{{01}}}}}{{{{A}_{{02}}}}} = \frac{{{{B}_{{01}}}}}{{{{B}_{{02}}}}} = \frac{{{{C}_{{01}}}}}{{{{C}_{{02}}}}} = \frac{{{{D}_{{01}}}}}{{{{D}_{{02}}}}} = \frac{{{{A}_{{01}}} + {{B}_{{01}}}}}{{{{A}_{{02}}} + {{B}_{{02}}}}} = \\ = \frac{{{{C}_{{01}}} + {{D}_{{01}}}}}{{{{C}_{{02}}} + {{D}_{{02}}}}} = \frac{{{{A}_{{01}}} + {{B}_{{01}}} + {{C}_{{01}}} + {{D}_{{01}}}}}{{{{A}_{{02}}} + {{B}_{{02}}} + {{C}_{{02}}} + {{D}_{{02}}}}}. \\ \end{gathered} $
Оптимальным условиям разделения при экстракции отвечают соотношения [12, 13, 16]:
$E_{1}^{0} = E_{2}^{0},$
${{(K_{{\text{B}}}^{*}K_{{\text{A}}}^{*})}^{{1/2}}} = {{(K_{{\text{D}}}^{*}K_{{\text{C}}}^{*})}^{{1/2}}} = 1{\text{/}}{{k}_{0}} = 1,$
${{\beta }_{{{\text{AB}}}}} = \alpha _{{{\text{AB}}}}^{{1/2}} = \frac{{{{A}_{2}}{\text{/}}{{A}_{0}}}}{{{{B}_{2}}{\text{/}}{{B}_{0}}}},$
${{\beta }_{{{\text{CD}}}}} = \alpha _{{{\text{CD}}}}^{{1/2}} = \frac{{{{C}_{2}}{\text{/}}{{C}_{0}}}}{{{{D}_{2}}{\text{/}}{{D}_{0}}}},$
${{B}_{1}}{\text{/}}{{B}_{0}} = \alpha _{{{\text{AB}}}}^{{1/2}}{\text{/}}(\alpha _{{{\text{AB}}}}^{{1/2}} + 1) = {{A}_{2}}{\text{/}}{{A}_{0}},$
${{B}_{2}}{\text{/}}{{B}_{0}} = 1{\text{/}}(\alpha _{{{\text{AB}}}}^{{1/2}} + 1) = {{A}_{1}}{\text{/}}{{A}_{0}},$
${{D}_{1}}{\text{/}}{{D}_{0}} = \alpha _{{{\text{CD}}}}^{{1/2}}{\text{/}}(\alpha _{{{\text{CD}}}}^{{1/2}} + 1) = {{C}_{2}}{\text{/}}{{C}_{0}},$
${{D}_{2}}{\text{/}}{{D}_{0}} = 1{\text{/}}(\alpha _{{{\text{CD}}}}^{{1/2}} + 1) = {{C}_{1}}{\text{/}}{{C}_{0}},$
$\Delta {{n}_{{\text{B}}}} = {{B}_{1}}--{{B}_{{01}}} = {{B}_{{02}}}--{{B}_{2}} = \frac{{{{B}_{0}}}}{2}\frac{{\alpha _{{{\text{AB}}}}^{{1/2}} - 1}}{{\alpha _{{{\text{AB}}}}^{{1/2}} + 1}},$
$\Delta {{n}_{{\text{A}}}} = {{A}_{2}}--{{A}_{{02}}} = {{A}_{{01}}}--{{A}_{1}} = \frac{{{{A}_{0}}}}{2}\frac{{\alpha _{{{\text{AB}}}}^{{1/2}} - 1}}{{\alpha _{{{\text{AB}}}}^{{1/2}} + 1}},$
$\Delta {{n}_{{\text{D}}}} = {{D}_{1}}--{{D}_{{01}}} = {{D}_{{02}}}--{{D}_{2}} = \frac{{{{D}_{0}}}}{2}\frac{{\alpha _{{{\text{CD}}}}^{{1/2}} - 1}}{{\alpha _{{{\text{CD}}}}^{{1/2}} + 1}},$
$\Delta {{n}_{{\text{C}}}} = {{C}_{2}}--{{C}_{{02}}} = {{C}_{{01}}}--{{C}_{1}} = \frac{{{{C}_{0}}}}{2}\frac{{\alpha _{{{\text{CD}}}}^{{1/2}} - 1}}{{\alpha _{{{\text{CD}}}}^{{1/2}} + 1}},$
которые характеризуют также симметричное разделение молекулярных и изотопных смесей в общей теории процессов разделения (βAB = $\alpha _{{{\text{AB}}}}^{{1/2}}$, βCD = $\alpha _{{{\text{CD}}}}^{{1/2}}$ [1, 16]).

Термодинамический эффект разделения заключается в изменении свободной энергии системы и энтропии смешения [4, 5, 19]. Энергия Гиббса системы (рис. 1) в начальном состоянии [10, 20]

$\begin{gathered} {{G}_{{\text{Н}}}} = {{A}_{{01}}}\mu _{{\text{A}}}^{{{{E}_{1}}}} + {{B}_{{01}}}\mu _{{\text{B}}}^{{{{E}_{1}}}} + {{C}_{{01}}}\mu _{{\text{C}}}^{{{{E}_{1}}}} + {{D}_{{01}}}\mu _{{\text{D}}}^{{{{E}_{1}}}} + E_{1}^{0}\mu _{{{{E}_{1}}}}^{{{{E}_{1}}}} + \\ + \;{{A}_{{02}}}\mu _{{\text{A}}}^{{{{E}_{2}}}} + {{B}_{{02}}}\mu _{{\text{B}}}^{{{{E}_{2}}}} + {{C}_{{02}}}\mu _{{\text{C}}}^{{{{E}_{2}}}} + {{D}_{{02}}}\mu _{{\text{D}}}^{{{{E}_{2}}}} + E_{2}^{0}\mu _{{{{E}_{2}}}}^{{{{E}_{2}}}}, \\ \end{gathered} $
где химический потенциал $\mu _{i}^{j}$ для бесконечно разбавленных растворов [11, 2124].
$\begin{gathered} \mu _{{\text{A}}}^{{{{E}_{1}}}} = \mu _{{\text{A}}}^{{*{{E}_{1}}}} + RT\ln [{{A}_{{01}}}{\text{/}}({{A}_{{01}}} + {{B}_{{01}}} + {{C}_{{01}}} + \\ + \;{{D}_{{01}}} + E_{1}^{0})], \\ \end{gathered} $
$\begin{gathered} \mu _{{\text{B}}}^{{{{E}_{1}}}} = \mu _{{\text{B}}}^{{*{{E}_{1}}}} + RT\ln [{{B}_{{01}}}{\text{/}}({{A}_{{01}}} + {{B}_{{01}}} + {{C}_{{01}}} + \\ + \;{{D}_{{01}}} + E_{1}^{0})], \\ \end{gathered} $
$\begin{gathered} \mu _{{\text{C}}}^{{{{E}_{1}}}} = \mu _{{\text{C}}}^{{*{{E}_{1}}}} + RT\ln [{{C}_{{01}}}{\text{/}}({{A}_{{01}}} + {{B}_{{01}}} + {{C}_{{01}}} + \\ + \;{{D}_{{01}}} + E_{1}^{0})], \\ \end{gathered} $
$\begin{gathered} \mu _{{\text{D}}}^{{{{E}_{1}}}} = \mu _{{\text{D}}}^{{*{{E}_{1}}}} + RT\ln [{{D}_{{01}}}{\text{/}}({{A}_{{01}}} + {{B}_{{01}}} + {{C}_{{01}}} + \\ + \;{{D}_{{01}}} + E_{1}^{0})], \\ \end{gathered} $
$\begin{gathered} \mu _{{{{E}_{1}}}}^{{{{E}_{1}}}} = \mu _{{{{E}_{1}}}}^{{0{{E}_{1}}}} + RT\ln [E_{1}^{0}{\text{/}}({{A}_{{01}}} + {{B}_{{01}}} + {{C}_{{01}}} + \\ + \;{{D}_{{01}}} + E_{1}^{0})], \\ \end{gathered} $
$\begin{gathered} \mu _{{\text{A}}}^{{{{E}_{2}}}} = \mu _{{\text{A}}}^{{*{{E}_{2}}}} + RT\ln [{{A}_{{02}}}{\text{/}}({{A}_{{02}}} + {{B}_{{02}}} + {{C}_{{02}}} + \\ + \;{{D}_{{02}}} + E_{2}^{0})], \\ \end{gathered} $
$\begin{gathered} \mu _{{\text{B}}}^{{{{E}_{2}}}} = \mu _{{\text{B}}}^{{*{{E}_{2}}}} + RT\ln [{{B}_{{02}}}{\text{/}}({{A}_{{02}}} + {{B}_{{02}}} + {{C}_{{02}}} + \\ + \;{{D}_{{02}}} + E_{2}^{0})], \\ \end{gathered} $
$\begin{gathered} \mu _{{\text{C}}}^{{{{E}_{2}}}} = \mu _{{\text{C}}}^{{*{{E}_{2}}}} + RT\ln [{{C}_{{02}}}{\text{/}}({{A}_{{02}}} + {{B}_{{02}}} + {{C}_{{02}}} + \\ + \;{{D}_{{02}}} + E_{2}^{0})], \\ \end{gathered} $
$\begin{gathered} \mu _{{\text{D}}}^{{{{E}_{2}}}} = \mu _{{\text{D}}}^{{*{{E}_{2}}}} + RT\ln [{{D}_{{02}}}{\text{/}}({{A}_{{02}}} + {{B}_{{02}}} + {{C}_{{02}}} + \\ + \;{{D}_{{02}}} + E_{2}^{0})], \\ \end{gathered} $
$\begin{gathered} \mu _{{{{E}_{2}}}}^{{{{E}_{2}}}} = \mu _{{{{E}_{2}}}}^{{0{{E}_{2}}}} + RT\ln [E_{2}^{0}{\text{/}}({{A}_{{02}}} + {{B}_{{02}}} + {{C}_{{02}}} + \\ + \;{{D}_{{02}}} + E_{2}^{0})], \\ \end{gathered} $
Здесь R – газовая постоянная, Т – абсолютная температура,
$\mu _{i}^{{*j}} = \mu _{i}^{{0j}} + RT\ln \gamma _{i}^{{0j}},$
где $\gamma _{i}^{{0j}}$ – коэффициенты активности при бесконечном разбавлении [11, 13, 22, 24].

Энергия Гиббса в конечном состоянии (состояние равновесия [20, с. 150])

$\begin{gathered} {{G}_{{\text{K}}}} = {{A}_{1}}\mu _{{\text{A}}}^{{{{E}_{1}}}} + {{B}_{1}}\mu _{{\text{B}}}^{{{{E}_{1}}}} + {{C}_{1}}\mu _{{\text{C}}}^{{{{E}_{1}}}} + {{D}_{1}}\mu _{{\text{D}}}^{{{{E}_{1}}}} + E_{1}^{0}\mu E_{{{{E}_{1}}}}^{{{{E}_{1}}}} + \\ + \;{{A}_{2}}\mu _{{\text{A}}}^{{{{E}_{2}}}} + {{B}_{2}}\mu _{{\text{B}}}^{{{{E}_{2}}}} + {{C}_{2}}\mu _{{\text{C}}}^{{{{E}_{2}}}} + {{D}_{2}}\mu _{{\text{D}}}^{{{{E}_{2}}}} + E_{2}^{0}\mu E_{{{{E}_{2}}}}^{{{{E}_{2}}}}. \\ \end{gathered} $
Здесь
$\mu _{{\text{A}}}^{{{{E}_{1}}}} = \mu _{{\text{A}}}^{{*{{E}_{1}}}} + RT\ln [{{A}_{1}}{\text{/}}({{A}_{1}} + {{B}_{1}} + {{C}_{1}} + {{D}_{1}} + E_{1}^{0})],$
$\mu _{{\text{B}}}^{{{{E}_{1}}}} = \mu _{{\text{B}}}^{{*{{E}_{1}}}} + RT\ln [{{B}_{1}}{\text{/}}({{A}_{1}} + {{B}_{1}} + {{C}_{1}} + {{D}_{1}} + E_{1}^{0})],$
$\mu _{{\text{C}}}^{{{{E}_{1}}}} = \mu _{{\text{C}}}^{{*{{E}_{1}}}} + RT\ln [{{C}_{1}}{\text{/}}({{A}_{1}} + {{B}_{1}} + {{C}_{1}} + {{D}_{1}} + E_{1}^{0})],$
$\mu _{{\text{D}}}^{{{{E}_{1}}}} = \mu _{{\text{D}}}^{{*{{E}_{1}}}} + RT\ln [{{D}_{1}}{\text{/}}({{A}_{1}} + {{B}_{1}} + {{C}_{1}} + {{D}_{1}} + E_{1}^{0})],$
$\mu _{{{{E}_{1}}}}^{{{{E}_{1}}}} = \mu _{{{{E}_{1}}}}^{{0{{E}_{1}}}} + RT\ln [E_{1}^{0}{\text{/}}({{A}_{1}} + {{B}_{1}} + {{C}_{1}} + {{D}_{1}} + E_{1}^{0})],$
$\mu _{{\text{A}}}^{{{{E}_{2}}}} = \mu _{{\text{A}}}^{{*{{E}_{2}}}} + RT\ln [{{A}_{2}}{\text{/}}({{A}_{2}} + {{B}_{2}} + {{C}_{2}} + {{D}_{2}} + E_{2}^{0})],$
$\mu _{{\text{B}}}^{{{{E}_{2}}}} = \mu _{{\text{B}}}^{{*{{E}_{2}}}} + RT\ln [{{B}_{2}}{\text{/}}({{A}_{2}} + {{B}_{2}} + {{C}_{2}} + {{D}_{2}} + E_{2}^{0})],$
$\mu _{{\text{C}}}^{{{{E}_{2}}}} = \mu _{{\text{C}}}^{{*{{E}_{2}}}} + RT\ln [{{C}_{2}}{\text{/}}({{A}_{2}} + {{B}_{2}} + {{C}_{2}} + {{D}_{2}} + E_{2}^{0})],$
$\mu _{{\text{D}}}^{{{{E}_{2}}}} = \mu _{{\text{D}}}^{{*{{E}_{2}}}} + RT\ln [{{D}_{2}}{\text{/}}({{A}_{2}} + {{B}_{2}} + {{C}_{2}} + {{D}_{2}} + E_{2}^{0})],$
$\mu _{{{{E}_{2}}}}^{{{{E}_{2}}}} = \mu _{{{{E}_{2}}}}^{{0{{E}_{2}}}} + RT\ln [E_{2}^{0}{\text{/}}({{A}_{2}} + {{B}_{2}} + {{C}_{2}} + {{D}_{2}} + E_{2}^{0})],$
Из равенства химических потенциалов при равновесии по данным [11, 13]
$\begin{gathered} \mu _{{\text{A}}}^{{*{{E}_{2}}}}\, - \,\mu _{{\text{A}}}^{{*{{E}_{1}}}}\, = \,RT\ln \frac{{{{A}_{1}}{\text{/}}({{A}_{1}}\, + \,{{B}_{1}}\, + \,{{C}_{1}}\, + \,{{D}_{1}}\, + \,E_{1}^{0})}}{{{{A}_{2}}{\text{/}}({{A}_{2}}\, + \,{{B}_{2}}\, + \,{{C}_{2}}\, + \,{{D}_{2}}\, + \,E_{2}^{0})}} = \\ = RT\ln (\gamma _{{\text{A}}}^{{0{{E}_{2}}}}{\text{/}}\gamma _{{\text{A}}}^{{0{{E}_{1}}}}) = RT\ln K_{{\text{A}}}^{*}, \\ \end{gathered} $
$\begin{gathered} \mu _{{\text{B}}}^{{*{{E}_{2}}}}\, - \,\mu _{{\text{B}}}^{{*{{E}_{1}}}}\, = \,RT\ln \frac{{{{B}_{1}}{\text{/}}({{A}_{1}}\, + \,{{B}_{1}}\, + \,{{C}_{1}}\, + \,{{D}_{1}}\, + \,E_{1}^{0})}}{{{{B}_{2}}{\text{/}}({{A}_{2}}\, + \,{{B}_{2}}\, + \,{{C}_{2}}\, + \,{{D}_{2}}\, + \,E_{2}^{0})}} = \\ = RT\ln (\gamma _{{\text{B}}}^{{0{{E}_{2}}}}{\text{/}}\gamma _{{\text{B}}}^{{0{{E}_{1}}}}) = RT\ln K_{{\text{B}}}^{*}, \\ \end{gathered} $
$\begin{gathered} \mu _{{\text{C}}}^{{*{{E}_{2}}}}\, - \,\mu _{{\text{C}}}^{{*{{E}_{1}}}}\, = \,RT\ln \frac{{{{C}_{1}}{\text{/}}({{A}_{1}}\, + \,{{B}_{1}}\, + \,{{C}_{1}}\, + \,{{D}_{1}}\, + \,E_{1}^{0})}}{{{{C}_{2}}{\text{/}}({{A}_{2}}\, + \,{{B}_{2}}\, + \,{{C}_{2}}\, + \,{{D}_{2}}\, + \,E_{2}^{0})}} = \\ = RT\ln (\gamma _{{\text{C}}}^{{0{{E}_{2}}}}{\text{/}}\gamma _{{\text{C}}}^{{0{{E}_{1}}}}) = RT\ln K_{{\text{C}}}^{*}, \\ \end{gathered} $
$\begin{gathered} \mu _{{\text{D}}}^{{*{{E}_{2}}}}\, - \,\mu _{{\text{D}}}^{{*{{E}_{1}}}}\, = \,RT\ln \frac{{{{D}_{1}}{\text{/}}({{A}_{1}}\, + \,{{B}_{1}}\, + \,{{C}_{1}}\, + \,{{D}_{1}}\, + \,E_{1}^{0})}}{{{{D}_{2}}{\text{/}}({{A}_{2}}\, + \,{{B}_{2}}\, + \,{{C}_{2}}\, + \,{{D}_{2}}\, + \,E_{2}^{0})}} = \\ = RT\ln (\gamma _{{\text{D}}}^{{0{{E}_{2}}}}{\text{/}}\gamma _{{\text{D}}}^{{0{{E}_{1}}}}) = RT\ln K_{{\text{D}}}^{*}. \\ \end{gathered} $
Изменение энергии Гиббса (конечная величина [4])
$\begin{gathered} --\Delta G = --({{G}_{{\text{K}}}}--{{G}_{{\text{Н}}}}) = \Delta {{n}_{{\text{B}}}}RT{\text{ln}}\frac{{{{B}_{1}}{\text{/}}E_{1}^{0}}}{{{{B}_{2}}{\text{/}}E_{2}^{0}}}-- \\ - \;\Delta {{n}_{{\text{A}}}}RT{\text{ln}}\frac{{{{A}_{1}}{\text{/}}E_{1}^{0}}}{{{{A}_{2}}{\text{/}}E_{2}^{0}}} + \Delta {{n}_{{\text{D}}}}RT{\text{ln}}\frac{{{{D}_{1}}{\text{/}}E_{1}^{0}}}{{{{D}_{2}}{\text{/}}E_{2}^{0}}}-- \\ \end{gathered} $
$\begin{gathered} - \;\Delta {{n}_{{\text{C}}}}RT{\text{ln}}\frac{{{{C}_{1}}{\text{/}}E_{1}^{0}}}{{{{C}_{2}}{\text{/}}E_{2}^{0}}} + RT(\Delta {{n}_{{\text{B}}}}--\Delta {{n}_{{\text{A}}}} + \Delta {{n}_{{\text{D}}}} - \\ --\;\Delta {{n}_{{\text{C}}}})\left( {\frac{{{{A}_{2}}\, + \,{{B}_{2}}\, + \,{{C}_{2}}\, + \,{{D}_{2}}}}{{E_{2}^{0}}} - \frac{{{{A}_{1}}\, + \,{{B}_{1}}\, + \,{{C}_{1}}\, + \,{{D}_{1}}}}{{E_{1}^{0}}}} \right) - \\ \end{gathered} $
$\begin{gathered} - \;RT[{{B}_{1}}\ln {{B}_{1}}{\text{/}}E_{1}^{0} + {{B}_{2}}\ln {{B}_{2}}{\text{/}}E_{2}^{0} + {{A}_{1}}\ln {{A}_{1}}{\text{/}}E_{1}^{0} + \\ + \;{{A}_{2}}\ln {{A}_{2}}{\text{/}}E_{2}^{0} + {{D}_{1}}\ln {{D}_{1}}{\text{/}}E_{1}^{0} + {{D}_{2}}\ln {{D}_{2}}{\text{/}}E_{2}^{0} + \\ \end{gathered} $
$\begin{gathered} + \;{{C}_{1}}\ln {{C}_{1}}{\text{/}}E_{1}^{0} + {{C}_{2}}\ln {{C}_{2}}{\text{/}}E_{2}^{0}--{{B}_{{01}}}\ln {{B}_{{01}}}{\text{/}}E_{1}^{0} - \\ --\;{{B}_{{02}}}\ln {{B}_{{02}}}{\text{/}}E_{2}^{0}--{{A}_{{01}}}\ln {{A}_{{01}}}{\text{/}}E_{1}^{0}--{{A}_{{02}}}\ln {{A}_{{02}}}{\text{/}}E_{2}^{0} - \\ \end{gathered} $
$\begin{gathered} --\;{{D}_{{01}}}\ln {{D}_{{01}}}{\text{/}}E_{1}^{0}--{{D}_{{02}}}\ln {{D}_{{02}}}{\text{/}}E_{2}^{0}--{{C}_{{01}}}\ln {{C}_{{01}}}{\text{/}}E_{1}^{0} - \\ --\;{{C}_{{02}}}\ln {{C}_{{02}}}{\text{/}}E_{2}^{0}]--RT(\Delta {{n}_{{\text{B}}}}--\Delta {{n}_{{\text{A}}}} + \Delta {{n}_{{\text{D}}}} - \\ \end{gathered} $
$\begin{gathered} --\;\Delta {{n}_{{\text{C}}}})\left( {\frac{{{{A}_{2}}\, + \,{{B}_{2}}\, + \,{{C}_{2}}\, + \,{{D}_{2}}}}{{E_{2}^{0}}} - \frac{{{{A}_{1}}\, + \,{{B}_{1}}\, + \,{{C}_{1}}\, + \,{{D}_{1}}}}{{E_{1}^{0}}}} \right) = \\ = \Delta {{n}_{{\text{B}}}}RT\ln K_{{\text{B}}}^{*}--\Delta {{n}_{{\text{A}}}}RT\ln K_{{\text{A}}}^{*} - T\Delta {{S}_{{{\text{AB}}}}} + \\ + \;\Delta {{n}_{{\text{D}}}}RT\ln K_{{\text{D}}}^{*}--\Delta {{n}_{{\text{C}}}}RT\ln K_{{\text{C}}}^{*}--T\Delta {{S}_{{{\text{CD}}}}}, \\ \end{gathered} $
где [23]
$\begin{gathered} \Delta {{S}_{{{\text{AB}}}}} = R\left\{ {{{B}_{0}}\left[ {\frac{{{{B}_{1}}}}{{{{B}_{0}}}}\ln \frac{{{{B}_{1}}{\text{/}}{{B}_{0}}}}{{E_{1}^{0}{\text{/}}(E_{1}^{0} + E_{2}^{0})}}} \right.} \right. + \\ + \;\left. {\frac{{{{B}_{2}}}}{{{{B}_{0}}}}\ln \frac{{{{B}_{2}}{\text{/}}{{B}_{0}}}}{{E_{2}^{0}{\text{/}}(E_{1}^{0} + E_{2}^{0})}}} \right] + {{A}_{0}}\left[ {\frac{{{{A}_{1}}}}{{{{A}_{0}}}}\ln \frac{{{{A}_{1}}{\text{/}}{{A}_{0}}}}{{E_{1}^{0}{\text{/}}(E_{1}^{0} + E_{2}^{0})}}} \right. + \\ + \;\left. {\left. {\frac{{{{A}_{2}}}}{{{{A}_{0}}}}\ln \frac{{{{A}_{2}}{\text{/}}{{A}_{0}}}}{{E_{2}^{0}{\text{/}}(E_{1}^{0} + E_{2}^{0})}}} \right]} \right\}, \\ \end{gathered} $
$\begin{gathered} \Delta {{S}_{{{\text{CD}}}}} = R\left\{ {{{D}_{0}}\left[ {\frac{{{{D}_{1}}}}{{{{D}_{0}}}}\ln \frac{{{{D}_{1}}{\text{/}}{{D}_{0}}}}{{E_{1}^{0}{\text{/}}(E_{1}^{0} + E_{2}^{0})}}} \right.} \right. + \\ + \;\left. {\frac{{{{D}_{2}}}}{{{{D}_{0}}}}\ln \frac{{{{D}_{2}}{\text{/}}{{D}_{0}}}}{{E_{2}^{0}{\text{/}}(E_{1}^{0} + E_{2}^{0})}}} \right] + {{C}_{0}}\left[ {\frac{{{{C}_{1}}}}{{{{C}_{0}}}}\ln \frac{{{{C}_{1}}{\text{/}}{{C}_{0}}}}{{E_{1}^{0}{\text{/}}(E_{1}^{0} + E_{2}^{0})}}} \right. + \\ + \;\left. {\left. {\frac{{{{C}_{2}}}}{{{{C}_{0}}}}\ln \frac{{{{C}_{2}}{\text{/}}{{C}_{0}}}}{{E_{2}^{0}{\text{/}}(E_{1}^{0} + E_{2}^{0})}}} \right]} \right\}. \\ \end{gathered} $
Для членов одного гомологического ряда [17]
$K_{{\text{B}}}^{*} = K_{{\text{F}}}^{*}\alpha _{{{\text{AB}}}}^{{1/2}},\quad K_{{\text{A}}}^{*} = K_{{\text{F}}}^{*}{\text{/}}\alpha _{{{\text{AB}}}}^{{1/2}},$
$K_{{\text{D}}}^{*} = K_{{\text{F}}}^{*}\alpha _{{{\text{CD}}}}^{{1/2}},\quad K_{{\text{C}}}^{*} = K_{{\text{F}}}^{*}{\text{/}}\alpha _{{{\text{CD}}}}^{{1/2}}$
и характерна линейная зависимость логарифма коэффициента активности γ0 от числа атомов углерода [24].

Записанные соотношения для ∆G можно получить по аналогии с данными [12], используя предыдущие уравнения и применяя для предельно разбавленных растворов формулы:

$\begin{gathered} \ln [A{\text{/}}(A + B + C + D + {{E}^{0}})] = \\ = \ln A{\text{/}}{{E}^{0}}--\frac{{A + B + C + D}}{{{{E}^{0}}}}, \\ \end{gathered} $
$\begin{gathered} \ln [B{\text{/}}(A + B + C + D + {{E}^{0}})] = \\ = \ln B{\text{/}}{{E}^{0}}--\frac{{A + B + C + D}}{{{{E}^{0}}}}, \\ \end{gathered} $
$\begin{gathered} \ln [C{\text{/}}(A + B + C + D + {{E}^{0}})] = \\ = \ln C{\text{/}}{{E}^{0}}--\frac{{A + B + C + D}}{{{{E}^{0}}}}, \\ \end{gathered} $
$\begin{gathered} \ln [D{\text{/}}(A + B + C + D + {{E}^{0}})] = \\ = \ln D{\text{/}}{{E}^{0}}--\frac{{A + B + C + D}}{{{{E}^{0}}}}, \\ \end{gathered} $
$\begin{gathered} \ln [{{E}^{0}}{\text{/}}(A + B + C + D + {{E}^{0}})] = \\ = \ln \left[ {1{\text{/}}\left( {1 + \frac{{A + B + C + D}}{{{{E}^{0}}}}} \right)} \right] = \\ --\frac{{A + B + C + D}}{{{{E}^{0}}}}. \\ \end{gathered} $
При отдельном разделении пар AB и CD имеем [12]:
$\begin{gathered} - \Delta {{G}_{{{\text{AB}}}}} = \Delta {{n}_{{\text{B}}}}RT\ln \frac{{{{B}_{1}}{\text{/}}E_{1}^{0}}}{{{{B}_{2}}{\text{/}}E_{2}^{0}}} - \Delta {{n}_{{\text{A}}}}RT\ln \frac{{{{A}_{1}}{\text{/}}E_{1}^{0}}}{{{{A}_{2}}{\text{/}}E_{2}^{0}}} + \\ + \;RT(\Delta {{n}_{{\text{B}}}} - \Delta {{n}_{{\text{A}}}})\left[ {\frac{{{{A}_{2}} + {{B}_{2}}}}{{E_{2}^{0}}} - \frac{{{{A}_{1}} + {{B}_{1}}}}{{E_{1}^{0}}}} \right] - \\ \end{gathered} $
$\begin{gathered} --RT[{{B}_{1}}\ln {{B}_{1}}{\text{/}}E_{1}^{0} + {{B}_{2}}\ln B{\text{/}}E_{2}^{0} + {{A}_{1}}\ln {{A}_{1}}{\text{/}}E_{1}^{0} + \\ + \;{{A}_{2}}\ln {{A}_{2}}{\text{/}}E_{2}^{0}--{{B}_{{01}}}\ln {{B}_{{01}}}{\text{/}}E_{1}^{0}--{{B}_{{02}}}\ln {{B}_{{02}}}{\text{/}}E_{2}^{0} - \\ --\;{{A}_{{01}}}\ln {{A}_{{01}}}{\text{/}}E_{1}^{0}--{{A}_{{02}}}\ln {{A}_{{02}}}{\text{/}}E_{2}^{0}] - \\ \end{gathered} $
$\begin{gathered} --\;RT\left( {\Delta {{n}_{{\text{B}}}}--\Delta {{n}_{{\text{A}}}}} \right)\left[ {\frac{{{{A}_{2}} + {{B}_{2}}}}{{E_{2}^{0}}} - \frac{{{{A}_{1}} + {{B}_{1}}}}{{E_{1}^{0}}}} \right] = \\ = \Delta {{n}_{{\text{B}}}}RT\ln K_{{\text{B}}}^{*}--\Delta {{n}_{{\text{A}}}}RT\ln K_{{\text{A}}}^{*}--T\Delta {{S}_{{{\text{AB}}}}}, \\ \end{gathered} $
$ - \Delta {{G}_{{{\text{CD}}}}} = \Delta {{n}_{{\text{D}}}}RT\ln K_{{\text{D}}}^{*}--\Delta {{n}_{{\text{C}}}}RT\ln K_{{\text{C}}}^{*}--T\Delta {{S}_{{{\text{CD}}}}},$
Убыль энергии Гиббса, отвечающая оптимальным условиям проведения процесса разделения, следует, соответственно [12], из уравнения
$\begin{gathered} - \Delta G = \frac{{{{q}_{{{\text{AB}}}}}}}{2}RT\frac{{\alpha _{{{\text{AB}}}}^{{1/2}} - 1}}{{\alpha _{{{\text{AB}}}}^{{1/2}} + 1}}\ln \alpha _{{{\text{AB}}}}^{{1/2}} - {{q}_{{{\text{AB}}}}}RT \times \\ \times \;\left[ {\frac{{\alpha _{{{\text{AB}}}}^{{1/2}}}}{{\alpha _{{{\text{AB}}}}^{{1/2}} + 1}}\ln \frac{{2\alpha _{{{\text{AB}}}}^{{1/2}}}}{{\alpha _{{{\text{AB}}}}^{{1/2}} + 1}} + \frac{1}{{\alpha _{{{\text{AB}}}}^{{1/2}} + 1}}\ln \frac{2}{{\alpha _{{{\text{AB}}}}^{{1/2}} + 1}}} \right] + \\ \end{gathered} $
$\begin{gathered} + \;\frac{{{{q}_{{{\text{CD}}}}}}}{2}RT\frac{{\alpha _{{{\text{CD}}}}^{{1/2}} - 1}}{{\alpha _{{{\text{CD}}}}^{{1/2}} + 1}}\ln \alpha _{{{\text{CD}}}}^{{1/2}} - {{q}_{{{\text{CD}}}}}RT \times \\ \times \;\left[ {\frac{{\alpha _{{{\text{CD}}}}^{{1/2}}}}{{\alpha _{{{\text{CD}}}}^{{1/2}} + 1}}\ln \frac{{2\alpha _{{{\text{CD}}}}^{{1/2}}}}{{\alpha _{{{\text{CD}}}}^{{1/2}} + 1}} + \frac{1}{{\alpha _{{{\text{CD}}}}^{{1/2}} + 1}}\ln \frac{2}{{\alpha _{{{\text{CD}}}}^{{1/2}} + 1}}} \right] = \\ - (\Delta {{G}_{{{\text{AB}}}}} + \Delta {{G}_{{{\text{CD}}}}}), \\ \end{gathered} $
где qAB = A0+ B0, qCD = C0+ D0,
$\begin{gathered} \left[ {\frac{{{{\alpha }^{{1/2}}}}}{{{{\alpha }^{{1/2}}} + 1}}\ln \frac{{2{{\alpha }^{{1/2}}}}}{{{{\alpha }^{{1/2}}} + 1}} + \frac{1}{{{{\alpha }^{{1/2}}} + 1}}\ln \frac{2}{{{{\alpha }^{{1/2}}} + 1}}} \right]\; \times \\ \times \;{{\left[ {\frac{{{{\alpha }^{{1/2}}} - 1}}{{{{\alpha }^{{1/2}}} + 1}}\ln {{\alpha }^{{1/2}}}} \right]}^{{ - 1}}} = f(\alpha ). \\ \end{gathered} $
При α < 1.5 имеем f(α) ≈ 0.25 [12], поэтому для тонкого разделения [1], когда α 1 ≪ 1, можно принять
$\begin{gathered} - \Delta G = \frac{{RT}}{2}\left[ {{{q}_{{{\text{AB}}}}}\frac{{\alpha _{{{\text{AB}}}}^{{1/2}} - 1}}{{\alpha _{{{\text{AB}}}}^{{1/2}} + 1}}\ln \alpha _{{{\text{AB}}}}^{{1/2}}} \right. - \\ - \;\left. {{{q}_{{{\text{CD}}}}}\frac{{\alpha _{{{\text{CD}}}}^{{1/2}} - 1}}{{\alpha _{{{\text{CD}}}}^{{1/2}} + 1}}\ln \alpha _{{{\text{CD}}}}^{{1/2}}} \right] = \\ \end{gathered} $
$\begin{gathered} = \frac{{RT}}{2}[{{q}_{{{\text{AB}}}}}({{\alpha }_{{{\text{AB}}}}}--1){\text{/}}8 + {{q}_{{{\text{CD}}}}}({{\alpha }_{{{\text{CD}}}}}--1){\text{/}}8] \\ = \frac{{RT}}{2}({{W}_{{{\text{AB}}}}} + {{W}_{{{\text{CD}}}}}). \\ \end{gathered} $
Здесь W характеризует приращение ценности смеси и максимальную разделительную способность системы при отдельном разделении пар AB и CD [1, 5, 6, 25]. Название “работа разделения” для W = q 1)/ 8 предложено Юри в 1939 г. [1, 5]. Таким образом, установлена связь между изменением энергии Гиббса экстракционной системы и фундаментальными понятиями общей теории процессов разделения. Эти понятия имеют большое значение [6], так как позволяют дать сравнительную оценку эффективности процессов разделения [2, 5, 7].

Из проведенного расчета и полученных уравнений следует, что ∆G = GK– GН < 0. Если ∆G < 0, то с точки зрения химической термодинамики процесс является самопроизвольным и способным [8, 9] совершить обратимую работу WT = = ‒∆G, где ∆G – изменение энергии Гиббса, вычисленное при полной неравновесности процесса. Эта функция обладает свойством термодинамической аддитивности и в рассмотренном примере разделения многокомпонентной смеси, состоящей из выбранных пар AB и CD, представляет собой сумму работ WT,AB = GAB и WT,CD = = GCD отдельного разделения указанных пар [7]. Таким образом, подтверждено предполагаемое ранее [7] свойство термодинамической аддитивности работы разделения.

ОСНОВНЫЕ ВЫВОДЫ

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

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

  1. Розен А.М. Теория разделения изотопов в колоннах. М.: Атомиздат, 1960. 440 с.

  2. Кузнецов П.Г., Соколов В.А., Седин И.П. // Разделение и анализ углеводородных газов. М.: Изд-во АН СССР, 1963. С. 114.

  3. Жуховицкий А.А., Туркельтауб Н.М. Газовая хроматография. М.: Гостоптехиздат, 1962. С. 7, 9, 10.

  4. Кириллин В.А., Шейндлин А.Е., Шпильрайн Э.Э. Термодинамика растворов. М.: Энергия, 1980. С. 225.

  5. Изотопы: свойства, получение, применение / Под ред. В.Ю. Баранова. М.: ИздАТ, 2000. 703 с.

  6. Сулаберидзе Г.А., Борисевич В.Д. // Журн. физ. химии. 2008. Т. 82. № 7. С. 1392.

  7. Чижков В.П., Бойцов В.Н. // Там же. 2019. Т. 93. № 4. С. 632.

  8. Еремин Е.И. Основы химической термодинамики. М.: Высшая школа, 1974. 340 с.

  9. Лейтес И.П., Сосна М.Х., Семенов В.П. Теория и практика химической энерготехнологии. М.: Химия, 1988. 280 с.

  10. Эверет Д. Введение в химическую термодинамику. М.: Изд-во иностр. лит., 1963. 299 с.

  11. Жидкостная колоночная хроматография / Под ред. З. Дейла, К. Мацека, Я. Янака. М.: Мир, 1978. Т. 1. 554 с.

  12. Чижков В.П., Бойцов В.Н. // Журн. физ. химии. 2020. Т. 94. № 2. С. 301.

  13. Трейбал Р. Жидкостная экстракция. М.: Химия, 1966, 722 с.

  14. Чижков В.П., Бойцов В.Н. // Журн. физ. химии. 2012. Т. 86. № 1. С. 139.

  15. Чижков В.П., Бойцов В.Н., Демин А.В. // Там же. 2011. Т. 85. № 3. С. 565.

  16. Бенедикт М., Пигфорд Т. Химическая технология ядерных материалов. М.: Атомиздат, 1960. 560 с.

  17. Чижков В.П., Бойцов В.Н. // Журн. физ. химии. 2014. Т. 88. № 9. С. 1440.

  18. Чижков В.П., Варивончик Э.А., Горячко Ю.В. и др. // Там же. 1991. Т. 65. № 1. С. 1.

  19. Бреслер С.Е. Радиоактивные элементы. М.: Гос. изд-во техн.-теор. лит., 1957. 545 с.

  20. Даниэлс Ф., Олберти Р. Физическая химия. М.: Мир, 1978. 645 с.

  21. Кричевский И.Р. Понятия и основы термодинамики. М.: Химия, 1970, с. 358, 360

  22. Герасимов Я.И., Дервинг В.П., Еремин Е.Н. и др. Курс физической химии. Т. 1. М.: ГНТИХЛ. 1963, с. 216.

  23. Чижков В.П., Бойцов В.Н. // Журн. физ. химии. 2005. Т. 79. № 3. С. 391.

  24. Руководство по газовой хроматографии / Под ред. Э. Лейбница, Х.Г. Штруппе. М.: Мир, 1988. Т. 1. 479 с.

  25. Чижков В.П., Бойцов В.Н. // Журн. физ. химии. 2015. Т. 89. № 4. С. 700.

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