SciELO - Scientific Electronic Library Online

 
vol.21 número4Búsqueda de Catalizadores Sólidos Básicos para la Producción de BiodieselAdsorción de Cadmio en Solución utilizando como Adsorbente Material Tobáceo Modificado índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

Compartir


Información tecnológica

versión On-line ISSN 0718-0764

Inf. tecnol. v.21 n.4 La Serena  2010

http://dx.doi.org/10.4067/S0718-07642010000400009 

Información Tecnológica Vol. 21(4), 67-73 (2010) doi:10.1612/inf.tecnol.4370it.09

 

QUIMICA Y APLICACIONES

 

Estudio Teórico de la Hidrólisis de Ésteres Alifáticos y Aromáticos

 

Theoretical Study of Aliphatic and Aromatic Esters Hydrolysis

 

Silvana C. Caglieri y Héctor R. Macaño

Universidad Tecnológica Nacional, Facultad Regional Córdoba, Departamento de Ingeniería Química, Av. Cruz Roja Argentina esquina Maestro López, Ciudad Universitaria (5016), Córdoba - Argentina (e-mail :scaglieri@quimica.frc.utn.edu.ar ; hmacano@quimica.frc.utn.edu.ar)


Resumen

Se realizó un estudio teórico comparativo de reactividad del acetato de etilo, éster alifático y del benzoato de etilo, éster aromático, frente a la reacción de hidrólisis catalizada por ácido, a través del análisis de los intermediarios de ambas reacciones. Se modelaron y determinaron las geometrías de todas las especies que intervienen en ambas hidrólisis y se calcularon los parámetros geométricos y las energías de los reactivos y productos, empleando dos métodos de estructura electrónica, DFT y AM1. Los valores obtenidos se compararon con datos experimentales y de la literatura observándose una gran concordancia. Siguiendo el mismo procedimiento se determinaron los parámetros geométricos y energías de los intermediarios de reacción, observándose una mayor estabilidad para el intermediario alifático, que reportó una energía inferior comparada con la del aromático. La diferencia de energía entre ambos fue de 34.11 kcal /mol.

Palabras claves: hidrólisis ácida, éster alifático, éster aromático, reactividad


Abstract

A comparative theoretical study of reactivity of ethyl acetate hydrolysis, aliphatic ester, and ethyl benzoate hydrolysis, aromatic ester, catalyzed by acid, from the analysis of intermediates in both reactions was carried out. The geometry of all species involved in both hydrolysis was determined and the energies of all reagents and products were calculated, using two methods of electronic structure, DFT and AM1. The calculated values were compared with experimental and literature data finding good agreement. Following the same procedure, the geometric parameters and energies of reaction intermediates were determined, observing greater stability for the aliphatic intermediate, with lower energy compared with the aromatic. The energy difference between them was 34.11 kcal/mol.

Keywords: acid hydrolysis, aliphatic ester, aromatic ester, reactivity


 

INTRODUCCIÓN

El estudio teórico de la hidrólisis ácida de los ésteres alifáticos y aromáticos, cuyos productos son alcoholes y ácidos orgánicos, es de gran interés, no sólo por la utilidad de los productos de reacción dentro de la industria química, sino que desde el punto de vista ambiental dicha reacción constituye una vía para la degradación de los poliésteres, polímeros sintéticos conformados por monómeros ésteres. Han sido reportados estudios sobre síntesis y caracterización de poliésteres alifáticos biodegradables (Pepic et al., 2008), investigaciones sobre biodegradabilidad de copoliésteres alifático-aromático (Chen et al., 2008) y sobre el mecanismo de reacción de la hidrólisis de poliésteres alifáticos (Partini y Pantani, 2007). Todos estos estudios coinciden en que los poliésteres formados por monómeros aromáticos poseen mejores propiedades mecánicas y físicas, pero son menos biodegradables y más resistentes a la degradación por hidrólisis que aquellos poliésteres que sólo contienen monómeros alifáticos, esto lleva a la necesidad de sintetizar copoliésteres con unidades alifáticas y aromáticas para poder satisfacer ambos requerimientos.

