Документ взят из кэша поисковой машины. Адрес оригинального документа : http://chem.msu.ru/rus/vmgu/116/443.pdf Дата изменения: Tue Oct 4 14:49:16 2011 Дата индексирования: Sun Apr 10 11:06:34 2016 Кодировка: Windows-1251 Поисковые слова: великий аттрактор

ВЕСТН. МОСК. УН-ТА. СЕР. 2. ХИМИЯ. 2011. Т. 52. ? 6

443

УДК 541.128.7

ОБРАЗОВАНИЕ ДИССИПАТИВНЫХ СТРУКТУР ПРИ ГОМОГЕННОМ ОКИСЛЕНИИ БИОСУБСТРАТОВ: ОЦЕНКА ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК И ПАРАМЕТРИЗАЦИЯ ВРЕМЕННЫХ РЯДОВ
У.Г. Магомедбеков, У.Г. Гасангаджиева, Х.М. Гасанова, Н.Х. Магомедбеков (Дагестанский государственный университет, кафедра общей и неорганической химии; e-mail: ukhgmag@mail.ru)
В реакции жидкофазного окисления цистеина показана возможность образования диссипативных структур в системе биосубстрат-оксигенированные комплексы железа (II). Приведены результаты по оценке динамических характеристик исследуемых процессов (тип динамики, размерности фазового пространства и аттрактора, показатели Ляпунова и энтропии Колмогорова-Синая). На основе подходов фликкер-шумовой спектроскопии проведена параметризация полученных временных рядов.

Ключевые слова: диссипативные структуры, фазовое пространство, аттрактор, показатели Ляпунова, энтропия Колмогорова-Синая, фликкер-шумовая спектроскопия, биосубстрат, цистеин, окисление.

В исследовании колебательных химических реакций важная роль принадлежит установлению типа динамики протекающих процессов с определением его основных параметров [1]. Принципиальным является получение ответа на вопросы: какого типа эволюцией порожден тот или иной аналитический сигнал и каким образом можно идентифицировать экспериментально полученные временные последовательности [2]. В статье представлены результаты количественной оценки динамических характеристик и параметризации экспериментальных временных рядов, получаемых при образовании диссипативных структур в гомогенной системе био субстрат-оксигенированные комплексы переходных металлов, на примере окисления цистеина (Cys, R) в присутствии координационных соединений железа (II) с диметилглиоксимом (DMG), гистидином (His) и молекулярным кислородом. Экспериментальная часть Методики приготовления растворов и проведения измерений описаны нами в [3]. В отличие от методики, приведенной в [3], наблюдение за динамическим поведением системы проводили путем синхронной регистрации изменения потенциала (с точностью +0,01 мВ) тремя точечными платиновыми электродами (Pt1, Pt2 и Pt3) относительно хлорсеребряного электрода во времени без перемешивания раствора (при перемешивании реакционной смеси осцилляции потенциа ла не наблюдают ся). Схема экспериментальной установки представлена на рис. 1 (расстоя-

ние между соседними индикаторными электродами составляло 18 мм). Выбор метода исследования обусловлен тем, что значения потенциала определяются соотношением окисленных и во сстановленных форм компонентов реакционной сме си. Наблюдая изменение этих величин, можно отслеживать функцию состояния системы по изменению концент рации ре агентов во времени. Характерные кривые зависимости изменения потенциала во времени представлены на рис. 2. Обнаружено, что флуктуации потенциала наблюдаются при следующих условиях проведения эксперимента: СCys = (1,25-7,50)ћ10-3 моль/л, Сcat = (2,5- 10,5)ћ10-5 моль/л, Т = 40-70њС и pH 7,30-7,50 (СCys и Сcat - концентрации реагента и оксигенированного комплекса соответственно). Полученные результаты свидетельствуют о том, что жидкофазное окисление цистеина в указанных условиях протекает в колебательном режиме, причем наблюдаемые о сцилляции носят как временной, так и пространственно-временной характер, т.е. имеет место образование диссипативных ст руктур. На это указывает тот факт, что значения потенциала в каждой фиксированной временной точке трех зависимостей не совпадают (рис. 2). Наряду с этим экспериментально получено, что флуктуации потенциала наблюдаются и при проведении измерений в случае, когда индикаторные электроды и электрод сравнения находятся в одной ячейке, а также при проведении измерений с использованием в каче ств е индикаторного элект рода и элект рода


444

ВЕСТН. МОСК. УН-ТА. СЕР. 2. ХИМИЯ. 2011. Т. 52. ? 6

