Statistical thermodynamics in the framework of the lattice fluid model, 2. Binary mixture of polydisperse polymers of special distribution

In this paper, the Gibbs free energy, the equation of state and the chemical potentials of polydisperse multicomponent polymer mixtures are derived. For general binary mixtures of polydisperse polymers, we also give the Gibbs free energy, the equation of state and the chemical potentials and derive the stability criteria and spinodal. Furthermore, binary polydisperse polymer mixtures of special distribution, i.e., Flory distribution, uniform distribution and Schulz distribution, are discussed and the influence of polydispersity on the interaction energy parameter is considered. For these special-distribution systems, the spinodal curves are simulated and the influence of chain length and polydispersity on the spinodal curves is discussed. The results suggest that the spinodal temperature of the mixture with a given volume fraction of one component decreases with increasing polydispersity and the extent of the shift decreases with increasing degree of polymerization when η = M̄_w/M̄_n is given. In addition, the variations of the spinodal curves with polydispersity and chain length are shown and they are qualitatively compared with the experimental results.     

Statistical thermodynamics of polydisperse polymer blends with strong interaction

A statistical thermodynamics theory of polydisperse polymer mixtures with strong interaction between dissimilar components based on a lattice fluid model is formulated. Expressions for the free energy, equation of state, phase stability and spinodal for a polydisperse, binary polymer mixture with strong interaction are derived.     

Modeling the flow of fluid/particle mixtures in microchannels: encapsulating nanoparticles within monodisperse droplets

We develop a hybrid computational approach for simulating mixtures of binary fluids and mobile, submicron particles. The model couples a lattice Boltzmann method for the binary fluid with a Brownian dynamics model for the particles. The particles can exhibit preferential wetting interactions with the different components of the fluid. As an illustration of the method, we carry out simulations in two dimensions to compare the spinodal decomposition of a pure binary fluid with the phase separation of binary blends that contain either fixed or mobile particles. We then isolate conditions where the flow of a binary fluid/particle mixture past surfaces with well-defined asperities leads to the formation of monodisperse droplets, which encapsulate the nanoparticles. The findings provide guidelines for creating multiphase emulsions with well-controlled morphologies.     

Calculation of the stability and of the phase equilibrium in the system polystyrene+cyclohexane+carbon dioxide based on equations of state

In order to calculate spinodals for polymer systems with an equation of state (EOS), we developed a stability theory using continuous thermodynamics. Here, the mixture considered consists of a polydisperse polymer and two monodisperse components as for example a solvent and a gas. We derived the spinodal equation on the base of the segment-molar Helmholtz energy of the mixture. As a result, a determinant similar to that of the monodisperse case is obtained, but the polydisperse polymer is identified by its weight average of the molecular weight. Furthermore, our paper provides the equations for the cloud-point curve derived with the aid of continuous thermodynamics. The final equations are applied to the system polystyrene+cyclohexane+carbon dioxide using the EOS of Sako, Wu and Prausnitz (SWP-EOS). For parameter fit and to prove the accuracy of the treatment, experimental data of the high-pressure equilibrium of the binary subsystems and of the ternary system were taken from literature.     

Phase coexistence in a polydisperse charged hard-sphere fluid: polymer mean spherical approximation

We have reconsidered the phase behavior of a polydisperse mixture of charged hard spheres (CHSs) introducing the concept of minimal size neutral clusters. We thus take into account ionic association effects observed in charged systems close to the phase boundary where the properties of the system are dominated by the presence of neutral clusters while the amount of free ions or charged clusters is negligible. With this concept we clearly pass beyond the simple level of the mean spherical approximation (MSA) that we have presented in our recent study of a polydisperse mixture of CHS [Yu. V. Kalyuzhnyi, G. Kahl, and P. T. Cummings, J. Chem. Phys. 120, 10133 (2004)]. Restricting ourselves to a 1:1 and possibly size-asymmetric model we treat the resulting polydisperse mixture of neutral, polar dimers within the framework of the polymer MSA, i.e., a concept that--similar as the MSA--readily can be generalized from the case of a mixture with a finite number of components to the polydisperse case: again, the model belongs to the class of truncatable free-energy models so that we can map the formally infinitely many coexistence equations onto a finite set of coupled, nonlinear equations in the generalized moments of the distribution function that characterizes the system. This allows us to determine the full phase diagram (in terms of binodals as well as cloud and shadow curves), we can study fractionation effects on the level of the distribution functions of the coexisting daughter phases, and we propose estimates on how the location of the critical point might vary in a polydisperse mixture with an increasing size asymmetry and polydispersity.     

Phase behavior of weakly polydisperse sticky hard spheres: perturbation theory for the Percus-Yevick solution

