Zhao-Yan Sun

A novel self-consistent-field lattice model for block copolymers

We develop a self-consistent-field lattice model for block copolymers and propose a novel and general method to solve the self-consistent-field equations. The approach involves describing the polymer chains in a lattice and employing a two-stage relaxation procedure to evolve a system as rapidly as possible to a free-energy minimum. In order to test the validity of this approach, we use the method to study the microphases of rod-coil diblock copolymers. In addition to the lamellar and cylindrical morphologies, micellar, perforated lamellar, gyroid, and zigzag structures have been identified without any prior assumption of the microphase symmetry. Furthermore, this approach can also give the possible orientation of the rods in different structures.

Self-assembly of rod-coil-rod ABA-type triblock copolymers

Self-assembled behavior of symmetric ABA rod-coil-rod triblock copolymer melts is studied by applying self-consistent-field lattice techniques in three-dimensional space. The phase diagram is constructed to understand the effects of the chain architecture on the self-assembled behavior. Four stable structures are observed for the ABA rod-coil-rod triblock, i.e., spherelike, lamellar, gyroidlike, and cylindrical structures. Different from AB rod-coil diblock and BAB coil-rod-coil triblock copolymers, the lamellar structure observed in ABA rod-coil-rod triblock copolymer melts is not stable for high volume fraction of the rod component (f(rod)=0.8), which is attributed to the intramolecular interactions between the two rod blocks of the polymer chain. When 0.3<f(rod)<0.7, the intramolecular interactions between the two rigid blocks of the polymer chain are decrease, which results in the occurrence of some interesting metastable mixed structures. These structures have not been observed in polymers containing only one rod block, such as a lamellar-alt-lamellar structure. The results are expected to provide guidance for the design of microstructures in experiments.

Model, self-assembly structures, and phase diagram of soft Janus particles

Janus particles exhibit interesting self-assembly behavior and functional performances. In particular, soft and deformable Janus particles, as diverse as Janus micelles, Janus microgels, and Janus dendrimers, should receive more attention due to their unique chemical and physical properties and enormous potential applications. Gaining control over precise and predictable self-assembled structures and understanding the fundamental details of self-assembly remain a formidable challenge. Here we present a novel mesoscale model for soft Janus particles, which successfully reflects their physical nature by directly mapping onto experimentally measurable particle properties. By properly tuning Janus balance and the strength of attraction between attractive patches, soft Janus particles can reversibly self-assemble into a number of fascinating hierarchical superstructures in dilute solutions, such as micelles, wormlike strings, single helices, double helices, bilayers, tetragonal bilayers, and complex supermicelles. Our work demonstrates that soft Janus particles with deformable and non-centrosymmetric characteristics hide many surprises in the design and fabrication of hierarchically self-assembled superstructures.

Study of self-assembly of symmetric coil-rod-coil ABA-type triblock copolymers by self-consistent field lattice method

The self-assembly of symmetric coil-rod-coil ABA-type triblock copolymer melts is studied by applying self-consistent field lattice techniques in a three-dimensional space. The self-assembled ordered structures differ significantly with the variation of the volume fraction of the rod component, which include lamellar, wave lamellar, gyroid, perforated lamellar, cylindrical, and spherical-like phases. To understand the physical essence of these phases and the regimes of occurrence, we construct the phase diagram, which matches qualitatively with the existing experimental results. Compared with the coil-rod AB diblock copolymer, our results revealed that the interfacial grafting density of the separating rod and coil segments shows important influence on the self-assembly behaviors of symmetric coil-rod-coil ABA triblock copolymer melts. We found that the order-disorder transition point changes from f(rod)=0.5 for AB diblock copolymers to f(rod)=0.6 for ABA triblock copolymers. Our results also show that the spherical-like and cylindrical phases occupy most of the region in the phase diagram, and the lamellar phase is found stable only at the high volume fraction of the rod.

A simulation model for soft triblock Janus particles and their ordered packing

We present a mesoscale simulation model that is suitable for describing the deformable and non-centrosymmetric characteristics of soft triblock Janus particles. The model parameters are readily mapped onto experimental systems under different ambient conditions. We examine the influence of Janus balance and the flexibility of Janus particle aggregates on the packing structures. Some ordered structures, such as the hexagonal columnar structure and the body-centered tetragonal structure, are observed in our simulations. Our study demonstrates that the Janus balance and the flexibility of Janus particle aggregates can be tuned to obtain various ordered packing structures. The soft Janus particles with soft and deformable characteristics may bring new excitement to materials science.

A molecular-dynamics simulation study on the dependence of Lennard-Jones gas-liquid phase diagram on the long-range part of the interactions

