Modeling of electrolyte sorption – from phase equilibria to dynamic separation systems

In this thesis, general approach is devised to model electrolyte sorption from aqueous solutions on solid materials. Electrolyte sorption is often considered as unwanted phenomenon in ion exchange and its potential as an independent separation method has not been fully explored. The solid sorbents studied here are porous and non-porous organic or inorganic materials with or without specific functional groups attached on the solid matrix. Accordingly, the sorption mechanisms include physical adsorption, chemisorption on the functional groups and partition restricted by electrostatic or steric factors. The model is tested in four Cases Studies dealing with chelating adsorption of transition metal mixtures, physical adsorption of metal and metalloid complexes from chloride solutions, size exclusion of electrolytes in nano-porous materials and electrolyte exclusion of electrolyte/non-electrolyte mixtures. The model parameters are estimated using experimental data from equilibrium and batch kinetic measurements, and they are used to simulate actual single-column fixed-bed separations. Phase equilibrium between the solution and solid phases is described using thermodynamic Gibbs-Donnan model and various adsorption models depending on the properties of the sorbent. The 3-dimensional thermodynamic approach is used for volume sorption in gel-type ion exchangers and in nano-porous adsorbents, and satisfactory correlation is obtained provided that both mixing and exclusion effects are adequately taken into account. 2-Dimensional surface adsorption models are successfully applied to physical adsorption of complex species and to chelating adsorption of transition metal salts. In the latter case, comparison is also made with complex formation models. Results of the mass transport studies show that uptake rates even in a competitive high-affinity system can be described by constant diffusion coefficients, when the adsorbent structure and the phase equilibrium conditions are adequately included in the model. Furthermore, a simplified solution based on the linear driving force approximation and the shrinking-core model is developed for very non-linear adsorption systems. In each Case Study, the actual separation is carried out batch-wise in fixed-beds and the experimental data are simulated/correlated using the parameters derived from equilibrium and kinetic data. Good agreement between the calculated and experimental break-through curves is usually obtained indicating that the proposed approach is useful in systems, which at first sight are very different. For example, the important improvement in copper separation from concentrated zinc sulfate solution at elevated temperatures can be correctly predicted by the model. In some cases, however, re-adjustment of model parameters is needed due to e.g. high solution viscosity.

Do We Understand Hydrodynamic Dispersion in Sedimenting Suspensions?

Microstructural changes, that is an important feature for the understanding of the velocity variance in sedimentation is investigated with numerical simulations. The simulations are used to describe velocity fluctuations and hydrodynamic dispersion in a suspension of interacting point-particles sedimenting in a rectangular box with periodic sides and impenetrable bottom and top. It is observed how the positions of the particles evolve in a finite container. The suspension that was initially random in the gravity direction only, tends to be fully randomized as a result of the relative arrangements of the particles and the hydrodynamic interactions between them. The computer simulations, based on statistics over a significant number of particle configurations, suggest velocity variances and diffusivities dependent on the size of the simulated system but with anisotropy in velocity fluctuations and diffusion coefficients nearly independent of the box size.

Saline distribution during multicomponent salting in pre-cooked quail eggs

The relationship of NaCl with problems of arterial hypertension has led to a reduction in the levels of this salt in food production. KCl has been used as a partial substitute for NaCl since it cannot be completely substituted without affecting the acceptability of the end product. In this study, the diffusion that occurs during quail egg salting in static and stirred brine was simulated. The mathematical model used was based on a generalization of the Fick's 2nd law, and the COMSOL Multiphysics software was used to simulate the diffusion in the NaCl-KCl-water system. The deviations in the simulated data and experimental data were 2.50% for NaCl and 6.98% for KCl in static brine, while in the stirred brine they were 3.48% for NaCl and 4.72% for KCl. The simulation results presented good agreement with the experimental values and validated the predictive capacity of the model.

Le rendu en demi-ton avec sensibilité à la structure

