Conformations and orientational ordering of semiflexible polymers in spherical confinement.

Semiflexible polymers in lyotropic solution confined inside spherical nanoscopic “containers” with repulsive walls are studied by molecular dynamics simulations and density functional theory, as a first step to model confinement effects on stiff polymers inside of miniemulsions, vesicles, and cells. It is shown that the depletion effects caused by the monomer-wall repulsion depend distinctly on the radius R of the sphere. Further, nontrivial orientational effects occur when R, the persistence length ℓp, and the contour length L of the polymers are of similar magnitude. At intermediate densities, a “shell” of wall-attached chains is forming, such that the monomers belonging to those chains a…

Necessary Catheter Diameters for Mechanical Thrombectomy with ADAPT

BACKGROUND AND PURPOSE: Large-bore catheters allow mechanical thrombectomy in ischemic stroke by engaging and retrieving clots without additional devices (direct aspiration first-pass technique [ADAPT]). The purpose of this study was to establish a model for minimal catheter diameters needed for ADAPT. MATERIALS AND METHODS: We established a theoretic model for the calculation of minimal catheter diameters needed for ADAPT. We then verified its validity in 28 ADAPT maneuvers in a porcine in vivo model. To account for different mechanical thrombectomy techniques, we factored in ADAPT with/without a hypothetic 0.021-inch microcatheter or 0.014-inch microwire inside the lumen of the aspiration…

Ordering, phase behavior, and correlations of semiflexible polymers in confinement.

Semiflexible polymers are ubiquitous in biological systems, e.g., as building blocks of the cytoskeleton, and they also play an important role in various materials due to their ability to form liquid-crystalline order. These rigid macromolecules are characterized by numerous (hierarchical) length-scales that define their static and dynamic properties. Confinement can promote uniform order, e.g., through capillary nematization in narrow slits, but it can also introduce long-ranged disruptions of the nematic ordering field through (unavoidable) topological defects in spherical containers. This Perspective concentrates on the theoretical description and computational modeling of such confined …

Understanding the properties of liquid-crystalline polymers by computational modeling

Abstract A topical review of recent theoretical work on the properties of lyotropic solutions and melts containing semiflexible polymers in thermal equilibrium is given, with a focus on the liquid-crystalline and smectic order of these systems in the bulk and under confinement. Starting with a discussion of single chain properties in terms of the Kratky-Porod worm-like chain model and its limitations, extensions along the lines of Onsager’s theory for the isotropic-nematic transition of solutions of hard rods are briefly reviewed. This discussion is followed by a review of recent Molecular Dynamics simulations and classical Density Functional Theory calculations. It is argued that, even in …

Directed Assembly of Soft Colloids through Rapid Solvent Exchange

We studied the directed assembly of soft nanoparticles through rapid micromixing of polymers in solution with a nonsolvent. Both experiments and computer simulations were performed to elucidate the underlying physics and to investigate the role of various process parameters. In particular, we discovered that no external stabilizing agents or charged end groups are required to keep the colloids separated from each other when water is used as the nonsolvent. Furthermore, the size of the nanoparticles can be reliably tuned through the mixing rate and the ratio between polymer solution and nonsolvent. Our results demonstrate that this mechanism is highly promising for the mass fabrication of un…

Optimizing endovascular stroke treatment: removing the microcatheter before clot retrieval with stent-retrievers increases aspiration flow

BackgroundFlow control during endovascular stroke treatment with stent-retrievers is crucial for successful revascularization. The standard technique recommended by stent-retriever manufacturers implies obstruction of the respective access catheter by the microcatheter, through which the stent-retriever is delivered. This, in turn, results in reduced aspiration during thrombectomy. In order to maximize aspiration, we fully retract the microcatheter out of the access catheter before thrombectomy—an approach we term the ‘bare wire thrombectomy’ (BWT) technique. We verified the improved throughput with systematic in vitro studies and assessed the clinical effectiveness and safety of this metho…

Janus Nanostructures from ABC/B Triblock Terpolymer Blends

Lamella-forming ABC triblock terpolymers are convenient building blocks for the synthesis of soft Janus nanoparticles (JNPs) by crosslinking the B domain that is &ldquo

The smectic phase in semiflexible polymer materials: A large scale Molecular Dynamics study