We study the effects of size polydispersity on the gas-liquid phase behavior of mixtures of sticky hard spheres. To achieve this, the system of coupled quadratic equations for the contact values of the partial cavity functions of the Percus-Yevick solution [R. J. Baxter, J. Chem. Phys. 49, 2770 (1968)] is solved within a perturbation expansion in the polydispersity, i.e., the normalized width of the size distribution. This allows us to make predictions for various thermodynamic quantities which can be tested against numerical simulations and experiments. In particular, we determine the leading order effects of size polydispersity on the cloud curve delimiting the region of two-phase coexistence and on the associated shadow curve; we also study the extent of size fractionation between the coexisting phases. Different choices for the size dependence of the adhesion strengths are examined carefully; the Asakura-Oosawa model [J. Chem. Phys. 22, 1255 (1954)] of a mixture of polydisperse colloids and small polymers is studied as a specific example.     

Investigation of particle interactions in concentrated colloidal suspensions using frequency domain photon migration: monodisperse systems

Frequency domain photon migration (FDPM) measurements were conducted to assess particle interactions of concentrated, monodisperse, polystyrene samples obtained directly from industry by using multiple scattering light. The angle-integrated static structure factor, S(q), and static structure factor at small wave vector q, S(0), were obtained from FDPM measurements at high volume fractions ranging from 0.05 to 0.3, and were compared with those obtained from the monodisperse hard sphere Percus-Yevick (HSPY) model. Effects of different colloid sizes on structure factor evaluated at two different wavelengths were also investigated. Results show that the monodisperse HSPY model is suitable for accounting for particle interactions and local microstructures in these colloidal suspensions. Upon using the HSPY model, particle sizes of polystyrene suspensions were recovered from FDPM measurements at high volume fractions (up to 0.3), which agree well with the DLS measurement of diluted sample ( approximately 0.001). The study of polydispersity effect shows that the FDPM method can be successfully used for recovering the mean particle size of polydisperse colloidal suspension with low polydispersity (<16%) under the assumption of monodisperse hard sphere systems.     

Scattering function and the dynamics of phase separation in polymer mixtures under shear flow

The phenomenological mean-field theory describing concentration fluctuations and spinodal decomposition of binary mixtures of long flexible macromolecules is generalized to mixtures under steady shear flow. This shear flow leads to a partial orientation and stretching of the coils, as well as to an anisotropic deformation of concentration fluctuations. Generalizing the approach of Onuki and Kawasaki, we obtain the collective scattering function describing these concentration fluctuations in the mixture under shear flow. Both the steady-state situation in the one-phase region and the initial stages of spinodal decomposition for concentrations inside of the spinodal curve are considered.Contributed paper delivered at the Tagung der Deutschen Physikalischen Gesellschaft, Fachausschuß Polymerphysik, Berlin, March 30–April 3, 1987.     

Multiphase coexistence in polydisperse colloidal mixtures

The authors study the phase behavior of mixtures of monodisperse colloidal spheres with a depletion agent which can have arbitrary shape and can possess a polydisperse size or shape distribution. In the low concentration limit considered here, the authors can employ the free-volume theory and take the geometry of particles of the depletion agent into account within the framework of fundamental measure theory. The authors apply their approach to study the phase diagram of a mixture of (monodisperse) colloidal spheres and two polydisperse polymer components. By fine tuning the distribution of the polymer, it is possible to construct a complex phase diagram which exhibits two stable critical points.     

Phase coexistence in polydisperse athermal polymer-colloidal mixture

A theoretical scheme developed earlier [Y. V. Kalyuzhnyi et al., Chem. Phys. Lett. 443, 243 (2007)] is used to calculate the full phase diagram of polydisperse athermal polymer-colloidal mixture with polydispersity in both colloidal and polymeric components. In the limiting case of bidisperse polymer-colloidal mixture, theoretical results are compared against computer simulation results. We present the cloud and shadow curves, critical binodals, and distribution functions of the coexisting phases and discuss the effects of polydispersity on their behavior. According to our analysis polydispersity extends the region of the phase instability, shifting the critical point to the lower values of the pressure and density. For the high values of the pressure polydispersity causes strong fractionation effects, with the large size colloidal particles preferring the low-density shadow phase and long chain length polymeric particles preferring the high-density shadow phase.     

Size fractionation in a phase-separated colloidal fluid

Phase separation of a polydisperse colloidal dispersion implies size fractionation. An application of this effect is given by size-selective purification procedures associated with the colloidal synthesis of so-called monodisperse nanoparticles. We used electron microscopy to determine detailed particle size distributions of coexisting colloidal fluid phases containing highly polydisperse iron oxide nanoparticles with a log-normal distribution (sigma = 0.54 for the total system). Analysis of N approximately 10000 particles per phase yields the first five statistical moments of the distributions. Within experimental error, the interdependence of the statistical moments is in quantitative agreement with the "universal law of fractionation" proposed by Evans, Fairhurst, and Poon [Phys. Rev. Lett. 1998, 81, 1326], even though the theory was derived in the limit of slight polydispersity.     

