ISSN 2007-1957 EQUILIBRIO QUÍMICO EN SISTEMAS HOMOGÉNEOS LÍQUIDOS Salvador Pérez Cárdenas Instituto Politécnico Nacional [email protected] Adelina Rosas Mercado Instituto Politécnico Nacional [email protected] Adelina Pérez Rosas Instituto Politécnico Nacional [email protected] Resumen El análisis termodinámico del equilibrio químico en sistemas líquidos puede implicar serias complicaciones, que suelen resolverse haciendo simplificaciones que no siempre son del todo válidas. Por ejemplo, puede tenerse inmiscibilidad parcial o total de los componentes del sistema reaccionante o bien, aún en el caso de miscibilidad total, una desviación del comportamiento de una solución ideal. La simplificación más frecuentemente utilizada es precisamente suponer idealidad de la solución lo que puede dar lugar a cálculos muy alejados del comportamiento real del sistema reaccionante. En este trabajo se tratará el caso de un sistema reaccionante multicomponente, totalmente miscible, pero no considerado como una solución ideal, introduciendo para ello el modelo de solución de Wilson (Grant M. Wilson, 1964)1 que para sistemas totalmente miscibles suele dar muy buenos resultados. Palabras clave: equilibrio químico, sistemas homogéneos líquidos, multicomponente 1 Ejemplar 15. Julio-diciembre de 2016 ISSN 2007-1957 En el estado líquido, la actividad de cada componente de una mezcla homogénea está dada por: fˆi i xi Pi satisat PFi ai 0 i xiisat PFi sat f1 Pi (1) Donde: a i : actividad del componente i en la mezcla líquida, i : coeficiente de actividad de i en la mezcla, x i : fracción mol de i en la fase líquida, isat : coeficiente de fugacidad de i en condiciones de saturación, PFi : factor de Poynting de i. Si la fase vapor puede considerarse una mezcla de gases ideales: isat = 1 (2) Y como el factor de Poynting es, a presiones bajas o moderadas, muy cercano a 1, entonces puede establecerse que: a i i xi (3) que es la forma más usual de calcular la actividad de un líquido en una mezcla. Si además, la mezcla líquida puede considerarse una solución ideal: i 1 (4) a i xi (5) Y, en consecuencia, La constante termodinámica de de equilibrio químico se define como: P a i K R a (6) i donde i es el coeficiente estequiométrico de cada especie química, producto o reactivo, que interviene en la reacción. Usando (3): 2 Ejemplar 15. Julio-diciembre de 2016 ISSN 2007-1957 P x i K i (7) x R i Pero, para una mezcla líquida ideal, según (4), la expresión (7) se reduce a: P x i K R x (8) i Supongamos un caso específico; sea la reacción en fase líquida C2H5OH + CH3COOH <=====> CH3COOC2H5 + H2O que ocurre a 25°C. Supongamos que la mezcla líquida de cuatro componentes es homogénea, es decir, que no hay problemas de inmiscibilidad parcial o total y que además la solución es ideal (lo cual no necesariamente es el caso); en esta circunstancia, los coeficientes de actividad serán 1 para todos los componentes y entonces será válida la ecuación (8). Para tener los valores de fracción mol, habrá que hacer un balance de materia. Supongamos que la proporción de alimentación de los reactivos es de 1 mol de etanol por cada mol de ácido acético; entonces el balance de materia quedara de la siguiente manera: Líquido Moles en equilibrio Fracción mol en equilibrio Etanol 1− (1 − Acido acético 1− (1 − Acetato de etilo Agua Total 2 1 Y entonces: 2 K 1 2 (9) 3 Ejemplar 15. Julio-diciembre de 2016 ISSN 2007-1957 Por otro lado, de tablas de propiedades termodinámicas, puede encontrarse que, a 25°C: H 0f (cal/mol) Líquido G f0 (cal/mol) Etanol − 66,356 − 41,770 Acido acético − 116,400 − 93,800 Acetato de etilo − -112,570 − 77,600 − 68,317 − 56,690 Agua De donde, GR0 P 298 i G f0 R i G f0 298 1,280 cal/mol H2O (10) 298 Y como, GR0 T RT ln K entonces: G R0 K exp RT T 0.1153 (11) Con (9) y (11): 0.1153 2 1 2 (12) de donde, = 0.2535 Con este valor, al sustituir en la última columna de la tabla de balance de materia: Líquido Etanol Acido acético Acetato de etilo Agua xi 0.3732 0.3732 0.1268 0.1268 4 Ejemplar 15. Julio-diciembre de 2016 ISSN 2007-1957 ¿Y si la temperatura no fuera 25°C? Supongamos que la reacción ocurre a 75°C (348.15 K), temperatura inferior a la de ebullición normal (1 atm), de todos los compuestos presentes, lo que garantiza que la fase líquida sea la única presente; Líquido Teb0 (°C) Etanol 78.3 Acido acético 117.9 Acetato de etilo Agua 77.1 100.0 El efecto de la temperatura sobre la constante de equilibrio químico puede evaluarse de la siguiente manera: G R0 H R0 T G R0 T (13) A la expresión (13) se le conoce como ecuación de Gibbs-Helmholtz. Por otro lado, T G R0 G R0 T T2 T G R0 T P T G R0 T H 0 2 R T (14) Con (13) en (14): Que puede desde luego escribirse en la siguiente forma: T G R0 RT Y como, GR0 entonces: T H R0 RT 2 (15) RT ln K G R0 T ln K RT (16) 5 Ejemplar 15. Julio-diciembre de 2016 ISSN 2007-1957 que al sustituirse en (15), lleva a: H R0 ln K T RT 2 Como K es sólo función de T, la ecuación anterior puede escribirse: d ln K H R0 dT R T2 (17) Integrando entre límites, desde una temperatura T1 hasta otra temperatura T2, para cada de las cuales la constante de equilibrio vale K1 y K2, suponiendo que H R0 es constante (procedimiento análogo al que conduce a la ecuación de Clausius-Clapeyron): ln K 2 H R0 K1 R Para nuestro caso, si H R0 H R0 ln 1 1 T1 T2 298 (18) 1869 cal/mol H2O, la ecuación (18) podría escribirse 1869 1 1 0.1153 1.987 298.15 348.15 K De donde, a 75°C, K = 0.1810, que del balance de materia implícito en la ecuación (12), conduce a: 2 0.1810 1 2 De donde, = 0.5798 ¿Y si la temperatura es de 25°C y se considera la mezcla líquida como solución no ideal? Para una mezcla multicomponente, totalmente miscible, el modelo de solución que suele funcionar mejor es el llamado de Wilson (1964)1: x ln i 1 ln x j ij k ki j k x j kj (19) j Donde los valores de ij se denominan parámetros adimensionales, que están relacionados con los parámetros ij , llamados energéticos y que usualmente se dan en cal/mol, mediante: 6 Ejemplar 15. Julio-diciembre de 2016 ISSN 2007-1957 ij ij exp R T Vi L V jL (20) Los valores de ij son los que suelen aparecer en la literatura especializada; cuando i = j, ij = 0. Vi L son los volúmenes molares de cada componente puro en la fase líquida. Para la reacción que se está considerando, si denominamos al etanol como componente 1, al ácido acético como componente 2, al acetato de etilo como 3 y al agua como 4, tenemos, de la literatura especializada, que, ij (cal/mol) j\i 1 2 3 4 1 0 500.00 844.69 382.30 2 600.00 0 400.00 300.00 3 −178.81 700.00 0 550.00 4 955.45 200.00 450.00 0 Los valores en cursiva son estimados, a falta de de datos conocidos; en consecuencia los resultados obtenidos adelante son meramente ilustrativos aunque no muy lejos de la realidad. Los volúmenes molares de los componentes puros en fase líquida en función de la temperatura, están dados por las siguientes expresiones (temperatura en kelvin, volumen en cm3/mol): V1 L 53.7004 0.03111T 0.00016T 2 (21-a) V2L 174.2891 0.66642T 0.000958T 2 (21-b) V3L 82.8668 0.03096T 0.00028T 2 (21-c) V4L 22.8860 0.03641T 0.000069T 2 (21-d) A 25°C: V1 L = 58.65, V2L = 57.58, V3L = 98.53 y V4L = 18.16 cm3/mol. Una forma de estimar la miscibilidad o no de un par de especies químicas, puede ser a través de los parámetros de solubilidad de Hildebrand2, definidos como: U ebi i L Vi 0.5 0.5 Hebi RT Vi L 25C 25C (22) Las Hebi pueden obtenerse fácilmente de las ecuaciones para la presión de saturación de cada componente en función de la temperatura, Antoine por ejemplo: 7 Ejemplar 15. Julio-diciembre de 2016 ISSN 2007-1957 ln Pi sat Ai Bi Tebi Ci (23) dado que, como es bien sabido, d ln Pi sat Hebi dTeb RTeb 2 (24) De esta manera pueden obtenerse los parámetros de solubilidad de cada componente: Comp 1 2 3 4 0.5 cal 3 cm 25C i 12.9174 8.9744 8.9695 23.3746 La miscibilidad de cualquier par puede estimarse observando los siguientes criterios3: Si: | i j | 11: miscibilidad total, 11 | i j | 14.5 : rango de incertidumbre, | i j | 14.5 : inmiscibilidad parcial o total Para el caso considerado: Etanol(1)-ácido acético(2): |12.9174 – 8.9744| = 3.9430 Etanol(1)-acetato de etilo(3): |12.9174 – 8.9695| = 3.9473 Etanol(1)-agua(4): |12.9174 – 23.3746| = 10.4574 Ácido acético(2)-acetato de etilo(3): |8.9744 – 8.9695| Ácido acético(2) agua(4): |8.9744 – 23.3746| = 14.4002 Acetato de etilo(3)-agua(4): |8.9695 – 23.3746| = 14.4051 = 0.0049 8 Ejemplar 15. Julio-diciembre de 2016 ISSN 2007-1957 Aunque los dos últimos pares están dentro del rango de incertidumbre, consultando otras fuentes4 puede constatarse que en ambos casos se da miscibilidad total. Como consecuencia de lo anterior, puede entonces utilizarse el modelo de solución de Wilson para el cálculo de los coeficientes de actividad que, con las ecuaciones (19) y (20) conduce a los siguientes resultados: ij j\i 1 2 3 4 1 1 0.7224 0.2426 0.1625 2 0.2162 1 0.3056 0.1111 3 1.3414 0.5111 1 0.1214 4 0.6441 3.8700 4.5234 1 Con lo que finalmente pueden calcularse los coeficientes de actividad y con ellos las actividades de los compuestos que intervienen en la reacción. Los resultados se resumen en la siguiente tabla: cm 3 Vi L mol i ai 58.65 1.5270 0.5702 0.3734 98.53 1.2617 0.4711 3 0.1266 59.15 1.7830 0.2257 4 0.1266 18.16 1.0550 1.1336 Comp. xi 1 0.3734 2 Y entonces, como conclusión: P a i K R a a3 a 4 0.2257 1.1336 0.1122 a1 a 2 0.5702 0.4711 i Valor de K que coincide con el obtenido a 25°C a partir de las energías libres de Gibbs de formación a la misma temperatura, que es K = 0.1153 [(11)]. 9 Ejemplar 15. Julio-diciembre de 2016 ISSN 2007-1957 ____________ Notas: 1. G. M. Wilson. J. Am. Chem. Soc., vol. 86. pp. 127-130, 1964. 2. J. H. Hildebrand and R. L. Scott, Regular Solutions, Prentice-Hall, USA, 1962. Hay otro parámetro de solubilidad, el de Hansen, mucho mejor para predecir miscibilidad, pero bastante más complicado de calcular que el de Hildebrand; para nuestro propósito este último es suficiente. 3. S. Pérez C., XVII Congreso Nacional de Termodinámica. México, 2002. 4. Handbook of Chemistry and Physics, 90th edition. Chemical Rubber Co. (CRC); Solubility Chart. 10 Ejemplar 15. Julio-diciembre de 2016
© Copyright 2024