Isotermas dinámicas de sorción y análisis termodinámico de granos de quinua
Dynamic sorption isotherms and thermodynamic analysis of quinoa grains
Brislyn Osmayo Huacan1; Ciro Suca-Colana1; Jhony Mayta-Hancco1
1 Universidad Nacional de Moquegua - Escuela Profesional de Ingeniería Agroindustrial, Jr. Ancash s/n, Mariscal Nieto, Código Postal 1800, Moquegua, Peru.
ORCID de los autores:
C. Suca-Colana: https://orcid.org/0000-0003-3135-6715
J. Mayta-Hancco: https://orcid.org/0000-0002-5014-4681
RESUMEN
Las isotermas de sorción de dos variedades de granos de quinua (Chenopodium quinoa) fueron determinados usando los métodos DVS y DDI a tres temperaturas (20, 30 y 40 °C) en un rango de actividad de agua de 0,1 – 0,9. Las isotermas DDI fueron usadas para determinar la actividad de agua crítica aplicando de la segunda derivada Savistky-Golay. Siete modelos matemáticos fueron usados para ajustar los datos experimentales de sorción. Calor isostérico de sorción, entropía de sorción, energía libre de Gibbs y compensación entalpía-entropía fueron evaluadas. Los resultados indicaron que el contenido de humedad en equilibrio de los granos de quinua (PK y BJ) tienden a incrementar junto a la actividad de agua y la disminución de la temperatura. Las isotermas de los granos de quinua (PK y BJ) presentaron una forma sigmoidal tipo II. Los modelos Chung-Pfost y Smith fueron los mejores ajustados a los datos experimentales. La región de la actividad de agua crítica para los granos de quinua fue 0,75~0,85. El calor isostérico de sorción, entropía de sorción disminuyeron al incremento del contenido de humedad. Así mismo el incremento de contenido de humedad y temperatura influyeron en el decrecimiento de los valores de energía libre de Gibbs. La compensación entalpía-entropía fue corroborada (Tβ ≠ Thm), los procesos fueron llevados a cabo por la entalpía (Tβ > Thm) y son espontáneos (∆G > 0).
Palabras clave: Quinua, isotermas de sorción; actividad de agua crítica; propiedades termodinámicas; DDI; DVS.
ABSTRACT
The sorption isotherms of two varieties of quinoa grains (Chenopodium quinoa) were determined using the DVS and DDI methods at three temperatures (20, 30, and 40 °C) in a water activity range of 0,1–0,9. The DDI isotherms were used to determine the critical water activity by applying the Savistky-Golay second derivative. Seven mathematical models were used to fit the experimental sorption data. Isosteric heat of sorption, entropy of sorption, Gibbs free energy, and enthalpy-entropy compensation were evaluated. The results indicated that the equilibrium moisture content of quinoa grains (PK and BJ) tends to increase with water activity and decreasing temperature. The isotherms of quinoa grains (PK and BJ) presented a type II sigmoidal shape. The Chung-Pfost and Smith models were the best fit to the experimental data. The critical water activity region for quinoa grains was 0,75~0,85. The isothermal heat of sorption and sorption entropy decreased with increasing moisture content. Likewise, the increase in moisture content and temperature influenced the decrease in Gibbs free energy values. The enthalpy-entropy compensation was corroborated (Tβ ≠ Thm), the processes were carried out by enthalpy (Tβ > Thm) and are spontaneous (∆G > 0).
Keywords: Quinoa, sorption isotherms; critical water activity; thermodynamic properties; DDI; DVS.
1. Introducción
La quinua (Chenopodium quinoa) es originario de la región andina de Bolivia y Perú desde hace 5000 – 7000 años y distribuida a alrededor del mundo e incluye 250 especies (Vega-Gálvez et al., 2010). Dado que la semilla de la quinua no es un grano, sino un fruto, se le denomina pseudo – cereal e incluso pseudo – oleaginosa, esto debido a su inusual composición y a su excepcional equilibrio entre aceite, proteína y grasa (Tolaba et al., 2004). Los granos de quinua tienen un alto contenido de proteína (15-17,5%), rico en aminoácidos esenciales como la histidina y lisina, que no son encontrados en otros cereales (Miranda et al., 2011; Bustos-Vanegas et al., 2018). Los granos de quinua crecen en regiones entre 2500-4000 metros sobre el nivel del mar, causando que el grano contenga entre 9 y 10 % de humedad al momento de su cosecha (Vega-Gálvez et al., 2010; Pumacahua-Ramos et al., 2016). Sin embargo, las regiones donde la quinua es expendida cuentan con una humedad relativa alta, así mismo, la capa exterior de la quinua es altamente higroscópica, permitiendo que la humedad relativa del aire influya en el grano incrementando su humedad, disminuyendo rápidamente su estabilidad de almacenamiento (Arslan-Tontul, 2021).
El control del contenido de humedad durante el procesamiento de alimentos es un método antiguo de conservación, y se consigue eliminando el agua o aglutinándola de manera que el alimento se vuelva estable al deterioro microbiano y químico, siendo de mucha importancia conocer las características de sorción de humedad del alimento (Al-Muhtaseb et al., 2002). Un camino para presentar las características de sorción es mediante las isotermas de sorción de humedad, que presentan la relación entre el agua adsorbida o desorbida en el alimento y la actividad de agua a temperatura constante bajo condiciones de equilibrio (Perdomo et al., 2009).
Para determinar las isotermas de sorción de alimentos usan diversos métodos, como el método gravimétrico, sorción dinámica de vapor (DVS) e isoterma dinámica de punto de rocío (DDI). El método gravimétrico o tradicional, usa soluciones de sales saturadas, que presentan varios inconvenientes, como largos tiempos para alcanzar el equilibrio, dificultad para obtener datos precisos y grandes cantidades de muestra (Oduola et al., 2022). Estos inconvenientes fueron resueltos por el uso de analizadores dinámicos automatizados de sorción de vapor de agua (DVS y DDI), siendo el método DVS con más similitudes con el método convencional, ya que determina la humedad en equilibrio a actividades de agua conocidas (programadas 0,05-0,95), sin embargo el método DDI no alcanza el estado de equilibrio, ya que en este método somete a la muestra a un flujo de aire (aire húmedo o aire seco, dependiendo del proceso de sorción) durante un cierto intervalo de tiempo a continuación, el flujo se detiene y se registra la actividad del agua y del cambio de peso (contenido de agua), logrando obtener isotermas con mayor resolución (mayor cantidad de puntos) (Roos, 2016; Oduola el al., 2022).
Existen una gran cantidad modelos matemáticos que describen las isotermas de sorción en alimentos, como el modelo de GAB, BET, Smith, Oswin, entre otros, estos modelos son ecuaciones teóricas, semi empíricas y empíricas, y algunas de estas son dependientes de la temperatura. No existe un modelo matemático único que se ajuste a todos los alimentos, es por ello, que el modelo de sorción más apropiado puede seleccionarse teniendo en cuenta el grado de ajuste a los datos experimentales y la simplicidad del modelo (Arslan-Tontul, 2021).
La alta resolución del método DDI hace posible caracterizar cambios abruptos en las propiedades de sorción resultantes de los cambios inducidos por la humedad; la actividad de agua en el punto de la inflexión de la isoterma de sorción es llamada actividad de agua crítica, y es identificada así porque el alimento pierde su estabilidad (Carter et al., 2015).
Las propiedades termodinámicas, como el calor isostérico de sorción, entropía de sorción, energía libre de Gibbs y teoría de compensación entalpía-entropía pueden ser estimados a partir de las isotermas de sorción, y pueden ser usadas para el análisis y diseño de procesos alimentarios y estudiar la interacción humedad-sólido (Miranda et al., 2012).
Un estudio previo sobre las propiedades termodinámicas (calor isostérico de sorción, entropía de sorción y energía libre de Gibbs) a partir de isotermas de sorción de humedad en granos de quinua fue reportado por Arslan-Tontul (2021), usando el método gravimétrico para la construcción de isotermas, sin embargo, este método no permite determinar la actividad de agua crítica, como se logra usando métodos de sorción dinámicos (DDI), así mismo, mostrar la relación entre el calor isostérico y entropía de sorción (teoría isocinética) permite determinar si el proceso de sorción es un proceso ordenado (relacionado a la entalpía) o un proceso desordenado (relacionado a la entropía), y consecuentemente determinar si el proceso de sorción es espontáneo o no espontáneo.
El objetivo de este trabajo fue determinar y modelar isotermas de sorción de dos variedades de granos de quinua a diferentes temperaturas (20, 30 y 40 °C). Siete modelos matemáticos dependientes de la temperatura fueron selec-cionados para describir el comportamiento de sorción. Determinar las actividades de agua crítica e investigar las propiedades termodinámicas (calor isostérico de sorción, entropía de sorción, energía libre de Gibbs, compensación entalpía-entropía) del proceso de sorción.
2. Metodología
2.1. Obtención de muestra y tamizado
Los granos de quinua utilizada en este estudio fueron las variedades Blanca de Juli (BJ) y Pasankalla (PK) que fueron adquiridas de la empresa Agroindustrias Andinas Altiplano SRL de la ciudad de Juliaca, Perú. Se utilizó granos enteros y de tamaño grande según los requisitos de calidad de los granos de quinua INACAL (2022) mediante el uso de un tamiz 12 ASTM en un tamizador vibracional Analysette 18 (Fritsch GmbH, Alemania) durante un tiempo de tamizado de 5 min, con intervalos de 5 s y una amplitud de 2 mm, respectivamente. Previo al análisis, las muestras fueron almacenadas en frascos de vidrio herméticos a 4 ± 1 ºC, para el mantenimiento y uniformidad del contenido de humedad de los granos de quinua de ambas variedades.
2.2. Determinación de humedad y actividad de agua
El contenido de humedad se obtuvo por gravimetría utilizando un horno con circulación de aire forzado a 105 ± 1 ºC durante 24 h (Brasil, 2009). La actividad de agua (aw) se determinó utilizando un higrómetro (Aqualab, Decagon Series 3, USA) a 25 ± 1 °C y la calibración se realizó con un estándar de 8,5 mol L-1 de LiCl. En ambos casos los análisis se realizaron por triplicado.
2.3. Determinación de isotermas de sorción
Las isotermas de sorción de los granos de quinua fueron determinados por el método de punto de rocío utilizando un analizador de vapor de sorción (Aqualab VSA, 1122, Decagon, WA, USA) a 20, 30 y 40 °C. Se realizó una prueba estándar de isoterma usando una muestra de celulosa microcristalina proporcionada por Decagon (datos no mostrados). La calibración se realizó de acuerdo con las recomendaciones realizadas por el fabricante del equipo usando patrones desechables de aw: un patrón de 6,0 mol L-1 de NaCl, para una aw = 0,76 y; uno de 13,41 mol L-1 de LiCl, para una aw = 0,25 (ambos a 25 ± 1 °C). En cada prueba se utilizó cerca de 1500 mg de muestra, o lo suficiente para cubrir todo el fondo de la cápsula de acero inoxidable. Cada isoterma se generó en el intervalo de 0,10 ≥ aw ≤ 0,95 a temperatura ambiente.
Se utilizó el método DVS (del inglés dynamic vapor sorption), con paso de aw de 0,1, tolerancia de 0.1% de variación de masa (weight trigger) y sin tiempo de espera. Se utilizaron los puntos de equilibrio para la generación de isotermas, siendo los demás descartados.
También se utilizó el método DDI (del inglés dynamic dew point isotherm) para la generación de isotermas en el intervalo de aw = 0,1 hasta aw = 0,9, caudal de 100 mL min-1, resolución de 0,01 y sin tiempo de espera. Estos ajustes permitieron la generación de isotermas con resoluciones lo suficientemente altas como para permitir la detección de puntos de inflexión. Los puntos de inflexión agudos en la isoterma indican una actividad de agua crítica (awc) para la transición de fase. La awc de los granos de quinua a temperatura constante se determinó mediante estrategias de suavizado de curvas de la segunda derivada desarrollada por la empresa Meter Group, que se basa en el método de suavizado y diferenciación de Savitzky and Golay (1964). Las isotermas generadas por ambos métodos (DVS y DDI) fueron por duplicado.
2.4. Análisis de datos de sorción
Los datos experimentales del contenido de Xe/aw se usaron para construir curvas isotérmicas en una hoja de cálculo de Excel, con datos de aw y Xe ingresados los ejes x e y, respectivamente. Se utilizaron siete ecuaciones de isotermas recomendadas para granos y cereales, para ajustar los datos de Xe de los granos de quinua (Tabla 1). El ajuste se realizó mediante análisis de regresión no lineal con el programa Statistica v12 para Windows. El análisis de regresión no lineal se realizó con el algoritmo de Levenberg-Marquardt; este es un método de búsqueda para minimizar la suma de cuadrados de las diferencias entre los datos predichos y experimentales.
Tabla 1
Ecuaciones de Xe/aw utilizadas en el estudio
Modelo | Ecuación |
X_e=[(ln(1-a_w))/(-A(t+B))]^(1/C) | |
X_e=A-B ln[-(t+C)*ln(a_w)] | |
X_e=[-exp(A+Bt) /ln (a_w)]^(1/C) | |
X_e=(A+Bt)(a_w/(1-a_w)) ^(1/C) | |
X_e=(AB(C/t) a_w)/((1-Ba_w) [1-Ba_w+B(C/t)] a_w) | |
X_e=exp(A-Bt)/(C-ln(a_w)) | |
Smith | X_e=A-Bt-C*ln(1-a_w) |
Donde Xe es el contenido de humedad de equilibrio (%, base seca); aw es la actividad de agua (decimal); t es temperatura (°C) y A, B y C son constantes empíricas específicas de las ecuaciones.
Los criterios utilizados para determinar la mejor ecuación para los datos Xe/aw fueron el coeficiente de determinación (R2), error relativo medio (MRE) (Ecuación 1) y el error estándar estimado (SEE) (Ecuación 2). Los valores de MRE ≤ 10% es indicativo de que el modelo es bueno en la predicción de las humedades de equilibrio (Aguerre et al., 1989). Este criterio fue utilizado para la selección del modelo adecuado en la determinación de las propiedades termodiná-micas de los granos de quinua. Los valores mínimos de SEE y el valor máximo de R2 reportan el mejor modelo.
Por otro lado, las diferencias entre el valor pronosticado por la ecuación dada y el valor experimental del Xe se representaron para producir gráficos de los residuales. Estos gráficos se evaluaron por su tendencia aleatoria (A) o tendenciosa (T).
|
| (1) |
|
| (2) |
Donde Xexp es el valor experimental del contenido de humedad en equilibrio (%, b.s.), Xpre es el valor predicho del contenido de humedad en equilibrio a través del modelo (%, b.s.), N el número de mediciones experimentales y DF grados de libertad de los residuales (número de datos observados menos número de parámetros del modelo).
2.5. Propiedades termodinámicas
Las propiedades termodinámicas de los granos de quinua fueron obtenidos por el método descrito por Corrêa et al. (2012) como se expresan en las ecuaciones 3 al 7.
| (3) |
| (4) |
| (5) |
| (6) |
| (7) |
|
|
| (8) |
Donde Ru es la constante universal de los gases (8,314 J mol‒1 K‒1), Ta es temperatura absoluta (K), Δhst es calor isostérico de sorción líquido (kJ kg‒1), L es calor latente de vaporización del agua libre (kJ kg‒1), Qst es calor isostérico total (kJ kg‒1), ΔS es entropía diferencial de sorción (kJ kg‒1 K‒1), ΔG es energía libre de Gibbs (kJ kg‒1), TB es temperatura isocinética (K), Thm es temperatura armónica media (K), n es número de temperaturas utilizadas y, m es número de pares de datos de entalpía y entropía.
3. Resultados y discusión
3.1. Contenido de humedad y actividad de agua
El resultado del análisis del contenido de hume-dad y actividad de agua de las dos variedades de quinua se muestran en la Tabla 2. Ambas variedades de quinua presentaron valores similares de humedad y actividad de agua. Los valores de humedad fueron similares a los obtenidos por Miranda et al. (2012), Pumacahua-Ramos et al. (2016) y Ramos Diaz et al. (2017) en diferentes muestras de granos de quinua. Por otro lado, nuestros resultados de humedad de las dos variedades de quinua se ubican por debajo de las recomendaciones de humedad (13%) establecido en la norma técnica peruana NTP 205.062:2021 para granos andinos (INACAL, 2022) y norma para la quinua del Codex Alimentarius (CXS 333-2019, 2020).
Tabla 2
Contenido de humedad y actividad de agua de dos variedades de quinua entera
| Variedad1 | |
BJ | PK | |
Contenido de humedad (g/100 g) | 11,03a ± 0,62 | 10,25a ± 1,08 |
Actividad de agua (adimensional) | 0,39a ± 0,01 | 0,38a ± 0,01 |
Los valores representan promedios ± desviación estándar. Distintos superíndices en la misma fila indican diferencias entre variedades de quinua (p < 0,05). Prueba de rango múltiple de Duncan.
1 BJ: quinua variedad Blanca de Juli; PK: quinua variedad Pasankalla.
Los valores de actividad de agua (aw) de los granos de quinua en estudio variaron de 0,38 a 0,39. Estos valores se ubican por debajo de 0,7 y no favorecen el crecimiento de hongos durante el almacenamiento de los granos y sus productos (Abdullah et al., 2000; Aldred et al., 2004; Rosentrater, 2022).
3.2. Isotermas de sorción de agua
Los datos experimentales del contenido de humedad en equilibrio (Xe) vs actividad de agua (aw) de las dos variedades de quinua (BJ y PK) a las temperaturas de trabajo (20, 30 y 40 °C) se muestran en la Figura 1. En general, se observó una correlación positiva entre Xe y aw. para las dos variedades de quinua. Arslan-Tontul (2021) también reportó incrementos de Xe y aw en el proceso de adsorción de agua de granos de quinua blanca y atribuido a la alta higroscopicidad de los granos de quinua y el aumento de la presión de vapor circundante en el alimento.
Por otro lado, la Xe de ambas variedades de quinua se correlacionan negativa con el incremento de la temperatura. Al respecto, el incremento de temperaturas puede disminuir la capacidad de adsorción de agua en los granos debido al aumento de la energía cinética asociada a las moléculas de agua presente en los alimentos que a su vez disminuye de las fuerzas de atracción entre las moléculas de agua y en consecuencia su evaporación (Arslan-Tontul, 2021; Miranda et al., 2012). Varios autores (Arslan-Tontul, 2021; Bustos-Vanegas et al., 2018; Miranda et al., 2012; Pumacahua-Ramos et al., 2016) observaron estas correlaciones positivas y negativas en las muestras de granos de quinua. Asimismo, el tipo de isoterma corresponder a las isotermas tipo II según la clasificación de Brunauer et al. (1938) y es característica en cereales, granos y semillas (Miranda et al., 2012).
Los valores del Xe obtenidos en el proceso de desorción siempre fueron superiores a los valores obtenidos por adsorción en las diferentes condiciones de temperatura y aw, evidenciando el fenómeno de histéresis (Figura 1). Estos resultados se corresponden a los encontrados en semillas de quinua (Bustos-Vanegas et al., 2018; Miranda et al., 2012). La diferencia entre la desorción y adsorción es conocida como histéresis y es un proceso irreversible (Wolf et al., 1972) que permite evaluar el deterioro químico, microbiológico y la estabilidad del producto durante el almacenamiento (da Silva et al., 2018). Lahsasni et al. (2004) explican el fenómeno de histéresis considerando un poro de estructura rígida conectado a su entorno por un pequeño capilar, durante la adsorción, el capilar comienza a llenarse como resultado del aumento de la humedad relativa. Cuando la presión parcial del vapor en el aire se vuelve mayor que la presión del vapor del líquido en el capilar, la humedad se moverá hacia el poro. En la desorción, el poro está inicialmente lleno de líquido en la saturación. Este líquido puede evaporarse solo cuando la presión del aire circundante se vuelve más baja que la presión de vapor del líquido dentro del capilar. Adicionalmente, los mismos autores afirman que las diferencias entre adsorción y desorción se dan debido a la gran variedad de diámetros capilares en el sistema poroso (Bustos-Vanegas et al., 2018; Lahsasni et al., 2004).

