Физикохимия поверхности и защита материалов, 2021, T. 57, № 6, стр. 574-580

К условиям расслоения двухкомпонентного адсорбционного слоя при адсорбции по изотерме Фрумкина–Дамаскина

Э. М. Подгаецкий *

Федеральное государственное бюджетное учреждение науки, Институт прикладной механики Российской академии наук
119991 Ленинский просп., 32а, Москва, Россия

* E-mail: Podgaetsky@mail.ru

Поступила в редакцию 03.03.2021
После доработки 12.07.2021
Принята к публикации 13.07.2021

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

Аннотация

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

Ключевые слова: совместная адсорбция, изотерма Фрумкина–Дамаскина, условия расслоения

ВВЕДЕНИЕ

В [1] рассматривался равновесный двухкомпонентный адсорбционный слой – частиц “а” с поверхностной концентрацией ${{\tilde {\Gamma }}_{1}}$ и частиц “b” с поверхностной концентрацией ${{{{\tilde {Г}}}}_{2}}$, образованный из жидкого раствора таких же частиц “а” и “b” на электронейтральной поверхности твердого тела, но с учетом ее деформации $\vartheta $. Выведенные в [1] уравнения изотерм совместной адсорбции базируются на уравнениях совместной адсорбции на твердой недеформированной поверхности

(1)
${{В}_{1}}{{с}_{1}}\left| {_{{\vartheta = 0}} = {{А}_{1}}({{{\text{Г}}}_{1}}{\text{,}}{{{\text{Г}}}_{2}})} \right.,$
(2)
${{В}_{2}}{{с}_{2}}\left| {_{{\vartheta = 0}}} \right. = {{А}_{2}}({{{\text{Г}}}_{1}},{{{\text{Г}}}_{2}}),$
где ${{В}_{1}},\,\,{{В}_{2}}$ – безразмерные константы, ${{{\text{Г}}}_{1}},\,\,{{{\text{Г}}}_{2}}$ – безразмерные значения ${{{{\tilde {Г}}}}_{1}}$, ${{{{\tilde {Г}}}}_{2}}$ частиц “а” и “b”, адсорбируемых единичной площадью твердой поверхности
(3)
${{{\text{Г}}}_{1}} = \frac{{{{{{{\tilde {Г}}}}}_{1}}}}{{{\text{Г*}}}},\,\,\,\,{{{\text{Г}}}_{2}} = \frac{{{{{{{\tilde {Г}}}}}_{2}}}}{{{\text{Г*}}}},$
а ${\text{Г*}} = {{{\text{Г}}}_{\infty }}$ – максимальное суммарное значение ${{{{\tilde {Г}}}}_{1}}$ и ${{{{\tilde {Г}}}}_{2}}$, адсорбируемое единичной площадью твердой поверхности, ${{с}_{1}},{{с}_{2}}$ – соответствующие безразмерные объемные концентрации.

Система уравнений (1), (2) неизбежно с определенного этапа ведет к поиску условий, при которых обеспечивается единственность значений ${{{\text{Г}}}_{1}},{{{\text{Г}}}_{2}}$ при заданных ${{с}_{1}},{{с}_{2}}$, т.е. решений уравнений (1), (2) относительно ${{{\text{Г}}}_{1}},{{{\text{Г}}}_{2}}$. Такое условие при однокомпонентной адсорбции по изотерме, например Фрумкина, [2]

(4)
${{В}_{1}}{{с}_{1}}_{{\left| \begin{subarray}{l} \vartheta = 0 \\ {{{\text{Г}}}_{2}} = 0 \end{subarray} \right.}} = \frac{{{{{\text{Г}}}_{1}}\exp ( - 2{{a}_{1}}{{{\text{Г}}}_{1}})}}{{1 - {{{\text{Г}}}_{1}}}} = {{А}_{F}}({{{\text{Г}}}_{1}}),$
где ${{а}_{1}}$ – аттракционная постоянная Фрумкина; имеет простой вид

${{а}_{1}} < 2.$

Значения ${{а}_{1}} > 2$ традиционно трактуются как возможность образования при определенных значениях объемной концентрации ${{с}_{1}}$ двух поверхностных фаз [3] по аналогии с уравнением Ван-дер-Ваальса. Двухфазность в адсорбционном слое также может быть вызвана физико-химически одинаковыми молекулами “а” и “b”, но с разной ориентацией к твердой поверхности при их адсорбции [4]. В этом случае описание такой ситуации математически не отличается от общей постановки в системе (1), (2).

ИСХОДНЫЕ УРАВНЕНИЯ

В [5, 6] исследование уравнений (1), (2) проводилось в частном случае совместной адсорбции – по изотерме Фрумкина–Дамаскина [7]