Los ésteres se hidrolizan lentamente ya que poseen grupos salientes muy básicos y además el agua es un nucleófilo débil. El aumento de la velocidad de la hidrólisis se puede conseguir con el agregado de un ácido como catalizador. Para que un catalizador aumente la velocidad de una reacción debe aumentar la velocidad del paso lento o paso determinante de la misma. La hidrólisis es una reacción de sustitución nucleofílica sobre carbono insaturado, el ácido que actúa como catalizador, primero protona el oxígeno del carbonilo y posteriormente se produce el ataque nucleofílico del agua para formar el intermediario tetraédrico, etapa determinante de la velocidad de reacción. El ácido aumenta la velocidad de formación del intermediario tetraédrico al protonar el oxígeno del carbonilo, porque un grupo carbonilo protonado es más susceptible de ser atacado por un nucleófilo que un grupo carbonilo no protonado.

Se han llevado a cabo estudios teóricos sobre la hidrólisis de ésteres alifáticos catalizada con ácidos (Hori et al., 2007), con bases (Haeffner et al., 1999) y aplicación de diversos métodos computacionales para estudiar la velocidad de dicha reacción mediante catálisis alcalina (Zhang et al., 2005). Tanto estos estudios, como trabajos experimentales y estudios cinéticos (Calvar et al., 2006) concuerdan en que dicha reacción transcurre a través de la formación de un intermediario tetraédrico. En este artículo se presenta un estudio teórico comparativo de la hidrólisis ácida del acetato de etilo y del benzoato de etilo, cuyos productos de reacción son el etanol y los ácidos acético y benzoico respectivamente, empleando la Química Computacional como herramienta para estudiar las estructuras y propiedades de compuestos orgánicos y predecir comportamientos de reacciones químicas.

METODOLOGÌA

Para poder comparar la reactividad de ambos ésteres frente a la hidrólisis, se diseñaron y optimizaron las estructuras de todas las especies que intervienen en las reacciones, determinándose los parámetros geométricos óptimos correspondientes, calculando además las energías mínimas de todos los compuestos, reactivos y productos, que participan en las reacciones como así también la energía de los respectivos intermediarios. Para realizar estos cálculos se emplearon métodos de estructura electrónica. Estos métodos dentro de la química computacional, utilizan las leyes de la mecánica cuántica como base para los cálculos, mediante la resolución de la ecuación de Schrödinger, que independiente del tiempo toma la forma de la siguiente ecuación:

H Ψ(r) = E Ψ(r)                                                                                                                                        (1)

Siendo H una función diferencial que incluye la energía cinética y potencial de núcleos y electrones, denominada operador hamiltoniano, Ψ la función de onda molecular y E energía. La solución que obtiene el menor valor de energía es la que corresponde al estado basal.

Dentro de los métodos de estructura electrónica, se empleó el método DFT basado en la teoría de funcionales de densidad (Parr y Yang, 1989), el cual incluye efectos de correlación electrónica y dentro de este se empleó el B3LYP (Becke, 1997). Este método denominado ab-initio, considera toda la estructura atómica sin introducir parámetros empíricos en la resolución de la ecuación de Schrödinger. Se optó por la función de base 3-21G* que indica que se utilizan tres funciones gaussianas para representar los orbitales atómicos internos, mientras que los externos, capa de valencia, se representan mediante dos funciones para la parte contraída y una para la parte difusa, más una función de polarización para átomos pesados, indicada mediante el *, que representa el efecto de los orbitales tipo d.

También para calcular las energías mínimas de todas las moléculas que intervienen en ambas reacciones se empleó el método semiempírico AM1, Austin Model 1 (Dewart et al., 1985), que utiliza parámetros derivados de resultados experimentales para simplificar los cálculos cuánticos y reducir su coste computacional. Este método omite los cálculos relativos a los electrones de capa interna, que los trata como si formasen parte del núcleo y se centra únicamente en los electrones externos, considerando a éstos electrones los responsables de la reactividad. Todos los cálculos se realizaron empleando el programa Gaussian03 (Frisch et al., 2004).