Figura 1. Contenido de humedad de equilibrio para granos de quinua variedad Blanca de Juli (BJ) y Pasankalla (PK) a diferentes temperaturas. Las líneas (---) corresponden a ajustes de los modelos de Chung-Pfost modificado y (―) Smith.
3.3. Evaluación de la capacidad de ajuste de los modelos de sorción
Los valores de los parámetros de los siete modelos de equilibrio higroscópico ajustados a los datos experimentales se muestran en la Tabla 3. En el presente estudio, todos los modelos de sorción para los granos de quinua BJ y PK utilizados mostraron una buena calidad de ajuste (R2 > 95%), excepto el modelo modificado de Henderson, GAB y Harkins-Jura, para la quinua PK. Se obtuvieron valores de R2 > 99,3% para el proceso de adsorción de los granos de quinua BJ y PK con el modelo de Smith, mientras que, en el proceso de desorción, los modelos de Chung-Pfost modificado y Smith son los que presentan mejor ajuste.
Según el error relativo medio (MRE) en el proceso de adsorción de la quinua BJ (Tabla 3), los modelos modificados de Chung-Pfost y Oswin presentaron valores de ~ 9%; mientras que con el modelo de Smith se obtuvo 3,3%. En el proceso de desorción, los modelos modificados de Henderson, Oswin y GAB presentaron valores de ~ 8 a ~ 9%, mientras que el modelo modificado de Chung-Pfost se obtuvo un 3,7%. En el proceso de adsorción de los granos de quinua PK, todos los modelos presentaron valores de MRE > 10%, excepto el modelo de Smith que obtuvo 4,1%. Sin embargo, en el proceso de desorción, los valores más bajos de MRE se obtuvieron con el modelo modificado Chung-Pfost y Smith con 7,0% y 6,7%, respectivamente.
Tabla 3
Parámetros del modelo matemático ajustado al contenido de humedad de equilibrio de granos de quinua variedad Blanca de Juli (BJ), Pasankalla (PK) y sus respectivos coeficientes de determinación (R2), error relativo medio (MRE), error estándar estimado (SEE) y distribución residual (RD)
Proceso | Pará- metros | Modelo | ||||||
Modificado de Henderson | Modificado de Chung-Pfost | Modificado de Halsey | Modificado de Oswin | Modificado de GAB | Harkins-Jura | Smith | ||
BJ | ||||||||
Adsorción | A B C R2 (%) MRE (%) SEE RD | 0,00025 9,84 1,77 96,28 12,23 1,56 T | 132,82 ‒0,89 0,17 98,06 9,48 1,13 T | 5,68 ‒0,033 2,16 98,61 10,09 0,96 A | 15,80 ‒0,17 2,77 98,35 9,32 1,05 A | 7,11 0,79 466,61 95,74 14,30 1,67 T | 2,92 0,015 0,41 97,23 10,46 1,35 T | 11,10 0,21 8,45 99,68 3,27 0,46 A |
Desorción | A B C R2 (%) MRE (%) SEE RD | 0,00013 3,03 2,04 97,85 9,05 1,18 A | 145,70 ‒3,07 0,17 99,40 3,69 0,63 A | 6,77 ‒0,04 2,44 97,69 11,34 1,15 T | 17,88 ‒0,19 3,14 98,78 8,07 0,89 T | 12,54 0,56 203,02 95,98 11,59 1,61 T | 3,22 0,02 0,58 97,84 8,77 1,18 T | 12,89 0,24 8,23 99,14 6,97 0,75 A |
PK | ||||||||
Adsorción | A B C R2 (%) MRE (%) SEE RD | 0,00018 ‒1,98 2,00 93,55 16,07 2,08 T | 115,05 ‒8,46 0,18 97,03 10,95 1,42 T | 6,76 ‒0,05 2,37 97,39 12,48 1,33 T | 18,72 ‒0,24 3,08 96,60 12,82 1,52 T | 7,57 0,77 594,83 90,69 21,33 2,48 T | 3,23 0,02 0,52 94,62 14,31 1,90 T | 14,56 0,30 8,09 99,36 4,13 0,66 A |
Desorción | A B C R2 (%) MRE (%) SEE RD | 0,00011 ‒6,08 2,21 96,01 14,98 1,69 A | 116,47 ‒10,01 0,17 98,85 7,09 0,91 A | 7,87 ‒0,06 2,62 96,99 14,79 1,47 T | 21,16 ‒0,28 3,40 97,58 12,28 1,32 T | 26,70 0,23 183,74 92,90 15,96 2,23 T | 3,52 0,02 0,68 96,38 13,40 1,61 T | 16,67 0,34 8,09 99,27 6,67 0,73 A |
Donde: A, aleatoria; T, tendenciosa.
Todos los modelos de sorción de agua de los granos de quinua BJ y PK presentaron valores del error estándar estimado (SEE) menores que 2 (Tabla 3), excepto los modelos modificados de Henderson y GAB para la quinua PK. En el proceso de adsorción se obtuvieron valores de 0,5 y 0,7 con el modelo de Smith para de la quinua BJ y PK, respectivamente. Sin embargo, en el proceso de desorción de la quinua BJ y PK se obtuvo valores de SEE de 0,6 y 0,7 con los modelos de Chung-Pfost modificado y Smith, respectivamente. Se observa que los modelos de Smith y el modificado de Chung-Pfost, fueron los más adecuados para describir el fenómeno de adsorción y desorción para la quinua BJ, respectivamente. Por otro lado, el modelo de Smith fue el que mejor describe los datos del Xe de la quinua PK. En general, los parámetros estadísticos están de acuerdo con lo señalado por varios autores (Aguerre et al., 1989; Bustos-Vanegas et al., 2018; Zeymer et al., 2022), tales como un alto coeficiente de determinación (R2), bajo error estándar estimado (SEE), error relativo medio (MRE) menor al 10% y distribución residual reducida y aleatoria. Bustos-Vanegas et al. (2018) reportaron al modelo modificado de Halsey como el más adecuado para describir el fenómeno de adsorción y desorción de semillas de quinua a temperaturas de 15, 25, 35 y 50 °C utilizando el método estático-gravimétrico. Estas diferencias en comparación a nuestros resultados podrían ser debido al método utilizado para la determinación del Xe en los granos de quinua, puesto que el método de punto de rocío fue utilizado en la presente investigación. Por otro lado, existen reportes del uso de modelos teóricos para describir el proceso de adsorción y desorción de agua de granos de quinua utilizando los modelos de Peleg (Pumacahua-Ramos et al., 2016), GAB de tres parámetros y Halsey (Arslan-Tontul, 2021), Iglesias-Chirife y Oswin (Miranda et al., 2012), GAB de cuatro parámetros (Rohvein et al., 2004) y GAB de tres parámetros (Tolaba et al., 2004). Estos reportes son diferentes a nuestros resultados debido a que los modelos matemáticos descritos en la Tabla 1 contienen tres parámetros, incluyendo la temperatura y son consideradas como ecuaciones estándar para cereales y oleaginosas por la Sociedad Americana de Ingenieros Agrícolas (Yu et al., 1999). Otros factores que pueden influir en la selección del modelo matemático puede ser la variedad, el proceso de lavado previo al estudio y la composición química de la quinua (Rohvein et al., 2004; Tolaba et al., 2004).
3.4. Actividad de agua crítica
La actividad de agua asociada con el punto de inflexión en la curva DDI es caracterizado como la actividad de agua crítica a la temperatura de la isoterma, este cambio abrupto, indica una activi-dad de agua crítica por cambios en las propie-dades de difusión de la humedad del alimento (Carter et al., 2016).
La Figura 2 muestra las curvas DDI de los granos de quinua a 20, 30 y 40 °C. Se observa un punto de inflexión entre 0,75~0,85 en ambas varieda-des, siendo verificado mediante la aplicación de la segunda derivada de Savistky-Golay. Un comportamiento similar fue encontrado por Carter et al. (2016) en granos de trigo blanco, indicando que el punto de inflexión indica un cambio repentino en las propiedades de sorción asociados con la consecución de una actividad de agua específica en la que la pendiente se acentúa y se absorbe más humedad a medida que aumenta la actividad de agua, este fenómeno sería similar a una transición vítrea en la que la movilidad y el volumen libre aumenta repentinamente, haciendo que los materiales polares como el azúcar y el almidón estén disponibles para el agua que en el estado previo de baja movilidad.
3.5. Calor isostérico de sorción (Qst)
El Qs es un parámetro termodinámico calculado a partir del Xe a diferentes temperaturas (Arslan-Tontul, 2021) y se define como el calor total de sorción menos el calor latente de vaporización del agua a una temperatura dada y usada comúnmente para estudiar la humedad y diseñar los procesos de secado de alimentos (Miranda et al., 2012; Tolaba et al., 2004).

