A dissipative particle dynamics simulation study on phase diagrams for the self-assembly of amphiphilic hyperbranched multiarm copolymers in various solvents

Single core and multicore aggregates from a polymer mixture: A dissipative particle dynamics study

HYPOTHESIS

Multicore block copolymer aggregates correspond to self-assembly such that the polymer system spontaneously phase separates to multiple, droplet-like cores differing in the composition from the polymer surroundings. Such multiple core aggregates are highly useful capsules for different applications, e.g., drug transport, catalysis, controlled solvation, and chemical reactions platforms. We postulate that polymer system composition provides a direct means for designing polymer systems that self-assemble to such morphologies and controlling the assembly response.

SIMULATIONS

Using dissipative particle dynamics (DPD) simulations, we examine the self-assembly of a mixture of highly and weakly solvophobic homopolymers and an amphiphilic block copolymer in the presence of solvent. We map the multicore vs single core (core-shell particles) assembly response and aggregate structure in terms of block copolymer concentration, polymer component ratios, and chain length of the weakly solvophobic homopolymer.

FINDINGS

For fixed components and polymer chemistries, the amount of block copolymer is the key to controlling single core vs multicore aggregation. We find a polymer system dependent critical copolymer concentration for the multicore aggregation and that a minimum level of incompatibility between the solvent and the weakly solvophobic component is required for multicore assembly. We discuss the implications for polymer system design for multicore assemblies. In summary, the study presents guidelines to produce multicore aggregates and to tune the assembly from multicore aggregation to single core core-shell particles.

Implicit-Solvent Coarse-Grained Simulations of Linear-Dendritic Block Copolymer Micelles

The design of nanoassemblies can be conveniently achieved by tuning the strength of the hydrophobic interactions of block copolymers in selective solvents. These block copolymer micelles form supramolecular aggregates, which have attracted great attention in the area of drug delivery and imaging in biomedicine due to their easy-to-tune properties and straightforward large-scale production. In the present work, we have investigated the micellization process of linear-dendritic block copolymers in order to elucidate the effect of branching on the micellar properties. We focus on block copolymers formed by linear hydrophobic blocks attached to either dendritic neutral or charged hydrophilic blocks. We have implemented a simple protocol for determining the equilibrium micellar size, which permits the study of linear-dendritic block copolymers in a wide range of block morphologies in an efficient and parallelizable manner. We have explored the impact of different topological and charge properties of the hydrophilic blocks on the equilibrium micellar properties and compared them to predictions from self-consistent field theory and scaling theory. We have found that, at higher degrees of branching in the corona and for short polymer chains, excluded volume interactions strongly influence the micellar aggregation as well as their effective charge.

Hybrid Molecules Consisting of Lysine Dendrons with Several Hydrophobic Tails: A SCF Study of Self-Assembling

In this article, we used the numerical self-consistent field method of Scheutjens-Fleer to study the micellization of hybrid molecules consisting of one polylysine dendron with charged end groups and several linear hydrophobic tails attached to its root. The main attention was paid to spherical micelles and the determination of the range of parameters at which they can appear. A relationship has been established between the size and internal structure of the resulting spherical micelles and the length and number of hydrophobic tails, as well as the number of dendron generations. It is shown that the splitting of the same number of hydrophobic monomers from one long tail into several short tails leads to a decrease in the aggregation number and, accordingly, the number of terminal charges in micelles. At the same time, it was shown that the surface area per dendron does not depend on the number of hydrophobic monomers or tails in the hybrid molecule. The relationship between the structure of hybrid molecules and the electrostatic properties of the resulting micelles has also been studied. It is found that the charge distribution in the corona depends on the number of dendron generations G in the hybrid molecule. For a small number of generations (up to G=3), a standard double electric layer is observed. For a larger number of generations (G=4), the charges of dendrons in the corona are divided into two populations: in the first population, the charges are in the spherical layer near the boundary between the micelle core and shell, and in the second population, the charges are near the periphery of the spherical shell. As a result, a part of the counterions is localized in the wide region between them. These results are of potential interest for the use of spherical dendromicelles as nanocontainers for drug delivery.

Mesoscale Modeling of Agglomeration of Molecular Bottlebrushes: Focus on Conformations and Clustering Criteria

