Dissipative particle dynamics simulation study on the binary mixture phase separation coupled with polymerization

Coarse-Grained Modeling of Liquid-Liquid Phase Separation in Cells: Challenges and Opportunities

Liquid-liquid phase separation (LLPS) within cells gives rise to membraneless organelles, which play pivotal roles in numerous cellular functions. A comprehensive understanding of the functional aspects of intrinsically disordered protein (IDP) condensates necessitates elucidating their inherent structures and establishing correlations with biological functions. Coarse-grained (CG) molecular dynamics (MD) simulations present a promising avenue for gaining insights into LLPS mechanisms of biomacromolecules. Essential to this endeavor is the development of tailored CG force fields for MD simulations, incorporating the full spectrum of biomolecules involved in the formation of condensates and accounting for real-time biochemical reactions coupled to the LLPS. Moreover, developing accurate theoretical frameworks and establishing links between condensate structure and its function are imperative for a thorough comprehension of LLPS of biological systems.

Structural heterogeneity in tetra-armed gels revealed by computer simulation: Evidence from a graph theory assisted characterization

Designing homogeneous networks is considered one typical strategy for solving the problem of strength and toughness conflict of polymer network materials. Experimentalists have proposed the hypothesis of obtaining a structurally homogeneous hydrogel by crosslinking tetra-armed polymers, whose homogeneity was claimed to be verified by scattering characterization and other methods. Nevertheless, it is highly desirable to further evaluate this issue from other perspectives. In this study, a coarse-grained molecular dynamics simulation coupled with a stochastic reaction model is applied to reveal the topological structure of a polymer network synthesized by tetra-armed monomers as precursors. Two different scenarios, distinguished by whether internal cross-linking is allowed, are considered. We introduce the Dijkstra algorithm from graph theory to precisely characterize the network structure. The microscopic features of the network structure, e.g., loop size, dispersity, and size distribution, are obtained via the Dijkstra algorithm. By comparing the two reaction scenarios, Scenario II exhibits an overall more idealized structure. Our results demonstrate the feasibility of the Dijkstra algorithm for precisely characterizing the polymer network structure. We expect this work will provide a new insight for the evaluation and description of gel networks and further help to reveal the dynamic process of network formation.

Dissipative Particle Dynamics Simulation for Reaction-Induced Phase Separation of Thermoset/Thermoplastic Blends

Reaction-induced phase separation occurs during the curing reaction when a thermoplastic resin is dissolved in a thermoset resin, which enables toughening of the thermoset resin. As resin properties vary significantly depending on the morphology of the phase-separated structure, controlling the morphology formation is of critical importance. Reaction-induced phase separation is a phenomenon that ranges from the chemical reaction scale to the mesoscale dynamics of polymer molecules. In this study, we performed curing simulations using dissipative particle dynamics (DPD) coupled with a reaction model to reproduce reaction-induced phase separation. The curing reaction properties of the thermoset resin were determined by ab initio quantum chemical calculations, and the DPD parameters were determined by all-atom molecular dynamics simulations. This enabled mesoscopic simulations, including reactions that reflect the intrinsic material properties. The effects of the thermoplastic resin concentration, molecular weight, and curing conditions on the phase-separation morphology were evaluated, and the cure shrinkage and stiffness of each cured resin were confirmed to be consistent with the experimental trends. Furthermore, the local strain field under tensile deformation was visualized, and the inhomogeneous strain field caused by the phase-separated structures of two resins with different stiffnesses was revealed. These results can aid in understanding the toughening properties of thermoplastic additives at the molecular level.

Reaction-induced morphological transitions in a blend of diblock copolymers and reactive monomers: dissipative particle dynamics simulation