Figura 2. Isotermas dinámicas de desorción para granos de quinua variedad Blanca de Juli (BJ) y Pasankalla (PK) a 20, 30 y 40 °C. Los gráficos de la segunda derivada de Savistky-Golay también se muestran en el eje secundario y debajo de la curva original. El valor máximo en la curva de la 2da derivada identifica el punto de inflexión en la curva original.
En la Figura 3 se observa una clara tendencia decreciente del Qst con el incremento de Xe para las dos variedades de quinua. Este fenómeno fue similar a los obtenidos por Pumacahua-Ramos et al. (2016), Bustos-Vanegas et al. (2018), Miranda et al. (2012) y Arslan-Tontul (2021) en granos de quinua y lo atribuyeron a la disminución de la interacción agua-sólido de los sitios polares en la superficie del producto que están cubiertos por moléculas de agua formando la monocapa, mientras el contenido de humedad aumenta.
Los valores de Qst en el proceso de adsorción de la quinua BJ fue mayor que la desorción, mientras que en la quinua PK los valores tanto de adsorción como desorción fueron menores a ~ 10000 kJ kg-1 (Figura 3). En general, los valores de Qst para las dos variedades de quinua se encuentran entre 5 y 18% de humedad. Cuando los niveles de humedad fueron superiores a ~18%, los calores de sorción fueron similares para ambas variedades de quinua, incluido el calor latente de vaporización de 2430 kJ kg-1 a una temperatura media de 30 °C.

