Гидродинамические характеристики стандартов полистиролсульфоната в растворах различной ионной силы

Гидродинамические характеристики стандартов полистиролсульфоната в растворах различной ионной силы

Санкт-Петербургский государственный университет

Физический факультет

Кафедра физики полимеров

Гидродинамические характеристики стандартов полистиролсульфоната в растворах различной ионной силы

Выпускная квалификационная работа на соискание академической степени магистра физики

Выполнена студентом Смирновым Артемом Михайловичем

Научный руководитель: доктор физ.-мат. наук Павлов Г.М

Рецензент: Мельников Александр Борисович

Санкт-Петербург

год.

Введение

Молекулярные свойства полимеров определяют физико-химические характеристики материалов, изготовленных на их основе. Конформация и размеры молекул ДНК, РНК и белков оказывают прямое влияние на способность этих молекул выполнять свои биологические функции. В связи с этим, изучение молекулярных свойств полимеров имеет принципиальное значение для решения прикладных и фундаментальных задач физикохимии высокомолекулярных соединений, биофизики, химии, фармакологии и др.

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

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

1. Основы статистики линейных полимерных цепей

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

.1 Идеальная полимерная цепь

Модель идеальной макромолекулы играет в физике полимеров такую же роль, как представление об идеальном газе в обычной молекулярной физике. Эта модель представляет собой цепочку из бестелесных звеньев; каждое из звеньев идеальной цепи соединено с двумя ближайшими по цепи соседями, но не взаимодействует ни с молекулами растворителя, ни с другими звеньями этой же или других макромолекул. Так же как существует много идеальных газов (одноатомный, двухатомный и т.д. - важно лишь, чтобы молекулы не взаимодействовали друг с другом), так есть и целый ряд моделей идеальных цепей; они различаются структурой звеньев и устройством связей между ближайшими соседями, но «идеальность» во всех случаях состоит в отсутствии объемных взаимодействий. Круг реальных условий, при которых молекулы ведут себя как реальные, не очень широк - в основной это разбавленные растворы полимеров в так называемых θ-растворителях, а также полимерные расплавы. Тем не менее, идеальные модели очень полезны, так как позволяют составить правильное представление о характере теплового движения макромолекул, другими словами - об энтропийных свойствах полимерного вещества.

.2 Свободно-сочлененная цепь

Для описания конформационных свойств полимерных молекул важное значение имеет модель свободно-сочлененной цепи, введенная в рассмотрение и разработанная Куном и Марком. В этой модели реальная полимерная цепь заменяется эквивалентной, состоящей из N прямолинейных сегментов длиной A, пространственные ориентации которых взаимно независимы. Полная длина эквивалентной цепи L принимается равной длине полностью вытянутой (без деформации валентных углов) реальной цепи (контурная длина):

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

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

Развитие статистики полимерных цепей на основе свободно-сочлененной модели привело к установлению важной закономерности - гауссова распределения расстояний h в ансамбле длинных цепных молекул (каждая длиной L):

Здесь  - вероятность того, что для произвольно выбранной из ансамбля цепи (состоящей из N сегментов длиной A) расстояние между ее концами лежит в пределах от h до dh. Цепи, удовлетворяющие этому распределению, называются гауссовыми. Из этой формулы может быть получен любой из моментов распределения / В частности, для среднего квадрата h получается величина

Отсюда следует, что средний квадрат расстояния между концами гауссовой цепи равен

т.е. пропорционален ее контурной длине L. Это свойство является основным для гауссовых цепей. Из последнего равенства следует, что линейные статистические размеры гауссовой цепи (характеризуемые величиной ) при заданной контурной длине L пропорциональны корню квадратному из длины сегмента Куна A. Поэтому эта длина статистического сегмента может служить мерой равновесной жесткости полимерной цепи.

.3 Гибкость полимерной цепи

Любая полимерная макромолекула обладает гибкостью, но механизм гибкости у разных полимеров разный.

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

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