The dissipative particle dynamics (DPD) method is applied to the morphological transitions of microphase-separated domains in a mixture of symmetric AB-diblock copolymers and reactive C-monomers, where polymerization and cross-linking reactions take place among C-monomers. The initial structure for the DPD simulation is an equilibrated cylindrical domain structure prepared by the density-biased Monte Carlo method with density profiles obtained from the self-consistent field theory. By introducing a cross-linking reaction among reactive C-monomers, we confirmed that the DPD simulation reproduces the morphological transitions observed in experiments, where the domain morphology changes due to segregation between A-blocks of diblock copolymers and cross-linking networks of C-monomers. When the cross-linking reaction of C-monomers is sufficiently fast compared to the deformation of the domains, the initial cylindrical domains are preserved, while the distance between the domains increases. On the other hand, when the formation of the cross-linking network is slow, the domains can deform and reconnect with each other in the developing cross-linking network. In this case, we observe morphological transitions from the initial domain morphology with a large-curvature interface to another domain morphology with a smaller-curvature interface, such as the transition from the cylindrical phase to the lamellar phase. We calculated the spatial correlations in the microphase-separated domains and found that such correlations are affected by the speed of the formation of the cross-linking network depending on whether the bridging between microphase-separated domains occurs in a nucleation and growth process or in a spinodal decomposition process.

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.

Understanding of supramolecular solution polymerization and interfacial polymerization via forming multiple hydrogen bonds: a computer simulation study

By employing dissipative particle dynamics (DPD) simulations combined with stochastic polymerization models, we have conducted a detailed simulation study of supramolecular solution polymerization as well as interfacial polymerization employing a coarse-grained model which is closer to the real monomer structure. By adding bending angle potentials to coarse-grained models representing supramolecular reactive monomers, we achieved monomer model simulations for different kinds of multiple hydrogen bonds. Our simulation results indicated that for the interfacial polymerization system, the volume of the monomer caused a strong steric hindrance effect, which in turn led to a low average degree of polymerization of the product. Therefore, by appropriately reducing the volume of the reaction monomer (corresponding to different confinement ascribed to the multiple hydrogen bonds), the average polymerization degree, the degree of reaction and the polymerization rate of the monomer can be effectively improved. For the solution polymerization system and the interfacial polymerization system, a certain proportion of rigid monomers and flexible monomers (60% rigid monomers and 40% flexible monomers) are mixed. High molecular weight products can thus be obtained via the polymerization reaction. The simulation strategy proposed in this study can not only provide theoretical guidance for better design of new supramolecular systems, but also provide ideas for the further synthesis of higher molecular weight supramolecular polymers.

Understanding of supramolecular emulsion interfacial polymerization in silico

The composition and structure of a membrane determine its functionality and practical application. We study the supramolecular polymeric membrane prepared by supramolecular emulsion interfacial polymerization (SEIP) on the oil-in-water droplet via the computer simulation method. The factors that may influence its structure and properties are investigated, such as the degree of polymerization and molecular weight distribution (MWD) of products in the polymeric membranes. We find that the SEIP can lead to a higher total degree of polymerization as compared to the supramolecular interfacial polymerization (SIP). However, the average chain length of products in the SEIP is lower than that of the SIP due to its obvious interface curvature. The stoichiometric ratio of reactants in two phases will affect the MWD of the products, which further affects the performance of the membranes in practical applications, such as drug release rate and permeability. Besides, the MWD of the product by SEIP obviously deviates from the Flory distribution as a consequence of the curvature of reaction interface. In addition, we obtain the MWD for the emulsions whose size distribution conforms to the Gaussian distribution so that the MWD may be predicted according to the corresponding emulsion size distribution. This study helps us to better understand the controlling factors that may affect the structure and properties of supramolecular polymeric membranes by SEIP.

Polymerization-Induced Reassembly of Gemini Molecules toward Generating Porous Two-Dimensional Polymers