Using dissipative particle dynamics, we characterize dynamics of aggregation of molecular bottlebrushes in solvents of various qualities by tracking the number of clusters, the size of the largest cluster, and an average aggregation number. We focus on a low volume fraction of bottlebrushes in a range of solvents and probe three different cutoff criteria to identify bottlebrushes belonging to the same cluster. We demonstrate that the cutoff criteria which depend on both the coordination number and the length of the side chain allows one to correlate the agglomeration status with the structural characteristics of bottlebrushes in solvents of various qualities. We characterize conformational changes of the bottlebrush within the agglomerates with respect to those of an isolated bottlebrush in the same solvents. The characterization of bottlebrush conformations within the agglomerates is an important step in understanding the relationship between the bottlebrush architecture and material properties. An analysis of three distinct cutoff criteria to identify bottlebrushes belonging to the same cluster introduces a framework to identify both short-lived transient and long-lived agglomerates; the same approach could be further extended to characterize agglomerates of various macromolecules with complex architectures beyond the specific bottlebrush architecture considered herein.

Passive and active tracer dynamics in polymer solutions with isotropic-to-nematic phase transition

Macromolecular crowding plays a crucial role in determining the dynamics in a living cell. We adopt Langevin dynamics simulations to investigate the anomalous diffusion dynamics of passive and active particles in a solution of polymer chains with tunable stiffness. The solution's anisotropic feature is modulated by changing both the polymer stiffness and volume fraction, where isotropic-to-nematic phase transition is involved. Our results demonstrate the significant impact of polymer flexibility on the dynamics of both passive and active probes. The distinct diffusion mechanism for an active particle is clarified by the interplay between polymer stiffness, crowdedness and activity. Polymer stiffness leads to a global inhibition effect on passive particle diffusion. The diffusion coefficient exhibits an intriguing non-monotonic variation at increasing polymer stiffness, which is due to the fact that the alignment of polymer chains is beneficial for diffusion along the nematic direction but unfavorable for that in the direction perpendicular to it. In sharp contrast, polymer stiffness plays a dominant role in facilitating active particle diffusion. Self-propulsion of the particle can utilize stiffness-induced elastic interactions more efficiently, which promotes its mobility in both directions. Meanwhile, an active particle might have a stronger ability to take advantage of the polymer alignment, contributing substantially enhanced diffusivity. In addition, the diffusion coefficient of an active particle is subject to a tendency of degeneration against varying volume fraction. This counter-intuitive behavior is due to the contrasting factors that increasing crowdedness induces a lower particle speed but a longer persistent motion time.

DPD Modelling of the Self- and Co-Assembly of Polymers and Polyelectrolytes in Aqueous Media: Impact on Polymer Science

This review article is addressed to a broad community of polymer scientists. We outline and analyse the fundamentals of the dissipative particle dynamics (DPD) simulation method from the point of view of polymer physics and review the articles on polymer systems published in approximately the last two decades, focusing on their impact on macromolecular science. Special attention is devoted to polymer and polyelectrolyte self- and co-assembly and self-organisation and to the problems connected with the implementation of explicit electrostatics in DPD numerical machinery. Critical analysis of the results of a number of successful DPD studies of complex polymer systems published recently documents the importance and suitability of this coarse-grained method for studying polymer systems.

Self-assembly of miktoarm palm tree-like star copolymers in a selective solvent

Amphiphilic miktoarm star copolymers with one long solvophobic arm (a "stem") and several short solvophilic arms (the "leaves") were studied in a selective solvent using mesoscopic computer simulations. The conventional morphologies (spherical, cylindrical and vesicular) as well as the mixed ones were obtained. However, the resulting diagram of states appeared to be different from the diagram of the linear diblock copolymer with the analogous composition. Namely, the increase of the number of leaves at fixed solvophobic-solvophilic ratio leads to the transition from the vesicles to the cylinders, while the latter ones eventually transform into spherical micelles in the case of highly branched copolymers. The observed effect appears due to the increase of the interfacial area between the collapsed and swollen blocks per single macromolecule. In turn, the increase of the solvent selectivity shifts the stability region of the cylindrical micelles to the region of more symmetric copolymer composition. Meanwhile, the compatibility between the blocks has a weak effect on the resulting morphology. Finally, it was found that the increase in the number of leaves and the simultaneous decrease in their length results in the localization of higher amount of solvophilic segments near the core-solvent interface, which in the case of cylindrical micelles significantly affects the shape of the aggregates making them thinner and longer.

