Quantitative description of the self-assembly of patchy particles into chains and rings

We numerically study a simple fluid composed of particles having a hard-core repulsion complemented by two patchy attractive sites on the particle poles. An appropriate choice of the patch angular width allows for the formation of ring structures which, at low temperatures and low densities, compete with the growth of linear aggregates. The simplicity of the model makes it possible to compare simulation results and theoretical predictions based on the Wertheim perturbation theory, specialized to the case in which ring formation is allowed. Such a comparison offers a unique framework for establishing the quality of the analytic predictions. We find that the Wertheim theory describes remarkably well the simulation results.

Scaling behavior and equilibrium lengths of knotted polymers

We use numerical simulations to investigate how the chain length and topology of freely fluctuating knotted polymer rings affect their various spatial characteristics such as the radius of the smallest sphere enclosing momentary configurations of simulated polymer chains. We describe how the average value of a characteristic changes with the chain size and how this change depends on the topology of the modeled polymers. Although the scaling profiles of a spatial characteristic for distinct knot types do not intersect (at least, in the range of our data), the profiles for nontrivial knots intersect the corresponding profile obtained for phantom polymers, i.e., those that are free to explore all available topological states. For each knot type, this point of intersection defines its equilibrium length with respect to the spatial characteristic. At this chain length, a polymer forming a given knot type will not tend to increase or decrease. on average, the value of the spatial characteristic when the polymer is released from its topological constraint. We show interrelations between equilibrium lengths defined with respect to spatial characteristics of different character and observe that they are related to the lengths of ideal geometric configurations of the corresponding knot types.

Some examples of modeling in chemical engineering: thermal treatment of polymers, phase equilibrium calculations and process design

Presentation submitted to PSE Seminar, Chemical Engineering Department, Center for Advanced Process Design-making (CAPD), Carnegie Mellon University, Pittsburgh (USA), October 2012.

Equilibrium swelling measurements of network and semicrystalline polymers in supercritical carbon dioxide using high-pressure NMR

The extent of swelling of cross-linked poly(dimethylsiloxane) and linear low-density poly(ethylene) in supercritical CO2 has been investigated using high-pressure NMR spectroscopy and microscopy. Poly(dimethylsiloxane) was cross-linked to four different cross-link densities and swollen in supercritical CO2. The Flory-Huggins interaction parameter, x, was found to be 0.62 at 300 bar and 45 degrees C, indicating that supercritical CO2 is a relatively poor solvent compared to toluene or benzene. Linear low-density poly(ethylene) was shown to exhibit negligible swelling upon exposure to supercritical CO2 up to 300 bar. The effect Of CO2 pressure on the amorphous region of the poly(ethylene) was investigated by observing changes in the H-1 T-2 relaxation times of the polymer. These relaxation times decreased with increasing pressure, which was attributed to a decrease in mobility of the polymer chains as a result of compressive pressure.

AN ASSESSMENT ON THE USE OF THE DEBYE-HUCKEL EQUATION FOR THE THERMODYNAMIC MODELING OF AQUEOUS SYSTEMS CONTAINING POLYMERS AND SALTS

In this work, a study on the role of the long-range term of excess Gibbs energy models in the modeling of aqueous systems containing polymers and salts is presented. Four different approaches on how to account for the presence of polymer in the long-range term were considered, and simulations were conducted considering aqueous solutions of three different salts. The analysis of water activity curves showed that, in all cases, a liquid-phase separation may be introduced by the sole presence of the polymer in the long-range term, regardless of how it is taken into account. The results lead to the conclusion that there is no single exact solution for this problem, and that any kind of approach may introduce inconsistencies.

Liquid-Liquid Equilibrium of Aqueous Two-Phase Systems Containing Some Synthetic Polyelectrolytes and Polyethylene Glycol

Experimental results are presented for the liquid-liquid equilibrium of aqueous two-phase systems containing a synthetic polyelectrolyte (polysodium acrylate, polysodium methacrylate, and polysodium ethylene sulfonate) and polyethylene glycol at (298.2 and 323.2) K. A total of 40 phase diagrams were obtained, comprising data both of the binodal curve (obtained through cloud-point measurements) and of equilibrium compositions. The influences of temperature, the nature of the polyelectrolyte monomer unit, and the chain length of both types of polymers are analyzed and discussed.

H-1 NMR study of the states of water in equilibrium Poly(HEMA-co-THFMA) hydrogels

The spin-spin relaxation times, T-2, of hydrated samples of poly(hydroxymethyl methacrylate), PHEMA, poly(tetrahydrofurfuryl methacrylate),PTHFMA, and the,corresponding HEMA-THFMA copolymers have been examined to probe the states of,the imbibed water in these polymers. The decay in the transverse magnetization of water. in fully hydrated samples of PHEMA, PTHFMA, and copolymers of HEMA and THFMA was described by a multiexponential function. The short component of T-2 was interpreted as water molecules that were strongly interacting with the polymer chains. The intermediate component of T-2 was assigned to water residing in the porous structure of the samples. The long component of T-2 was believed to arise from water residing in the remnants of cracks formed in the polymer network during water sorption.

