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/N1/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.