Журнал физической химии, 2019, T. 93, № 8, стр. 1277-1280

Кластерная модель структуры 1,3-пропандиола и ее описание на основе диэлектрических данных

В. И. Журавлев a*

a Московский государственный университет имени М.В. Ломоносова, Химический факультет
Москва, Россия

* E-mail: zhura-061@yandex.ru

Поступила в редакцию 20.09.2018
После доработки 20.09.2018
Принята к публикации 09.10.2018

Полный текст (PDF)

Аннотация

Продолжено развитие подходов, связывающих значения статической (εs) и комплексной (ε*) диэлектрической проницаемости с особенностями строения исследуемых веществ, в частности, с параметрами кластеров. В качестве объекта изучения впервые выбран 1,3-пропандиол (1,3-ПД). Приведены значения статической диэлектрической проницаемости εs; высокочастотной диэлектрической проницаемости (ε∞DH); дипольного момента кластера μс(D); параметры внутри- (nDH) и межкластерной (mDH) корреляции, времена релаксации τDHexp = 1/ωp, τDHmid и τVFT, кажущаяся энергия активации $\Delta H_{{{\text{DHexp}}}}^{\# }$, $\Delta H_{{{\text{DHth}}}}^{\# }$, $\Delta H_{{{\text{VFT}}}}^{\# }$ 1,3-пропандиола, ZHB = $H_{{\exp }}^{\# }$(T)/23.84 (423 K) – число Н-связей, разрываемых в результате диэлектрической релаксации.

Ключевые слова: статическая и комплексная диэлектрическая проницаемость, параметры кластеров, Н-связи, диэлектрическая релаксация

Изучение структуры и свойств жидких систем относится к фундаментальным задачам современной науки. Понимание молекулярной природы макроскопических свойств жидкостей и растворов важно для разработки научных основ создания жидкофазных материалов с заданным комплексом свойств. Наиболее сложными для теоретического изучения и весьма важными для практического применения являются молекулярные системы с сильно направленными взаимодействиями, к числу которых относятся многоатомные спирты [1]. Решение этих важных вопросов тормозится отсутствием строгой количественной теории жидкого состояния вещества. Последняя, для своего развития, нуждается в целом ряде экспериментальных данных о структуре жидкостей [17]. Современное развитие экспериментальной методики определения ε*(ω) и модельных представлений для интерпретации экспериментальных данных дает возможность приблизиться к пониманию деталей молекулярной структуры жидкостей.

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

(1)
$\begin{gathered} F(\tau ) = \frac{{\sin \beta \pi }}{\pi }{{\left( {\frac{\tau }{{{{\tau }_{\beta }} - \tau }}} \right)}^{\beta }},\quad {\text{е с л и }}\quad \tau \leqslant {{\tau }_{\beta }}, \\ F(t) = 0,\quad {\text{е с л и }}\quad t > {{t}_{\beta }}. \\ \end{gathered} $
Здесь τβ – максимальное время релаксации, β – параметр, характеризующий ширину спектра времен релаксации. Частотная зависимость комплексной диэлектрической проницаемости, соответствующая функции распределения (1) имеет вид:

(2)
$\varepsilon {\text{*}}(\omega ) = {{\varepsilon }_{\infty }} + \frac{{{{\varepsilon }_{s}} - {{\varepsilon }_{\infty }}}}{{{{{(1 + i\omega {{\tau }_{\beta }})}}^{\beta }}}}.$

График зависимости ε'' = f(ε') представляет собой скошенную дугу (рис. 1).

Рис. 1.

График зависимости ε'' = f(ε'), характерный для уравнения (2).

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

В диэлектрической спектроскопии основным параметром, характеризующим изучаемый объект, является эффективное время диэлектрической релаксации (τ = 1/ωр). Величина τ связана со структурой веществ и растворов, но сделать определенные выводы о строении только на основе τ представляется весьма затруднительным. В наших предыдущих работах [917] был осуществлен анализ диэлектрических спектров ряда веществ на основе уравнений, предложенных Диссадо и Хиллом. Интерес к их подходу обусловлен тем, что он позволяет определить на основе релаксационных данных не один параметр (τ), а несколько параметров, зависящих от строения жидкости (μс, m, n), где μс – дипольный момент кластера, m, n – параметры межкластерной и внутрикластерной корреляции [11, 16, 17]. В настоящей работе продолжено развитие подходов, связывающих значения статической (εs) и комплексной (ε*) диэлектрической проницаемости с особенностями строения исследуемых веществ, в частности с параметрами кластеров. В качестве объекта изучения впервые был выбран 1,3-пропандиол (1,3-ПД).

Ранее нами было обнаружено, что в большинстве случаев значения эффективного времени диэлектрической релаксации (τ), рассчитываемые с помощью уравнения Диссадо–Хилла (ДХ) [2, 3] практически совпадают с (τ), определяемыми по традиционному уравнению Дэвидсона–Коула (ДК) [4]. Но уравнение (ДХ), и в этом его преимущество, дает возможность определить параметры, которые можно связать с характеристиками кластеров (m, n, ${{\bar {\mu }}_{{\text{c}}}}$) [11].

Полученные в [1, 6] экспериментальные данные о дисперсии диэлектрической проницаемости ε*(ω) = ε'(ω) + iε"(ω) 1,3-ПД описывались уравнением (ДХ) [2, 3]:

(3)
$\begin{gathered} \frac{{\varepsilon {\text{*}}(\omega ) - {{\varepsilon }_{{\infty {\text{DH}}}}}}}{{{{\varepsilon }_{s}} - {{\varepsilon }_{{\infty {\text{DH}}}}}}} = {{\left( {1 + i\frac{\omega }{{{{\omega }_{p}}}}} \right)}^{{n - 1}}} \times \\ \times \;\frac{{_{2}{{F}_{1}}[1 - n,1 - m;2 - n;{{{(1 + i\omega {\text{/}}{{\omega }_{p}})}}^{{ - 1}}}]}}{{_{2}{{F}_{1}}(1 - n,1 - m;2 - n;1)}}, \\ \end{gathered} $
где ε∞DH – высокочастотный предел области дисперсии ε*(ω) = ε'(ω) + iε"(ω), частота релаксации ωDH = 1/τDH, где τDH – время релаксации, n, m – параметры внутри- и межкластерной корреляции, 2F1[…] гауссовские гипергеометрические функции.

Среднестатистический квадрат дипольного момента кластера $\left\langle {\mu _{{\text{c}}}^{2}} \right\rangle $ 1,3-ПД рассчитывался по уравнению {ДХ} (4) [2]

(4)
$\begin{gathered} {{\varepsilon }_{{\text{s}}}}\, - \,{{\varepsilon }_{{\infty {\text{DH}}}}}\, = \,{{\left( {\frac{{{{\omega }_{{{\text{DH}}}}}}}{\zeta }} \right)}^{n}}\cos \left( {\frac{{n\pi }}{2}} \right)\frac{{_{{{{N}_{{\text{A}}}}}}}}{{{{V}_{{\text{m}}}}}}\left\langle {\mu _{{\text{c}}}^{2}} \right\rangle \frac{{(1 - M_{{\text{e}}}^{2})}}{{kT}} \times \\ \times \;{{\left\{ {1 - \frac{{{{T}_{{\text{c}}}}}}{T}(1 - M_{{\text{e}}}^{2})} \right\}}^{{ - 1}}}\frac{{\Gamma ({\text{1}} - n)}}{m}, \\ \end{gathered} $
где Г(1 – n) – гамма-функция, Tc= Tgd, Me= tanh[(BDH+ kTcMe)/kT] – единичный вектор продольной составляющей дипольного момента кластера, Vm – молярный объем, ζ – частота колебаний молекул внутри кластера.

Параметры Tc= Tgd (Тс – температура, которую обычно выбирают равной температуре фазового перехода, например, равной температуре стеклования Тg) и BDH находились при анализе температурных зависимостей времен релаксации τDHth в рамках модели (ДХ) отдельно для каждого интервала температур [11]:

(5)
$\begin{gathered} {{\omega }_{{{\text{DHth}}}}} = 1{\text{/}}2\pi {{\tau }_{{^{{{\text{DH}}}}}}} = {{v}_{{0{\text{DH}}}}}\exp ( - {{U}_{{{\text{DH}}}}}{\text{/}}kT) \times \\ \times \;\cosh \left( {\frac{{{{B}_{{{\text{DH}}}}} + k{{T}_{{\text{c}}}}{{M}_{{\text{e}}}}}}{{kT}}} \right)\left\{ {1 - \frac{{{{T}_{{\text{c}}}}}}{T}(1 - M_{{\text{e}}}^{2})} \right\}, \\ \end{gathered} $
где ${{v}_{{0{\text{DH}}}}}$ – максимальная скорость перехода через минимальный средний активационный барьер Δ = $E_{{{\text{DH}}}}^{\# }$ + BDH+ кТcMe, k – постоянная Больцмана. Основной параметр – характеристическая температура Tc, полагалась равной температуре стеклования Tgd = 180.28 K, при которой τDH(Tgd) = = τgWLF(Tgd) = 100(c) в уравнения Вильямса–Лэндела–Ферри (ВЛФ) [911], используемого для описания температурной зависимости интегрального времени релаксации τDH. Методика выбора параметра BDH = 2.248 кДж/моль, не зависящего от температуры, рассмотрена в [911]. Cогласно модели (ДХ) параметр BDH= ћζBNA = = hζBNA/(2π), где h = 6.62 × 10–34 Дж с – постоянная Планка, связан с циклической частотой ζB = = 2πBDH/(hNA) = 35.41 × 1012 Гц. Следовательно, частоту fB можно рассчитать с помощью соотношения fB= BDH/(hNA) = 5.64 ТГц.

Полученные результаты представлены в таблице 1 и на рис. 2.

Таблица 1.  

Экспериментальные характеристики

T, K εs ε∞DH mDH nDH τDhеxp, пс μс, D (Tge = = 180.3 K) μс, D (Tg = = 166.4 K) τVFT, пс $\Delta H_{{{\text{DHexp}}}}^{\# }$, $\Delta H_{{{\text{DHth}}}}^{\# }$, $\Delta H_{{{\text{VFT}}}}^{\# }$, ZHB
кДж/моль
293 35.10 3.20 1.000 0.260 495 90.03 97.81 544 37.28 37.56 40.62 1.56
303 33.00 3.19 1.000 0.248 332 75.70 82.12 319 35.84 36.11 38.18 1.50
313 31.40 3.21 1.000 0.237 220 64.71 70.10 200 34.49 34.73 36.08 1.45
323 29.90 3.33 1.000 0.218 135 53.41 57.75 131 33.22 33.44 34.27 1.39
333 28.50 3.31 1.000 0.195 84.5 44.18 47.67 90.4 32.03 32.21 32.69 1.34
343 27.00 3.41 1.000 0.169 62.0 36.90 39.73 64.57 30.90 31.05 31.30 1.30
353 26.06 3.40 1.000 0.150 47.8 32.82 35.27 47.62 29.84 29.94 30.07 1.25
363 24.20 3.32 0.906 0.130 35.7 26.62 28.56 36.10 28.84 28.89 28.97 1.21
383 22.40 3.30 1.000 0.087 23.9 22.09 23.60 22.24 27.00 26.91 27.10 1.13
403 20.90 3.39 0.931 0.068 15.1 18.48 19.70 14.76 25.34 25.11 25.56 1.06
423 19.10 3.46 0.900 0.050 9.25 15.84 16.83 1.04 23.84 23.44 24.28 1.00

Обозначения: εs – статическая диэлектрическая проницаемость; ε∞DH – высокочастотная диэлектрическая проницаемость, μс(D) – дипольный момент кластера, параметры внутри – nDH и mDH межкластерной корреляции, времена релаксации τDHexp = 1/ωp, τDHmid и τVFT, кажущаяся энергия активации $\Delta H_{{{\text{DHexp}}}}^{\# }$, $\Delta H_{{{\text{DHth}}}}^{\# }$, $\Delta H_{{{\text{VFT}}}}^{\# }$ 1,3-пропандиола, ZHB = = $H_{{\exp }}^{\# }$(T)/23.84 (423 K) – число Н-связей, разрываемых в результате диэлектрической релаксации.

Рис. 2.

Графики: а – диэлектрической проницаемости 1,3-пропандиола, б – диэлектрических потерь 1,3-пропандиола, в – статической диэлектрической проницаемости, полученной в работах [1, 6]; 1 – эксперимент, 2 – расчет.

На рис. 3 представлены температурные зависимости εs, τ и μс, полученные из наших измерений [1] и с привлечением данных работы [6] при Tge = 180.3 K. Для сравнения на графике для μс приведен расчет при Tge = 166.4 K, полученной в работе [7] калориметрическим методом.

Рис. 3.

Температурные зависимости τ и μс; а – расчет времени релаксации в рамках модели (DX), б – расчет дипольного момента кластера при разных значениях Tge.

В работе [18] авторы повторили работу известного шведского ученого Гестблома [19, 20], но не сослались на его работы. Так же авторы не сослались на работу [1], а сослались на совершенно другую работу этого автора. Их результаты по фактору корреляции весьма сильно отличаются от работ [1, 19, 20]. Авторы не привели табличные данные по эксперименту, а из графиков получить данные по εs и другим экспериментальным величинам не представляется возможным.

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

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

  1. Журавлев В.И. // Журн. физ. химии. 1992. Т. 66. № 1. С. 225.

  2. Dissado L.A, Hill R.M. // Philos. Mag. B. 1980. V. 41. 625.

  3. Dissado L.A., Hill R.M. // Proc. R. Sos. Lond. A. 1983. V. 390. P. 131.

  4. Davidson D.W., Cole R.H. // J. Chem. Phys. 1951. V. 19. P. 1484.

  5. Ферри Дж. Вязкоупругие свойства полимеров. М.: Изд-во иностр. лит., 1963.

  6. Онищенко С.Г., Исаев И.Д. // Прикладная физика твердого тела. Сб. статей. Махачкала. 1973. С. 61.

  7. Kiyoshi Takeda, Osamu Yamamuro, Itaru Tsukushi et al. // J. Mol. Struct. 1999. V. 479. P. 227.

  8. Фрелих Г. Теория диэлектриков. М.: Изд-во иностр. лит., 1960. С. 247.

  9. Усачева Т.M., Лифанова H.B., Журавлев В.И., Матвеев В.К. // Журн. структур. химии. 2009. Т. 50. С. 967.

  10. Усачева Т.M., Лифанова H.B., Журавлев В.И., Матвеев В.К. // Там же. 2011. Т. 52. № 6. С. 1186.

  11. Журавлев В.И. // Журн. физ. химии. 2019. Т. 93. № 5. С. 707.

  12. Усачева Т.М., Журавлев В.И. // Там же. 2013. Т. 87. № 3. С. 443.

  13. Журавлев В.И. // Там же. 2015. Т. 89. № 12. С. 53.

  14. Журавлев В.И. // Там же. 2016. Т. 90. № 8. С. 1.

  15. Журавлев В.И. // Там же. 2016. Т. 90. № 10. С. 1484.

  16. Усачева Т.M., Лифанова H.B., Журавлев В.И., Матвеев В.К. // Там же. 2017. Т. 91. № 6. С. 991.

  17. Усачева Т.M. // Там же. 2018. Т. 92. № 5. С. 760.

  18. Ganesh V., Mashalkar, Pankaj A. Chalikwar, Ashok C. Kumbharkhane // Physics and Chemistri of Liquids. 2015. V. 53. № 3. P. 307.

  19. Hanna F.F., Gestblom Bo, Soliman A. // Phys. Chem. Chem. Phys. 2000. V. 2. P. 5071.

  20. Hanna F.F., Gestblom Bo, Soliman A. // J. Mol. Liq. 2002. V. 95. P. 27.

Дополнительные материалы отсутствуют.