(5)
${{A}_{1}} = \frac{{{{{\text{Г}}}_{1}}\exp ( - 2{{a}_{1}}{{{\text{Г}}}_{1}} - 2{{a}_{3}}{{{\text{Г}}}_{2}})}}{{1 - {{{\text{Г}}}_{1}} - {{{\text{Г}}}_{2}}}},$
(6)
${{A}_{2}} = \frac{{{{{\text{Г}}}_{2}}\exp ( - 2{{a}_{2}}{{{\text{Г}}}_{2}} - 2{{a}_{3}}{{{\text{Г}}}_{1}})}}{{1 - {{{\text{Г}}}_{1}} - {{{\text{Г}}}_{2}}}},$
в области $\Delta $ значений переменных ${{{\text{Г}}}_{1}},{{{\text{Г}}}_{2}}$
(7)
$\begin{gathered} \Delta \,:\,0 < {{{\text{Г}}}_{1}} + {{{\text{Г}}}_{2}} < 1, \hfill \\ \,\,\,\,\,0 < {{{\text{Г}}}_{1}},\,\,0 < {{{\text{Г}}}_{2}}, \hfill \\ \end{gathered} $
используя сформулированный в [5] общий критерий единственности решения системы (5), (6) относительно ${{{\text{Г}}}_{1}},{{{\text{Г}}}_{2}}$
(8)
$J\left( {\frac{{{{A}_{1}},{{A}_{2}}}}{{{{{\text{Г}}}_{1}},{{{\text{Г}}}_{2}}}}} \right) \ne 0,\,\,\,\,({{{\text{Г}}}_{1}},{{{\text{Г}}}_{2}}) \in \Delta ,$
где J – якобиан преобразований ${{А}_{1}}({{{\text{Г}}}_{1}},{{{\text{Г}}}_{2}}),$ ${{А}_{2}}({{{\text{Г}}}_{1}},{{{\text{Г}}}_{2}})$
(9)
$J \equiv \frac{{\partial {{А}_{1}}}}{{\partial {{{\text{Г}}}_{1}}}}\frac{{\partial {{А}_{2}}}}{{\partial {{{\text{Г}}}_{2}}}} - \frac{{\partial {{А}_{1}}}}{{\partial {{{\text{Г}}}_{2}}}}\frac{{\partial {{А}_{2}}}}{{\partial {{{\text{Г}}}_{1}}}},$
но в частном случае значений ${{а}_{1}},{{а}_{2}},{{а}_{3}}$, когда ${{а}_{1}} = {{а}_{2}}$.

В последующем анализе это ограничение не обязательно и будут выводиться достаточные условия по параметрам ${{а}_{1}},{{а}_{2}},{{а}_{3}}$, при которых условие (8) выполняется.

Подставляя функции А1, А2 – (5),(6) – в левую часть (8), получим выражение функции J1, Г2)

(10)
$J = \frac{{{{e}^{g}}{{J}_{0}}}}{{{{{(1 - {\text{Г}})}}^{4}}}},$
где

(11)
$\begin{gathered} {{J}_{0}} \equiv [{{{\text{Г}}}_{1}} + (1 - 2{{а}_{1}}{{{\text{Г}}}_{1}})(1 - {\text{Г}})] \times \\ \times \,\,[{{{\text{Г}}}_{2}} + (1 - 2{{а}_{2}}{{{\text{Г}}}_{2}})(1 - {\text{Г}})] - {{{\text{Г}}}_{1}}{{{\text{Г}}}_{2}}{{[1 - 2{{а}_{3}}(1 - {\text{Г}})]}^{2}}, \\ g \equiv - 2{{а}_{1}}{{{\text{Г}}}_{1}} - 2{{а}_{2}}{{{\text{Г}}}_{2}} - 2{{а}_{3}}{{{\text{Г}}}_{2}} - 2{{а}_{3}}{{{\text{Г}}}_{1}}, \\ {\text{Г}} \equiv {{{\text{Г}}}_{1}} + {{{\text{Г}}}_{2}}. \\ \end{gathered} $

Знакоопределяющий множитель J0 в (11) можно представить в виде

(12)
$\begin{gathered} {{J}_{0}} = {{{\text{Г}}}_{1}}{{{\text{Г}}}_{2}} + {{{\text{Г}}}_{1}}(1 - 2{{а}_{2}}{{{\text{Г}}}_{2}})(1 - {\text{Г}}) + \\ + \,\,{{{\text{Г}}}_{2}}(1 - 2{{а}_{1}}{{{\text{Г}}}_{1}})(1 - {\text{Г}}) + {{(1 - {\text{Г}})}^{2}}(1 - 2{{а}_{1}}{{{\text{Г}}}_{1}}) \times \\ \times \,\,(1 - 2{{а}_{2}}{{{\text{Г}}}_{2}}) - {{{\text{Г}}}_{1}}{{{\text{Г}}}_{2}} + 4{{а}_{3}}{{{\text{Г}}}_{1}}{{{\text{Г}}}_{2}}(1 - {\text{Г}}) - \\ - \,\,4{{{\text{Г}}}_{1}}{{{\text{Г}}}_{2}}а_{3}^{2}{{(1 - {\text{Г}})}^{2}} = (1 - {\text{Г}})[{{{\text{Г}}}_{1}}(1 - 2{{а}_{2}}{{{\text{Г}}}_{2}}) + \\ + \,\,{{{\text{Г}}}_{2}}(1 - 2{{а}_{1}}{{{\text{Г}}}_{1}}) + (1 - {\text{Г}})(1 - 2{{а}_{1}}{{{\text{Г}}}_{1}})(1 - 2{{а}_{2}}{{{\text{Г}}}_{2}}) + \\ + \,\,4{{а}_{3}}{{{\text{Г}}}_{1}}{{{\text{Г}}}_{2}} - 4а_{3}^{2}{{{\text{Г}}}_{1}}{{{\text{Г}}}_{2}}(1 - {\text{Г}})] \equiv (1 - {\text{Г}}){{J}_{s}}, \\ \end{gathered} $
где