Toward enzyme-responsive polymersome drug delivery

In drug delivery, enzyme-responsive drug carriers are becoming increasingly relevant because of the growing association of disease pathology with enzyme overexpression. Polymersomes are of interest to such applications because of their tunable properties. While polymersomes open up a wide range of chemical and physical properties to explore, they also present a challenge in developing generalized rules for the synthesis of novel systems. Motivated by this issue, in this perspective, we summarize the existing knowledge on enzyme-responsive polymersomes and outline the main design choices. Then, we propose heuristics to guide the design of novel systems. Finally, we discuss the potential of an integrated approach using computer simulations and experimental studies to streamline this design process and close the existing knowledge gaps.

Advances in self-assembling of pH-sensitive polymers: A mini review on dissipative particle dynamics

Dissipative Particle Dynamics (DPD) is a mesoscopic simulation program used to simulate the behavior of complex fluids. This work systematically reviews the use of DPD to simulate the self-assembly process of pH-sensitive drug-loaded nanoparticles. pH-sensitive drug-loaded nanoparticles have the characteristics of good targeting and slow release in the body, which is an ideal method for treating cancer and other diseases. As an excellent simulation method, DPD can help people explore the loading and release laws of drugs with complex molecular structures and has extensive applications in other medical fields. This article reviews the self-assembly process of pH-sensitive polymers under neutral conditions and explores the factors that affect the self-assembly structure. It points out that different hydrophilic-hydrophobic ratios, molecular structures, and component distributions will affect the morphology, stability and drug carrying capacity of micelles. This article also introduces the release mechanism of the drug in detail and introduces the factors that affect the release. This article can help relevant researchers to follow the latest advances in the DPD simulation and pH-sensitive drug nano-carrier and insight people to investigate the further application of DPD simulation in biomedical science.

Effects of active crowder size and activity-crowding coupling on polymer translocation

Polymer translocation in complex environments is crucially important to many biological processes in life. In the present work, we adopted two-dimensional Langevin dynamics simulations to study the forced and unbiased polymer translocation dynamics in active and crowded media. The translocation time and probability are analyzed in terms of active force F_a, volume fraction φ and also the crowder size. The non-trivial active crowder size effect and activity-crowding coupling effect as well as the novel mechanism of unbiased translocation between two active environments with different active particle sizes are clarified. Firstly, for forced translocation, we reveal an intriguing non-monotonic dependence of the translocation time on the crowder size in the case of large activity. In particular, crowders of intermediate size similar to the polymer segment are proven to be the most favorable for translocation. Moreover, a facilitation-inhibition crossover of the translocation time with increasing volume fraction is observed, indicating a crucial activity-crowding coupling effect. Secondly, for unbiased translocation driven by different active crowder sizes, the translocation probability demonstrates a novel turnover phenomenon, implying the appearance of an opposite directional preference as the active force exceeds a critical value. The translocation time in both directions decreases monotonically with the active force. The asymmetric activity effect together with the entropic driving scenario provides a reasonable picture for the peculiar behavior observed in unbiased translocation.

Tunable depletion force in active and crowded environments

We adopt two-dimensional Langevin dynamics simulations to study the effective interactions between two passive colloids in a bath crowded with active particles. We mainly pay attention to the significant effects of active particle size, crowding-activity coupling, and chirality. First, a transition of depletion force from repulsion to attraction is revealed by varying particle size. Moreover, larger active crowders with sufficient activity can generate strong attractive force, which is in contrast to the cage effect in passive media. It is interesting that the attraction induced by large active crowders follows a linear scaling with the persistence length of active particles. Second, the effective force also experiences a transition from repulsion to attraction as volume fraction increases, as a consequence of the competition between the two contrastive factors of activity and crowding. As bath volume fraction is relatively small, activity generates a dominant repulsion force, while as the bath becomes concentrated, crowding-induced attraction becomes overwhelming. Lastly, in a chiral bath, we observe a very surprising oscillation phenomenon of active depletion force, showing an evident quasiperiodic variation with increasing chirality. Aggregation of active particles in the vicinity of the colloids is carefully examined, which serves as a reasonable picture for our observations. Our findings provide an inspiring strategy for the tunable active depletion force by crowding, activity, and chirality.

Molecular Modelling Guided Modulation of Molecular Shape and Charge for Design of Smart Self-Assembled Polymeric Drug Transporters