Abstract Semiflexible polymers in concentrated lyotropic solution are studied within a bead-spring model by molecular dynamics simulations, focusing on the emergence of a smectic A phase and its properties. We systematically vary the density of the monomeric units for several contour lengths that are taken smaller than the chain persistence length. The difficulties concerning the equilibration of such systems and the choice of appropriate ensemble (constant volume versus constant pressure, where all three linear dimensions of the simulation box can fluctuate independently) are carefully discussed. Using HOOMD-blue on graphics processing units, systems containing more than a million monomeri…

Ultra-coarse-graining of homopolymers in inhomogeneous systems

Abstract We develop coarse-grained (CG) models for simulating homopolymers in inhomogeneous systems, focusing on polymer films and droplets. If the CG polymers interact solely through two-body potentials, then the films and droplets either dissolve or collapse into small aggregates, depending on whether the effective polymer–polymer interactions have been determined from reference simulations in the bulk or at infinite dilution. To address this shortcoming, we include higher order interactions either through an additional three-body potential or a local density-dependent potential (LDP). We parameterize the two- and three-body potentials via force matching, and the LDP through relative entr…

Influence of polymer flexibility on nanoparticle dynamics in semidilute solutions

The hierarchical structure and dynamics of polymer solutions control the transport of nanoparticles (NPs) through them. Here, we perform multi-particle collision dynamics simulations of solutions of semiflexible polymer chains with tunable persistence length lp to investigate the effect of chain stiffness on NP transport. The NPs exhibit two distinct dynamical regimes - subdiffusion on short time scales and diffusion on long time scales. The long-time NP diffusivities are compared with predictions from the Stokes-Einstein relation (SER), mode-coupling theory (MCT), and a recent polymer coupling theory (PCT). Increasing deviations from the SER as the polymer chains become more rigid (i.e. as…

Rapid Production of Internally Structured Colloids by Flash Nanoprecipitation of Block Copolymer Blends.

Colloids with internally structured geometries have shown great promise in applications ranging from biosensors to optics to drug delivery, where the internal particle structure is paramount to performance. The growing demand for such nanomaterials necessitates the development of a scalable processing platform for their production. Flash nanoprecipitation (FNP), a rapid and inherently scalable colloid precipitation technology, is used to prepare internally structured colloids from blends of block copolymers and homopolymers. As revealed by a combination of experiments and simulations, colloids prepared from different molecular weight diblock copolymers adopt either an ordered lamellar morph…

Entropic Unmixing in Nematic Blends of Semiflexible Polymers

Binary mixtures of semiflexible polymers with the same chain length but different persistence lengths separate into two coexisting different nematic phases when the osmotic pressure of the lyotropic solution is varied. Molecular Dynamics simulations and Density Functional Theory predict phase diagrams either with a triple point, where the isotropic phase coexists with two nematic phases, or a critical point of unmixing within the nematic mixture. The difference in locally preferred bond angles between the constituents drives this unmixing without any attractive interactions between monomers.

Phase Separation and Nematic Order in Lyotropic Solutions: Two Types of Polymers with Different Stiffnesses in a Common Solvent

The interplay of the isotropic-nematic transition and phase separation in lyotropic solutions of two types of semiflexible macromolecules with pronounced difference in chain stiffness is studied by Density Functional Theory and Molecular Dynamics simulations. While the width of the isotropic-nematic two-phase coexistence region is narrow for solutions with a single type of semiflexible chain, the two-phase coexistence region widens for solutions containing two types of chains with rather disparate stiffness. In the nematic phase, both types of chains contribute to the nematic order, with intermediate values of the order parameter compared to the corresponding single component solutions. As …

Structure and Shear Response of Janus Colloid–Polymer Mixtures in Solution

We investigate the structure and rheological properties of dilute colloid-polymer mixtures at rest and under shear via molecular simulations that take into account hydrodynamic interactions. Mixtures of amphiphilic Janus colloids (JCs) and hydrophobic/amphiphilic polymers are considered for various solvent qualities and polymer concentrations. Free polymers, small polymer droplets, and hybrid aggregates coexist in mixtures with slightly hydrophobic homopolymers. As the solvent quality worsens, all polymers aggregate into small droplets, covered and stabilized by the JCs. In mixtures with amphiphilic polymers, we observe the coexistence of free polymers, purely polymeric micelles, and hybrid…

Self-assembly of colloidal micelles in microfluidic channels.

The self-assembly of amphiphilic Janus colloids in microfluidic channels is studied using hybrid molecular dynamics simulations with fully resolved hydrodynamic interactions incorporated through the multi-particle collision dynamics algorithm. The simulations are conducted at a density and temperature where the Janus particles spontaneously self-assemble into spherical micelles to minimize the interface between the solvophobic caps and the surrounding solvent. In confined systems, this contact area can also be reduced by aggregation at the channel walls. Indeed, a sizable fraction of free particles and small clusters with three and four members are found at the walls when the microfluidic c…

Gas Transport in Interacting Planar Brushes

Recent experiments on melts of spherical nanoparticles (NPs) densely grafted with polymer chains show enhanced gas transport relative to the neat polymer (without NPs). As a means of understanding this unexpected behavior, we consider here the simpler case of two interacting planar brushes, under conditions representing a polymer melt far below its critical point (i.e., where the "free volume" or holes act akin to a poor solvent). Computer simulations illustrate, in agreement with mean-field ideas, that the density profile far away from the walls is flat but with a value that is marginally larger than the corresponding polymer melt under identical state conditions. We find that tracer parti…

Semiflexible Polymers in Spherical Confinement: Bipolar Orientational Order Versus Tennis Ball States

Densely packed semiflexible polymers with contour length L confined in spheres with radius R of the same order as L cannot exhibit uniform nematic order. Depending on the chain stiffness (which we vary over a wide range), highly distorted structures form with topological defects on the sphere surface. These structures are completely different from previously observed ones of very long chains winding around the inner surface of spheres and from nematic droplets. At high densities, a thin shell of polymers close to the sphere surface exhibits a tennis ball texture due to the confinement-induced gradual bending of polymer bonds. In contrast, when the contour length of the chains is significant…

In Silico Design Enables the Rapid Production of Surface-Active Colloidal Amphiphiles

A new technology platform built on the integration of theory and experiments to enable the design of Janus colloids with precision control of surface anisotropy and amphiphilicity could lead to a disruptive transformation in the next generation of surfactants, photonic or phononic materials, and coatings. Here, we exploit molecular dynamics (MD) simulations to guide the rational design of amphiphilic polymer Janus colloids by Flash NanoPrecipitation (FNP), a method capable of the production of colloids with complex structure without the compromise of reduced scalability. Aided by in silico design, we show in experiments that amphiphilic Janus colloids can be produced using a unique blend of…

On the Stability of Polymeric Nanoparticles Fabricated through Rapid Solvent Mixing.

We study the stability of polymeric nanoparticles fabricated through the rapid mixing of polymers in a good solvent with a poor solvent that is miscible with the good solvent. In previous experiments where water was used as the poor solvent, a negative surface charge was measured on the precipitated nanoparticles, which led to the long-time stability of the dispersion. It was argued that these charges originate presumably from either water or hydroxide adsorption at the hydrophobic nanoparticle surface or from impurities in the feed streams that preferentially adsorb on the precipitated nanoparticles. To elucidate the origin of this stabilization mechanism, we performed experiments wherein …

Structure of Polymer-Grafted Nanoparticle Melts.

The structure of neat melts of polymer-grafted nanoparticles (GNPs) is studied via coarse-grained molecular dynamics simulations. We systematically vary the degree of polymerization and grafting density at fixed nanoparticle (NP) radius and study in detail the shape and size of the GNP coronas. For sufficiently high grafting density, chain sections close to the NP core are extended and form a dry layer. Further away from the NP, there is an interpenetration layer, where the polymer coronas of neighboring GNPs overlap and the chain sections have almost unperturbed conformations. To better understand this partitioning, we develop a two-layer model, representing the grafted polymer around an N…

Self-assembly of semiflexible polymers confined to thin spherical shells

Confinement effects are critical for stiff macromolecules in biological cells, vesicles, and other systems in soft matter. For these molecules, the competition between the packing entropy and the enthalpic cost of bending is further shaped by strong confinement effects. Through coarse-grained molecular dynamics simulations, we explore the self-assembly of semiflexible polymers confined in thin spherical shells for various chain lengths, chain stiffnesses, and shell thicknesses. Here, we focus on the case where the contour and persistence length of the polymers are comparable to the radius of the confining cavity. The range of ordered structures is analyzed using several order parameters to …

Blends of Semiflexible Polymers: Interplay of Nematic Order and Phase Separation

Mixtures of semiflexible polymers with a mismatch in either their persistence lengths or their contour lengths are studied by Density Functional Theory and Molecular Dynamics simulation. Considering lyotropic solutions under good solvent conditions, the mole fraction and pressure is systematically varied for several cases of bending stiffness κ (the normalized persistence length) and chain length N. For binary mixtures with different chain length (i.e., NA=16, NB=32 or 64) but the same stiffness, isotropic-nematic phase coexistence is studied. For mixtures with the same chain length (N=32) and large stiffness disparity (κB/κA=4.9 to 8), both isotropic-nematic and nematic-nematic unmixing oc…

Cell-sized confinements alter molecular diffusion in concentrated polymer solutions due to length-dependent wetting of polymers

Living cells are characterized by the micrometric confinement of various macromolecules at high concentrations. Using droplets containing binary polymer blends as artificial cells, we previously showed that cell-sized confinement causes phase separation of the binary polymer solutions because of the length-dependent wetting of the polymers. Here we demonstrate that the wetting-induced heterogeneity of polymers also emerges in single-component polymer solutions. The resulting heterogeneity leads to a slower transport of small molecules at the center of cell-sized droplets than that in bulk solutions. This heterogeneous distribution is observed when longer polymers with lower wettability are …

Surface Activity of Soft Polymer Colloids

We investigate the behavior of polymer colloids at the interface between two immiscible liquids using molecular dynamics simulations. We study several colloid morphologies with various degrees of amphiphilicity, that is, purely solvophobic homogeneous and Janus particles and amphiphilic Janus and core-shell particles. Regardless of the specific morphology, the polymer colloids irreversibly anchor at the liquid-liquid interface, accompanied by a marked reduction of the interfacial tension, γ. Purely solvophobic particles lower γ because they reduce the interfacial area shared by the two immiscible liquids, whereas amphiphilic colloids have an additional enthalpic contribution. At the liquid-…

Coupling of Nanoparticle Dynamics to Polymer Center-of-Mass Motion in Semidilute Polymer Solutions

We investigate the dynamics of nanoparticles in semidilute polymer solutions when the nanoparticles are comparably sized to the polymer coils using explicit- and implicit-solvent simulation methods. The nanoparticle dynamics are subdiffusive on short time scales before transitioning to diffusive motion on long time scales. The long-time diffusivities scale according to theoretical predictions based on full dynamic coupling to the polymer segmental relaxations. In agreement with our recent experiments, however, we observe that the nanoparticle subdiffusive exponents are significantly larger than predicted by the coupling theory over a broad range of polymer concentrations. We attribute this …

Modeling hydrodynamic interactions in soft materials with multiparticle collision dynamics

Multiparticle collision dynamics (MPCD) is a flexible and robust mesoscale computational technique for simulating solvent-mediated hydrodynamic interactions in soft materials. Here, we provide a critical overview of the MPCD method and summarize its current strengths and limitations. The capabilities of the method are highlighted by reviewing its recent applications to simulate diverse phenomena, ranging from the flow of complex fluids and thermo-osmotic transport to bacterial swimming and active particle self-assembly. We also discuss outstanding challenges and emerging methodological developments that are expected to greatly expand the applicability of MPCD to other systems of technologic…

Nanoparticle dynamics in semidilute polymer solutions: Rings versus linear chains

We study the dynamics of nanoparticles in semidilute solutions of ring and linear polymers using hybrid molecular dynamics–multiparticle collision dynamics simulations. The dynamics of the monomers, the polymer centers-of-mass, and the nanoparticles coincide for these two architectures for solutions of the same monomer concentration. The long time diffusivities of the nanoparticles follow the predictions of a polymer coupling theory [Cai et al., Macromolecules 44, 7853–7863 (2011)], suggesting that nanoparticle dynamics are coupled to segmental relaxations for both polymer architectures examined here. At intermediate time scales, the nanoparticle dynamics are characterized by subdiffusive e…

Cross-stream migration of a Brownian droplet in a polymer solution under Poiseuille flow

The migration of a Brownian fluid droplet in a parallel-plate microchannel was investigated using dissipative particle dynamics computer simulations. In a Newtonian solvent, the droplet migrated toward the channel walls due to inertial effects at the studied flow conditions, in agreement with theoretical predictions and recent simulations. However, the droplet focused onto the channel centerline when polymer chains were added to the solvent. Focusing was typically enhanced for longer polymers and higher polymer concentrations with a nontrivial flow-rate dependence due to droplet and polymer deformability. Brownian motion caused the droplet position to fluctuate with a distribution that prim…

Disentangling the Role of Chain Conformation on the Mechanics of Polymer Tethered Particle Materials

[Image: see text] The linear elastic properties of isotropic materials of polymer tethered nanoparticles (NPs) are evaluated using noncontact Brillouin light spectroscopy. While the mechanical properties of dense brush materials follow predicted trends with NP composition, a surprising increase in elastic moduli is observed in the case of sparsely grafted particle systems at approximately equal NP filling ratio. Complementary molecular dynamics simulations reveal that the stiffening is caused by the coil-like conformations of the grafted chains, which lead to stronger polymer–polymer interactions compared to densely grafted NPs with short chains. Our results point to novel opportunities to …

Stratification of polymer mixtures in drying droplets: Hydrodynamics and diffusion

We study the evaporation-induced stratification of a mixture of short and long polymer chains in a drying droplet using molecular simulations. We systematically investigate the effects of hydrodynamic interactions (HI) on this process by comparing hybrid simulations accounting for HI between polymers through the multiparticle collision dynamics technique with free-draining Langevin dynamics simulations neglecting the same. We find that the dried supraparticle morphologies are homogeneous when HI are included but are stratified in core--shell structures (with the short polymers forming the shell) when HI are neglected. The simulation methodology unambiguously attributes this difference to th…

Spatial Demixing of Ring and Chain Polymers in Pressure-Driven Flow

We investigate mixtures of ring and linear polymers in solution at various number ratios, ranging from pure chains to pure rings, and at densities around the overlap concentration. In bulk and at r...

Segregation in Drying Binary Colloidal Droplets

When a colloidal suspension droplet evaporates from a solid surface, it leaves a characteristic deposit in the contact region. These deposits are common and important for many applications in printing, coating, or washing. By the use of superamphiphobic surfaces as a substrate, the contact area can be reduced so that evaporation is almost radially symmetric. While drying, the droplets maintain a nearly perfect spherical shape. Here, we exploit this phenomenon to fabricate supraparticles from bidisperse colloidal aqueous suspensions. The supraparticles have a core-shell morphology. The outer region is predominantly occupied by small colloids, forming a close-packed crystalline structure. Tow…

Structure and dynamics of amphiphilic Janus spheres and spherocylinders under shear.

We study the structure formation and flow properties of colloidal dispersions comprised of Janus spheres, Janus spherocylinders, and their mixtures, using hybrid molecular dynamics simulations that take into account hydrodynamic interactions. We systematically vary the Janus balance and the shape anisotropy of the particles, and explore a range of colloid volume fractions in the liquid regime of the phase diagram. At rest, Janus spheres with small hydrophobic patches form spherical micelles for all investigated colloid concentrations. In contrast, Janus spheres with an entirely hydrophobic hemisphere aggregate to larger worm-like micelles and network-like structures. Janus spherocylinders e…

BoltzmaNN: Predicting effective pair potentials and equations of state using neural networks

Neural networks (NNs) are employed to predict equations of state from a given isotropic pair potential using the virial expansion of the pressure. The NNs are trained with data from molecular dynamics simulations of monoatomic gases and liquids, sampled in the NVT ensemble at various densities. We find that the NNs provide much more accurate results compared to the analytic low-density limit estimate of the second virial coefficient and the Carnahan-Starling equation of state for hard sphere liquids. Furthermore, we design and train NNs for computing (effective) pair potentials from radial pair distribution functions, g(r), a task that is often performed for inverse design and coarse-graini…

Stratification Dynamics in Drying Colloidal Mixtures

Stratification in binary colloidal mixtures was investigated using implicit-solvent molecular dynamics simulations. For large particle size ratios and film Péclet numbers greater than unity, smaller colloids migrated to the top of the film, while big colloids were pushed to the bottom, creating an "inverted" stratification. This peculiar behavior was observed in recent simulations and experiments conducted by Fortini et al. [ Phys. Rev. Lett. 2016 , 116 , 118301 ]. To rationalize this behavior, particle size ratios and drying rates spanning qualitatively different Péclet number regimes were systematically studied, and the dynamics of the inverted stratification were quantified in detail. Th…

Evaporation-induced assembly of colloidal crystals

Colloidal crystals are often prepared by evaporation from solution, and there is considerable interest to link the processing conditions to the crystal morphology and quality. Here, we study the evaporation-induced assembly of colloidal crystals using massive-scale nonequilibrium molecular dynamics simulations. We apply a recently developed machine-learning technique to characterize the assembling crystal structures with unprecedented microscopic detail. In agreement with previous experiments and simulations, faster evaporation rates lead to earlier onset of crystallization and more disordered surface structures. Surprisingly, we find that collective rearrangements of the bulk crystal durin…

Transport coefficients of self-propelled particles: Reverse perturbations and transverse current correlations

The reverse perturbation method [Phys. Rev. E 59, 4894 (1999)] for shearing simple liquids and measuring their viscosity is extended to the Vicsek model (VM) of active particles [Phys. Rev. Lett. 75, 1226 (1995)] and its metric-free version. The sheared systems exhibit a phenomenon that is similar to the skin effect of an alternating electric current: Momentum that is fed into the boundaries of a layer decays mostly exponentially toward the center of the layer. It is shown how two transport coefficients, i.e., the shear viscosity $\ensuremath{\nu}$ and the momentum amplification coefficient $\ensuremath{\lambda}$, can be obtained by fitting this decay with an analytical solution of the hydr…

Equilibrium Dynamics and Shear Rheology of Semiflexible Polymers in Solution

We study the structure and dynamics of semidilute solutions of semiflexible polymers at rest and under shear using hybrid molecular dynamics simulations that take hydrodynamic interactions into account. We show that the polymer center-of-mass diffusion coefficient significantly decreases with increasing chain stiffness at fixed monomer density. The zero-shear viscosity shows a corresponding increase due to the intermolecular interactions of stiffer chains. We apply steady shear flow to the polymer solutions and show that at high shear rates the flow properties become almost independent of polymer stiffness. We characterize the polymer conformations under shear and find that in this regime p…

Topology-Sensitive Microfluidic Filter for Polymers of Varying Stiffness

The separation of polymers based on their size, rigidity, and topology is an essential but also highly challenging task for nanoscience and engineering. Using hybrid molecular dynamics simulations that correctly take into account hydrodynamics, we have designed microfluidic channels for separating linear from ring polymers in dilute solutions. We establish that the transport velocity of the polymers is independent of their topology and rigidity when the channel walls are smooth and repulsive. However, when the walls are decorated with attractive spots arranged on lines parallel to the flow, ring polymers exhibit an order of magnitude higher transport velocity compared to linear chains. The …

On the applicability of density dependent effective interactions in cluster-forming systems

We systematically studied the validity and transferability of effective, coarse-grained, pair potentials in ultrasoft colloidal systems. We focused on amphiphilic dendrimers, macromolecules which can aggregate into clusters of overlapping particles to minimize the contact area with the surrounding (implicit) solvent. Simulations were performed for both the monomeric and coarse-grained model in the liquid phase at densities ranging from infinite dilution up to values close to the freezing point. For every state point, each macromolecule was mapped onto a single interaction site and the effective pair potential was computed using a coarse-graining technique based on force-matching. We found e…

Densely packed semiflexible macromolecules in a rigid spherical capsule

The ordering of semiflexible polymers with persistence length lp and contour length L confined in a sphere of radius R is studied by molecular dynamics simulations of a coarse-grained model. Monomer densities are chosen where the corresponding bulk lyotropic solution or melt is a well-ordered nematic, and purely repulsive walls of the rigid confining sphere are considered. It is found that polymers close to the walls are bent according to the curvature of the confining spheres with all their monomers in a few layers parallel to the sphere surface, whereas the remaining macromolecules closer to the sphere center have one chain end and their center of mass far from the surface. The latter cha…

Dynamics of single semiflexible polymers in dilute solution

We study the dynamics of a single semiflexible chain in solution using computer simulations, where we systematically investigate the effect of excluded volume, chain stiffness, and hydrodynamic interactions. We achieve excellent agreement with previous theoretical considerations, but find that the crossover from the time τb, up to which free ballistic motion of the monomers describes the chain dynamics, to the times W−1 or τ0, where anomalous monomer diffusion described by Rouse-type and Zimm-type models sets in, requires two decades of time. While in the limit of fully flexible chains the visibility of the anomalous diffusion behavior is thus rather restricted, the t3/4 power law predicted…

Controlling Janus Nanodisc Topology through ABC Triblock Terpolymer/Homopolymer Blending in 3D Confinement

Janus particles have drawn considerable interest as colloidal surfactants, microswimmers, and building blocks for colloidal lattices. So far, research primarily focused on spherical Janus particles for which a number of fabrication methods are well established. Janus particles with geometric anisotropy offer shape-dependent properties in addition to surface anisotropy, but their synthesis is more challenging. Here, we report a variety of polymeric Janus nanoparticles synthesized from ABC triblock terpolymer microphases in microemulsion droplets. Evaporation-induced assembly of the ABC triblock terpolymers led to prolate microparticles with A/C lamellae stacked along the particle’s major axi…

Adsorption and structure formation of semiflexible polymers on spherical surfaces

Abstract Rigid spheres with a short-range attractive potential are taken as a coarse-grained model of vesicles, which contain a solution of semiflexible polymers in their interior. Assuming good solvent conditions with an implicit description of the solvent, effective monomers experience bond-length and bond-angle potentials as well as excluded-volume interaction. Due to the attractive vesicle surface, phase separation occurs between a thin shell of adsorbed monomers at the surface and a rather dilute, and therefore, disordered polymer solution in the sphere interior. While at a planar attractive surface the wormlike chains would exhibit liquid crystalline (nematic and smectic) order, the c…

Entropic Unmixing in Nematic Blends of Semiflexible Polymers.

Binary mixtures of semiflexible polymers with the same chain length, but different persistence lengths, separate into two coexisting different nematic phases when the osmotic pressure of the lyotropic solution is varied. Molecular Dynamics simulations and Density Functional Theory predict phase diagrams either with a triple point, where the isotropic phase coexists with two nematic phases or a critical point of unmixing within the nematic mixture. The difference in locally preferred bond angles between the constituents drives this unmixing without any attractive interactions between monomers.

Axial dispersion of Brownian colloids in microfluidic channels

Controlled production of patchy particles from the combined effects of nanoprecipitation and vitrification

Using molecular dynamics simulations, we study a simple and scalable method for fabricating patchy nanoparticles via the assembly of binary polymer blends under a rapid solvent exchange. Patchiness can be achieved by incorporating a glassy component, which kinetically traps the particle morphology along the path to the equilibrium configuration. Our simulations reveal that the number of surface patches increases for larger nanoparticles and for more asymmetric blend ratios, while the size distribution of the patches remains rather uniform. Other than multi-patch nanoparticles, Janus structures have been obtained for small nanoparticles. Further, ribbon structures with elongated surface doma…

Structured Nanoparticles from the Self-Assembly of Polymer Blends through Rapid Solvent Exchange

Molecular dynamics simulations were performed to study systematically the rapid mixing of a polymer blend in solution with a miscible nonsolvent. In agreement with experiments, we observe that polymers self-assemble into complex nanoparticles, such as Janus and core-shell particles, when the good solvent is displaced by the poor solvent. The emerging structures can be predicted on the basis of the surface tensions between the polymers as well as between the polymers and the surrounding liquid. Furthermore, the size of the nanoparticles can be independently tuned through the mixing rate and the polymer concentration in the feed stream; meanwhile, the composition of the nanoparticles can be c…

Tuning Selectivities in Gas Separation Membranes Based on Polymer-Grafted Nanoparticles

Polymer membranes are critical to many sustainability applications that require the size-based separation of gas mixtures. Despite their ubiquity, there is a continuing need to selectively affect the transport of different mixture components while enhancing mechanical strength and hindering aging. Polymer-grafted nanoparticles (GNPs) have recently been explored in the context of gas separations. Membranes made from pure GNPs have higher gas permeability and lower selectivity relative to the neat polymer because they have increased mean free volume. Going beyond this ability to manipulate the mean free volume by grafting chains to a nanoparticle, the conceptual advance of the present work is…

Efficient mesoscale hydrodynamics: Multiparticle collision dynamics with massively parallel GPU acceleration

Abstract We present an efficient open-source implementation of the multiparticle collision dynamics (MPCD) algorithm that scales to run on hundreds of graphics processing units (GPUs). We especially focus on optimizations for modern GPU architectures and communication patterns between multiple GPUs. We show that a mixed-precision computing model can improve performance compared to a fully double-precision model while still providing good numerical accuracy. We report weak and strong scaling benchmarks of a reference MPCD solvent and a benchmark of a polymer solution with research-relevant interactions and system size. Our MPCD software enables simulations of mesoscale hydrodynamics at lengt…

Multi-scale simulations of polymeric nanoparticle aggregation during rapid solvent exchange.

Using a multi-scale approach which combines both molecular dynamics (MD) and kinetic Monte Carlo (KMC) simulations, we study a simple and scalable method for fabricating charge-stabilized nanoparticles through a rapid solvent exchange, i.e., Flash NanoPrecipitation (FNP). This multi-scale approach is based on microscopic information from MD simulations and uses a KMC algorithm to access macroscopic length- and time scales, which allows direct comparison with experiments and quantitative predictions. We find good agreement of our simulation results with the experiments. In addition, the model allows us to understand the aggregation mechanism on both microscopic and macroscopic levels and det…

Erratum: “Evaporation-induced assembly of colloidal crystals” [J. Chem. Phys. 149, 094901 (2018)]

Phase behavior of flexible and semiflexible polymers in solvents of varying quality.

The interplay of nematic order and phase separation in solutions of semiflexible polymers in solvents of variable quality is investigated by density functional theory (DFT) and molecular dynamics (MD) simulations. We studied coarse-grained models, with a bond-angle potential to control chain stiffness, for chain lengths comparable to the persistence length of the chains. We varied both the density of the monomeric units and the effective temperature that controls the quality of the implicit solvent. For very stiff chains, only a single transition from an isotropic fluid to a nematic is found, with a phase diagram of "swan-neck" topology. For less stiff chains, however, also unmixing between…

Coil-Globule Collapse of Polystyrene Chains in Tetrahydrofuran-Water Mixtures.

We study the coil and globule states of a single polymer chain in solution by performing molecular dynamics simulations with a united atom model. Specifically, we characterize the structural properties of atactic polystyrene chains with N = 20–150 monomers in tetrahydrofuran–water mixtures at varying mixing ratios. We find that the hydrophobic polymers form rather open coils when the mole fraction of water, XW, is roughly below 0.25, whereas the chains collapse into globules when XW ≳ 0.75. We confirm the theoretically expected scaling laws for the radius of gyration, Rg, in these regimes, i.e., Rg ∝ N3/5 and Rg ∝ N1/3 for good and poor solvent conditions, respectively. For poor solvent con…

Nematic order in solutions of semiflexible polymers: Hairpins, elastic constants, and the nematic-smectic transition

Coarse-grained models of lyotropic solutions of semiflexible polymers are studied by both molecular dynamics simulations and density functional theory calculations, using an implicit solvent bead-spring model with a bond-angle potential. We systematically vary the monomer density, persistence length, and contour length over a wide range and explore the full range from the isotropic-nematic transition to the nematic-smectic transition. In the nematic regime, we span the entire regime from rigid-rod like polymers to thin wormlike chains, confined in effective straight tubes caused by the collective nematic effective ordering field. We show that the distribution of bond angles relative to the …

Stratification in Drying Polymer–Polymer and Colloid–Polymer Mixtures

Drying polymer-polymer and colloid-polymer mixtures were studied using Langevin dynamics computer simulations. Polymer-polymer mixtures vertically stratified into layers, with the shorter polymers enriched near the drying interface and the longer polymers pushed down toward the substrate. Colloid-polymer mixtures stratified into a polymer-on-top structure when the polymer radius of gyration was comparable to or smaller than the colloid diameter, and a colloid-on-top structure otherwise. We also developed a theoretical model for the drying mixtures based on dynamical density functional theory, which gave excellent quantitative agreement with the simulations for the polymer-polymer mixtures a…

Employing artificial neural networks to find reaction coordinates and pathways for self-assembly

Capturing the autonomous self-assembly of molecular building blocks in computer simulations is a persistent challenge, requiring to model complex interactions and to access long time scales. Advanced sampling methods allow to bridge these time scales but typically require to construct accurate low-dimensional representations of the transition pathways. In this work, we demonstrate for the self-assembly of two single-stranded DNA fragments into a ring-like structure how autoencoder architectures based on unsupervised neural networks can be employed to reliably expose transition pathways and to provide a suitable low-dimensional representation. The assembly occurs as a two-step process throug…

Self-Assembly of Polymer Blends and Nanoparticles through Rapid Solvent Exchange.

Molecular dynamics simulations were performed to study the fabrication of polymeric colloids containing inorganic nanoparticles (NPs) via the flash nanoprecipitation (FNP) technique. During this process, a binary polymer blend, initially in a good solvent for the polymers, is rapidly mixed with NPs and a poor solvent for the polymers that is miscible with the good solvent. The simulations reveal that the polymers formed Janus particles with NPs distributed either on the surface of the aggregates, throughout their interior, or aligned at the interface between the two polymer domains, depending on the NP-polymer and NP-solvent interactions. The loading and surface density of NPs can be contro…