(12a)
$\begin{gathered} {{J}_{s}} \equiv {{G}_{0}} + {{G}_{1}}, \\ {{G}_{0}} \equiv {{{\text{Г}}}_{1}}(1 - 2{{а}_{2}}{{{\text{Г}}}_{2}}) + {{{\text{Г}}}_{2}}(1 - 2{{а}_{1}}{{{\text{Г}}}_{1}}) + \\ + \,\,(1 - {\text{Г}})(1 - 2{{а}_{1}}{{{\text{Г}}}_{1}})(1 - 2{{а}_{2}}{{{\text{Г}}}_{2}}) = \\ = 2{{а}_{2}}(1 - 2{{а}_{1}}{{{\text{Г}}}_{1}}){\text{Г}}_{2}^{2} - 2{{а}_{2}}{{{\text{Г}}}_{2}} \times \\ \times \,\,[1 - 2{{а}_{1}}{{{\text{Г}}}_{1}}(1 - {{{\text{Г}}}_{1}})] + 1 - 2{{а}_{1}}{{{\text{Г}}}_{1}}(1 - {{{\text{Г}}}_{1}}), \\ {{G}_{1}} \equiv 4{{а}_{3}}{{{\text{Г}}}_{1}}{{{\text{Г}}}_{2}}[1 - {{а}_{3}}(1 - {\text{Г}})]. \\ \end{gathered} $

Якобиан (9) с учетом представлений (10)–(12), (12а) примет вид

(13)
$J = \frac{{{{e}^{g}}{{J}_{0}}}}{{{{{(1 - {\text{Г}})}}^{4}}}} = \frac{{{{e}^{g}}{{J}_{s}}}}{{{{{(1 - {\text{Г}})}}^{3}}}}.$

Знакоопределяющим множителем для якобиана J в области $\Delta $ согласно равенства (13) является функция Js, а условие (8) можно заменить неравенством

(14)
${{J}_{s}} > 0,\,\,\,\,({{{\text{Г}}}_{1}},{{{\text{Г}}}_{2}}) \in \Delta $
вследствие равенства

(14a)
$J\left| {_{{{{{\text{Г}}}_{1}} = {{{\text{Г}}}_{2}} = 0}}} \right. = {{J}_{s}}\left| {_{{{{{\text{Г}}}_{1}} = {{{\text{Г}}}_{2}} = 0}}} \right. = 1 > 0.$

Найдем вначале условия при которых ${{J}_{s}} > 0$ на всей границе области $\Delta $.

Пусть Г2 = 0 . Тогда границей $\Delta $ является отрезок 0 < Г1 < 1, а необходимым условием для неравенства (14) на этом отрезке является

(15)
$\begin{gathered} {{J}_{s}}\left| {_{{{{{\text{Г}}}_{2}} = 0}}} \right. = {{{\text{Г}}}_{1}} + (1 - {{{\text{Г}}}_{1}})(1 - 2{{а}_{1}}{{{\text{Г}}}_{1}}) = \\ = 1 - 2{{а}_{1}}{{{\text{Г}}}_{1}} + 2{{а}_{1}}{\text{Г}}_{1}^{2} \equiv {{f}_{1}} > 0, \\ 0 < {{{\text{Г}}}_{1}} < 1. \\ \end{gathered} $

Решением неравенства (15) будет (Приложение 1)

(16)
${{а}_{1}} < 2.$

Полагая теперь Г1 = 0, придем к условию

(17)
$\begin{gathered} {{J}_{s}}\left| {_{{{{{\text{Г}}}_{1}} = 0}}} \right. = {{{\text{Г}}}_{2}} + (1 - {{{\text{Г}}}_{2}})(1 - 2{{а}_{2}}{{{\text{Г}}}_{2}}) = \\ = 1 - 2{{а}_{2}}{{{\text{Г}}}_{2}} + 2{{а}_{2}}{\text{Г}}_{2}^{2} > 0,\,\,\,\,0 < {{{\text{Г}}}_{2}} < 1. \\ \end{gathered} $

Аналогично неравенству (15) найдем решение неравенства (17)

(18)
${{а}_{2}} < 2.$

Пусть теперь ${\text{Г}} \equiv {{{\text{Г}}}_{1}} + {{{\text{Г}}}_{2}} = 1$. Тогда, полагая ${{{\text{Г}}}_{2}} = 1 - {{{\text{Г}}}_{1}}$ из (14) и (12а), придем к необходимому условию

(19)
$\begin{gathered} {{J}_{s}}\left| {_{{{\text{Г}} = 1}}} \right. = 1 - 2({{а}_{1}} + {{а}_{2}}){{{\text{Г}}}_{1}}(1 - {{{\text{Г}}}_{1}}) + 4{{а}_{3}}{{{\text{Г}}}_{1}}(1 - {{{\text{Г}}}_{1}}) = \\ = 2({{а}_{1}} + {{а}_{2}} - 2{{а}_{3}}){\text{Г}}_{1}^{2} - 2({{а}_{1}} + {{а}_{2}} - 2{{а}_{3}}){{{\text{Г}}}_{1}} + 1 > 0, \\ 0 < {{{\text{Г}}}_{1}} < 1 - {{{\text{Г}}}_{2}}. \\ \end{gathered} $

