Lattice model of equilibrium polymerization. V. Scattering properties and the width of the critical regime for phase separation

Theory of thermoreversible gelation and anomalous concentration fluctuations in polyzwitterion solutions

We present a theoretical framework to investigate thermoreversible phase transitions within polyzwitterion systems, encompassing macrophase separations (MPS) and gelation. In addition, we explore concentration fluctuations near critical points associated with MPS, as well as tricritical and bicritical points at the intersection of MPS and gelation. By utilizing mean-field percolation theory and field theory formalism, we derive the Landau free energy in terms of polyzwitterion concentration with fixed dipole strengths and other experimental variables, such as temperatures and salt concentrations. As the temperature decreases, the dipoles can form cross-links, resulting in polyzwitterion associations. The associations can grow to a gel network and enhance the propensity for MPS, including liquid-liquid, liquid-gel, and gel-gel phase separations. Remarkably, the associations also impact critical behaviors. Using the renormalization group technique, we find that the critical exponents of the polyzwitterion concentration correlation functions significantly deviate from those in the Ising universality class due to the presence of polyzwitterion associations, leading to crossover critical behaviors.

Applicability of the Generalized Stokes-Einstein Equation of Mode-Coupling Theory to Near-Critical Polyelectrolyte Complex Solutions

We examine whether the mode-coupling theory of Kawasaki and Ferrell (KF) [Kawasaki, K. Kinetic Equations and Time Correlation Functions of Critical Fluctuations. Ann. Phys. 1970, 61 (1), 1-56; Ferrell, R. A. Decoupled-Mode Dynamical Scaling Theory of the Binary-Liquid Phase Transition. Phys. Rev. Lett. 1970, 24 (21), 1169-1172] can describe dynamic light scattering (DLS) measurements of the dynamic structure factor of near-critical polyelectrolyte complex (PC) solutions that have been previously shown to exhibit a theoretically unanticipated lower critical solution temperature type phase behavior, i.e., phase separation upon heating, and a conventional pattern of static critical properties (low angle scattering intensity and static correlation, ξ_s) as a function of reduced temperature. Good qualitative accord is observed between our DLS measurements and the KF theory. In particular, we observe that the collective diffusion coefficient D_c of the PC solutions obeys the generalized Stokes-Einstein equation (GSE), D_c = k_BT/6πηξ_s, where ξ_s is specified from our previous measurements and where η is measured by capillary rheometry under the same thermodynamic conditions as in our previous study of these solutions, allowing for a no-free-parameter test of the GSE. We also find that even the wavevector (q)-dependent collective diffusion coefficient D_c(q), measured by varying the scattering angle in the DLS measurements over a large range, is also well-described by the mean-field version of the KF theory. We find it remarkable that the KF theory provides such a robust description of collective diffusion in these complex charged polyelectrolyte blends under near-critical conditions given that charge fluctuations and association of the polymers might be expected to lead to physical complications that would invalidate the standard model of uncharged fluid mixtures.

Aggregation kinetics of irreversible patches coupled with reversible isotropic interaction leading to chains, bundles and globules

In the present study we are performing simulation of simple model of two patch colloidal particles undergoing irreversible diffusion limited cluster aggregation using patchy Brownian cluster dynamics. In addition to the irreversible aggregation of patches, the spheres are coupled with isotropic reversible aggregation through the Kern–Frenkel potential. Due to the presence of anisotropic and isotropic potential we have also defined three different kinds of clusters formed due to anisotropic potential and isotropic potential only as well as both the potentials together. We have investigated the effect of patch size on self-assembly under different solvent qualities for various volume fractions. We will show that at low volume fractions during aggregation process, we end up in a chain conformation for smaller patch size while in a globular conformation for bigger patch size. We also observed a chain to bundle transformation depending on the attractive interaction strength between the chains or in other words depending on the quality of the solvent. We will also show that bundling process is very similar to nucleation and growth phenomena observed in colloidal system with short range attraction. We have also studied the bond angle distribution for this system, where for small patches only two angles are more probable indicating chain formation, while for bundling at very low volume fraction a tail is developed in the distribution. While for the case of higher patch angle this distribution is broad compared to the case of low patch angles showing we have a more globular conformation. We are also proposing a model for the formation of bundles which are similar to amyloid fibers using two patch colloidal particles.