Figura 3. Efecto del contenido de humedad de equilibrio en el calor isostérico de sorción de granos de quinua variedad Blanca de Juli (BJ) y Pasankalla (PK). Las líneas (---) corresponden a los valores de calor latente de vaporización (HL).
3.6. Entropía de sorción
La ΔS es proporcional al número de sitios de sorción disponibles correspondientes a un nivel de energía específico, que describe el grado de desorden y aleatoriedad (disposición espacial) de movimiento de las moléculas de agua (Cano-Higuita et al., 2015). La entropía, por tanto, caracteriza o define el grado de orden o aleatoriedad existente en el sistema agua-sorbente, lo que puede ayudar a interpretar procesos como la disolución, la cristalización y el hinchamiento (Hassini et al., 2015).
La Figura 4 muestra la dependencia de la entropía de sorción al contenido de humedad, se observó que ΔS muestra valores máximos en valore más bajos de Xe, donde se produce la sorción de la monocapa.

Figura 4. Efecto del contenido de humedad de equilibrio en la entropía de sorción de granos de quinua variedad Blanca de Juli (BJ) y Pasankalla (PK).
Comportamientos similares fueron mostrados por Bustos-Vanegas et al. (2018) y Arrslan-Tontul (2020), debido a que en un el alto contenido de humedad, las moléculas de agua presentes en la superficie disminuyen la capacidad de unión de más moléculas a los sitios activos de macro-moléculas. Así, la tasa de sorción y la entropía disminuyen. Por el contrario, los sitios de unión activos disponibles son más numerosos con un contenido de humedad más bajo, lo que aumenta la entropía.
Los valores de ΔS en el proceso de adsorción de la quinua BJ fueron mayores al proceso de desorción en humedades menores a ~8%, muy por el contrario, la quinua PK, mostró valores altos en el proceso de desorción que en el proceso de adsorción ~25 kJ kg-1 K-1 en humedades inferiores a ~15%. Los valores de ΔS fueron similares a humedades mayores a ~20% en ambas variedades.
3.7. Energía libre de Gibbs (ΔG)
La ΔG es una función termodinámica que representa la cantidad de energía máxima liberada en un proceso a una temperatura y presión constantes (equilibrio entre entalpía y entropía), que está disponible para ser usada (Bustos-Vanegas et al., 2018).
La variación de valores de ΔG con el contenido de humedad y la temperatura son mostrados en la Figura 5. La figura indica que la energía libre disminuyó desde ~30 kJ kg-1 K-1 a valores cercanos a cero con el incremento del contenido de humedad en ambas variedades de quinua. Resultados similares fueron obtenidos Bustos-Vanegas et al. (2018), y Arslan-Tontul (2021), y atribuyendo este resultado a los alimentos hidrófilos que presentan más niveles de energía debido a la alta afinidad de adsorción de agua. Generalmente, si ΔG es negativo, el proceso de sorción se considera espontáneo (el proceso ocurre por sí mismo), por el contrario, el proceso es considerado como no espóntaneo (Zhang et al., 2022). La Figura 5 muestra que los procesos de adsorción y desorción en ambas variedades de quinua son espontáneos.

