Search

Designing complex cluster crystals with a specific function using simple colloidal building blocks remains a challenge in materials science. Herein, we propose a conceptually new design strategy for constructing complex cluster crystals via hierarchical self-assembly of simple soft Janus colloids. A novel and previously unreported colloidal cluster-χ (χc) phase, which resembles the essential structural features of α-manganese but at a larger length scale, is obtained through molecular dynamics simulations. The formation of the χc phase undergoes a remarkable two-step self-assembly process, that is, the self-assembly of clusters with specific size dispersity from Janus colloids, followed by the highly ordered organization of these clusters. More importantly, the dynamic exchange of particles between these clusters plays a critical role in stabilizing the χc phase. Such a conceptual design framework based on intercluster exchange has the potential to effectively construct novel complex cluster crystals by hierarchical self-assembly of colloidal building blocks.

The capability of polyelectrolyte brushes to spontaneously clean oil fouling via water is determined by factors including water wettability and the self-assembled structures of hydrated polyelectrolytes. Although the charged groups of polyelectrolytes provide the original source of water wettability, the self-assembled structures play a significant role in the self-cleaning performances. Here, we employ coarse-grained molecular dynamics simulations to study the general self-cleaning characteristics of two types of surface-grafted polyelectrolyte brushes (i.e., zwitterionic and anionic polyelectrolytes). It has been found that the high grafting density is favorable to fouling reduction for both polyzwitterions and polyanions. To be specific, the hydrated polyzwitterions form an intermolecular cross-linked network via zwitterionic complexes, resulting in better self-cleaning capabilities than the polyanions at lower grafting densities. However, polyanions form bundles with each consisting of several chains via hydrophobic interactions and electrostatic repulsions presenting better self-cleaning performances than the polyzwitterions at higher grafting densities.

Topological defects in liquid crystals under confined geometries have attracted extensive research interests. Here, we perform molecular dynamics simulations to investigate the formation and transition of defect patterns in two-dimensional smectic Gay-Berne liquid crystals with a simple rectangular confinement boundary. Two typical types of defect patterns, bridge and diagonal defect patterns, are observed, which can be transformable continuously between each other over time. The transition usually starts from the line or point defect regions, and the competition between neighboring and opposite boundary effects induces the continuous realignments of the smectic layers to connect the neighboring or opposite walls. The relative stability of these two defect patterns can be controlled by changing the confinement conditions. These results deepen our understanding of transition kinetics of defect patterns in confined liquid crystals.

Reversible chemistries have been extensively explored to construct highly crystalline covalent organic frameworks (COFs) via defect correction. However, the mechanisms of defect correction that can explain the formation of products as single crystals, polycrystal/crystallites, or amorphous solids remain unknown. Herein, we employed molecular dynamics simulations combined with a polymerization model to investigate the growth kinetics of two-dimensional COFs. By virtue of the Arrhenius two-state model describing reversible reactions, we figured out the conditions in terms of active energy and binding energy for different products. Specifically, the ultraslow growth of COFs under high reversibility of reactions corresponding to low binding energies resulted in a single crystal by inhibiting the emergence of nuclei as well as correcting defects through continually dropping small defective fragments off at crystal boundaries. High bonding energies responsible for the high nucleation rate and rapid growth that incorporated defects in crystals and caused the division of crystals through defect correcting processes led to small crystallites or polycrystals. The insights into the mechanisms help us to understand and further control the growth kinetics by exploiting reversible conditions to synthesize COFs of higher quality.

Synthesis of covalent organic frameworks with long-range molecular ordering is an outstanding challenge due to the fact that defects against predesigned topological symmetries are prone to form and break crystallization. The physical origins and controlling parameters of topological defects remain scarcely understood. By virtue of molecular dynamics simulations, we found that pentagons for combination [C4 + C4] and [C4 + C2] and heptagons for [C3 + C3] and [C3 + C2] were initial defects for growth dynamics with both uncontrolled and suppressed nucleation, further inducing more complex defects. The defects can be significantly reduced by achieving the growth with monomers added to a single nucleus, agreeing well with previous simulations and experiments. To understand the nature of defects, we proposed a parameter φ to describe the range of biased rotational angle between two monomers, within which chemical reactions are allowed. The parameter φ shows a monotonic relationship with defect population, which is demonstrated to be highly computable by using density functional theory calculations. When φ < 20, we can even observe defect-free growth for the four combinations, irrespective of growth dynamics. The results are essential for screening and designing condensation reactions for the synthesis of single crystals of high quality.