Из (19) получим решение этого неравенства относительно величины $({{а}_{1}} + {{а}_{2}} - 2{{а}_{3}})$, аналогичное решениям неравенств (15) и (17)

(20)
${{а}_{1}} + {{а}_{2}} - 2{{а}_{3}} < 2.$

Неравенство (20), как и (16), (18), являются необходимыми для выполнения условия (14). Поэтому нарушение любого из них ведет к неединственности решений уравнений (1), (2), но не при любых значений с1, с2.

Выразим теперь производную функции Js по а2 с учетом равенств (12) и (12а)

(21)
$\begin{gathered} \frac{{\partial {{J}_{s}}}}{{\partial {{а}_{2}}}} = \frac{{\partial ({{G}_{0}} + {{G}_{1}})}}{{\partial {{а}_{2}}}} = \frac{{\partial {{G}_{0}}}}{{\partial {{а}_{2}}}} = \\ = - 2{{{\text{Г}}}_{1}}{{{\text{Г}}}_{2}} - 2{{{\text{Г}}}_{2}}(1 - {\text{Г}})(1 - 2{{а}_{1}}{{{\text{Г}}}_{1}}) = \\ = - 2{{{\text{Г}}}_{2}}[{{{\text{Г}}}_{1}} + (1 - {\text{Г}})(1 - 2{{а}_{1}}{{{\text{Г}}}_{1}})] = - 2{{Г}_{2}} \times \\ \times \,\,[1 - {{{\text{Г}}}_{2}} - 2{{а}_{1}}{{{\text{Г}}}_{1}}(1 - {{{\text{Г}}}_{2}}) + 2{{а}_{1}}{\text{Г}}_{1}^{2}] = - 2{{{\text{Г}}}_{2}}{{F}_{1}}({{{\text{Г}}}_{1}}), \\ \end{gathered} $
(21a)
${{F}_{1}} \equiv 2{{а}_{1}}{\text{Г}}_{1}^{2} - {{v}_{{_{1}}}}{{{\text{Г}}}_{1}} + 1 - {{{\text{Г}}}_{2}},\,\,\,\,{{v}_{1}} \equiv 2{{а}_{1}}(1 - {{{\text{Г}}}_{2}}).$

Аналогично найдем с учетом (12), (12а) производную функции Js по а1

(22)
$\begin{gathered} \frac{{\partial {{J}_{s}}}}{{\partial {{а}_{1}}}} = \frac{{\partial ({{G}_{0}} + {{G}_{1}})}}{{\partial {{а}_{1}}}} = \frac{{\partial {{G}_{0}}}}{{\partial {{а}_{1}}}} = \\ = - 2{{{\text{Г}}}_{1}}[1 - {{{\text{Г}}}_{1}} - 2{{а}_{2}}{{{\text{Г}}}_{2}}(1 - {{{\text{Г}}}_{1}}) + 2{{а}_{2}}{\text{Г}}_{2}^{2}] = \\ = - 2{{{\text{Г}}}_{1}}{{F}_{2}}({{{\text{Г}}}_{2}}), \\ \end{gathered} $
(22a)
${{F}_{2}} \equiv 2{{а}_{2}}{\text{Г}}_{2}^{2} - {{v}_{2}}{{{\text{Г}}}_{2}} + 1 - {{{\text{Г}}}_{1}},\,\,\,{{v}_{2}} \equiv 2{{а}_{2}}(1 - {{{\text{Г}}}_{1}}).$

В Приложении 2 показано, что при условиях (16), (18) выполняются неравенства

(23)
$\begin{gathered} {{F}_{1}}({{{\text{Г}}}_{1}}) > 0,\,\,\,\,0 < {{{\text{Г}}}_{1}} < 1 - {{{\text{Г}}}_{2}}, \\ {{F}_{2}}({{{\text{Г}}}_{2}}) > 0,\,\,\,\,0 < {{{\text{Г}}}_{2}} < 1 - {{{\text{Г}}}_{1}}. \\ \end{gathered} $

Из (23) с учетом (21), (22) вытекают условия монотонной зависимости функции Js по параметрам а1 и а2

(23a)
$\frac{{\partial {{J}_{s}}}}{{\partial {{а}_{1}}}} < 0,$
(23b)
$\frac{{\partial {{J}_{s}}}}{{\partial {{а}_{2}}}} < 0.$

В Приложении 3 выведены с учетом (16), (18), (23а), (23b) условия (24)–(27)

(24)
${{а}_{1}} < 1 + {{а}_{3}},$
(25)
${{а}_{2}} < 1 + {{а}_{3}},$
(26)
${{а}_{1}} + {{а}_{3}} < 4,$
(27)
${{а}_{2}} + {{а}_{3}} < 4,$
при которых выполняется неравенство (14),

При этом можно отметить, что условие (20), необходимое для выполнения неравенства (14) следует из (24) и (25). Действительно, складывая левые части этих неравенств и соответственно их правые части, получим неравенство (20).

Из (14), (13) получим в итоге

