Parallel Algorithm for Cell Dynamics Simulation of Block Copolymers

Hybrid Time-Dependent Ginzburg–Landau Simulations of Block Copolymer Nanocomposites: Nanoparticle Anisotropy

Block copolymer melts are perfect candidates to template the position of colloidal nanoparticles in the nanoscale, on top of their well-known suitability for lithography applications. This is due to their ability to self-assemble into periodic ordered structures, in which nanoparticles can segregate depending on the polymer–particle interactions, size and shape. The resulting coassembled structure can be highly ordered as a combination of both the polymeric and colloidal properties. The time-dependent Ginzburg–Landau model for the block copolymer was combined with Brownian dynamics for nanoparticles, resulting in an efficient mesoscopic model to study the complex behaviour of block copolymer nanocomposites. This review covers recent developments of the time-dependent Ginzburg–Landau/Brownian dynamics scheme. This includes efforts to parallelise the numerical scheme and applications of the model. The validity of the model is studied by comparing simulation and experimental results for isotropic nanoparticles. Extensions to simulate nonspherical and inhomogeneous nanoparticles are discussed and simulation results are discussed. The time-dependent Ginzburg–Landau/Brownian dynamics scheme is shown to be a flexible method which can account for the relatively large system sizes required to study block copolymer nanocomposite systems, while being easily extensible to simulate nonspherical nanoparticles.

Large scale three dimensional simulations of hybrid block copolymer/nanoparticle systems

Block copolymer melts self-assemble in the bulk into a variety of nanostructures, making them perfect candidates to template the position of nanoparticles. The morphological changes of block copolymers are studied in the presence of a considerable filling fraction of colloids. Furthermore, colloids can be found to assemble into ordered hexagonally close-packed structures in a defined number of layers when softly confined within the phase-separated block copolymer. A high concentration of interface-compatible nanoparticles leads to complex long-lived block copolymer morphologies depending on the polymeric composition. Macrophase separation between the colloids and the block copolymer can be induced if colloids are unsolvable within the matrix. This leads to the formation of ellipsoid-shaped polymer-rich domains elongated along the direction perpendicular to the interface between block copolymer domains.

Nucleation of the BCC phase from disorder in a diblock copolymer melt: Testing approximate theories through simulation

We examine nucleation of the stable body-centred-cubic (BCC) phase from the metastable uniform disordered phase in an asymmetric diblock copolymer melt. Our comprehensive, large-scale simulations of the time-dependent, mean-field Landau-Brazovskii model find that spherical droplets of the BCC phase nucleate directly from disorder. Near the order-disorder transition, the critical nucleus is large and has a classical profile, attaining the bulk BCC phase in an interior that is separated from disorder by a sharp interface. At greater undercooling, the amplitude of BCC order in the interior decreases and the nucleus interface broadens, leading to a diffuse critical nucleus. This diffuse nucleus becomes large as the simulation approaches the disordered phase spinodal. We show that our simulation follows the same nucleation pathway that Cahn and Hilliard found for an incompressible two-component fluid, across the entire metastable region. In contrast, a classical nucleation theory calculation based on the free energy of a planar interface between coexisting BCC and disordered phases agrees with simulation only in the limit of very small undercooling; we can expand this region of validity somewhat by accounting for the curvature of the droplet interface. A nucleation pathway involving a classical droplet persists, however, to deep undercooling in our simulation, but this pathway is energetically unfavourable. As a droplet grows in the simulation, its interface moves with a constant speed, and this speed is approximately proportional to the undercooling.

Block copolymer nanocontainers

Using cell dynamics computer simulation, we perform a systematic study of thin block copolymer films around a nanoparticle. Lamellar-, cylinder-, and sphere-forming block copolymers are investigated with respect to different film thicknesses, particle radii, and boundary conditions at the film interfaces. The obtained structures include standing lamellae and cylinders, "onions", cylinder "knitting balls", "golf ball", layered spherical, "virus"-like and mixed morphologies with T-junctions and U-type defects. The kinetics of the structure formation and difference with planar thin films are discussed. Our simulations suggest that novel porous nanocontainers can be formed by the coating of a sacrificial nanobead by a block copolymer layer with a well-controlled nanostructure. In addition, first scanning force microscopy experiments on a model system reveal surface structures similar to those predicted by our simulations.

Kinetic pathways of gyroid-to-cylinder transitions in diblock copolymers under external fields: cell dynamics simulation

Using cell dynamics simulation we investigate the cubic gyroid morphology of block copolymer melts under simple shear flow and electric field. The electric field should be stronger than a certain critical value to induce transition to a cylindrical phase. In the case of simple steady shear the gyroid-to-cylinder transition was observed even for a very weak shear. Quantitative analysis of pathways of gyroid-to-cylinder transition is performed by means of Minkowski functionals. We found that the kinetics of the gyroid-to-cylinder transition are different under electric field and shear flow. Moreover, the gyroid structure under different strengths of electric field shows different pathways. Different types of intermediates such as five-, four-fold connections and "winding" cylinders are found for different pathways.