Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 0,940

MOLECULAR ASSOCIATION IN BINARY MIXES OF ALIPHATIC ETHERS AND 1-ALCOHOLS

Rudakov A.M. 1 Glagoleva M.A. 1 Lipanova N.V. 1 Khmelevskaya N.D. 1 Skorobogatko D.S. 1
1 National Research Nuclear University "Moscow Engineering Physics Institute”
In frameworks the cluster model of molecular association and the ideal associated solutions theory are proved the equations of the dependences on concentration of the thermodynamic properties (a component coefficient activities, the excess thermodynamic functions, vapour pressure) for binary mixes nonelectrolytes, showing positive deviations from the Raoult’s law). The key parameters of the model equations are the association number A1 and D dispersion of cluster distributions on stoihiometrical numbers which the pair above solutions pays off on the data on the general(common) pressure. The equations of model it is used for the description of the thermodynamic properties binary mixes aliphatic ethers with alcohols. The conclusion is made, that nonideality the investigated binary mixes is caused, mainly, by molecular association of alcohols. It is established, that in a dimethyl ether the values А1 and D decrease with growth of molecular weight of alcohol.
molecular association
1-alcohols
aliphatic ether
thermodynamic function

Введение

Добавки эфиров и спиртов к моторному топливу широко используются для повышения октанового числа. Поэтому значительное внимание в литературе уделяется изучению равновесного давления пара тройных систем спирт – эфир - алкан и соответствующих бинарных подсистем [1,4–9]. Смеси низкомолекулярных спиртов и эфиров неограниченно смешиваются друг с другом и проявляют положительные отклонения от закона Рауля. Можно полагать, что неидеальность этих смесей определяется ассоциацией спирта.

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

Теоретическая часть

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

, i = 1, 2,… (1)

Число молекул, присоединившееся к произвольно выбранной молекуле в кластере, определим выражением:

, (2)

где Аi – число стехиометрической ассоциации. Присутствие в растворе кластеров с различным стехиометрическим составом равносильно утверждению о существовании распределения кластеров по стехиометрическим числам ассоциации. Математическое ожидание этого распределения по определению определяется соотношением:

, , (3)

где хi и ni – мольная доля и количество растворенного вещества, содержащего i молекул в кластере, – суммарное количество растворенного вещества, – количество чистого растворителя в одном килограмме, моль.

Можно полагать, что рассматриваемый раствор является идеальным ассоциированным [1]. В этом случае выражение термодинамической константы равновесия ki для процессов (1) имеет вид:

. (4)

Из (3) и (4) следует, что:

, (5)

где х0 – мольная доля мономеров растворенного вещества.

После логарифмирования и дифференцирования (5) получим:

, (6)

где D – дисперсия распределения кластеров по числам ассоциации. Разложение правой части (6) в ряд по степеням x0 при ограничении нулевым приближением дает:

, , (7)

где D1 = D(x = 1) и A1 = A(x = 1) – дисперсия и среднее распределение кластеров по числам ассоциации в стандартном состоянии при x = 1. Интегрирование (7) в интервале (, x0), где – мольная доля мономеров растворенного вещества при x = 1, приводит к

. (8)

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

и . (9)

Из (9) следует, что при x ≈ 1 . Тогда из (8) имеем:

. (10)

Пусть для приготовления раствора было взято ns моль растворителя и такое количество N, моль, мономерных молекул растворенного вещества, что после наступления равновесия в растворе оказалось n моль различных кластеров растворенного неэлектролита. Очевидно, что N = n(1+A). В этом случае для коэффициента активности растворителя fs будет справедливо:

. (11)

После простых преобразований имеем:

, (12)

или, учитывая (10) и малую по сравнению с 1 величину Ах, находим:

. (13)

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

. (14)

В свою очередь определение избыточной мольной энергии Гиббса и соотношения (13) и (14) позволяют записать:

(15)

Подставляя (15) в уравнение Гиббса – Гельмгольца, для избыточной мольной энтальпии находим:

. (16)

В случае, когда D1 зависит от Т, с использованием определения дисперсии дискретного распределения получим:

, (17)

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

Для избыточного мольного объема в соответствии с его определением получено следующее выражение:

, (18)

– изменение стандартного мольного объема, сопровождающее присоединение молекулы неэлектролита к кластеру. В предположении об идеальности паровой фазы с использованием (13) и (14) можно записать следующее выражение для давления пара над раствором:

, (19)

где и – давление насыщенного пара над чистыми – растворителем и растворенным веществом, соответственно. Полученные уравнения использованы для описания литературных экспериментальных данных [6,7,10] по общему давлению пара и избыточным свойствам бинарных смесей спирт – эфир. При описании этих данных минимизацию проводили с использованием стандартных отклонений σ:

, (20)

где YE – избыточное термодинамическое свойство или давление пара, v – количество экспериментальных точек, b – число эмпирических параметров в уравнении.

Результаты и обсуждение

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

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

Таблица 1

Описание экспериментальных данных по равновесию жидкость – пар

уравнением (19) (растворитель эфир)

T, K

v

Источник

Параметры уравнения

А1

D1

, кПа

Диметиловый эфир – метанол

323.15

34

[5]

0.71

0.55

7.66

Диметиловый эфир – этанол

323.15

37

[5]

0.55

0.46

7.98

Диметиловый эфир – пропан-1-ол

323.15

37

[5]

0.37

0.36

9.37

Диметиловый эфир – бутан-1-ол

323.15

39

[5]

0.27

0.27

9.96

Дипропиловый эфир – пропан-1-ол

298.15

15

[2]

1.14

0.82

0.06

Дипропиловый эфир – пропан-2-ол

298.15

14

[2]

1.18

0.90

0.05

Избыточные энтальпия и энтропия. В литературе имеются сведения об энтальпиях смешения при различных составах рассматриваемых систем [7,10]. Результаты, представлены в табл. 2.

Таблица 2

Описание экспериментальных данных по энтальпии уравнением (17) (растворитель эфир)

Источник

v

T, K

, кДж/моль

, Дж/моль

Диметиловый эфир – метанол

[5]

20

323.15

-1.75

-0.88

5.3 (6.3%)

Диметиловый эфир – этанол

[5]

21

323.15

-4.59

-0.04

4.1 (0.9%)

Диметиловый эфир – пропан-1-ол

[5]

21

323.15

-7.44

0.10

5.9 (1.3%)

Диметиловый эфир – бутан-1-ол

[5]

21

323.15

-10.63

0.11

5.1 (1.1%)

Дипропиловый эфир – пропан-1-ол

[2]

13

298.15

-2.46

-0.41

13.2 (2.1%)

Дипропиловый эфир – пропан-2-ол

[2]

15

298.15

-3.45

2.71*10-3

16.5 (2.5%)

Можно видеть, что уравнение (17) кластерной модели удовлетворительно описывает экспериментальные данные при выборе в качестве растворителя эфира.

Наличие данных по GE и HE позволяет провести расчет ТSE по определению ТSE = HE - GE в двух вариантах. В первом варианте ТSE рассчитывали с использованием значений GE и HE, вычисленных по уравнениям (15) и (17). Эти результаты на рис. 1 и 2 представлены сплошными линиями. Во втором варианте использовали экспериментальные данные по HE, а значения GE рассчитывали по уравнению (15). Результаты на рис. 1,2 представлены точками.

Описание: 1

Рис. 1. Зависимости от мольной доли спирта избыточных термодинамических функций для системы диметиловый эфир – метанол при 323.15 К: 1 – HE, 2 - GE, 3 – TSЕ

Рис. 2. Зависимости от мольной доли спирта избыточных термодинамических функций для системы диметиловый эфир – этанол при 323.15 К: 1 – HE, 2 - GE, 3 – TSЕ

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

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

Из данных табл. 2 по значениям (энергия, высвобождающаяся при распаде 1 моля ассоциатов спиртов), учитывая значение RT, можно сделать вывод об энергетике ассоциатов. Величина для большинства систем находится в пределах от 1.7 до 10.6 кДж/моль (табл. 2). Принято считать, что энергия водородной связи находится в пределах от 5 до 30 кДж/моль. Следовательно, ассоциаты спиртов в рассмотренных эфирах образуются за счет водородных связей и ван-дер-ваальсовских взаимодействий.

Рецензенты:

Соловьев Сергей Николаевич, д.х.н., профессор, зав. каф. общей и неорганической химии Федерального государственного бюджетного образовательного учреждения высшего профессионального образования «Российский химико-технологический университет им. Д. И. Менделеева», г. Москва.

Аллахвердов Грант Рантович, д.х.н., профессор, Институт фармацевтических реактивов РЕФАРМ, г. Москва.