Методика обробки експериментальних даних

[ виправити ] текст може містити помилки, будь ласка перевіряйте перш ніж використовувати.

скачати

Методика обробки експериментальних даних

Вся процедура обробки експериментальних даних може бути розділена на два етапи. На першому проводиться первинна обробка відомостей, отриманих при проведенні експерименту за хімічним рівноваги, з метою визначення значень констант рівноваги. До цього ж етапу відноситься і статистична обробка даних, що дозволяє провести вибракування помилкових відомостей.
На другому етапі отримані значення констант рівноваги піддаються аналізу c метою безпосереднього визначення ролі ентальпійного і ентропійного факторів у рівновазі позиційної ізомеризації досліджуваних речовин і вироблення підходів до прогнозування рівноваги перетворень споріднених структур.
4.1. Первинна обробка експериментальних даних
Первинна обробка результатів дослідження хімічного рівноваги заснована на рекомендаціях [38, 101-104] і включала в себе наступне.
Розрахунок відносин концентрацій компонентів, що характеризують константи рівноваги досліджуваних реакцій.
Аналіз вивчених відносин з метою встановлення моменту досягнення рівноваги і включення в обробку тільки рівноважних даних.
2. Виняток грубих помилок всередині серії визначень константи рівноваги. Для кожної температури дослідження серіями вважали досліди, що розрізняються між собою або складом вихідної суміші, або кількістю каталізатора. Відбраковування промахів у серіях проводили з використанням критерію - найбільшого за абсолютною величиною нормованого вибіркового відхилення:
, Де - Значення константи рівноваги; - Середнє арифметичне значення константи рівноваги в серії j; - Зміщена дисперсія ; - Число результативних визначень константи рівноваги в серії j.
Процентні точки найбільшого за абсолютною величиною нормованого вибіркового відхилення запозичені з роботи [102].
3. Розрахунок середнього арифметичного значення константи рівноваги і дисперсії відтворюваності в серіях після виключення грубих помилок:


Зіставлення дисперсій відтворюваності констант рівноваги в серіях при одній температурі. Ця стадія дисперсійного аналізу є вельми корисною, оскільки дозволяє контролювати помилки відтворюваності, що виникають на всіх етапах отримання експериментальної інформації. Перевірку рівності дисперсій відтворюваності в серіях виконували за двома критеріями: Фішера - якщо число серій дорівнювало двом і Бартлетта - коли кількість серій перевищувало два [38, 102]. Якщо нуль-гіпотеза виконувалася, то дисперсію відтворюваності обчислювали за наступною формулою:
, Де .
Для всіх вивчених в даній роботі перетворень дисперсії відтворюваності констант рівноваги в серіях були однорідні.
Розрахунок середнього значення константи рівноваги для кожної температури дослідження:

Перевірку значущості розбіжності середніх значень констант рівноваги в серіях. Для цього обчислювали дисперсію :

де m - число серій, w j - вага серії j, дорівнює числу визначень n j.
Величина характеризувала розсіювання значень щодо . Число ступенів свободи дисперсії дорівнювало . Перевірку гіпотези рівності середніх значень констант рівноваги в серіях проводили за допомогою розподілу Фішера. Якщо нуль-гіпотеза виконувалася, то обчислювали зведену дисперсію

