Reaction-induced morphological transitions in a blend of diblock copolymers and reactive monomers: dissipative particle dynamics simulation

The dissipative particle dynamics (DPD) method is applied to the morphological transitions of microphase-separated domains in a mixture of symmetric AB-diblock copolymers and reactive C-monomers, where polymerization and cross-linking reactions take place among C-monomers. The initial structure for the DPD simulation is an equilibrated cylindrical domain structure prepared by the density-biased Monte Carlo method with density profiles obtained from the self-consistent field theory. By introducing a cross-linking reaction among reactive C-monomers, we confirmed that the DPD simulation reproduces the morphological transitions observed in experiments, where the domain morphology changes due to segregation between A-blocks of diblock copolymers and cross-linking networks of C-monomers. When the cross-linking reaction of C-monomers is sufficiently fast compared to the deformation of the domains, the initial cylindrical domains are preserved, while the distance between the domains increases. On the other hand, when the formation of the cross-linking network is slow, the domains can deform and reconnect with each other in the developing cross-linking network. In this case, we observe morphological transitions from the initial domain morphology with a large-curvature interface to another domain morphology with a smaller-curvature interface, such as the transition from the cylindrical phase to the lamellar phase. We calculated the spatial correlations in the microphase-separated domains and found that such correlations are affected by the speed of the formation of the cross-linking network depending on whether the bridging between microphase-separated domains occurs in a nucleation and growth process or in a spinodal decomposition process.

Simulation of the RAFT polymerization in 3D: steric restrictions and incompatibility between species

In this work we developed a RAFT polymerization model taking into account the main reactions of the experimental RAFT process and implemented that model in dissipative particle dynamics (DPD). With...

Simulation of polymerization induced phase separation in model thermosets

Polymerization induced phase separation (PIPS) in a three component thermoset is studied using molecular dynamics simulations of a new coarse-grained thermoset model. The system includes two crosslinker molecules, which differ in their glass transition temperatures (T_g) and chain length and thus have the potential for phase separation. One crosslinker has a high T_g corresponding to a rubbery behavior, and simulations were performed for a short length (4 beads) and a long length (33 beads). The resin and other crosslinker have low T_g. A coarse-grained model is developed with these features and with interaction parameters determined so that for either rubbery crosslinker length, the system is in the liquid state at the cure temperature. For sufficiently slow reaction rates, the long rubbery molecule exhibits PIPS into a bicontinuous array of nanoscale domains, but the short one does not, reproducing recent experimental results. The simulations demonstrate that the reaction rates must be slow enough to allow diffusion to yield phase separation. Particularly, the reaction rate corresponding to the secondary amine must be very slow, else the structure of crosslinked clusters and the substantially increased diffusion time will prevent PIPS.

Phase Diagram for Ideal Diblock-Copolymer Micelles Compared to Polymerization-Induced Self Assembly

In this work we constructed a detailed phase diagram for the solutions of ideal diblock-copolymers and compared such diagram with that obtained during polymerization-induced self-assembly (PISA); a wide range of polymer concentrations as well as chain compositions was studied. As the length of the solvophobic block nB increases (the length of the solvophilic block nA was fixed), the transition from spherical micelles to cylinders and further to vesicles (lamellae) occurs. We observed a rather wide transition region between the spherical and cylindrical morphology in which the system contains a mixture of spheres and short cylinders, which appear to be in dynamic equilibrium; the transition between the cylinders and vesicles was found to be rather sharp. Next, upon increasing the polymer concentration in the system, the transition region between the spheres and cylinders shifts towards lower nB/nA values; a similar shift but with less magnitude was observed for the transition between the cylinders and vesicles. Such behavior was attributed to the increased number of contacts between the micelles at higher polymer volume concentrations. We also found that the width of the stability region of the cylindrical micelles for small polymer volume concentrations is in good quantitative agreement with the predictions of analytical theory. The obtained phase diagram for PISA was similar to the case of presynthesized diblock copolymer; however, the positions of the transition lines for PISA are slightly shifted towards higher nB/nA values in comparison to the presynthesized diblock copolymers, which is more pronounced for the case of the cylinders-to-vesicles transition. We believe that the reason for such behavior is the polydispersity of the core-forming blocks: The presence of the short and long blocks being located at the micelle interface and in its center, respectively, helps to reduce the entropy losses due to the insoluble block stretching, which leads to the increased stability of more curved micelles.