.4 Персистентная (червеобразная) модель Порода

Для описания конформационных свойств жесткоцепных молекул наиболее подходящей является модель персистентной, или червеобразной цепи Кратки-Порода, в которой в отличие от свободно-сочлененной цепи Куна учитывается ориентационное близкодействие элементов, составляющих цепь. В основе модели Порода (так же как и в модели Куна) лежит цепь длиной L, состоящая из n прямолинейных сегментов длиной ΔL, так что . Однако в отличие от свободно-сочлененной цепи пространственные ориентации соседних элементов здесь не вполне взаимно независимы - направление первого элемента в определенной мере передается по цепи. Корреляция между элементами выражается в том, что среднее (по всем конформациям) значение  угла между соседними элементами не равно нулю и одинаково для всех элементов цепи (в случае свободно-сочлененной цепи k=0). Таким образом, совокупность величин ΔL и k является мерой ориентационного близкодействия, т.е. корреляции для персистенции в цепи. В итоге, угол  между элементами 1 и n определяется соотношением

здесь a - длина персистенции.

Если, оставляя постоянными значения L и a, перейти к пределу  (тогда ), то изломанная персистентная цепь превращается в непрерывную червеобразную цепь, которая определяется соотношением

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

.5 Размер идеальной полимерной цепи

Простейшая характеристика пространственного размера полимерной цепи - среднеквадратичной расстояние между ее концами; сравнение размера с контурной длиной характеризует степень запутанности цепи в пространстве.

Гибкость полимерной цепи ведет к тому, что макромолекулы никогда не имеют прямолинейной формы - напротив, любая достаточно длинная цепь, извиваясь, в любой момент имеет форму случайного запутанного клубка. Размер клубка в пространстве не характеризуется контурной длиной цепи, различие этих величин определяется тем, насколько запутана цепь. Именно поэтому необходимо обсуждать вопрос о размерах клубкообразных конформаций цепей. Для характеристики степени гибкости макромолекулы можно наряду с персистентной длиной использовать величину сегмента Куна.

.6 Структура и жесткость полимерной цепи

Основным механизмом, определяющим жесткость типичных гибкоцепных полимеров, является взаимодействие боковых групп, которое и учитывается в конформационной статистике полимерных цепей. Однако этот механизм не приводит к столь большим значениям , которые наблюдаются у жесткоцепных полимеров. Систематические данные по этому вопросу были получены при исследовании конформационных свойств гребнеобразных молекул - полиалкилакрилатов полиалкилметакрилатов с различной длиной боковых цепей. С удлинением боковой алкильной группы жесткость основной цепи несколько возрастает и для высших гомологов в 2-3 раза превосходит жесткость полиметилакрилата. Однако взаимодействие боковых цепей приводит к более значительному увеличению их равновесной жесткости и ориентационной упорядоченности по сравнению с основной цепью. Этот эффект значительно резче проявляется в молекулах, боковые группы которых, имея цепное строение, кроме того, имеют мезогенную структуру.

.7 Влияние эффектов исключенного объема

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

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

Однако, поскольку всякая реальная цепная молекула имеет бόльшую или меньшую гибкость, при ее тепловом движении всегда возможны случайные сближения атомов и групп, значительно удаленных друг от друга по цепи. При таких сближениях неизбежно возникновение взаимодействия между сблизившимися элементами цепи, имеющего характер их взаимного отталкивания, тем большего, чем бόльший эффективный объем занимает взаимодействующая пара элементов («исключенный объем»). Эти взаимодействия, являющиеся взаимодействиями дальнего порядка, принято называть эффектами исключенного объема, поскольку в их основе лежит невозможность для двух элементов цепи одновременно занимать в пространстве один и тот же элемент объема. Эффекты исключенного объема возмущают конформацию клубкообразной молекулы, приводя (в силу возникающих отталкиваний) к увеличению средних расстояний между ее элементами, в том числе и к увеличению  и . Количественно эти возмущения характеризуют коэффициентами αh и αR линейного увеличения размеров молекулярного клубка, определяемыми соотношениями