In situ polymerization of preorganized amphiphilic monomers on various substrates provides a flexible synthetic route to construct high-quality two-dimensional polymers (2DPs) with designed functionalities. However, the detailed polymerization kinetics of these monomers in 2D confinement and their impact on the structural features of 2DPs have not been efficiently explored. Here, using dissipative particle dynamics (DPD) simulations, we unveil the similarity of the polymerization kinetics of the amphiphilic Gemini molecules in both a 2D-confined space and solution and emphasize the key role of the initiator concentration in modifying the morphology of 2DPs. More interestingly, introducing a spacer group into the Gemini monomer facilitates the formation of porous 2DPs. The size and periodic arrangement of pores in these 2DPs could be directly controlled by the Gemini molecular geometries and polymerization kinetics. The insights based on our DPD simulations provide valuable guidelines for the rational design and synthesis of 2DPs from a wider range of amphiphilic molecules.

Microscopic characteristics of Janus nanoparticles prepared via a grafting-from reaction at the immiscible liquid interface

The dynamic process of synthesizing Janus nanoparticles (JNPs) at a water/oil two-phase interface using a grafting-from reaction is investigated via dissipative particle dynamics simulations. We find that the interfacial tension, the initial monomer concentration, and the reaction probability can greatly influence the microscopic characteristics of JNP structure. It is difficult to synthesize a symmetric JNP with an equal volume ratio between hydrophilic and hydrophobic parts by grafting-from methods unless the physical chemical conditions in the two phases are strictly symmetric, and there is always a disordered domain on the JNP at a two immiscible solvents interface. Interestingly, for certain routes for synthesizing JNPs with a grafting-from method, the higher interfacial tension between the water and oil phases may enhance the degree of disorder of the grafted chains. The asymmetric initial monomer concentration in solution and the reaction probability can be used to control the syntheses of asymmetric JNPs.

Tethering solvophilic blocks to the ends of polymer brushes: an effective method for adjusting surface patterns

The effect of different lengths of solvophilic A and C blocks on the assembled configuration of intermediate solvophobic B-blocks in both ABA and ABC polymer brush systems is investigated via dissipative particle dynamics simulations. For the AB diblock copolymer brush with solvophilic A-blocks being grafted to the surface, B-blocks self-assemble into spherical micelle structures that are immersed in a layer formed by the A-blocks. Tethering a very small solvophilic block A(C) at the free end of the polymer brush pulls the B-blocks toward the polymer brush/solvent interface and increases their local density which can significantly change the B-block self-assembled structure from spherical micelles to ripples. By increasing the length of the outermost solvophilic blocks, the lateral density distribution of B-blocks can be further changed, resulting in the domain size of the ripple structure first decreasing and then increasing. Compared to the ABA system, the incompatibility between the A and C blocks can effectively reduce the vertical domain separation caused by the fusion of the upper and lower A blocks. Then, based on an AB diblock copolymer brush system with self-assembled spherical micelles, we introduce extremely short free solvophilic blocks A(C) in dilute solution that can be tethered to the free ends of the polymer brush by using a reaction model [Liu et al., J. Chem. Phys., 2007, 127, 144903]. We find that the micelles' coalescence is mainly affected by the content of tethered reactive solvophilic blocks, and only weakly affected by the reaction rate of the reversible reactions. This strategy of tethering solvophilic blocks to the ends of polymer brushes can be an effective way for the fabrication of stimuli-responsive surfaces and for adjusting nanoscopic surface patterns.

Simultaneous polymer chain growth with the coexistence of bulk and surface initiators: insight from computer simulations

By Brownian dynamics simulations we study the simultaneous polymer chain growth process with the coexistence of bulk and surface initiators. We find that when the surface initiator density is low enough, the practical experimental way to estimate the dispersity (Đ) of surface-initiated chains on the basis of the dispersity of bulk-initiated chains remains valid as long as the conformations of grafted chains remain within the mushroom regime (i.e., the grafted chains are sparsely distributed). On the other hand, although the average chain lengths of surface and bulk polymers could be equivalent when certain conditions are met, their mass distributions are still different. We also find that increasing the fraction of surface initiators leads to an enlarged disparity in Đ and average length between surface and bulk chains, which is inconsistent with previous studies. This study helps in better understanding the cooperative competition and suppressing effect of bulk chains on surface grown chains, as well as the cause of the dispersity of the surface grown chains as compared to their bulk counterparts with the coexistence of bulk and surface initiators.

