ПРАВИТЕЛЬСТВО РОССИЙСКОЙ ФЕДЕРАЦИИ
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ
УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ
«САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ»
(СПбГУ)
Кафедра Молекулярной биофизики и физики полимеров
Направление «Прикладные математика и физика»
Конформация сверхразветвленных полимеров:
экспериментальные данные и компьютерные расчеты
диплом на соискание академической степени бакалавра физики
Бакалаврская работа студентки
дневного отделения
___________________ Лазаренко Дарьи Сергеевны
___________________
Научный руководитель:
к.ф.-м.н., Е. В. Лебедева
___________________
Рецензент:
к.ф.-м.н. М. А. Безрукова
Санкт-Петербург
2016
СОДЕРЖАНИЕ
1. Введение
2
Цели
7
Объекты исследования
7
2. Основные методы изучения полимеров
9
1.1.
Вискозиметрия
9
1.2.
Прибор для измерения вязкости
10
2.1
Денситометрия
11
2.2
Прибор для изучения плостности
11
3.1
Динамическое рассеяние света
12
3.2
Прибор для измерения динамического рассеяния света
13
4.
Другие методы
14
3. Эксперимент
16
Обработка экспериментальных данных
16
4. Моделирование
24
5. Выводы
32
6. Благодарности
34
7. Список литературы
35
1
ВВЕДЕНИЕ
Одним из самый востребованных направлений в науке о современных
материалах, в том числе и о полимерах, является исследование
высокофункциональных макромолекул сложной структуры. К таким объектам
относятся высокомолекулярные соединения с большим количеством концевых
функциональных групп, а именно дендримеры, гребнеобразные полимеры,
сверхразветвленные полимеры и др. В данной работе исследуются образцы
сверхразветвленных полимеров. В настоящее время большой акцент делается
на их синтезе, изучению свойств и поиску новых путей применения в жизни и
промышленности. Это связано с тем, что сверхразветвленные полимеры
обладают уникальной топологией по сравнению с линейными или сшитыми
аналогами идентичного состава. Главной особенностью сверхразветвленных
полимеров является "прогрессивный" характер ветвления. Структура таких
полимеров не является регулярной благодаря наличию циклов, разному
расстоянию между узлами ветвления и наличию "засохших" ветвей (т.е. ветвей,
на которых не возможен дальнейший рост).
Развиваясь из точечного разветвляющего центра путем непрерывных
ветвлений, подобные объекты уже на ранних стадиях роста приобретают
глобулярную форму, которая очень мало зависит от внешних условий.
Основные их качества – это наноразмерная структура, близкая к сферической
форма, высокое количество терминальных групп и разветвленная структура,
которые определяют следующие физико-химические свойства: высокая
растворимость и термическая стабильность, низкая вязкость растворов,
способность к образованию композитов, и применение в методах разделения и
концентрирования [1]. Наличие в своем составе большого количества концевых
функциональных групп делает возможным дальнейшие изменения
сверхразветвленных полимеров для придания им свойств, требуемых при
создании различного рода материалов [2]. Концевые группы обладают
2
достаточной свободой локализации, легко модифицируются при воздействии
окружающей среды. Сложность исследования данных полимеров заключается в
корректности изучения изменения структуры этих объектов под действием
внешних факторов и взаимодействия функциональных групп с этими
условиями среды. Корректный учет взаимодействий групп друг с другом и
молекулами раствора, их влияние на подвижность и структуру исследуемых
объектов вс е еще являет ся нерешенной проблемой. Подробный
экспериментальный анализ и построение теорий для описания таких
соединений является сложной процедурой из-за нетривиальной архитектуры
молекул. Из-за древообразной структуры сверхразветвленных полимеров,
концевые группы могут составлять половину от всей молекулы. Степень их
ионизации в водных растворах зависит от pH, в результате чего, при
определенноых условиях, сверхразветвленные полимеры могут образовывать
интерполиэлектролитные комплексы с противоположно заряженными
молекулами. С помощью этого свойства структуры такого типа можно
использовать как молекулярные наноконтейнеры для доставки лекарств и генов
в клетки организма, а также в качестве пищевых добавок и в полиграфической
промышленности [3]. В общей сложности, сверхразветвленные полимеры
нашли широкое применение во многих областях, а именно полимерном
материаловедении, технике, медицине, биологии, микроэлектроники и др.
Благодаря преобразованию концевых групп можно контролировать
реакционную способность, растворимость, адгезию к поверхности,
люминесцентные и электрохимические свойства древообразных структур [4].
Использование дендритных полимеров в хроматографиче ских и
электрофоретических методах на данный момент позволило решить ряд
аналитических задач на уменьшение сорбции [5], [6].
Сверхразветвленные полимеры различного строения обычно синтезируют
«традиционными» методами полимерной химии - поликонденсацией и
полимеризацией. Широко используются также различные типы реакций
3
полиприсоединения, полимеризации с раскрытием кольца. Но, как показал
анализ имеющихся публикаций, несмотря на то, что при помощи реакций
полициклоприсоединения можно синтезировать сверхразветвленные полимеры
с самыми разнообразными и необычными свойствами и с такой структурой,
которую практически невозможно получить при помощи других методов [2].
Сверхразветвленные молекулы, по сравнению с дендримерами, заметно
сокращают затраты ресурсов и временя синтеза. Тем не менее, структура таких
полимеров не контролируема как у дендримеров, и функциональные группы не
расположены в упорядоченном положении. Это происходит из-за
статистического характера стадий связывания, стерических затруднений в росте
цепей, и способности к реагированию функциональных групп, так как
распространение происходит только с двух сторон между блоками ветвления,
что производит нерегулярные и линейные сегменты. Таким образом, степень
ветвления для дендримеров и разветвленного полимера сильно отличается, хотя
по всей композиции они аналогичны. Из-за несовершенной конструкции,
полидисперсности молекулярной массы, отсутствия упорядоченного
расположения функциональных групп и неконтролируемой формы
сверхразветвленных полимеров, они могут быть использованы в качестве
функциональных материалов в развивающихся областях, таких как материалы
для покрытий и модифицирования. В связи с этим, сверхразветвленные
полимеры могут иметь промежуточные свойства между теми, которые имеют
линейные и дендритные полимеры [7]. Часто эти молекулы синтезируются в
одном сосуде без длительных стадий пошаговой реакции и очистки,
необходимых для получения «традиционных» дендримеров. Несмотря на
большую полидисперсность и дефекты их химической структуры, вызванные
внутренней циклизацией и побочными реакциями, сверхразветвленные
полимеры обладают, в значительной степени, аналогичными характеристиками
компактных структур со значительной долей концевых групп, расположенных
на внешней стороне молекулы [8]. Несмотря на их нерегулярную химическую
4
структуру, множественные слабые межмолекулярные взаимодействия могут
способствовать сверхразветвленной полимерной сборке хорошо упорядоченных
одномерных микроскопических или даже макроскопических волокон [9], [10].
Синтетические, модифицированные сверхразветвленные полифенелены могут
быть использованы в качестве макроинициаторов (катализаторов) для синтеза
звездообразных полимеров путем анионной полимеризации [11].
Степень разветвленности сверхразветвленных полимеров принимает
значения от 0 до 1, где предельные значения соответствуют образованию
линейного полимера и дендримера соответственно. Как было описано выше,
контролировать степень разветвленности сложно, этот процесс нерегулярен.
Нерегулярность приводит к образованию большого числа геометрических
изомеров, которые растут с увеличением молекулярного веса. [12]. Оказывается,
что такие изомеры также влияют на свойства сверхразветвленных полимеров.
Одним из наиболее интересных аспектов в сверхразветвленных полимерах
является получение сверхразветвленных полимеров с узким молекулярномассовым распределением и степенью разветвленности близкой к 1.
Полимеры имеют переменный гидродинамический радиус в зависимости
от свойств растворителя. Гидродинамические радиусы сверхразветвленных
полимеров имеют меньшие размеры, чем у их линейных аналогов с той же
молярной массой. Это указывает на то, что молекулярные массы дендримеров и
сверхразветвленных полимеров, определенных с помощью гельпроникающей
хроматографии следует рассматривать с некоторым скептицизмом.
Гидродинамические радиусы были также чувствительны к полярности
функциональных групп на периферии [7].
Влияние концевых групп сверхразветвленных полимеров на тепловые
свойства было детально изучено Хоукеру и Чу [13], которые подготовили поли
(эфиры кетонов) с одинаковым строительным блоком, но с разными степенями
разветвленности и конечными группами. Они продемонстрировали, что
5
тепловые свойства сверхразветвленных полимеров с различной степенью
разветвленности не зависят от макромолекулярного строения, но сильно зависят
от свойств концевых групп. Ким и Вебстер [14] предположили, что стеклование
сверхразветвленных полимеров обусловлено поступательным движением
молекулы, а не посегментным движением цепи. Они также отметили, что
свойства в большей степени зависят от длины цепи между точками ветвления, а
не от молекулярной массы, иными словами, сверхразветвленные полимеры с
более короткой длиной цепи ведут себя так же, как олигомеры, даже если имеют
большие молекулярные веса.
Следует отметить, что для успешного применения сверхразветвленных
полимеров с низкой полидисперсностью в тех или иных областях, необходимо
постоянно совершенствовать методы их синтеза, не смотря на существующий
уровень синтетических технологий, достигнутый за последнее время. Также,
это позволит провести систематическое и фундаментальное исследование
взаимосвязи между структурой и функциональностью полимера. [7].
Данная работа вносит вклад в исследование сверхразветвленных
полимеров и дает возможность посмотреть еще с одной стороны на изучение
этих структур, а именно, при помощи компьютерных ресурсов.
6
Целями данной работы являются:
- детальное исследование свойств сверхразветвленных пиридилфениленовых
полимеров в растворе методами молекулярной гидродинамики (денситометрия,
вискозиметрия, скоростная седиментация)
и оптики (динамическое
светорассеяние и двойное лучепреломление в потоке),
- сопо ставление гидродинамиче ских характеристик, полученных
экспериментально и из компьютерных расчетов.
Объекты исследования
Сверхразветвленные пиридилфениленовые полимеры, исследуемые в
данной работе, имеют структуру, содержащую пиридиновые и бензольные
кольца. Химические структуры исследуемых объектов представлены на рис. 1.
F-551
F-536
F-547
Рис. 1. Химические структуры исследуемых пиридилфениленовых полимеров.
Объекты были синтезированы в Институте элементоорганических
соединений им. А. Н. Несмеянова РАН, г. Москва [15].
В работе изучены физико-химические характеристики сверхразветвленных
пиридилфениленовых полимеров с различными фокальными группами в
растворе. В качестве растворителя был выбран тетрахлорэтан (ТХЭ, вязкость
растворителя η0 = 1.609 cП; плотность ρ0 = 1.587 г/см3; показатель преломления
7
n0 = 1.4916). Изучение гидродинамических и оптических свойств исследуемых
растворов полимеров выполняли при температуре Т = 250С.
8
ОСНОВНЫЕ МЕТОДЫ ИЗУЧЕНИЯ ПОЛИМЕРОВ
Существует ряд экспериментальных методов, которые позволяют
определить физические характеристики полимеров.
В данной работе, для изучения объектов, были использованы следующие
методы: вискозиметрия, денситометрия, динамическое рассеяние света,
скоростная седиментация и двойное лучепреломление в потоке.
1.1.
Вискозиметрия
Вискозиметрия является наиболее широко применяемым транспортным
методом определения гидродинамических свойств макромолекул.
В ходе эксперимента по определению характеристической вязкости
измеряют вязкости растворителя
η0
[ η]
и ряда концентраций раствора η .
Теоретически зависимость [] от с описывается уравнением [16]:
η−η0
=[ η ]+k c [ η ]2
cη0
(1)
г д е k` - константа Хаггинса, [η] – характеристическая вязкость полимера.
Характеристическая вязкость полимера в соответствии с уравнением (1)
определяется соотношением:
η −η o
¿
ηo c
(2)
Величина характеристической вязкости [η] имеет размерность удельного
объема и позволяет определить величину молекулярной массы М полимера.
9
1.2. Прибор для измерения вязкости
Измерение вязкости проводили на микровискозиметре Lovis 2000 M/ME
(рис. 2), действие которого основано на принципе «катящегося шарика».
Данный прибор применяется для измерения динамической и кинематической
вязкости жидкостей.
Образец заливают в стеклянны капилляр, который затем помещают в
капиллярный блок с возможностью термостатирования. Этот блок может
менять угол наклона, чем расширяет диапазон измерений и подвергает образец
различным скоростям сдвига.[Руководство по эксплуатации Lovis 200 M/ME
Микровискозиметр] Для измерений использовали стеклянный капилляр с
внутренним диаметром 1.59 мм и золотой шарик 1.50 мм в диаметре.
Погрешность измерений не превышает 5%.
Рис. 2. Микровискозиметр Lovis 200 M/ME.
2.1. Денситометрия
10
Денситометрия, как метод измерения плотности, используется для
определения удельного парциального объема.
2.2. Прибор для измерения плотности
Измерение плотности выполняли на денситометре RM40. Прибор
изображен на рис. 3. RM40 – автоматический прибор для измерения плотности
растворов полимеров. Принцип действия RM40 основан на измерении частоты
ко л е б а н и й и з м е р и т е л ь н о й к ю в е т ы , в ы з ы в а е м ы х п е р е м е н н ы м
электромагнитным полем. Под действием возбуждающего поля пустая кювета
будет колебаться с собственной частотой, а при заполнении измерительной
кюветы веществом будет меняться её масса и, соответственно, частота
колебаний кюветы. Подобно маятнику, чем больше плотность образца, а значит
и его масса, заключенная в кювете, тем ниже частота колебаний. Погрешность
измерений составляет порядка 1%.
Рис.3. Денситометр RM40.
3.1. Динамическое рассеяние света
11
Метод динамического или квазиупругого рассеяния света, который часто
называют фотонной корреляционной спектроскопией или спектроскопией
о п т и ч е с ко го с м е ш е н и я и с п ол ь зуе т с я дл я и з м е р е н и я р а з м е р о в
высокомолекулярных объектов. Данный метод позволяет определить
коэффициент диффузии дисперсных частиц в жидкости путем анализа
корреляционной функции флуктуаций интенсивности рассеянного света. Далее,
из коэффициента диффузии рассчитывается гидродинамический радиус частиц.
В методе квазиупругого рассеяния света измеряется частотный спектр
рассеянного света. Используется монохроматический исходный пучок света от
лазера с интенсивностью �0, частотой �0 и длиной волны �0. Поскольку
рассеивающие объекты – макромолекулы движутся, свет, рассеянный под углом
� не монохроматичен. Регистрируется интенсивность света I, рассеянного под
этим углом.
В разбавленном растворе макромолекулы движутся как точечные
рассеиватели с коэффициентом диффузии �.
Зависимость I от частоты рассеянного света � представляет собой спектр
рассеянного света. Характеристическая ширина этой кривой Δ� ~ Dq2 где q –
волновой вектор. Таким образом, измерение спектра рассеянного света
позволяет определить коэффициент поступательной диффузии полимерных
клубков D.
Автокорреляционная функция представляет собой убывающую
экспоненту. Зная ее постоянную времени, которая называется временем
релаксации τ или временем корреляции, коэффициент поступательной
диффузии Dt можно определить по формуле:
1
2
= Dt q
τ
(3)
3.2. Прибор для измерения динамического рассеяния света
12
Универсальный цифровой коррелятор является основным измерительным
у с т р о й с т в ом ф ото н н о й ко р р е л я ц и о н н о й с п е кт р о с ко п и и ( Ф КС ) ,
обеспечивающим необходимую точность измерений. Коррелятор работает в
реальном масштабе времени.
Измерения динамического рассеяния света были выполнены на установке
Photocor Instruments. с коррелятором реального времени PhotoCor-FC (рис. 4),
содержащим 288 каналов. Образец освещали линейно поляризованным светом
𝐻� − 𝑁� лазера с длиной волны �0 = 632.8нм и мощностью 25 мВт.
Использовали стеклянные цилиндрические кюветы с диаметром 1 и 0.5 см,
которые погружали в иммерсионную жидкость – декан, имеющую близкий
показатель преломления со стеклом кюветы. Рассеянный свет регистрировали в
диапазоне углов � о т 25o д о 150o с помощью фотоэлектронного умножителя,
работающего в режиме счета фотонов [17]. Погрешность измерений составляет
около 5%.
Рис. 4. Спектрометр Photocor Complex.
4. Другие методы
13
- Метод скоростной седиментации позволяет определять отношение
молекулярной массы к коэффициенту поступательного трения (M/f), который, в
свою очередь, связан с размером макромолекул. На величину этого отношения
определенное влияние оказывает также полидисперсность полимерных
образцов. Определение молекулярного веса полимера с помощью метода
скоростной седиментации основано на установлении седиментационного
равновесия в растворах полимера. Частицы с большей массой движутся в
направлении донной части кюветы в установке со скоростью более высокой,
чем макромолекулы с меньшей массой. Распределение скоростей седиментации
зависит от градиента концентраций, который устанавливается в граничной
области, и его изменения во времени.
При седиментации поступательное движение молекул происходит под
действием центробежной силы Fc. В условиях скоростной седиментации она
уравновешивается силой трения. Записав скорость движения молекул как dx/dt,
коэффициент седиментации можно определить из соотношения:
dx
dt
M ( 1−v ́ ρ0 )
s≡ 2 =
∙
,
f
ω x NA
(4)
где x – расстояние между молекулой и осью вращения, ω – угловая скорость
вращения ротора, N A
ρ0
– число Авогадро, ( 1− v ´ ρ0 )
– фактор «плавучести»,
– плотность растворителя, v ̅ - парциальный удельный объем полимера и f –
коэффициент поступательного трения [16].
По данным наблюдений строят график изменения скорости седиментации
при различной концентрации и определяют по этому графику константу
седиментации s0 данного полимера при бесконечном разбавлении его раствора.
Погрешность метода скоростной седиментации может достигать 10 %.
14
- Двойное лучепреломление в потоке (эффект Максвелла, ДЛП) –
динамооптический метод используемый для изучения оптических,
конформационных и гидродинамических свойств молекул полимеров в
растворах, позволяющий непосредственно получать сведения о размерах и
кон ф иг урац и и растворенных макром олекул [ 16] , [ 18] . Д войно е
лучепреломление в потоке в изучаемом растворе возникает под действием
механических сдвиговых напряжений, которые создаются в ламинарном потоке
в динамооптиметре. Необходимым условием возникновения ДЛП в растворах
является асферичность формы жестких частиц или молекул, поскольку только
при выполнении этого условия возможна кинематическая ориентация частиц в
сдвиговом потоке. Величина наблюдаемого двойного лучепреломления Δn
определяется разностью показателей преломления для обыкновенного и
необыкновенного лучей, которая зависит как от характера функции
распределения частиц в потоке по углам ориентации, так и от оптических
свойств растворенных частиц.
Погрешность эксперимента не превышает 10%.
15
ЭКСПЕРИМЕНТ
Одним из важнейших вопросов, возникающих при исследовании
высокомолекулярных соединений, является определение молекулярной массы
изучаемого полимера. Обычно эта задач может быть успешно решена одним из
двух «абсолютных» методов — светорассеяния или седиментационнодиффузионного анализа полимера в растворе [16]. Однако, как метод
светорассеяния, так и метод седиментации оказываются неэффективными в
применении ко многим широко используемым в настоящее время полимерам.
Поэтому для определения молекулярных масс таких полимеров приходится
использовать некоторые гидродинамические характеристики молекул, которые
могут быть определены экспериментально с достаточной надежностью,
например, характеристическую вязкость [η].
Обработка экспериментальных данных
В результате экспериментов различными методами для исследуемых
объектов были получены значения характеристической вязкости, плотности,
коэффициентов диффузии D и констант седиментации s0, рассчитаны величины
удельного парциального объема
v ̅, гидродинамического радиуса
молекулярной массы МsD, оптических коэффициентов сдвига
Rh,
n
.
❑
характеристическая вязкость [η]
В ходе эксперимента были измерены времена прохождения шариком
полной длины капилляра в растворе t и растворителе t0, по которым рассчитали
значения относительной вязкости η/η 0. Для определения величины
характеристической вязкости построили зависимость приведенной вязкости
16
ηпр =
ηsp
c
=
(η−η0 )
cη0
от концентрации раствора с (рис. 5). Видно, что все
зависимости хорошо аппроксимируются прямыми линиями. При экстраполяции
этих зависимостей к бесконечному разбавлению были получены значения
характеристической вязкости пиридилфениленовых полимеров.
Рис. 5. Концентрационные зависимости приведенной вязкости для исследуемых
образцов 1 - 536, 2 - 547, 3 - 551 в ТХЭ.
коэффициент диффузии D и гидродинамический радиус Rh
Коэффициент поступательной диффузии был определен по результатам
эксперимента по динамическому светорассеянию.
Обработка автокорреляционных функций рассеянного света, выполненная
в программе DynaLS, позволила для разных концентраций раствора с
определить времена релаксации τ для разных углов рассеяния �. Строится
17
зависимость обратного времени релаксации 1/τ от квадрата волнового вектора
q2, которая описывается прямолинейной зависимостью, проходящей через
начало координат. Величина тангенса угла наклона этой зависимости задает
коэффициент поступательной диффузии Dс. Коэффициенты диффузии D для
каждого из объектов были получены путем экстраполяции зависимости Dc от с
к нулевой концентрации ( рис. 6).
Рис. 6 . Зависимость коэффициента диффузии Dc от концентрации раствора с
для сверхразветвленных пиридилфениленов в ТХЭ.
Если рассматривать полимерные макромолекулы как непротекаемые
сферы, то их гидродинамический радиус Rh можно определить по формуле
Стокса – Эйнштейна:
D
T
k
6 0 Rh
,
(5)
где k – постоянная Больцмана.
удельный парциальный объем v ̅
18
Удельный парциальный объем находят из значений инкремента плотности
по формуле v ̅ =
1−dρ/dc
, где
ρ0
dρ /dc
– инкремент плотности (тангенс угла
наклона зависимости ( ρ−ρ 0 ) от концентрации раствора с, где ρ – плотность
раствора)
Полученные методом денситометрии экспериментальные данные хорошо
аппроксимировались прямолинейными зависимостями, проходящими через
начало координат. Значение удельного парциального объема v̅ оказалось
одинаковым для всех образцов и составило 0,754 см3/г.
гидродинамический инвариант A0
Используя полученные значения D
и
[η]
, был рассчитан
гидродинамический инвариант:
η0 D ( M [ η ] )1/3
A0=
T
.
(6)
Для каждой системы полимер-растворитель во всей практически
доступной области молекулярных весов параметр А0 действительно является
инвариантом. Практически он не зависит от длины молекулярной цепи и
качества растворителя, а для умеренно полидисперсных образцов и от степени
их полидисперсности. Диапазон значений А0: 2.88 - 3.84 ×1010 эрг×К-1мол-1/3,
где наименьшее значение соответствует конформации макромолекулы
«палочка», а наибольшее значение – конформация «гауссов клубок» [16].
молекулярная масса МsD
По значениям коэффициентов диффузии D и седиментации s0 были
рассчитаны молекулярные массы полимеров:
19
S
МsD
RT
= D 1−v ̅ ρ .
0
(7)
Молекулярная масса, являющаяся одной из основных характеристик
любого химического соединения, приобретает в случае высокомолекулярных
соединений особую роль, т.к. она также служит мерой длины цепной молекулы.
оптический коэффициент сдвига
n
❑
Из экспериментов по двойному лучепреломлению в потоке были
получены значения оптических коэффициентов сдвига.
В ходе эксперимента измерялись величины
и N – величина углов в
равновесном и неравновесном состоянии системы и число оборотов ротора в
секунду соответственно, которые связаны с двулучепреломлением n и
градиентом скорости потока g соотношением:
Δn sin 2ϕ
=y
g
N ,
(8)
где y – постоянная, определяемая из геометрии прибора.
Оптический коэффициент сдвига растворенного полимера определяется
по формуле:
Δn ( Δn/g ) р− ра −( Δn/g ) р− ля
=
Δτ
η0 ( ηr−1)
,
где ηr = η/ η0 – относительная вязкость раствора и растворителя.
(9)
По экспериментальным данным были построены зависимости двойного
лучепреломления n от градиента скорости потока g (рис. 7). А также, Δn от
напряжения сдвига Δτ (рис. 8).
20
Рис. 7. Зависимость двойного лучепреломленияn от градиента скорости g для
образцов 1 - 536 (с=0,96 г/дл), 2 - 547(с=0,78 г/дл), 3 - 551 (с=0,79 г/дл) в ТХЭ
Рис. 8. Зависимость двойного лучепреломления Δn от напряжения сдвига Δ для
образцов 1 - 536 (с=0,96 г/дл), 2 - 547(с=0,78 г/дл), 3 - 551 (с=0,79 г/дл) в ТХЭ
21
Видно, что экспериментальные зависимости хорошо экстраполируются
прямыми линиями, проходящими через начало координат. Линейный характер
данных зависимостей позволяет сделать вывод об отсутствии в системе
большого количества крупных ассоциатов.
Все экспериментальные данные, полученные вышеперечисленными
методами гидродинамики и оптики, представлены в таблице 1.
В таблице представлены значения n
= MsD/M0
– количество
повторяющихся звеньев в структуре, где M0 – молекулярная масса одного
повторяющегося звена.
Следует напомнить, что эксперименты по ДЛП характеризуют степень
асимметрии формы изучаемых объектов. Т.е. наблюдается эффект ДЛП в том
случае, когда исследуемая макромолекула по форме отлична от сферы.
Таблица 1. Характеристическая вязкость - [η] (дл/г), коэффициент флотации
–s0 (Св), коэффициент диффузии D (см2/c), гидродинамический радиус Rh (нм),
молекулярная масса MSd (Да), молекулярная масса Mw (определенная методом
ГПХ), количество звеньев n, гидродинамический инвариант A0 (г·см2·с-2·град1
·моль-1/3), оптический коэффициент сдвига n/
∆n
(г-1·см·с2) для образцов
∆τ
сверхразветвленных пиридилфениленовых полимеров в тетрахлорэтане.
Образе
[η]
ц
s0×1 D× Rh
MSD× Mw*× n
А0×
013
10-3
1010
107
10-3
n
❑
1010
536
0.2
-1.9 2.9
4.5 46
49.6
8
3.3
45
551
0.32 -2.1 2.0
6.5 74
77
11 3.3
35
4
3.7
46
547
0.23 -1.9 3.3
3.9 41
70.2
6
5
22
* Значения масс получены из экспериментов по гель проникающей хроматографии (Институт
химии СПБГУ).
Из таблицы 1 видно, что значения оптических коэффициентов сдвига
и зу ч а е м ы х н а м и с т ру кту р отл и ч н ы от н ул я , т. е . ко н ф о рм а ц и я
сверхразветвленных пиридилфениленовых полимеров отлична от сферической.
Чтобы понять к какой конформации изучаемые объекты больше всего
приближены, было решено сопо ставить экспериментальные
гидродинамические параметры, с характеристиками, полученными из
компьютерных расчетов.
23
МОДЕЛИРОВАНИЕ
Компьютерное моделирование сложных полимерных систем позволяет
наиболее эффективно проводить теоретические исследования. Все еще
остаются открытыми вопросы, касающиеся изучения структуры и подвижности
сверхразветвленных полимеров в растворах методами компьютерного
моделирования. Например, нет ответа на вопрос, когда необходим явный учет
противоионов в таких системах, а когда возможно применение приближения
Дебая-Хюккеля для описания электростатических взаимодействий,
существенно увеличивающего масштабы исследуемых времен [3].
Исследование комплексов сверхразветвленных полимеров является важным
шагом в анализе возможности использования их в качестве молекул, действие
которых может контролироваться изменением внешней среды. Таким образом,
теоретическое исследование комплексов, образованных сверхразветвленными
полимерами, методами компьютерного моделирования является актуальным как
для создания высокофункциональных полимерных систем, так и для
дальнейшего развития фундаментальной науки о полимерах нового поколения.
Ко м п ь ю т е р н ы е р а с ч е т ы г и д р о д и н а м и ч е с к и х п а р а м е т р о в
сверхразветвленных пиридилфениленовых полимеров были получены двумя
способами.
В первом случае расчеты были выполнены для двух предельных случаев
конформаций моделей, состоящих из шарообразных бусин: «палочка» и
«сфера». Для этого, в структуре полимеров было выделено повторяющееся
элементарное звено, приближенное эквивалентной сферой, из которого в
последствии были собраны модели макромолекул.
Во втором случае моделировалась и обсчитывалась целая макромолекула,
где каждый ее атом был приближен сферой радиуса Ван-дер-Ваальса. Были
рассмотрены три разных конформации: «сфера», «палочка» и «крест».
Остановимся на методах более подробно.
24
- В первом методе моделирование происходило в вакууме. Структуру
повторяющегося звена «прорисовывали» в программе HyperChem и проводили
оптимизацию геометрии моделей звена (рис. 9).
Рис. 9. 3D модель структурного звена F-547.
Цель моделирования состояла в том, чтобы изучить всевозможные
конформации данного структурного звена. Для этого варьировались торсионные
углы вблизи центральной точки ветвления. В результате был получен набор
упакованных звеньев, что позволило рассчитать радиусы инерции
R g=
√
∑ M i r2i
∑ Mi
для повторяющегося звена каждой из исследуемых структур.
В ряде работ было показано, что отношение гидродинамического радиуса
к радиусу инерции для такого рода структур Rg/Rh ~1 [19]. Значения
гидродинамического радиуса для исследуемых объектов менялось не
значительно во всех рассмотренных конформациях.
25
Впоследствии для каждого из рассмотренных объектов повторяющийся
элемент был заменен сферой с эквивалентным коэффициентом трения
=60Rh.
В качестве моделей для расчета гидродинамических характеристик были
выбраны два предельных случая конформации полимерных структур, поскольку
их свойства уже предсказаны теорией. За основу были взяты сферическая
модель системы (формула Стокса для коэффициента поступательного трения) и
модель цепного ожерелья повторяющихся бус (теории Кирквуда и Хирста –
Штокмайера [18]). На рис. 10 представлены две предельные модели
полимерной цепи с повторяющимся звеном в виде «бусинок».
Рис. 10. Модели предельных конформаций для расчета гидродинамических
характеристик, в качестве повторяющегося элементарного звена
представлена структура F-547.
Далее для этих двух моделей были определены величины коэффициентов
диффузии, характеристической вязкости и гидродинамического инварианта, с
использованием следующих соотношений:
«Палочка»
«Сфера»
26
[
D=
[η] =
1 1
1
= 1+ (ln ( n )−1)
f n
n
kT
f
]
---kT
D = 6 π R ’ η0
2 π ( nl ) N A ln ( l / d )
45 M
η0 D ( M [ η ] )1/3
A0=
T
A 0=
2.5 V N A
M
η0 D ( M [ η ] )1/3
A0=
T
[η]=
η0 DM [ η ]
T
1
3
A 0=
η0 DM [ η ]
T
1
3
T - температура, K; f – коэффициент вращательного трения; n – количество
элементов в цепи; D – коэффициент диффузии, см 2/ с ; d = 2Rh – диаметр
элемента цепи, см; [η] – характеристическая вязкость, дл/г; l = d·n – длина
цепи, см; R’(зависит от количества звеньев в сфере и их упаковки) – радиус
«большой сферы», см; M – молекулярная масса, Дa; A0 – гидродинамический
инвариант, эрг×К-1мол-1/3;
η0
- вязкость растворителя, Пс; V – объем
«большой сферы», см3; k, N a , π - константы.
Радиус «большой сферы» определяется геометрически исходя из
количества составных элементов структуры.
Результаты компьютерных расчетов для двух моделей и результаты
эксперимента представлены в таблице 2.
В результате сопоставления характеристик, полученных путем
компьютерных расчетов и экспериментальными методами, видно, что
конформация объекта F-536, в структуре которого отсутствуют кислородные
мостики, ближе к сферической. Что касается двух других объектов, то они
имеют промежуточную конформацию между предельными случаями.
Следует отметить, что в этом методе параметры растворителя не учитываются
до момента расчета гидродинамического инварианта.
Табл. 2. Сводная таблица коэффициентов диффузии, характеристической
вязкости и гидродинамического инварианта для моделей «палочка», «сфера»
27
полученных по результатам компьютерных расчетов и из эксперимента для
всех типов исследованных структур.
Сфера
2.35
2,90
2.27
2,00
3,30
2,10
2.40
2,00
1.60
0.33
0,2
0.44
0,44
0,23
0.57
0.61
0,32
0.69
3.14
3,30
3.34
2,78
3,75
3.14
2.82
3,34
3.00
палочка
палочка
ент Эксперим
Сфера
ент Эксперим
F-551
палочка
D×107,
F-547
ент Эксперим
Сфера
F-536
см2/с
[η],
дл/г
А0×101
0, эрг·К1
·
мол-1/3
Исходя из результатов, полученных из программы HyperChem первым
спо собом, было решено прове сти изучение сверхразветвленных
пиридилфениленовых полимеров в специализированной программе HYDRO
[20], которая предназначена для поатомного обсчета структуры целиком.
Отличия второго метода заключаются в учете растворителя на этапе
моделирования структуры и возможности подбирать абсолютно любые
конформации.
- Во втором методе все расчеты интересующих нас гидродинамических
параметров, а именно коэффициента диффузии, гидродинамического радиуса,
константы седиментации и характеристической вязкости проходили
непосредственно в программе HYDRO.
28
Предварительно, в программе HyperChem, для исследованных структур
были подобраны три конформации, а именно «сфера», «палочка» и «крест» и
проведена оптимизация геометрии структуры. Прикладывая поле с
использованием метода молекулярной динамики в программе HyperChem,
молекулы выстраивались в наиболее «удобную» конформацию. Каждый
образец обрабатывался в программе 50 пс. Эта вспомогательная процедура дает
возможность узнать координаты каждого атома в молекуле.
29
Рис. 11. Типы моделей на примере F-536 в интерфейсе программы HYDRO:
А. - «палочка», Б. - «сфера» и В. - «крест».
Для получения характеристик в программе HYDRO необходимо
учитывать такие параметры как: плотность и вязкость растворителя;
температуру при которой выполняется эксперимент и удельный парциальный
объем полимера. Все эти величины были найдены ранее из эксперимента и
оказались одинаковыми для всех исследуемых структур.
Зная координаты всех элементов структуры, их количество и их Ван-дерВаальсовые радиусы [22] были рассчитаны искомые гидродинамические
параметры. HYDRO строит наиболее вероятную конформацию по теории
Кирквуда [23] и рассчитывает все интересующие константы. На рис. 11
представлены конформации для пиридилфениленового полимера F-536,
обсчитываемые в программе.
Полученные с помощью программы HYDRO гидродинамические
параметры представлены в таблице 3.
Также, с использованием найденных величин были рассчитаны
гидродинамические инварианты по формуле (6).
Н е с м от ря н а з н ач ит е льны й ра з б ро с д а нны х, получ е нны х
экспериментально и из программы HYDRO, можно оценить близость к той или
иной конформации для каждой из исследованных структур сверхразветвленных
пиридилфениленовых полимеров.
Анализ полученных данных показывает, что наиболее близкой
конформацией для объекта F-536 можно считать «сферу» (из совпадения
констант седиментации и диффузии). Если же учесть, что полученный методом
ДЛП оптический коэффициент сдвига отличен от нуля, тогда можно говорить
что при дальнейшей интерпретации данных для этой структуры стоит
использовать модель сплюснутого или вытянутого эллипсоида с относительно
небольшой асиммет рией формы. Из т аблицы т акже видно, что
30
экспериментальные характеристики для структур F-547 и F-551 лежат между
параметрами для конформаций «палочка» и «крест», вычисленными в HYDRO,
что также хорошо согласуется с экспериментальными данными метода ДЛП.
Т а б л и ц а 3 . Гидродинамические параметры сверхразветвл енных
пиридилфениленовых полимеров, полученные в программе HYDRO
экспериментально.
Наимено
Параметры «палочка»
«сфера»
«крест»
эксперимент
А0×1010
s×1013
3.20
2.92
3.04
3.75
-1.0
-1.1
-1.1
-1.9
[η]
0.24
2.9
0.15
3.0
0.19
2.9
0.23
3.3
А0×1010
s×1013
3.19
2.96
2.99
3.30
-1.1
-1.3
-1.2
-1.9
[η]
0.27
2.5
0.14
2.9
0.18
2.7
0.2
2.9
А0×1010
s×1013
3.32
2.87
3.06
3.34
-0.96
-1.3
-1.4
-2.1
[η]
0.79
1.7
0.20
2.3
0.20
2.5
0.32
2.0
вание
объекта
F-547
n=6
D×107
F-536
n=8
D×107
F-551
n = 11
D×107
31
и
ВЫВОДЫ
В работе изучены гидродинамические, оптические и конформационные
свойства сверхразветвленных пиридилфениленовых полимеров с различными
фокальными группами в растворе. Проведено всестороннее исследование
характеристик полимеров в разбавленных растворах методами молекулярной
гидродинамики (вискозиметрия, денситометрия, скоростная седиментация),
оп ти к и ( д и н ам и ч е ско е расс ея ние свет а, д инам ич е ско е д войно е
лучепреломление в потоке). Получены основные гидродинамические
характеристики этих «макромолекул-наночастиц» такие, как величина
характеристической вязкости, удельный парциальный объем, коэффициенты
диффузии и седиментации, охарактеризованы молекулярные массы полимеров
по седиментационно-диффузионному анализу для каждой из систем.
Определена величина гидродинамического инварианта и эффективного
гидродинамического радиуса макромолекул для исследуемых объектов.
Полученные значения величин оптического коэффициента сдвига для
изучаемых полимеров оказались достаточно большими и характерны для
молекул, не обладающих сферической симметрией.
Выполнено моделирование и рассчитаны гидродинамические параметры
двумя различными методами.
Во-первых, проведены компьютерные расчеты для исследуемых систем с
использованием программы HyperChem. Рассчитаны величины коэффициентов
диффузии, характеристической вязкости и гидродинамического инварианта для
двух предельных случаев конформации полимерных структур: «палочка» и
«сфера». При сопоставлении экспериментальных данных и данных,
полученных компьютерным расчетом, выявлено удовлетворительное
соответствие. Сделанные расчеты также указывают на то, что форма
исследуемых систем не является сферической.
32
Во-вторых, гидродинамические параметры сверхразветвленных
пиридилфениленовых полимеров (коэффициент диффузии, константа
седиментации, характеристическая вязкость) получены с использованием
специализированной программы HYDRO. Сопоставление характеристик,
полученных по результатам моделирования и эксперимента, даёт основания
полагать, что конформация объекта F-536, в структуре которого отсутствуют
кислородные мостики, ближе к сферической. При дальнейшей интерпретации
данных для этой структуры стоит использовать модель сплюснутого или
вытянутого эллипсоида с относительно небольшой асимметрией формы.
Для всех изученных сверхразветвленных пиридилфениленовых
полимеров экспериментальные гидродинамиче ские и оптиче ские
характеристики хорошо согласуются с результатами, полученными из
компьютерных расчетов .
33
БЛАГОДАРНОСТИ
Работа выполнена при финансовой поддержке РФФИ (проект № 15-33-50978).
Коэффициенты седиментации были определены на оборудовании Центра
диагностики функциональных материалов для медицины, фармакологии и
наноэлектроники Научного парка СПбГУ.
34
СПИСОК ЛИТЕРАТУРЫ
[1] Berzelius, J.J. Annual Report of Advances in Physics and Chemistry. 1847 – 26.
[2] Малков Г. В. Синтез и исследование свойств сверхразветвленных
полимеров, полученных по реакциям полициклоприсоединения. // Диссертация
на соискание степени кандидата химических наук : 02.00.06, 02.00.04; [Место
защиты: Ин-т хим. физики им. Н.Н. Семенова РАН]. - Москва, 2008. – 190 с.
[3] Люлин С. В. Структура и динамика сверхразветвленных полимеров и
комплексов на их основене. Компьютерное моделирование. // Диссертации на
соискание ученой степени доктора физико-математических наук: 02.00.06;
Санкт-Петербург, 2010 - 315 с.
[4] Jiang G., Chen W., Xia W. Chromophore-Containing Polymers for Trace
Explosive Sensors // Designed Monomers & Polymers. 2008. V.11. P.105.
[5] Chong-Qi Shou , Chang-Li Zhou, Chun-Bin Zhao , Zhi-Liang Zhang , Guan-Bin
Li, Li-Ren Chen. Preparation and evaluation of non-bonded hyperbranched polymercoated columns for capillary electrophoresis. // Talanta 63 (2004) 887–891.
[6] Карцова и др. Сорбционные и хроматографические процессы. //СанктПетербург, 2013. Т. 13. Вып. 5 – 6 с.
[ 7 ] К. Inoue. Functional dendrimers, hyperbranched and star polymers // Prog.
Polym. Sci. 25 (2000) 453 ± 571 465.
[8] Yan, D.; Frey, H.; Gao, C., Eds.; Hyperbranched Polymers // Wiley, Weinheim,
2011.
[9] Yan, D.; Zhou, Y.; Hou, J. Science 2004, 303, 65–67.
[10] S. Peleshanko, V. V. Tsukruk. Assembling Hyperbranched Polymerics // Journal
of Polymer Physics: Part B: Polymer Physics 2012, 50, 83 – 100.
[11] Kim, Y. H.; Webster, O. W. Macromolecules 1992, 25, 5561–5572/
[12] Flory PJ. J Am Chem. Soc. 1952; 74:2718.
[13] Hawker CJ, Chu F. Macromolecules 1996; 29: 4370.
[14] Kim YH, Webster OW. Macromolecules 1992; 25: 5561.
[15] N.V. Kuchkina, M.S. Zinatullina, E.S. Sercova, P.S. Vlasov, A.S. Peregulov, Z.B.
Shifrina. Hyperbranched pyridylphenelene polymers based on the first-generation
dendrimer as multifunctional monomer // The Royal Society of Chemistry: RSC
Adv., 2015 – 6 c.
[16] Цветков В. Н. Жесткоцепные полимерные молекулы. - Л.:Наука, 1986. –
380 с.
35
[17] Bruce J. Berne, Robert Pecora Dynamic Light Scattering. // Courier Corporation,
1976 – 376 c.
[18] В. Н. Цветков, В. Е. Эскин, С. Я. Френкель, «Структура макромолекул в
растворах». М., Главная редакция физ.-мат. литературы, 1964, 720с.
[19] W. Burchard Solution properties of branched macromolecules // Advanced in
polymer science, 1999, V. 43, p. 113-134
[20] J. Garc´ıa de la Torre, S. Navarro, M.C. L´opez Mart´ınez, F.G. D´ıaz, J.J. L
´opez Cascales Biophys. J. 1994, 67(2), 530-531.
[22] Ю. В. Зефиров, П.М. Зоркий. Ван – дер – Ваальсовы радиусы атомов в
кристаллохимии (исторический очерк).
[23] J. Garc´ıa de la Torre. Dynamic electro-optic properties of macromolecules and
nanoparticles in solution: A review of computational and simulation methodologies //
Colloids and Surfaces B: Biointerfaces 56 (2007) 4 – 15
36
Отзывы:
Авторизуйтесь, чтобы оставить отзыв