Здесь  и  - средний квадрат расстояния между концами цепи и ее радиус инерции в отсутствие объемных эффектов;  и  - те же величины, возмущенные объемными эффектами.

Конформации реальных полимерных молекул изучаются в разбавленных растворах, где объемные эффекты существенно зависят от взаимодействий молекул полимера с молекулами растворителя и для одного и того же полимера могут быть весьма различны в различных растворителях. Подбором достаточно «плохого» растворителя и соответствующей температуры («θ-температуры») влияние конечного объема мономерной единицы можно скомпенсировать взаимным притяжением единиц цепи. В этих условиях эффекты исключенного объема отсутствуют, в равенствах () коэффициенты . С улучшением термодинамического качества растворителя и соответствующим усилением взаимодействий полимер-растворитель притяжение между элементами цепи не в состоянии компенсировать их отталкивания, эффект исключенного объема увеличивается и α становится больше единицы. Причем возмущенные размеры молекул  и  растут быстрее, чем пропорционально длине цепи L. В итоге получаем в грубом приближении

где ε>0.

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

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

2. Тория вращательного трения для моделей персистентного ожерелья и червеобразного цилиндра

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

2.1 Ожерелье

Вращательное трение червеобразной цепи было рассмотрено Хирстом, использовавшим и развившим формализм, разработанный Кирквудом, в его теории вращательного трения палочкообразного ожерелья.

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

Эти выражения в области малых  переходят в соотношения

Откуда, при  следует  и , что соответствует конформации прямой тонкой палочки. В области  следует , что соответствует гауссову клубку с распределением сегментов, в средней сферически-симметричным относительно средней точки цепи.

Окончательные выражения для коэффициента вращательного трения  червеобразной цепи при ее вращении вокруг оси  или  Хирст получает для двух предельных случаев.

Для длинных цепей, когда  (червеобразный клубок),

Если выразить  через диаметр  эквивалентной стоксовой сферы  и принять модель червеобразного ожерелья с соприкасающимися бусами , то последние выражения трансформируются в выражения

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

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

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

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

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

.2 Червеобразный цилиндр

Под червеобразным цилиндром понимают цилиндр, изогнутый таким образом, что форма его осевой линии описывается уравнением червеобразной цепи

и соответственно расстояние между любыми двумя его точками на осевой линии определяется по формуле

Гидродинамическое сопротивление, испытываемое таким телом при его движении в вязкой жидкости, вычисляется методом Озеена - Бюргерса.

Теория характеристической вязкости раствора жесткоцепных молекул на основе модели червеобразного цилиндра была разработана Ямакавой, применившим для этой цели метод Озеена - Бюргерса.

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

Пропуская все вычисления, можно сказать, что для червеобразного сфероцилиндра (при ) с хорошим приближением может быть использовано выражение

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

3. Гидродинамика растворов полимеров, характеристическая вязкость, взаимодействия в цепи

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

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

3.1 Поступательное трение

Количественной характеристикой трения при поступательном движении тела в окружающей его жидкости является коэффициент поступательного трения , определяемый выражением

где  - скорость движения теля, вызванного действием силы . Направление  совпадает с направлением

Строгое решение уравнений гидродинамики Навье-Стокса для простейшего случая движения шара в вязкой жидкости приводит к формуле

где  - коэффициент вязкости окружающей жидкости;  - диаметр шара.

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

.2 Вращательное трение

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

Коэффициент пропорциональности  - коэффициент вращательного трения шара. Он пропорционален объему шара  и вязкости растворителя , согласно равенству

где  - диаметр шара.

Для вытянутых эллипсоидов вращения ()

3.3 Модель червеобразного цилиндра

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

и соответственно расстояние между любыми двумя точками на осевой линии определяется по формуле

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

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

.4 Характеристическая вязкость

