Исследование структуры и дипольной подвижности водородосвязанных растворов методом временной диэлектрической спектроскопии тема автореферата и диссертации по физике, 01.04.03 ВАК РФ

На правах рукописи

ЛУНЁВ ИВАН ВЛАДИМИРОВИЧ

ИССЛЕДОВАНИЕ СТРУКТУРЫ И ДИПОЛЬНОЙ подвижности ВОДОРОДОСВЯЗАННЫХ РАСТВОРОВ МЕТОДОМ ВРЕМЕННОЙ ДИЭЛЕКТРИЧЕСКОЙ СПЕКТРОСКОПИИ

Специальность 01.04.03 - радиофизика

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

Работа выполнена на кафедре радиоэлектроники Казанского государственного университета.

Научный руководитель: кандидат физико-математических наук,

доцент Ю.А. Гусев;

Научный консультант: кандидат физико-математических наук,

с.н.с. В.И. Архипов;

Официальные оппоненты:

Ведущая организация:

доктор физико-математических наук,

профессор Н.К. Гайсин;

доктор физико-математических наук,

профессор Г.Ф. Новиков;

Казанский физико-технический институт

им. Е.К. Завойского КазНЦ РАН

Защита диссертации состоится 12 ноября 2007 г. в 14— часов на заседании диссертационного совета Д 212.081.18 при Казанском государственном университете по адресу: 420008, Казань, ул. Кремлевская 18, физический корпус.

С диссертацией можно ознакомиться в библиотеке Казанского государственного университета.

Автореферат разослан «11» октября 2007 г.

Ученый секретарь Диссертационного совета Д 212.081.18,

Казанский государственный университет

доктор физико-математических наук, профессор ----— A.B. Карпов

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

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

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

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

Некоторые из перечисленных ограничений снимаются при использовании методов временной диэлектрической спекгроскопии(ВДС), успешно развивающейся в последнее время. Метод В ДС позволяет получать непрерывный диэлектрический спектр исследуемого объекта. Существенным преимуществом метода ВДС является возможность получения релаксационных характеристик непосредственно во временной области К сожалению, метод ВДС обладает недостаточно широким частотным диапазоном Ю6-Ю10 Гц, что приводит к трудностям в определении диэлектрических параметров вещества при измерениях в области температур, ниже 0°С.

Цель работы Экспериментальное и теоретическое исследование концентрационной зависимости диэлектрической релаксации водного раствора изопропилового спирта в диапазоне концентраций от 0,1 до 0,95 мольных долей воды в растворе, в диапазоне температур от -30°С до 20°С, методом ВДС Сформулированная цель потребовала решения следующих задач Задачи исследования

Экспериментальное исследование диэлектрической релаксации водных растворов изопропилового спирта методом ВДС, в диапазоне температур от -30°С, до + 20°С

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

Расширение рабочего частотного диапазона ВДС спектрометра путем применения высокоимпедансного активного пробника

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

Методы исследования При решении поставленных задач использовался метод сосредоточенной емкости ВДС [8] При обработке результатов применялись основные положения математической статистики, численные методы аппроксимации и решения систем алгебраических уравнений [9] На защиту выносятся следующие положения

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

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

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

Достоверность результатов обеспечивается-

Выполнением калибровочных измерений диэлектрической проницаемости эталонных жидкостей с известными свойствами и совпадением полученных результатов с литературными данными [1] с погрешностью 3%-5%,

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

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

Благодаря спроектированному активному пробнику нижняя граница частотного диапазона временного диэлектрического спектрометра доведена до 10"1 Гц

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

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

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

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

Апробация работы. Основные результаты работы были доложены и обсуждались на следующих конференциях

6th Liquid Matter Conference of the European Physical Society, Utrecht, the Netherlands, 2005,

4th Conference of the International Dielectric Society & 9th International Conference on Dielectric & Related Phenomena IDS & DRP 2006,Poznan, Poland

Молодежная научно-практическая конференция, посвященная 10-летию филиала КГУ в г Зеленодольске, 2006,

XIII Всероссийская конференция «Структура и динамика молекулярных систем», Яльчик, 2007,

Публикации. По материалам диссертации опубликовано 10 работ Структура и объем работы Диссертация состоит из введения и пяти глав, работа изложена на 132 страницах, содержит 27 рисунков и 2 таблицы Список цитируемой литературы содержит 90 наименований

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

Содержание работы

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

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

Измерения комплексной диэлектрической проницаемости в широком диапазоне частот и температур показали, что диэлектрические

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

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

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

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

Вторая глава посвящена практическому применению методов временной диэлектрической спектроскопии Рассмотрены метод многократных отражений и метод сосредоточенной емкости

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

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

Рис.1. Блок схема метода ВДС [10].

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

Рис. 2. Хщхжтерная форма сигналов, регистрируемого в эксперимента ВДС экспериментах. У0(1) -импульс отраженный от открытого конца линии; КО) — импульс отраженный от ячейки с образцом.

Напряжение на измерительной ячейке есть У(1у=ЩугУ<^, где У^) и Щ) - падающий и отраженный сигналы соответственно (рис. 2.). Тогда ток текущий через ячейку

Д0= , 0)

где 20 характеристическое сопротивление передающей линии. Акшвное сопротивление ячейки с образцом по постоянном)' току есть

и ток провод имости

Таким образом, полный ток через ячейку с образцом

Заряд <2(1) конденсатора, заполненного образцом, связан с диэлектрической функцией отклика <рф и приложенным к конденсатору напряжением Уф соотношением

V о у

г

где С0 - емкость пустой ячейки. Поскольку (¿)= > и комплексная

о

оницаемость есть

СО

£ * (со) = Е^ + - ехр(- 1й)()с11, (6)

диэлектрическая проницаемость есть

1 Х(Г0(/)-2ф))

(7)

где Ь - оператор преобразования Лапласа.

В методе многократных отражений исследуемым образцом заполняют отрезок (длиной /) коаксиальной линии Сигнал претерпевает многократные отражения на

Связь между элементом матрицы рассеяния 8ц(со) и е (со) образца задается

1 -р2(со) е'1М

(В)