RESULTADOS Y DISCUSIÓN

Ambas hidrólisis, tanto la del acetato de etilo, éster alifático como la del benzoato de etilo, éster aromático, transcurren por el mismo mecanismo de reacción. En las Fig. 1 y Fig. 2 se describen las etapas a través de las cuales se llevan a cabo las hidrólisis ácidas de ambos ésteres. Ambas reacciones son reversibles. La hidrólisis ácida de los ésteres es la reacción inversa a la esterificación. En la Fig. 3 se observan las estructuras de los intermediarios de reacción, alifático y aromático de ambas hidrólisis.

La Tabla 1 reporta los resultados obtenidos de la aplicación de los métodos de cálculo mencionados, es decir en la misma se pueden observar las longitudes de enlace (r), los ángulos de enlace (9), los momentos dipolares (µ) y la energía (E) de los reactivos y productos que intervienen en cada una de las reacciones motivo de estudio. Se observa una marcada concordancia entre los valores de energía que resultan de los cálculos teóricos indicados con E en kcal/mol y los valores obtenidos de bibliografía indicados con AHÍ en kcal/mol en fase gaseosa (Yaws, 1999).

Fig. 1: Mecanismo de la hidrólisis ácida del acetato de etilo

Fig. 2: Mecanismo de la hidrólisis ácida del benzoato de etilo

Fig. 3: Estructura intermediarios alifático y aromático

Tabla 1: Parámetros geométricos y energías de reactivos y productos

En las Tabla 2 y Tabla 3 se muestran los parámetros geométricos, valores de longitudes y ángulos de enlaces de los intermediarios alifáticos y aromáticos respectivamente.

Tabla 2: Parámetros geométricos del intermediario alifático

En ambos intermediarios el carbono central tetraédrico está unido a tres oxígenos y a un carbono. La longitud de enlace r (C1 -C4) en el intermediario alifático es de 1.521 Å, superior a la longitud r (C1 - C4) de 1.508 Å que presenta el intermediario aromático. El anillo bencénico es un grupo más voluminoso que el metilo, pero es plano y proporciona un carbono hibridado sp2 para enlazarse al carbono central tetraédrico, esto demuestra su mayor proximidad al carbono central tetraédrico comparado con el metilo que si bien es un grupo más pequeño, proporciona un carbono hibridado sp3 para enlazarse con dicho carbono, de allí que se encuentre más alejado. La existencia y proximidad de un grupo más voluminoso ocasiona un mayor impedimento estérico y los cálculos de energía a nivel teórico indican que el intermediario alifático es 34.11 kcal/mol más estable que el intermediario aromático.

Tabla 3: Parámetros geométricos del intermediario aromático

Si analizamos el mecanismo de la reacción de sustitución nucleofílica catalizada por ácido, luego de la protonación del oxígeno del carbonilo debida al medio, el carbono del carbonilo experimenta una deficiencia electrónica que en el caso del acetato de etilo está compensada por el efecto inductivo sp3→ sp2 ejercido por el grupo metilo y en el caso del benzoato de etilo, la densidad de carga positiva puede deslocalizarse en el anillo bencénico. Esta deslocalización de la carga compensa más la deficiencia electrónica que el efecto debido a la diferencia de hibridación, por lo tanto en el caso del acetato de etilo, el carbono del carbonilo está más susceptible a ser atacado por el agente nucleofílico, siendo por lo tanto mayor la reactividad de este éster frente a la hidrólisis catalizada por ácidos. Este análisis se convalida al comparar las respectivas energías de activación (Ea), obteniéndose las mismas a partir de las especies protonadas. Para la hidrólisis ácida del acetato de etilo se obtuvo una Ea= 16.39 kcal/mol, valor que se aproxima al de 15.6 kcal/mol obtenido experimental mente (Kirby, 1972) y para la hidrólisis ácida del benzoato de etilo se obtuvo una Ea= 24.45 kcal/mol, valor que se aproxima al de 23.82 kcal/mol calculado a partir del valor experimental de 19.26 kcal/mol, reportado para la reacción inversa de esterificación (Pipus et al., 2000).