We describe the algorithm of employing multi-GPU power on the basis of Message Passing Interface (MPI) domain decomposition in a molecular dynamics code, GALAMOST, which is designed for the coarse-...

Hollow polyhedral cages hold great potential for application in nanotechnological and biomedical fields. Understanding the formation mechanism of these self-assembled structures could provide guidance for the rational design of the desired polyhedral cages. Here, by constructing kinetic network models from extensive coarse-grained molecular dynamics simulations, we elucidated the formation mechanism of the dodecahedral cage, which is formed by the self-assembly of patchy particles. We found that the dodecahedral cage is formed through increasing the aggregate size followed by structure rearrangement. Based on this mechanistic understanding, we improved the productivity of the dodecahedral cage through the rational design of the patch arrangement of patchy particles, which promotes the structural rearrangement process. Our results demonstrate that it should be a feasible strategy to achieve the rational design of the desired nanostructures via the kinetic analysis. We anticipate that this methodology could be extended to other self-assembly systems for the fabrication of functional nanomaterials.

Molecular dynamics simulations are useful tools to unveil molecular mechanisms of polymer phase separation, self-assembly, adsorption, and so on. Due to large molecular size and slow relaxation of the polymer chains, a great amount of issues related to large-distance chain displacement cannot be tackled easily with conventional molecular dynamic simulations. Systematic coarse-graining and enhanced sampling methods are two types of improvements that can boost spatiotemporal scales in polymer simulations. We present two typical ways to obtain the coarse-graining potential either by fitting to correct liquid structures or by fitting to available thermodynamic properties of polymer systems. The newly proposed anisotropic coarse-grained particle model can be used to describe aggregation and assembly of polymeric building blocks from disk-like micelles to Janus particles. We also present a stochastic polymerization model combined with coarse-grained simulations to investigate the problems strongly influenced by the coupling of polymerization and excluded volume effects. Finally, a facile implementation of integrated tempering sampling method is illustrated to be very efficient on bypassing local energy minima and having access to true equilibrium polymer structures.

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.

Molecular dynamics simulations are adopted to calculate the equation of state characteristic parameters P*, rho*, and T* of isotactic polypropylene (iPP) and poly(ethylene-co-octene) (PEOC), which can be further used in the Sanchez-Lacombe lattice fluid theory (SLLFT) to describe the respective physical properties. The calculated T* is a function of the temperature, which was also found in the literature. To solve this problem, we propose a Boltzmann fitting of the data and obtain T* at the high-temperature limit. With these characteristic parameters, the pressure-volume-temperature (PVT) data of iPP and PEOC are predicted by the SLLFT equation of state. To justify the correctness of our results, we also obtain the PVT data for iPP and PEOC by experiments. Good agreement is found between the two sets of data. By integrating the Euler-Lagrange equation and the Cahn-Hilliard relation, we predict the density profiles and the surface tensions for iPP and PEOC, respectively. Furthermore, a recursive method is proposed to obtain the characteristic interaction energy parameter between iPP and PEOC. This method, which does not require fitting to the experimental phase equilibrium data, suggests an alternative way to predict the phase diagrams that are not easily obtained in experiments. As an example, in the framework of SLLFT, the spinodal curve for the iPP/PEOC blend is predicted at the low molecular weights that are used in the simulations.

We develop a mesoscale nonequilibrium simulation model to study the effect of steady shear on the hierarchical self-assembly of soft disklike particles in dilute solutions. By properly tuning shear rates and solvent conditions, soft disklike particles can self-assemble into flexible threads and bundle-like structures along the flow direction. Shear flow facilitates the self-assembly of soft disklike particles into one-dimensional long threads along the flow direction; however, it suppresses the formation of flexible bundles from the threads while decreasing the solvent quality. The relatively well-defined bundle structures along the flow direction can only be obtained when the solvent condition becomes even worse. Our study elucidates how the solvent condition and shear rate can be utilized to control the shear-induced self-assembled structures, which would enable designed nanofabrication.

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.