Dans ce mémoire nous allons présenter une méthode de diffusion d’erreur originale qui peut reconstruire des images en demi-ton qui plaisent à l’œil. Cette méthode préserve des détails fins et des structures visuellement identifiables présentes dans l’image originale. Nous allons tout d’abord présenter et analyser quelques travaux précédents afin de montrer certains problèmes principaux du rendu en demi-ton, et nous allons expliquer pourquoi nous avons décidé d’utiliser un algorithme de diffusion d’erreur pour résoudre ces problèmes. Puis nous allons présenter la méthode proposée qui est conceptuellement simple et efficace. L’image originale est analysée, et son contenu fréquentiel est détecté. Les composantes principales du contenu fréquentiel (la fréquence, l’orientation et le contraste) sont utilisées comme des indices dans un tableau de recherche afin de modifier la méthode de diffusion d’erreur standard. Le tableau de recherche est établi dans un étape de pré-calcul et la modification est composée par la modulation de seuil et la variation des coefficients de diffusion. Ensuite le système en entier est calibré de façon à ce que ces images reconstruites soient visuellement proches d’images originales (des aplats d’intensité constante, des aplats contenant des ondes sinusoïdales avec des fréquences, des orientations et des constrastes différents). Finalement nous allons comparer et analyser des résultats obtenus par la méthode proposée et des travaux précédents, et démontrer que la méthode proposée est capable de reconstruire des images en demi-ton de haute qualité (qui préservent des structures) avec un traitement de temps très faible.

Nanostructure des particules polymériques : aspects physiques, chimiques et biologiques

Les nanotechnologies appliquées aux sciences pharmaceutiques ont pour but d’améliorer l’administration de molécules actives par l’intermédiaire de transporteurs nanométriques. Parmi les différents types de véhicules proposés pour atteindre ce but, on retrouve les nanoparticules polymériques (NP) constituées de copolymères “en bloc”. Ces copolymères permettent à la fois l’encapsulation de molécules actives et confèrent à la particule certaines propriétés de surface (dont l’hydrophilicité) nécessaires à ses interactions avec les milieux biologiques. L’architecture retenue pour ces copolymères est une structure constituée le plus fréquemment de blocs hydrophiles de poly(éthylène glycol) (PEG) associés de façon linéaire à des blocs hydrophobes de type polyesters. Le PEG est le polymère de choix pour conférer une couronne hydrophile aux NPs et son l’efficacité est directement liée à son organisation et sa densité de surface. Néanmoins, malgré les succès limités en clinique de ces copolymères linéaires, peu de travaux se sont attardés à explorer les effets sur la structure des NPs d’architectures alternatives, tels que les copolymères en peigne ou en brosse. Durant ce travail, plusieurs stratégies ont été mises au point pour la synthèse de copolymères en peigne, possédant un squelette polymérique polyesters-co-éther et des chaines de PEG liées sur les groupes pendants disponibles (groupement hydroxyle ou alcyne). Dans la première partie de ce travail, des réactions d’estérification par acylation et de couplage sur des groupes pendants alcool ont permis le greffage de chaîne de PEG. Cette méthode génère des copolymères en peigne (PEG-g-PLA) possédant de 5 à 50% en poids de PEG, en faisant varier le nombre de chaînes branchées sur un squelette de poly(lactique) (PLA). Les propriétés structurales des NPs produites ont été étudiées par DLS, mesure de charge et MET. Une transition critique se situant autour de 15% de PEG (poids/poids) est observée avec un changement de morphologie, d’une particule solide à une particule molle (“nanoagrégat polymére”). La méthode de greffage ainsi que l’addition probable de chaine de PEG en bout de chaîne principale semblent également avoir un rôle dans les changements observés. L’organisation des chaînes de PEG-g-PLA à la surface a été étudiée par RMN et XPS, méthodes permettant de quantifier la densité de surface en chaînes de PEG. Ainsi deux propriétés clés que sont la résistance à l’agrégation en conditions saline ainsi que la résistance à la liaison aux protéines (étudiée par isothermes d’adsorption et microcalorimétrie) ont été reliées à la densité de surface de PEG et à l’architecture des polymères. Dans une seconde partie de ce travail, le greffage des chaînes de PEG a été réalisé de façon directe par cyclo-adition catalysée par le cuivre de mPEG-N3 sur les groupes pendants alcyne. Cette nouvelle stratégie a été pensée dans le but de comprendre la contribution possible des chaines de PEG greffées à l’extrémité de la chaine de PLA. Cette librairie de PEG-g-PLA, en plus d’être composée de PEG-g-PLA avec différentes densités de greffage, comporte des PEG-g-PLA avec des PEG de différent poids moléculaire (750, 2000 et 5000). Les chaines de PEG sont seulement greffées sur les groupes pendants. Les NPs ont été produites par différentes méthodes de nanoprécipitation, incluant la nanoprécipitation « flash » et une méthode en microfluidique. Plusieurs variables de formulation telles que la concentration du polymère et la vitesse de mélange ont été étudiées afin d’observer leur effet sur les caractéristiques structurales et de surface des NPs. Les tailles et les potentiels de charges sont peu affectés par le contenu en PEG (% poids/poids) et la longueur des chaînes de PEG. Les images de MET montrent des objets sphériques solides et l'on n’observe pas d’objets de type agrégat polymériques, malgré des contenus en PEG comparable à la première bibliothèque de polymère. Une explication possible est l’absence sur ces copolymères en peigne de chaine de PEG greffée en bout de la chaîne principale. Comme attendu, les tailles diminuent avec la concentration du polymère dans la phase organique et avec la diminution du temps de mélange des deux phases, pour les différentes méthodes de préparation. Finalement, la densité de surface des chaînes de PEG a été quantifiée par RMN du proton et XPS et ne dépendent pas de la méthode de préparation. Dans la troisième partie de ce travail, nous avons étudié le rôle de l’architecture du polymère sur les propriétés d’encapsulation et de libération de la curcumine. La curcumine a été choisie comme modèle dans le but de développer une plateforme de livraison de molécules actives pour traiter les maladies du système nerveux central impliquant le stress oxydatif. Les NPs chargées en curcumine, montrent la même transition de taille et de morphologie lorsque le contenu en PEG dépasse 15% (poids/poids). Le taux de chargement en molécule active, l’efficacité de changement et les cinétiques de libérations ainsi que les coefficients de diffusion de la curcumine montrent une dépendance à l’architecture des polymères. Les NPs ne présentent pas de toxicité et n’induisent pas de stress oxydatif lorsque testés in vitro sur une lignée cellulaire neuronale. En revanche, les NPs chargées en curcumine préviennent le stress oxydatif induit dans ces cellules neuronales. La magnitude de cet effet est reliée à l’architecture du polymère et à l’organisation de la NP. En résumé, ce travail a permis de mettre en évidence quelques propriétés intéressantes des copolymères en peigne et la relation intime entre l’architecture des polymères et les propriétés physico-chimiques des NPs. De plus les résultats obtenus permettent de proposer de nouvelles approches pour le design des nanotransporteurs polymériques de molécules actives.