з числом ступенів свободи . На цьому обробка закінчувалася.
Для всіх досліджених реакцій розбіжності між серійними константами рівноваги не перевищували дисперсії відтворюваності, що вказувало на відсутність систематичних відхилень.
4.2. Аналіз рівноважних даних
Отримані значення констант рівноваги аналізувалися за допомогою підходу, заснованого на послідовному виключенні з вихідної рідиннофазної константи рівноваги вкладів, обумовлених:
· Міжмолекулярними взаємодіями,
· Симетрією зовнішнього обертання молекул,
· Обертанням молекули як цілого,
· Змішанням конформерів,
· Коливальним рухом,
· Обертанням окремих груп у молекулах реагентів.
Результатом подібного виключення є перехід до газофазної константі рівноваги , Що розраховується наступним чином:
= Exp (ln K p s - D S (or) / R - D S (mix) / R - D S (vib) / R - D S (ir) / R),
де - Бессімметрійная газофазних константа рівноваги реакції, D S (or) = S o II, g, (or) - S o I, g, (or); D S (mix) = S o II, g, (mix) - S o I, g, (mix); D S (vib) = S o II, g, (vib) - S o I, g, (vib); D S (ir) = S o II, g, (ir) - S o I, g, (ir) .- відповідно, ентропійні вклади, обумовлені обертанням молекули як цілого, змішанням конформерів, коливальним рухом і загальмованим обертанням груп в молекулі, що розраховуються як різниця відповідних ентропійних складових кінцевих (II) і вихідних речовин (I ).
Аналіз констант рівноваги проводився в такий спосіб.
Шляхом виключення вкладу на симетрію молекул перебувала бессімметрійная жидкофазная константа реакції , Тут числа симетрії s формуються на підставі як симетрії зовнішнього обертання молекули, так і симетрії внутрішнього обертання тих груп, симетрія яких не може бути врахована при аналізі ентропійного внеску, обумовленого внутрішнім обертанням. Для систем, розглянутих в даній роботі, такими групами були феніленовие фрагменти.
Шляхом зняття вкладу на міжмолекулярні взаємодії розраховувалася бессімметрійная газофазних константа рівноваги реакції . Тиску насиченої пари розраховувалися методом Лі-Кеслера [50] або за експериментальними даними. Застосування до розрахунку тисків насиченої пари методики, описаною в розділі 2.1, дозволяє забезпечити похибка розрахунку не більше 10% відн. для всіх тисків, наведених у цій роботі. Питання розрахунку критичних температур і тисків (T c і P c) викладені в главі 2, ацентріческіе фактори (w) розраховувалися за рівнянням Лі-Кеслера [50].
Решта вклади вимагали залучення інформації про геометрію, енергетичні характеристики молекул і частоти коливального спектра. Для одержання такої інформації нами використовувались різні розрахункові методи. Остаточна обробка інформації та обчислення ентропійних вкладів виконувалася за допомогою програми Entropy, опис якої буде наведено у п. 4.3.
Геометрія молекули оптимізувати методи молекулярної механіки (силове поле MMX на базі силового поля Аллінджера MM2) програмою PCModel 3.2. Для оптимізації молекул біфенілів використовувалася PCModel 4.0, що володіє великими можливостями при розрахунках в p-електронних системах. Вихідною інформацією були оптимізована геометрія молекули для найбільш стійкого конформера та інформація про зміну енергії молекули при обертанні кожного з дзиг, зберігаються в окремих файлах. Для формування потенційної кривої бар'єру обертання кожного з дзиг використовувались значення потенційної енергії молекули при зміні двогранного кута між обраними зв'язками дзиги і кістяка від 0 о до 360 о з кроком 10 о, при цьому на кожному фіксованому значенні кута проводилася оптимізація геометрії молекули.
На підставі відомостей про геометрію молекули розраховувалося твір головних центральних моментів інерції I A I B I C, що є вільним членом кубічного рівняння
,

де , , , , , - Моменти інерції молекули (тут n - число атомів у молекулі; m i - маса i-го атома; x i, y i, z i - координати i-го атома в системі координат з центром, що знаходиться в центрі інерції молекули). Звідси

Надалі розраховувалася сума станів жорсткого ротатора

і внесок в ентропію, обумовлений обертанням молекули як цілого

де s - число симетрії молекули, h - постійна Планка, k - постійна Больцмана.
На підставі отриманих раніше відомостей про зміну енергії молекули при обертанні кожної з її груп розраховувався внесок в ентропію від змішування конформерів

де m - загальна кількість конформерів (у нашому випадку враховувалися всі стани, отримані при повороті дзиги від 0 до 350 про про с кроком 10 о, тобто m = 36 ∙ n, де n - число груп, що обертаються в молекулі), x i - мольна частка кожної конформації