CONCLUSIONES

De acuerdo a los resultados y discusión presentados en este trabajo, el intermediario tetraédrico alifático resulta ser más estable que el aromático, por lo tanto el acetato de etilo presenta una mayor reactividad frente a la hidrólisis ácida, comparado con el benzoato de etilo, presentando este último una menor velocidad de reacción. La facilidad con que los ésteres alifáticos se hidrolizan coincide con la aplicación de esta reacción en la degradación de poliésteres alifáticos.

La existencia de una marcada concordancia entre los resultados teóricos obtenidos y los datos bibliográficos, valida el empleo de los métodos teóricos DFT y AM1 como herramientas para el estudio de la estructura de los reactivos y productos que participan del sistema químico planteado. La aproximación entre los resultados teóricos obtenidos y datos experimentales, valida el empleo de dichos métodos para estimar la cinética del sistema químico planteado.

REFERENCIAS

Becke, A.D.; Density-functional thermochemistry V. Systematic optimization of exchange-correlation functionals, J. Chem. Phys.: 107, 8554-8560 (1997).        [ Links ]

Calvar, N., B. González y A. Domínguez, Esterification of acetic acid with ethanol: Reaction kinetics and operation in a packed bed reactive distillation column,Chemical Engineering and Processing: 46, 1317-1323 (2007).        [ Links ]

Chen, Y., L. Tan, L. Chen, Y.Yang y X. Wang, Study on Biodegradable aromatic/aliphatic copolyesters, Brazilian Journal of Chemical Engineering: 25, 321-335 (2008).        [ Links ]

Dewar, M.J.S., E.G.Zoebisch, E.F.Healy y J.J.P.Stewart; AM1: A New General Purpose Quantum Mechanical Molecular Model, J. Am. Chem.Soc.:107, 3902-3909 (1985).        [ Links ]

Frisch, M.J., G. W. Trucks, y otros 80 autores; Gaussian 03, Revision E.01; Gaussian, Inc., Wallingford CT, Pittsburgh, USA (2004).        [ Links ]

Haeffner, F., C.H. Hu, T. Brinck y T. Nori; The catalytic effect of water in basis hydrolysis of methyl acetate a theoretical study, J. Mol. Struct.Theochem: 459, 85-93 (1999).        [ Links ]

Hori, K., Y. Ikenaga, K. Arata, T. Takahashi, K. Kasai, Y. Noguchi, M. Sumimoto y H. Yamamoto; Theoretical study on the reaction mechanism for the hydrolysis of esters and amides under acidic conditions, Tetrahedron: 63, 1264-1269 (2007).        [ Links ]

Kirby, A.J., Hydrolysis and Formation of Esters of Organic Acids, In Comprehensive Chemical Kinetics, pp 57-161 Elsevier Publishing Company, Amsterdam, The Netherlands (1972).        [ Links ]

Parr, R.G. y W. Yang; Density-functional theory of atoms and molecules, Oxford University Press.-Nueva York (1989).        [ Links ]

Partini M. y R. Pantani, FTIR analysis of hydrolysis in aliphatic polyesters, Polymer Degradation and Stability 92, 1491-1497 (2007).        [ Links ]

Pepic, D., E. Zagaz, M.Zigon, A.Krzan, M. Kunaver y J. Djonlagic; Synthesis and characterization of biodegradable aliphatic copolyesters with poly (ethylene oxide) soft segments, European Polymer Journal: 44,904-917 (2008)        [ Links ]

Pipus,G., I. Plazl y T. Koloini; Esterification of benzoic acid in microwave tubular flor reactor, Chemical Engineering Journal: 76,239-245 (2000).        [ Links ]

Yaws , C.L., Chemical Properties Handbook, Mc Graw Hill, Nueva York (1999).        [ Links ]

Zhang, H., X. Qu y H. Ando; A simple method for reaction rate prediction of ester hydrolysis, J. Mol. Struct.Theochem: 725, 31-37 (2005).        [ Links ]

Creative Commons License Todo el contenido de esta revista, excepto dónde está identificado, está bajo una Licencia Creative Commons