PFG-NMR measurements of the self-diffusion coefficients of water in equilibrium poly(HEMA-co-THFMA) hydrogels

The self-diffusion coefficients for water in a series of copolymers of 2-hydroxyethyl methacrylate, HEMA, and tetrahydrofurfuryl methacrylate, THFMA, swollen with water to their equilibrium states have been studied at 310 K using PFG-NMR. The self-diffusion coefficients calculated from the Stejskal-Tanner equation, D-obs, for all of the hydrated polymers were found to be dependent on the NMR storage time, as a result of spin exchange between the proton reservoirs of the water and the polymers, reaching an equilibrium plateau value at long storage times. The true values of the diffusion coefficients were calculated from the values of D-obs, in the plateau regions by applying a correction for the fraction of water protons present, obtained from the equilibrium water contents of the gels. The true self-diffusion coefficient for water in polyHEMA obtained at 310 K by this method was 5.5 x 10(-10) m(2) s(-1). For the copolymers containing 20% HEMA or more a single value of the self-diffusion coefficient was found, which was somewhat larger than the corresponding values obtained for the macroscopic diffusion coefficient from sorption measurements. For polyTHFMA and copolymers containing less than 20% HEMA, the PFG-NMR stimulated echo attenuation decay curves and the log-attenuation plots were characteristic of the presence of two diffusing water species. The self-diffusion coefficients of water in the equilibrium-hydrated copolymers were found to be dependent on the copolymer composition, decreasing with increasing THFMA content.

Equilibrium self-assembly of colloids with distinct interaction sites: Thermodynamics, percolation, and cluster distribution functions

We calculate the equilibrium thermodynamic properties, percolation threshold, and cluster distribution functions for a model of associating colloids, which consists of hard spherical particles having on their surfaces three short-ranged attractive sites (sticky spots) of two different types, A and B. The thermodynamic properties are calculated using Wertheim's perturbation theory of associating fluids. This also allows us to find the onset of self-assembly, which can be quantified by the maxima of the specific heat at constant volume. The percolation threshold is derived, under the no-loop assumption, for the correlated bond model: In all cases it is two percolated phases that become identical at a critical point, when one exists. Finally, the cluster size distributions are calculated by mapping the model onto an effective model, characterized by a-state-dependent-functionality (f) over bar and unique bonding probability (p) over bar. The mapping is based on the asymptotic limit of the cluster distributions functions of the generic model and the effective parameters are defined through the requirement that the equilibrium cluster distributions of the true and effective models have the same number-averaged and weight-averaged sizes at all densities and temperatures. We also study the model numerically in the case where BB interactions are missing. In this limit, AB bonds either provide branching between A-chains (Y-junctions) if epsilon(AB)/epsilon(AA) is small, or drive the formation of a hyperbranched polymer if epsilon(AB)/epsilon(AA) is large. We find that the theoretical predictions describe quite accurately the numerical data, especially in the region where Y-junctions are present. There is fairly good agreement between theoretical and numerical results both for the thermodynamic (number of bonds and phase coexistence) and the connectivity properties of the model (cluster size distributions and percolation locus).

Criticality of colloids with distinct interaction patches: The limits of linear chains, hyperbranched polymers, and dimers

We use a simple model of associating fluids which consists of spherical particles having a hard-core repulsion, complemented by three short-ranged attractive sites on the surface (sticky spots). Two of the spots are of type A and one is of type B; the bonding interactions between each pair of spots have strengths epsilon(AA), epsilon(BB), and epsilon(AB). The theory is applied over the whole range of bonding strengths and the results are interpreted in terms of the equilibrium cluster structures of the phases. In addition to our numerical results, we derive asymptotic expansions for the free energy in the limits for which there is no liquid-vapor critical point: linear chains (epsilon(AA)not equal 0, epsilon(AB)=epsilon(BB)=0), hyperbranched polymers (epsilon(AB)not equal 0, epsilon(AA)=epsilon(BB)=0), and dimers (epsilon(BB)not equal 0, epsilon(AA)=epsilon(AB)=0). These expansions also allow us to calculate the structure of the critical fluid by perturbing around the above limits, yielding three different types of condensation: of linear chains (AA clusters connected by a few AB or BB bonds); of hyperbranched polymers (AB clusters connected by AA bonds); or of dimers (BB clusters connected by AA bonds). Interestingly, there is no critical point when epsilon(AA) vanishes despite the fact that AA bonds alone cannot drive condensation.

Effect of Knotting on the Shape of Polymers

Momentary configurations of long polymers at thermal equilibrium usually deviate from spherical symmetry and can be better described, on average, by a prolate ellipsoid. The asphericity and nature of asphericity (or prolateness) that describe these momentary ellipsoidal shapes of a polymer are determined by specific expressions involving the three principal moments of inertia calculated for configurations of the polymer. Earlier theoretical studies and numerical simulations have established that as the length of the polymer increases, the average shape for the statistical ensemble of random configurations asymptotically approaches a characteristic universal shape that depends on the solvent quality. It has been established, however, that these universal shapes differ for linear, circular, and branched chains. We investigate here the effect of knotting on the shape of cyclic polymers modeled as random isosegmental polygons. We observe that random polygons forming different knot types reach asymptotic shapes that are distinct from the ensemble average shape. For the same chain length, more complex knots are, on average, more spherical than less complex knots.