Nanomedicine employs molecular materials for prevention and treatment of disease. Recently, smart nanoparticle (NP)-based drug delivery systems were developed for the advanced transport of drug molecules. Rationally engineered organic and inorganic NP platforms hold the promise of improving drug targeting, solubility, prolonged circulation, and tissue penetration. However, despite great progress in the synthesis of NP building blocks, more interdisciplinary research is needed to understand their self-assembly and optimize their performance as smart nanocarriers. Multi-scale modeling and simulations provide a valuable ally to experiment by mapping the potential energy landscape of self-assembly, translocation, and delivery of smart drug-loaded NPs. Here, we highlight key recent advances to illustrate the concepts, methods, and applications of smart polymer-based NP drug delivery. We summarize the key design principles emerging for advanced multifunctional polymer topologies, illustrating how the unusual architecture and chemistry of dendritic polymers, self-assembling polyelectrolytes and cyclic polymers can provide exceptional drug delivery platforms. We provide a roadmap outlining the opportunities and challenges for the effective use of predictive multiscale molecular modeling techniques to accelerate the development of smart polymer-based drug delivery systems.

Using reactive dissipative particle dynamics to understand local shape manipulation of polymer vesicles

Biological cells have long been of interest to researchers due to their capacity to actively control their shape. Accordingly, there is significant interest in generating simplified synthetic protocells that can alter their shape based on an externally or internally generated stimulus. To date, most progress has been made towards controlling the global shape of a protocell, whereas less is known about generating a local shape change. Here, we seek to better understand the possible mechanisms for producing local morphological changes in a popular protocell system, the block copolymer vesicle. Accordingly, we have combined Dissipative Particle Dynamics (DPD) and the Split Reactive Brownian Dynamics algorithm (SRBD) to produce a simulation tool that is capable of modeling the dynamics of self-assembled polymer structures as they undergo chemical reactions. Using this Reactive DPD or RDPD method, we investigate local morphological change driven by either the microinjection of a stimulus or an enzymatically-produced stimulus. We find that sub-vesicle-scale morphological change can be induced by either a solvent stimulus that swells the vesicle membrane, or by a reactant stimulus that alters the chemistry of the block polymer in the membrane corona. Notably, the latter method results in a more persistent local deformation than the former, which we attribute to the slower diffusion of polymer chains relative to the solvent. We quantify this deformation and show that it can be modulated by altering the interaction parameter of the parts of the polymer chain that are affected by the stimulus.

Predicting asymmetric phospholipid microstructures in solutions

Asymmetric phospholipid microstructures, such as asymmetric phospholipid membranes, have potential applications in biological and medicinal processes. Here, we used the dissipative particle dynamics simulation method to predict the asymmetric phospholipid microstructures in aqueous solutions. The asymmetric phospholipid membranes, tubes and vesicles are determined and characterized by the chain density distributions and order parameters. The phase diagrams are constructed to evaluate the effects of the chain length on the asymmetric structure formations at equilibrium states, while the average radius of gyration and shape factors are calculated to analyze the asymmetric structure formations in the non-equilibrium processes. Meanwhile, we predicted the mechanical properties of the asymmetric membranes by analyzing the spatial distributions of the interface tensions and osmotic pressures in solutions.

Amphiphilic hyperbranched polymers: synthesis, characterization and self-assembly performance

Abstract

A series of amphiphilic hyperbranched polymers (AHP-s, the “s” refers to the algebra of AHP) were synthesized by the reaction between hydroxyl-terminated hyperbranched polymers (HBP-s, the “s” refers to the algebra of HBP) and palmitoyl chloride. FTIR, NMR and GPC were used to determine the structure of AHP-s, the results showed that AHP-s exhibits core-shell structure. The thermal properties of polymers were investigated by DSC and TGA. It was found that AHP-2, AHP-3 and AHP-4 display higher thermal stability than AHP-1 (AHP-1, AHP-2, AHP-3 and AHP-4 represent the first, second, third and fourth generation AHP, respectively). Furthermore, the self-assembly performance of AHP-s in THF solvent was investigated by TEM and SEM. Finally, the encapsulation capacity of the AHP-s for methyl orange (MO) was explored at different concentrations of AHP-s and pH conditions. It was found that AHP-s is capable of accommodating hydrophilic guest MO. Moreover, the higher generation of AHP-s, the stronger encapsulation capacity obtained. And the encapsulation capacity closely associated with the pH of encapsulation system.