Phase coexistence in polydisperse multi-Yukawa hard-sphere fluid: high temperature approximation

High temperature approximation (HTA) is used to describe the phase behavior of polydisperse multi-Yukawa hard-sphere fluid mixtures. It is demonstrated that in the frames of the HTA the model belongs to the class of "truncatable free energy models," i.e., the models with thermodynamical properties (Helmholtz free energy, chemical potential, and pressure) defined by the finite number of generalized moments. Using this property we were able to calculate the complete phase diagram (i.e., cloud and shadow curves as well as binodals) and size distribution functions of the coexisting phases of several different models of polydisperse fluids. In particular, we consider polydisperse one-Yukawa hard-sphere mixture with factorizable Yukawa coefficients and polydisperse Lennard-Jones (LJ) mixture with interaction energy parameter and/or size polydispersity. To validate the accuracy of the HTA we compare theoretical results with previously published results of more advanced mean spherical approximation (MSA) for the one-Yukawa model and with the Monte Carlo (MC) computer simulation results of [Wilding et al. J. Chem. Phys. 121, 6887 (2004); Phys. Rev. Lett. 95, 155701 (2005)] for the LJ model. We find that overall predictions of the HTA are in reasonable agreement with predictions of the MSA and MC, with the accuracy range from semiquantitative (for the phase diagram) to quantitative (for the size distribution functions).     

Spinodal decomposition in asymmetric polymer mixtures

We present a preliminary numerical study of spinodal decomposition in an asymmetric polymer mixture, i.e., of polymer with different chain lengths, in three dimensions with full Flory-Huggins-de Gennes free energy, numerically integrating the time evolution equations for the conserved order parameter. For the sake of comparison, we also present a numerical study of the symmetric polymer mixture. The results indicate that the scaled structure factor for the asymmetric polymer mixture is much broader than that of a symmetric polymer mixture. It is interesting that the growth exponents are not symmetric around the critical quench, i.e., growth exponents on the two sides of the critical composition are different. In addition to that, the magnitudes of the pair correlation functions of asymmetric mixtures are very small for x larger than the characteristic domain size r_g and the oscillations seen in the symmetric mixture are almost absent. We have attributed this finding to the rough interfaces and broader domain size distribution in the phase separated asymmetric polymer mixtures. Therefore, the simulation reveals that the asymmetry plays an important role for the spinodal decomposition dynamics of polymer mixtures.     

Size distribution of aggregates in micellar solutions

A general expression for the equilibrium size distribution of polydisperse ionic aggregates in micellar solutions is proposed. This expression accounts for the interactions between the micelles. The interaction is modeled via screened electrostatic potential. Asymptotic formulae for nearly monodisperse case and for high polydispersity of the micelles are presented. The results show that the interactions lead to the stabilization of the monodispersity in the first case. In the second case the interactions cause an increase in the polydispersity of the micellar solution.     

Effects of colloid polydispersity on the phase behavior of colloid-polymer mixtures

We study theoretically the equilibrium phase behavior of a mixture of polydisperse hard-sphere colloids and monodisperse polymers, modeled using the Asakura-Oosawa model [S. Asakura and F. Oosawa, J. Chem. Phys. 22, 1255 (1954)] within the free volume approximation of H. N. W. Lekkerkerker, W. C. K. Poon, P. N. Pusey, A. Stroobants, and P. B. Warren [Europhys. Lett. 20, 559 (1992)]. We compute full phase diagrams in the plane of colloid and polymer volume fractions, using the moment free energy method. The intricate features of phase separation in pure polydisperse colloids combine with the appearance of polymer-induced gas-liquid coexistence to give a rich variety of phase diagram topologies as the polymer-colloid size ratio xi and the colloid polydispersity delta are varied. Quantitatively, we find that polydispersity disfavors fluid-solid against gas-liquid separation, causing a substantial lowering of the threshold value xi(c) above which stable two-phase gas-liquid coexistence appears. Phase splits involving two or more solids can occur already at low colloid concentration, where they may be kinetically accessible. We also analyze the strength of colloidal size fractionation. When a solid phase separates from a fluid, its polydispersity is reduced most strongly if the phase separation takes place at low colloid concentration and high polymer concentration, in agreement with experimental observations. For fractionation in gas-liquid coexistence we likewise find good agreement with experiment, as well as with perturbative theories for near-monodisperse systems.     