Total curvature and total torsion of knotted polymers

Previous work on radius of gyration and average crossing number has demonstrated that polymers with fixed topology show a different scaling behavior with respect to these characteristics than polymers with unrestricted topology. Using numerical simulations, we show here that the difference in the scaling behavior between polymers with restricted and unrestricted topology also applies to the total curvature and total torsion. For each knot type, the equilibrium length with respect to a given spatial characteristic is the number of edges at which the value of the characteristic is the same as the average for all polygons. This number appears to be correlated to physical properties of macromolecules, for example gel mobility as measured by the separation between distinct knot types. We also find that, on average, closed polymers require slightly more total curvature and slightly less total torsion than open polymers with the corresponding number of monomers.

The effects of additives and chemical modification on the solution properties of thermo-sensitive polymers

Cette thèse concerne l’étude de phase de séparation de deux polymères thermosensibles connus-poly(N-isopropylacylamide) (PNIPAM) et poly(2-isopropyl-2-oxazoline) (PIPOZ). Parmi des études variées sur ces deux polymères, il y a encore deux parties de leurs propriétés thermiques inexplicites à être étudiées. Une partie concerne l’effet de consolvant de PNIPAM dans l’eau et un autre solvant hydromiscible. L’autre est l’effet de propriétés de groupes terminaux de chaînes sur la séparation de phase de PIPOZ. Pour ce faire, nous avons d’abord étudié l’effet de l’architecture de chaînes sur l’effet de cosolvant de PNIPAMs dans le mélange de méthanol/eau en utilisant un PNIPAM en étoile avec 4 branches et un PNIPAM cyclique comme modèles. Avec PNIPAM en étoile, l’adhérence de branches PNIPAM de à un cœur hydrophobique provoque une réduction de Tc (la température du point de turbidité) et une enthalpie plus faible de la transition de phase. En revanche, la Tc de PNIPAM en étoile dépend de la masse molaire de polymère. La coopérativité de déhydratation diminue pour PNIPAM en étoile et PNIPAM cyclique à cause de la limite topologique. Une étude sur l’influence de concentration en polymère sur l’effet de cosolvant de PNIPAM dans le mélange méthanol/eau a montré qu’une séparation de phase liquide-liquide macroscopique (MLLPS) a lieu pour une solution de PNIPAM dans le mélange méthanol/eau avec la fraction molaire de méthanol entre 0.127 et 0.421 et la concentration en PNIPAM est constante à 10 g.L-1. Après deux jours d’équilibration à température ambiante, la suspension turbide de PNIPAM dans le mélange méthanol/eau se sépare en deux phases dont une phase possède beaucoup plus de PNIPAM que l’autre. Un diagramme de phase qui montre la MLLPS pour le mélange PNIPAM/eau/méthanol a été établi à base de données expérimentales. La taille et la morphologie de gouttelettes dans la phase riche en polymère condensée dépendent de la fraction molaire de méthanol. Parce que la présence de méthanol influence la tension de surface des gouttelettes liquides, un équilibre lent de la séparation de phase pour PNIPAM/eau/méthanol système a été accéléré et une séparation de phase liquide-liquide macroscopique apparait. Afin d’étudier l’effet de groupes terminaux sur les propriétés de solution de PIPOZ, deux PIPOZs téléchéliques avec groupe perfluorodécanyle (FPIPOZ) ou groupe octadécyle (C18PIPOZ) comme extrémités de chaîne ont été synthétisés. Les valeurs de Tc des polymères téléchéliques ont beaucoup diminué par rapport à celle de PIPOZ. Des micelles stables se forment dans des solutions aqueuses de polymères téléchéliques. La micellization et la séparation de phase de ces polymères dans l’eau ont été étudiées. La séparation de phase de PIPOZs téléchéliques suit le mécanisme de MLLPS. Des différences en tailles de gouttelettes formées à l’intérieur de solutions de deux polymères ont été observées. Pour étudier profondément les différences dans le comportement d’association entre deux polymères téléchéliques, les intensités des signaux de polymères correspondants et les temps de relaxation T1, T2 ont été mesurés. Des valeurs de T2 de protons correspondants aux IPOZs sont plus hautes.

Microwave-assisted hydrogel synthesis: a new method for crosslinking polymers in aqueous solutions

It has been found that hydrogels may be formed by microwave irradiation of aqueous solutions containing appropriate combinations of polymers. This new method of hydrogel synthesis yields sterile hydrogels without the use of monomers, eliminating the need for the removal of unreacted species from the final product. Results for two particularly successful combinations, poly(vinyl alcohol) with either poly(acrylic acid) or poly(methylvinylether-alt-maleic anhydride), are presented. Irradiation using temperatures of 100–150 °C was found to yield hydrogels with large equilibrium swelling degrees of 500–1000 g g−1. Material leached from both types of hydrogel shows little cytotoxicity towards HT29 cells.