En el presente trabajo de tesis se elabora un modelo matemático que predice el coeficiente de partición de proteínas en sistemas de dos fases acuosas (ATPS), a partir del efecto de hidrofobicidad de las proteínas y la diferencia de energía de solvatación electrostática entre las fases . Para esto, se utilizan los datos experimentales de partición de 11 proteínas, en los sistemas PEG+fosfato, PEG+sulfato, PEG+citrato y PEG+dextrano, cada uno con concentración alta (8.8% p/p), intermedia (0.6% p/p) y nula (0% p/p) de NaCl, definiendo un total de 12 sistemas de características diferentes. Se propusieron 26 modelos en el presente trabajo, obteniendo los mejores resultados con el modelo , el cual disminuye el error de predicción en 6 de los 12 sistemas estudiados, con respecto a los resultados obtenidos por Salgado et al. (2008), donde sólo se utiliza la hidrofobicidad de las proteínas como variable. Los sistemas más favorecidos con la inclusión simultánea de ambas variables corresponden a aquellos con alta (8.8% p/p) concentración de NaCl, donde disminuyó el error de predicción entre un 4.0% y un 34.1% con respecto a los resultados de Salgado et al. Por lo tanto, la adición de una alta concentración de NaCl constituye una fácil alternativa para aumentar la calidad predictiva del modelo, siempre que el diseño lo permita. Debido a que los fenómenos fisicoquímicos involucrados en la partición de proteínas en un sistema PEG+Dextrano (formado por dos polímeros) son distintos a los involucrados en los sistemas PEG+sal, se elaboró un modelo aplicable sólo al sistema PEG+Dextrano, con el cual disminuyó el error de predicción entre un 1.7% y un 29.5% con respecto al de Salgado et al. En este ATPS, no sólo se logró mejorar la predicción del coeficiente de partición para todas las concentraciones de NaCl consideradas, sino que fue posible disminuir el error de predicción con mayor frecuencia que en los demás ATPS, considerando los 26 modelos propuestos. La complejidad del fenómeno involucrado en la partición de proteínas en los sistemas PEG+sal, con la consecuente dificultad de mejorar modelos predictivos aplicados a estos sistemas, explica porqué los ATPS de soluciones poliméricas no han sido reemplazados a nivel industrial por los económicos sistemas salinos. El aporte del presente trabajo de tesis no sólo consiste en presentar una herramienta matemática para predecir el coeficiente de partición de proteínas en ATPS, sino que además, el modelo indica en qué proporción se debe modificar la hidrofobicidad de las proteínas y la diferencia de energía de solvatación electrostática para optimizar la separación.

El argumento fundamental de esta tesis es el estudio de la solvatación iónica por medio de cálculos con ordenador. Tres lineas de investigación han sido seguidas:(i) Solvatación y mobilidad ionica. Las características principales del processo de intercambio entre la primera y la segunda capa de hidratación iónica para Li+ en agua se ha encontrado ser independiente del estado termodinamico en gran medida. Ha sido demostrado que el desplazamiento cuadrático medio de moléculas pertenecientes a complejos inertes está caracterizado por un largo transitorio debido a la lenta relajación rotacional del complejo. El incremento del coeficiente de difusión iónico debido a los intecambios en la capa de solvatación ha sido calculado por primera vez en el caso de Li+ y Na+. Finalmente, se han derivado leyes de probabilidad que ponen en relación la estereoquímica y la velocidad iónica instantanea.(ii) Plastificantes. Se propone un nuevo procedimiento para el desarrollo de campos de fuerza intramoleculares, que funciona satisfactoriamente en el caso de dos moleculas de interés en las Batterias a Iones de Litio: carbonato de etileno y -butirolactona. Respecto a la solvatación de Li+ en los dos solventes, el ión está coordenado por 4 moleculas a través del oxigeno del carbonilo con pequeñas distorsiones de la geometría molecular. La nueva asignación de los modos vibracionales hecha para las dos moléculas ha permitido calcular los cambios inducidos por el ión litio, explicando varias caracteristicas de los espectros esperimentales.(iii) Polarización. Se ha estudiado la eficacia de los metodos de polarización más comunes para simulaciones de Dinámica Molecular en dímeros ión-molécula, usando calculos ab initio como referencia. En lugar de centrarnos en la superficie de energía potencial completa (procedimiento típico), se ha considerado solo la parte electrostática. Se han desarrollado nuevos modelos polarizables para agua y tetracloruro de carbono, que reproducen el comportamiento de sistemas carga-molécula. Ha sido encontrado que, en el caso de dímeros ión-molécula, se requiere una corrección de amortiguamento de la polarización a cortas distancias. El método de los dipolos puntuales junto al método de amortiguamento de Thole reproduce satisfactoriamente las características principales para cationes y aniones atómicos. / The underlying topic of this thesis is the study of ion solvation by means of computer calculations. Three lines of investigation have been followed:(i) Solvation and Ionic Mobility. The main features of the exchange process between first and second ionic hydration shells for Li+ in water have been found to be largely independent of the thermodynamic state. It has been shown that the mean square displacement of molecules belonging to inert complexes is characterized by a long transient due to the slow rotational relaxation of the complex. The increase of the ionic diffusion coefficient due to solvation shell exchanges has been computed for the first time in the case of Li+ and Na+. Finally, probability laws have been derived which relate the stereochemistry and the instantaneous ionic velocity.(ii) Plasticizers. A new approach for the development of intramolecular force fields is proposed, which performs satisfactorily in the case of two molecules of interest for Lithium Ion Batteries: ethylene carbonate and -butyrolactone. Concerning the solvation of Li+ in both solvents, it is coordinated by 4 molecules through the carbonyl oxygen with slight distorsions of the molecular geometry. The new vibrational mode assignment performed for both molecules has allowed to compute the vibrational shifts induced by the lithium ion, explaining a number of features present in the experimental spectra.(iii) Polarization. The performance of the most commonly used polarization methods for Molecular Dynamics simulation is studied for ion-molecule dimers, using ab initio calculations as benchmark. Instead of focusing on the full potential energy surface (the standard approach), only the electrostatic part is considered. New polarizable models have been developed for water and carbon tetrachloride, which reproduce the behaviour of charge-molecule systems. In the case of ion-molecule dimers it has been found that a polarization damping correction is required at short distances. The point dipole method in conjunction with the Thole damping scheme reproduces rather satisfactorily the main features both for atomic cations and anions.

