Polymer brush relaxation during and after polymerization – Monte Carlo simulation study

The Kinetics of Polymer Brush Growth in the Frame of the Reaction Diffusion Front Formalism

We studied the properties of a reaction front that forms in irreversible reaction-diffusion systems with concentration-dependent diffusivities during the synthesis of polymer brushes. A coarse-grained model of the polymerization process during the formation of polymer brushes was designed and investigated for this purpose. In this model, a certain amount of initiator was placed on an impenetrable surface, and the "grafted from" procedure of polymerization was carried out. The system consisted of monomer molecules and growing chains. The obtained brush consisted of linear chains embedded in nodes of a face-centered cubic lattice with excluded volume interactions only. The simulations were carried out for high rafting densities of 0.1, 0.3, and 0.6 and for reaction probabilities of 0.02, 0.002, and 0.0002. Simulations were performed by means of the Monte Carlo method while employing the Dynamic Lattice Liquid model. Some universal behavior was found, i.e., irrespective of reaction rate and grafting density, the width of the reaction front as well as the height of the front show for long times the same scaling with respect to time. During the formation of the polymer layer despite the observed difference in dispersion of chain lengths for different grafting densities and reaction rates at a given layer height, the quality of the polymer layer does not seem to depend on these parameters.

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.

Polymerization and Structure of Opposing Polymer Brushes Studied by Computer Simulations

A model of the polymerization process during the formation of a pair of polymer brushes was designed and investigated. The obtained system consisted of two impenetrable parallel surfaces with the same number of chains grafted on both surfaces. Coarse-grained chains embedded in nodes of a face-centered cubic lattice with excluded volume interactions were obtained by a ‘grafted from’ procedure. The structure of synthesized macromolecular systems was also studied. Monte Carlo simulations using the dynamic lattice liquid model were employed using dedicated parallel machine ARUZ in a large size and time scale. The parameters of the polymerization process were found to be crucial for the proper structure of the brush. It was found that for high grafting densities, chains were increasingly compressed, and there is surprisingly little interpenetration of chains from opposite surfaces. It was predicted and confirmed that in a polydisperse sample, the longer chains have unique configurations consisting of a stretched stem and a coiled crown.

Dynamics of Opposing Polymer Brushes: A Computer Simulation Study

Opposing polymer brush systems were synthesized and investigated by molecular modeling. Chains were restricted to a face-centered cubic lattice with the excluded volume interactions only. The system was confined between two parallel impenetrable walls, with the same number of chains grafted to each surface. The dynamic properties of such systems were studied by Monte Carlo simulations based on the dynamic lattice liquid model and using a highly efficient parallel machine ARUZ, which enabled the study of large systems and long timescales. The influence of the surface density and mean polymer length on the system dynamic was discussed. The self-diffusion coefficient of the solvent depended strongly on the degree of polymerization and on the polymer concentration. It was also shown that it is possible to capture changes in solvent mobility that can be attributed to the regions of different polymer densities.

The Competition of Termination and Shielding to Evaluate the Success of Surface-Initiated Reversible Deactivation Radical Polymerization

One of the challenges for brush synthesis for advanced bioinspired applications using surface-initiated reversible deactivation radical polymerization (SI-RDRP) is the understanding of the relevance of confinement on the reaction probabilities and specifically the role of termination reactions. The present work puts forward a new matrix-based kinetic Monte Carlo platform with an implicit reaction scheme capable of evaluating the growth pattern of individual free and tethered chains in three-dimensional format during SI-RDRP. For illustration purposes, emphasis is on normal SI-atom transfer radical polymerization, introducing concepts such as the apparent livingness and the molecular height distribution (MHD). The former is determined based on the combination of the disturbing impact of termination (related to conventional livingness) and shielding of deactivated species (additional correction due to hindrance), and the latter allows structure-property relationships to be identified, starting at the molecular level in view of future brush characterization. It is shown that under well-defined SI-RDRP conditions the contribution of (shorter) hindered dormant chains is relevant and more pronounced for higher average initiator coverages, despite the fraction of dead chains being less. A dominance of surface-solution termination is also put forward, considering two extreme diffusion modes, i.e., translational and segmental. With the translational mode termination is largely suppressed and the living limit is mimicked, whereas with the segmental mode termination occurs more and the termination front moves upward alongside the polymer layer growth. In any case, bimodalities are established for the tethered chains both on the level of the chain length distribution and the MHD.