Monte Carlo Simulation of Surface-Initiated Polymerization: Heterogeneous Reaction Environment

Modular addition strategy-regulated polymerization-induced self-assembly: an in silico experiment

We propose a modular addition strategy-regulated polymerization-induced self-assembly (PISA) system to effectively control the reaction kinetics and self-assembly morphologies. We validated this strategy by performing in silico experiments on a well-established PISA system. Two categories of modular addition strategies, i.e., the multistep addition strategy and the constant rate addition strategy, were investigated. Results showed that the modular addition operation of macromolecular chain transfer agents (macro-CTAs) effectively regulated the width of the molecular weight distribution for the hydrophobic PSt block, which further led to an assembly of vesicle structures with irregular aspherical cavities. Besides, we found a new transition pathway for the formation of vesicles, which involved generation of small vesicles in the early stage followed by a gradual growth in the intermediate and late stages. In the constant rate addition strategy, with the increase in the addition rate of macro-CTA, we found that the morphology basically tended to change from a micellar structure to a vesicle structure. This study holds potential to inspire future work toward the improvement of experimental techniques in PISA-relevant systems.

Surface-Initiated Polymerization with an Initiator Gradient: A Monte Carlo Simulation

Due to the difficulty of accurately characterizing properties such as the molecular weight (Mn) and grafting density (σ) of gradient brushes (GBs), these properties are traditionally assumed to be uniform in space to simplify analysis. Applying a stochastic reaction model (SRM) developed for heterogeneous polymerizations, we explored surface-initiated polymerizations (SIPs) with initiator gradients in lattice Monte Carlo simulations to examine this assumption. An initial exploration of SIPs with 'homogeneously' distributed initiators revealed that increasing σ slows down the polymerization process, resulting in polymers with lower molecular weight and larger dispersity (Đ) for a given reaction time. In SIPs with an initiator gradient, we observed that the properties of the polymers are position-dependent, with lower Mn and larger Đ in regions of higher σ, indicating the non-uniform properties of polymers in GBs. The results reveal a significant deviation in the scaling behavior of brush height with σ compared to experimental data and theoretical predictions, and this deviation is attributed to the non-uniform Mn and Đ.

Molecular Dynamics Simulations of Ideal Living Polymerization: Terminal Model and Kinetic Aspects

Living polymerization is an important synthetic approach to achieving precise control of synthesized polymers, which is crucial for their applications. The molecular weight distribution (MWD) prescribes the macroscopic properties of polymers and hence is a key feature to characterize polymerization. In this work, we present a systematic molecular dynamics simulation study of ideal living polymerization in bulk and surface-initiated systems based on a terminal stochastic reaction model. The evolution of polymer dispersity and MWD along with the polymerization process is examined. We demonstrate that MWD is generally well captured by the Schulz-Zimm distribution for bulk and surface-initiated systems with low grafting densities. However, as the grafting density in the surface-initiated case increases, heterogeneity in chain growth emerges due to the kinetic trapping of reactive sites, which causes the starving of short chains and the thriving of minority long chains such that a shoulder region shows up in MWD. This effect can be enhanced by kinetic compressing induced by polymerization. In addition, the interplay of bonding reaction kinetics and other kinetic properties (e.g., mass transfer and polymer relaxation) is further explored, alongside the influences of bonding probability and reactant concentration. We expect that this investigation will aid in our understanding of typical kinetic aspects of living polymerization.

A Simple Stochastic Reaction Model for Heterogeneous Polymerizations

The stochastic reaction model (SRM) treats polymerization as a pure probability‐based issue, which is widely applied to simulate various polymerization processes. However, in many studies, active centers were assumed to react with the same probability, which cannot reflect the heterogeneous reaction microenvironment in heterogeneous polymerizations. Recently, we have proposed a simple SRM, in which the reaction probability of an active center is directly determined by the local reaction microenvironment. In this paper, we compared this simple SRM with other SRMs by examining living polymerizations with randomly dispersed and spatially localized initiators. The results confirmed that the reaction microenvironment plays an important role in heterogeneous polymerizations. This simple SRM provides a good choice to simulate various polymerizations.