(28)
$J > 0,\,\,\,\,({{{\text{Г}}}_{1}},{{{\text{Г}}}_{2}}) \in \Delta .$

То есть выполняется условие единственности (8) решения уравнений (5), (6) относительно ${{{\text{Г}}}_{1}},{{{\text{Г}}}_{2}}$, реализуемое при выполнении неравенств (16), (18), (24)–(27).

Неединственность решений уравнений (1), (2) возникает при нарушении условия (8), т.е. условия (28), но, как отмечено в [6], не при любых значениях с1 > 0, с2 > 0, а в определенных их областях.

ПРИМЕРЫ ЧИСЛЕННОГО РЕШЕНИЯ УРАВНЕНИЙ (5), (6)

Рассмотрим частный случай уравнений (5), (6)

(29)
${{а}_{1}} = {{а}_{2}} = 0,$
(30)
${{\beta }_{1}} = \frac{{{{{\text{Г}}}_{1}}\exp ( - 2{{a}_{3}}{{{\text{Г}}}_{2}})}}{{1 - {{{\text{Г}}}_{1}} - {{{\text{Г}}}_{2}}}},$
(31)
${{\beta }_{2}} = \frac{{{{{\text{Г}}}_{2}}\exp ( - 2{{a}_{3}}{{{\text{Г}}}_{1}})}}{{1 - {{{\text{Г}}}_{1}} - {{{\text{Г}}}_{2}}}},$
(31а)
${{\beta }_{1}} \equiv {{В}_{1}}{{с}_{1}},\,\,\,\,{{\beta }_{2}} \equiv {{В}_{2}}{{с}_{2}}.$

Исследования системы уравнений (5), (6) при ${{а}_{3}} \ne 0$ выпали из поля зрения работ [7, 8], поэтому рассматриваемый теперь случай (29), но при ${{а}_{3}} \ne 0$, представляет интерес. Разделив уравнение (30) на уравнение (31), получим уравнение

(32)
$\frac{{{{{\text{Г}}}_{1}}}}{{{{{\text{Г}}}_{2}}}}\exp (2{{а}_{3}}{{{\text{Г}}}_{1}} - 2{{а}_{3}}{{{\text{Г}}}_{2}}) = {{\beta }_{{12}}},$
(32а)
${{\beta }_{{12}}} \equiv \frac{{{{\beta }_{1}}}}{{{{\beta }_{2}}}}.$

Используя условие (29), найдем решение уравнения (31) относительно Г2

(33)
${{{\text{Г}}}_{2}} = \frac{{{{\beta }_{2}}(1 - {{{\text{Г}}}_{1}})}}{{{{\beta }_{2}} + \exp ( - 2{{а}_{3}}{{{\text{Г}}}_{1}})}} \equiv {{f}_{2}}({{{\text{Г}}}_{1}}).$

Заменяя уравнение (31) его решением (33) относительно Г2, получим два уравнения

(34)
${{\beta }_{1}} = \frac{{{{{\text{Г}}}_{1}}\exp ( - 2{{a}_{3}}{{f}_{2}})}}{{1 - {{{\text{Г}}}_{1}} - {{f}_{2}}}},$
(35)
$\frac{{{{{\text{Г}}}_{1}}}}{{{{f}_{2}}}}\exp (2{{а}_{3}}{{{\text{Г}}}_{1}} - 2{{а}_{3}}{{f}_{2}}) = {{\beta }_{{12}}},$
и каждое из них может служить для определения Г1. Логарифмируя левую и правую части (35), придем к уравнению
(36)
$ - {{а}_{3}} = \frac{{\ln ({{\beta }_{{12}}}y)}}{{2{{{\text{Г}}}_{1}}(y - 1)}},$
где

(36а)
$y \equiv \frac{{{{f}_{2}}}}{{{{{\text{Г}}}_{1}}}}.$

Заменяя уравнение (35) уравнением (36), получим в итоге вместо (34), (35) два уравнения

(37)
$1 = \frac{{{{{\text{Г}}}_{1}}\exp ( - 2{{a}_{3}}{{f}_{2}})}}{{(1 - {{{\text{Г}}}_{1}} - {{f}_{2}}){{\beta }_{1}}}} = \frac{{\exp ( - 2{{а}_{3}}{{f}_{2}})}}{d} \equiv {{{\text{Ф}}}_{0}}({{{\text{Г}}}_{1}}),$
(38)
$ - {{а}_{3}} = \frac{{\ln ({{\beta }_{{12}}}y)}}{{2{{{\text{Г}}}_{1}}(y - 1)}} \equiv {{{\text{Ф}}}_{3}}({{{\text{Г}}}_{1}}),$
где

(38а)
$d \equiv \frac{{{{\beta }_{1}}(1 - {{{\text{Г}}}_{1}} - {{f}_{2}})}}{{{{{\text{Г}}}_{1}}}}.$

Каждое из уравнений (37), (38) с учетом определения y в (36а), d в (38а) и f2 в (33) является уравнением относительно единственного неизвестного Г1 при заданных всех других параметрах, входящих в эти уравнения. Значения Г2 в качестве решения системы (5), (6) и соответствующие найденному Г1 следует определять по выражению функции f2 (33).

Далее ограничимся отрицательными значениями ${{а}_{3}}$