Рис. 1. Схема потенциометриче ской установки: 1 - жидкостный термостат, 2 - термометр, 3 - электролитический ключ, 4 - точечные платиновые электроды (Pt1, Pt2 и Pt3), 5 - мультитест (иономер) ИПЛ-111-3, 6 - компьютер

сравнения точечных платиновых электродов, помещенных в исследуемый раствор (известно, что, если пространственные критические явления не имеют места, то изменение потенциала в указанных случаях должно быть равно нулю). Обсуждение результ атов Анализ динамики протекающих процессов При интерпретации экспериментальных результатов, полученных в виде временных рядов, основными задачами являются количественная оценка характеристик динамических систем, а также установление возможности динамического прогнозирования и идентификации ('паспортизации') временных рядов. В настоящей работе для определения количественных параметров динамики протекающих химических процессов мы проводили дискретное преобразование Фурье (ДПФ), реконструкцию динамики по временным последовательностям с определением размерностей фазового пространства и аттрактора, вычисление характеристических показателей Ляпунова (i) и энтропии Колмогорова-Синая (КС-энтропия, h). Идентификацию временных рядов осуществляли при помощи подходов фликкер-шумовой спектроскопии (ФШС) - метода анализа динамических систем, успешно развивающегося в последнее время и позволяющего провести параметризацию реальных хаотических сигналов, продуцируемых открытыми системами [4]. Дискретное преобразование Фурье Результаты ДПФ показали, что при разных условиях эксперимента число частот до статочно велико, и выделить из них специфические не удается, что сви-

детельствует о реализации хаотического режима колебаний. Следует обратить внимание на различимость спектров Фурье каждого из трех синхронно регистрируемых временных рядов, полученных при регистрации тремя электродами (Pt1, Pt2 и Pt3), что подтверждает возможность реализации в исследуемых системах детерминированного хаоса пространственновременного характера. Построение фазовых портретов При анализе полученных результатов использована идея реконструкции аттрактора, с помощью которой можно восстановить фазовое пространство и иметь представление об особенностях динамики системы по изменению одной переменной [4]. В соответствии с

Рис. 2. Изменение потенциала (Е) во времени (t) (СCys = 6,25ћ10-3 моль/л, Сcat = 2,5ћ10-5 моль/л, Т = 50њС, pH 7,40): 1 - Pt1, 2 - Pt2, 3 - Pt3


ВЕСТН. МОСК. УН-ТА. СЕР. 2. ХИМИЯ. 2011. Т. 52. ? 6

445

этим, на основе экспериментальных временных рядов построены фазовые портреты в трехмерном пространстве в координатах Еt+2ф - Еt+ф - Еt, где ф = 1 с - временная задержка (рис. 3). Вид полученных кривых (рис. 3) свидетельствует о том, что траектории движения динамических систем ограничены областью фазового пространства (аттрактором), поэтому процессы, протекающие в исследуемых динамических системах, носят детерминированный характер. Сравнение вида фазовых портретов показывает различие траекторий движения динамических систем, что подтверждает факт реализации диссипативных структур в исследуемой системе. Определение размерностей фазового пространства и аттрактора Анализ особенностей динамики изучаемых процессов проведен с помощью интегральной корреляционной функции аттрактора [5], представленной в виде где (x) - функция Хевисайда ((x) = 1 при x > 0; (x) = 0 при x 0). Эта функция позволяет на основе определения наклона линейной области зависимости ln C(r) от ln r вычислить размерности аттрактора d (нижняя граница) для разных значений размерности фазового пространства n. Зависимости в координатах ln C(r)-ln r для временного ряда, полученного при регистрации потенциала электродом Pt1 (n = 2-8), представлены на рис. 4. Размерность фазового пространства и аттрактора определяли на основании зависимости d от n (рис. 5).
Рис. 4. Зависимость ln C(r) от ln r при n: 1- 2, 2 - 3, 3 - 4, 4 - 5, 5 - 6, 7 - 8

Рис. 5. Зависимость размерности аттрактора d от размерностей пространства n

Рис. 3. Фазовые портреты в координатах Еt+2- Еt+- Еt (здесь и далее приведены зависимости для электрода Pt1)

Получено, что при одинаковых значениях размерностей фазового пространства величины размерностей аттрактора для трех синхронно регистрируемых сигналов электродами Pt1, Pt2 и Pt3 (табл. 1) различаются. Тот факт, что размерности аттракторов имеют дробные значения, свидетельствует о детерминированном характере исследуемых процессов, а соответствующие проявляемой динамике аттракторы являются странными, т.е. реализуется динамиче ский (детерминированный) хаос. Величины размерностей фазовых пространств соответствуют пяти, следовательно, при моделировании кинетических закономерностей протекающих процессов в исследуемых системах необходимо учитывать число компонентов, равное пяти, т.е. для описания кинетических закономерностей исследуемых процессов необходима система из пяти дифференциальных уравнений.