Figura 5. Energía libre de Gibbs del proceso de sorción de granos de quinua.
3.8. Teoría de compensación entalpía-entropía
La teoría de compensación entalpía-entropía o teoría isocinética es usada para evaluar fenó-menos físicos y químicos como las reacciones de sorción y proponer una relación lineal entre ΔS y Δh, y permite comprobar si habrá mayor intera-cción molecular debido a la reducción de la libertad o al enlace de las moléculas del alimento, generando mayor organización u orden (relacio-nado con la entalpía) frente a la desorganización y mayor libertad de las moléculas del alimento (relacionado con la entropía) (Spada et al., 2013). La teoría solo es aplicable si la temperatura iso-cinética (Tβ) es diferente de la temperatura media harmónica (Thm) (Pavón-García et al., 2015).
La Figura 6 muestra la relación entalpía-entropía con el proceso de sorción de granos de quinua, mostrando una fuerte correlación entre las dos propiedades termodinámicas, las Tβ fueron desiguales a Thm en ambas variedades de quinua (Tabla 4). En todos los casos Tβ fue mayor que Thm, indicando que los procesos de sorción fueron llevados a cabo por la entalpía, indicando ser un proceso ordenado, así mismo los valores positivos de ΔG, confirma que el proceso de sorción es un proceso espontáneo. Estudios previos, muestran resultados similares para la compensación ental-pía-entropía obtenidos por Bustos-Vanegas et al. (2018), Hassini et al. (2015) y Corrêa et al. (2015).
4. Conclusiones
La actividad de agua y contenido de humedad de los granos de quinua BJ y PK fueron de 0,39 y 11,03%, 0,38 y 10,25% respectivamente, cuyos valores evitan el crecimiento microbiano.
Los datos de sorción para los granos de quinua BJ fueron ajustados al modelo de Smith y Chung-Pfost en los procesos de adsorción y desorción respectivamente, a diferencia de los datos granos de quinua PK que fueron ajustados al modelo de Smith en ambos procesos y las isotermas obtenidas fueron de tipo II.
La actividad de agua crítica de los granos de quinua fluctuó entre 0,75~0,85, identificando que la permeabilidad de la humedad de la superficie al grano cambia significativamente cuando la actividad de agua del grano aumenta hasta valores superiores a de actividad de agua crítica. Los valores de entalpía, entropía y energía libre de Gibbs decrecieron con el incremento del contenido de humedad en ambos granos de quinua. La teoría de compensación entalpía-entropía fue confirmado y se puede concluir que el mecanismo de sorción es controlado por la entalpía para ambas variedades de quinua.
La determinación de las propiedades termodiná-micas y actividad de agua crítica en granos de quinua abre la puerta a posteriores investiga-ciones sobre el modelado de estabilidad y vida útil, diseño de procesos de almacenamiento defi-niendo empaques que permitan una mayor esta-bilidad del producto, relación entre las propie-dades termodinámicas y composición química.

