DOI: 10.6060/tcct.20186103.5662 УДК: 544.3.01+544.32
МОДЕЛИРОВАНИЕ ТЕРМОДИНАМИЧЕСКИХ СВОЙСТВ И ФАЗОВЫХ РАВНОВЕСИЙ В СИСТЕМЕ LiCl - H2O В РАМКАХ МОДЕЛИ ПИТЦЕРА
С.Е. Працкова
Светлана Евгеньевна Працкова
Кафедра аналитической и физической химии, Челябинский государственный университет, ул. Братьев Кашириных, 129, Челябинск, Российская Федерация, 454001 E-mail: [email protected]
В статье приведены результаты термодинамического моделирования системы LiCl - H2O. Для описания термодинамических свойств водных растворов электролита наиболее широко используется модель Питцера. Модельные уравнения просты и содержат небольшое количество переменных параметров, имеющих определенный физический смысл. Это линейные комбинации вириальных коэффициентов в разложении избыточной свободной энергии Гиббсараствора по числу молей ионов. Выведены уравнения для расчета среднего ионного коэффициента активности для компонента системы LiCl и осмотического коэффициента воды. Рассчитаны энергетические параметры модели с учетом опытных данных по осмотическому коэффициенту воды Ф от состава раствора при 25 °С и от температуры для системы LiCl - H2O. Для проверки адекватности модели Питцера была решена обратная задача, то есть были рассчитаны значения осмотического коэффициента воды и среднего ионного коэффициента от состава раствора. Вычисленные значения хорошо согласуются с экспериментальными данными. Определены значения стандартных энергий Гиббса образования кристаллогидратов LiClH2O, LiCl 2H2O. Построена двойная фазовая диаграмма LiCl - H2O. На расчетной диаграмме имеется эвтектическая точка с координатами 25,3 мас. % LiCl, t = -80 °C и две перитектики - 48,2 мас. % LiCl, t = 50 °C; 56,5 мас. % LiCl, t = 102 °C. Произведен расчет избыточных парциальных молярных энтальпий и энтропий компонентов раствора. Как следует из полученных данных, водные растворы хлорида лития образуются с поглощением тепла, система стабилизируется за счет увеличения энтропии. Раствор испытывает положительные отклонения от идеальности.
Ключевые слова: модель Питцера, осмотический коэффициент, средний ионный коэффициент активности, фазовые диаграммы, избыточные парциальные молярные энтальпии и энтропии
MODELING OF THERMODYNAMIC PROPERTIES AND PHASE EQUILIBRIUM IN SYSTEM OF LiCl - H2O WITHIN FRAME OF PITZER MODEL
S.E. Pratskova
Svetlana E. Pratskova
Department of Analytical and Physical Chemistry, Chelyabinsk State University, Bratiev Kashirinykh st., 129, Chelyabinsk, 454001, Russia E-mail: [email protected]
Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol. 2018. V. 61. N 3
Lithium chloride is used as a solid electrolyte in chemical sources of current, organic synthesis, pyrotechnics and as a desiccant. To describe the thermodynamic properties of aqueous electrolyte solutions, the Pitzer's model is most widely used. Model equations are simple and contain a small number of variable parameters that have a certain physical meaning. They are linear combinations of virial coefficients in the expansion of the excess free Gibbs energy of the solution of the number of moles of ions. The results of thermodynamic modeling of the LiCl - H2O system within the frame of this model are presented in the article. Equations for calculating the average ionic coefficient of activity for the component of the LiCl system and the osmotic coefficient of water are derived. The energy parameters of the model are calculated taking into account the experimental data on the osmotic coefficient of water Ф on the composition of the solution at 25 °C and on the temperature for the LiCl - H2O system. To verify the adequacy of the Pitzer's model, the inverse problem was solved, i.e. the values of the osmotic coefficient of water and the average ionic coefficient of the solution composition were calculated. The calculated values are in good agreement with the experimental data. The values of the standard Gibbs energies for the formation of crystalline hydrates LiCH2O, LiCl2H2O are determined. A double phase diagram of LiCl - H2O is constructed. On the calculated diagram, there is a eutecticpoint with coordinates of 25.3 wt.% LiCl, t = -80 °C and two peritectics of 48.2 wt. % LiCl, t = 50 °C; 56.5 wt. % LiCl, t = 102 °C. The excess partial molar enthalpies and entropies of the solution components were calculated. As follows from the data obtained, aqueous solutions of lithium chloride are formed with heat absorption, the system is stabilized by increasing entropy. The solution experiences positive deviations from ideality.
Key words: Pitzer's model, osmotic coefficient, average ionic activity coefficient, phase diagrams, excess partial molar enthalpies and entropies
Для цитирования:
Працкова С.Е. Моделирование термодинамических свойств и фазовых равновесий в системе LiCl - H2O в рамках модели Питцера. Изв. вузов. Химия и хим. технология. 2018. Т. 61. Вып. 3. С. 21-25 For citation:
Pratskova S.E. Modeling of thermodynamic properties and phase equilibrium in system of LiCl - H2O within frame of Pitzer model. Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol. 2018. V. 61. N 3. P. 21-25
ВВЕДЕНИЕ
Хлорид лития используется как твердый электролит в химических источниках тока, в органическом синтезе, пиротехнике и как осушитель. Термодинамические свойства водного раствора хлорида лития являются предметом изучения многих ученых [1-4].
Экспериментальная диаграмма ЫС1 - Н20 содержит кристаллогидраты: ЫСЬШО, ЫО^ШО. Точка эвтектики имеет координаты 25,3 мас. % ЫС1, - 80 °С и перитектические точки - 48,2 мас. % ЫС1, 50 °С; 56,5 мас. % ЫС1, 100,5 °С [1]. В настоящем исследовании рассматриваются справочные данные по термодинамическим свойствам водного раствора хлорида лития [5].
Для описания термодинамических свойств водных растворов электролитов самой широко применяемой является модель Питцера [6,7]. Модельные уравнения просты и содержат небольшое число варьируемых параметров, которые имеют определенный физический смысл. Они представ-
ляют собой линейные комбинации вириальных коэффициентов в разложении избыточной свободной энергии Гиббса раствора по числу молей ионов.
РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ
Выведено уравнение для расчета среднего ионного коэффициента активности ЫС1:
Ыу±=-АФ[т12/(1 + Ьт+ 21п(1 + Ьт-12)/Ь] + 2р°1 а-т +
' (1)
+ р ,С - [(т§'(21 /2) + 2mg (2112)] + т2Сф+,С -,
где АФ - предельный угол наклона Дебая-Хюк-келя для осмотического коэффициента, равный 0,391 [7]; ти+= т = т - моляльность соли;
Ь - постоянная для всех систем, равная 1,2 по Питцеру; с1- ,РР;1+СГ, Сф + сг - энергетические
параметры теории, функции g, g' определяются как g = 2[1-(1+х) • ехр-х)]/х2, / = -2[1 - (1+х+—) • ехр(-х)]/х2,
х = а
41 , где а (а = 2, а = 0) - масштабный множитель для зависимости в от ионной силы.
Уравнение для расчета осмотического коэффициента воды:
Ф = 1 + [-АФт/(1 + Ьт12) + В(о)+ т +
^Ы+,С1-
C ф
Li + ,Cl -
= 0,2587-8,8-10"4-Т (Т = 273 - 298, К),
+ b (1)+ m exp(-2m /2 ) + m 2 С5 +
^ Li+ ,СГ ^ Li + ,СГ
где mL.+ = m = m - моляльность соли.
Расчет энергетических параметров теории Pl* сг, + cr, Сф+ сг проводился методом МНК, для
этого была составлена система уравнений выражения (1), где варьировались справочные данные по моляльности и соответствующего ей среднего коэффициента активности у± при 25 °С [5]. Параметры теории получились следующими:
Рйс - = 0,2°44, P +,С1 - = -0,0291, С и +сг = -0,0039.
Также были рассчитаны энергетические параметры, как функции от температуры, для этого использовались уравнение (1) и данные о значениях коэффициента активности электролита при различных температурах табл. 1.
Таблица1
Средние коэффициенты активности электролита
при различных температурах Table 1. The average coefficients of activity of electro-
2r ф
(2)
Значения среднего коэффициента
m активности у±
T=0°C T=50°C T=75°C T=100 °C
0,1 0,801 0,783 - 0,756
0,2 0,800 0,747 - 0,714
0,5 0,754 0,722 - 0,676
1,0 0,792 0,746 0,715 0,681
2,0 0,959 0,871 0,818 0,763
3,0 1,230 1,071 0,984 0,899
4,0 1,636 1,355 1,217 1,087
6,0 3,117 2,308 1,961 1,665
8,0 6,248 4,051 3,230 2,584
10,0 12,51 6,972 5,190 3,889
Параметры теории получились следующими:
b
(о)
Li + ,Cl-
= -1,9796+0,0073-Т (Т = 273 - 298, К),
в(о)+ = 663,7625-6,1704-Т+0,0191 -Т2-2-10-5
" Li + ,Cl- ' ' '
(Т = 298 - 373 К);
b
bn-ar = 10,2678-0,0346-Т (Т = 273 - 298, К), (1+ = -18,4307+0,1134-Т-1,9-10"4-Т2+4,5-10"8-Т3
Li + ,Cl~ 5 5 5 5
(Т = 298 - 373 К);
= -0,0717+1,5144-10_4-Т+8,5185-Т2-0,2-10_8-Т3
Li + ,Cl-
(Т = 298 - 373 К).
Для проверки адекватности модели Пит-цера решали обратную задачу, т.е. находили значения осмотического коэффициента воды и среднего ионного коэффициента от состава раствора. Для этого подставляли найденные параметры теории в уравнения (1), (2). Полученные данные представлены в табл. 2.
Таблица2
Справочные [5] и расчетные данные осмотического коэффициента воды и среднего коэффициента активности для системы LiCl - H2O при 25 °С Table 2. The reference [5] and calculated data of the osmotic water coefficient and average activity coefficient
m, моль/1000 г H2O фэксп фрасч у^эксп у±расч
0,1 0,942 0,929 0,793 0,769
0,2 0,941 0,925 0,761 0,727
0,3 0,946 0,929 0,747 0,708
0,4 0,953 0,937 0,742 0,700
0,5 0,962 0,948 0,741 0,698
0,6 0,972 0,961 0,744 0,701
0,8 0,993 0,989 0,755 0,715
1,0 1,015 1,019 0,773 0,737
2,0 1,142 1,185 0,921 0,924
3,0 1,285 1,356 1,156 1,217
4,0 1,443 1,524 1,501 1,630
5,0 1,611 1,687 1,995 2,193
6,0 1,786 1,843 2,697 2,945
7,0 1,963 1,993 3,683 3,938
8,0 2,138 2,136 5,047 5,234
9,0 2,306 2,271 6,895 6,905
10,0 2,463 2,399 9,333 9,033
11,0 2,605 2,520 12,45 11,71
12,0 2,728 2,634 16,29 15,04
13,0 2,881 2,740 20,82 19,14
14,0 2,912 2,838 25,93 24,11
15,0 2,972 2,929 31,48 30,06
16,0 3,012 3,012 37,25 37,09
17,0 3,036 3,087 43,13 45,31
18,0 3,049 3,155 49,14 54,73
Как следует из полученных данных, расчетные величины хорошо согласуются с экспериментальными данными, относительная погрешность укладывается в 5%, следовательно, параметры модели подобраны корректно.
Для моделирования диаграммы состояния ЫС1 - Н2О необходимы данные об энергии Гиббса образования кристаллогидратов Ь1С1Н2О, Ь1С12Н2О. Для этого рассмотрено равновесие:
ЫС1(ж)+Н2О(ж) = ЫС1пН2О(т) (3)
Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol. 20l8. V. б!. N 3
Выражение ДG°т для кристаллогидрата: АО°т = ЯТ 1п аис1 (^) + ЯТ 1п ан г0{ ж) (4)
Подставляя выражения (1) в (4) и задаваясь составом системы и температурой, найдены выражения для энергии Гиббса образования кристаллогидратов:
Да°т(ПС1-Н20) = -87000+29,97-Т, Дж/моль, Да°т(ПС1-2Н20) = -64000+23,38-Т, Дж/моль.
Результаты расчета фазовой диаграммы ПС1 - Н20 представлены на рисунке, на который нанесены экспериментальные точки [1]. На рассчитанной диаграмме есть точка эвтектики с координатами 25,3 мас.% ЫС1, 1 = -80 °С и две перитектики 48,2 мас. % ЫС1, 1 = 50 °С; 56,5 мас. % ЫС1, г = 102 °С.
AH2 =-uRT2 • (Э lng± / ЭT)
(5)
m, p v '
Аналогичное выражение для растворителя: Щ = -RT2 • (ЭlnajЭТ)m,p. (6)
Для относительной парциальной молярной энтропии растворенного вещества и растворителя можно записать следующие выражения:
S°m" = S2 - S° = AH¡ / T -uR • ln(Q • m • g± ),(7)
8отн = _ 3о =ан1/ Т - Я ■ 1по1. (8)
Относительные молярные энтальпия и энтропия раствора рассчитаны по уравнениям:
(9)
AH = AHl • xl + AH 2 • x2,
S = S°m" • x, + S°mH ■ x2. (10)
Результаты расчетов термодинамических функций представлены в табл. 3.
Таблица 3
Избыточные относительные парциальные молярные энтальпии и энтропии компонентов раствора и
раствора в целом при 25 °С Table 3. Excess relative partial molar enthalpies and entropies of the components of the solution and solution as a whole at 25 °C
1ЛС12НгОп4 ЫСШгОю
Рис. Расчетная (1) и экспериментальная (2) диаграмма LiCl -
H2O
Fig. The calculated (1) and experimental (2) diagram of LiCl -H2O
Так как определить абсолютное значение энтальпии раствора невозможно, можно установить, насколько энтальпия вещества, находящегося в данном растворе, отличается от энтальпии этого вещества в другом состоянии или в другом растворе. Разность между парциальной молярной энтальпией компонента раствора и парциальной молярной энтальпией компонента в стандартном состоянии называют относительной парциальной молярной энтальпией компонента в данном растворе [8].
Коэффициенты активности связаны с другими парциальными молярными характеристиками растворенного вещества и растворителя фундаментальными термодинамическими соотношениями.
Выражение для расчета относительной парциальной молярной энтальпии растворенного вещества по данным температурной зависимости коэффициента активности [8]:
m, моль/ l000 г H2O Hl, кДж/ моль S нтн Д ж/' моль-К H 2, кДж/ моль S2 нтн Д ж/' моль-К АН, кДж/ моль <^отн Дж/ моль-К
l,0 0,019 0,22 12,82 48,12 0,25 1,06
l,2 0,041 0,32 11,99 41,68 0,29 1,20
l,6 0,10 0,60 9,59 27,35 0,37 1,35
2,0 0,18 0,95 6,71 12,32 0,40 1,34
2,5 0,29 1,45 3,01 -6,05 0,41 1,13
3,0 0,42 2,00 -0,44 -23,0 0,37 0,72
4,0 0,65 3,10 -5,74 -50,4 0,22 -0,50
5,0 0,81 3,99 -8,17 -67,2 0,07 -1,89
6,0 0,83 4,41 -7,27 -72,1 0,03 -3,05
7,0 0,61 4,12 -2,81 -64,5 0,22 -3,57
8,0 0,09 2,87 5,33 -44,2 0,75 -3,05
9,0 -0,78 0,42 17,23 -10,8 1,07 -1,15
l0,0 -2,10 -3,46 32,93 35,6 3,24 2,50
ll,0 -3,92 -9,02 52,46 95,27 5,39 8,22
l2,0 -6,33 -16,5 75,85 168,1 8,27 16,2
l3,0 -9,37 -26,1 103,0 254,2 11,9 27,0
l4,0 -13,1 -38,1 134,2 353,5 16,5 40,6
15,0 -17,7 -52,8 169,2 466,1 22,0 57,5
l6,0 -23,1 -70,3 208,0 592,0 28,5 77,7
l7,0 -29,4 -90,9 250,8 731,2 36,2 101
18,0 -36,8 -114 297,5 883,7 45,0 129
Как следует из полученных данных, водные растворы хлорида лития образуются с поглощением тепла, система стабилизируется за счет увеличения энтропии. Раствор испытывает положительные отклонения от идеальности.
ВЫВОДЫ
В рамках модели Питцера выведены уравнения для расчета среднего ионного коэффициента активности для компонента системы ЫС1 и осмотического коэффициента воды. Рассчитаны энергетические параметры модели с учетом опытных данных по осмотическому коэффициенту воды Ф от состава раствора при 25 °С и от температуры для
системы ЫС1 - Н2О. Рассчитанные значения осмотического коэффициента воды и среднего ионного коэффициента хорошо согласуются с экспериментальными данными. Определены значения стандартных энергий Гиббса образования кристаллогидратов ЬЮШО, ЬЮ^ШО. Построена двойная фазовая диаграмма ЫС1 - Н2О. Произведен расчет избыточных парциальных молярных энтальпий и энтропий компонентов раствора.
ЛИТЕРАТУРА
1. Monnin C., Dubois M., Papaiconomou N., Simonin J.P.
Thermodynamics of the LiCl + H2O System. Chem. Eng. Data. 2002. V. 47. No. 6. P. 1331-1336. DOI: 10.1021/je0200618.
2. Simonin J.P., Blum L., Turq P. Real ionic solutions in the Mean Spherical Approximation. 1. Simple salts in the primitive model. Phys. Chem. 1996. V. 100. P. 7704-7709. DOI: 10.1021/jp953567o.
3. Simonin J.P. Real ionic solutions in the Mean Spherical Approximation. 2. Pure strong electrolytes up to very high concentrations and mixtures, in the primitive model. Phys. Chem. 1997. V. 101. P. 4313-4320.
4. Simonin J.P., Bernard O., Blum L. Real ionic solutions in the Mean Spherical Approximation. 3. Osmotic and activity coefficients for associating electrolytes in the primitive model. Phys. Chem. 1998. V. 102. P. 4411-4417.
5. Cohen-Adad R., Lorimer J.W. Solubility of LiCl in water. In Solubility Data Series. 1991.
6. Працкова С.Е. Фазовые диаграммы и моделирование термодинамических свойств системы H2O - KOH - KCl. Вестн. Башкир. ун-та. 2017. Т. 22. № 3. C. 713-716.
7. Працкова С.Е. Моделирование термодинамических свойств системы H2O-UOH методом Питцера. Научн. Ведом. Белгород. гос. ун-та. Сер.: Естеств. науки. 2017. Т. 267. № 18. Вып. 40. C. 11-14.
8. Семиохин И.А. Физическая химия. М.: МГУ. 2001. 272 с.
REFERENCES
1. Monnin C., Dubois M., Papaiconomou N., Simonin J.P.
Thermodynamics of the LiCl + H2O System. Chem. Eng. Data. 2002. V. 47. No. 6. P. 1331-1336. DOI: 10.1021/je0200618.
2. Simonin J.P., Blum L., Turq P. Real ionic solutions in the Mean Spherical Approximation. 1. Simple salts in the primitive model. Phys. Chem. 1996. V. 100. P. 7704-7709. DOI: 10.1021/jp953567o.
3. Simonin J.P. Real ionic solutions in the Mean Spherical Approximation. 2. Pure strong electrolytes up to very high concentrations and mixtures, in the primitive model. Phys. Chem. 1997. V. 101. P. 4313-4320.
4. Simonin J.P., Bernard O., Blum L. Real ionic solutions in the Mean Spherical Approximation. 3. Osmotic and activity coefficients for associating electrolytes in the primitive model. Phys. Chem. 1998. V. 102. P. 4411-4417.
5. Cohen-Adad R., Lorimer J.W. Solubility of LiCl in water. In Solubility Data Series. 1991.
6. Pratskova S.E. Phase diagrams and modeling of the thermodynamic properties of the H2O-KOH-KCl system. Vest. Bashkir. Univ. 2017. V. 22. N 3. P. 713-716 (in Russian).
7. Pratskova S.E. Modeling of the thermodynamic properties of the H2O-LiOH system by the Pitzer method. Nauch. Ve-dom. Belgorod. Gos. Univ. Ser.: Estesstv. Nauki. 2017. V. 267. N 18. Iss. 40. P. 11-14 (in Russian).
8. Semiokhin I.A. Physical Chemistry. M.: MSU. 2001. 272 p. (in Russian).
Поступила в редакцию 24.07.2017 Принята к опубликованию 29.01.2018
Received 24.07.2017 Accepted 29.01.2018