446

ВЕСТН. МОСК. УН-ТА. СЕР. 2. ХИМИЯ. 2011. Т. 52. ? 6

Таблица 1 Характеристики динамики процесса окисления цистеина, полученные разными методами анализа временных рядов
Электрод Pt1 Pt2 Pt3 Размерность d 2,60 2,65 2,52 n 5 5 5 1 0,18 0,15 0,17 Показатели Ляпунова 2 0 0 0 3 -0,35 -0,38 -0,41 h 0,18 0,15 0,17 t, c 5,6 6,7 5,9

Вычисление характеристических показателей Ляпунова и энтропии Колмогорова-Синая Важная о собенно сть хаотиче ского движения - чрезвычайная чувствительность траектории к начальным условиям. Экспоненциальную расходимость-сходимость фазовых траекторий системы оценивают при помощи вычисления характеристических показателей Ляпунова ( i), являющихся количественной мерой при определении неустойчивости [6]. Результаты расчетов с использованием некоммерческой программы TISEAN 2 [7] по экспериментальным данным (рис. 2) показателей Ляпунова (1, 2, 3), энтропии Колмогорова-Синая (КС-энтропии, h) и времени (t = 1/h), на которое можно предсказать поведение системы, приведены в табл. 1, а график зависимости показателей Ляпунова от длины временного ряда - на рис. 6. П олученные ре зульт аты показывают ( табл . 1, рис. 6), что в трехмерном фазовом пространстве при значениях 1 > 0, 2 = 0 и 3 < 0 траектории движения динамической системы сходятся к странному аттрактору, т.е. наблюдается детерминированный хаос.

Существенно, что сумма всех ляпуновских показателей (1, 2, 3) отрицательна, что, в свою очередь, является одним из о сновных признаков проявления диссипативности. Подтверждением этого является и то, что для всех случаев h > 0. Данные табл. 1 отражают и то обстоятельство, что описание динамики протекающих процессов путем обработки временных рядов указанными методами приводит к сходным результатам. Это подтверждает правомочно сть обсуждаемых подходов для анализа полученных экспериментальных данных. Анализ динамики на основе фликкер-шумовой спектроскопии В последнее время получил развитие метод фликкер-шумовой спектроскопии, который позволяет придать информационную значимо сть нерегулярно стям анализируемых сигналов - всплескам, скачкам, изломам производных различных порядков на каждом про странственном, временном или энергетиче ском уровнях иерархической организации исследуемых систем [2, 8]. При описании свойств каждого из указанных типов нерегулярно стей анализируют спектры мощности S(f) (f - частота) временных рядов V(t) ( = 0) в виде:

и вычисляют пере ходные разно стные мо менты Ф(p)(): Ф(2)() = [V(t) - V(t + )2] . Следует подчеркнуть, что многопараметрические выражения анализируемых процессов для спектров мощности и структурных функций оказываются инвариантными для каждого из пространственно-временных уровней рассматриваемой системы [2, 8].

Рис. 6. Зависимость показателей Ляпунова от длины временного ряда: 1 - l1, 2 - l2, 1 - l3


ВЕСТН. МОСК. УН-ТА. СЕР. 2. ХИМИЯ. 2011. Т. 52. ? 6

447

На рис. 7 и в т а бл. 2 приве дены ре зул ьт ат ы ФШС-анализа динамики фиксированных изменений потенциала Е(t) трех точечных платиновых электродов (Pt1, Pt2 и Pt2) относительно хлорсеребряного (см. рис. 2). Программы для расчета любезно предоставлены профессором С.Ф. Тимашевым. Для описания 'экспериментальных' зависимостей Ф(2)() ис-

пользовали свободные параметры , H1 и T1, а для описания зависимости S(f) - параметры S(0), Т0 и n0. Результаты вычисления параметров ФШС для исследуемых временных последовательностей представлены в табл. 2. Видно, что каждый временной ряд имеет свой набор указанных параметров, и, поэтому представленные зависимости позволяют провести одно-

