Математическая модель химических осцилляций, возникающих в гомогенной системе цистеин–оксигенированные комплексы железа (II)

УДК 541.128.7

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ХИМИЧЕСКИХ ОСЦИЛЛЯЦИЙ, ВОЗНИКАЮЩИХ В ГОМОГЕННОЙ СИСТЕМЕ ЦИСТЕИН-ОКСИГЕНИРОВАННЫЕ КОМПЛЕКСЫ ЖЕЛЕЗА (II) У.Г. Магомедбеков, У.Г. Гасангаджиева, Х.М. Гасанова

Проведен анализ математической модели кинетики гомогенного окисления цистеина в присутствии оксигенированных комплексов железа (II) с диметилглиоксимом и цитози-ном, представленной в виде системы из трех нелинейных дифференциальных уравнений. Проведено упрощение модели при помощи стехиометрического анализа предложенной кинетической схемы. На основе качественного анализа системы дифференциальных уравнений установлено, что реализуется одно стационарное состояние с особой точкой типа «пространственное седло-фокус», из которого возможна бифуркация типа Андронова-Хопфа. При численном решении системы дифференциальных уравнений для разных начальных условий получено, что математическая модель имеет решение в виде предельного цикла, и процесс протекает в колебательном режиме при исходных концентрациях цистеина и катализатора, сопоставимых с таким их содержанием, при котором экспериментально наблюдаются химические осцилляции.

Ключевые слова: цистеин, окисление, колебательные процессы, математическое моделирование, кинетический анализ, стехиометрический анализ.

При интерпретации экспериментальных результатов по колебательным химическим реакциям наряду с привлечением методов нелинейной динамики и термодинамики открытых систем существенную роль играет математическое моделирование - установление качественных и количественных изменений в поведении нелинейных динамических систем при изменении описывающих их параметров [1]. В рамках данного подхода, как правило, проводят качественный анализ и численное решение кинетической модели, представленной в виде системы автономных нелинейных дифференциальных уравнений [2].

В настоящей статье приведены результаты по математическому моделированию кинетических закономерностей химических осцилляций, проявляющихся при гомогенном окислении цистеина (Cys, RSH2) в присутствии оксигенированных ком -плексов железа (II) с диметилглиоксимом (ДМГ) и цитозином (Cyt) (Cat). Основное внимание, наряду с анализом адекватности модели, обращено на установление возможности проявления кинетической неустойчивости при использовании в качестве управляющих параметров концентрации реагента и катализатора.

Экспериментальная часть

Исследования проводили на основе методического подхода, апробированного ранее при изучении колебательных окислительно-восстановительных превращений разных биосубстратов; особенности проведения эксперимента описаны нами в [2]. Особое внимание обращалось на чистоту реагентов. Все реагенты, использованные при выполнении работы (цистеин (Cys, RSH2), диметилглиоксим (ДМГ), цитозин (Cyt), трис-(гидроксиметил) аминометан, сульфат железа(11) (FeSO47H2O)) соответствовали квалификации «х.ч.».

Катализатор (cat) в виде комплексного соединеня железа (II) с ДМГ и Cyt получали смешением растворов соли железа и соответствующих реагентов при молярном соотношении Fе(II):ДМГ:Cyt = 1:2:2. В процесе приготовления катализатора соблюдали следующий порядок смешивания реагентов:

раствор железа (II) - ДМГ - Cyt.

Для каждого опыта использовали свежеприготовленный раствор катализатора. Для приготовления реакционной смеси вносили раствор реагента (Cys) и катализатора определенной концентрации в стеклянный реактор, а затем доводили объем до 50 мл буфер-

ной смесью. В начале каждого эксперимента для насыщения раствора кислородом через реакционную

з

среду пропускали ~20 см воздуха.

Для наблюдения за динамикой поведения системы проводили регистрацию изменения потенциала (Де, мВ, точность ±0,01 мВ) точечного платинового электрода (£ = 1 мм ) относительно хлорсеребряного во времени (г, с) без перемешивания раствора.