The particle-transfer molecular-dynamics technique is adopted to construct the Lennard-Jones fluid gas-liquid phase diagram. Detailed study of the dependence of the simulation results on the system size and the cutoff distance is performed to test the validity of the simulation technique. Both the traditional cutoff plus long-range correction (CPC) and Ewald summation methods are used in the simulations to calculate the interactions. In the intermediate range of temperatures, the results with the Ewald summation method are almost the same as those with the CPC method. However, in the range close to the critical point, the results with the CPC method deviate from those with the Ewald summation. Compared with the results obtained via the Ewald summation in a smaller system, simply increasing the system size in the CPC scheme may not give better results.

A possible route to fabricate patchy nanoparticles via self-assembly of a multiblock copolymer chain in one step

By using computer simulations, we propose a simple route to fabricate 7–17 nm particles with controllable patch symmetry viaself-assembly of a polymer chain in one step. A single chain of polystyrene-polymethylmethacrylate (PS-PMMA) multiblock copolymer, which is intended to be a generic representative for common hydrophobic multiblock copolymers, is used to fabricate the patchy particles in a solvent that is poor for both components. Various kinds of patchy particles, such as one-patch (with C∞v symmetry), two-patch (with D∞h symmetry), three-patch (with D3hsymmetry), four-patch (with Td symmetry), and cross-ribbon patchy particles, have been obtained. Our work demonstrates that a rational bottom-up design of patchy nanoparticles with controllable symmetry is possible by manipulating the block copolymer chain length and solvent quality.

Heterogeneous crystallization of hard spheres on patterned substrates

We report a numerical investigation of the crystallization of monodisperse hard spheres on different patterned substrates. We find that the duration of a metastable fluid state, which may last for relatively long time in the homogeneous crystallization, can be deeply reduced or almost eliminated when a substrate is used as the seed of crystallization. In the presence of the square patterned substrate, which has the basic character of the body-centered cubic (bcc) crystal structure, a transient bcc crystal phase was observed, suggesting that the bcc structure could be possible to be stabilized by the square patterned substrates. The process of crystallization becomes complicated when the patterned substrate is incommensurate with the bulk crystal. Furthermore, a purely face-centered cubic or hexagonal close packed crystal phase with stacking faults can be generated by using a certain patterned substrate without other factors such as gravity.

Monte Carlo simulation of a single ring among linear chains: Structural and dynamic heterogeneity

We perform lattice Monte Carlo simulation using the bond-fluctuation model to examine the conformation and dynamic properties of a single small flexible ring polymer in the matrix of linear chains as functions of the degree of polymerization of the linear chains. The average conformation properties as gauged by the mean-square radius of gyration and asphericity parameter are insensitive to the chain length for all the chain lengths examined (30, 100, 300, and 1000). However, in the longer chain (300 and 1000) samples, there is an increased spread in the distribution of the value of these quantities, suggesting structural heterogeneity. The center-of-mass diffusion of the ring shows a rapid decrease with increasing chain length followed by a more gradual change for the two longer chain systems. In these longer chain systems, a wide spread in the value of the apparent self-diffusion coefficient is also observed, as well as qualitatively different square displacement trajectories among the different samples, suggesting heterogeneity in the dynamics. A primitive path analysis reveals that in these long chain systems, the ring can exist in topologically distinct states with respect to threading by the linear chains. Threading by the linear chain can dramatically slow down and in some cases stall the diffusive motion of the ring. We argue that the life times for these topological conformers can be longer than the disentanglement time of the linear chain matrix, so that the ring exhibits nonergodic behavior on time scales less or comparable to the life time of these conformers. Our results suggest a picture of the ring diffusion as one where the diffusion path consists of distinctive segments, each corresponding to a different conformer, with slow interconversion between the different conformers.

Relationship between Structural Gel and Mechanical Gel for ABA Triblock Copolymer in Solutions: A Molecular Dynamics Simulation

Polymer gel exists ubiquitously in our daily life, as in food, cosmetics, drugs, and so on. From the structural point of view, the 3D network can be found in a structural gel. In most experimental work, the gel is identified by the sharp increase in modules; that is, the gel should have similar properties as those of a solid, which is named as mechanical gel. However, not all structural gels have strong mechanical responses. Therefore, studying the relationship between structural gel and mechanical gel is very important. In this work, we investigate the structure and mechanical properties of symmetric ABA copolymers with solvophobic end blocks during the sol-gel transition. Three typical systems with weak, middle, and strong solvophobicities are simulated. It is found that the gelation concentration, gel structure, and mechanical response of structural gel are strongly affected by the solvophobicity of ABA block copolymer. We also find that only the gel formed in strong solvophobic systems has a strong mechanical response. Furthermore, the influence of solvophobicity of A-block on the static and dynamic properties of ABA block copolymers in solutions is also studied to give a molecular understanding of physical gelation.