Теория характеристической вязкости  раствора жесткоцепных молекул на основе модели червеобразного цилиндра была разработана Ямакавой с сотрудниками, применившими для этой цели метод Озеена-Бюргерса.

Зависимость характеристической вязкости от молекулярных параметров Μ, L, d и A для червеобразной цепи можно представить в виде уравнения Флори

,

Где Φ является функцией этих параметров.

.5 Эффект исключенного объема

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

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

Поступательное трение жесткоцепного полимера, моделируемого червеобразным ожерельем с учетом эффектов исключенного объема, вычислялось в работе Хирста [Gray H.B., Bloomfield V.A., Hearst J.E., J. Chem. Phys., 1967, vol. 46, N 4, p. 1493-1498]. Вычисления проводились с использованием параметра ε, характеризующего отклонение статистических размеров цепи от гауссовых свойств за счет объемных эффектов. Второй момент распределения  по расстояниям  в гауссовых цепях, возмущенных объемными эффектами, связан с невозмущенной величиной  (в θ-условиях) соотношением

Для коэффициентов поступательного трения в области больших , в итоге, было получено выражение

где функция  табулирована в работе [Gray H.B., Bloomfield V.A., Hearst J.E., J. Chem. Phys., 1967, vol. 46, N 4, p. 1493-1498].

Таким образом, при наличии эффектов исключенного объема в области больших  соотношение  линейно относительно . Так же было показано, что в этих условиях  линейно относительно .

4. Полиэлектролиты

.1 Полианионы, поликатионы, полиамфолиты

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

Классические представители полиэлектролитов - полиакриловая и полиметакриловая кислоты.

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

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

К рассматриваемому классу поликислот относятся также многие полимеры биологического происхождения. Здесь стоит назвать в первую очередь нуклеиновые кислоты - ДНК и РНК, передающие генетическую информацию.

Другой класс полиэлектролитов - полиоснования.

Конфигурационные свойства полиоснований аналогичны свойствам поликислот. В частности, полиоснования, как и поликислоты, сильнее ионизованы в солевой форме.

Сочетание кислотных и основных групп в одной цепи приводит к образованию полиамфолитов, составляющих третий класс полиэлектролитов. Для каждого полиамфолита существует определенное, зависящее от его состава, значение pH, при котором количества положительных и отрицательных зарядов в цепи равны. Иными словами, суммарный заряд полиамфолита в этой изоэлектрической точке (ИЭТ) равен нулю. При pH ниже ИЭТ в цепи начинают доминировать положительные заряды. При достаточно низком pH ионизация всех кислотных групп оказывается подавленной и полиамфолит превращается в полиоснование. Наоборот, по мере повышения pH над ИЭТ полиамфолит постепенно превращается в поликислоту. Полиамфолиты играют огромную роль в природе: все белки относятся к полиамфолитам.

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

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

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

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

Флори построил теорию на основе объемных эффектов. Электростатическое отталкивание приводит к набуханию клубка, зависящему от ионной силы. Флори предполагал, что клубок вместе с окружающим его растворителем электрически нейтрален. Расчет показывает, что электростатические взаимодействия не могут превратить клубок в вытянутую молекулу, происходит лишь набухание.

Птицын развил более строгую теорию. Основное предположение 0 из-за большего экранирующего действия контрионов заряженные группы макромолекулы, расположенные далеко друг от друга по цепи, взаимодействуют лишь при случайном сближении в результате флуктуационного изгибания цепи. Также, из этой теории следует, что конформационные свойства заряженных макромолекул промежуточны между свойствами ненабухших клубков и жестких стержней:

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

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

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

.2 Полиэлектролитное набухание

Далее следует рассмотреть исследованное как экспериментально, так и теоретически так называемое полиэлектролитное набухание. Причина полиэлектролитного набухания - электростатическое отталкивание одноименно заряженных звеньев цепи, приводящее к развертыванию клубков и увеличению их линейных размеров. Поскольку о размерах макромолекул можно судить по характеристической вязкости [η], пропорциональной объему клубков, первые оценки полиэлектролитного набухания были произведены по измерениям [η] в зависимости от степени ионизации. Было показано, в частности, что при полной ионизации полиметакриловой кислоты [η] может возрасти на два порядка, чему соответствует увеличение линейных размеров клубков в 5-6 раз. Однако при этом сразу возникает вопрос, в какой степени полиэлектролитное набухание можно считать изотропным?