Polymer-grafted nanoparticles prepared via a grafting-from strategy: a computer simulation study

Nanoparticles (NPs) grafted with polymer chains prepared via a grafting-from strategy are studied through coarse-grained molecular dynamics simulations combined with our stochastic reaction model. A system involving multiple individual NPs, with grafting-from processes for all the NPs induced simultaneously, is simulated, so that chain growth competition on the same NP, as well as between neighbouring NPs, are both naturally considered. Our results imply that there should be an optimized range of NP sizes, as compared to monomer size, in which initiator sites are most easily induced. Besides, when the initiator density is high, a shielding effect from the sparse long chains on the most short chains or initiators evidently yields an extremely unbiased distribution of chains. We also adopt a representative polymer-tethered NP prepared via a grafting-from strategy to study the potential of mean force between NPs, so that the dispersion and stabilization abilities of such polymer-grafted NPs in a polymer matrix can be generally predicted during the preparation of polymer nanocomposite materials. Our study helps to elucidate the cause of chain dispersity during the grafting-from process and could act as a guide for better design and to improve the performance of polymer nanocomposites.

Micellar polymerization: Computer simulations by dissipative particle dynamics

Nowadays, micellar polymerization is widely used in different fields of industry and research, including modern living polymerization technique. However, this process has many variables and there is no comprehensive model to describe all features. This research presents simulation methodology which describes key properties of such reactions to take a guide through a variety of their modifications. Dissipative particle dynamics is used in addition to Monte Carlo scheme to simulate initiation, propagation, and termination events. Influence of initiation probability and different termination processes on final conversion and molecular-weight distribution are presented. We demonstrate that prolonged initiation leads to increasing in polymer average molecular weight, and surface termination events play major role in conversion limitation, in comparison with recombination. © 2018 Wiley Periodicals, Inc.

Polymerization-induced polymer aggregation or polymer aggregation-enhanced polymerization? A computer simulation study

In this study, using dissipative particle dynamics simulations coupled with the stochastic reaction model, we investigate the polymerization-induced polymer aggregation process and the polymer aggregation-enhanced polymerization process in a binary solution.

Dissipative particle dynamics study of phase separation in binary fluid mixtures in periodic and confined domains

We investigate the phase separation behavior of binary mixtures in two-dimensional periodic and confined domains using dissipative particle dynamics. Two canonical problems of fluid mechanics are considered for the confined domains: square cavity with no-slip walls and lid-driven cavity with one driven wall. The dynamics is studied for both weakly and strongly separating mixtures and different area fractions. The phase separation process is analyzed using the structure factor and the total interface length. The dynamics of phase separation in the square cavity and lid-driven cavity are observed to be significantly slower when compared to the dynamics in the periodic domain. The presence of the no-slip walls and the inertial effects significantly influences the separation dynamics. Finally, we show that the growth exponent for the strongly separating case is invariant to changes in the inter-species repulsion parameter.

A general strategy to synthesize chemically and topologically anisotropic Janus particles

Emulsion polymerization is the most widely used synthetic technique for fabricating polymeric particles. The interfacial tension generated with this technique limits the ability to tune the topology and chemistry of the resultant particles. We demonstrate a general emulsion interfacial polymerization approach that involves introduction of additional anchoring molecules surrounding the microdroplets to synthesize a large variety of Janus particles with controllable topological and chemical anisotropy. This strategy is based on interfacial polymerization mediated by an anchoring effect at the interface of microdroplets. Along the interface of the microdroplets, the diverse topology and surface chemistry features of the Janus particles can be precisely tuned by regulating the monomer type and concentration as well as polymerization time. This method is applicable to a wide variety of monomers, including positively charged, neutrally charged, and negatively charged monomers, thereby enriching the community of Janus particles.