это уравнение не имеет решения в явном ввде Вычислить в (т) при помощи этого уравнения можно лишь методом итераций [9] Однако в раде случаев нахождение точного решения затруднено - чаще всего, по причинам экспериментального характера, например из-за низкого отношения сигнал/шум при исследовании образца с высокой проводимостью

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

Основным условием, накладываемым на метод, является условие сосредоточенности емкости С1«с/(е (со^у,,), где с1- наибольший линейный размер измерительной ячейки, ув= а>/2ж - верхняя граница интересующего нас частотного диапазона, е (со) -комплекснаядиотегарическая проницаемость исследуемого образца, с - скорость света в вакууме Метод сосредоточенной емкости позволяет измерять параметры диэлектриков в диапазоне 106-10№ Гц Этот диапазон можно перекрыть двумя-тремя измерениями, варьируя емкость измерительной ячейки. Из прочих преимуществ можно отметить малое количество вещества, требуемое для измерений, возможность интерпретации результатов непосредственно во временной области.

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

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

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

Одно ю основных требований к сигналу при вычислениях, для всех кривых значения амплитуд сигналов должны быть определены для одних и тех же моментов времени. Если сигналы не согласованы по временной оси, их можно сдвигать относительно друг друга на промежутки времени, кратные Ли В программе реализована возможность единичных временных сдвигов на А1 для любой кривой, которые производятся по команде пользователя, но такой подход не учитывает объективности человека и может быть использован только для кошроля погрешности согласования. Более правильным является подход, разработанный для общего случая и подходящий как д ля временного сдвига, так и для сдвига по уровню На графике пользователем вручную выбирается участок на котором кривые должны совпадать и по которому необходимо произвести их совмещение. Дальнейший алгоритм построен на основе метода наименьших квадратов и на минимизации разности между выбранными участками кривых В случае если нужно совмести, кривые по амплитуде, сдвинув их к какой-то одной, для каждой из них поочфедно минимизируется сумма квадратов разностей

1=к

соответствующих точек (лежащих одна над другой) обрабатываемой у, и неподвижной У, 1фивой (на выбранном участке п-кточек). Если рассмотреть уравнение

5„к=±(У,-(у,+Х)\2 (Ю)

где X - величина сдвига и найти его минимум можно показать, что такой

способ минимизации в результате даст 1

Х=-

п — к

и 1 п

п-к

(Н)

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

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

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

О 5 10 15 20

t, ns

0 S 10 15 20

t, ns

a) 6)

Pu:. 3. Коррекция экжршкнтальных кривых no уровню амплитуды и временному сдвигу а) кривые до коррекгрщ б) кривые после коррекции.

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

t

u](t)= jut(f - t')h(t')cù' =ul(t)*h(t), (12)

—со

или в частотной области:

¡У;(й>)=Е/,.(©)я(<у), (13)

где u,(t) - сигнал, проходящий через фильтр. Видно, что если применил» фильтр ко всем экспериментальным кривым, то в выражениях (3), (4), (7) ничего не изменится т.к. функция Н(со) сократится. Это позволяет применять линейные цифровые фильтры, не алияя на результат вычислений. И хотя на первый взгляд их применение ничего не дает, они позволяют сильно уменьшить погрешность в вычислениях и что самое главное позволяют более корректно решать уравнение (5), (6), решение которого очень чувствительно к любым помехам в сигнале и тем более к случайным флукгуациям точек кривых.

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

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

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

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

> £

! ]

& щж

жмШ ]

IV * Йу

ж

№

148.5 148 0 147.5 147 0 146.5 Э 146.0 145.5 145.0 144 5 144.0

0.70 0.75 0.»

0.85 0,90 0.95 1.00 1.05 1, ПЭ

У,

'Ч/

^ \\

I

ч/

0.70 0.75 0.8

0.85 0.90 0.95 1.00 1.05 1, ПЭ

а) б)

Рис. 4. Применение цифровых линейных фильтров, а) кривые до обработки , б) кривые после обработки

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

В основных вычислениях метода ВДС трудности в численной реализации представляет решение уравнения

I

(14)

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

фЮ

4000000 300000020000001000000 о

-1000000 -2000000 -3000000 -4000000

I

Рис. 5. Расходящееся решение интегрального уравнения (14).

I, пэ

Рис. б. Сходящееся решение уравнения (14).

При разработке программы были протестированы разные алгоритмы, позволяющие решить поставленную задачу. № них был выбран метод, предложенный В.А. Гончаровым [10]. Дня данного класса функций, а именно монотонных функций, начальные значения которых нулевые, а значения при больших Г стремятся к наклонной или горизонтальной асимптоте, данный метод является одним из наиболее простых и приводит к хорошим результатам. В частности он легко реализуется численно.

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

1. Первый из использованных алгоритмов является полностью автоматизированным и не требует от пользователя указания, каких либо начальных условий для запуска и не требует дополнительной информации о виде функции. Результатом его работы, как и всех остальных алгоритмов, является вычисление значения угла наклона касательной, проведенной к функции отклика в области где последняя выходит на асимптоту. Основная задача, которую необходимо решить это и есть поиск области асимптоты. В данном случае она ищется методом последующего «перебора асимптот». Для этого сначала вся функция отклика <р(г) методом наименьших квадратов аппроксимируется прямой линией и запоминается дисперсия ф(1) относительно прямой. Далее убирается первая точка на и снова проводится аппроксимация прямой, вычисляется дисперсия. Значение полученных дисперсий сравнивается и если последняя дисперсия меньше первой, то теперь запоминается последняя и процедура по удалению точки в начале повторяется. Если же вновь

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

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

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

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

I, пв

Рис. 7. Функции диэлектрического отклика без учета проводимости и с учетам проводимости.

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

О)

®2 £Й I 2 )

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

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

В настоящее время среди импульсных методов наиболее широкое распространение получил метод сосредоточенной емкости, как наиболее простой в реализации. Однако для изучения сложных систем, таких как жидкие кристаллы, нефтяные дисперсные системы, биологические системы, его рабочий частотный диапазон 106Гц-^ Ю!0Гц недостаточно широк.

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