Одна из первых теорий эффекта полиэлектролитного набухания принадлежит А. Качальскому и Лифсону. Они полагали, что функция распределения расстояний между концами заряженной цепи имеет вид

где индекс «0» соответствует незаряженной цепи, а  - электростатическая энергия цепи, вычисляемая как сумма энергий отталкивания  всех пар заряженных звеньев. В первоначальном варианте теории, принадлежащем Качальскому, Кюнцле и В. Куну, экранирование электростатических взаимодействий, обусловленное образованием дебай-хюккелевской атмосферы противоионов вокруг заряженных групп цепи, не принималось во внимание. Основанием для подобного пренебрежения было следующее неправильное (как было установлено позже) допущение. Плотность атмосферы противоионов, характеризуемая параметром χ в теории Дебая-Хюккеля, определяется ионной силой раствора μ:

где

 - заряд электрона,  - валентность иона сорта ,  - число таких ионов в 1 см3, ε - диэлектрическая проницаемость раствора.

Если потенциальная энергия кулонова взаимодействия двух полностью изолированных ионов равна

то в среде, содержащей ионы противоположного знака (противоионы), благодаря образованию облака противоионов вокруг каждого иона энергия взаимодействия ослабевает и определяется выражением

Здесь  - сумма вандерваальсовых радиусов ионов. Параметр χ имеет геометрический смысл обратной величины эффективного радиуса ионной атмосферы, определяя расстояние от иона, за пределом которого осуществляется полное экранирование взаимодействий.

Качальский, Кюнцле и Кун предположили, что при очень большом разбавлении множитель Дебая-Хюккеля  можно принять равным единице, так как противоионы равномерно распределены в объеме раствора, тогда как заряды сосредоточены в малых дискретных областях, занятых макромолекулами. При этом экранирование действительно практически не должно иметь места, и при расчете  можно пользоваться кулоновым потенциалом. Полагая, что  - гауссова функция, они получили:

где  - число ионногенных групп в макромолекуле, т.е. в случае гомополимеров степень полимеризации, а  - степень ионизации, определяемая значением pH среды.

Это соотношение предсказывает анизотропное развертывание цепочек с переходом к практически полностью вытянутым конфигурациям. Предположение χ=0, сделанное при выводе этого уравнения, оказалось неверным потому, что на самом деле противоионы не распределены равномерно по объему раствора, а удерживаются электростатическим полем клубка. Кроме того, следует учитывать хотя и слабую, но все же конечную диссоциацию воды.

В теории полиэлектролитов рассматриваются две характеристические длины - радиус Бьеррума  и радиус экранирования Дебая-Хюккеля . Радиус Бьеррума  характеризует экранирующее действие растворителя. Радиус экранирования Дебая - Хюккеля  - расстояние, на котором распространяется действие электрического поля отельного заряда, помещенного в среду, содержащую другие заряды. Здесь  - элементарный заряд,  - электрическая постоянная,  - диэлектрическая проницаемость среды,  - ионная сила раствора,  - число -тых ионов в единице объема,  - заряд -того иона в единицах .

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

.3 Вискозиметрия полиэлектролитов

Характеристическая вязкость является одной из самых широко используемых молекулярно-гидродинамических характеристик полимеров. Это объясняется ее информативностью, а также сравнительной простотой и доступностью метода.

.4 Сущность явления

Характеристическая вязкость раствора высокомолекулярного вещества [η] имеет размерность удельного объема и служит мерой дополнительных потерь энергии, связанных с вращением макромолекул в потоке.