Graphical abstract

Shear-induced microstructures and dynamics processes of phospholipid cylinders in solutions

Shear-induced microstructures and their corresponding dynamic processes are investigated for phospholipid cylinders in aqueous solution by dissipative particle dynamic simulation. Various phospholipid cylinders with cross-sections, which are formed under shear-free flow, are selected to examine the effects of shear flow on their structures and dynamic processes. Shear flow induces the transition from cylinders into vesicles at weak rate and the transition into vesicle–lamella mixtures with increased shear rate and lamella structures at the strong shear rate. Then, the average radius of gyration and shape factors of the polymer chains in the dynamic processes are discussed in detail. Results show that shear flow causes the structure of the polymer chains to be elongated along the shear direction, and the configuration of the polymer chain can be rapidly transformed into an ellipsoid structure under strong shear.

Shape Transformations of Vesicles Self-Assembled from Amphiphilic Hyperbranched Multiarm Copolymers via Simulation

The understanding of shape transformations of vesicles is of fundamental importance in biological and clinical sciences. Hyperbranched polymer vesicles (branched polymersomes) are newly emerging polymer vesicles with special structure and property. They have also been regarded as a good model for biomembranes. However, the shape transformations of hyperbranched polymer vesicles have not been studied from either an experimental or theoretical level. Herein, the shape transformations of vesicles self-assembled from amphiphilic hyperbranched multiarm copolymers (HMCs) in response to the interaction parameters between the hydrophobic core and hydrophilic arms and the polymer concentrations are investigated carefully through dissipative particle dynamics (DPD) simulations. In the morphological phase diagram, two types of vesicles are obtained: one type corresponds to vesicles without holes formed at low concentrations including unilamellar vesicles, double-lamellar vesicles, discocyte-shaped vesicles, and tubular vesicles, and the other type corresponds to vesicles with holes formed at high concentrations including stomatocyte-shaped vesicles, toroidal vesicles, genus-3 (G-3) toroidal vesicles with three holes, and genus-4 (G-4) toroidal vesicles with four holes. In addition, both the self-assembly mechanisms and the dynamics for the formation of these vesicles have been systematically studied. The current work will offer theoretical support for fabricating novel vesicles with various shapes from hyperbranched polymers.

pH-Responsive Zwitterionic Copolymer DHA-PBLG-PCB for Targeted Drug Delivery: A Computer Simulation Study

In this work, the self-assembled behaviors of zwitterionic copolymer docosahexaenoic acid- b-poly(γ-benzyl-l-glutamate)- b-poly(carboxybetaine methacrylate) (DHA-PBLG-PCB) and the loading and release mechanism of the anticancer drug doxorubicin (DOX) was investigated via computer simulations. The effects of polymer concentration, drug content, and pH on polymeric micelles were explored by dissipative particle dynamics (DPD) simulations. Simulation results show that DHA-PBLG₁₅-PCB₁₀ can self-assemble into core-shell micelles; in addition, the drug-loaded micelles have a pH-responsive feature. DOX can be encapsulated into the core-shell micelle under normal physiological pH conditions, whereas it can be released under acidic pH conditions. The self-assembled behaviors of copolymer DHA-PBLG-PEG were also studied to have a comparison with those of DHA-PBLG-PCB. The DHA-PBLG₁₅-PCB₁₀ system has a stable structure and it has a great potential to serve as drug delivery vehicles for targeted drug delivery.

Self-assembly of phospholipid molecules in solutions under shear flows: Microstructures and phase diagrams

Shear-induced microstructures and their phase diagrams were investigated for phospholipid molecules in aqueous solution by dissipative particle dynamic simulation. Self-assembled microstructures, including spherical and cylindrical micelles, spherical vesicles, lamellae, undulated lamellae, perforated lamellae, and continuous networks, were observed under various shear flows and phospholipid concentrations, where the spatial inhomogeneity and symmetry were analysed. A series of phase diagrams were constructed based on the chain lengths under various phospholipid concentrations. The phase distributions showed that the structures with spherical symmetry could be shear-induced to structures with cylindrical symmetry in the dilute solutions. In the semi-concentrated solutions, the lamellae were located in most spaces under zero shear flows, which could be shear-induced into undulated lamellae and then into cylindrical micelles. For the concentrated solutions, the strong shear flows oriented the directions of multilayer lamellae and phase transitions appeared between several cylindrical network structures. These observations on shear-induced microstructures and their distributions revealed a promising approach that could be used to design bio-microstructures based on phospholipid molecules under shear flows.