Предлагаемое решение позволяет существенно расширить частотный диапазон ВДС спектрометра. Главной особенностью прибора является использование активного пробника с большим 10 МОм и более) входным сопротивлением, то позволяет свести к минимуму шунтирующее влияние коаксиального тракта на исследуемый образец. Прибор имеет весьма широкий рабочий диапазон 10"' Гц -н 10ю Гц, что позволяет изучать релаксационные процессы, обусловленные как коллективными, так и индивидуальными молекулярными движениями.

Рис. 8. Функциональная схема диэлектрического спектрометра с активным пробником. 1 -генератор скачка напряжения; 2 — смеситель; 3 - блок синхронизации и формирования аналогового сигнала; 4-аналого-цифровойпреобразователь; 5-компьютер.

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

Время, с

Рис. 9. Типичная форма регистрируемых кривых.

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

0.0 5.0x10'" 1.0x10* 1.5x10"' 2.0x10'

Время, с

Рис. 10. Функция отклика для глицерина, Т=28'С. Черная линия-экспериментальные точки, круглые точки-результат фитинга распределением КВВ.

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

' ' - х<г\Л

ф(0 = (*,-0

1-ехр

(16)

где

//

£ ^ -статическая диэлектрическая проницаемость, 6Х -диэлектрическая

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

рассчитаны параметры £s, , Т , (X . Как видно из таблицы 1, рассчитанные

параметры находятся в хорошем согласии с данными из справочной литературы [1].

е*, т, с а

Значения 39,91 3,99 11,1е-10 0,936

Погрешность определения ±0,009 ±0,09 ±4,6е-12 ±0,004

Табличные значения [1] 41,9 4,0 10,4е-10

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

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

Эксперимент проводился методом временной диэлектрической спектроскопии.

В ходе проведения эксперимента были получены функции диэлектрического отклика и рассчитаны диэлектрические спектры водных растворов изопропилового спирта в диапазоне концентраций от 0,1 до 0,95 мольных долей воды в раствора в диапазоне температур от -30°С до +20°С.

ФИКЦИЯ атлечфнчвошго от

Коиплекная диа,еюри1«с«ая проницаемости, Х = 0.5

Время, с

Рис.11. Экспериментально полученные функции диэлектрического отклика, концентрация X» = 0,5 мсшьных долей воды в растворе.

Частота, Гц

Рис. ^.Экспериментально полученные

спектры комплексной диэлектрической проницаемости, концентрация Ху, - 0,5 мольных с

Экспериментально полученные функции диэлектрического отклика и диэлектрические спектры апроксимировались уравнением КВВ и распределением Коула-Девидсона, соответственно. Дтя сравнения на рис.13, приведено время диэлектрической релаксации Т рассчтатое по функциям даэлекфического отклика и ча

1Е-9

0.2 о,ч 0,6 о,а 1,0 Концентрация, Хв

Рис. 13. Зависимость времени релаксации т от концентрации водного раствора изопропилового спирта при температуре Г = 20 С. Х„- мольная доля содержания воды Квадратные символы 1 я-времена релаксации, полуденные во временной области при помощи формулы (1). Круглые символы 2 о -времена релаксации, полученные в частотной области при помощи формулы Ксула-Дэвидсона

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

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

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

1Е-3-,

Т 1E-1Q:

1Е-11 - ■

0,0 0,2 0,4 0,6 0,8 1,0

Концентрация,X

Рис. 14. Зависимость времени диэлектрической релаксации г от концентрациях логарифмическом масштабе.

Концентрационная зависимость времени диэлектрической релаксации, представленная рис. 14 показывает, что в области концентраций 0.8 - 0.85 (в единицах мольной доли воды) происходит смена характера релаксационного процесса В области концентраций от 0 до 0.8 релаксация подобна той, что происходит в чистом изопропиловом спирте. Известно [11], что в одноатомных спиртах ассопцаты представляют собой линейные цепочки. Молекулы воды обладают на порядок большей подвижностью по сравнению с подвижностью молекул спирта. Поэтому увеличение концентрации воды, приводит к укорачиванию характерного времени жизни ассоциата В области концентраций от 0.8 до 1 характер релаксации становится таким же, как в чистой воде. Перестройка структуры ассоциагов происходит при концентрации, когда на одну молекулу спирта приходится около 9 молекул воды. Мы предполагаем, что это концентрация, при которой могут существовать водоподобные ассоциаты. Под водоподобными ассоциатами мы понимаем их пространственную протяжённость (возможно, близкую к тегграэдрической структуре), в отличие от цепочечных, спиртоподобных. Это согласуется с выводами работы [12], где на основе диэлектрических данных для чистой воды показано, что минимальная пространственно коррелированная область с характерным временем жизни около 10 пикосекунд включаегг в себя 9 молекул.

Концентрация, Х^

Рис.15. Концентрационная зависимость параметра а, при температурах Т = 2№С, Т = ¡(ТС, Т = (УС, Т=-1(ГС.

Параметр а меняется от 1 до 0.8, с минимумом в точке с концентрацией 0.83. Это означает, что разброс времён элементарных релаксационных процессов при критических концентрациях наиболее широкий.

В работах [2]-[4] отсутствует температурная зависимость диэлектрических параметров. Все эксперименты проводились при комнатной температуре. Однако, в работах [5]-[7] предполагается активационный характер переориентации диполей. В этом случае температурная зависимость времени диэлектрической релаксации описывается формулой Эйринга [13]

г(т)=1^ех{Щ о7)

где И - постоянная Планка, кв - постоянная Больцмана, К - газовая постоянная, ДО -свободная энергия активации. Частным случаем формулы Эйринга, когда можно пренебречь температурной зависимостью ~1/Г по сравнению с температурной зависимостью ~ щ^АСШ) является аррениусовская температурная зависимость времени релаксации

г(Т) = т0 ехр[~] (18)

где т0 - постоянная, которая по порядку величины — 10~13 с совпадает с периодом врашэтельных колебаний диполя.

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

Рис. 16. Зависимость времени диэлектрической релаксации водного раствора юопропилового спирта от обратной температуры в логарифмическом масштабе для концентрации 0.9.

Из рис. 16. видно, что в диапазоне примерно от -10°С до комнатной температуры действительно наблюдается температурная зависимость времени релаксации, близкая к аррениусовской. При температурах ниже -10°С эта зависимость становится нелинейной, что характерно для стеклообразующих жидкостей. Для описания стеклообразующих жидкостей используют эмпирии-ческую формулу Фогеля-Фалчера-Таммана [14]

<*>,«* (Т) = Аехр

( В ^

V

(19)

где Л(Т) = т'1 " частота пика потерь; А, В, Т0 - эмпирические параметры.

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

Температура, К Температура, К

Рис. 17. Температурная зависимость частоты пика потерь водного раствора изопропилового спирта для концентрации ХЛУ= 0.4, Х„,=0.7, Хш= 0.8, 0.95.

230 240 2SO 260 270 280 290 300

Температура, К

Ln(io„)

2D

240 2S0 Z6Q 270 280 200 300

Температура, К

На рис.17 представлены температурные зависимости частоты пика потерь изопропилового спирта, для концентрации Xw= 0.4, Xw= 0.7, Xw^ 0.8, Xw= 0.95.

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

Таблица 2. Температура структурного перехода водного раствора изопропилового стрта для

различных концентраций

Концентрация, Xw Температура структурного перехода, Т0 Погрешность определения т„,%

0,4 94,6 10,2

0,5 147,35 8,1

0,6 176,7 4,7

0,7 201,7 4

0,8 211,4 5,3

0,85 214,1 4,9

0,9 218,5 6,11

Основные результаты и выводы

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

2 Методом временной диэлектрической спеюроскопии получены функции диэлектрического отклика водных растворов изопропилового спирта в диапазоне концентраций от 01 до 0 95 мольных долей воды в растворе, в диапазоне температур от -30°С до +20°С. Обнаружена одна область релаксации, свидетельствующая о кооперативном характере дипальной поляризации раствора.

3 На базе разработанного активного пробника расширен частотный диапазон временного диэлектрического спектрометра 10"'-10|0Гц.

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

По материалам диссертации опубликованы следующие работы

1 Лунев ИВ Структурные переходы в дипольной ориешационной подвижности водородосвязанных растворов Диэлектрическая релаксация водного раствора изопропилового спирта / И.В Лунев, ВЖАрхипов, Гончаров ВА, МАКазымова II Известия вузов Поволжский регион -2007. -Выпуск.2,-С 121-130

2 Гончаров В А. Исследование диэлектрических характеристик конденсированных сред с использованием высокоимпедансного акшвного пробника / ВА Гончаров, Р В Картин, ИВ Лун®, В ИАрхипов // Известия вузов Поволжский регион — 2007 Выпуск.4, (в печати)

3 Лунев ИВ Влияние морфологии шин на релаксацию воды I ИВ. Лунев, Р Р Нишатуплин, АЮ Заввдонов, ЮА Гусев, ИР Манюров // Журнал технической физики -2001 -Т 71 -Выпуск. 11,-С. 127-129

4 Мамин РФ Фогостимулированные явления в релаксорах / РФ Мамин, С А Мигачев, МФ Садыков, ИВ Луне», В В Изотов, Ю А Гусев // Письма в ЖЭТФ -2003,-Т 11 -С 1232-1236

5 Ivan Lounev Analysis of dielectric relaxation data m water-saturatedsands and clays /Ivan Lounev, Raoul Nigmatullm, Andrey Zavidonov, Yu Gusev, Ibragim Manurov, Renat Muslimov // Journal of Non-Ciystalhne Solids — 2002 - V 305, -P 255-260

6 Лунев И В Диэлектрическая релаксация водного раствора изопропилового спирта /ИВ Лунев, В И Архипов, Д В Сараев, Ю А Гусев, MA Казымова И "Структура и динамика молекулярных систем", Сборник тезисов ВыпускХП -2005 -С 123

7 Лунев И В Диэлектрическая релаксация водного раствора изопропилового спирта / И.В Лунев, В И Архипов, Д В Сараев, Ю А Гусев, M А Казымова // «Структура и динамика молекулярных систем», Сборник статей Выпуск ХП -2005 -Т 2 -С 18-21

8 Lunev L Dielectnc relaxation water solution of isoamil alcohol / L Lunev, V Arkhipov,

V Gusev, M Vasilieva // 6th Liquid Matter Conference of the European Physical Society, Book of Abstract—2005 -P 85

9 Ivan Lounev Dielectric relaxation water solution of isopropyi alcohol / Ivan Lounev, Vladimir Aifchipov and Denis Saraev // 9th International Conference on Dielectric & Related Phenomena IDS & DRP, Abstract book.—2006 -P 231-232

10 Картин P В Оптимизация алгоритма обработай ситнолов ВДС спектрометра дан корректного расчета диэлектрических спектров / Р.В Картин, ИВ Лун®, В А Гончаров, Ю А Гусев // Молодежная научно-практическая конференция, посвященная 10-лешю филиала КГУ в г Зеленодольске, Сборник материалов конференцииДООб - С 21-24

Цитируемая литература

1 Ахадов Я.Ю Диэлектрические параметры чистых жидкостей/ Ахадов Я.Ю // Москва, МАИ, -1999 -С 854

2 Mashimo S The dielectnc relaxation of mixtures of water and primary alcohol / S Mashimo, S Kuwabara, S Yagihara, К Higasi // Journal of Chemical Physics, - 1989 -

V 90 -№6,- P 3292-3294

3 Mashimo S Structures of water and primary alcohol studied by microwave dielectnc analyses / S Mashimo, T Umehara, H. Redlin // Journal of Chemical Physics, - 1991 -

V 95-№9,-P 6257-6260

4 Bao J.Z Microwave dielectnc characterization of binary mixtures of water, methanol, and ethanol) J.Z Bao, M.L Swicord, С С. Davis II Journal of Chemical Physics, -1996 - V104, -№12.—P 4441-4450

5 Sato T Dynamical aspects of mixing schemes in ethanol-water mixtures in terms of excess partial molar activation free energy, enthalpy, and entropy of the dielectnc relaxation process / T Sato, A Chiba, R Nazaki // Journal of Chemical Physics , - 1999 -V110 - №5, - P 2508-2521

6 Sato T Hydrophobic hydration and molecular association in methanol-water mixtures studied by microwave dielectnc analysis / T. Sato, A Chiba, R Nozaki // Journal of Chemical Physics,-2000 —№6,—V112,—P 2924-2932

7 Sudo S Dielectnc relaxation tone and relaxation time distribution of alcohol-water mixtures / S Sudo, N Shmyashiki, Y Kitsuki, S Yagihara // Journal of Chemical Physics,-2002 - V106,-P 458-464

8 Goncharov VA Treatment of time-domain spectroscopy data in the lumped-capacitance method / V A Goncharov, Yu.D Feldman // Further considerations Chem. Phys Letters,-1980,~v71,-N3.-P 513-518

9 Корн Г Справочник по математике для научных работников и инженеров' Г Корн,Т Корн//М.«Наука»,-1978,-С381

10 Гончаров В А Флуюуационная и импульсная диэлекгромегрия и молекулярное движение в жидких кристаллах / Д иссертация на соискание ученой степени кандидата физико-математических наукДСФТИ им Е.К. Завойского КФ АН СССР, —1986

11 Крестов ГА Термодинамика ионных процессов в растворах / ГА Крестов // Ленинград,Химия -1984 -С.272

12. Arkhipov VI, Agmon N Relation between Macroscopic and Microscopic Dielectnc Relaxation Times m Water Dynamics. Israel Journal of Chemistiy, 2003, - V 43, pp 363-371

13 Hill N.E Dielectric properties and molecular behavior / NE H511, W.E Vaughan, AJH. Price, M. Davis//London, Remhold, — 1969 -P 69-71

14 Nigmatullin R.R New approach in fee description of dielectric relaxation phenomenon, correct deduction and interpretation of the Vogel-Fulcher-Tamman equation / RR. Nigmatullin, S L Osokin, GJ Smith//PhysicaC Condensed Matter,-2003 -V15 -P 3481-3516

Отпечатано с готового оригинал-макета в типографии Издательства Казанского государственного университета им В И Ульянова-Ленина Тираж 100 экз Заказ 39/10

420008, ул Профессора Нужина, 1/37 тел 231-53-59,292-65-60

Введение.

Глава I

Краткое изложение основных положений теории диэлектрической проницаемости.

1.1 Поляризация диэлектриков.

1.2 Статическая диэлектрическая проницаемость.

1.3. Временные и частотные свойства диэлектриков.

1.4 Характеристика современных представлений о молекулярном движении в водородосвязанных растворах.

Глава II

Методы временной диэлектрической спектроскопии.

2.1. Общие принципы методов ВДС.

2.2 Метод многократных отражений.

2.3. Метод шунтирующей сосредоточенной емкости.

Глава III

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

3.1 Особенности обработки экспериментальных данных.

3.2. Предварительная обработка экспериментальных данных.

3.3 Основные вычисления методов ВДС.

3.4 Реализация программного обеспечения.

Глава IV

Диэлектрический спектрометр на базе высокоимпедансного активного пробника.

4.1 Функциональная схема спектрометра, основные соотношения, предварительные результаты.

4.2Реализация вычислений.

ГЛАВА V

Структурные переходы в дипольной ориентационной подвижности водородосвязанных растворов.

5.1 Эксперимент.

5.2 Анализ погрешностей.

5.3 Анализ температурной зависимости времён диэлектрической релаксации.

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

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

В качестве объекта исследования нами выбраны изопропиловый спирт и его водные растворы. Этот выбор обусловлен следующими причинами. Во-первых, чистые жидкости могут быть использованы для тестирования создаваемых измерительных средств, поскольку диэлектрические характеристики их представлены во многих литературных источниках [1]. Систематического исследования этих веществ в высокочастотной области при изменении температуры в широком диапазоне не проводилось. В работах [2-7] были измерены диэлектрические параметры водных растворов некоторых первичных спиртов (пропиловый, этиловый, метиловый).

Предполагалось, что процесс диэлектрической релаксации имеет активационный характер, хотя измерения проводились лишь при комнатной температуре. Поэтому выводы работ [0-7] вызывают сомнение.

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

Некоторые из перечисленных ограничений снимаются при использовании методов временной диэлектрической спектроскопии(ВДС), успешно развивающейся в последнее время. Метод ВДС позволяет получать непрерывный диэлектрический спектр исследуемого объекта. Существенным преимуществом метода ВДС является возможность получения релаксационных характеристик непосредственно во временной области. К сожалению, метод ВДС обладает недостаточно широким частотным диапазоном 10б-Ю10 Гц, что приводит к трудностям в определении диэлектрических параметров вещества при измерениях в области температур, ниже 0°С.

Цель работы. Экспериментальное и теоретическое исследование концентрационной зависимости диэлектрической релаксации водного раствора изопропилового спирта в диапазоне концентраций от ОД до 0,95 мольных долей воды в растворе, в диапазоне температур от -30°С до 20°С, методом ВДС. Сформулированная цель потребовала решения следующих задач.

Задачи исследования.

Экспериментальное исследование диэлектрической релаксации водных растворов изопропилового спирта методом ВДС, в диапазоне температур от -30°С, до + 20°С.

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

Расширение рабочего частотного диапазона ВДС спектрометра путем применения высокоимпедансного активного пробника.

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

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

На защиту выносятся следующие положения:

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

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

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

Достоверность результатов обеспечивается:

Выполнением калибровочных измерений диэлектрической проницаемости эталонных жидкостей с известными свойствами и совпадением полученных результатов с литературными данными [1] с погрешностью 3%-5%;

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

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

Благодаря спроектированному активному пробнику нижняя граница частотного диапазона временного диэлектрического спектрометра доведена до 10"' Гц. .

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

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

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

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

Апробация работы. Основные результаты работы были доложены и обсуждались на следующих конференциях:

6th Liquid Matter Conference of the European Physical Society, Utrecht, the Netherlands, 2005;

4th Conference of the International Dielectric Society & 9th International Conference on Dielectric & Related Phenomena IDS & DRP 2006,Poznan, Poland.

Молодежная научно-практическая конференция, посвященная 10-летию филиала КГУ в г. Зеленодольске, 2006;

XIII Всероссийская конференция «Структура и динамика молекулярных систем», Яльчик, 2007;

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

Основные результаты и выводы

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

2. Методом временной диэлектрической спектроскопии получены функции диэлектрического отклика водных растворов изопропилового спирта в диапазоне концентраций от 0.1 до 0.95 мольных долей воды в растворе, в диапазоне температур от -30°С до +20°С. Обнаружена одна область релаксации, свидетельствующая о кооперативном характере дипольной поляризации раствора.

3. На базе разработанного активного пробника расширен частотный диапазон временного диэлектрического спектрометра Ю'МО10 Гц.

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

Заключение

Концентрационная зависимость времени диэлектрической релаксации, представленная рисунке 5.4 (кривая 1), показывает, что в области концентраций 0.8 - 0.85 (в единицах мольной доли воды) происходит смена характера релаксационного процесса. В области концентраций от 0 до 0.8 релаксация подобна той, что происходит в чистом изопропиловом спирте.

Известно [89], что в одноатомных спиртах ассоциаты представляют собой линейные цепочки. Молекулы воды обладают на порядок большей подвижностью по сравнению с подвижностью молекул спирта. Поэтому увеличение концентрации воды, приводит к укорачиванию характерного времени жизни ассоциата. В области концентраций от 0.8 до 1 характер релаксации становится таким же, как в чистой воде. Перестройка структуры ассоциатов происходит при концентрации, когда на одну молекулу спирта приходится около 9 молекул воды. Мы предполагаем, что это концентрация, 1 при которой могут существовать водоподобные ассоциаты. Под водоподобными ассоциатами мы понимаем их пространственную протяжённость (возможно, близкую к тетраэдрической структуре), в отличие от цепочечных, спиртоподобных. Это согласуется с выводами работы [90], где на основе диэлектрических данных для чистой воды показано, что минимальная пространственно коррелированная область с характерным временем жизни около 10 пикосекунд включает в себя 9 молекул.

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

Наши экспериментальные исследования показали, что результаты экспериментов проводимых во временной области (метод ВДС) нужно описывать временными функциями, минуя дискретное преобразование Фурье

1. Ахадов Я.Ю. Диэлектрические параметры чистых жидкостей/ Ахадов Я.Ю. // Москва, МАИ, 1999. -С. 854 .

2. Mashimo S. The dielectric relaxation of mixtures of water and primary alcohol / S. Mashimo, S. Kuwabara, S. Yagihara, K. Higasi // Journal of Chemical Physics, 1989. - V. 90. -№6, - P. 3292-3294.

3. Mashimo S. Structures of water and primary alcohol studied by microwave dielectric analyses / S. Mashimo, T. Umehara, H. Redlin // Journal of Chemical Physics, 1991. -V.95.- №9, -P.6257-6260.

4. Bao J.Z. Microwave dielectric characterization of binary mixtures of water, methanol, and ethanol / J.Z. Bao, M.L. Swicord, C.C. Davis // Journal of Chemical Physics, 1996. - V.104, - №12. - P. 4441-4450.

5. Sato T. Hydrophobic hydration and molecular association in methanol-water mixtures studied by microwave dielectric analysis / T. Sato, A. Chiba, R. Nozaki //. Journal of Chemical Physics, 2000. - №6, - V.l 12, 2924-2932.

6. Sudo S. Dielectric relaxation time and relaxation time distribution of alcohol-water mixtures / S. Sudo, N. Shinyashiki, Y. Kitsuki, S. Yagihara // Journal of Chemical Physics, 2002. - V.l06, P. 458-464.

7. Kremer F., Broadband Dielectric Spectroscopy / F. Kremer, A. Schonhals // Springer, 2002, - P. 729.

8. Resibois P. Classical kinetic theory of fluids /Р. Resibois, M. De Leener // N.Y.: Wiley,- 1997.-P. 422.

9. Мамин Р.Ф. Фотостимулированные явления в релаксорах / Р.Ф. Мамин, С.А. Мигачев, М.Ф. Садыков, И.В. Лунёв, В.В. Изотов, IO.A. Гусев // Письма в ЖЭТФ. -2003,-Т. 11.-С. 1232-1236.

10. Дебай П. Полярные молекулы / ГНТИ, М. Л., 1931, С. - 183.

11. Ансельм А. ЖЭТФ 12,274 (1942).

12. Михайлов Г.П. Современные теории дипольной поляризации молекулярных конденсированных систем / Г.П. Михайлов, Л.Л. Бурштейн // Успехи физических наук, т. LXXIV, вып. 1, -1961. С. 3-30.

13. Фрелих Г. Теория диэлектриков / Г. Фрелих, М: Иностранная литература, 1960. с. 252.

14. Vayghan W.E. Dielectric Relaxation/ Ann.Rev.Phys.Chem.- 1979.-V.30.i1. P. 103-124.

15. Hill N.E. The theory of Complex Permittivity// J.Phys.C. -1972.-V.5, No.4.-P.415-424.

16. Коффи В. Молекулярная диффузия и спектры / М. Ивенс, П.Григолини // Москва: Мир, 1987. - С. 370.

17. Дебай П. Теория электрических свойств молекул / П. Дебай, Г. Закк // Пер. с нем. М. Л., 1936,-С. 144.

18. A.Schallamach. Trans. Far. Soc., 42A (1946) 180-185.The dielectric relaxation of mixtures of dipolar liquids.

19. Alexandrov and Dzhian. Tech. Phys.U.S.S.R., 5 (1938) 836.

20. Hassion F.X. Dielectric properties of liquid ethanol and 2-propanol / F.X. HassionR.H.Cole// J.Chem.Phys., 23(10) 1955.-P. 1756-1762.

21. Denney D.J. Dielectric properties of methanol and methanol- 1-propanol solutions / DJ.Denney, R.H.Cole J. // Chem.Phys., 23(10) 1955 - P. 1767-1772.

22. McDuffie G.E Dielectric Properties of Glycerol-Water Mixtures / G.E.McDuffie, Jr., R.G.Quinn, T.A.Litovitz II J.Chem.Phys., 37(2) 1962. - P. 239-242.

23. J.B.Segar, H.E. Oberstar, Ind. Eng.Chem. 43 (Î951) 2117

24. Matsuoka S. Entropy, Free Volume, and Cooperative Relaxation / J. Res. Natl. Inst. Stand.Technol., 102 (1997) P. 213-228.

25. Bertolini D. Dielectric relaxation of water-alcohol and alcohol-alcohol solutions / D.Bertolini, M.Cassetari,G.Salvetti, E. Tombari, S.Veronesi // J.Non-CrystallineSolids. 131-133 (1991)-P.l 169-1173.

26. Bertolini D. Dielectric relaxation in binary solutions: Theory and experimental results / D.Bertolini, M.Cassetari, E. Tombari // J.Chem.Phys. 108(15) (1998)-P. 6416-6430.

27. Фельдман Ю.Д. Временная спектроскопия диэлектриков / Ю.Д. Фельдман, Ю.Ф. Зуев, В.М. Валитов // Приборы и Техника Эксперимента, №3,1979, с. 5-20.

28. Терлецкая JI. А. Методы временной диэлектрической спектроскопии для исследования динамических свойств диэлектриков / Л.А. Терлецкая, Ю.М. Кесслер, Ю.В. Подгорный // Известия Вузов, Химия и хим. технология. -1978. Т. XXI вып. 10. - С. 5-20.

29. M.J.C.van Gemert, Philips Res. Repts, 1973,28,530

30. A.T. Стар, Радиотехника и радиолокация, 1960, с. 85, «Советское радио».

31. Cole R.H., Mashimo S., Winsor P., Evaluation of dielectric behavior by time domain spectroscopy. 3. Precision difference method / J.Phys.Chem. 1980. - V. 84.-P. 786-793.

32. Chahine R., Bose Т.К. Comparative studies of various methods in time domain spectroscopy // J. Chem. Phys.- 1980. V. 72(2). - P. 808-815.

33. Gestblom B. On the method of reflection and short in dielectric time domain spectroscopy / Chem. Phys. Letters. -1980. v.74 № 2. - P. 333-336.

34. Goncharov V.A. The combined total-reflection method in dielectric time domain spectroscopy / Chem. Phys. Letters. -1985. v. 117 № 1. - P. 52-56.

35. Гончаров B.A. Методика экспериментального изучения диэлектрической релаксации путем измерений во временной области / Казанский физ.-тех. ин-т. Казань, 1983. С. 19. - Деп. в ВИНИТИ 01.07.83,-№3579.

36. Гастон М.А.Р. Справочник по волновым сопротивлениям фидерных линий СВЧ. -Мю.: Связь, 1976.-С. 152.

37. Goncharov V.A. Treatment of time-domain spectroscopy data in the lumped-capacitance method. Further considerations / V.A. Goncharov, Yu.D. Feldman // Chem. Phys. Letters, 1980, - v.71, -.N 3, - P. 513 - 518.

38. FeldmanY.D. Time domain treatment of TDS data for the lumped-capacitance method / Y.D. Feldman, V.A. Goncharov, Y.F. Zuev, V.M. Valitov // Chem. Phys. Letters, 1978, - v.58, -.N 2, - P. 304 - 306.

39. Hall D.G. Dielectric polarization of sulfuric acid solution / D.G. Hall, R.H. Cole // J. Chem. Phys.- 1981. V. 85. - P. 1065-1069.

40. Iskander M.F. A time domain technique for measurement of dielectric properties of biological substances / M.F. Iskander, S.S. Stuchly// IEEE Trans. Instr. Meas. -1972.-IM-21. -N4.-P. 425-429.

41. R.H. Cole, J. Phys. Chem., 1975, 79,1459.

42. Хиппель A. P., Диэлектрики и волны, пер. с англ., М., -1960. С. 440.

43. Брандт А.А. Исследования диэлектриков на сверхвысоких, частотах. М. (1963). 403 с.

44. Валитов В.М. Автоматический временной диэлектрический спектрометр / В.М. Валитов, И.В. Ермолина, Ю.Ф. Зуев, Ю.Д. Фельдман // Журнал физической химии., 1987, -С. 564-569.

45. Фельдман Ю.Д. Разностная методика анализа диэлектрических данных во временной области / Ю.Д. Фельдман, Ю.Ф. Зуев, И.В. Ермолина, В.А. Гончаров // Журнал физической химии., 1988, т. 62, -С. 557-561.

46. Бойко Б. П. Основы радиоэлектроники. Часть 1. Сигналы. Учебное пособие для студентов специальности Радиофизика и электроника. / Физический факультет Казанского государственного университета, 2001.

47. Баскаков С.И. Радиотехнические цепи и сигналы. / М.: высшая школа, 1988,448 с.

48. Румянцев К.Е. Радиотехнические цепи и сигналы: учебник для сред.проф. образования / К.Е. Румянцев, П.А. Землянухин, А.И. Окорочков // М.:

49. Издательский центр «Академия», 2005 384 с.

50. С.Г. Михлин, X.JI. Смолицкий. Приближенные методы решения дифференциальных и интегральных уравнений. / Издательство «Наука», главная редакция физико математической литературы, Москва, 1965.

51. А.Ф. Верлань, B.C. Сизиков. Интегральные уравнения: методы, алгоритмы, программы. Справочное пособие. / Киев, Наукова Думка, 1986.

52. Богнер Р. Введение в цифровую фильтрацию / Р. Богнер, А.1

53. Константинидис // Мир. 1976.

54. Рабинер JL, Гоулд Б. Теория и применение цифровой обработки сигналов. М.: Мир, 1978. - 846 с.

55. Д. Худсон. "Статистика для физиков". Мир, 1970.

56. Банди Б. Методы оптимизации, вводный курс. / Москва, «Радио и связь», 1988 г.

57. Тихонов А.Н. Методы решения некорректных задач / А.Н. Тихонов, В.Я. Арсенин //№: Наука, -1979, -С. 288.

58. Гончаров В.А. Флуктуационная и импульсная диэлектрометрия и молекулярное движение в жидких кристаллах / Диссертация на соискание ученой степени кандидата физико-математических наук, КФТИ им. Е.К. Завойского КФ АН СССР, 1986.

59. Верлань А.Ф. Методы решения интегральных уравнений с программами для ЭВМ / А.Ф. Верлань, B.C. Сизиков // Справочное пособие, 1978,-С. 292.

60. Калиткин Н.Н. Численные методы / М.: Наука, 1978. -С. 512.

61. Васильев А.А. Методы решения задач линейного программирования / А.А. Васильев, В.А. Никитенков, Т.М. Никитенкова // Сыктывкар.: Изд. СГУ. -1990.-С. 73.

62. Dev S.B. Computational techniques in analysis of dielectric relaxation measurements // Adv. Relax. Processes. -1972.-V.4-N 2. -P. 159-191.

63. Gestblom B. About Fourier transform of TDS data / B. Gestblom, E. Noreland // J. Phys. Chem. -1976. V.80.-P. 1631-1634.

64. Balcou Y. A method increase the accuracy of fast Fourier transform calculations for rapidly varying functions // IEEE Transactions Instr. Meas-1981. -V.IM-30- P. 38-40.

65. А.Я. Архангельский. Программирование в Delphi 7 / M.: ООО «Бином-Пресс», 2005 г.-1152 с.

66. А .Я. Архангельский. Delphi 7. Справочное пособие / М.: ООО «Бином11. Пресс», 2004 г. 1024 с.

67. Фаронов В.В. Система программирования Delphi / СПб.: БХВ-Петербург, 2004.

68. Feldman Y.D. Treatment of TDS data for the lumped-capacitance method. Frequency domain / Y.D. Feldman, V.A. Goncharov, Y.F. Zuev, V.M. Valitov // Chem. Phys. Letters-1979-V.65.-N1.-P.69-70.

69. Goncharov V.A. On the study of the dynamical behavior of dielectrics by thermal-noise correlation analysis // Chem. Phys. Letters.-1984.-V.lll.-N6-P.521-525.

70. M. Шур. Физика полупроводниковых приборов / М., Мир, 1992.

71. Щедрин А.И. Схемотехника аналоговых и аналого-цифровых электронных устройств / 2005. -С. 530.

72. Гончаров В.А. Исследование диэлектрических характеристик конденсированных сред с использованием высокоимпедансного активного пробника / В.А. Гончаров, Р.В. Каргин, И.В.Лунёв, В.ИЛрхипов // Известия вузов. Поволжский регион. 2007. Выпуск.4, (в печати).

73. Эйзенберг Д. Структура и свойства воды / Д. Эйзенберг, В.Кауцман // Л,. Гидрометеоиздат, 1975.-С. 181.

74. Лунёв И.В. Влияние морфологии глин на релаксацию воды / И.В. Лунёв, Р.Р.i

75. Нигматуллин, А.Ю. Завидонов, Ю.А Гусев, И.Р. Майоров // Журнал технической физики.-2001.-Т. 71.-Выпуск. 11,-С. 127-129.

76. Ivan Lounev Analysis of dielectric relaxation data in water-saturatedsands and clays /Ivan Lounev, Raoul Nigmatullin, Andrey Zavidonov, Yu. Gusev, Ibragim Manurov, Renat Muslimov // Journal of Non-Crystalline Solids. 2002. - V. 305, - P. 255-260.

77. И.В. .Лунёв Диэлектрические свойства врды в катион-замещенном монтмориллоните / И.В. Лунёв, P.P. Нигматуллин, Ю.А. Гусев, Д. В. Сараев // «Структура и динамика молекулярных систем», Сборник статей. Выпуск IX,-2002.-Т. 1.-С. 317-320.

78. Nigmatullin R.R. Theory of Dielectric Relaxation in Non-Crystalline Solids: From a Set of Micromotions to the Averaged Colective Motion in the Mesoscale Region. Physica В: Physics of Condensed Matter, V.358, 2005, - P. 201-215.

79. Лунёв И.В. Диэлектрическая релаксация водного раствора изопропилового спирта / И.В. Лунёв, В.И. Архипов, Д.В. Сараев, Ю.А. Гусев, М.А. Казымова // "Структура и динамика молекулярных систем", Сборник тезисов. Выпуск XII. -2005.-С. 123.

80. Лунёв И.В. Диэлектрическая релаксация водного раствора изопропилового спирта / И.В. Лунёв, В.И. Архипов, Д.В. Сараев, Ю.А. Гусев, М.А. Казымова //

81. Структура и динамика молекулярных систем», Сборник статей. Выпуск XII. -2005.-Т. 2.-С. 18-21.

82. Lunev I. Dielectric relaxation water solution of isoamil alcohol /1. Lunev, V.th

83. Arkhipov, Y. Gusev, M. Vasilieva // 6 Liquid Matter Conference of the European Physical Society, Book of Abstract. 2005. - P. 85.

84. Ivan Lounev Dielectric relaxation water solution of isopropyl alcohol / Ivan Lounev, Vladimir Arkhipov and Denis Saraev // 9th International Conference on Dielectric & Related Phenomena IDS & DRP, Abstract book. 2006. - P. 231 -232.

85. Романычев Г.Д. Анализ источников погрешностей метода временной диэлектрической спектроскопии / Г.Д. Романычев, И.В. Ермолина, Е.А. Полыгалов, Ю.Ф. Зуев, Д.В. Упшинский, Ю.Д. Фельдман // Измерительная Техника, 1992, №8, с. 61-63.

86. Yu. Feldman. Time-domain dielectric spectroscopy: An advanced measuring system / Yu. Feldman, A. Andrianov, E. Polygalov, I. Ermolina, G. Romanychev, Yu. Zuev, B.Milgotin // Rev. Sci. Instrum. 67 (9), September, 1996.

87. Garg K., Smyth C.P., Journal of Chemical Physics, 1965, V.69, P. 1294.

88. Bendler J.T A New Volgel-Like Law: Ionic Conductivity, Dielectric Relaxation, and Viscosity near the Glass Transition / J.T. Bendler, J.J. Fontanella, M.F. Shlesinger // Physical Review Letters, 2001, V.87, №5, P. 1951 - 1954.

89. Крестов Г.А. Термодинамика ионных процессов в растворах / Г.А. Крестов // Ленинград, Химия. -1984. С. 272.

90. Arkhipov V.I., Agmon N. Relation between Macroscopic and Microscopic Dielectric Relaxation Times in Water Dynamics. Israel Journal of Chemistry, 2003, V.43, pp.363-371.