де n - число груп, що обертаються в молекулі, m - кількість конформерів, E i - енергія молекули в даному стані рівна , Де - Початкове значення енергії, - Найменша енергія молекули, отримана при обертанні всіх можливих вовчків.
Для знаходження вкладу в ентропію, обумовленого коливальним рухом були використані розрахункові значення частот коливального спектра, отримані Напівемпіричним методом PM3 (HyperChem 5.0) для оптимізованої тим же методом геометрії молекули. Критерієм якості оптимізації служило відсутність в спектрі негативних значень частот.
Розрахунок внеску в ентропію, обумовленого коливальним рухом проводився наступним чином.

де ν i - частота з прийнятого до розрахунку набору, m - кількість частот у наборі. З повного набору частот коливального спектра виключені крутильні коливання, відповідні обертанню груп, що у розрахунку внеску в ентропію від загальмованого обертання, таким чином , Де n - число атомів у молекулі, n top - число вовчків. За відсутності надійних методик визначення крутильних коливань в спектрі нами застосовувалася наближена оцінка типів коливань з використанням режиму Animate програми HyperChem 5.0.
Інформація про геометрію молекули і потенційних кривих бар'єрів обертання дзиг використовувалася для розрахунку внеску в ентропію, обумовленого внутрішнім обертанням груп у молекулі. Ентропійний внесок визначався як

тут n - число максимумів потенційної кривої бар'єру обертання групи, s-число симетрії групи (підходи до визначення чисел симетрії груп, що обертаються були розглянуті у розділі 1), S fr - ентропія вільного обертання дзиги, - Різниця між ентропія вільного і загальмованого обертання, що визначається за таблицями Пітцер і Гуїнна [21] як функція і , Де V o - ефективний бар'єр обертання дзиги, Q fr - статистична сума по станах вільного внутрішнього обертання.
Величина ефективного бар'єру обертання приймалася рівною , Де - Залежність зміни потенційної енергії молекули від кута повороту дзиги φ. Для розрахунку V o отримані методом молекулярної механіки значення потенційної енергії молекули при заданих значеннях кута повороту дзиги описувалися за допомогою кубічного сплайна, потім отриманий сплайн інтегрувався за методом Сімпсона.
Статистична сума по станах вільного внутрішнього обертання розраховувалася як , Де I пр - приведений момент інерції дзиги, який розраховувався у відповідності з наступною процедурою [105].
Для групи, що обертається вводиться координатна система з осями x, y, z, розташованими наступним чином: вісь z збігається з віссю обертання дзиги, вісь x проходить через центр мас дзиги і перпендикулярна осі z, вісь y проходить через точку перетину осей x, z і перпендикулярна до них. Атоми дзиги, що лежать на осі z з подальшого розгляду виключаються. Далі проводиться розрахунок наступних величин: - Момент інерції дзиги щодо осі z, і - Твори моментів інерції, - Фактор незбалансованості дзиги.
Потім знаходяться направляючі косинуси осей x, y, z щодо головних центральних осей 1, 2, 3 інерції молекули. Напрямок осей вибирається таким чином, щоб обидві системи координат були або правими або лівими. При цьому має дотримуватися умова рівності одиниці визначника матриці напрямних косинусів, тобто

що може використовуватися для перевірки правильності визначення напрямних косинусів.
Наведений момент інерції розраховується наступним чином:

де . Тут r (i) - проекції на головні осі інерції молекули вектора, спрямованого з центру ваги молекули в центр координат дзиги, індекс i приймає значення 1, 2, 3 в циклічному порядку, тобто при i = 1 індекс i -1 дорівнює 3, а індекс i +1 при i = 3 дорівнює 1.
У ході аналізу для вивчення впливу кожного з ентропійних вкладів на константу рівноваги розглядався наступний ряд величин, одержуваних при послідовному виключенні з вихідної константи рівноваги K x кожного з позначених раніше ентропійних вкладів:
1. Бессімметрійная константа рівноваги в рідкій фазі -
2. Бессімметрійная константа рівноваги в газовій фазі -
3. з виключеним внеском на обертання молекули як цілого -
4. "Істотна" газофазних константа рівноваги, отримана після виключення з вкладу на змішання конформерів -
5. з виключеним внеском від коливального руху -
6. після виключення вкладу на внутрішнє обертання -
Значення ентропійних вкладів і констант рівноваги для вивчених перетворень приводиться в табл. 5.1.
Однією з цілей даної роботи є вдосконалення підходу до прогнозування хімічної рівноваги в реакціях позиційної ізомеризації алкілароматичних вуглеводнів та їх функціональних похідних. При прогнозуванні констант рівноваги з використанням рівняння
(4.1)
необхідна як можна більш достовірна інформація про ентальпії та зміну ентропії реакції. При цьому експериментальні дані за вказаними властивостями не охоплюють весь набір речовин, що представляють практичний і теоретичний інтерес. Таким чином, при прогнозуванні доводиться вдаватися до розрахункових методів.
Розрахунок ентальпії реакції, як правило, проводиться за допомогою адитивних методів з урахуванням ефектів взаємодії заступників, отриманих на основі експериментальних калориметричних даних.
Використання адитивних методів при прогнозуванні ентропії речовин більш проблематично через особливої ​​природи цієї властивості, незважаючи на те, що такі методи розвинені і застосовуються при масових розрахунках термодинамічних властивостей органічних речовин [50, 111]. Кращим, на наш погляд, підходом є розрахунок зміни ентропії реакцій у відповідності з методами статистичної термодинаміки [21, 105, 106, 112], де властивість представляється у вигляді суми внесків різних видів руху молекул:
. (4.2)
У зв'язку з тим, що в даній роботі розглядаються реакції позиційної ізомеризації, внесок поступального руху молекул ( ) І внесок електронних станів молекул ( ) Приймаються рівними нулю внаслідок рівності значень відповідних складових ентропії в ізомерних молекул. Методики розрахунку інших складових ентропії описані раніше в цьому розділі.
На закінчення слід зупинитися ще на одному питанні.
Очевидним є факт існування температурних залежностей у внесків, обумовлених змішанням конформерів, коливальним рухом і внутрішнім обертанням, причому вид і параметри цих залежностей визначаються структурою молекул. Таким чином, при проведенні аналізу виключення зазначених ентропійних вкладів призводить до зміни ентальпійного ефекту реакції при переході від до . Тобто, в загальному випадку, формування додаткової ентальпійного шкали [12, 144] з самостійним набором ефектів взаємодії заступників в молекулі видається цілком обгрунтованим.
Ці ефекти повинні коректно описувати рівновагу перетворень за умови, що для останніх виконується умова = 0.
Джерелом значень подібних ефектів служить експеримент за хімічним рівноваги. Значення цих ефектів є залежними від набору методів, за допомогою яких розраховуються ентропійні вклади в константу рівноваги, особливо вклади, обумовлені внутрішнім обертанням груп у молекулах і змішанням конформерів. Тому результативним при прогнозуванні може бути тільки спільне використання зазначених ентальпійним ефектів і всієї сукупності рекомендованих розрахункових методів.
Враховуючи те, що для позиційних ізомерів в порівняно вузькому (100-150 К) температурному інтервалі зміни з температурою невелика [21], слід очікувати невеликих змін з температурою окремих внесків і ентальпійним ефектів при переході від до . Це в ідеалі передбачає відсутність необхідності формування додаткової шкали для ефектів взаємодії заступників в молекулі.
Додати в блог або на сайт

Цей текст може містити помилки.

Хімія | Реферат
48.3кб. | скачати


Схожі роботи:
Методика обробки експериментальних даних 2
Статистичні способи обробки експериментальних даних
Аналіз пакетів обробки експериментальних даних SABR і BOOTSTRAP
Аналіз експериментальних даних
Застосування систем лінійних рівнянь для апроксимації експериментальних даних
Методи аналізу та обробки даних
Технічні засоби обробки даних
Структури та алгоритми обробки даних
Системи реєстрації та обробки даних
© Усі права захищені
написати до нас