Activity–crowding coupling effect on the diffusion dynamics of a self-propelled particle in polymer solutions

The anomalous diffusion dynamics of an active particle in polymer solutions is studied based on a Langevin Brownian dynamics simulation.

Self-assembly and stimuli-responsive behaviours of side-chain liquid crystalline copolymers: a dissipative particle dynamics simulation approach

Side-chain liquid crystalline copolymers are able to self-assemble into various aggregates in selective solvents, in particular, deformed structures.

Comparative study of the crowding-induced collapse effect in hard-sphere, flexible polymer and rod-like polymer systems

A systematic Langevin simulation is performed to study the crowding-induced collapse effect on a probed chain in three typical systems: hard sphere (HS), flexible polymer and rod-like polymer.

Relationship between water permeation and flip-flop motion in a bilayer membrane

The lipid bilayer membrane facilitates various biological reactions and is thus an essential structure that sustains all higher forms of life. The unique local environment of the lipid bilayer plays critical roles for the diffusion of biomolecules as well as water molecules in biological reactions. Although fluctuation of the cell membrane is expected to allow for the transport of some water molecules, the flip-flop of lipid molecules corresponds to lipid transport between membrane leaflets, and is considered to be an important process to regulate the lipid composition of biological membranes. However, the relationship between these flip-flop phenomena and surrounding water molecules remains poorly understood. We hypothesized that the flip-flop is caused by water molecules permeating through the cell membrane. To test this hypothesis, we used millisecond-order coarse-grained molecular simulations (dissipative particle dynamics) to investigate the distance between water molecules and lipid molecules depending on the position of the lipid molecule. The results clearly showed that water molecules affect the flip-flop motion in the early stage, but have minimal contribution to the subsequent behavior. Moreover, based on the results of dissipative particle dynamics simulation, we computed several first-passage-time (FPT) quantities to describe the detailed dynamics of water permeation. We modeled arrangements in the middle of the flip-flop process, which were compared with the arrangement without lipid molecules. Overall, our results indicate that lipid molecules located both in perpendicular and parallel arrangements largely affect water permeation. These findings provide new insight into the detailed relationship between water permeation and the flip-flop motion.

Self-Assembly of Lysine-Based Dendritic Surfactants Modeled by the Self-Consistent Field Approach

Implementing a united atom model, we apply self-consistent field theory to study structure and thermodynamic properties of spherical micelles composed of surfactants that combine an alkyl tail with a charged lysine-based dendritic headgroup. Following experiments, the focus was on dendron surfactants with varying tail length and dendron generations G0, G1, G2. The heads are subject to acetylation modification which reduces the charge and hydrophilicity. We establish a reasonable parameter set which results in semiquantitative agreement with the available experiments. The critical micellization concentration, aggregation number, and micelle size are discussed. The strongly charged dendronic surfactants micelles are stable for generation numbers G0 and G1, for progressively higher ionic strengths. Associates of G2 surfactants are very small and can only be found at extreme surfactant concentration and salt strengths. Micelles of corresponding weaker charged acetylated variants exist up to G2, tolerate significantly lower salt concentrations, but lose the spherical micelle topology for G0 at high ionic strengths.

Scission Free Energies for Wormlike Surfactant Micelles: Development of a Simulation Protocol, Application, and Validation for Personal Care Formulations

We present a scheme to calculate wormlike micelle scission free energies from a potential of mean force (PMF) derived from a weighted histogram analysis method (WHAM) applied to coarse grained dissipative particle dynamics (DPD) simulations. In contrast to previous related work, we use a specially chosen external potential based on a reaction coordinate that reversibly drives surfactants out of the nascent scission location. For the application to a model body wash formulation, we predict how addition of NaCl and small molecules such as perfume raw materials (PRMs) affect scission energies. The results show qualitative agreement and correct trends compared to recently determined scission energies for the same system; however, a more rigorous parametrization of the underlying DPD potential is required for quantitative agreement.

One-pot production of porous assemblies by PISA of star architecture copolymers: a simulation study

Porous vesicles can be produced in one-pot by the PISA of star architecture copolymers.