Optimización de operaciones unitarias utilizando líquidos iónicos como agentes materiales de separación

En este trabajo se estudia la modelización y optimización de procesos industriales de separación mediante el empleo de mezclas de líquidos iónicos como disolventes. Los disolventes habitualmente empleados en procesos de absorción o extracción suelen ser componentes orgánicos muy volátiles y dañinos para la salud humana. Las innovadoras propiedades que presentan los líquidos iónicos, los convierten en alternativas adecuadas para solucionar estos problemas. La presión de vapor de estos compuestos es muy baja y apenas varía con la temperatura. Por tanto, estos compuestos apenas se evaporan incluso a temperaturas altas. Esto supone una gran ventaja en cuanto al empleo de estos compuestos como disolventes industriales ya que permite el reciclaje continuo del disolvente al final del proceso sin necesidad de introducir disolvente fresco debido a la evaporación del mismo. Además, al no evaporarse, estos compuestos no suponen un peligro para la salud humana por inhalación; al contrario que otros disolventes como el benceno. El único peligro para la salud que tienen estos compuestos es por tanto el de contacto directo o ingesta, aunque de hecho muchos Líquidos Iónicos son inocuos con lo cual no existe peligro para la salud ni siquiera a través de estas vías. Los procesos de separación estudiados en este trabajo, se rigen por la termodinámica de fases, concretamente el equilibrio líquido-vapor. Para la predicción de los equilibrios se ha optado por el empleo de modelos COSMO (COnductor-like Screening MOdel). Estos modelos tienen su origen en el empleo de la termodinámica de solvatación y en la mecánica cuántica. En el desarrollo de procesos y productos, químicos e ingenieros frecuentemente precisan de la realización de cálculos de predicción de equilibrios de fase. Previamente al desarrollo de los modelos COSMO, se usaban métodos de contribución de grupos como UNIFAC o modelos de coeficientes de actividad como NRTL.La desventaja de estos métodos, es que requieren parámetros de interacción binaria que únicamente pueden obtenerse mediante ajustes por regresión a partir de resultados experimentales. Debido a esto, estos métodos apenas tienen aplicabilidad para compuestos con grupos funcionales novedosos debido a que no se dispone de datos experimentales para llevar a cabo los ajustes por regresión correspondientes. Una alternativa a estos métodos, es el empleo de modelos de solvatación basados en la química cuántica para caracterizar las interacciones moleculares y tener en cuenta la no idealidad de la fase líquida. Los modelos COSMO, permiten la predicción de equilibrios sin la necesidad de ajustes por regresión a partir de resultados experimentales. Debido a la falta de resultados experimentales de equilibrios líquido-vapor de mezclas en las que se ven involucrados los líquidos iónicos, el empleo de modelos COSMO es una buena alternativa para la predicción de equilibrios de mezclas con este tipo de materiales. Los modelos COSMO emplean las distribuciones superficiales de carga polarizada (sigma profiles) de los compuestos involucrados en la mezcla estudiada para la predicción de los coeficientes de actividad de la misma, definiéndose el sigma profile de una molécula como la distribución de probabilidad de densidad de carga superficial de dicha molécula. Dos de estos modelos son COSMO-RS (Realistic Solvation) y COSMO-SAC (Segment Activity Coefficient). El modelo COSMO-RS fue la primera extensión de los modelos de solvatación basados en continuos dieléctricos a la termodinámica de fases líquidas mientras que el modelo COSMO-SAC es una variación de este modelo, tal y como se explicará posteriormente. Concretamente en este trabajo se ha empleado el modelo COSMO-SAC para el cálculo de los coeficientes de actividad de las mezclas estudiadas. Los sigma profiles de los líquidos iónicos se han obtenido mediante el empleo del software de química computacional Turbomole y el paquete químico-cuántico COSMOtherm. El software Turbomole permite optimizar la geometría de la molécula para hallar la configuración más estable mientras que el paquete COSMOtherm permite la obtención del perfil sigma del compuesto mediante el empleo de los datos proporcionados por Turbomole. Por otra parte, los sigma profiles del resto de componentes se han obtenido de la base de datos Virginia Tech-2005 Sigma Profile Database. Para la predicción del equilibrio a partir de los coeficientes de actividad se ha empleado la Ley de Raoult modificada. Se ha supuesto por tanto que la fracción de cada componente en el vapor es proporcional a la fracción del mismo componente en el líquido, dónde la constante de proporcionalidad es el coeficiente de actividad del componente en la mezcla multiplicado por la presión de vapor del componente y dividido por la presión del sistema. Las presiones de vapor de los componentes se han obtenido aplicando la Ley de Antoine. Esta ecuación describe la relación entre la temperatura y la presión de vapor y se deduce a partir de la ecuación de Clausius-Clapeyron. Todos estos datos se han empleado para la modelización de una separación flash usando el algoritmo de Rachford-Rice. El valor de este modelo reside en la deducción de una función que relaciona las constantes de equilibrio, composición total y fracción de vapor. Para llevar a cabo la implementación del modelado matemático descrito, se ha programado un código empleando el software MATLAB de análisis numérico. Para comprobar la fiabilidad del código programado, se compararon los resultados obtenidos en la predicción de equilibrios de mezclas mediante el código con los resultados obtenidos mediante el simulador ASPEN PLUS de procesos químicos. Debido a la falta de datos relativos a líquidos iónicos en la base de datos de ASPEN PLUS, se han introducido estos componentes como pseudocomponentes, de manera que se han introducido únicamente los datos necesarios de estos componentes para realizar las simulaciones. El modelo COSMO-SAC se encuentra implementado en ASPEN PLUS, de manera que introduciendo los sigma profiles, los volúmenes de la cavidad y las presiones de vapor de los líquidos iónicos, es posible predecir equilibrios líquido-vapor en los que se ven implicados este tipo de materiales. De esta manera pueden compararse los resultados obtenidos con ASPEN PLUS y como el código programado en MATLAB y comprobar la fiabilidad del mismo. El objetivo principal del presente Trabajo Fin de Máster es la optimización de mezclas multicomponente de líquidos iónicos para maximizar la eficiencia de procesos de separación y minimizar los costes de los mismos. La estructura de este problema es la de un problema de optimización no lineal con variables discretas y continuas, es decir, un problema de optimización MINLP (Mixed Integer Non-Linear Programming). Tal y como se verá posteriormente, el modelo matemático de este problema es no lineal. Por otra parte, las variables del mismo son tanto continuas como binarias. Las variables continuas se corresponden con las fracciones molares de los líquidos iónicos presentes en las mezclas y con el caudal de la mezcla de líquidos iónicos. Por otra parte, también se ha introducido un número de variables binarias igual al número de líquidos iónicos presentes en la mezcla. Cada una de estas variables multiplican a las fracciones molares de sus correspondientes líquidos iónicos, de manera que cuando dicha variable es igual a 1, el líquido se encuentra en la mezcla mientras que cuando dicha variable es igual a 0, el líquido iónico no se encuentra presente en dicha mezcla. El empleo de este tipo de variables obliga por tanto a emplear algoritmos para la resolución de problemas de optimización MINLP ya que si todas las variables fueran continuas, bastaría con el empleo de algoritmos para la resolución de problemas de optimización NLP (Non-Linear Programming). Se han probado por tanto diversos algoritmos presentes en el paquete OPTI Toolbox de MATLAB para comprobar cuál es el más adecuado para abordar este problema. Finalmente, una vez validado el código programado, se han optimizado diversas mezclas de líquidos iónicos para lograr la máxima recuperación de compuestos aromáticos en un proceso de absorción de mezclas orgánicas. También se ha usado este código para la minimización del coste correspondiente a la compra de los líquidos iónicos de la mezcla de disolventes empleada en la operación de absorción. En este caso ha sido necesaria la introducción de restricciones relativas a la recuperación de aromáticos en la fase líquida o a la pureza de la mezcla obtenida una vez separada la mezcla de líquidos iónicos. Se han modelizado los dos problemas descritos previamente (maximización de la recuperación de Benceno y minimización del coste de operación) empleando tanto únicamente variables continuas (correspondientes a las fracciones o cantidades molares de los líquidos iónicos) como variables continuas y binarias (correspondientes a cada uno de los líquidos iónicos implicados en las mezclas).