(39)
${{а}_{3}} < 0.$

В случае (39) для функции f21) получим условие ее монотонного убывания по переменной Г1

(40)
$\frac{{\partial {{f}_{2}}}}{{\partial {{{\text{Г}}}_{1}}}} = {{\beta }_{2}}\frac{{2{{а}_{3}}(1 - {{{\text{Г}}}_{1}}){{e}^{{ - 2{{а}_{3}}{{{\text{Г}}}_{1}}}}} - ({{\beta }_{2}} + {{e}^{{ - 2{{а}_{3}}{{{\text{Г}}}_{{\text{1}}}}}}})}}{{{{{({{\beta }_{2}} + {{e}^{{ - 2{{а}_{3}}{{{\text{Г}}}_{1}}}}})}}^{2}}}} < 0,$
так как ${{\beta }_{2}} > 0,\,\,0 \leqslant {\text{Г}} < 1$.

Из (40) следует также монотонная зависимость функции y1)

(41)
$\frac{{\partial y}}{{\partial {{{\text{Г}}}_{1}}}} = \frac{\partial }{{\partial {{{\text{Г}}}_{1}}}}\left( {\frac{{{{f}_{2}}}}{{{{{\text{Г}}}_{1}}}}} \right) = \frac{{\frac{{\partial {{f}_{2}}}}{{\partial {{{\text{Г}}}_{1}}}}{{{\text{Г}}}_{1}} - {{f}_{2}}}}{{{\text{Г}}_{1}^{2}}} < 0,\,\,\,\,0 < {{{\text{Г}}}_{1}} < 1.$

Введем теперь точку ${{{\text{Г}}}_{z}}$ в интервале $0 < {{{\text{Г}}}_{z}} < 1$ из условия

(42)
$y\left| {_{{{{{\text{Г}}}_{1}} = {{{\text{Г}}}_{z}}}}} \right. = 1.$

Подставляя в (42) выражение y из (36а) с учетом (33), получим трансцендентное уравнение для определения ${{{\text{Г}}}_{z}}$

(43)
${{\beta }_{2}} = \frac{{{{{\text{Г}}}_{z}}{{e}^{{ - 2{{а}_{3}}{{{\text{Г}}}_{z}}}}}}}{{1 - 2{{{\text{Г}}}_{z}}}},\,\,\,0 < {{{\text{Г}}}_{z}} < 0.5.$

Функция y1) с учетом (35) и (42) удовлетворяет краевым условиям

(44)
$y\left| {_{{{{{\text{Г}}}_{1}} = 0}}} \right. = + \infty ,\,\,\,\,y\left| {_{{{{{\text{Г}}}_{1}} = {{{\text{Г}}}_{z}}}}} \right. = 1.$

Полагая вначале ${{\beta }_{1}} \ne {{\beta }_{2}}$, будем считать ${{\beta }_{1}} > {{\beta }_{2}}$, т.е.

(45)
${{\beta }_{{12}}} > 1.$

Из (41) и (42) находим

(46)
$y > 1,\,\,\,\,0 < {{{\text{Г}}}_{1}} < {{{\text{Г}}}_{z}}.$

Из (45) и (46) следует

(47)
${{{\text{Ф}}}_{3}}({{{\text{Г}}}_{1}}) > 0,\,\,\,\,0 < {{{\text{Г}}}_{1}} < {{{\text{Г}}}_{z}}.$

При этом из (44) получим асимптотику функции Ф31) на концах интервала $0 < {{{\text{Г}}}_{1}} < {{{\text{Г}}}_{z}}$

(48)
${{{\text{Ф}}}_{3}}\left| {_{{{{{\text{Г}}}_{1}} = + 0}}} \right. = + \infty ,\,\,\,\,{{{\text{Ф}}}_{3}}\left| {_{{{{{\text{Г}}}_{1}} = {{{\text{Г}}}_{z}} - 0}}} \right. = + \infty .$

Вид зависимости функции Ф31) в соседнем интервале ${{{\text{Г}}}_{z}} < {{{\text{Г}}}_{1}} < 1$ определяется неравенством (Приложение 4)

(49)
$\frac{{\partial {{{\text{Ф}}}_{3}}}}{{\partial {{{\text{Г}}}_{1}}}} > 0,\,\,\,\,{{{\text{Г}}}_{z}} < {{{\text{Г}}}_{1}} < 1.$

Введем также единственную с учетом (49) точку Г0 в интервале ${{{\text{Г}}}_{z}} < {{{\text{Г}}}_{0}} < 1$, где

(50)
${{{\text{Ф}}}_{{\text{3}}}}\left| {_{{{{{\text{Г}}}_{1}} = {{{\text{Г}}}_{0}}}}} \right. = {{\left. {\left[ {\frac{{\ln ({{\beta }_{{12}}}y)}}{{2{{{\text{Г}}}_{1}}(y - 1)}}} \right]} \right|}_{{{{{\text{Г}}}_{1}} = {{{\text{Г}}}_{0}}}}} = 0,\,\,\,{{\left. {\,{{\beta }_{{12}}}y} \right|}_{{{{{\text{Г}}}_{{\text{1}}}}{\text{ = }}{{{\text{Г}}}_{{\text{0}}}}}}} = 1.$