Copolymerization on Selective Substrates: Experimental Test and Computer Simulations

We explore the influence of a selective substrate on the composition and sequence statistics during the free radical copolymerization. In particular, we study the radical copolymerization of styrene and acrylic acid in bulk and in silica pores of different sizes. We show that the substrate affects both polymer composition and sequence statistics. We use dissipative particle dynamics simulations to study the polymerization process in detail, trying to pinpoint the parameters responsible for the observed differences in the polymer chain composition and sequences. The magnitude of the observed effect depends on the fraction of adsorbed monomer units, which cannot be described in the framework of the copolymerization theories based on the terminal unit model.

Effect of Cross-Linking on the Structure and Growth of Polymer Films Prepared by Interfacial Polymerization

Interfacial polymerization of tri- and bifunctional monomers (A3B2 polymerization) is investigated by dissipative particle dynamics to reveal an effect of cross-linking on the reaction kinetics and structure of the growing polymer film. Regardless of the comonomer reactivity and miscibility, the kinetics in an initially bilayer melt passes from the reaction to diffusion control. Within the crossover period, branched macromolecules undergo gelation, which drastically changes the scenario of the polymerization process. Comparison with the previously studied linear interfacial polymerization (Berezkin, A. V.; Kudryavtsev, Y. V. Linear Interfacial Polymerization: Theory and Simulations with Dissipative Particle Dynamics J. Chem. Phys. 2014, 141, 194906) shows similar conversion rates but very different product characteristics. Cross-linked polymer films are markedly heterogeneous in density, their average polymerization degree grows with the comonomer miscibility, and end groups are mostly trapped deeply in the film core. Products of linear interfacial polymerization demonstrate opposite trends as they are spontaneously homogenized by a convective flow of macromolecules expelled from the reactive zone to the film periphery, which we call the reactive extrusion effect and which is hampered in branched polymerization. Influence of the comonomer architecture on the polymer film characteristics could be used in various practical applications of interfacial polymerization, such as fabrication of membranes, micro- and nanocapsules and 3D printing.

Mesoscale modelling of environmentally responsive hydrogels: emerging applications

We review recent advances in mesoscale computational modeling, focusing on dissipative particle dynamics, used to probe stimuli-sensitive behavior of hydrogels.

Phase diagrams of block copolymer melts by dissipative particle dynamics simulations

Phase diagrams for monodisperse and polydisperse diblock copolymer melts and a random multiblock copolymer melt are constructed using dissipative particle dynamics simulations. A thorough visual analysis and calculation of the static structure factor in several hundreds of points at each of the diagrams prove the ability of mesoscopic molecular dynamics to predict the phase behavior of polymer systems as effectively as the self-consistent field-theory and Monte Carlo simulations do. It is demonstrated that the order-disorder transition (ODT) curve for monodisperse diblocks can be precisely located by a spike in the dependence of the mean square pressure fluctuation on χN, where χ is the Flory-Huggins parameter and N is the chain length. For two other copolymer types, the continuous ODTs are observed. Large polydispersity of both blocks obeying the Flory distribution in length does not shift the ODT curve but considerably narrows the domains of the cylindrical and lamellar phases partially replacing them with the wormlike micelle and perforated lamellar phases, respectively. Instead of the pure 3d-bicontinuous phase in monodisperse diblocks, which could be identified as the gyroid, a coexistence of the 3d phase and cylindrical micelles is detected in polydisperse diblocks. The lamellar domain spacing D in monodisperse diblocks follows the strong-segregation theory prediction, D∕N(1∕2) ~ (χN)(1∕6), whereas in polydisperse diblocks it is almost independent of χN at χN < 100. Completely random multiblock copolymers cannot form ordered microstructures other than lamellas at any composition.