Modulação da diferenciação neural de células troco embrionárias por transientes de cálcio intracelulares: papéis dos receptores purinérgicos e de canais de cálcio voltagem-dependentes

Receptores purinérgicos e canais de cálcio voltagem-dependentes estão envolvidos em diversos processos biológicos como na gastrulação, durante o desenvolvimento embrionário, e na diferenciação neural. Quando ativados, canais de cálcio voltagem-dependentes e receptores purinérgicos do tipo P2, ativados por nucleotídeos, desencadeiam transientes de cálcio intracelulares controlando diversos processos biológicos. Neste trabalho, nós estudamos a participação de canais de cálcio voltagem-dependentes e receptores do tipo P2 na geração de transientes de cálcio espontâneos e sua regulação na expressão de fatores de transcrição relacionados com a neurogênese utilizando como modelo células tronco (CTE) induzidas à diferenciação em células tronco neurais (NSC) com ácido retinóico. Descrevemos que CTE indiferenciadas podem ter a proliferação acelerada pela ativação de receptores P2X7, enquanto que a expressão e a atividade desse receptor precisam ser inibidas para o progresso da diferenciação em neuroblasto. Além disso, ao longo da diferenciação neural, por análise em tempo real dos níveis de cálcio intracelular livre identificamos 3 padrões de oscilações espontâneas de cálcio (onda, pico e unique), e mostramos que ondas e picos tiveram a frequência e amplitude aumentadas conforme o andamento da diferenciação. Células tratadas com o inibidor do receptor de inositol 1,4,5-trifosfato (IP3R), Xestospongin C, apresentaram picos mas não ondas, indicando que ondas dependem exclusivamente de cálcio oriundo do retículo endoplasmático pela ativação de IP3R. NSC de telencéfalo de embrião de camundongos transgênicos ou pré-diferenciadas de CTE tratadas com Bz-ATP, o agonista do receptor P2X7, e com 2SUTP, agonista de P2Y2 e P2Y4, aumentaram a frequência e a amplitude das oscilações espontâneas de cálcio do tipo pico. Dados, obtidos por microscopia de luminescência, da expressão em tempo real de gene repórter luciferase fusionado à Mash1 e Ngn2 revelou que a ativação dos receptores P2Y2/P2Y4 aumentou a expressão estável de Mash1 enquanto que ativação do receptor P2X7 levou ao aumento de Ngn2. Além disso, células na presença do quelante de cálcio extracelular (EGTA) ou do depletor dos estoques intracelulares de cálcio do retículo endoplasmático (thapsigargin) apresentaram redução na expressão de Mash1 e Ngn2, indicando que ambos são regulados pela sinalização de cálcio. A investigação dos canais de cálcio voltagem-dependentes demonstrou que o influxo de cálcio gerado por despolarização da membrana de NSC diferenciadas de CTE é decorrente da ativação de canais de cálcio voltagem-dependentes do tipo L. Além disso, esse influxo pode controlar o destino celular por estabilizar expressão de Mash1 e induzir a diferenciação neuronal por fosforilação e translocação do fator de transcrição CREB. Esses dados sugerem que os receptores P2X7, P2Y2, P2Y4 e canais de cálcio voltagem-dependentes do tipo L podem modular as oscilações espontâneas de cálcio durante a diferenciação neural e consequentemente alteram o padrão de expressão de Mash1 e Ngn2 favorecendo a decisão do destino celular neuronal.