Рис. 7. ФШС-анализ зависимости Е(t): а - анализируемая флуктуационная область зависимости Е(t); б - спектр S(f) в области низких частот; в - зависимости: 1 - 'экспериментальная' Ф(2)(), 2 - расчетная Ф(2)(), 3 - ре зонансная составляющая Фr(2)(); г - 'экспериментальная' зависимость Ф(2)(), получаемая за вычетом резонансной составляющей Фr(2)(), расчетная хаотическая составляющая Фс(2)() Таблица 2 'Паспортные' параметры флуктуационной динамики процесса окисления цистеина

Электрод Pt1 Pt2 Pt3

S(0)10-5, мВ2ћс 84,5 68,2 32,5

T0, c 1737 2170 1039

n

0

, мВ 41,63 40,66 26,42

1 0,59 0,53 0,68

T1, c 5852 7823 1511

2,1 1,9 2,1

Примечания. S(0) - параметр, характеризующий низкочастотный предел спектра мощности, формируемой нерегулярностями-вспле сками, n0 - параметр, характеризующий скоро сть потери корреляционных связей в последовательности нерегулярностейвспле сков, происходящих на временных интервалах t0, H1 - константа Херста, характеризующая скорость, с которой динамиче ская переменная теряет 'память' о своей величине на временных интервалах, меньших времени t1, - среднеквадратичное отклонение измеряемой динамической переменной.


448

ВЕСТН. МОСК. УН-ТА. СЕР. 2. ХИМИЯ. 2011. Т. 52. ? 6

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

(электродами Pt1, Pt2 и Pt3 в одной ячейке) сигналов, что также подтверждает заключение об образовании диссипативных структур в исследуемой системе. Авторы выражают благодарность профессору С.Ф. Тимашеву за консультации и помощь при подготовке данной статьи.

Работа выполнена при финансовой поддержке РФФИ (проект ? 09-03-96526 р_юг_а).
СПИСОК ЛИТЕРАТУРЫ 1. Магомедбеков У.Г., Гасангаджиева У.Г., Гасанова Х.М., Магомедбеков Н.Х. // Рос. хим. ж. (Ж. Рос. хим. об-ва им. Д.И. Менделеева). 2009. 53. ? 6. С. 74. 2. Тимашев С.Ф. Фликкер-шумовая спектроскопия: информация в хаотических сигналах. М., 2007. 3. Магомедбеков У.Г. // ЖФХ. 2002. 76. ? 4. С. 676. 4. Takens F. // Lecture Notes in Math. Berlin, 1985. 1025. Р. 99. 5. Grasberger P., Procaccia I. //Physica. D. 1983. 9. N 1. Р. 189. 6. Ляпунов А.М. Общая задача об устойчивости движения. Черновцы, 2000. 7. Программы для обработки временных рядов TISEAN 2.1. http://www.mpipks-dresden.mpg.de/~tisean 8. Тимашев С. Ф., Поляков Ю. С., Лакеев С. Г., Мисуркин П. И., Данилов А. И. // ЖФХ. 2010. 84. ? 10. С. 1980.
Поступила в редакцию 10.06.2011

FORMATION OF DISSIPATIVE STRUCTURES AT HOMOGENEOUS OXIDATION OF BIOSUBSTRATA: EVALUATION OF DYNAMIC CHARACTERISTICS AND PARAMETERIZATION OF TIME ROWS
U.G. Magomedbekov, Kh.M. Gasanova, U.G. Gasangadzhieva, N.Kh. Magomedbekov (Dagestan state university)
On an example of liquid-phase oxidation of cysteine the possibility of formation of dissipative structures in the system of biosubstrate-oxigenated complexes of iron (II) is shown. The results of evaluation of dynamic characteristics of investigated processes (type of dynamics, dimensions of phase space and attractor, indicators of Lyapunov and the entropy of Kolmogorov-Sinai) are displayed. The results of parameterization of the received time rows based on approaches of flicker-noise spectroscopy are presented.

Key words: dissipative structures, phase space, attractor, indicators of Lyapunov, the entropy of Kolmogorov-Sinai, flicker-noise spectroscopy, biosubstratе, cysteine, oxidation.
С в е дения об ав т орах: М агом е дбеков Ух у маали Г аджиевич - док т. хим. наук, профе ссор, прорек т ор по кадровой политике и филиа лам, зав. кафедрой общей и неорганиче ской химии ДГУ (ukhgmag@mail.ru); Гасанова Хадижат Магомедовна - доцент кафедры общей и неорганиче ской химии ДГУ, канд. хим. наук; Гасангаджиева Умукусум Гусейновна- доцент кафедры общей и неорганиче ской химии ДГУ, канд. хим. наук; Магомедбеков Нажму дин Хайдарбегович - докторант кафедры общей и неорганиче ской химии ДГУ, канд. хим. наук.