La presente investigación se centra en dos ejes principales. El primero está relacionado con metodologías de muestreo para la simulación de moléculas modeladas como cuerpos rígidos. Utilizando dichas metodologías fue posible desarrollar un segundo eje de investigación, el cual tiene como objeto de estudio a los cálculos de energía libre de solvatación y de propiedades termodinámicas relacionadas. El Capítulo 3 tiene como objeto de estudio al método de Monte Carlo Híbrido (MCH) aplicado a la simulación de cuerpos rígidos. Se realizó un estudio sistemático de la eficiencia de MCH, y se desarrolló un procedimiento de optimización en línea que permite obtener parámetros óptimos o cuasi-óptimos de MCH para los sistemas en fase líquida investigados: propano, benceno y agua. La implementación de MCH, realizada en LAMMPS, fue motivada por aspectos objetables del método del ensamble expandido: usualmente se combinan pasos de DM y MC violando la condición del balance detallado de probabilidad y, además, no se realiza una equilibración luego de las transiciones entre los macroestados. En el Capítulo 4 se realiza una validación empírica del método del ensamble expandido, para lo cual se determinó un perfil de energía libre “correcto”, que se obtuvo a partir de muestras independientes de cada macroestado empleando MCH. Asimismo, se evaluaron diferentes frecuencias de transición entre los macroestados, y en todos los casos el perfil de energía libre coincide con el resultado de referencia. De este análisis se desprende que los aspectos objetables no introducen errores sistemáticos importantes, y se concluye que el muestreo secuencial del ensamble expandido con DM y MC arroja resultados confiables. La validación del método del ensamble expandido permitió abordar el desafío de predecir la solvatación en líquidos iónicos y calcular propiedades macroscópicas requeridas en el diseño de procesos químicos, tal como se detalla en el Capítulo 5. El CO2 proveniente de procesos de combustión constituye una de las fuentes principales de gases de efecto invernadero de generación antropogénica y el potencial de [emim][B(CN)4] para su captura ha sido recientemente reportado. Con fines comparativos, se analizaron los solventes [emim][B(CN)4] y [emim][Tf2N], tanto en la condición de dilución infinita, como en altas concentraciones. Las simulaciones confirmaron lo observado experimentalmente, en tanto se requiere de una presión considerablemente menor en [emim][B(CN)4] para obtener una misma molalidad de CO2. Por otra parte, se predijeron los coeficientes de actividad a dilución infinita de los solutos hexano, benceno, ciclohexano, etanol y agua en [emim][B(CN)4], y se compararon con valores experimentales reportados recientemente en la literatura. Finalmente, en el Capítulo 6 se exponen las conclusiones y consideraciones finales de la presente investigación, como también algunas sugerencias para trabajos futuros. En el Capítulo 1 se introducen conceptos de Mecánica Estadística y Simulación Molecular (SM) que han sido fundamentales para el desarrollo de esta tesis. A esto antecede una revisión bibliográfica en la cual se destacan factores determinantes para que la Simulación Molecular constituya actualmente una herramienta de investigación en la Ingeniería Química. En el Capítulo 2 se analizan y desarrollan métodos de Dinámica Molecular (DM) de cuerpos rígidos en los ensambles NVE y NVT, motivado principalmente por resultados insatisfactorios obtenidos con un integrador NVT implementado en softwares muy utilizados por la comunidad de SM, y que acopla termostatos independientes a los grados de libertad traslacionales y rotacionales. En particular, para el sistema de agua líquida, dicha estrategia no logra reproducir la temperatura especificada, inclusive para un paso de tiempo pequeño de 1 fs. Ante esto, es natural preguntarse si el problema podría estar vinculado con el integrador NVE, que es la base del método NVT. Una anomalía detectada fue la violación del teorema de equipartición de la energía cinética, lo cual tiene origen en la discretización de la integración numérica. Esto ya ha sido reportado en la literatura, pero sólo se mostró el efecto de esos errores sobre diversas propiedades termodinámicas y, hasta donde se sabe, no existen trabajos que hayan revisado el método NVT mencionado. Esto motivó la derivación de un integrador NVE con una formulación matemática simplificada de las ecuaciones de movimiento. En la presente investigación se demuestra que el origen de las anomalías está relacionado con la disparidad de los momentos de inercia. Además, se presenta un procedimiento simple para disminuir el efecto de los errores de discretización, que viene dado por la redistribución de las masas atómicas. / The present investigation focuses on two main topics. The first is related to sampling strategies for the simulation of molecules modeled as rigid bodies. This trigger a second line of research, which aims to calculate solvation free energies and related thermodynamic properties. Chapter 1 introduces concepts of Statistical Mechanics and Molecular Simulation (MS) needed for the development of this thesis. This is preceded by a literature review that highlights the major factors that drove Molecular Simulation to become a research tool in Chemical Engineering. Chapter 2 is devoted to the analysis and development of integrators for Molecular Dynamics (MD) with rigid bodies in the NVE and NVT ensembles. This was mainly motivated by unsatisfactory results obtained with an NVT integrator that couples independent thermostats to the translational and rotational degrees of freedom, which is implemented in software packages widely used by the MS community. In the case of liquid water, such strategy fails to reproduce the specified temperature, even for a time step as small as 1 fs. Thus, it is natural to ask whether the problem could be linked to the NVE integrator, which is the basis of the NVT method. An anomalous violation of the kinetic-energy equipartition theorem was detected, originated by the time discretization employed in the numerical integration. This has already been reported in the literature, but the discussion was limited to the effect of these errors on various thermodynamic properties and, as far as is known, no one has reviewed the NVT method mentioned above. This motivated the derivation of an NVE integrator with a simplified mathematical formulation of the equations of motion. In the present investigation it is shown that the origin of the anomalies is related to a large difference in the moments of inertia of a rigid molecule. In addition, a simple procedure for reducing the effect of discretization errors is presented, which consists in redistributing the atomic masses. In Chapter 3, the Hybrid Monte Carlo (HMC) method is applied to the simulation of rigid bodies. A systematic study on the efficiency of HMC was carried out. Then, an on-the-fly optimization procedure was developed, which allows obtaining optimal or quasi-optimal parameters for HMC. The liquid-phase systems assessed were propane, benzene, and water. The implementation of HMC, carried out in LAMMPS, was motivated by questionable aspects of the expanded ensemble method: MD and MC steps are usually combined in a way that violates the condition of detailed balance of probabilities and, in addition, equilibration is not performed after a transition between macrostates takes place. In Chapter 4, an empirical validation of the expanded ensemble method is shown, for which a “correct” free energy profile was determined from independent samples drawn at each macrostate using HMC. Likewise, different transition frequencies between macrostates were evaluated, and in all cases the free energy profile agrees with the reference result. From this analysis, it seems that the questionable aspects do not introduce important systematic errors and the sequential sampling of an expanded ensemble alternating MD and MC yields reliable results. The validation of the expanded ensemble method allowed to tackle, in Chapter 5, the challenge of predicting solvation in ionic liquids and macroscopic properties required in the design of chemical processes. Carbon dioxide produced in combustion processes is one of the main sources of anthropogenic generation of greenhouse gases. On the other hand, the potential of [emim][B(CN)4] for its capture has recently been reported. For comparison, the solvents [emim][B(CN)4] and [emim][Tf2N] were analyzed, both in the condition of infinite dilution and at high concentrations. The simulations confirmed what was observed experimentally, that is, a considerably lower pressure is required in [emim][B(CN)4] to obtain the same molality of CO2. The infinite-dilution activity coefficients of solutes hexane, benzene, cyclohexane, ethanol, and water in [emim][B(CN)4] were also calculated, and a comparison was made with experimental values recently reported in the literature. Finally, Chapter 6 presents the conclusions and final considerations of the present investigation, as well as some suggestions for future work.