Effect of polymer architecture on Curcumin 1 encapsulation and release from PEGylated polymer nanoparticles: toward a drug delivery nano-platform to the CNS

We developed a nanoparticles (NPs) library from poly(ethylene glycol)–poly lactic acid comb-like polymers with variable amount of PEG. Curcumin was encapsulated in the NPs with a view to develop a delivery platform to treat diseases involving oxidative stress affecting the CNS. We observed a sharp decrease in size between 15 and 20% w/w of PEG which corresponds to a transition from a large solid particle structure to a “micelle-like” or “polymer nano-aggregate” structure. Drug loading, loading efficacy and release kinetics were determined. The diffusion coefficients of curcumin in NPs were determined using a mathematical modeling. The higher diffusion was observed for solid particles compared to “polymer nano-aggregate” particles. NPs did not present any significant toxicity when tested in vitro on a neuronal cell line. Moreover, the ability of NPs carrying curcumin to prevent oxidative stress was evidenced and linked to polymer architecture and NPs organization. Our study showed the intimate relationship between the polymer architecture and the biophysical properties of the resulting NPs and sheds light on new approaches to design efficient NP-based drug carriers.

Influence of the pseudopotential on the properties of liquid rubidium at 315 K

The influence of the pseudopotential on both the structure and the self-diffusion of liquid rubidium at the melting point has been investigated by means of molecular-dynamics calculations. The model potential considered has been computed from the pseudopotential of Ashcroft, the dielectric function of Geldart and Vosko, and a Born-Mayer term. Four different values for the core radius which enters as input in the pseudopotential have been considered. In this way we have been able to observe and interpret the effect of this contribution on the properties of the liquid.

From a discrete to a continuum model of cell dynamics in one dimension

Multiscale modeling is emerging as one of the key challenges in mathematical biology. However, the recent rapid increase in the number of modeling methodologies being used to describe cell populations has raised a number of interesting questions. For example, at the cellular scale, how can the appropriate discrete cell-level model be identified in a given context? Additionally, how can the many phenomenological assumptions used in the derivation of models at the continuum scale be related to individual cell behavior? In order to begin to address such questions, we consider a discrete one-dimensional cell-based model in which cells are assumed to interact via linear springs. From the discrete equations of motion, the continuous Rouse [P. E. Rouse, J. Chem. Phys. 21, 1272 (1953)] model is obtained. This formalism readily allows the definition of a cell number density for which a nonlinear "fast" diffusion equation is derived. Excellent agreement is demonstrated between the continuum and discrete models. Subsequently, via the incorporation of cell division, we demonstrate that the derived nonlinear diffusion model is robust to the inclusion of more realistic biological detail. In the limit of stiff springs, where cells can be considered to be incompressible, we show that cell velocity can be directly related to cell production. This assumption is frequently made in the literature but our derivation places limits on its validity. Finally, the model is compared with a model of a similar form recently derived for a different discrete cell-based model and it is shown how the different diffusion coefficients can be understood in terms of the underlying assumptions about cell behavior in the respective discrete models.