Generating minimal living systems from non-living materials and increasing their evolutionary abilities

We review lessons learned about evolutionary transitions from a bottom-up construction of minimal life. We use a particular systemic protocell design process as a starting point for exploring two fundamental questions: (i) how may minimal living systems emerge from non-living materials? and (ii) how may minimal living systems support increasingly more evolutionary richness? Under (i), we present what has been accomplished so far and discuss the remaining open challenges and their possible solutions. Under (ii), we present a design principle we have used successfully both for our computational and experimental protocellular investigations, and we conjecture how this design principle can be extended for enhancing the evolutionary potential for a wide range of systems.This article is part of the themed issue 'The major synthetic evolutionary transitions'.

GALAMOST: GPU-accelerated large-scale molecular simulation toolkit

GALAMOST [graphics processing unit (GPU)-accelerated large-scale molecular simulation toolkit] is a molecular simulation package designed to utilize the computational power of GPUs. Besides the common features of molecular dynamics (MD) packages, it is developed specially for the studies of self-assembly, phase transition, and other properties of polymeric systems at mesoscopic scale by using some lately developed simulation techniques. To accelerate the simulations, GALAMOST contains a hybrid particle-field MD technique where particle–particle interactions are replaced by interactions of particles with density fields. Moreover, the numerical potential obtained by bottom-up coarse-graining methods can be implemented in simulations with GALAMOST. By combining these force fields and particle-density coupling method in GALAMOST, the simulations for polymers can be performed with very large system sizes over long simulation time. In addition, GALAMOST encompasses two specific models, that is, a soft anisotropic particle model and a chain-growth polymerization model, by which the hierarchical self-assembly of soft anisotropic particles and the problems related to polymerization can be studied, respectively. The optimized algorithms implemented on the GPU, package characteristics, and benchmarks of GALAMOST are reported in detail.

DISSIPATIVE PARTICLE DYNAMICS SIMULATION STUDY ON CONTROLLING MOLECULAR WEIGHT DISTRIBUTION IN EMULSION POLYMERIZATION

The controlling factors on molecular weight distribution in emulsion polymerization are investigated using dissipative particle dynamics simulations. The propagation and bi-radical termination kinetic steps are taken into account in the simulations by coupling a Monte Carlo type reaction model. We find that monomer concentration can be very efficient on controlling molecular weight distribution and plays a decisive role on the formation of high molecular weight polymers. Increasing initiator concentration can effectively reduce the polymer molecular weight while increasing the polydispersity index. Moreover, increasing polymerization rate can slightly narrow the molecular weight distribution. We also find that, by suitably tuning the surfactant chain length, it may be possible to obtain an optimal molecular weight distribution in emulsion polymerization.

DISSIPATIVE PARTICLE DYNAMICS SIMULATION ON THE PREPARATION PROCESS OF MACROPOROUS STYRENE-DIVINYLBENZENE COPOLYMER

Dissipative particle dynamics (DPD) was carried out to study the copolymerization process of halogen-substituted styrene (ST-X), divinylbenzene (DVB), and dodecane (DO) in aqueous solution. The results showed that properly selecting the sizes of DO and water beads, and choosing the appropriate interaction parameters between beads are crucial in the simulation system. The influence of polymerization parameters, such as polymerization temperature, ST-X/DVB/DO/water ratio, and halogen substitution on ST on the copolymer morphology, were investigated in detail and confirmed that the simulation results can reproduce the macroscopic experimental phase separation. In addition, the pore volume of copolymers formed at different polymerization conditions is estimated, which consisted with the data measured in the experiments. Our work has demonstrated that DPD methods can be applied to study such copolymerization process.