Figura 6. Compensación entalpía-entropía del proceso de sorción de los granos de quinua.
Tabla 4
Teoría de compensación entalpía-entropía de granos de quinua
Proceso | Thm (K) | Tβ (K) | ∆G (kJ/kg) | Tipo de proceso | Control del proceso | |
Blanca de Juli | Adsorción | 302,93 | 315,95 | 44,55 | Espontaneo | Entalpía |
Desorción | 302,93 | 317,35 | 32,24 | Espontaneo | Entalpía | |
Pasankalla | Adsorción | 302,93 | 326,32 | 61,16 | Espontaneo | Entalpía |
Desorción | 302,93 | 335,97 | 30,98 | Espontaneo | Entalpía |
Agradecimientos
La presente investigación forma parte del proyecto “Diagramas de estado de semillas de quinua”.
Referencias bibliográficas
Abdullah, N., Nawawi, A., & Othman, I. (2000). Fungal spoilage of starch-based foods in relation to its water activity (aw). Journal of Stored Products Research, 36(1), 47-54. https://doi.org/10.1016/S0022-474X(99)00026-0
Aguerre, R., Suarez, C., & Viollaz, P. E. (1989). New BET type multilayer sorption isotherms. Part II: Modeling water sorption in foods. Lebensmittel-Wissenschaft und-Technologie, 22(4), 192–195.
Aldred, D., Magan, N., & Olsen, M. (2004). 7 - The use of HACCP in the control of mycotoxins: the case of cereals. In N. Magan & M. Olsen (Eds.), Mycotoxins in Food (pp. 139-173). Woodhead Publishing. https://doi.org/10.1533/9781855739086.2.139
Al-Muhtaseb, A. H., McMinn, W. A. M., & Magee, T. R. A. (2002). Moisture Sorption Isotherm Characteristics of Food Products: A Review. Food and Bioproducts Processing, 80(2), 118-128. https://doi.org/10.1205/09603080252938753
Arslan-Tontul, S. (2021). Moisture sorption isotherm and thermodynamic analysis of quinoa grains. Heat and Mass Transfer, 57(3), 543-550. https://doi.org/10.1007/s00231-020-02978-8
Brasil. (2009). Regras para análise de sementes (1a ed.). Ministério da Agricultura, Pecuária e Abastecimento. Brasília, DF.
Brunauer, S., Emmett, P. H., & Teller, E. (1938). Adsorption of Gases in Multimolecular Layers. Journal of the American Chemical Society, 60(2), 309-319. https://doi.org/10.1021/ja01269a023
Bustos-Vanegas, J. D., Corrêa, P. C., Zeymer, J. S., Baptestini, F. M., & Campos, R. C. (2018). Moisture Sorption Isotherms of Quinoa Seeds: Thermodynamic Analysis. Engenharia Agrícola, 38(6), 941-950. https://doi.org/10.1590/1809-4430
Cano-Higuita, D., Villa-Vélez, M., Harvey, A., Telis-Romero, J., Váquiro, H., Alexander, T., & Nicoletti, V. (2015). Influence of alternative drying aids on water sorption of spray dried mango mix powders: A thermodynamic approach. Food and Bioproducts Processing, 93, 19-28. https://doi.org/10.1016/j.fbp.2013.10.005.
Carter, B. P., Galloway, M. T., Campbell, G. S. et al. (2015). The critical water activity from dynamic dewpoint isotherms as an indicator of crispness in low moisture cookies. Food Measure 9, 463–470. https://doi.org/10.1007/s11694-015-9254-3.
Carter, Brady, P., Galloway, Mary, T., Campbel, Gaylon, S., & Carter, Arron, H. (2016). Changes in the Moisture Permeability of Grain at the Critical Water Activity from Dynamic Dewpoint Isotherms. American Society of Agricultural and Biological Engineers. 59(3), 1023-1028. doi: 10.13031/trans.59.11379.
Chen, C.-C., & Morey, R. V. (1989). Comparison of Four EMC/ERH Equations. Transactions of the ASAE, 32(3), 983-0990. https://doi.org/10.13031/2013.31103
Corrêa, P. C., Oliveira, G. H. H., & Santos, H. S. (2012). Thermody-namic Properties of Agricultural Products Processes. In I. Arana (Ed.), Physical Properties of Foods: Novel Measurement Techniques and Applications (1st ed., pp. 131–142). CRC Press. https://doi.org/10.1201/b11542
CXS 333-2019. (2020). Standard for Quinoa. In: Codex Alimentarius. Internacional Food Standards.
da Silva, H. W., Rodovalho, R. S., & Silva, I. L. (2018). Hysteresis and thermodynamic properties of water sorption in ‘Malagueta’ pepper seeds. Revista Brasileira de Engenharia Agrícola e Ambiental, 22(9), 658-663. https://doi.org/10.1590/1807-1929/agriambi.v22n9p658-663
Harkins, W. D., & Jura, G. (1944). Surfaces of Solids. XIII. A Vapor Adsorption Method for the Determination of the Area of a Solid without the Assumption of a Molecular Area, and the Areas Occupied by Nitrogen and Other Molecules on the Surface of a Solid. Journal of the American Chemical Society, 66(8), 1366-1373. https://doi.org/10.1021/ja01236a048
Hassini, L., Bettaieb, E., Desmorieux, H., Torres Sandoval, S., Touil, A. (2015). Desorption isotherms and thermodynamic properties of prickly pear seeds. Industrial Crops and Products, 67, 457-465. https://doi.org/10.1016/j.indcrop.2015.01.078.
Iglesias, H. A., & Chirife, J. (1976). Prediction of the effect of temperature on water sorption isotherms of food material. International Journal of Food Science & Technology, 11(2), 109-116. https://doi.org/10.1111/j.1365-2621.1976.tb00707.x
INACAL. (2022). Granos Andinos. Quinua en grano. Requisitos. In Normas Técnicas Peruanas. NTP 205.062:2021 (Vol. 3): Instituto Nacional de Calidad (INACAL). Perú.
Jayas, D. S., & Mazza, G. (1993). Comparison of Five, Three-parameter Equations for the Description of Adsorption Data of Oats. Transactions of the ASAE, 36(1), 119-125. https://doi.org/10.13031/2013.28322
Lahsasni, S., Kouhila, M., & Mahrouz, M. (2004). Adsorption–desorption isotherms and heat of sorption of prickly pear fruit (Opuntia ficus indica). Energy Conversion and Management, 45(2), 249-261. https://doi.org/10.1016/S0196-8904(03)00133-X
Miranda, M., Vega-Gálvez, A., Sanders, M., López, J., Lemus-Mondaca, R., Martínez, E., & Di Scala, K. (2012). Modelling the Water Sorption Isotherms of Quinoa Seeds (Chenopodium quinoa Willd.) and Determination of Sorption Heats. Food and Bioprocess Technology, 5(5), 1686-1693. https://doi.org/10.1007/s11947-011-0610-y
Oduola, A., Luthra, K., & Atungulu, G. (2022). Determination of moisture sorption isotherms of industrial hemp (Cannabis sativa L.) flower and leaf composite powders. Industrial Crops and Products, 168, 115201. https://doi.org/10.1016/j.indcrop.2022.115201
Pavon-García, L. M. A., Gallardo-Rivera, R., Román-Guerrero, A., Carrillo-Navas, H., Rodríguez-Huezo, M. E., Guadarrama-Lezama, A. Y., & Pérez-Alonso, C. (2015). Moisture sorption properties and storage stability conditions. Revista Mexicana de Ingeniería Química, 14(3), 601-613.
Perdomo, J., Cova, A., Sandoval, A. J., Laredo, E., Müller, A. J., & García Rodríguez, H. L. (2009). Glass transition temperatures and water sorption isotherms of cassava starch. Carbohydrate Polymers, 76(2), 305-313. https://doi.org/10.1016/j.carbpol.2008.10.023
Pfost, H. B., Maurer, S. G., Chung, D. S., & Milliken, G. A. (1976). Summarizing and reporting equilibrium moisture data for grains. ASAE(ASAE paper), 76–3520.
Pumacahua-Ramos, A., Gomez Vieira, J. A., Telis- Romero, J., Villa-Vélez, H. A., & Lopes Filho, J. F. (2016). Isotherms and isosteric heat of sorption of two varieties of Peruvian quinoa. Scientia Agropecuaria, 7(4), 409-417. https://doi.org/10.17268/sci.agropecu.2016.04.06
Ramos Diaz, J. M., Sundarrajan, L., Kariluoto, S., Lampi, A.-M., Tenitz, S., & Jouppila, K. (2017). Effect of Extrusion Cooking on Physical Properties and Chemical Composition of Corn-Based Snacks Containing Amaranth and Quinoa: Application of Partial Least Squares Regression. Journal of Food Process Engineering, 40(1), e12320. https://doi.org/10.1111/jfpe.12320
Rohvein, C., Santalla, E., & Gely, M. C. (2004). Estimation of sorption isotherm and the heat of sorption of quinoa (Chenopodium quinoa Willd.) seeds. Food Science and Technology International, 10(6), 409-413. https://doi.org/10.1177/1082013204049513
Rosentrater, K. A. (2022). Chapter 17 - Biochemical, functional, and nutritive changes during storage. In K. A. Rosentrater (Ed.), Storage of Cereal Grains and Their Products (Fifth Edition) (pp. 443-501). Woodhead Publishing. https://doi.org/10.1016/B978-0-12-812758-2.00010-6
Roos, Y. H., & Bambang, N. (2016). Dynamic water sorption for the study of amorphous content of vacuum-dried honey powder. Powder Technology, 301, 981-988. https://doi.org/10.1016/j.powtec.2016.07.055.
Spada, C., Zapata, J., Caciano, N., Ferreira, P., Ligia, M., Tessaro, D., & Cristina, I. (2013). Water adsorption isotherms of microcapsules with hydrolyzed pinhão (Araucaria angustifolia seeds) starch as wall material. Journal of Food Engineering, 114(1), 64-69. https://doi.org/10.1016/j.jfoodeng.2012.07.019
Savitzky, A., & Golay, M. J. E. (1964). Smoothing and Differentiation of Data by Simplified Least Squares Procedures. Analytical Chemistry, 36, 1627-1639. https://doi.org/10.1021/ac60214a047
Smith, S. E. (1947). The Sorption of Water Vapor by High Polymers. Journal of the American Chemical Society, 69(3), 646-651. https://doi.org/10.1021/ja01195a053
Thompson, T. L., Peart, R. M., & Foster, G. H. (1968). Mathematical Simulation of Corn Drying — A New Model. Transactions of the ASAE, 11(4), 582–586.
Tolaba, M. P., Peltzer, M., Enriquez, N., & Lucı́a Pollio, M. (2004). Grain sorption equilibria of quinoa grains. Journal of Food Engineering, 61(3), 365-371. https://doi.org/10.1016/S0260-8774(03)00143-2
Vega-Gálvez, A., Miranda, M., Vergara, J., Uribe, E., Puente, L., & Martinez, E. (2010). Nutrition facts and functional potential of quinoa (Chenopodium quinoa willd.), an ancient Andean grain: a review. Journal Science Food Agricultural (90), 2541-2547. https://doi.org/10.1002/jsfa.4158
Wolf, M., Walker, J. E., & Kapsalis, J. G. (1972). Water vapor sorption hysteresis in dehydrated food. Journal of Agricultural and Food Chemistry, 20(5), 1073-1077. https://doi.org/10.1021/jf60183a021
Yu, L., Mazza, G., & Jayas, D. S. (1999). Moisture sorption characteristics of freeze-dried, osmofreeze-dried, and osmo-air-dried cherries and blueberries. Transactions of the ASAE, 42(1), 141-147. https://doi.org/10.13031/2013.13189
Zhang, Zhen-shan, Li, Xiao-dan, Jía, Hui-jie, Liu, Yu-lan (2022). Moisture sorption isotherms and thermodynamic properties of tiger nuts: An oil-rich tuber. LWT, 167, 113866. https://doi.org/10.1016/j.lwt.2022.113866.
Zeymer, J. S., Corrêa, P. C., Oliveira, G. H. H. d., Araujo, M. E. V. d., Guzzo, F., & Baptestini, F. M. (2022). Moisture sorption isotherms and hysteresis of soybean grains. Acta Scientiarum. Agronomy, 45(1). https://doi.org/10.4025/actasciagron.v45i1.56615