Научная статья на тему 'Применение факторного анализа для оценки влияния геометрических параметров на приведенный момент инерции молекул 2-нитро-3-метилфенола'

Применение факторного анализа для оценки влияния геометрических параметров на приведенный момент инерции молекул 2-нитро-3-метилфенола Текст научной статьи по специальности «Математика»

CC BY
34
11
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ФАКТОРНЫЙ АНАЛИЗ / МОЛЕКУЛА / ГЕОМЕТРИЧЕСКИЕ ПАРАМЕТРЫ / ПРИВЕДЕННЫЙ МОМЕНТ ИНЕРЦИИ

Аннотация научной статьи по математике, автор научной работы — Куркчи Э.У., Шейх-Заде М.И.

Метод факторного анализа использован для оценки влияния геометрических параметров молекулы на величину приведенного момента инерции I пр на примере 2-нитро-3-метилфенола в цис -конформации. Получена математическая модель такого влияния в виде уравнения регрессии, из анализа которого сделан вывод, что доминирующее влияние на значение I пр этой молекулы оказывает геометрия волчка (О-Н-группа), а влияние геометрии заместителей в бензольном кольце оказывается заметно менее существенным.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Применение факторного анализа для оценки влияния геометрических параметров на приведенный момент инерции молекул 2-нитро-3-метилфенола»

Ученые записки Таврического национального университета им. В. И. Вернадского Серия «Биология, химия». Том 22 (61). 2009. № 4. С. 317-321.

УДК 539.194

ПРИМЕНЕНИЕ ФАКТОРНОГО АНАЛИЗА ДЛЯ ОЦЕНКИ ВЛИЯНИЯ ГЕОМЕТРИЧЕСКИХ ПАРАМЕТРОВ НА ПРИВЕДЕННЫЙ МОМЕНТ ИНЕРЦИИ МОЛЕКУЛ 2-НИТРО-3-МЕТИЛФЕНОЛА

Куркчи Э.У., Шейх-Заде М.И.

1РВУЗ «Крымский инженерно- педагогический университет, Симферополь, Украина,

e-mai: csepu@gala. net

Метод факторного анализа использован для оценки влияния геометрических параметров молекулы на величину приведенного момента инерции 1пр на примере 2-нитро-3-метилфенола в цис-конформации. Получена математическая модель такого влияния в виде уравнения регрессии, из анализа которого сделан вывод, что доминирующее влияние на значение 1пр этой молекулы оказывает геометрия волчка (О-Н-группа), а влияние геометрии заместителей в бензольном кольце оказывается заметно менее существенным. Ключевые слова: факторный анализ, молекула, геометрические параметры, приведенный момент инерции

ВВЕДЕНИЕ

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

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

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

Таблица 1.

Геометрические параметры молекулы фенола

г(С-С), А г(С-Н), А г(С-О), А г(О-Н), А ZССС, град ZССН, град ZСОН, град вч град литерату ра

1,3933 1,083 1,374 0,957 119,2120,9 119,2121,6 108,8 2,52 [3]

1,397 1,0761,084 1,364 0,956 120 120 109,0 2,2 [4]

- - 1,369 0,941 - - 107,5 4,833 [5]

1,311,44 - 1,41 - 117-124 - - 0 [6]

- - - - - - - 5 [7]

1,3954 1,082 1,379 0,944 - - 105,4 - [8]

в ) - угол между связью С-О и осью Z, проходящей через атомы С1 и С4.

МАТЕРИАЛЫ И МЕТОДЫ

Для решения поставленной задачи удобно воспользоваться методами планирования эксперимента (факторного анализа) [9], с помощью которых можно получить математическую модель влияния геометрических параметров на значение 1пр. В качестве объекта исследования выбран 2-КО2-3-СН3- фенол, так как полученные в данной работе результаты предполагается использовать в дальнейшем при определении вида потенциальной функции внутреннего вращения этой молекулы. Данные, имеющиеся в литературе (например, [2]), указывают на то, что значение 1пр существенно зависит от геометрии волчка. Кроме того, априори можно предположить, что значение 1пр будет зависеть и от угла 0 выхода плоскости нитрогруппы из плоскости бензольного кольца. Согласно [10] для 2-метил-нитробензола угол 0 = 30 Исходя из сказанного, в качестве параметров варьирования выбраны следующие величины: Zp, ZСОН, ze, Z0N0, г(С-О), г(О-Н), г(С-^).

РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ

Введем обозначения: х1= zp, х2= ZСОН, х3= г(С-О), х4=г(0-Н), х5=ze, х6= г(С-К), х7=Z0N0, ,у=1пр. Величины X; будем называть факторами. Ограничимся рассмотрением линейной модели, которая может быть представлена уравнением

7

регрессии у = Ь0 + ^ bi - хг (1)

г=1

