Monte Carlo study of the microphase separation of cross-linked polymer blends

Phase diagram of selectively cross-linked block copolymers shows chemically microstructured gel

We study analytically the intricate phase behavior of cross-linked AB diblock copolymer melts, which can undergo two main phase transitions due to quenched random constraints. Gelation, i.e., spatially random localisation of polymers forming a system-spanning cluster, is driven by increasing the number parameter μ of irreversible, type-selective cross-links between random pairs of A blocks. Self-assembly into a periodic pattern of A/B-rich microdomains (microphase separation) is controlled by the AB incompatibility χ inversely proportional to temperature. Our model aims to capture the system's essential microscopic features, including an ensemble of random networks that reflects spatial correlations at the instant of cross-linking. We identify suitable order parameters and derive a free-energy functional in the spirit of Landau theory that allows us to trace a phase diagram in the plane of μ and χ. Selective cross-links promote microphase separation at higher critical temperatures than in uncross-linked diblock copolymer melts. Microphase separation in the liquid state facilitates gelation, giving rise to a novel gel state whose chemical composition density mirrors the periodic AB pattern.

Self-consistent field approach for cross-linked copolymer materials

A generalized self-consistent field approach for polymer networks with a fixed topology is developed. It is shown that the theory reproduces the localization of cross-links, which is characteristic for gels. The theory is then used to study the order-disorder transition in regular networks of end-linked diblock copolymers. Compared to diblock copolymer melts, the transition is shifted towards lower values of the incompatibility parameter χ (the Flory-Huggins parameter). Moreover, the transition becomes strongly first order already at the mean-field level. If stress is applied, the transition is further shifted and finally vanishes in a critical point.

Self-organized stiffness in regular fractal polymer structures

We investigated membrane-like polymer structures of fractal connectivity such as Sierpinski gaskets and Sierpinski carpets applying the bond fluctuation model in three dimensions. Without excluded volume (phantom), both polymeric fractals obey Gaussian elasticity on larger scales determined by their spectral dimension. On the other hand, the swelling effect due to excluded volume is rather distinct between the two polymeric fractals: Self-avoiding Sierpinski gaskets can be described using a Flory-type mean-field argument. Sierpinski carpets having a spectral dimension closer to perfect membranes are significantly more strongly swollen than predicted. Based on our simulation results it cannot be excluded that Sierpinski carpets in athermal solvent show a flat phase on larger scales. We tested the self-consistency of Flory predictions using a virial expansion to higher orders. From this we conclude that the third virial coefficient contributes marginally to Sierpinski gaskets, but higher order virial coefficients are relevant for Sierpinski carpets.

Polymer-decorated tethered membranes under good- and poor-solvent conditions

We study tethered membranes grafted by polymer chains on one side. Mean-field and scaling arguments predicting a spontaneous curvature are compared to the results of lattice-based Monte Carlo simulations using the Bond Fluctuation Model, which are carried out for various grafting densities and chain lengths. We show that already slightly overlapping chains bend the membrane significantly. This proves the entropic origin for the bending stiffness, which is of order kT. To understand the membrane curvature under conditions of very small bending stiffness we apply a geometrical model which takes into account the state of chains at the overlap threshold. Applying a thermal solvent model for the grafted chains, we demonstrate that the bending direction of the membrane can be triggered by variation of the solvent quality. This indicates that polymer-decorated membranes may serve as switchable nanoscale devices.

Effect of network topology on phase separation in two-dimensional Lennard-Jones networks

We generate two-dimensional Lennard-Jones networks with random topology by preparing a perfect four-functional network of identical harmonic springs and randomly cutting some of the springs. Using molecular-dynamics simulations we find that the fraction p of active springs affects both the temperature of phase separation and the type of structures observed below this temperature, from networklike high-density patterns at p>0.5 ("gel") to dropletlike structures at p<0.5 ("sol"). In the gel domain, these patterns are determined by the interplay between free energy and network topology, with the former dominant as p-->1 and the latter as p-->0.5 .