Внутреннее трение или вязкость всякой жидкости проявляется в тех случаях, когда она находится в состоянии потока с отличным от нуля градиентом скорости. Простейший пример такого потока - ламинарный поток с постоянным градиентом скорости , направление которого нормально направлению скорости. Скорость жидкости при этом определяется выражениями

Чем больше внутреннее трение в жидкости, тем большее напряжение сдвига τ нужно приложить, чтобы поддерживать поток с заданным градиентом скорости . Последнее выражается формулой Ньютона:

Коэффициент пропорциональности η называют коэффициентом вязкости или просто вязкостью жидкости. Жидкости, для которых η не зависит от τ, называют ньютоновыми.

Определение характеристической вязкости вытекает из выражения

где  - удельная вязкость раствора.

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

Более или менее строгие теории вязкости растворов цепных молекул (Кирквуда - Райзмана, Зимма) приводят при достаточно большом числе сегментов в макромолекуле к совпадающему результату - соотношению

со значением коэффициента 1/моль.

Для гауссовых клубков  и последнее уравнение переходит в  с =0.5. Такая зависимость характеристической вязкости  действительно наблюдается в идеальных растворителях.

Использование данных вискозиметрии для определения молекулярного веса и размеров макромолекул

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

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

и определение  требует независимого определения . Даже если эта величина известна, препятствием к непосредственному использованию этой формулы служит еще некоторая неопределенность параметра , который монотонно убывает при улучшении качества растворителя по формуле:

где  - термодинамический параметр, характеризующий качество растворителя.

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

( - вязкость растворителя) зависит от  только через , что же касается параметра , то он по сравнению с  слабо зависит от качества растворителя.

Характеристическая вязкость  связана с молекулярными характеристиками линейной цепной макромолекулы известным соотношением Флори

где  - среднеквадратичное расстояние между концами полимерной цепи,  - молекулярная масса исследуемого образца,  - гидродинамический параметр Флори.

Широкое применение в полимерной науке и технологии имеет соотношение Куна-Марка-Хаувинка-Сакурады:

где  и  - параметры, определенные и табулированные для многих систем полимер-растворитель. Известен ряд экстраполяционных процедур, позволяющих получать величину  незаряженных макромолекул []. Наиболее известным является уравнение Хаггинса:

В паре с уравнением Хаггинса рассматривают уравнение Крэмера

где  - безразмерный параметр Крэмера.

Сравнительно недавно Федорс [] предложил построение

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

В то же время в реологии известно, что динамическая вязкость раствора полимера в широком интервале концентраций может быть представлена как

Из этого соотношения следует, что  в каждой точке зависимости  от  может рассматриваться как текущее значение характеристической вязкости:

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

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

5. Определение характеристической вязкости линейных полиэлектролитов

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

.1 Диффузия макромолекул в растворах

.1.1 Основы Явления

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

Преимущество метода диффузии по сравнению с другими гидродинамическими методами (например, вискозиметрией) заключается в том, что в диффузионных измерениях исследуемые макромолекулы изучаются практически в их «естественном» состоянии, не возмущенном воздействиями внешних сил (механических, электрических и т.п.).

Макроскопическая направленная диффузия в растворе возникает при наличии градиента концентрации  растворенного вещества  в направлении . Она описывается двумя уравнениями Фика, из которых первое служит также определением коэффициента диффузии :

Здесь  - поток вещества, диффундирующего в направлении  (масса растворенного вещества, прошедшего в единицу времени через единицу площади сечения, нормального к x).

Второе уравнение Фика, характеризующее нестационарный процесс диффузии, выражает изменение концентрации раствора во времени  в результате такой диффузии. Для одномерного случая (диффузия в одном направлении) оно вытекает из предыдущего уравнения, если учесть, что . Тогда получим:

Если в рассматриваемой системе коэффициент диффузии  есть величина, не зависящая от , то тогда получаем

Если в диффузионном опыте распределение  гауссово, следует

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

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

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

где υ - скорость, с которой движется частица в растворителе под действием силы  вдоль оси .