Для нахождения коэффициентов Ь0, Ъ (г=1, 2, ..., 7) была построена в кодированных величинах матрица планирования типа 27-4 (1/16-реплика) с генерирующими соотношениями х4= - хг х2, х5= - ху х3, х6= - х2- х3, х7=х1- х2- х3, которая представлена в табл. 2.

Нулевой уровень для каждого фактора выбирался таким образом, чтобы натуральное значение фактора на нижнем уровне равнялось значению этого фактора для равновесной конфигурации молекулы. Интервал варьирования каждого фактора выбирался таким, чтобы охватить разброс значений этого фактора, представленный в доступной нам литературе. Расчеты проводились для цис-формы 2-К02-3-СН3-фенола. Геометрия равновесной конфигурации принималась как для фенола [3], расстояния С-К и геометрия N02- группы принимались как для нитробензола [11], геометрия СН3-группы и расстояние С3-С принимались как в толуоле [12]. Значения у находились по формулам, приведенным в [1] с использованием рабочей матрицы в именованных значениях факторов. Полученные значения у приведены в табл.2.

Таблица 2.

Матрица планирования №1 типа 27-4 в кодированных величинах

Хо Х1 Х2 Х3 Х4 Хз Хб Х7 у-103 а.е.м.- А

+ - - - - - - - 752,449

+ + - - + + - + 940,12б

+ - + - + - + + 80б,533

+ + + - - + + - б82,444

+ - - + - + + + 755,189

+ + - + + - + - 975,3бб

+ - + + + + - - б98,280

+ + + + - - - + 72б,291

Рассмотрим какие из факторов х^ оказывают доминирующее влияние на значение у. Для этого необходимо вычислить значения коэффициентов Ъ в уравнении (1). Матрица планирования, представленная в табл.1, позволяет получить совместные оценки коэффициентов регрессии

Ъ' = А - #24 - Лз - Рб7, Ь2 = #2 - #4 - Аб - #57, Ь = Л - Аз - Аб - #47,

Ь'4 = #4 -А2 -#37 - Лб, Ь5 = А -Л13 - #27 - Дб, Ь'б = Л -#17 -#23 -#45,

Ь7 = #7 - #16 - #25 - #34.

Для того, чтобы оценить линейные члены отдельно от парных взаимодействий, была реализована еще одна матрица планирования типа 27-4 (1/16-реплика) с генерирующими соотношениями: х4=хг х2, х5= х^ х3, хб= х2- х3, х7=хг х2- х3. Такая матрица позволяет получить совместную оценку для коэффициентов Ъ", в которых все парные взаимодействия имеют знаки плюс. Тогда раздельные оценки для линейных членов можно получить путем усреднения первой и второй системы оценок: Ъ1 = (Ъ/+ Ъ")/2. Полученное таким образом уравнение регрессии для рассматриваемой задачи имеет вид:

1б2

y = (798,267 + 25,443 • x - 61,405 • x2 + 2,706 • x3 + 69,790 • x4 - 7,744 • x5 + 6,266 • x6 + + 7,023 • x7) •lO-3 (2)

Для оценки значимости коэффициентов уравнения (2) нужно знать величину дисперсии воспроизводимости S2^. Данная работа представляет по существу вычислительный эксперимент [13]. Специфика вычислительного эксперимента такова, что значение S2^ оказывается равной нулю. Это связано с тем, что в каждой строке матрицы планирования уи = уи]-, где и- номер строки матрицы, уи - среднее арифметическое значение у в и-той строке, уи]- - результат отдельного вычисления в

и-той строке,]- номер отдельного вычисления в и-той строке (]=1, 2, ..., п), п- число вычислений у в и-той строке. Поэтому ошибки АЬг- оказываются равными нулю и все коэффициенты в уравнении (2) нужно признать значимыми. Так как использованные реплики являются насыщенными, то число степеней свободы f=0 и нет возможности строгой проверки адекватности полученной модели с помощью критерия Фишера [2]. Оценку адекватности можно провести косвенным образом путем сравнения значений уи, полученных по формулам, приведенным в [1] и значений уи, вычисленных по уравнению [2]. Например, значения у1 и y1, т.е.для равновесной конфигурации рассматриваемой молекулы, оказываются равными 752,44940-3 а.е.м.- Ä2 и 756,189 10-3 а.е.м.- Ä2 соответственно (разница составляет

0.5.%), что является вполне приемлемым.

Из уравнения (2) видно, что доминирующее влияние на значение у оказывают факторы х4, х2, х1, т.е. геометрия волчка (ОН-группа), причем увеличение значений факторов х1, х4 приводит к росту величины у, а увеличение значения фактора х2 приводит к уменьшению величины у. Влияние остальных факторов оказывается заметно менее существенным.

ВЫВОДЫ

1. Получена математическая модель влияния геометрических параметров на значение 1пр в ^ис-форме 2-нитро-3-метилфенола.

2. Показано, что доминирующее влияние на 1пр указанной молекулы оказывает геометрия волчка (ОН-группа), а влияние геометрии заместителей в бензольном кольце оказывается заметно менее существенным

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