Classifying general nonlinear force laws in cell-based models via the continuum limit

though discrete cell-based frameworks are now commonly used to simulate a whole range of biological phenomena, it is typically not obvious how the numerous different types of model are related to one another, nor which one is most appropriate in a given context. Here we demonstrate how individual cell movement on the discrete scale modeled using nonlinear force laws can be described by nonlinear diffusion coefficients on the continuum scale. A general relationship between nonlinear force laws and their respective diffusion coefficients is derived in one spatial dimension and, subsequently, a range of particular examples is considered. For each case excellent agreement is observed between numerical solutions of the discrete and corresponding continuum models. Three case studies are considered in which we demonstrate how the derived nonlinear diffusion coefficients can be used to (a) relate different discrete models of cell behavior; (b) derive discrete, intercell force laws from previously posed diffusion coefficients, and (c) describe aggregative behavior in discrete simulations.

Determination of eddy-diffusivity in the lowermost stratosphere

We present a 2D-advection-diffusion model that simulates the main transport pathways influencing tracer distributions in the lowermost stratosphere (LMS). The model describes slow diabatic descent of aged stratospheric air, vertical (cross-isentropic) and horizontal (along isentropes) diffusion within the LMS and across the tropopause using equivalent latitude and potential temperature coordinates. Eddy diffusion coefficients parameterize the integral effect of dynamical processes leading to small scale turbulence and mixing. They were specified by matching model simulations to observed CO distributions. Interestingly, the model suggests mixing across isentropes to be more important than horizontal mixing across surfaces of constant equivalent latitude, shining new light on the interplay between various transport mechanisms in the LMS. The model achieves a good description of the small scale tracer features at the tropopause with squared correlation coefficients R2 = 0.72…0.94.

Nuclear Magnetic Resonance Study of Hydrated Bentonite

Impedance spectroscopy and nuclear magnetic resonance (NMR) were used to investigate the mobility of water molecules located in the interlayer space of H(+) - exchanged bentonite clay. The conductivity obtained by ac measurements was 1.25 x 10(-4) S/cm at 298 K. Proton ((1)H) lineshapes and spin-lattice relaxation times were measured as a function of temperature over the temperature range 130-320 K. The NMR experiments exhibit the qualitative features associated with the proton motion, namely the presence of a (1)H NMR line narrowing and a well-defined spin-lattice relaxation rate maximum. The temperature dependence of the proton spin-lattice relaxation rates was analyzed with the spectral density function appropriate for proton dynamics in a two-dimensional system. The self-diffusion coefficient estimated from our NMR data, D similar to 2 x 10(-7) cm(2)/s at 300 K, is consistent with those reported for exchanged montmorillonite clay hydrates studied by NMR and quasi-elastic neutron scattering (QNS).

Application of Microelectrode Voltammetry to Study the Properties of Surfactant Solutions: Alkyltrimethylammonium Bromides

Microelectrode cyclic voltammetry (MV) has been employed to investigate the micellar properties of solutions of homologous alkyltrimethylammonium bromides, RMe(3)ABr, R = C(10), C(12), and C(14), in water and in the presence of added NaBr. The micellar self-diffusion coefficient was calculated from the limiting current for the reversible electron transfer of micelle-bound ferrocene. From the values of this property, other parameters were calculated, including the micellar hydrodynamic radius, RH, and aggregation number, N(agg); the latter was also theoretically calculated. We determined the values of the diffusion coefficient as a function of various experimental variables and observed the following trends: The diffusion coefficient decreases as a function of increasing surfactant concentration (no additional electrolyte added); it decreases as a function of increasing surfactant concentration at fixed NaBr concentration; and it shows a complex dependence (increase then decrease) on the NaBr concentration at a fixed RMe(3)ABr concentration. The value of the intermicellar interaction parameter decreases and then increases as a function of increasing NaBr concentration. These results are discussed in terms of intermicellar,interactions and the effect of NaBr on the micellar surface charge density and sphere-to-rod geometry change. The NaBr concentration required to induce the latter change increases rapidly as a function of decreasing the length of R: no geometry change was detected for C(10)Me(3)ABr. Values of N(agg) increase as I function of increasing the length of R and are in good agreement with both literature values and values that were calculated theoretically. Thus, MV is a convenient and simple technique for obtaining fundamental properties of surfactant solutions, including additive-induced changes of micellar parameters (N(agg)) and morphology changes.