Согласно Эйнштейну, при стационарном процессе диффузии, вызванном наличием градиента концентрации , «движущей силой» служит градиент  осмотического давления . При этом сила, действующая на частицы в 1 см3 раствора, равна - , а на одну частицу соответственно

Под действием силы  частицы движутся со скоростью υ и создают диффузионный поток вещества

Сопоставляя предыдущие уравнения, получаем:

где- осмотическое давление раствора.

В случае идеального раствора невзаимодействующих частиц осмотическое давление определяется законом Вант-Гоффа

В этом случае «движущая сила» диффузии равна

Отсюда получаем уравнение Эйнштейна

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

6. Схема диффузометра

.1 Седиментация

.1.1 Основные понятия

В отличие от диффузии, где поступательное движение молекул происходит под действием сил, вызванных градиентом осмотического давления, при седиментации - это центробежная сила , вызванная вращением ротора ультрацентрифуги:

Здесь  - масса молекулы;  - расстояние ее от оси вращения;  - угловая скорость ротора;  - множитель плавучести, зависящий от плотности растворителя  и парциального удельного объема  полимера.

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

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

Эта формула показывает, что скорость седиментации  в точке  пропорциональна центростремительному ускорению  в этой точке, а коэффициент пропорциональности  называют коэффициентом седиментации. Единицей измерения в системе СГС является см/сек·дин. Практической мерой является единица Сведберга (S)= СГС.

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

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

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

Коэффициент седиментации макромолекул определяется экспериментально при скоростном центрифугировании путем фиксирования положения (координата ) седиментационной границы в соответствующие моменты времени . При этом используется выражение

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

Другими формами этого выражения также являются

где  - координата седиментационной границы в начальный момент времени ().

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

Зависимость коэффициента седиментации  от концентрации  исследуемого полимерного раствора обычно выражена значительно заметнее, чем в случае диффузии;  уменьшается с возрастанием . Поэтому для определения константы седиментации  необходима экстраполяция экспериментальных значений  к условиям предельно разбавленного раствора.

Последнее выполняется с использованием эмпирического линейного соотношения

Константа седиментации для монодисперсного полимера, пропорциональна отношению :

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

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

Зависимость константы седиментации от молекулярного веса.

В соответствии с гидродинамическими теориями коэффициента поступательного трения, зависимость  от молекулярного веса  может быть представлена формулами

Это эмпирическое соотношение Марка - Куна - Хаувинка,  - тот же параметр, что в формулах

. Экспериментальные данные и их обсуждение

Было изучено гидродинамическое поведение стандартов полистирол-4-сульфоната натрия

Рис. 1 Структурная формула полистирол-4-сульфоната натрия (ПС).

в чистой воде, 0.2 M и 4.17 М NaCl.

Таблица 1. Гидродинамические характеристики образцов полистиролсульфоната натрия в 0.2M NaCl, вязкость в 4,17 М NaCl и в воде при 250С

* - коэффициенты седиментации измерены при 200С

,

Исследование изотермическое поступательной диффузии проводилось в диффузометре В.Н. Цветкова. Результаты представлены на графике

рис.2. Зависимость дисперсии диффузионной границы s2 (см2) от времени процесса диффузии t (с) для полистиролсульфонатов натрия в 0.2M NaCl при 250С.

Экспериментальные точки с хорошей точностью ложатся на прямые.

Седиментация проводилась на ультрацентрифуге

Рис. 3. Концентрационная зависимость коэффициента седиментации для образца 4 в 0.2 М NaCl при 20°С

На основе диффузионно-седиментационного анализа были определены молекулярные массы (табл.1)

Исследования вискозиметрии проводились с применением капиллярного вискозиметра Освальда. Характеристическая вязкость рассчитывалась с помощью уравнений Хаггинса, Крэмера и Федорса. Сравнение результатов применения всех трех уравнений представлено на рис.

Сводный график результатов исследования вискозиметрии в 0,2 М NaCl и в 4,17 М NaCl по Хаггинсу на рис

