Scaling properties of the morphological measures at the early and intermediate stages of the spinodal decomposition in homopolymer blends

Mass-transfer driven spinodal decomposition in a ternary polymer solution

Nonsolvent induced phase separation (NIPS) is a widely occuring process used in industrial membrane production, nanotechnology and Nature to produce microstructured polymer materials. A variety of process-dependent morphologies are produced when a polymer solution is exposed to a nonsolvent that, following a period where mass is exchanged, precipitates and solidifies the polymer. Despite years of investigation, both experimental and theoretical, many questions surround the pathways to the microstructures that NIPS can produce. Here, we provide simulation results from a model that simultaneously captures both the processess of solvent/nonsolvent exchange and phase separation. We show that the time it takes the nonsolvent to diffuse to the bottom of the film is an important timescale, and that phase separation is possible at times both much smaller and much larger than this scale. Our results include both one-dimensional simulations of the mass transfer kinetics and two- and three-dimensional simulations of morphologies at both short and long times. We find good qualitative agreement with experimental heuristics, but we conclude that an additional model for the vitrification process will be key for fully explaining experimental observations of microstructure formation.

A multi-fluid model for microstructure formation in polymer membranes

We develop a multi-fluid model for a ternary polymer solution using the Rayleighian formalism of Doi and Onuki, and give an efficient pseudo-spectral method for solving both the diffusion and momentum equations that result. Subsequently, we find that the numerical simulation is capable of describing systems at the micron length-scale and easily reaches millisecond time-scales. In addition, we characterize the model thermodynamics and kinetics including the (i) phase behavior, (ii) structure of the interfaces, (iii) mutual diffusion coefficients, (iv) bulk spinodal decomposition kinetics with and without hydrodynamics and (v) spinodal decomposition in the presence of an interface with a non-solvent bath. We obtain good qualitative agreement with the expected thermodynamic and kinetic behavior. We also show that a linear stability analysis of the diffusion equation quantitatively predicts the fastest growing mode obtained from simulation and gives insight into the phase separation process relevant for the evolution of microstructure in phase-separating ternary polymer solutions.

Phase separation in thermal systems: a lattice Boltzmann study and morphological characterization

We investigate thermal and isothermal symmetric liquid-vapor separations via a fast Fourier transform thermal lattice Boltzmann (FFT-TLB) model. Structure factor, domain size, and Minkowski functionals are employed to characterize the density and velocity fields, as well as to understand the configurations and the kinetic processes. Compared with the isothermal phase separation, the freedom in temperature prolongs the spinodal decomposition (SD) stage and induces different rheological and morphological behaviors in the thermal system. After the transient procedure, both the thermal and isothermal separations show power-law scalings in domain growth, while the exponent for thermal system is lower than that for isothermal system. With respect to the density field, the isothermal system presents more likely bicontinuous configurations with narrower interfaces, while the thermal system presents more likely configurations with scattered bubbles. Heat creation, conduction, and lower interfacial stresses are the main reasons for the differences in thermal system. Different from the isothermal case, the release of latent heat causes the changing of local temperature, which results in new local mechanical balance. When the Prandtl number becomes smaller, the system approaches thermodynamical equilibrium much more quickly. The increasing of mean temperature makes the interfacial stress lower in the following way: σ=σ(0)[(T(c)-T)/(T(c)-T(0))](3/2), where T(c) is the critical temperature and σ(0) is the interfacial stress at a reference temperature T(0), which is the main reason for the prolonged SD stage and the lower growth exponent in the thermal case. Besides thermodynamics, we probe how the local viscosities influence the morphology of the phase separating system. We find that, for both the isothermal and thermal cases, the growth exponents and local flow velocities are inversely proportional to the corresponding viscosities. Compared with the isothermal case, the local flow velocity depends not only on viscosity but also on temperature.

Coarse-grained model for mechanosensitive ion channels