Используя определение y в (36а) и выражение функции f2 (33), из (50) найдем трансцендентное уравнение для определения Г0

(51)
${{\beta }_{1}} = \frac{{{{{\text{Г}}}_{0}}({{\beta }_{2}} + {{e}^{{ - 2{{а}_{3}}{{{\text{Г}}}_{0}}}}})}}{{1 - {{{\text{Г}}}_{0}}}}.$

Из (41), (42) и (45) придем к неравенствам

(52)
${{\beta }_{{12}}}y > 1,\,\,\,\,{{{\text{Г}}}_{z}} < {{{\text{Г}}}_{1}} < {{{\text{Г}}}_{0}},\,\,\,\,y\left| {_{{{{{\text{Г}}}_{1}} = {{{\text{Г}}}_{0}}}}} \right. = \frac{1}{{{{\beta }_{{12}}}}} < 1.$

Условия (52) и (41) с учетом определения функции Ф3 в (36) приводят к неравенству

(53)
${{{\text{Ф}}}_{3}}({{{\text{Г}}}_{1}}) < 0,\,\,\,\,0 < {{{\text{Г}}}_{z}} < {{{\text{Г}}}_{1}} < {{{\text{Г}}}_{0}}.$

Асимптотика функции Ф31) при ${{{\text{Г}}}_{1}} \to 1$ вытекает с учетом (33) из равенства

$y\left| {_{{{{{\text{Г}}}_{1}} = 1}}} \right. = 0,$
то есть

(54)
${{{\text{Ф}}}_{{\text{3}}}}\left| {_{{{{{\text{Г}}}_{1}} = 1}}} \right. = + \infty .$

Поэтому область значений Г1, где возможно наличие корней уравнения (38) при а3 < 0, с учетом (47) и (53) сужается к двум интервалам

(54а)
$0 < {{{\text{Г}}}_{1}} < {{{\text{Г}}}_{z}},\,\,\,\,{{{\text{Г}}}_{0}} < {{{\text{Г}}}_{1}} < 1.$

При этом вследствие монотонного роста функции Ф31) в области (49) корень уравнения (38) в интервале ${{{\text{Г}}}_{0}} < {{{\text{Г}}}_{1}} < 1$ является единственным.

Таким образом уравнение (38) при условии (45) и с учетом асимптотик функции Ф31) (48) может иметь в области $0 < {{{\text{Г}}}_{1}} < {{{\text{Г}}}_{z}}$ неединственное решение, а в области ${{{\text{Г}}}_{0}} < {{{\text{Г}}}_{1}} < 1$ единственное.

Пусть теперь

(55)
${{\beta }_{{12}}} = 1.$

В случае (55) точки Гz и Г0 согласно их определению в (43) и в (51) совпадают. При этом сингулярность функции Ф31) в этой точке исчезает

(56)
$\begin{gathered} \mathop {\lim }\limits_{{{{\text{Г}}}_{1}} \to {{{\text{Г}}}_{z}}} {{{\text{Ф}}}_{3}}({{{\text{Г}}}_{1}}) = \mathop {\lim }\limits_{{{{\text{Г}}}_{1}} \to {{{\text{Г}}}_{z}}} \frac{1}{{2{{{\text{Г}}}_{1}}}}\mathop {\lim }\limits_{y \to 1} \frac{{\ln y}}{{y - 1}} = \\ = \frac{1}{{2{{{\text{Г}}}_{z}}}}\mathop {\lim }\limits_{y \to 1} \frac{{{{y}^{{ - 1}}}}}{1} = \frac{1}{{2{{{\text{Г}}}_{z}}}} > 0. \\ \end{gathered} $

Равенство Г0 = Гz и условие (56) означают, что область ${{{\text{Г}}}_{z}} < {{{\text{Г}}}_{{\text{1}}}} < {{{\text{Г}}}_{{\text{0}}}}$, где Ф31) < 0, исчезает, т.е. при условии (55)

(57)
${{{\text{Ф}}}_{3}}({{{\text{Г}}}_{1}}) > 0,\,\,\,\,0 < {{{\text{Г}}}_{1}} < 1.$

При этом функция Ф31) на концах интервала $0 < {{{\text{Г}}}_{1}} < 1$ с учетом ее определения (38) имеет асимптотики

(58)
${{{\text{Ф}}}_{3}}\left| {_{{{{{\text{Г}}}_{1}} = 0 + 0}}} \right. = + \infty ,\,\,\,\,{{{\text{Ф}}}_{3}}\left| {_{{{{{\text{Г}}}_{1}} = 1 - 0}}} \right. = + \infty .$

Возможное отсутствие решений уравнения (38), например, при –а3 < 1 (а3 < 0), будет означать лишь, что решение системы (30), (31) следует искать среди корней уравнения (37), т.к. краевые значения функции Ф01) в интервале $0 < {{{\text{Г}}}_{1}} < 1$

(59)
${{{\text{Ф}}}_{0}}\left| {_{{{{{\text{Г}}}_{1}} = 0}}} \right. = 0,\,\,\,\,{{{\text{Ф}}}_{0}}\left| {_{{{{{\text{Г}}}_{1}} = 1 - 0}}} \right. = + \infty $
обеспечивают существование корня в этом интервале.