Expressão dos genes codificadores de canais de sódio Nav 1.7, Nav 1.8 e Nav 1.9 em portadores da Síndrome de Ardência Bucal

A síndrome de ardência bucal (SAB) é uma condição caracterizada pelo sintoma de ardência na mucosa oral, na ausência de qualquer sinal clínico. Sua etiologia é desconhecida e, até o momento, não dispõe de tratamento efetivo. Há entretanto características de doença neuropática que justificam investigações nesse sentido. O objetivo desse estudo foi mensurar a expressão gênica dos receptores de canais de sódio, Nav 1.7, Nav 1.8 e Nav 1.9, nos pacientes portadores de SAB. A casuística foi composta por dois grupos sendo o grupo de estudo composto por 12 pacientes portadores de SAB, selecionados através do critério estabelecido pela International Headache Society, em 2013 e o grupo controle composto por 4 pacientes não portadores de SAB. As amostras analisadas foram coletadas do dorso lingual, por meio de biópsia realizada com punch de 3 mm e profundidade de 3 mm, estas foram submetidas ao método de análise RT-PCR em tempo real. A expressividade dos genes de canais de sódio foi avaliada nos indivíduos portadores de SAB em relação aos do grupo controle, sendo esta calculada a partir da normalização dos dados da quantificação destes com os da expressão do gene constitutivo (GAPDH), pelo método de Cicle Threshold comparativo e analisados estatisticamente por meio do teste estatístico Mann-Whitnney. Observou-se o aumento da expressão gênica do Nav 1.7 (fold-change = 38.70) e diminuição da expressão gênica do Nav 1.9 (fold-change = 0.89), porém sem diferenças estatisticamente significativas entre os grupos analisados. O gene Nav 1.8 não foi expresso em nenhuma das amostras analisadas. O Nav 1.7 expressa-se tanto em neurônios nociceptivos quanto no sistema nervoso autônomo e mutações no Nav 1.9 tem sido associada a perda de percepção dolorosa. Os resultados obtidos embora não estatisticamente significativos são compatíveis com as características da doença, justificando a extensão dos estudos na linha expressão de genes codificadores dos canais de sódio em pacientes com SAB.

Dinâmica molecular de sistemas iônicos poliatômicos: modelos polarizável e não-polarizável

Simulações de sais de carbonato fundidos pelo método de Dinâmica Molecular (MD) foram efetuadas com o modelo polarizável de cargas flutuantes (FC). O modelo de cargas flutuantes implementa os efeitos de polarização pelo método de Lagrangiano estendido, onde as variáveis extras são as próprias cargas parciais do íon poliatômico. O modelo FC foi parametrizado por meio de cálculos ab inito, aplicado ao ânion carbonato. Cálculos de Química Quântica ab initio foram utilizados para corroborar o modelo proposto para o ânion carbonato. Os sistemas investigados consistem em misturas de carbonatos alcalinos fundidos, Li2CO3/K2CO3, os quais são utilizados como eletrólitos em células a combustível. As simulações MD foram utilizadas para verificar o efeito da polarização dos ânions sobre a estrutura e dinâmica do líquido. Estudamos o efeito da inclusão de polarização sobre a condutividade do eletrólito.