1. Марголин Л.Н. Вычисление приведенных моментов инерции для внутреннего вращения в симметричных молекулах / Л.Н. Марголин, Ю.А. Пентин, В.И. Тюлин // Опт. и спектр. - 1973. -Т. 35.- № 5.-С. 824-827.

2. Quade C.R. Internal rotation in completely asymmetric molecules. III Theory for molecules with twofold potential barriers / C.R. Quade // J.Chem^s. - 1967. -V. 47. - № 3. - P. 1073-1090.

3. Larsen N.W. Microwave spectra of the six mono- 13C- substituted phenols and some monodeuterated species of phenol. Complete substitution structure and absolute dipole moment / N.W. Larsen // J.Mol.Struct. -1979. - V. 51. - № 2. - Р. 175-190.

4. Pedersen T. Microwave spectra of the six monodeuteriophenols. Molecular structure, dilope moment and barrier to internal rotation of phenol / T. Pedersen, N.W. Larsen, L. Nygaard // J.Mol.Struct. -1969. -V. 4. - № 1.-Р. 59-77.

5. Quade C.R. Contributions of internal rotation to the rotational coefficient in phenol / C.R. Quade // J.Chem.Phys. - 1968. - V. 48. -№ 2. - P. 5490-5493.

6. Китайгородский А.И. Строение органического вещества. Данные структурных исследований. 1929 -1970 / А.И. Китайгородский, П.М. Зорький, В.К. Бельский. - М.: Наука, 1980. -648 с.

7. Bist H.D. The vibrational spectrum and torsion of phenol / H.D. Bist, J.C.D. Brand, D.R. Williams // J.Mol.Spectrosc. -1967. - V. 24. -№ 4. - Р. 402-412.

8. Forest H. Microwave spectra of some isotopically substituted phenol / H. Forest, B.P. Dailey // J.Chem.Phys. - 1966. - V. 45. - № 5. - P. 1736-1746.

9. Тихомиров В.Б. Планирование и анализ эксперимента / В.Б. Тихомиров. - М.: Легкая индустрия, 1974.- 264 с.

10. Шорыгин П.П. О зависимости спектров ароматических нитросоединений от угла поворота нитрогруппы вокруг связи С-N / П.П. Шорыгин, З.Ф. Ильичева // Изв.АН СССР. Сер.физ. -1958. -Т. 22.- № 9.- С. 1058-1062.

11. Вилков Л.В. Определение геометрического строения свободных молекул / Л.В. Вилков, В.С. Мастрюков, Н.И. Садова. -Л.: Химия, 1978.-224 с.

12. Свердлов Л.М. Колебательные спектры многоатомных молекул / Л.М. Свердлов, М.А. Ковнер, Е.П. Крайнов. - М.: Наука, 1970. - 560 с.

13. Основы научных исследований/ Под ред. Крутова В.И. и Попова В.В. - М.: Высшая школа. -1989. - 400 с.

Куркчи Э.У. Застосування факторного анаизу щодо оцшки впливу геометричних napaMeTpiB на приведений момент шерцй молекул 2-штро-3-метилфенолу / Э.У. Куркчи, М.И. Шейх-Заде //

Вчеш записки Тавршського нащонального ушверситету iм. В.1. Вернадського. Сeрiя „Бюлопя, хiмiя". - 2009. - Т. 22 (61). - № 4. - С. 317-321.

Метод факторного аналiзу застосовано для оцшки впливу геометричних параметрiв молекули на величину приведеного моменту шерцп 1пр на прикладi 2-штро-3-метилфенолу в цг'с-конформаци. Одержана математична модель такого впливу у виглядi рiвняння регреси, з аналiзу котрого зроблено висновок, що домшуючий вплив на значення 1пр ще! молекули оказуе геометрiя дзиТи (О-Н-група), а вплив геометрп замюниюв у бензеновому кшьщ е помино менш впливовим.

Ключовi слова: факторний аналiз, молекула, геометричш параметри, приведений момент iнерцi!..

Kurchi E. The application of factor analysis for the estimation of the influence of geometric parameters on the reduced moment of inertia in 2-nitro-3-methylphenol molecules / E. Kurchi, M.I. Sheikh-Zade //

Scientific Notes of Taurida V.Vernadsky National University. - Series: Biology, chemistry. - 2009. - V.22 (61). - №4. -Р. 317-321.

The method of factor analysis has been used for the estimation of the influence of geometric parameters on the value of the reduced moment of inertia Ir on the example 2-nitro-3-methylphenol in cis- conformation. The mathematical model of such influence in the form of regression equation has been obtained, from the analysis of which a conclusion has been made that a dominant influence on Ir value of this molecule makes the top geometry (OH-group), but the influence of the geometry substitutes in benzene ring is noticeably less essential.

Keywords: factor analysis, molecule, geometric parameters, reduced moment of inertia.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.