Учитывая краевые значения (58) функции Ф3, парное (неединственное) решение уравнения (38) находится при заданном значении а3 < 0 (как пересечение графика функции Ф31) с горизонтальной прямой –а3) в интервале $0 < {{{\text{Г}}}_{1}} < {{{\text{Г}}}_{z}}$, если ${{\beta }_{{12}}} > 1$. График этой функции легко рассчитывается. Если же ${{\beta }_{{12}}} = 1$, то решение системы (30), (31) следует искать среди корней обоих уравнений (38) и (37) в зависимости от величины (–а3), т.е. от критерия единственности в данных условиях (24).

В табл. 1–5 приведены рассчитанные указанным способом значения Г1 и Г2 и соответствующие значения функций Ф3 и Ф0 при фиксированных ${{\beta }_{1}},{{\beta }_{2}},{{а}_{3}}$, и значения ${{{\text{Г}}}_{z}}$ и ${{{\text{Г}}}_{0}}$, рассчитанные из уравнений (43) и (51), когда ${{{\text{Г}}}_{z}} \ne {{{\text{Г}}}_{0}}$.

Таблица 1.  

а3 = –0.9, ${{\beta }_{1}} = 4,\,\,{{\beta }_{2}} = 2$

Г1 Г2 Ф3 Ф0 Гz Г0
0.665 0.12617 0.8991 0.9991 0.3418 0.478
Таблица 2.  

а3 = –0.9, ${{\beta }_{1}} = 4,\,\,{{\beta }_{2}} = 4$

Г1 Г2 Ф0
0.398 0.398 0.9998
Таблица 3.  

а3 = –1.5, ${{\beta }_{1}} = 4,\,\,{{\beta }_{2}} = 2$

Г1 Г2 Ф3 Ф0 Гz Г0
0.741 0.0461 1.499 0.9992 0.3071 0.364
Таблица 4.  

а3 = –1.5, ${{\beta }_{1}} = {{\beta }_{2}} = 4$

Г1 Г2 Ф3 Ф0
0.1622 0.5956 1.5006 0.9994
0.596 0.1619 1.5009 1.0007
Таблица 5.  

а3 = –3, ${{\beta }_{1}} = 4,\,\,{{\beta }_{2}} = 2,$ Гz = 0.2419, Г0 = 0.316

Г1 Г2 Ф3 Ф0
0.0946 0.481 3.0008 0.9993
0.054 0.559 2.9989 1.0001
0.794 0.00345 3.0003 1.0005

Полученные численно результаты полностью соответствуют условиям единственности (16), (18), (24)–(27) (–а3 < 1 при условии (29)) в табл. 1, 2 и демонстрируют возможность неединственности (при –а3 > 1), но при определенных значениях ${{\beta }_{1}},\,{{\beta }_{2}}$табл. 3–5.

В случае однокомпонентной адсорбции по изотерме Фрумкина (4) при а1 > 2 (притяжение между молекулами адсорбата) как известно, образуется S-образность ее формы. При а1 < 0 (отталкивание между молекулами) S-образность исчезает. Для двухкомпонентной изотермы (5), (6) Фрумкина–Дамаскина полученные результаты уже в частном случае (29) (табл. 4, 5) демонстрируют обратный случаю однокомпонентной адсорбции эффект – неединственность решения системы (30), (31) при а3 < 0 (отталкивание между молекулами “а” и “b”) с ростом величины (–а3).

ЗАКЛЮЧЕНИЕ

1. Исследовалась система нелинейных уравнений двухкомпонентной адсорбции в задаче об условиях расслоения в таком адсорбционном слое при адсорбции по изотерме Фрумкина–Дамаскина.

2. Выведены аналитически условия единственности решения нелинейных уравнений двухкомпонентной адсорбции относительно поверхностных концентраций каждого компонента при заданных их объемных концентрациях в случае трехпараметрического описания совместной адсорбции по изотерме Фрумкина–Дамаскина.

3. На частном примере параметров двухкомпонентной адсорбции по изотерме Фрумкина–Дамаскина показано, что большое значение параметра взаимодействия между молекулами разных адсорбатов (но при отталкивании – отрицателное) приводит к неединственности решения нелинейных уравнений совместной адсорбции, в отличие от случая однокомпонентной адсорбции при отрицательных значениях аттракционного параметра Фрумкина.

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

  1. Подгаецкий Э.М. // Физикохимия поверхности и защита материалов. 2012. Т. 48. № 4. С. 155.

  2. Фрумкин А.Н. // Тр. хим. ин-та им. Л.Я. Карпова. 1925. № 4. С. 56.

  3. Оуара К., Лифшиц В.Г., Саранин А.А. и др. // Введение в физику поверхности. М.: Наука, 2006. 490 с.

  4. Дамаскин Б.Б. // Электрохимия. 1977. Т. 12. № 6. С. 816.

  5. Подгаецкий Э.М. // Электрохимия. 1974. Т. 10. № 4. С. 666.

  6. Подгаецкий Э.М. // Электрохимия. 1975. Т. 11. № 11. С. 1759.

  7. Дамаскин Б.Б. // Электрохимия. 1969. Т. 5. № 2. С. 346.

  8. Харкац Ю.И. // Электрохимия. 1980. Т. 16. № 12. С. 1820.

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