Mixtures of two self- and mutually-associating liquids: Phase behavior, second virial coefficients, and entropy-enthalpy compensation in the free energy of mixing

The theoretical description of the phase behavior of polymers dissolved in binary mixtures of water and other miscible solvents is greatly complicated by the self- and mutual-association of the solvent molecules. As a first step in treating these complex and widely encountered solutions, we have developed an extension of Flory-Huggins theory to describe mixtures of two self- and mutually-associating fluids comprised of small molecules. Analytic expressions are derived here for basic thermodynamic properties of these fluid mixtures, including the spinodal phase boundaries, the second osmotic virial coefficients, and the enthalpy and entropy of mixing these associating solvents. Mixtures of this kind are found to exhibit characteristic closed loop phase boundaries and entropy-enthalpy compensation for the free energy of mixing in the low temperature regime where the liquid components are miscible. As discussed by Widom et al. [Phys. Chem. Chem. Phys. 5, 3085 (2003)], these basic miscibility trends, quite distinct from those observed in non-associating solvents, are defining phenomenological characteristics of the "hydrophobic effect." We find that our theory of mixtures of associating fluids captures at least some of the thermodynamic features of real aqueous mixtures.

Coupling of isotropic and directional interactions and its effect on phase separation and self-assembly

The interactions of molecules and particles in solution often involve an interplay between isotropic and highly directional interactions that lead to a mutual coupling of phase separation and self-assembly. This situation arises, for example, in proteins interacting through hydrophobic and charged patch regions on their surface and in nanoparticles with grafted polymer chains, such as DNA. As a minimal model of complex fluids exhibiting this interaction coupling, we investigate spherical particles having an isotropic interaction and a constellation of five attractive patches on the particle's surface. Monte Carlo simulations and mean-field calculations of the phase boundaries of this model depend strongly on the relative strength of the isotropic and patch potentials, where we surprisingly find that analytic mean-field predictions become increasingly accurate as the directional interactions become increasingly predominant. We quantitatively account for this effect by noting that the effective interaction range increases with increasing relative directional to isotropic interaction strength. We also identify thermodynamic transition lines associated with self-assembly, extract the entropy and energy of association, and characterize the resulting cluster properties obtained from simulations using percolation scaling theory and Flory-Stockmayer mean-field theory. We find that the fractal dimension and cluster size distribution are consistent with those of lattice animals, i.e., randomly branched polymers swollen by excluded volume interactions. We also identify a universal functional form for the average molecular weight and a nearly universal functional form for a scaling parameter characterizing the cluster size distribution. Since the formation of branched clusters at equilibrium is a common phenomenon in nature, we detail how our analysis can be used in experimental characterization of such associating fluids.

Phase field method for nonequilibrium dynamics of reversible self-assembly systems

Phase field methods are extended to describe the nonequilibrium dynamics of reversible self-assembly systems, an extension that is complicated by the mutual coupling of many non-conserved order parameters into a set of highly nonlinear partial differential equations. Further complications arise because the sum of all non-conserved order parameters equals a conserved order parameter. The theory is developed for the simplest model of reversible self-assembly in which no additional constraints are imposed on the self-assembly process since the extension to treat more complex self-assembly models is straightforward. Specific calculations focus on the time evolution of the cluster size distribution for a free association system that is rapidly dropped from one ordered state to a more ordered state within the one-phase region. The dynamics proceed as expected, thereby providing validation of the theory which is also capable of treating systems with spatial inhomogeneities.

Lattice cluster theory of associating telechelic polymers. III. Order parameter and average degree of self-assembly, transition temperature, and specific heat

The lattice cluster theory of strongly interacting, structured polymer fluids is applied to determine the thermodynamic properties of solutions of telechelic polymers that may associate through bifunctional end groups. Hence, this model represents a significant albeit natural extension of a diverse array of prior popular equilibrium polymerization models in which structureless "bead" monomers associate into chain-like clusters under equilibrium conditions. In particular, the thermodynamic description of the self-assembly of linear telechelic chains in small molecule solvents (initiated in Paper II) is systematically extended through calculations of the order parameter Φ and average degree <N> of self-assembly, the self-assembly transition temperature T(p), and the specific heat C(V) of solutions of telechelic molecules. Special focus is placed on examining how molecular and thermodynamic parameters, such as the solution composition φ, temperature T, microscopic interaction energies (ε(s) and ε), and length M of individual telechelic chains, influence the computed thermodynamic quantities that are commonly used to characterize self-assembling systems.