Alkoxy Chain Effect on the Viscosity of a Quaternary Ammonium Ionic Liquid: Molecular Dynamics Simulations

The viscosity of ionic liquids based on quaternary ammonium cations is reduced when one of the alkyl chains is replaced by an alkoxy chain (Zhou et al. Chem. Eur. J. 2005, 11, 752.). A microscopic picture of the role played by the ether function in decreasing the viscosity of quaternary ammonium ionic liquids is provided here by molecular dynamics (MD) simulations. A model for the ionic liquid N-ethyl-N,N-dimethyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide, MOENM(2)E TFSI, is compared to the tetraalky-lammonium counterpart. The alkoxy derivative has lower viscosity, higher ionic diffusion coefficients, and higher conductivity than the tetraalkyl system at the same density and temperature. A clear signature of the ether function on the liquid structure is observed in cation-cation correlations, but not in anion-anion or anion-cation correlations. In both the alkyl and the alkoxy ionic liquids, there is aggregation of long chains of neighboring cations within micelle-like structures. The MD simulations indicate that the less effective assembly between the more flexible alkoxy chains, in comparison to alkyl chains, is the structural reason for higher ionic mobility in MOENM(2)E TFSI.

Metal transport parameters of a gneiss saprolitic silty soil for liner design

Diffusion coefficients and retardation factors of two metal cations (Cd2+ and Pb2+) were measured for a compacted Brazilian saprolitic soil derived from gneiss, aiming to assess its geoenvironmental performance as a liner for waste disposal sites. This soil occurs extensively all over the country in very thick layers, but has not been used in liners because of its hydraulic conductivity, higher than 10(-9) m/s when compacted at optimum water content of standard Proctor energy, but which can be reduced by means of appropriate compaction techniques or additives. Batch, column, and diffusion tests were carried out with monospecies synthetic solutions at pH 1, 3, and 5.5. Measured diffusion coefficients varied between 0.5 and 4 X 10(-10) m(2)/s. Retardation factors show that cadmium, a very mobile cation, is not adsorbed at pH I but is significantly retained at pH 3 and pH 5.5, whereas lead is retained at all tested pH values though slightly at pH 1. Estimated retardation factors from batch tests were 1.3-2.3 times those resulting from column tests and at its highest when obtained by diffusion tests; whereas batch tests allow a more complete exposure of the soil grains to the solution, time-dependent nonspecific adsorption may take longer to occur. The importance of contact time was observed and should be considered in further investigations. Its significant retention of metals suggests a promising utilization of this soil as a bottom liner for wastes landfills.

On the molecular mass of the extracellular hemoglobin of Glossoscolex paulistus: Analytical ultracentrifugation reexamination

The giant extracellular hemoglobin of Glossoscolex paulistus (HbGp) is constituted by Subunits containing heme groups with molecular masses (M) in the range of 15 to 19 kDa, monomers of 16 kDa (d), and trimers of 51 to 52 kDa (abc) linked by nonheme structures named linkers of 24 to 32 kDa (L). HbGp is homologous to Lumbricus terrestris hemoglobin (HbLt). Several reports propose M of HbLt in the range of 3.6 to 4.4 MDa. Based on subunits M determined by mass spectrometry and assuming HbGp stoichiometry of 12(abcd)(3)L(3) (Vinogradov model) plus 144 heme groups, a Value of M for HbGp oligomer of 3560 kDa can be predicted. This Value is nearly 500 kDa higher than the unique HbGp M Value reported in the literature. In the current work, sedimentation velocity analytical ultracentrifugation (AUC) experiments were performed to obtain M for HbGp in oxy and cyano-met forms. s(20,w)(0), values of 58.1 +/- 0.2 S and 59.6 +/- 0.2 S, respectively, for the two oxidation forms were obtained. The ratio between sedimentation and diffusion coefficients supplied values for M of approximately 3600 100 and 3700 100 kDa for oxy and cyano-met HbGp forms, respectively. An independent determination of the partial specific volume, V(bar), for HbGp was performed based on density measurements, providing a value of 0.764 +/- 0.008, in excellent agreement with the estimates from SEDFIT software. Our results show total consistency between M obtained by AUC and recent partial characterization by mass spectrometry. Therefore, HbGp possesses M very close to that of HbLt, suggesting an oligomeric assembly in agreement with the Vinogradov model. (c) 2008 Elsevier Inc. All rights reserved.