scieee AI-readable full text Open interactive document viewer

Desarrollo de nanoemulgeles con aceite esencial de mandarina ecológica y goma guar

Antúnez Domínguez, Jesús Manuel

Full text

Índice general 1. Objetivos Generales 4 2. Fundamento teórico 5 2.1. Química Verde. Aceite esencial de Mandarina . . . . . . . . . . . . . . . . . . . 5 2.2. Emulsiones ....................................... 6 2.2.1. Termodinámica de las emulsiones . . . . . . . . . . . . . . . . . . . . . . 7 2.2.2. Parámetros de estabilización . . . . . . . . . . . . . . . . . . . . . . . . . 8 2.2.2.1. Tamañodegota........................... 8 2.2.2.2. Estabilizantes. El sistema HLB . . . . . . . . . . . . . . . . . . 10 2.2.2.3. Modificadores. Geles . . . . . . . . . . . . . . . . . . . . . . . 13 2.3. Mecanismos de desestabilización . . . . . . . . . . . . . . . . . . . . . . . . . . . 14 2.3.1. Separación gravitacional . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 2.3.2. Floculación................................... 16 2.3.3. Coalescencia.................................. 16 2.3.4. Maduración de Ostwald . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17 2.3.5. Inversióndefase................................ 17 2.4. Reología......................................... 18 2.4.1. Comportamientos ideales . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 2.4.2. Comportamientos reales . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 2.4.2.1. Dependencia con la velocidad de cizalla . . . . . . . . . . . . 20 2.4.2.2. Dependencia con el tiempo de cizalla . . . . . . . . . . . . . . 21 2.4.2.3. Comportamiento viscoelástico . . . . . . . . . . . . . . . . . . 21 2.4.3. Ensayosreológicos .............................. 22 3. Desarrollo Experimental 25 3.1. Equipamiento...................................... 25 3.1.1. Materiales ................................... 25 3.1.2. Dispositivos .................................. 25 1 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales 3.1.3. Medidasdeseguridad............................. 27 3.2. Procedimientos ..................................... 27 3.2.1. Formulación .................................. 27 3.2.2. Preparación de nanoemulgeles . . . . . . . . . . . . . . . . . . . . . . . . 28 3.2.3. Ensayos..................................... 29 3.2.3.1. Medidas de la DTG . . . . . . . . . . . . . . . . . . . . . . . . 29 3.2.3.2. Ensayos reológicos . . . . . . . . . . . . . . . . . . . . . . . . . 29 3.2.3.2.1. Barrido de esfuerzo: . . . . . . . . . . . . . . . . . . 29 3.2.3.2.2. Barrido de frecuencia: . . . . . . . . . . . . . . . . . 29 3.2.3.2.3. Curva de flujo: . . . . . . . . . . . . . . . . . . . . . . 29 3.2.3.3. Medidas de estabilidad . . . . . . . . . . . . . . . . . . . . . . 30 4. Resultados y Discusión 31 4.1. MedidasdelaDTG................................... 31 4.2. Ensayosreológicos................................... 35 4.2.1. Barridos de esfuerzo. Cálculo del esfuerzo crítico . . . . . . . . . . . . . 35 4.2.2. Barrido de frecuencias . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38 4.2.3. Curvas de flujo. Ajuste al modelo de Cross . . . . . . . . . . . . . . . . . 38 4.2.4. Medidas de estabilidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40 5. Conclusión 45 6. Nomenclatura 46 7. Referencias 48 Referencias........................................... 48 2 Resumen Mediante la síntesis de nanoemulgeles, se consiguió mejorar la utilidad y la seguridad del aceite esencial de mandarina ecológica para su uso en la nueva generación de productos de la química verde. Varios aspectos se tuvieron en cuenta para asegurar la optimización de la estabilidad del formato en el tiempo: • Se tomó el tamaño de gota como parámetro fundamental de control, alcanzando el tamaño nanómetrico en las gotas de la fase dispersa de aceite esencial usando el homogeneizador a presión Micofluidizer. • Se valoró el efecto de dos agentes tensioactivos en conjunto, el polisorbato 80, también conocido como Tween 80, y el monooleato de sorbitano conocido como Span 80, de tal manera que se compararon las propiedades resultantes para diferentes proporciones relativas de cada uno respecto a la formulación general. • Se reforzó el sistema a través del uso de goma guar de grado alimentario, que no solo proporciona estructura a la emulsión, sino que también altera las propiedades reológicas del conjunto para hacerlas más atractivas. Finalmente, una vez preparados los nanoemulgeles, fueron caracterizados para determinar el tamaño de gota, las propiedades reológicas a través de barridos de deformación, de frecuencia y curvas de flujo; y se estudiaron los mecanismos de desestabilización mediante la técnica de Multiple Light Scattering (MLS). Palabras clave:aceite esencial de mandarina, emulsiones, emulgeles, reología, tensioactivos. 3 Sección 1 Objetivos Generales Los objetivos perseguidos con la realización de este trabajo son: 1. Obtención de nanoemulgeles estables, a partir del aceite esencial de mandarina ecológica y goma guar. 2. Optimización del tamaño de gota y estabilidad según la proporción de tensioactivos Tween/Span empleados en la emulsión. 3. Estudiar la reología de los nanoemulgeles resultantes para la caracterización de sus propiedades y estructura. 4. Identificar los mecanismos de desestabilización para cada una de las formulaciones sintetizadas. 4 Sección 2 Fundamento teórico 2.1. Química Verde. Aceite esencial de Mandarina Una de las motivaciones esenciales en este trabajo está destinadao a su aplicación en el desarrollo de productos renovables y más respetuosos con el medio ambiente. Este es uno de los cometidos principales de la química verde, rama en pleno desarrollo de la química en la que la prioridad de la investigación y producción reside en la sostenibilidad medioambiental. La base de la química verde reside en principios relacionados no sólo con aspectos ecológicos, sino también sociales y económicos. ( Anastas y Warner, 1998). Favoreciendo la prevención de la generación de residuos en lugar de una posterior gestión, siendo además cualquiera de los componentes involucrados sustancias de toxicidad nula o reducida. En cuanto al origen de los componentes, debe ser preferiblemente renovable y su explotación sostenible. En el procesado del producto, se debe fomentar el ahorro energético y material, evitando el uso de sustancias que no formen parte del resultado final y al término de su vida, los componentes deben degradarse en sustancias inocuas y cuyo impacto no persista en el medio ambiente. En la elaboración de las emulsiones de este trabajo se tienen en cuenta varios de estos preceptos. Al provenir del cultivo de la mandarina, el aceite esencial usado es renovable, es más, por su denominación de ecológico, se asegura la sostenibilidad en su producción. Se trata de un compuesto no tóxico y regulado como “sustancia generalmente reconocida como segura” (GRAS) ( FDA, 2017), aunque en laboratorio requiere ciertas medidas de seguridad para su manipulación a la vez que tiene riesgos de peligrosidad asociados a su inflamabilidad. Sus efectos en el medio ambiente son de corta duración, y aunque es tóxico para la vida acuática; su baja solubilidad, alta volatilidad y biodegradabilidad, lo hacen seguro para el medio ambiente (Nikfar y Behboudi, 2014). El aceite esencial de mandarina tiene reconocidas muchas cualidades de interés en la industria química: antifúngico y antibacteriano, insecticida y repelente, antiproliferativo, quimioprotector, antioxidante y por supuesto, posee excepcionales características organolépticas ( Dosoky y Setzer, 2018). Por ello, es un componente de gran interés en aplicaciones medicinales, industriales y alimentarias. Su componente mayoritario es el D-limoneno, seguido del γ -terpineno, disolventes 5 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales de gran interés en la química verde ; ahí su potencial uso en el desarrollo de nuevos productos (Ciriminna R, Lomeli-Rodriguez M, Demma Carà P, Lopez-Sanchez JA, 2014). El aceite esencial proviene del procesamiento de la cáscara de la mandarina (Citrus nobilis oCitrus reticulata) mediante la técnica de prensado en frío, que permite obtener los extractos a través de grandes compresiones sin elevar la temperatura. Por lo tanto, la ventaja de este proceso respecto a otros comunes, como la destilación, reside en una mejor conservación de los elementos volátiles y susceptibles de desnaturalizarse u oxidarse a altas temperaturas que están presentes en la cáscara. Estos fenómenos también ocurren, aunque en menor medida, a temperatura ambiente; siendo necesario garantizar las condiciones de conservación adecuadas para su correcto almacenamiento: un recipiente hermético, en lugar resguardado de la luz y a temperatura constante preferiblemente baja, incluso refrigerado. De lo contrario, sus propiedades se verán considerablemente mermadas. Por consiguiente, la estabilidad de sus componentes es un factor clave en la utilidad del aceite esencial. Por estas razones, la aplicación directa del aceite esencial se vuelve entonces inviable en la mayoría de los casos, y se requiere la síntesis de un formato lo que separe y conserve durante su uso de forma sencilla. Además, como se dijo anteriormente, su uso en estado puro puede resultar peligroso por su inflamabilidad y efectos irritantes al contacto. Por ello, se aprovecha su uso en emulsiones que permiten el mantenimiento de sus características durante mayor tiempo además de disminuir la cantidad necesaria del aceite, suponiendo una optimización del recurso y un aumento de la seguridad. 2.2. Emulsiones Una emulsión es un tipo de coloide, esto es, una mezcla homogénea formada por dos fases inmiscibles. Está formado por una fase continua que contiene las partículas de la fase dispersa distribuidas por todo su volumen. El que la mezcla sea homogénea conlleva que estas partículas o gotas de la fase dispersa tengan tamaño micrométrico o menor, siendo imposible diferenciar los dos componentes a simple vista. En el caso de las emulsiones, ambas fases se encuentran en estado líquido, por lo que una de ellas suele ser una fase orgánica y la otra acuosa. Dependiendo de cuál sea la fase continua, las emulsiones se clasifican en emulsiones O/W (oil in water) si es el agua o W/O (water in oil) si por el contrario lo es el aceite. El primer parámetro de importancia de una emulsión es entonces la fracción de fase dispersa respecto al total, que puede calcularse en volumen φ o en masa φm , estando ambas relacionadas por: 6 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales φm=φρ2 φρ2+(1−φ)ρ1 (2.1) Existen también sistemas más complejos en los cuales gotas de una emulsión son a su vez dispersadas en una fase continua formando emulsiones múltiples (Figura 2.1). En el caso de tres componentes se forman emulsiones O/W/O o W/O/W, en las que la fase dispersa original puede o no coincidir con la fase continua final.(Rosano, Gandolfo, y Hidrot, 1998) Figura 2.1: Clasificación de las emulsiones simples y múltiples según su composición. 2.2.1. Termodinámica de las emulsiones Una de las principales características de las emulsiones es que son termodinámicamente inestables y con el paso del tiempo las fases tienden a separarse. Sin embargo, si pueden llegar a ser cinéticamente estables y mantenerse por períodos prolongados dadas ciertas circunstancias, que serán desarrolladas más adelante en el apartado 2.2.2. La termodinámica de las emulsiones procede del cambio en la energía libre del sistema. El cambio de una única intercara entre los líquidos inmiscibles, que consideramos esférica, respecto a un sistema en el que la fase dispersa se divide en varias esferas o gotas, viene dado por: ∆Gemul =∆Aσ−T∆S(2.2) Donde la contribución que proviene del aumento de área por la tensión superficial σ que 7 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales soporta se compensa con el aumento de la entropía. En la mayoría de casos, la suma suele ser positiva, lo que conlleva que la formación de emulsiones sea un proceso no espontáneo que requiere el suministro de energía ( T. F. Tadros, 2013). Aun así, existen excepciones como es el caso del efecto Ouzo, en el que la fase dispersa se encuentra previamente disuelta en un disolvente que sí es soluble en la fase continua, provocando la saturación y separación del primero al mezclarlos y formando una emulsión espontánea (Vitale y Katz, 2003). Una emulsión tiende a descomponerse por uno de los varios procesos de desestabilización existentes en ausencia de la participación de ninguna forma de estabilización. Por el contrario, cuando están presentes, se forman valles de estabilidad en la energía libre que dificultan y retrasan su degradación (Figura 2.2). Figura 2.2: Transición de energía libre de un estado II (emulsión) a I (separación) sin y con mecanismos de estabilización respectivamente, donde cada curva representa una forma de desestabilización del sistema. En la segunda figura se aprecia la aparición de un estado intermedio estable III (T. F. Tadros, 2013). 2.2.2. Parámetros de estabilización 2.2.2.1. Tamaño de gota Uno de los parámetros principales el comportamiento de una emulsión es el tamaño de gota de la fase dispersa, cuya influencia afecta no solo a la reología y a la apariencia exterior, sino también en su estabilidad. Por este motivo, el tamaño de gota de las emulsiones desarrolladas va a ser una de las variables de control en el procedimiento de desarrollo. En base a su tamaño de gota, las emulsiones se suelen clasificar en macroemulsiones y nanoemulsiones. En las primeras, las gotas tienen un tamaño superior a una micra, mientras que en las otras siempre se encuentran en el intervalo comprendido entre un nanómetro y un micrómetro (Sjöblom, 2006). Un tamaño de gota menor supone un efecto más débil de la fuerza de gravedad sobre la partícula, y por lo tanto, las gotas seguirán recorridos diferentes entre sí por el movimiento 8 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales Figura 2.7: Estructura molecular de la goma guar empleada como modificador de la textura de los nanoemulgeles.(Mudgil y cols., 2014) rrirá rápidamente. Es por eso que la separación de las fases es un proceso dinámico, relacionado con el movimiento y la interacción de las gotas dispersas entre sí. Dependiendo de las condiciones en las que se produzca, se diferencian los siguientes procesos de desestabilización (Figura 2.8), aunque sus efectos están estrechamente conectados y la aparición de uno frecuentemente favorece el desarrollo de los otros, siendo bajo ninguna circunstancia excluyentes. 2.3.1. Separación gravitacional Los principales generadores de movimiento de las gotas dispersas en la emulsión en reposo son la acción de la fuerza de gravedad y el movimiento browniano. Mientras que la primera tiende a agruparlas por su diferencia de densidad respecto a la fase continua, la segunda es aleatoria y promueve una distribución uniforme por todo el volumen disponible. La fuerza gravitatoria tenderá a agrupar el líquido más denso en la parte inferior del contenedor y al más liviano en la superior. Si la fase dispersa es más densa se produce la sedimentación de las gotas en el fondo, mientras que si es menos densa se produce el cremado y se forma una capa de fase dispersa por encima de la emulsión. Este comportamiento viene descrito por la ley de Stokes: νstokes = − 2gr 2¡ρ1−ρ2¢ 9η1 (2.5) Esta velocidad límite, la velocidad de cremado, resulta de la atracción de la gravedad y del rozamiento con el medio disperso. El resultado de la confluencia de las gotas en una zona de la muestra genera la diferenciación de fases a simple vista: una capa de gotas unidas o emulsión rica en gotas conocida como sedimento o crema, una zona intermedia que desaparece 15 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales progresivamente de la emulsión inicial, y una parte mayoritariamente acuosa pobre en gotas dispersas llamada serum. 2.3.2. Floculación La floculación es un proceso de agregación de las gotas en cúmulos, y aunque esta unión no es total y puede ser reversible, sí contribuye a la separación gravitacional aumentando el tamaño efectivo de las gotas en emulsiones diluidas, esto es, de baja fracción dispersa φ . Sin embargo, si la emulsión es muy concentrada, las gotas pueden formar un entramado que impida el movimiento y consecuente unión entre ellas, aumentando la estabilidad general de la emulsión de forma similar a un gel. La floculación de una emulsión está dictada por las interacciones entre las gotas, de modo que las fuerzas de repulsión y atracción entre ellas se igualan en una distancia de equilibrio que mantiene sus superficies unidas. A este efecto hay que añadir la acción de fuerzas externas, que por un lado pueden aumentar la frecuencia con la que las gotas colisionan y se unen, o si son suficientemente elevadas, separar los agregados formados mediante esfuerzos de cizalla. 2.3.3. Coalescencia En este proceso la unión de las gotas es permanente y comprende su fusión en una sola. Al igual que en el caso de la floculación, la coalescencia supone un aumento en la distribución de tamaños de gota, pero además, también implica un descenso en la superficie de contacto entre los líquidos. Por consiguiente, este mecanismo es el que más acerca al sistema al estado termodinámicamente estable de una única interfase y suele estar ligado a la etapa final de otros mecanismos de desestabilización. Pueden diferenciarse dos tipos de coalescencia: la producida por la colisión de dos gotas de forma inmediata; y la que se produce por el contacto prolongado entre dos gotas. Ambas están relacionadas con debilidades en la interfase que producen la fusión de las gotas, no obstante, en el primer caso es evitable fácilmente con el uso de tensioactivos que provoquen la repulsión entre gotas, mientras que en el segundo aparece espontáneamente por fluctuaciones a pesar de dicha repulsión. El primer tipo da lugar a la separación gravitacional mientras que el segundo está condicionado por el contacto entre las gotas, producido tanto por la cremado (o sedimentación) como por la floculación o la acción de fuerzas externas. 16 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales 2.3.4. Maduración de Ostwald La maduración de Ostwald es un proceso por el cual se produce la transferencia de masa desde las gotas de menor tamaño a la de mayor tamaño, desencadenando el crecimiento de las últimas. Este proceso tiene su origen en las diferencias de presión de la intercara que viene dadas por la presión de Laplace (2.3). La separación entre las dos fases no es totalmente abrupta y existe cierto grado de progresión en la concentración de un medio a otro, que se refleja en cierta solubilidad del líquido disperso en los alrededores de la gota. La solubilidad dependerá de la presión sobre la superficie, que es mayor para gotas menores. Como consecuencia, una mayor parte de la masa de las gotas de menor tamaño se encuentra solubilizada en la intercara y es susceptible de su transporte y absorción por parte de una gota mayor. La velocidad a la que ocurre la difusión causante de maduración de Ostwald viene dada por la teoría LSW (Lifshitz, Slezov y Wagner), la cual aunque formulada para sólidos, describe adecuadamente el fenómeno en líquidos: ωostwald =dr 3 c dt =8 9 DC∞σV2 m RgT(2.6) Resulta llamativa la existencia de un radio crítico r c , al cual las gotas no experimentan ni reducción ni crecimiento, este valor fluctúa con el tiempo y experimentalmente ha demostrado aproximarse al radio medio de las gotas de la emulsión ( Taylor, 1995). A causa de la maduración de Ostwald se producen desigualdades en la distribución de tamaños de gota y se fomenta la aparición de distribuciones multimodales. 2.3.5. Inversión de fase La inversión de fase es un cambio en la estructura de la emulsión, intercambiando las posiciones entre la fase dispersa y el medio continuo, esto es, una emulsión O/W pasa a ser W/O o viceversa. Este tipo de cambios está estrechamente relacionado con el tensioactivo usado, y dependiendo de las condiciones en las que se produce el proceso de inversión de fase se pueden distinguir dos tipos (Salager, 2016): • Transicional: alteraciones en la temperatura provocan una variación de la afinidad de los grupos hidrófilos y lipófilos de un tensioactivo, resultando en un cambio en el índice HLB y la existencia de la conocida como la temperatura de inversión de fase (PIT), en la cual se produce el cambio reversible de una emulsión a la otra (McClements, 1999). • Catastrófica: la inversión se produce de forma dinámica, añadiendo de forma continua más cantidad de la fase dispersa para elevar su concentración súbitamente. 17 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales Este fenómeno también ha suscitado interés por su capacidad de aplicación en el proceso de emulsificación, permitiendo crear emulsiones a partir de otras más sencillas. Figura 2.8: Representación del aspecto de cada uno de los mecanismos de desestabilización de emulsiones (T. F. Tadros, 2013). 2.4. Reología La reología es la ciencia dedicada al estudio del flujo y la deformación de los materiales en relación a las tensiones ejercidas. En el entorno industrial, el conocimiento de estas propiedades es de vital importancia en, por un lado, el procesamiento y envasado del material y por el otro, en las características típicamente deseadas por los consumidores de las industrias alimentaria y cosmética. La deformación se define como la variación de la forma de un material sin afectar a su integridad tras aplicar un esfuerzo, mientras que el flujo es la continuidad de dicha deformación en el material durante el período que se mantiene la tensión. Existen dos formas en las que se puede producir el flujo ante una fuerza, cortante o extensional (Figura 2.9). En el primero, las partículas del cuerpo se desplazan en planos paralelos al esfuerzo, es decir, deslizan entre sí. En el segundo, el desplazamiento ocurre axialmente, y puede significar el acercamiento o 18 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales alejamiento de las partículas.(Barnes, 2000) Figura 2.9: Representación de las diferencias entre esfuerzos cortantes o de cizalla y esfuerzos extensionales en el flujo que producen (Barnes, 2000). 2.4.1. Comportamientos ideales Los materiales pueden exhibir dos comportamientos bien diferenciados respecto a las deformaciones que se ejercen sobre ellos (McClements, 1999): • Sólidos ideales: son capaces de almacenar la energía transmitida por el esfuerzo en los enlaces de sus partículas, de manera que una vez este cese, el material recuperará su forma original. Se conocen como sólidos elásticos y su proporcionalidad entre la tensión aplicada y la deformación resultante sigue la ley de Hooke: τ=Eγ(2.7) La constante de proporcionalidad es el módulo elástico, cuyo significado físico es la rigidez del material, que no es más que su capacidad para almacenar energía. La ecuación (2.7) es útil solo en el caso de un material homogéneo e isótropo, y, dependiendo el tipo de esfuerzo, el módulo elástico será diferente: módulo de Young para tensiones o compresiones simples, módulo de cizalla para tensiones de cizalla y módulo de compresibilidad para compresión uniforme. Sin embargo, en la realidad, los materiales solo son capaces de actuar de forma elástica por deformaciones reducidas, tras lo cual se produce la deformación plástica y flujo o, por el contrario, se pierde la integridad del material y se produce la ruptura. • Líquidos ideales: la energía aplicada por el esfuerzo se disipa e invierte en el desplazamiento de las partículas y en la ruptura de los enlaces de la estructura, que se rompen y crean con facilidad. Se produce entonces el flujo del material durante el periodo de aplicación de la fuerza y una vez esta cesa también lo hace el flujo. No existe recuperación alguna de la energía y ni de la forma inicial. Se conocen como fluidos newtonianos y siguen la relación: τ=η˙ γ(2.8) 19 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales La constante de proporcionalidad se conoce como viscosidad, y su valor da una idea de la resistencia que ejerce el líquido al flujo. Aquí, la ecuación (2.8) exige que el líquido sea incompresible, homogéneo e isótropo. Aunque existe la viscosidad extensional, este comportamiento se asocia habitualmente a esfuerzos de cizalla ( McClements, 1999). Como en el caso anterior, la viscosidad es difícilmente constante y varía con parámetros como la temperatura. 2.4.2. Comportamientos reales Más interesante desde el punto de vista reológico, resultan los cambios en la viscosidad relacionados con el esfuerzo producido en los fluidos no newtonianos: 2.4.2.1. Dependencia con la velocidad de cizalla La viscosidad puede variar dependiendo de la velocidad de cizalla (Figura 2.10), se diferencian los siguientes tipos de materiales: • Materiales pseudoplásticos: reducen su viscosidad conforme aumenta la velocidad a la que se produce la deformación. Este comportamiento es común en emulsiones por mecanismos como la floculación, ya que al unirse las gotas provocan un aumento aparente de la viscosidad, pero, al ejercer un esfuerzo estas se separan y el valor de la viscosidad disminuye. • Materiales dilatantes: su viscosidad se incrementa al aumentar la velocidad de deformación. Es común entre geles y polímeros con largas cadenas unidas puntos de enredo y entrecruzamientos que van perdiendo libertad de movimiento al fluir. • Plásticos de Bingham: la dependencia entre esfuerzo y velocidad de deformación es lineal, no obstante, no existe el flujo de materia hasta que no se supera cierto umbral de esfuerzo, asociado a la pérdida de estructura por parte del material. Figura 2.10: Representación de los comportamientos dependientes de la velocidad de cizalla para la viscosidad y el esfuerzo respectivamente. 20 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales 2.4.2.2. Dependencia con el tiempo de cizalla También existen diferentes comportamientos en relación con cómo afecta a la viscosidad la duración del esfuerzo. Aquí, la velocidad de cizalla es constante, pero, aun así, se aprecia un cambio en el valor de la viscosidad. Este cambio es reversible cuando desaparece la tensión, aunque igualmente requiere un tiempo de recuperación. • Tixotropía: la viscosidad disminuye cuanto mayor sea el tiempo en el que se ha ejercido el esfuerzo. • Antitixotrpía: también conocida como reopexia, consiste en el aumento de la viscosidad conforme mayor es la duración de la tensión. Figura 2.11: Representación de los comportamientos dependientes del tiempo de cizalla, formando ciclos de histéresis cuyo área es equivalente a la pérdida de energía. Características de estos comportamientos son curvas de histéresis al retroceder en el barrido de velocidades de cizalla (Figura 2.11), puesto que existen pérdidas de energía por cambios reversibles en la estructura del material (Grumezescu, 2016). 2.4.2.3. Comportamiento viscoelástico Los comportamientos de sólido y líquido ideales se alejan por lo tanto de lo que se encuentra en la naturaleza y se consideran más bien como casos extremos. En el caso de muchos sistemas materiales, se produce un comportamiento intermedio, con cierto grado de elasticidad y cierto grado de viscosidad. Esto se conoce como viscoelasticidad, y en numerosos casos, este comportamiento puede aproximarse por la suma de dos contribuciones, una perfectamente elástica, representada por un muelle, y otra perfectamente viscosa, representada por un amortiguador (Figura 2.12) (Grumezescu, 2016). Figura 2.12: Modelos simples basados en muelles y amortiguadores para la descripción de comportamiento viscoelástico. 21 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales 2.4.3. Ensayos reológicos La viscoelasticidad es el principal dominio de estudio de la reología, para su caracterización, se llevan a cabo ensayos reológicos. La caracterización reológica se realiza a través de los ensayos dinámicos de cizalla oscilatoria en un reómetro. En los barridos de esfuerzo, se ejerce una señal de tensión sinusoidal con frecuencia (2.9) y se recorren valores de baja amplitud (SAOS). La deformación del material será entonces el objeto de medida (2.10): τ=τ0cos(ωt)(2.9) γ=γ0sen(ωt+δ)(2.10) Existe un desfase δ entre la variable de control y su respuesta que da información sobre el comportamiento general del sistema. Para sólidos ideales elásticos su valor es 0 y para líquidos viscosos ideales es de 90 ° , el comportamiento viscoelástico comprende entonces todo el intervalo entre ambos valores (Figura 2.13). Figura 2.13: Comportamientos de los modelos ideales respecto al viscoelástico en la relación de esfuerzos y deformaciones oscilatorias. A partir de este valor, y debido a la forma sinusoidal, es posible la separación de las contribuciones elástica y viscosa a través de las funciones seno y coseno, desfasadas 90 ° al igual que las respuestas de un sólido y un líquido ideales: G0=µγ0 τ0¶cosδ(2.11) G"=µγ0 τ0¶senδ(2.12) tgδ=G" G0(2.13) 22 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales Estas magnitudes son los módulos viscoelásticos, siendo el primero el módulo elástico o de almacenamiento (2.11) y el segundo el módulo viscoso o de pérdidas (2.12). Ambos constituyen el módulo complejo (2.14). G∗=G0 +iG"(2.14) En una representación de los módulos frente a la deformación, para valores bajos, ambos serán constantes y paralelos en la conocida como región lineal viscoelástica (LVR). Dependiendo de cuál de los dos módulos sea mayor, se puede hablar de un sólido o un líquido viscoelástico, pues predomina el comportamiento elástico o viscoso respectivamente. Llegado el valor de tensión de cizalla conocido como tensión crítica τc , este comportamiento se pierde debido a la destrucción de la estructura interna, que se traduce en una caída súbita del módulo elástico G’. Posteriormente, también le ocurre a G” aunque es posible que antes se alcance un punto de corte entre ambos módulos, el punto de flujo, que indica el cambio a un comportamiento más fluido en sólidos viscoelásticos (Figura 2.14). Figura 2.14: Comportamiento típico de un sólido y un líquido viscoelásticos para un barrido de los esfuerzos, apreciándose en ambos casos la meseta de la LVR. Otra forma de distinguir un comportamiento más sólido o fluido se halla en otro ensayo dinámico, el barrido de frecuencias. La variable a recorrer es la frecuencia y se representa la forma de los módulos al respecto en la LVR, esto es, sin llegar a producir efectos destructivos en la muestra (Figura 2.15). Aún así, el comportamiento viscoelástico de un material puede no limitarse a un solo régimen, y cambiar con la frecuencia, de ahí el interés por este tipo de ensayos. Para conocer el comportamiento del viscoelástico del material, también es necesaria la caracterización de la viscosidad. Como se vio en la sección 2.4.2.1, la viscosidad puede depender de la velocidad de deformación, y su relación se muestra en las curvas de flujo. Se realizan a través de medidas del par de torsión sobre la muestra al incrementar la velocidad de cizalla, 23 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales Figura 2.15: Comportamiento de los módulos elástico y viscoso para diferentes sistemas viscoelásticos en función de la frecuencia. procurando mantener el estado estacionario para evitar efectos dependientes del tiempo. En el caso de los nanoemulgeles, lo habitual es un comportamiento pseudoplástico. Figura 2.16: Comportamiento de la viscosidad en algunos materiales pseudoplásticos comprendida entre dos valores límite. Este comportamiento, se distingue por partir de un valor inicial estable a bajas velocidades de cizalla, que se identifica como el intervalo en el que se consigue mantener un orden irregular en la estructura por la acción del movimiento browniano. Al aumentar la velocidad de cizalla se produce un descenso marcado por la separación de agregados y la ordenación de las gotas y partículas en dirección del esfuerzo, disminuyendo las interacciones entre estos elementos que contribuyen a la viscosidad. No obstante, este nivel de ordenación tiene un límite alcanzable, y la viscosidad tiende a estancarse entorno a un valor límite η∞ . Este comportamiento es posible describirlo a través del modelo de Cross (2.15) ( Santos, Trujillo-cayado, y Calero, 2014). η=η0 1+³˙ γ ˙ γC´1−n(2.15) Las constantes presentes en (2.15) tienen además significados físicos que aportan información del sistema; ˙ γc representa la velocidad crítica a la que se desvía el comportamiento pseudoplástico del material del aquel propio de un fluido newtoniano. La viscosidad inicial del sistema antes de que se produzca su disminución es η0 ynun parámetro llamado índice de flujo cuyo valor está comprendido entre cero y la unidad, correspondientes a los extremos del comportamiento viscoelástico, sólido ideal y líquido newtoniano respectivamente. 24 Sección 4 Resultados y Discusión 4.1. Medidas de la DTG Para empezar, la distribución del tamaño de gota sirvió como variable principal para distinguir el procesado óptimo con el Microfluidizer, para ello, se ensayaron muestras con una proporción de tensioactivos al 50/50. Figura 4.1: Influencia del procesado de las emulsiones en la DTG. La distribución de gotas menores se corresponde con una sola pasada por el homogeneizador Microfluidizer a una presión de 25000 psi (Figura 4.1). Además, este proceso no solo se 31 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales 15000psi 1p 15000psi 2p 25000psi 1p 25000psi 2p D3,2(µm) 0.227 0.866 0.192 0.246 D4,3(µm) 0.263 0.441 0.208 0.279 Span 1.179 1.320 0.723 1.066 Cuadro 4.1: Valores de D3,2, D4,3 y span de las muestras en dependiendo del procesado. corresponde al mínimo en el tamaño de las gotas, sino que también supone una distribución más estrecha, con tamaños de gotas más similares entre sí. Efectuar varias pasadas tuvo un efecto contraproducente, parecido al de aplicar menores presiones, se destaca entonces la susceptibilidad por parte de estas nanoemulsiones al sobreprocesamiento. La combinación de ambos factores produjo la peor emulsión, dispersa y con gotas más cercanas al límite de macroemulsiones. Los valores de D4,3 , D3,2 y span también concuerdan con el comportamiento descrito (Cuadro 4.1), sirviendo de apoyo para dilucidar que una sola pasada a 25000 psi es el procedimiento óptimo para la síntesis de nanoemulsiones. A continuación, se comparó el estado de las nanoemulsiones según su proporción de tensioactivos en función del tiempo transcurrido. Comparando las distribuciones, resulta fácil deducir que las proporciones extremas producen peores resultados y duran menos, es más, en el caso de estos nanoemulgeles, se perdió el tamaño nanométrico antes de los 7 días. Las demás muestras obtuvieron resultados favorables, en especial, la proporción a partes iguales, aunque al final su distribución tendió a distribuciones similares al resto. Se presentan además los diámetros medios de Sauter, de volumen y el valor de span obtenidos en el primer día del nanoemulgel (Cuadro 4.2). Figura 4.2: Comparación de la DTG de las emulsiones en los días 1 y 3 desde su creación. 32 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales Figura 4.3: Comparación de la DTG de las emulsiones en los días 7 y 14 desde su creación. Día 1 100/0 75/25 50/50 25/75 0/100 D3,2(µm) 0.367 0.338 0.232 0.308 1.143 D3,2(µm) 0.691 0.383 0.275 3.980 9.789 Span 3.063 0.949 1.342 1.752 23.259 Cuadro 4.2: Valores de D3,2, D4,3 y de las muestras según su proporción en el primer día. Estos parámetros permiten diferenciar las cualidades de las nanoemulsiones más parecidas en la distribución de tamaños de gota, así la emulsión 75/25 a pesar de ser parecida a la 25/75, posee en realidad mejores características en relación a los diámetros medios de las gotas y por ello, probablemente presente mejor estabilidad también. Finalmente, se utiliza la DTG para establecer diferencias en la calidad y estabilidad de los nanoemulgeles generados con distintas proporciones de tensioactivos Tween 80/Span 80. Se compara por un lado la evolución del tamaño de las gotas con el paso del tiempo para una misma proporción, que comprende cuatro medidas en el intervalo de 7 días. En el caso de la muestra que sólo contenía Span 80, el tamaño ya estaba fuera de lo aceptable para nanoemulsiones casi desde el principio, viéndose menos afectada por el paso del tiempo. Por el otro, se incluyen, además, medidas realizadas en los pasos intermedios para analizar la influencia de cada uno en el tamaño de gota. En general, resultaron emulsiones monomodales, contando con un único tamaño de gota predominante, salvo en el caso de la emulsión que contaba únicamente con Tween 80. Este fenómeno se vio acrecentado por la agitación asociada a la mezcla del gel con la nanoemulsión, y, con el tiempo degeneró en una amplia distribución de tamaños incluso por encima de la micra, a pesar de tener un tamaño deseable al inicio. En todos los casos el paso del tiempo supuso un aumento tanto del tamaño de las gotas como de la amplitud de la distribución en la que se encuentran, indicando que el mecanismo de desestabilización predominante para las 33 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales Figura 4.4: DTG de la emusión 100/0 a lo largo del tiempo y en el procesado. Figura 4.5: DTG de la emusión 75/25 a lo largo del tiempo y en el procesado. Figura 4.6: DTG de la emusión 50/50 a lo largo del tiempo y en el procesado. 34 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales Figura 4.7: DTG de la emusión 25/75 a lo largo del tiempo y en el procesado. Figura 4.8: DTG de la emusión 0/100 a lo largo del tiempo y en el procesado. nanoemulsiones es la coalescencia. 4.2. Ensayos reológicos 4.2.1. Barridos de esfuerzo. Cálculo del esfuerzo crítico La finalidad de este ensayo es la determinación de la región lineal viscoelástica del nanoemulgel. Este parámetro es útil para conocer los esfuerzos a los que se puede someter un producto de estos materiales sin destrucción de microestructura. Este ensayo también proporciona otra información interesante sobre el material, en este caso, al representar los módulos de almacenamiento y pérdida, si tiene un comportamiento más próximo al de los sólidos o al de los líquidos. Un valor del módulo de pérdidas mayor que el de almacenamiento, indica una consistencia más fluida en todos los emulgeles desarrollados. Esto se explica por su alto contenido en agua, 35 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales pues la fase dispersa no es lo suficientemente concentrada para dar consistencia del sistema y la acción modificadora de la goma guar afecta principalmente a la textura en concentraciones semidiluidas. Se denomina esfuerzo crítico al valor de esfuerzo para el cual el módulo elástico G’ pierde su valor constante. Este valor indica el final de la LVR, y se determinó aproximando visulamente el final de la región plana en la representación del módulo elástico (Cuadro 4.3). A partir del esfuerzo crítico, la tendencia es decreciente. Figura 4.9: Barrido de esfuerzos para los nanoemulgeles 0/100 y 25/75 donde la flecha indica el esfuerzo crítico. Figura 4.10: Barrido de esfuerzos para los nanoemulgeles 50/50 y 75/25 donde la flecha indica el esfuerzo crítico. Los resultados apuntan a una reducción en el esfuerzo crítico conforme aumenta la cantidad de Tween 80. Dadas las cantidades con las que se experimenta, se destaca un comportamiento 36 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales Figura 4.11: Barrido de esfuerzos para el nanoemulgel 100/0 donde la flecha indica el esfuerzo crítico, seguida de la representación de los esfuerzos críticos en función de la proporción de Tween 80. 100/0 75/25 50/50 25/75 0/100 τc(Pa) 0.14 0.5 1.1 1.1 1.1 Cuadro 4.3: Valores aproximados del esfuerzo crítico para cada uno de los nanoemulgeles. considerablemnte sensible, y, por la importancia del parámetro, surge la necesidad de poder aproximarse a un valor de referencia a partir de su formulación. En una representación de los valores de críticos de esfuerzo frente al contenido en Tween 80, es posible apreciar una dependencia lineal entre ambas con la condición de que este representa al menos la mitad del tensioactivo, y dado el alto grado de correlación, se ha establecido una función que permita predecir su valor para un nanoemulgel de proporción diferente (Cuadro 4.4). A partir de la mitad en proporción, al disminuir el contenido en Tween 80, el valor aproximado del esfuerzo crítico se estanca. Valor Error M (Pa) 1.9953 0.1928 N(Pa) -0.0189 0.0025 R20.98 - Cuadro 4.4: Resultados de la pendiente My ordenada en el origen Nde la regresión lineal que relaciona la proporción en Tween 80 y el esfuerzo crítico. 37 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales 4.2.2. Barrido de frecuencias Partiendo del ensayo anterior, se estudia el comportamiento de los emulgeles ante los esfuerzos de cizalla oscilatoria dependiendo de la frecuencia, siempre dentro de la LVR. Figura 4.12: Barrido de frecuencias para los nanoemulgeles El comportamiento de los emulgeles resulta el propio de un líquido viscoelástico (Figura 2.15). El módulo de pérdidas se encuentra por encima del módulo de almacenamiento, siendo ambos crecientes con la frecuencia. La tendencia del módulo elástico es superar al módulo viscoso, en cuyo caso abandonaría el comportamiento de líquido viscoelástico para mostrar el propio de geles (Figura 2.15). El módulo elástico apenas oscila con la variación de concentraciones de tensioactivo, y aunque si existe cierta variación en el módulo viscoso, el comportamiento general es el mismo. Se induce así, que la reología del nanoemulgel viene dada principalmente por la concentración de goma guar, invariante en este trabajo. 4.2.3. Curvas de flujo. Ajuste al modelo de Cross Como ensayo reológico, se encargan de caracterizar la viscosidad del material en función de la velocidad de cizalla. Para los emulgeles preparados, se comprobó el cumplimiento del 38 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales modelo de Cross, asociado a materiales viscoelásticos con comportamiento pseudoplástico. Figura 4.13: Curvas de flujo de los nanoemulgeles 0/100 y 25/75, cuyo comportamiento no coincide con el modelo de Cross. La viscosidad de las muestras pobres en Tween 80 presentaban un comportamiento pseudoplástico, pero el modelo de Cross no se ajusta a su tendencia. A diferencia de una viscosidad con valores límite (Figura 2.16), su viscosidad tiene una tendencia permanentemente decreciente por lo que se puede observar en el intervalo del ensayo. Figura 4.14: Curvas de flujo de los nanoemulgeles 50/50 y 75/25, cuyo comportamiento sí se adapta al modelo de Cross. En cambio, en el resto de las muestras se comprobó un ajuste satisfactorio del modelo que, permite obtener información del material dado el significado físico de sus constantes (Cuadro 4.5). Además de la viscosidad inicial del sistema ante la ausencia de esfuerzos externos, se destaca el valor de la velocidad de cizalla crítica necesaria para producir el descenso de su viscosidad, de gran interés en el procesado, transporte y aplicación de posibles productos basados en los nanoemulgeles. También, se corrobora una vez más la predominancia del comportamiento líquido frente al sólido por un valor del índice de flujo más cercano a la unidad que a cero. 39 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales 50/50 75/25 100/0 Valor Error Valor Error Valor Error η0(Pa·s) 0.02122 5.09492·10−50,02037 5.66179·10−50.01411 3.43844·10−5 ˙ γ(s−1)42.39856 3.88628 44.21342 5.05597 70.43609 7.84077 n 0.69905 0.0152 0.68379 0.02091 0.66525 0.02254 R20.99879 - 0.99831 - 0.99862 - Cuadro 4.5: Resultados del ajuste del modelo de Cross al comportamiento de los nanoemulgeles 50/50, 75/25 y 100/0 con sus errores asociados. Figura 4.15: Curva de flujo del nanoemulgel 0/100 y comparativa de todas las curvas de flujo juntas. Al representar todas las curvas de flujo juntas se observa como el comportamiento viscoso es altamente dependiente de la proporción de Tween 80. Todas las curvas que satisfacen el modelo de Cross muestran una variación similar, encontrándose la principal diferencia en la viscosidad inicial. Esto concuerda con los barridos de esfuerzo, que en mayor proporción de Tween 80, mostraban más cercanía en el valor de los módulos viscoelásticos, indicando un comportamiento menos viscoso. 4.2.4. Medidas de estabilidad A partir del método de MLS es posible identificar los mecanismos de desestabilización de las nanoemulsiones. En cualquier proporción de los tensioactivos es posible observar es la separación gravitacional, que en este caso, se trata de cremado. Prueba de ello es la disminución en el tiempo de la intensidad de retrodispersión en la zona inferior del contenedor debido a la menor concentración de gotas de aceite esencial de mandarina. Este mecanismo se puede destacar en la desestabilización de las emulsiones unida a la coalescencia, la cual ya se apreciaba en la evolución de la DTG. Este fenómeno se observa en la disminución progresiva de la retrodispersión en la zona central, salvo para el nanoemulgel 75/25, en el que es apenas 40 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales 47 Sección 7 Referencias Referencias Anastas, P., y Warner, J. (1998). Green Chemistry: Theory and Practice. (First ed.). Oxford, Reino Unido: Oxford University. Barnes, H. A. (2000). Handbook of Elementary Rheology. Aberystwyth: University of Wales, Institute of Non-Newtonian Fluid Mechanics. Ciriminna R, Lomeli-Rodriguez M, Demma Carà P, Lopez-Sanchez JA, P. M. (2014). Chem- Comm. Chem. Commun(50), 15288–15296. Dosoky, N. S., y Setzer, W. N. (2018). Biological activities and safety of citrus spp. Essential oils. International Journal of Molecular Sciences,19(7), 1–25. FDA. (2017). Sec. 182.20 Essential oils, oleoresins (solvent-free), and natural extractives (including distillates). Grumezescu, A. M. (2016). Emulsions (Vol. 3). Elsevier Science Ltd. ICI Americas Inc. (1980). The HLB SYSTEM. Launay, B., Cuvelier, G., y Martinez-reyes, S. (1997). Viscosity of locust bean, guar and xanthan gum solutions in the Newtonian domain. Carbohydrate Polymers,34(4), 385–395. McClements, D. J. (1999). Food Emulsions: Principles, Practiques and Tecniques. Washington D.C.: CRC Press. Miller, R. (2016). Emulsifiers : Types and Uses (First ed.). Elsevier Science Ltd. Mudgil, D., Barak, S., y Khatkar, B. (2014, 06). Guar gum: processing, properties and food applications- a review. Journal of Food Science and Technology -Mysore-,51, 409-418. Nikfar, S., y Behboudi, A. F. (2014). Limonene. Encyclopedia of Toxicology: Third, 78–82. Oliveira, D., Michel, R. C., Mcbride, A. J. A., Moreira, A. S., y Lomba, R. F. T. (2013). Concentration Regimes of Biopolymers Xanthan , Tara , and Clairana , Comparing Dynamic Light Scattering and Distribution of Relaxation Time. PLoS ONE,8(5), 2–8. doi: 10.1371/journal.pone.0062713 Rosano, H. L., Gandolfo, F. G., y Hidrot, J. D. P. (1998). Stability of W1/O/W2 multiple emulsions influence of ripening and interfacial interactions. Colloids and Surfaces A: 48 Antúnez Domínguez, Jesús Manuel Doble Grado en Física e Ing. de Materiales Physicochemical and Engineering Aspects,138(1), 109–121. Salager, J.-l. (2016). Emulsion Phase Inversion Phenomena. En Emulsions and emulsion stability (Second ed., cap. Chapter 4). London, United Kingdom: CRC Press Taylor & Francis. Santos, J., Trujillo-cayado, L. A., y Calero, N. (2014). Physical Characterization of Eco-Friendly O / W Emulsions Developed Through a Strategy Based on Product Engineering Principles. American Institute of Chemical Engineers,60(7). Sjöblom, J. (2006). Emulsions and Emulsion Stability (Second ed.). Bergen (Noruega): CRC Press. Stang, M., Schuchmann, H., & Schubert, H. (2001). Emulsification in High-Pressure Homogenizers. Egineering in Life Sciences,1, 151–157. Tadros, T., Izquierdo, P., Esquena, J., y Solans, C. (2004). Formation and stability of nanoemulsions. Advances in Colloid and Interface Science,108-109, 303–318. Tadros, T. F. (2013). Emulsion Formation , Stability , and Rheology. En Emulsion formation and stability (pp. 1–75). Weinheim, Germany: Wiley-VCH. Taylor, P. (1995). Ostwald ripening in emulsions. Colloids and Surfaces A: Physicochemical and Engineering Aspects,99, 175–185. Trujillo-Cayado, L. A., Santos, J., Alfaro, M. C., Calero, N., y Muñoz, J. (2016). A further step in the development of oil-in-water emulsions formulated with a mixture of green solvents. Industrial and Engineering Chemistry Research,55(27), 7259–7266. Villalobos-Castillejos, F., Granillo-Guerrero, V. G., Leyva-Daniel, D. E., Alamilla-Beltrán, L., Gutiérrez-López, G. F., Monroy-Villagrana, A., y Jafari, S. M. (2018). Fabrication of Nanoemulsions by Microfluidization. Nanoemulsions: Formulation, Applications, and Characterization, 207–232. Vitale, S. A., y Katz, J. L. (2003). Liquid Droplet Dispersions Formed by Homogeneous Liquid - Liquid Nucleation : “ The Ouzo Effect ”. Langmuir,19, 4105–4110. 49