In this work, a coarse-grained model for a kind of membrane protein, the mechanosensitive channel of small conductance (MscS), is proposed. The basic structure of the MscS is preserved when the protein is coarse-grained. For the coarse-grained model, the channels show two different states, namely the open and closed states, depending on the model parameters. Under the same membrane tension, the state of the ion channel is found to be critically determined by the protein structure, especially the length of the transmembrane alpha-helix. It is also found that for the protein with certain size, the gating transition occurs when the membrane tension is applied, resembling in a real mechanosensitive channel.

Characterization of the dynamics of block copolymer microdomains with local morphological measures

We investigate the structure formation in thin films of cylinder forming block copolymers. With in situ scanning probe microscopy image sequences can be recorded with high temporal (2 min per frame) and spatial (10 nm) resolution. We compare different image processing methods for quantitative analysis of the large amount of data. Computing local Minkowski functionals yields local geometrical and morphological information about the observed structures and enables us to track their evolution with time. An alternative characterization method is to reduce the gray scale images to their skeleton and to classify and count the branching points of the skeletonized structure. We tracked the temporal evolution of these measures and computed correlation functions.

Coarsening of bicontinuous structures via nonconserved and conserved dynamics

Coarsening subsequent to phase separations occurs in many two-phase mixtures. While unique scaled particle size distributions have been determined for highly asymmetric mixtures in which spherical particles form in a matrix, it is not known if a unique scaled structure exists for symmetric mixtures, which yield bicontinuous structures having intricately interpenetrating phase domains. Using large-scale simulations, we have established that unique scaled microstructures exist in symmetric mixtures evolving via nonconserved and conserved dynamics. We characterize their morphologies by the interfacial shape distribution, a counterpart to the particle size distribution, and their topologies by the genus. We find that the two dynamics result in unique, but different, scaled interfacial shape distributions, with conserved dynamics yielding a narrower distribution around zero mean curvature. In contrast, the two scaled structures are topologically similar, having nearly equal values of the scaled genus.

Viscoelastic phase separation of protein solutions

In addition to the known behavior of normal phase separation and gelation, we report novel phase-separation behavior of protein solutions as their intermediate case. A network structure of the protein-rich phase may be formed even if it is the minority phase, contrary to the conventional wisdom. This behavior is characteristic of viscoelastic phase separation found in polymer solutions. This kinetic pathway may play crucial roles in the complex phase ordering of protein solutions, in particular, protein network formation in biological systems and foods.

Kinetic pathways of sheared block copolymer systems derived from Minkowski functionals

We employ Minkowski functionals to analyze the kinetics of pattern formation under an applied external shear flow. The considered pattern formation model describes the dynamics of phase separating block copolymer systems. For our purpose, we have chosen two block copolymer systems (a melt and a solution) that exhibit a hexagonal cylindrical morphology as an equilibrium structure. Our main objective is the determination of efficient choices for the treshold values that are required for the calculation of the Minkowski functionals. We find that a minimal set of two treshold values (one from which should be equal to an average density value and another to a higher density value) is sufficient to unraffle the phase separation kinetics. Given these choices, we focus on the influence of the degree of phase separation, and the instance at which the shear is applied, on the kinetic pathways. We also found a remarkable similarity of the time evolution of Euler characteristic and the segregation parameter for the average density choice.

The unphysical pinning of the domain growth during the separation of homopolymer blends near the spinodal

We simulate model B of mesoscopic dynamic with the Flory-Huggins free energy for the homopolymer blend. We concentrate the study on the rescaling of the spatial coordinates in the model. We show that the commonly used rescaling of the spatial coordinates, by the function vanishing at the spinodal, leads to the unphysical freezing of the domains. We study the evolution of the system in two different processes: One is the growth process induced by the temperature quench to the metastable or unstable region of the phase diagram and the second is the quench-jump process in which we first allow the system to separate below the spinodal and next heat it up below or above the spinodal. The proper rescaling avoiding the unphysical pinning of the domain growth at the spinodal is proposed.

Morphological changes during the order-disorder transition in the two- and three-dimensional systems of scalar nonconserved order parameters