Выбор метода исследования обусловлен тем, что значения потенциала определяются соотношением окисленных и восстановленных форм компонентов реакционной смеси, следовательно, по изменению концентрации реагентов во времени можно отслеживать функцию состояния системы [1, 3].

Одна из кривых зависимости изменения потенциала во времени для условий, когда ссу8 = 6,25 10-3 моль/л, сса = 2,5 10 5 моль/л, г = 40СС и рН 7,65, представлена на рис. 1, а, эта же кривая после снятия тренда показана на рис. 1, б.

Результаты эксперимента показали, что гомогенное окисление цистеина в присутствии в качестве катализатора оксигенированных комплексов железа (II) с диметилглиоксимом и цитозином протекает в колебательном режиме, причем флуктуацион-ные явления реализуются в условиях, когда ССу8 = (2,5—7,5)10-3 моль/л, Сса1 = (0,2-1,0)104 моль/л, г = 40-60°С, рН 7,60-7,75. Появлению колебаний предшествует индукционный период, наличие которого

Рис. 1. Изменение потенциала (ДЕ) во времени (г): а - экспериментальная кривая, б - кривая после снятия тренда

указывает на автокаталитическии характер протекающих процессов.

Обсуждение результатов Обоснование кинетической схемы процесса

Анализ литературных данных [4-6] показывает, что механизм процесса окисления молекулярным кислородом серосодержащих биосубстратов в присутствии комплексных соединении переходных металлов сложен и до конца не установлен. Поэтому при составлении кинетической схемы мы использовали данные, полученные другими авторами [4-8].

Известно [4-8], что цистеин содержит три атома водорода, способных к ионизации (по сульфгидрильной, карбоксильной и аминогруппе с рК 1,71; 8,37 и 10,70 соответственно). В условиях опыта (рН 7,60-7,75) он находится в ионизированной форме по карбоксилу в виде Я8Ы,- и ионизированный реагент подвергается окислению молекулярным кислородом, связанным со смешаннолигандным комплексом железа (II) с диметилглиоксимом и азотистым основанием с образованием пероксида водорода и ион-радикала типа ЯБ * . Кроме того, согласно [7], соединения железа (II) способны активировать молекулярный кислород, восстанавливая его до О^- , Н202, Н02 или Н20. Взаимодействие сульфгидрильных группировок цистеина с перекисью водорода в водных растворах в интервале рН 6-12 [4-8] протекает по схеме:

Я8Ы- + Ы202 ^ Я80Ы- + Ы20.

На основании этих предпосылок в качестве модельного процесса гомогенного окисления цистеина выбрана схема, включающая следующие стадии:

В стадию 6 данной кинетической схемы введен коэффициент / (/ > 1), который учитывает возмож-

ность автокаталитического производства частиц НО2 при / > 1. В этой связи необходимо сделать одно пояснение. Коэффициент / по стехиометрии должен быть равен единице. Однако каждая из стадий, приведенных в рассматриваемой схеме, и стадия 6 в том числе, является сложным процессом, который может протекать по разным направлениям и с образованием разных продуктов. В определенных условиях взаимодействие между частицами Я80Ы- и Н2О2 с образованием НО2 может протекать в автокаталитическом режиме. Поэтому гипотеза о том, что/ может быть больше 1 не лишена смысла.

При обосновании предлагаемого механизма учитываются два обстоятельства:

1) с целью упрощения модели нами отобраны те стадии из общего механизма реакции окисления исследуемого реагента, которые учитывают особенности процесса после выхода его в колебательный режим;