Anomalous potentials from inverse analyses of interfacial polydisperse attractive colloidal fluids

This paper investigates effects of using monodisperse inverse analyses to extract particle-particle and particle-surface potentials from simulated interfacial colloidal fluids of polydisperse attractive particles. Effects of polydispersity are investigated as functions of particle concentration and attractive well depth and range for van der Waals and depletion potentials. Forward Monte Carlo simulations are used to generate particle distribution functions for polydisperse interfacial colloidal fluids from which inverted potentials are obtained using an inverse Ornstein-Zernike analysis and an inverse Monte Carlo simulation method. Attractive potentials are successfully recovered for monodisperse colloidal fluids, but polydispersity that is unaccounted for in inverse analyses produces (1) apparent softening of strong forces, (2) anomalous repulsive and attractive interactions, and (3) aphysical particle overlaps. This investigation provides insights into the role of polydispersity in altering the equilibrium structure and corresponding inverted potentials of attractive colloidal fluids near surfaces. These findings should assist the design and interpretation of optical microscopy experiments involving interfacial colloidal fluids similar to the simulated experiments reported here.     

Investigation of structure factors in dense colloidal suspensions using frequency domain photon migration: polydisperse systems

Frequency domain photon migration (FDPM) technique was employed to investigate the structure factors of dense, polydisperse colloidal suspensions. The angle-integrated structure factors, [S(q)], extracted from FDPM measurements of scattering properties at volume fractions ranging from 0.05 to 0.4, were compared with the values predicted from the polydisperse hard sphere Percus-Yevick (HSPY) model, as well as decoupling approximation (DA) and local monodisperse approximation (LMA) models that incorporated independently measured particle size information. Results show that the polydisperse HSPY model is the most suitable for accounting for particle interactions which predominantly arise from volume exclusion effects. Furthermore, the influence of size polydispersity upon [S(q)] is most significant at high volume fractions. The static structure factors at small wave vector q, S(0), were also assessed from dual wavelength FDPM measurements by using the small wave number approximation as well as the local monodisperse approximation. The measured S(0) agrees well with the values predicted by the polydisperse HSPY model.     

Isotropic-nematic phase transition of nonaqueous suspensions of natural clay rods

A novel model system for studying the behavior of hard colloidal rods is presented, consisting of sterically stabilized particles of natural sepiolite clay. Electron microscopy and scattering results confirmed that the organophilic clay particles were individual, rigid rods when dispersed in organic solvents. With a length-to-diameter ratio of approximately 27, the particles showed nematic ordering for volume fractions phi > 0.06. Polarizing microscopy revealed that the phase separation process involved nucleation, growth, and coalescence of nematic domains. The phase volumes and particle concentrations in the coexisting phases were determined. The dependence of these quantities on the total concentration of the suspension agrees well with Onsager's [Ann. N. Y. Acad. Sci. 51, 627 (1949)] isotropic-nematic phase transition theory extended to bidisperse and polydisperse rod systems, and with previous experimental results for rigid rodlike particles. Particle size distributions were obtained by analyzing transmission electron microscopy images. A significant fractionation with respect to rod length (but not diameter) was observed in the coexisting isotropic and nematic phases. The relative polydispersity of both daughter phases was distinctly smaller than that of the parent suspension. The phase behavior of these daughter fractions agrees well with the predictions for hard spherocylinders of corresponding aspect ratios. An isotropic-nematic-nematic phase equilibrium was seen to develop in phase separated samples after 1 month standing and is ascribed to the effect of polydispersity and possibly gravity. The second nematic phase appearing is dominated by very long rods.     

Phase behavior and interfacial properties of nonadditive mixtures of Onsager rods

Within a second virial theory, we study bulk phase diagrams as well as the free planar isotropic-nematic interface of binary mixtures of nonadditive thin and thick hard rods. For species of the same type, the excluded volume is determined only by the dimensions of the particles, whereas for dissimilar ones it is taken to be larger or smaller than that, giving rise to a nonadditivity that can be positive or negative. We argue that such a nonadditivity can result from modeling of soft interactions as effective hard-core interactions. The nonadditivity enhances or reduces the fractionation at isotropic-nematic (IN) coexistence and may induce or suppress a demixing of the high-density nematic phase into two nematic phases of different composition (N(1) and N(2)), depending on whether the nonadditivity is positive or negative. The interfacial tension between coexisting isotropic and nematic phases shows an increase with increasing fractionation at the IN interface, and complete wetting of the IN(2) interface by the N(1) phase upon approach of the triple-point coexistence. In all explored cases bulk and interfacial properties of the nonadditive mixtures exhibit a striking and quite unexpected similarity with the properties of additive mixtures of different diameter ratio.