The order-disorder transition is studied in a system of a scalar nonconserved order parameter. We use this well studied system to show that the application of the methods of topology and geometry reveals that our knowledge of the kinetic pathways by which the order-disorder transition proceeds is far from being complete. We show that in two-dimensional (2D) and 3D systems there are three dynamical regimes in the evolution of the system: early, intermediate, and late. In the intermediate regime two length scales govern the behavior of the system, whereas in the early and intermediate regime there is only one length scale. The size distribution of the domain area indicates the pathway by which the domains change their size. There are only two types of domains in a 2D system: circular and elongated with well defined characteristics (scaling of the area with the contour length) which in the late regime do not depend on time after rescaling by the average area and contour in the system. The elongated domains continuously change into circular domains reducing in this way the overall dissipation in the system. In order to reach a Lifshitz-Cahn-Allen (LCA) late stage regime the number of elongated domains must be strongly reduced. In the intermediate regime the number of elongated domains is large and simple LCA scaling does not hold. In a 3D symmetric system we always have a bicontinuous structure that evolves by cutting small connections. The late stage regime seems to be associated with the appearance of the preferred nonzero mean curvature. The early-intermediate regime crossover is associated with the saturation of the order parameter inside the domains, while the intermediate-late stage regime crossover is related to the global breaking of the +/- order parameter symmetry (marked by the appearance of the nonzero mean curvature but still zero average magnetization). The times for the occurrence of these crossovers do not depend on the size of the system.

Dynamics of spinodal decomposition in finite-lifetime systems: Nonlinear statistical theory based on a coarse-grained lattice-gas model

We study theoretically dynamics of the spinodal decomposition in finite-lifetime systems to clarify effects of the interparticle interactions beyond the Ginzburg-Landau-Wilson phenomenology. Our theory is based on the coarse-grained Hamiltonian derived from the interacting lattice-gas model with a finite lifetime. The information of a system is reduced to closed-form coupled integrodifferential equations for the single-point distribution function and the dynamical structure factor. These equations involve explicitly the interparticle interactions. The finite lifetime prevents the phase separation and the order formation in the cw creation case; domains cannot grow to be larger than an asymptotic characteristic size [k(max)(t --> infinity)](-1). Power-law dependence of k(max)(t --> infinity) on the interparticle interaction and the particle lifetime is also found numerically. The finite lifetime prevents the phase separation, i.e., the lower critical wave number k((1))(c) appears and domains of size larger than [k((1))(c)](-1) cannot grow.

Scaling of the Euler characteristic, surface area, and curvatures in the phase separating or ordering systems

We present robust scaling laws for the Euler characteristic and curvatures applicable to any symmetric system undergoing phase separating or ordering kinetics. We apply it to the phase ordering in a system of the nonconserved scalar order parameter and find three scaling regimes. The appearance of the preferred nonzero curvature of an interface separating +/- domains marks the crossover to the late stage regime characterized by the Lifshitz-Cahn-Allen scaling.

Influence of the free-energy functional form on simulated morphology of spinodally decomposing blends

The spinodal decomposition of a binary mixture has been studied within several mesoscopic models. It has been found that the form of the equilibrium free energy has a crucial effect on the morphological development in asymmetric blends. We have shown that the principal quantity that determines the topology of the interface (and type of morphology) is the equilibrium minority phase volume fraction, while the transition from bicontinuous to droplet morphology can be treated as a percolation. The concentration dependence of the square gradient coefficient attributed for the Flory-Huggins-de Gennes free energy has no significant influences on the average domain growth, but can be distinguished experimentally from its constant-coefficient alternative by measuring the maximum wave vector of the scattering intensity as a function of the minority phase volume fraction for spinodally decomposing asymmetric blends. The concentration dependence of the Onsager coefficient has the weak, systematic effect of slowing down the morphological development. The local shape of the interface is not affected considerably by the concentration dependence of the square gradient and Onsager coefficient.