CC BY
23
УДК 544.3.03
И. П. Анашкин, А. В. Клинов, Н. Р. Христолюбов
СУЩЕСТВОВАНИЕ ЛИНИИ ZENO ДЛЯ ИОННЫХ ЖИДКОСТЕЙ И ЕЕ СВЯЗЬ С КРИТИЧЕСКОЙ ТОЧКОЙ
Ключевые слова: ионная жидкость, линия Zeno.
Показано, что в области температур, где экспериментально измерены PVT данные ионных жидкостей, линия Zeno с высокой точностью может быть описана прямой. По положению линии Zeno были рассчитаны свойства в критической точке и проведено сравнение с имеющимися литературными данными.
Keywords: ionic liquid, Zeno line.
It is shown that ionic liquids Zeno line with high precision can be described by a straight line in the temperature region where the PVT data experimentally measured. Critical properties were calculated using the line Zeno and were compared with the available literature data.
Введение
Ионные жидкости обладают рядом специфических свойств, таких как низкое давление насыщенных паров, высокая растворяющая способность и др. Благодаря этим свойствам перспективным является их применение как растворителей, катализаторов и теплоносителей. Ионные жидкости обладают достаточной для практического применения температурой разложения. Так, ионные жидкости на основе имидазола разлагаются в диапазоне 240-520 °С [1].
Это не позволяет экспериментально изучить всю фазовую диаграмму. Несмотря на то, что на практике большинство процессов химической технологии проводятся при температурах ниже температуры разложения ионных жидкостей, информация о свойствах в данной области может быть использована, например, для расчета линии фазового равновесия [2,3] или нахождения параметров потенциалов межмолекулярного взаимодействия [4,5]. В данной работе было определено положение линии Zeno [6,7] на фазовой диаграмме для различных ионных жидкостей и проведена оценка свойств в критической точке.
Расчет положения критической точки
В работах [2,3] предложен метод для расчета парожидкостного равновесия, основанный на предположении, что линия фазового равновесия асимптотически стремится к линии Zeno и пересекает в точке с нулевой температурой. Выражение для бинодали жидкости было представлено в виде:
1
(
2
P(T) = Pc +- Pc "PB + 3pB
Tb
л
т +
3
+ —
2
(
Tc
PB - Pc - PBT^ TB
(1)
где т=1-Т/Тс , Тв - температура Бойля, Тс -критическая температура, рв - плотность Бойля, рс -критическая плотность. По данному выражению можно определить линию фазового равновесия по температуре и плотности Бойля, критической температуре и критической плотности. Также была предложена связь, позволяющая сократить
^ + = S1.
(2)
требуемые исходные данные:
Лс
Тв Рв
где Б1 - параметр, слабо зависящий от типа вещества. Так, на примере множества различных потенциалов с различными параметрами притяжения и отталкивания, было показано, что 8^0,67^0,71. Данное выражение позволяет сократить количество необходимых
экспериментальных данных, требуемых для нахождения положения линии фазового равновесия по уравнению (1).
Температуры разложения ионных жидкостей не позволяют проводить экспериментальное определение положения критической точки на фазовой диаграмме, поэтому прямое использование корреляций (1) и (2) для определения фазового равновесия затруднительно. Для определения критической точки было использовано уравнение:
Tc , 2Pc
= 1.
(3)
tb Pb
На фазовой диаграмме данная зависимость представляет из себя медиану треугольника, образованного линией Zeno и осями координат. Обработка экспериментальных данных показала, что данная линия проходит в области критической точки [8]. Исходя из предположения, что критическая точка находится на пересечении линий, описываемых уравнениями 2 и 3, критические значения температуры и давления можно вычислить по температуре и плотности Бойля:
Tc = 2S1Tb - Tb , (4)
Pc = (1 - Si)Pb . (5)
Значение параметра S1=0,71 было найдено минимизацией относительного отклонения от экспериментального значения ряда веществ. В таблице 1 представлено сравнение экспериментальных и расчетных данных. Отклонения от экспериментальных данных для большинства веществ не превышают 10 %. Отклонения выше среднего значения наблюдаются для металлов и веществ с выраженным дипольным моментом (такие как вода и аммиак).
2
т
Таблица 1 - Сравнение расчетных свойств в критической точке с экспериментальными значениями [2]
Вещество Экспериментальные значения Расчет Относительное отклонение
Tc, K pc, г/см3 Tb, K pB, г/см3 Tc, K pc, г/см3 ДТс,% Apo, %
Ar 150,69 0,536 392,84 1,87 164,99 0,542 9,49 1,18
Ne 44,49 0,482 118,85 1,63 49,92 0,473 12,20 1,93
Kr 209,48 0,909 537,98 3,24 225,95 0,940 7,86 3,37
Xe 289,73 1,1 740,02 3,95 310,81 1,146 7,28 4,14
NH3 405,4 0,225 935,92 0,95 393,09 0,276 3,04 22,44
CO2 304,13 0,468 741,37 1,8 311,38 0,522 2,38 11,54
этан 305,13 0,207 779,35 0,74 327,33 0,215 7,27 3,67
этилен 282,35 0,214 714,12 0,78 299,93 0,226 6,23 5,70
фтор 144,44 0,593 385,06 2,01 161,73 0,583 11,97 1,70
гексан 507,82 0,233 1235,84 0,9 519,05 0,261 2,21 12,02
метан 190,56 0,163 498,06 0,57 209,19 0,165 9,77 1,41
N2 126,19 0,313 327,62 1,1 137,60 0,319 9,04 1,92
O2 154,58 0,436 401,14 1,53 168,48 0,444 8,99 1,77
пропен 365,57 0,223 894,42 0,86 375,66 0,249 2,76 11,84
CF3Cl (R13) 302,3 0,583 763,31 2,13 320,59 0,618 6,05 5,95
CF2QH (R22) 369,3 0,524 907,89 1,99 381,31 0,577 3,25 10,13
CF2H2 (R32) 351,26 0,424 823,92 1,74 346,05 0,505 1,48 19,01
Cs 1924 0,39 4114,47 1,96 1728,08 0,568 10,18 45,74
Hg 1751 5,8 6552,6 14,4 2752,09 4,176 57,17 28,00
H2O 647,3 0,32 1268 1,2 532,56 0,348 17,73 8,75
Среднее 9,82 10,11
Линия Zeno и критические параметры ионных жидкостей
Для определения положения линии Zeno проведена обработка экспериментальных PVT данных различных ионных жидкостей. На рис. 1 представлены линии Zeno для некоторых из исследуемых жидкостей. Из рисунка видно, что в области, где имеются экспериментальные данные, линия Zeno хорошо аппроксимируется прямой. Температура и плотность Бойля определялись экстраполяцией полученных прямых. Значения критической температуры и критической плотности определялись по уравнениям 4 и 5 соответственно. Результаты расчетов представлены в таблице 2. Также в таблице 2 представлены значения из работы [9], полученные по уравнению состояния с использованием метода групповых вкладов.
500-1
450
400
350
300
250
о [C2mim]|Tf2N]
□ [C8mim][PFe]
❖ [C8mim]|Tf2N]
а [C3mpip][Tf2N]
о
♦ о
д
а Д а а
о о
2.6 2.8 3 3.2 3.4 3.6 3.8 4 4.2 4.4 р, моль/л
Рис. 1 - Линии Zeno ионных жидкостей
Таблица 2 - Температуры и плотности Бойля ионных жидкостей; свойства в критической точке определенные по корреляциям (4), (5) и по уравнению состояния [9]
Ионная жидкость Tb, К pB, моль/л ссылка Tc, К pc, моль/л Tc [9], К pc[9], моль/л
[C2mim] [Tf2N] 1860,5 4,648 [10] 781,4 1,348 - -
[C2eim] [Tf2N] 1769,5 4,523 [11] 743,2 1,312 - -
[C3mim] [Tf2N] 1824,0 4,376 [12] 766,1 1,269 - -
[C5mim] [Tf2N] 1841,8 3,887 [12] 773,5 1,127 - -
[C6mim] [Tf2N] 1894,9 3,660 [13] 795,9 1,061 - -
[C7mim] [Tf2N] 1839,2 3,499 [10] 772,5 1,015 - -
[C8mim] [Tf2N] 1815,8 3,337 [14] 762,7 0,968 - -
[C10mim] [Tf2N] 1767,6 3,395 [11] 742,4 0,984 - -
[C3mpy] [Tf2N] 1819,7 4,185 [15] 764,3 1,214 - -
[C3mpip] [Tf2N] 1858,5 4,000 [15] 780,6 1,160 - -
[C2mim] [BF4] 2257,4 8,516 [14] 948,1 2,470 596,2 1,849
[C4mim] [BF4] 2192,8 6,176 [16] 921,0 1,791 643,2 1,527
[C4mim] [BF4] 2143,3 6,061 [17] 900,2 1,758 - -
[C6mim] [BF4] 2162,3 5,243 [14] 908,1 1,521 690,0 1,300
[C2mim] [PF6] 2166,1 8,233 [14] 909,8 2,388 674,0 1,543
[C4mim] [PF6] 1992,6 5,681 [17] 836,9 1,647 719,4 1,311
[C4mim] [PF6] 2287,6 5,542 [16] 960,8 1,607 - -
[C8mim] [PF6] 2102,1 4,946 [14] 882,9 1,434 810,8 1,009
[C4mim] [Cl] 2196,0 7,630 [18] 922,3 2,213 789,0 1,758
[C2mim] [MeSO4] 2121,8 6,779 [11] 891,1 1,966 - -
[C2mim] [EtSO4] 2271,7 6,061 [19] 954,1 1,758 1067,5 1,516
[C2mim] [EtSO4] 2124,9 6,159 [11] 892,4 1,786 - -
[C2mim] [CF3SO4] 1931,4 5,594 [15] 811,2 1,622 - -
[C4mim] [C(CN)3] 2002,6 4,111 [10] 841,1 1,192 - -
Выводы
Как было определено в предыдущей работе [4], температура и плотность Бойля коррелирует с силой межмолекулярного взаимодействия и размером молекулы соответственно. Найденные по экспериментальным данным температуры и плотности Бойля позволяют качественно оценить влияние химического строения на межмолекулярное взаимодействие. Сила межмолекулярного взаимодействия в большей степени зависит от аниона ионной жидкости. Так, ионные жидкости анионом которых является
бис(трифторметилсульфонил)имид, вне зависимости от катиона обладают, меньшей температурой Бойля по сравнению с ионными жидкостями на основе других анионов. С увеличением размера катиона и аниона возрастает плотность Бойля. Таким образом, можно разработать процедуру, аналогичную работам [4,5], позволяющую определять параметры потенциалов межмолекулярного взаимодействия по линии Zeno.
Свойства в критической точке, определенные по линии Zeno и по уравнению состояния [9], значительно отличаются. Данное различие обусловлено несколькими факторами. Во-
первых, экспериментальное измерение PVT данных возможно в небольшом диапазоне температур, это не позволяет провести верификацию уравнений состояния при высоких температурах. Во-вторых, в данной работе для определения критических свойств использовалось несколько приближений, апробированных на широком классе веществ: одноатомных, двухатомных газах, углеводородах и т. д. В следующих работах планируется получить аналогичные корреляции, разработанные по данным молекулярного моделирования заряженных частиц.
Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 1403-00251 а.
Литература
1. S. Sowmiah, V. Srinivasadesikan, M.-C. Tseng, Y.-H. Chu, Molecules, 14, 9, 3780-3813 (2009)
2. E.M. Apfelbaum, V.S. Vorob'ev, J. Phys. Chem. B, 112, 41, 13064-13069 (2008)
3. E.M. Apfelbaum, V.S. Vorob'ev, J. Chem. Phys., 130, 21, 214111 (2009)
4. И.П. Анашкин, А.В. Клинов, Вестник Казанского Технологического Университета, 15, 8, 273-276 (2012)
5. И.П. Анашкин, А.В. Клинов, Журнал Физической Химии, 87, 11, 1805-1813 (2013)
6. D. Ben-Amotz, D.R. Herschbach, Isr. J. Chem., 30, 59-68 (1990)
7. И.П. Анашкин, А.В. Клинов, Вестник Казанского Технологического Университета, 20, 11-15 (2011)
8. E.M. Apfelbaum, V.S. Vorob'ev, G.A. Martynov, J. Phys. Chem. B, 110, 16, 8474-8480 (2006)
9. J.O. Valderrama, R.E. Rojas, Ind. Eng. Chem. Res., 48, 14, 6890-6900 (2009)
10. R.L. Gardas, M.G. Freire, P.J. Carvalho, I.M. Marrucho, I.M.A. Fonseca, A.G.M. Ferreira, J.A.P. Coutinho, J. Chem. Eng. Data, 52, 5, 1881-1888 (2007)
11. L.I.N. Tomé, P.J. Carvalho, M.G. Freire, I.M. Marrucho, I.M.A. Fonseca, A.G.M. Ferreira, J.A.P. Coutinho, R.L. Gardas, J. Chem. Eng. Data, 53, 8, 1914-1921 (2008)
12. J.M.S.S. Esperanza, Z.P. Visak, N.V. Plechkova, K.R. Seddon, H.J.R. Guedes, L.P.N. Rebelo, J. Chem. Eng. Data, 51, 6, 2009-2015 (2006)
13. J. Safarov, R. Hamidova, S. Zepik, H. Schmidt, I. Kul, A. Shahverdiyev, E. Hassel, J. Mol. Liq., 187, 137-156 (2013)
14. R. Taguchi, H. Machida, Y. Sato, R.L. Smith, J. Chem. Eng. Data, 54, 1, 22-27 (2009)
15. R.L. Gardas, H.F. Costa, M.G. Freire, P.J. Carvalho, I.M. Marrucho, I.M.A. Fonseca, A.G.M. Ferreira, J.A.P. Coutinho, J. Chem. Eng. Data, 53, 3, 805-811 (2008)
16. H. Machida, Y. Sato, R.L. Smith, Fluid Phase Equilibria, 264, 1-2, 147-155 (2008)
17. A. Tekin, J. Safarov, A. Shahverdiyev, E. Hassel, J. Mol. Liq., 136, 1-2, 177-182 (2007)
18. H. Machida, R. Taguchi, Y. Sato, R.L. Smith, Jr., J. Chem. Eng. Data, 56, 4, 923-928 (2011)
19. H. Schmidt, M. Stephan, J. Safarov, I. Kul, J. Nocke, I.M. Abdulagatov, E. Hassel, J. Chem. Thermodyn., 47, 68-75 (2012)
© И. П. Анашкин - программист, каф. процессов и аппаратов химической технологии КНИТУ, anashkin.ivan@gmail.com; А. В. Клинов - д-р техн. наук, проф., зав. каф. процессов и аппаратов химической технологии КНИТУ, alklin@kstu.ru; Н. Р. Христолюбов - магистрант каф. общей химической технологии, nikryd@mail.ru.
© I. P. Anashkin - programmer, department of chemical engineering KNRTU, anashkin.ivan@gmail.com; A. V. Klinov - doctor of technical science, professor, head of department, alklin@kstu.ru; N. R. Khristolubov - student KNRTU, nikryd@mail.ru.