Equilibrium polymerization models of re-entrant self-assembly

As is well known, liquid-liquid phase separation can occur either upon heating or cooling, corresponding to lower and upper critical solution phase boundaries, respectively. Likewise, self-assembly transitions from a monomeric state to an organized polymeric state can proceed either upon increasing or decreasing temperature, and the concentration dependent ordering temperature is correspondingly called the "floor" or "ceiling" temperature. Motivated by the fact that some phase separating systems exhibit closed loop phase boundaries with two critical points, the present paper analyzes self-assembly analogs of re-entrant phase separation, i.e., re-entrant self-assembly. In particular, re-entrant self-assembly transitions are demonstrated to arise in thermally activated equilibrium self-assembling systems, when thermal activation is more favorable than chain propagation, and in equilibrium self-assembly near an adsorbing boundary where strong competition exists between adsorption and self-assembly. Apparently, the competition between interactions or equilibria generally underlies re-entrant behavior in both liquid-liquid phase separation and self-assembly transitions.

Critical behavior on approaching a special critical point in a complex fluid

The critical behavior of osmotic susceptibility is investigated in the re-entrant complex mixture 1-propanol (P)+water (W)+potassium chloride (KCl) through light-scattering measurements. The measurements are performed on approaching a special critical point [i.e., the double critical point (DCP)] in this mixture, along the line of upper critical solution temperatures (T(U)'s), by varying t from the high temperature one-phase region. The light-scattering data analysis emphasizes the need for corrections to the asymptotic Ising behavior and yields very large magnitudes for the correction-to-scaling amplitudes A(1) and A(2), with the first-correction amplitude A(1) being negative, signifying a nonmonotonic crossover behavior of the susceptibility exponent in this mixture. For the T(U) closest to the DCP, the effective susceptibility exponent gamma(eff) displays a nonmonotonic crossover from its nearly doubled three dimensional (3D)-Ising value toward its nearly double mean-field value with an increase in t. While for that far away from the DCP, it displays a nonmonotonic crossover from its single-limit Ising value toward a value slightly lower than its mean-field value of 1 with an increase in t. This feature of the effective susceptibility exponent is interpreted in terms of the possibility of a nonmonotonic crossover to the mean-field value from lower values in the nonasymptotic high t region. The renormalized Ising regime extends over a larger t range for the sample (or T(U)) closest to the DCP when compared to that far away from it. The in-between T(U)'s display a trend toward shrinkage in the renormalized Ising regime as T(U) shifts away from the DCP. Nevertheless, the crossover to the mean-field behavior is completed only beyond t>10(-2) for the T(U)'s studied. The observed crossover behavior is attributed to the presence of strong ion-induced clustering in this mixture, as revealed by various structure probing techniques, while the observed unique trend in the crossover behavior is discussed in terms of the varying influence of the DCP on the critical behavior along the T(U) line. The crossover behavior for the T(U)'s is pronounced and more sharp compared to the T(L)'s (lower critical solution temperatures) [U. K. Pradeep, J. Chem. Phys. 129, 134506 (2008)] in this mixture, although there exists no difference in the growth of the mesoscale clusters in the lower and upper one-phase regions in this mixture. Our observations suggest the need to look at the crossover behavior probably from two perspectives, namely, the dielectric effect and the clustering effect. The effective susceptibility exponent as a function of the field variable t(UL), instead of the conventional variable t, displays a sharp nonmonotonic crossover from its asymptotic 3D-Ising value ( approximately 1.24) toward a value slightly lower than its nonasymptotic mean-field value of 1, as that observed in the t analysis for the T(U) far away from the influence of the DCP.

Lattice model of equilibrium polymerization. VII. Understanding the role of "cooperativity" in self-assembly

Cooperativity is an emergent many-body phenomenon related to the degree to which elementary entities (particles, molecules, organisms) collectively interact to form larger scale structures. From the standpoint of a formal mean field description of chemical reactions, the cooperativity index m, describing the number of elements involved in this structural self-organization, is the order of the reaction. Thus, m for molecular self-assembly is the number of molecules in the final organized structure, e.g., spherical micelles. Although cooperativity is crucial for regulating the thermodynamics and dynamics of self-assembly, there is a limited understanding of this aspect of self-assembly. We analyze the cooperativity by calculating essential thermodynamic properties of the classical mth order reaction model of self-assembly (FAm model), including universal scaling functions describing the temperature and concentration dependence of the order parameter and average cluster size. The competition between self-assembly and phase separation is also described. We demonstrate that a sequential model of thermally activated equilibrium polymerization can quantitatively be related to the FAm model. Our analysis indicates that the essential requirement for "cooperative" self-assembly is the introduction of constraints (often nonlocal) acting on the individual assembly events to regulate the thermodynamic free energy landscape and, thus, the thermodynamic sharpness of the assembly transition. An effective value of m is defined for general self-assembly transitions, and we find a general tendency for self-assembly to become a true phase transition as m-->infinity. Finally, various quantitative measures of self-assembly cooperativity are discussed in order to identify experimental signatures of cooperativity in self-assembling systems and to provide a reliable metric for the degree of transition cooperativity.

Does equilibrium polymerization describe the dynamic heterogeneity of glass-forming liquids?

A significant body of evidence indicates that particles with excessively high or low mobility relative to Brownian particles form in dynamic equilibrium in glass-forming liquids. We examine whether these "dynamic heterogeneities" can be identified with a kind of equilibrium polymerization. This correspondence is first checked by demonstrating the presence of a striking resemblance between the temperature dependences of the configurational entropy s(c) in both the theory of equilibrium polymerization and the generalized entropy theory of glass formation in polymer melts. Moreover, the multiple characteristic temperatures of glass formation are also shown to have analogs in the thermodynamics of equilibrium polymerization, supporting the contention that both processes are varieties of rounded thermodynamic transitions. We also find that the average cluster mass (or degree of polymerization) varies in nearly inverse proportionality to s(c). This inverse relation accords with the basic hypothesis of Adam-Gibbs that the number of particles in the cooperatively rearranging regions (CRR) of glass-forming liquids scales inversely to s(c) of the fluid. Our identification of the CRR with equilibrium polymers is further supported by simulations for a variety of glass-forming liquids that verify the existence of stringlike or polymeric clusters exhibiting collective particle motion. Moreover, these dynamical clusters have an exponential length distribution, and the average "string" length grows upon cooling according to the predictions of equilibrium polymerization theory. The observed scale of dynamic heterogeneity in glass-forming liquids is found to be consistent with this type of self-assembly process. Both experiments and simulations have revealed remarkable similarities between the dynamical properties of self-assembling and glass-forming liquids, suggesting that the development of a theory for the dynamics of self-assembling fluids will also enhance our understanding of relaxation in glass-forming liquids.

Simulation and theory of flexible equilibrium polymers under poor solvent conditions

Grand canonical Monte Carlo simulation and simple statistical thermodynamic theory are used to model the aggregation and phase separation of systems of reversibly polymerizing monomers, capable of forming chains with or without the ability to cyclize into rings, with isotropic square-well attractions between nonbonded pairs of monomers. The general trend observed in simulation of chain-only systems, as predicted in a number of published theoretical works, is that the critical temperature for phase separation increases and the critical monomer density decreases with rising polymer bond strength. Introduction of the equilibrium between chains and rings into the theory lowers the predicted critical temperature and increases the predicted critical density. While the chain-only theories predict a vanishing critical density in the limit of complete polymerization, when ring formation is taken into account the predicted critical density in the same limit approaches the density of the onset of the ring-chain transition. The theoretically predicted effect of cyclization on chemical potential is in good qualitative agreement with a subset of simulation results in which chain-only systems were compared with equilibrium mixtures of rings and chains. The influence of attractions on the aggregation number and radius of gyration of chains and rings observed in simulations is also discussed.