2) наличие автокаталитических стадий (реакции 3, 4 по отношению к частицам 118 ' и 3, 5 (при/> 1) по отношению к радикалам НО2), а также реализация обратных связей между отдельными стадиями процесса взаимодействия продуктов реакции с исходными веществами и интермедиатами (реакции 1, 2 и 2, 3); эти условия являются необходимыми для возникновения критических явлений.

Предложенная кинетическая схема включает девять реактантов (компонентов реакционной среды), следовательно, для объяснения особенностей кинетики протекающих процессов необходимо проанализировать систему из девяти автономных обыкновенных дифференциальных уравнений (ОДУ) нелинейного типа. Проведение качественного анализа такого типа систем, их численное интегрирование, а также интерпретация получаемых результатов являются довольно трудоемкими процессами и в большинстве случаев невозможны. Поэтому проведено упрощение модели с целью уменьшения числа дифференциальных уравнений в их системах на основе стехиометрического анализа предложенной кинетической схемы.

Стехиометрический анализ кинетической схемы

Стехиометрический анализ рассматриваемой кинетической схемы проведен с использованием подходов, разработанных в [9-10] с целью выявить линейные связи между реактантами и стадиями процесса. Согласно этим подходам, при нумерации компонентов

реакционной среды и стадий механизма по порядку

— 2+

и при введении обозначений (Ах - Я8Н ; А2-¥ек02 ; А3 - ; А4 - ЯО'2; А5 - А6 - Н202; А1 - Я801Г а8 - Н20; А9 - 02) для / = 2 система однородных стехиометрических уравнений, состоящая из восьми реакций и девяти компонентов реакционной среды, приводится к виду, представленному на схеме 1. Полученной системе уравнений соответствует стехио-метрическая матрица В (схема 2), где строки матрицы соответствуют реактантам, а столбцы - стадиям. Неравенство ранга матрицы числу стадий (матрица В имеет неполный столбцовый ранг) свидетельствует о том, что данная матрица содержит хотя бы один на-

бор линейно независимых вектор-столбцов. В нашем случае таковыми являются вектор-столбцы Ь (/ = 1, 2, 3, 5, 7), система которых соответствует базисным стадиям приведенной кинетической схемы, и которые, в свою очередь, составляют базис стехиометрического пространства. Из полученного базиса выстраивается базисная матрица В (схема 3), обладающая свойствами гкВ = гкВ, причем ёеК^ В)ф 0 [15, 16]. Данный результат показывает, что ранг матрицы В равен пяти (гкВ = 5), откуда следует, что количество ключевых реактантов, которое необходимо выбрать для описания кинетики рассматриваемого процесса, соответствует также пяти. Для подтверждения данного за-

С х е м а 1

-А! -А2 А3 а4 А5 0 0 0 0 = 0

-А! 0 А3 -а4 0 А6 0 0 0 = 0

0 0 -2а3 -2а4 0 0 2а7 0 а9 = 0

-А! 0 2а3 0 0 0 -а7 А8 0 = 0

0 0 А3 2а4 0 "А6 -а7 0 0 = 0

0 0 0 -2а4 0 А6 0 0 а9 = 0

0 А2 0 0 -А5 0 0 0 -А9 = 0

С х е м а 2

В =

-1 -1 0 -1 0 0 0

-1 0 0 0 0 0 1

1 1 -2 2 1 0 0

1 -1 -2 0 2 -2 0

1 0 0 0 0 0 -1

0 1 0 0 -1 1 0

0 0 2 -1 -1 0 0

0 0 0 1 0 0 0

0 0 1 0 0 1 -1

С х е м а 3

ключения необходимо проанализировать схему рассматриваемой реакции на правильность и полноту с использованием матрицы стехиометрических коэффициентов (элементной матрицы) Н, которая устанавливает связь между числами молей элементов и химических соединений.

Для описания состава реактантов в качестве первичных выбраны следующие элементы: G1 - Я8, G2 -Ы, G3 - Бек, G4 - О (Ог - столбик соответствующий первичному элементу). Тогда, элементная матрица Н, строки которой соответствуют компонентам, имеет вид, представленный на схеме 4.

Следует указать, что на пересечении столбца и строки матрицы Н можно получить состав вещества в целом по количеству первичных элементов.

Анализ матрицы Н с использованием ее сингулярного разложения и обращения квадратной матри-т т

цы Н Н (Н - транспонированная элементная матрица) показывает, что ранг ее соответствует четырем (ткН = 4).

Этот результат свидетельствует о том, что стехи-ометрическая и элементная матрицы представляют собой матрицы полного столбцового ранга. При этом сумма рангов этих матриц равна количеству компонентов реакционной среды, т.е. гкВ + ткН = 9. Отсюда следует, что система стехиометрических уравнений (схема реакции) является полной.

При проведении стехиометрического анализа в качестве ключевых нами выбраны следующие реак-

танты: Рер22+ (Л2% ЯБ ' (А3), Я0'2 (Л4% Б^ОН" (А7) и Ы20 (А8); тогда матрицы неключевых (Внк) и ключевых (Вк) компонентов соответственно будут иметь вид, представленный на схеме 5. Ранг матрицы соответствует ткВк= гкВ к = 5, следовательно, матрица Вк обращаема.

Известно [8], что размерность подсистемы дифференциальных уравнений в полной системе кинетических уравнений можно снизить до ранга стехиоме-трической матрицы ключевых реактантов ткВк (в рассматриваемом случае до пяти). Важно, что при этом дифференциальные уравнения записываются относительно концентраций только ключевых компонентов реакционной среды. В рассматриваемом случае подсистема дифференциальных уравнений еще не является замкнутой. Для преодоления данной некорректности необходимо подсистему дифференциальных уравнений относительно ключевых реактантов дополнить системой линейных связей между концентрациями ключевых и неключевых реактантов, после чего полученная система алгебро-дифференциаль-ных уравнений становится замкнутой.

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

С = С

— нк нк

+ анк(ск-сн°к),

С х е м а 4

С х е м а 5

где Онк - матрица линейного преобразования вектора изменения состава ключевых реактантов Ск в вектор изменения состава неключевых реактантов Сн; С°нк - вектор мольных концентраций реактантов в начальный момент времени (вектор начального состава). Матрица линейного преобразования Онк представлена на схеме 6. Тогда выражение зависимости вектора концентраций неключевых реактантов от концентрации ключевых записывается соотношением, представленным на схеме 7.

На основании проведенного анализа полную систему кинетических уравнений с учетом линейных связей между ключевыми и неключевыми компонентами реакционной среды можно представить в виде системы алгебро-дифференциальных уравнений (схема 8), где Сх - [ВНЯЛ; С2- |Те^022+]; С3- [Яв" ]; С4- [НО*]; С5-\¥ек2+]; С6 -[ Н202]; С7 -[ Я801Г]; С8 - [Н20]; С9 - [02]. Таким образом, проведенный стехиометрический анализ обсуждаемой кинетической схемы позволяет снизить размерность подсистемы дифференциальных уравнений в полной системе кинетических уравнений с 9 до 5.

Следует иметь в виду, что к некоторым ключевым реактантам можно применять принципы стационарной кинетики, а именно принцип квазистационарности Боденштейна-Темкина [11]. Так, концентрации разных форм катализатора значительно меньше концентрации устойчивых реактантов, следовательно, можно принять, что ^С2/Л = 0. По этой же причине можно пренебречь изменением концентрации воды, так как реакция проводится в водной среде и поэтому можно принять = 0. Таким образом,

применение принципов стационарной кинетики к рассматриваемой системе позволяет дифференциальные уравнения относительно С2 и С8 преобразовать к виду алгебраических, число дифференциальных уравнений в полной системе алгебро-дифференциальных уравнений уменьшается до трех, и в этом виде (дополненная четырьмя линейными связями относительно неключевых реактантов и двумя алгебраическими уравнениями, полученными при использовании принципов стационарной кинетики) она представлена на схеме 9. Введем обозначения: [Я8Ы-]исх = СА; С2 = Сь; с3 = Сх; с4 = с; с7 = С; с1 = СА- Сх- С2; тогда дифферен-

=

— НК ВК

-1

о -

1 о

0,5 -0,5

0 0

-0,5 -0,5

С х е м а 6

0 -2

0 0

0,5 -1

-0,5 0

С х е м а 7

или

С с = с

= с\ ~(с3 -с3°)-2(с8 -с8и) (с2

= с,

б

+ 0,5(с7

С

о>

2>

+ 0,5(с3-^3 о

С,0)-

0,5 (С4-С°4)-

-с7и)-(с8

С8°)

С9 =Сд -

-С,0)-

С х е м а 8

С х е м а 9

С х е м а 10

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

Качественный анализ математической модели

Основными задачами анализа динамики неравновесных процессов являются исследование существования стационарных состояний (единственность или множественность), определение устойчивости

Рис. 2. Зависимость изменения концентрации частицХ1 = ; у = ; 2 1 = от времени

1

/ 1

/ / / * \ 2105

3 10-5/ 1-ю-4/

Рис. 3. Зависимость изменения концентрации частицX. = ; У. = ; 2г; = от времени

этих состояний с помощью какого-либо подхода (например, по Ляпунову), исследование проведения реакции в различных пространствах - фазовом, динамическом и (или) параметрическом. Следует подчеркнуть, что проявление критических явлений в химической кинетике возможно только в том случае, когда стационарное состояние неустойчиво. Для получения ответа на эти вопросы необходимо использование методов качественного анализа дифференциальных уравнений и решение их численными методами [12-15].

Введем обозначения: х = юСх, у = цС^ 2 = уС% г = дг', а = СА, Ь = Сь, ю = к3/к2, п = к2/к1, у = к4/к1, д = к4, в = к2к4кк1, в к^/к^ [л к^/^, р ф

к5/к4, а = к6к1/к22 и запишем систему уравнений, приведенных на схеме 10, в безразмерном виде:

в(dx/dt) = (а - цх - р2)Ь + (а - цх - р2)у - ху +

+ 2(а - цх - р2)2 +

в,(dy/dt) = (а - цх - р2)Ь + (а - цх - р2)у - ху +

+ & - 2ау2;

dz/dt = ху -

Стационарные состояния определяются на основе решения системы алгебраических уравнений (2), полученных приравниванием нулю левых частей системы дифференциальных уравнений (1):

(а - цх - р2)Ь - ху + 2(а - цх - р2)г +

+ (а - цх - р2)у + ф2 = 0;

(а - цх - р2)Ь - ху - (а - цх - р2)у + - 2ау = 0;

ху -

Можно показать, что неподвижная (стационарная) точка является единственной с координатами

х = а/ц, у = 0, 2 = 0.

Матрица Раусса-Гурвица, полученная в результате линеаризации нелинейной системы дифференциальных уравнений, при помощи которой можно проанализировать устойчивость стационарного состояния, представляется для нашего случая в виде, показанном на схеме 11.

С х е м а 11

X3 + А1Х2 + А2Х + А3 = 0

с учетом (3) соответствуют:

А1 = [ф + Ьц + а/ц ]; А2 = = [Ьцф + а/ц(Ьр+(1-У)ф)]; А3 = -афЬ(/- 1).

Анализ устойчивости стационарного состояния проведен для двух случаев: 1) / = 1 (отсутствие автокаталитического производства частиц НО2 (реакции 4 и 6 в кинетической схеме)); 2) / > 1 (наличие автокатализа по отношению к указанным частицам).

В первом случае коэффициент А3 = 0, следовательно, корни характеристического многочлена (4) соответствуют: = 0, Х2 Ф 0 и Х3 Ф 0. В соответствии с [13-14], если одно из значений корней характеристического уравнения обращается в нуль, то в такой системе в окрестности стационарного состояния возможна смена качественных структур при переходе через границу области Раусса-Гурвица, и следовательно, может реализоваться неустойчивое состояние.

При / > 1 коэффициент А3 и все корни характеристического уравнения (Х1, Х2 и Х3) отличны от нуля, значит действительная часть хотя бы одного из кор -ней будет положительной. В такой ситуации, в соответствии с теоремой Ляпунова [12], особая точка является неустойчивой, и возможно возникновение химических нестабильностей.

Для установления типа стационарной точки, согласно [13], проведем анализ выражения

А*=-А12А2+4А13Л^

-4А3 -18А1А2А3+27А3 .

Отсюда получаем, что коэффициенты (А1, А2 и А3) характеристического уравнения

Анализ этого выражения, проведенный с помощью процедуры, рассмотренной в [13], показывает, что при А* > 0 (или | - 18 А, а2а3 + 27 А 32| > |- а2

2 33 3 I23 3 I

а2 +4 а1 А 3+ 4 А3 |)) рассматриваемое стационарное состояние становится неустойчивым, причем особая точка в этом случае относится к типу «пространственное седло-фокус». При указанных условиях выполнимы условия реализации бифуркации типа Андронова-Хопфа (бифуркация в предельный цикл (тор)), т.е. в исследуемой системе возможно возникновение критических явлений в виде проявления колебательного режима.

Численная реализация математической модели

Изменения, происходящие в динамических системах любого типа, в частности химических, и их эволюция проявляются главным образом при изменении определенных управляющих параметров, входящих

в систему дифференциальных уравнений, представленной в виде математической модели исследуемого процесса. В рассматриваемом случае в качестве управляющих параметров выбраны концентрации реагента (цистеин) и катализатора (оксигенированные комплексы железа (II)). Численное интегрирование проведено с помощью программы Mathcad с подпрограммой Як^ар1:, которая иллюстрирует определенное согласие с физическими явлениями в системах с сосредоточенными параметрами и позволяет представлять полученные результаты в терминах траекторий динамической системы и ее фазовых портретов.

В результате расчета констант скоростей на основании данных работ [4-8] получены следующие величины: к1 ~ 4,0102; к2 ~ 4,5106; к3 ~ 4,2104; к4 ~ 1,2106; к5 ~ 2,1106; к6 ~ 9,0105 (лмольЛ1). Тогда параметры е, е', ц, р, ф, о и 0 с учетом величин констант скоростей принимают значения: е ~ 3,2105; е' ~ 2,710-1; ц ~ 1,1102; р ~ 3,310-4; ф ~ 1,75; о ~ 1,8105.

На рис. 1-3 представлены результаты численного интегрирования математической модели в координатах «концентрация реагентов-время», а также фазовые портреты системы, полученные для условий: а = 1,010-4; Ь = 2,510-5; е ~ 3,2105; е' ~ 2,710-1; ц ~ 1,1102; р ~ 3,310-4; ф ~ 1,75; о ~ 1,8 105 при / = 2 (число точек равно 2000, а величина шага соответствует 1).

Полученные данные показывают (рис. 3, а), что при указанных концентрациях реагента и катализатора реализуется колебательный режим. Вместе с этим, результаты по расчетам фазовых портретов (рис. 3, б, в) свидетельствуют о реализации особой точки типа «пространственное седло-фокус», из которой возможна бифуркация в тор (бифуркация типа Ан-дронова-Хопфа); эти данные находятся в согласии с

результатами качественного анализа математической модели.

Для выяснения характера эволюции динамики поведения системы был проведен компьютерный эксперимент для различных значений концентраций реагента (а) и катализатора (Ь), выбранных в качестве управляющих параметров. В качестве примера на рис. 4 представлены результаты расчета при значениях а = 4,010-4 и Ь = 2,510-5.

В результате численного анализа математической модели, получено, что область существования химических колебаний в исследуемой системе соответствуют следующим концентрациям реагента и катализатора: Ссу = 1,010 4 -5 10-4, Ссай = 1,0 10-52,5 10-5 (моль/л). Полученные данные демонстрируют удовлетворительное согласие с экспериментально наблюдаемыми значениями концентраций реагента и катализатора, при которых реализуется флуктуацион-ный режим.

Наряду с этим, анализ рассматриваемой математической модели позволяет оценить важную характеристику колебательных химических реакций - индукционный период. Индукционный период, полученный на основе расчетов, соответствует 850-900 с, что сопоставимо с экспериментальными величинами (~900-1200 с, см. рис. 1, б).

Таким образом, на основе полученных результатов можно заключить, что предложенная математическая модель удовлетворительно описывает процесс гомогенного окисления цистеина в присутствии оксигени-рованных комплексов железа (II) с диметилглиокси-мом и цитозином в колебательном режиме.

Авторы выражают благодарность докт. хим. наук, профессору М.З. Зейналову за консультации в части стехиометрического анализа кинетической схемы.

Работа выполнена на базе НОЦ «Нелинейная химия», ЦКП «Аналитическая спектроскопия», ГК №16.552.11.7051 по ФЦП «Исследования и разработки по приоритетным направлениям развития научно-технического комплекса России на 2007 - 2012 годы» и при финансовой поддержке РФФИ

СПИСОК ЛИТЕРАТУРЫ

1.Магомедбеков У.Г. // Вестн. Моск. ун-та. Сер. 2. Химия. 2001. 42. С. 75.

2.Магомедбеков У.Г., Гасангаджиева У.Г., Гасанова Х.М., Магомедбеков Н.Х. // Рос. хим. ж. (Ж. Рос. хим. об-ва им. Д.И. Менделеева), 2009. 56. № 6. С. 74.

ЗМагомедбеков У.Г., Гасангаджиева У.Г., Гасанова Х.М., Магомедбеков Н.Х. Гаджибалаева З.М., Магомедова Р.А. // Изв. вузов. Сев.-Кавк. регион. Естественные науки. 2008. № 6 (146). С. 48.

4.Общая органическая химия. Т. 5. Соединения фосфора и серы / Под ред. И.О. Сазерленда и Д.Н. Джонса. М., 1983.

5. Торчинский Ю.М. Сульфгидрильные и дисульфидные группы белков. М., 1971.

6.Harman L.S., Mottley C., Mason R.P. // J. Biol. Chem. 1984. 259. № 9. Р. 5606.

7.СычевА.Я., ИсакВ.Г. // Усп. Хим. 1995. 64. № 12. С. 1183.

8.Radil R., Beckman J.S., Bush K.M., Freeman B.A. // J. Biol. Chem. 1991. 266. № 7. Р. 4244.

9. Горский В.Г., Зейналов М.З. Физико-химические и математические основы феноменологической кинетики сложных реакций. Махачкала, 1997.

10. Горский В.Г., ЗейналовМ.З., Гаджибалаева З.М. // ТОХТ. 2006. 40. № 6. C. 659.

11. Темкин О.Н. Гомогенный металлокомплексный катализ. Кинетические аспекты. М., 2008.

12. Ляпунов А.М. Общая задача об устойчивости движения. Черновцы, 2000.

13.Баутин Н.Н. Поведение динамических систем вблизи границ области устойчивости. М., 1984.

14. Малинецкий Г.Г. Хаос. Структуры. Вычислительный эксперимент: Введение в нелинейную динамику. М., 2000.

15. Быков В.И., Цыбенова С.Б. Нелинейные модели химической кинетики. М., 2011.

Поступила в редакцию 20.04.13

MATHEMATICAL MODEL OF CHEMICAL OSCILLATIONS ARISING IN THE HOMOGENEOUS SYSTEM OF CYSTEINE - OXIGENATED COMPLEXES OF IRON

U.G. Magomedbekov, U.G. Gasangadzhieva, Kh.M. Gasanova.

The analysis of the mathematical model of the kinetics of homogeneous oxidation of cysteine in the presence of oxygenated complexes of iron (II) with dimethylglyoxime and cytosine represented as a system of three nonlinear differential equations was carried out. Simplification of the model on the basis of stoichiometric analysis of the offered kinetic scheme was conducted. On the basis of qualitative analysis of the system of differential equations it was determined that one stationary state with a singular point of the "spatial saddle-focus" type from which the bifurcation of Andronov-Hopf type is realized. Numerically solving the system of differential equations for different initial conditions it was found that the mathematical model has a solution in the form of a limit cycle and the process runs in an oscillatory mode at initial concentrations of cysteine and catalyst comparable with its amounts at which chemical oscillations can be experimentally observed. It was made a conclusion that the mathematical model satisfactorily describes processes running in the observed system.

Key words: oscillatory processes, mathematical modeling, kinetic analysis, stoichiometric analysis, cysteine, oxidation.

Сведения об авторах: Магомедбеков Ухумаали Гаджиевич - зав. кафедрой общей и неорганической химии Дагестанского государственного университета, профессор, докт. хим. наук (ukhgmag@mail.ru); Гасангаджиева Умукусум Гусейновна - доцент кафедры общей и неорганической химии Дагестанского государственного университета, канд. хим. наук; ГасановаХадижат Магомедовна - доцент кафедры общей и неорганической химии Дагестанского государственного университета, канд. хим. наук.