Рис.. Построения Хаггинса hsp/с (см3×г-1)-с (г×см-3) для определения значений характеристической вязкости для полистиролсульфонатов натрия в 0.2M NaCl и в 4,17M NaCl при 250С.

Кроме того для определения характеристической вязкости было использовано уравнение Крэмера для двух образцов (№1 и №5) (рис. )

А также для бессолевых растворов полиэлектролитов

Рис. Зависимость lnhr - c для водных бессолевых растворов полистиролсульфонатов натрия при 250С.

На основе проведенных измерений была посчитана характеристическая вязкость исследованных образцов.

Также были построены зависимости Куна-Марка-Хаувинка-Сакурады, на основе чего получены скейлинговые соотношения, связывающие молекулярную массу с гидродинамическими характеристиками. Оказалось, что скейлинговый индекс, характеризующий зависимость характеристической вязкости от молекулярной массы при переходе от 0.2 M к 4.17М NaCl уменьшается. Таким образом, с увеличением ионной силы растворителя в макромолекуле происходят конформационные изменения, что приводит к уменьшению объема макромолекулярного клубка. При переходе от 4.17 М NaCl к чистой воде объем молекул увеличивается в 200-230 раз в зависимости от молекулярной массы.

Рис.6 - зависимости [h]~M в 0.2M NaCl ([h]~M0.9), в 4,17 M NaCl ([h]~M0.7) и в чистой воде([h]~M1.7).

В результате анализа зависимости между  и ,  и  в двойном логарифмическом масштабе (рис) получаем соотношения типа Куна-Марка-Хаувинка-Сакурады для молекул полистиролсульфоната в 0,2 М растворе NaCl при 25 .

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

дает среднее значение

которое использовали при дальнейшей интерпретации гидродинамических характеристик.

Рис.7. Зависимости Куна-Марка-Хаувинка-Сакурады для полистиролсульфонатов натрия в 0.2M растворе NaCl.

На основе теории Грея-Блюмфельда-Хирста, описывающей гидродинамическое поведение червеобразных цепей с учетом, как эффектов протекания, так и эффектов исключенного объема, проведены оценки длины статистического сегмента и гидродинамического поперечника молекул полистирол-4-сульфоната натрия.

Для определения длины статистического сегмента Куна и гидродинамического поперечника макромолекулы использовали соотношение

где  моль-1 - вязкостный параметр Флори,  - гидродинамический параметр Флори,  - масса единицы длины полимерной цепи,  - длина статистического сегмента,  - гидродинамический диаметр, .

Соотношение описывает зависимость коэффициента седиментации червеобразного ожерелья с учетом эффектов дальнодействия и внутримолекулярного протекания (Теория Грея-Блюмфельда-Хирста)[].

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

Рис. 8. Построение (M2/[h])1/3 - M(1-e)/2, используемое для оценки длины сегмента Куна A молекул полистиролсульфоната натрия

Далее оценивалась электростатическая составляющая, вносящая свой вклад в величину A

где  - радиус Дебая-Хюккеля

 - радиус Бьеррума

Таблица 2. Сравнение общей и электростатистической компоненты жесткости ПС для различных ионных сил.

Список использованной литературы

1.       Цветков В.Н. Жесткоцепные полимеры. М.: Наука, 1986.

.        Павлов Г.М., Зайцева И.И., Губарев А.А., Гаврилова И.И., Панарин, Е.Ф.// Журнал прикладной химии. 2006. т. 79. вып. 9, С. 1506-1509.

.        Павлов Г.М., Губарев А.А., Зайцева И.И., Сибилева М.А.// Журнал прикладной химии. 2006. т. 79. вып. 9, С. 1423-1428

.        Цветков В.Н., Эскин В.Е., Френкель С.Я. Структура макромолекул в растворах. М.: Наука, 1964.

5.       Gray G., Bloomfield V., Hearst J.// J. Chem. Phys. 1967. V. 46