Microphase separation in cross-linked polymer blends. Efficient replica RPA post-processing of simulation data for homopolymer networks

We investigate the behaviour of randomly cross-linked (co)polymer blends using a combination of replica theory and large-scale molecular dynamics simulations. In particular, we derive the analogue of the random phase approximation for systems with quenched disorder and show how the required correlation functions can be calculated efficiently. By post-processing simulation data for homopolymer networks we are able to describe neutron scattering measurements in heterogeneous systems without resorting to microscopic detail and otherwise unphysical assumptions. We obtain structure function data which illustrate the expected microphase separation and contain system-specific information relating to the intrinsic length scales of our networks.

Glassy correlations and microstructures in randomly cross-linked homopolymer blends

We consider a microscopic model of a polymer blend that is prone to phase separation. Permanent cross-links are introduced between randomly chosen pairs of monomers, drawn from the Deam-Edwards distribution. Thereby, not only density but also concentration fluctuations of the melt are quenched-in in the gel state, which emerge upon sufficient cross-linking. We derive a Landau expansion in terms of the order parameters for gelation and phase separation, and analyze it on the mean-field level, including Gaussian fluctuations. The mixed gel is characterized by thermal as well as time-persistent (glassy) concentration fluctuations. Whereas the former are independent of the preparation state, the latter reflect the concentration fluctuations at the instant of cross-linking, provided the mesh size is smaller than the correlation length of phase separation. The mixed gel becomes unstable to microphase separation upon lowering the temperature in the gel phase. Whereas the length scale of microphase separation is given by the mesh size, at least close to the transition, the emergent microstructure depends on the composition and compressibility of the melt. Hexagonal structures, as well as lamellas or random structures with a unique wavelength, can be energetically favorable.

Computer simulation of polymer networks: Swelling by binary Lennard-Jones mixtures

The swelling of regular, tightly meshed model networks is investigated by a molecular-dynamics-Monte Carlo hybrid technique. The chemical equilibrium between two simulation boxes representing the gel phase and a solvent bath, respectively, is obtained by subjecting the Lennard-Jones particles of a binary mixture, serving as explicit solvent, to the particle transfer step of Gibbs ensemble-Monte Carlo. The swelling behavior, especially preferential absorption of a single component, whose dependence on temperature, pressure, and fluid composition is studied, also depends significantly on the size of the central simulation cell. These finite-size effects correlate well with those exhibited by the density of solvent-free (dry) networks. A theoretical expression, whose derivation is based on network elasticity (of dry networks) yields finite-size scaling behavior in good accord with simulation results for both dry networks and gels in contact with solvent baths. This expression can be used to extrapolate the swelling behavior of simulated finite systems to infinite system size.

Finite size effects in tightly meshed polymer networks

Molecular dynamics computer simulations on regular, tightly meshed model networks exhibit variations of the network density with system size. We show that these variations are due mainly to network elasticity. A theoretical expression derived on the basis of the self-consistent-field approach yields finite size scaling behavior in good accord with the simulation for a wide range of thermodynamic conditions.

Swelling of model polymer networks with different cross-link densities: a computer simulation study

The swelling of model polymer networks with different cross-link densities is studied via molecular dynamics simulation. During the simulation, the solvent particles, consisting of one interaction center or six interaction centers, respectively, are transferred between two coupled simulation boxes. The gel box includes both network and solvent particles, whereas the solvent box contains solvent only. The particle transfer is controlled by the solvent chemical potential difference in the two boxes, which is calculated via the Widom test particle method for the one-site solvent and via Rosenbluth sampling for the chainlike solvent. The equilibrium swelling ratio of the network as well as the solvent diffusion coefficient under subcritical and supercritical conditions are computed as functions of the network cross-link density for a wide range of temperatures and pressures. In addition, the simulated swelling behavior is compared to a Flory-Huggins-type theory, which yields qualitative agreement for the systems studied here.