Mesoscale simulations of diffusion and sedimentation in shape-anisotropic nanoparticle suspensions

We determine the long-time self-diffusion coefficient and sedimentation coefficient for suspensions of nanoparticles with anisotropic shapes (octahedra, cubes, tetrahedra, and spherocylinders) as a function of nanoparticle concentration using mesoscale simulations. We use a discrete particle model for the nanoparticles, and we account for solvent-mediated hydrodynamic interactions between nanoparticles using the multiparticle collision dynamics method. Our simulations are compared to theoretical predictions and experimental data from existing literature, demonstrating good agreement in the majority of cases. Further, we find that the self-diffusion coefficient of the regular polyhedral shapes can be estimated from that of a sphere whose diameter is the average of their inscribed and circumscribed sphere diameters.

Sequence-dependent material properties of biomolecular condensates and their relation to dilute phase conformations

Material properties of phase-separated biomolecular condensates, enriched with disordered proteins, dictate many cellular functions. Contrary to the progress made in understanding the sequence-dependent phase separation of proteins, little is known about the sequence determinants of condensate material properties. Using the hydropathy scale and Martini models, we computationally decipher these relationships for charge-rich disordered protein condensates. Our computations yield dynamical, rheological, and interfacial properties of condensates that are quantitatively comparable with experimentally characterized condensates. Interestingly, we find that the material properties of model and natural proteins respond similarly to charge segregation, despite different sequence compositions. Molecular interactions within the condensates closely resemble those within the single-chain ensembles. Consequently, the material properties strongly correlate with molecular contact dynamics and single-chain structural properties. We demonstrate the potential to harness the sequence characteristics of disordered proteins for predicting and engineering the material properties of functional condensates, with insights from the dilute phase properties.

Sequence-Dependent Material Properties of Biomolecular Codensates and their Relation to Dilute Phase Conformations

Material properties of phase-separated biomolecular assemblies, enriched with disordered proteins, dictate their ability to participate in many cellular functions. Despite the significant effort dedicated to understanding how the sequence of the disordered protein drives its phase separation to form condensates, little is known about the sequence determinants of condensate material properties. Here, we computationally decipher these relationships for charged disordered proteins using model sequences comprised of glutamic acid and lysine residues as well as naturally occurring sequences of LAF1's RGG domain and DDX4's N-terminal domain. We do so by delineating how the arrangement of oppositely charged residues within these sequences influences the dynamical, rheological, and interfacial properties of the condensed phase through equilibrium and non-equilibrium molecular simulations using the hydropathy scale and Martini models. Our computations yield material properties that are quantitatively comparable with experimentally characterized condensate systems. Interestingly, we find that the material properties of both the model and natural proteins respond similarly to the segregation of charges, despite their very different sequence compositions. Condensates of the highly charge-segregated sequences exhibit slower dynamics than the uniformly charge-patterned sequences, because of their comparatively long-lived molecular contacts between oppositely charged residues. Surprisingly, the molecular interactions within the condensate are highly similar to those within a single-chain for all sequences. Consequently, the condensate material properties of charged disordered proteins are strongly correlated with their dense phase contact dynamics and their single-chain structural properties. Our findings demonstrate the potential to harness the sequence characteristics of disordered proteins for predicting and engineering the material properties of functional condensates, with insights from the dilute phase properties.

Sequence-Dependent Conformational Transitions of Disordered Proteins During Condensation

Intrinsically disordered proteins (IDPs) can form biomolecular condensates through phase separation. It is recognized that the conformation of IDPs in the dense and dilute phases as well as at the interfaces of condensates can critically impact the resulting properties associated with their functionality. However, a comprehensive understanding of the conformational transitions of IDPs during condensation remains elusive. In this study, we employ a coarse-grained polyampholyte model, comprising an equal number of oppositely charged residues-glutamic acid and lysine-whereby conformations and phase behavior can be readily tuned by altering the protein sequence. By manipulating the sequence patterns from perfectly alternating to block-like, we obtain chains with ideal-like conformations to semi-compact structures in the dilute phase, while in the dense phase, the chain conformation is approximately that of an ideal chain, irrespective of the protein sequence. By performing simulations at different concentrations, we find that the chains assemble from the dilute phase through small oligomeric clusters to the dense phase, accompanied by a gradual swelling of the individual chains. We further demonstrate that these findings are applicable to several naturally occurring proteins involved in the formation of biological condensates. Concurrently, we delve deeper into the chain conformations within the condensate, revealing that chains at the interface show a strong sequence dependence, but remain more collapsed than those in the bulk-like dense phase. This study addresses critical gaps in our knowledge of IDP conformations within condensates as a function of protein sequence.

Structure Formation in Supraparticles Composed of Spherical and Elongated Particles

We studied the evaporation-induced formation of supraparticles from dispersions of elongated colloidal particles using experiments and computer simulations. Aqueous droplets containing a dispersion of ellipsoidal and spherical polystyrene particles were dried on superamphiphobic surfaces at different humidity values that led to varying evaporation rates. Supraparticles made from only ellipsoidal particles showed short-range lateral ordering at the supraparticle surface and random orientations in the interior regardless of the evaporation rate. Particle-based simulations corroborated the experimental observations in the evaporation-limited regime and showed an increase in the local nematic ordering as the diffusion-limited regime was reached. A thin shell of ellipsoids was observed at the surface when supraparticles were made from binary mixtures of ellipsoids and spheres. Image analysis revealed that the supraparticle porosity increased with an increasing aspect ratio of the ellipsoids.

Confined Dynamics in Spherical Polymer Brushes

We investigate the dynamics of polymers grafted to spherical nanoparticles in solution using hybrid molecular dynamics simulations with a coarse-grained solvent modeled via the multiparticle collision dynamics algorithm. The mean-square displacements of monomers near the surface of the nanoparticle exhibit a plateau on intermediate time scales, indicating confined dynamics reminiscent of those reported in neutron spin-echo experiments. The confined dynamics vanish beyond a specific radial distance from the nanoparticle surface that depends on the polymer grafting density. We show that this dynamical confinement transition follows theoretical predictions for the critical distance associated with the structural transition from confined to semidilute brush regimes. These findings suggest the existence of a hitherto unreported dynamic length scale connected with theoretically predicted static fluctuations in spherical polymer brushes and provide new insights into recent experimental observations.

Conformational Properties of Polymers at Droplet Interfaces as Model Systems for Disordered Proteins

Polymer models serve as useful tools for studying the formation and physical properties of biomolecular condensates. In recent years, the interface dividing the dense and dilute phases of condensates has been discovered to be closely related to their functionality, but the conformational preferences of the constituent proteins remain unclear. To elucidate this, we perform molecular simulations of a droplet formed by phase separation of homopolymers as a surrogate model for the prion-like low-complexity domains. By systematically analyzing the polymer conformations at different locations in the droplet, we find that the chains become compact at the droplet interface compared with the droplet interior. Further, segmental analysis revealed that the end sections of the chains are enriched at the interface to maximize conformational entropy and are more expanded than the middle sections of the chains. We find that the majority of chain segments lie tangential to the droplet surface, and only the chain ends tend to align perpendicular to the interface. These trends also hold for the natural proteins FUS LC and LAF-1 RGG, which exhibit more compact chain conformations at the interface compared to the droplet interior. Our findings provide important insights into the interfacial properties of biomolecular condensates and highlight the value of using simple polymer physics models to understand the underlying mechanisms.

Aggregation of amphiphilic nanocubes in equilibrium and under shear

We investigate the self-assembly of amphiphilic nanocubes into finite-sized aggregates in equilibrium and under shear, using molecular dynamics (MD) simulations and kinetic Monte Carlo (KMC) calculations. These patchy nanoparticles combine both interaction and shape anisotropy, making them valuable models for studying folded proteins and DNA-functionalized nanoparticles. The nanocubes can self-assemble into various finite-sized aggregates ranging from rods to self-avoiding random walks, depending on the number and placement of the hydrophobic faces. Our study focuses on suspensions containing multi- and one-patch cubes, with their ratio systematically varied. When the binding energy is comparable to the thermal energy, the aggregates consist of only few cubes that spontaneously associate/dissociate. However, highly stable aggregates emerge when the binding energy exceeds the thermal energy. Generally, the mean aggregation number of the self-assembled clusters increases with the number of hydrophobic faces and decreases with increasing fraction of one-patch cubes. In sheared suspensions, the more frequent collisions between nanocube clusters lead to faster aggregation dynamics but also to smaller terminal steady-state mean cluster sizes. The results from the MD and KMC simulations are in excellent agreement for all investigated two-patch cases, whereas the three-patch cubes form systematically smaller clusters in the MD simulations compared to the KMC calculations due to finite-size effects and slow aggregation kinetics. By analyzing the rate kernels, we are able to identify the primary mechanisms responsible for (shear-induced) cluster growth and breakup. This understanding allows us to tune nanoparticle and process parameters to achieve desired cluster sizes and shapes.

Conformational Properties of Polymers at Droplet Interfaces as Model Systems for Disordered Proteins

Polymer models serve as useful tools for studying the formation and physical properties of biomolecular condensates. In recent years, the interface dividing the dense and dilute phases of condensates has been discovered to be closely related to their functionality, but the conformational preferences of the constituent proteins remain unclear. To elucidate this, we perform molecular simulations of a droplet formed by liquid‒liquid phase separation of homopolymers, as a surrogate model for the prion-like low-complexity domains. By systematically analyzing the polymer conformations at different locations in the droplet, we find that the chains become compact at the droplet interface compared to the droplet interior. Further, segmental analysis revealed that the end sections of the chains are enriched at the interface to maximize conformational entropy, and are more expanded than the middle sections of the chains. We find that the majority of chain segments lie tangential to the droplet surface and only the chain ends tend to align perpendicular to the interface. These trends also hold for the natural proteins FUC LC and LAF-1 RGG, which exhibit more compact chain conformations at the interface compared with the droplet interior. Our findings provide important insights into the interfacial properties of biomolecular condensates and highlight the value of using simple polymer physics models to understand the underlying mechanisms.

Role of Strong Localized vs Weak Distributed Interactions in Disordered Protein Phase Separation

Interaction strength and localization are critical parameters controlling the single-chain and condensed-state properties of intrinsically disordered proteins (IDPs). Here, we decipher these relationships using coarse-grained heteropolymers comprised of hydrophobic (H) and polar (P) monomers as model IDPs. We systematically vary the fraction of P monomers XP and employ two distinct particle-based models that include either strong localized attractions between only H-H pairs (HP model) or weak distributed attractions between both H-H and H-P pairs (HP+ model). To compare different sequences and models, we first carefully tune the attraction strength for all sequences to match the single-chain radius of gyration. Interestingly, we find that this procedure produces similar conformational ensembles, nonbonded potential energies, and chain-level dynamics for single chains of almost all sequences in both models, with some deviations for the HP model at large XP. However, we observe a surprisingly rich phase behavior for the sequences in both models that deviates from the expectation that similarity at the single-chain level will translate to a similar phase-separation propensity. Coexistence between dilute and dense phases is only observed up to a model-dependent XP, despite the presence of favorable interchain interactions, which we quantify using the second virial coefficient. Instead, the limited number of attractive sites (H monomers) leads to the self-assembly of finite-sized clusters of different sizes depending on XP. Our findings strongly suggest that models with distributed interactions favor the formation of liquid-like condensates over a much larger range of sequence compositions compared to models with localized interactions.

Effect of Surfactant Selectivity on Shape and Inner Morphology of Triblock Terpolymer Microparticles

The confined assembly of block copolymers has become a useful tool to prepare microparticles with controlled anisotropy and inner structure. While a solid understanding of the behavior of AB diblock copolymers exists, knowledge on the parameters affecting ABC triblock terpolymer assembly is much more limited. In this work, we analyze the effect of block-selective surfactants, sodium-4-vinylbenzenesulfonate (VBS) and sodium dodecylsulfate (SDS), in the evaporation-induced confined assembly (EICA) of a polystyrene-block-polybutadiene-block-poly(methyl methacrylate) triblock terpolymer (SBM). Despite using the same terpolymer and emulsification process, SDS results in ellipsoidal microparticles with axially stacked lamellae while VBS results in spherical microparticles with concentric lamellae or 3D spiral morphology. This change in morphology upon switching the surfactant was further substantiated by molecular simulations and enhances our understanding of terpolymer microphase separation in confinement. This article is protected by copyright. All rights reserved.

Rapid Precipitation of Ionomers for Stabilization of Polymeric Colloids

Polymeric colloids have shown potential as "building blocks" in applications ranging from formulations of Pickering emulsions and drug delivery systems to advanced materials, including colloidal crystals and composites. However, for applications requiring tunable properties of charged colloids, obstacles in fabrication can arise through limitations in process scalability and chemical versatility. In this work, the capabilities of flash nanoprecipitation (FNP), a scalable nanoparticle (NP) fabrication technology, are expanded to produce charged polystyrene colloids using sulfonated polystyrene ionomers as a new class of NP stabilizers. Through experimental exploration of formulation parameters, increases in the ionomer content are shown to reduce the particle size, mitigating a significant trade-off between the final particle size and inlet concentration; thus, expanding the processable material throughput of FNP. Further, the degree of sulfonation is found to impact stabilization with optimal performance achieved by selecting ionomers with intermediate (2.45-5.2 mol %) sulfonation. Simulations of single ionomer chains and their arrangement in multicomponent NPs provide molecular insights into the assembly and structure of NPs wherein the partitioning of ionomers to the particle surface depends on the polymer molecular weight and degree of sulfonation. By combining the insights from simulations with diffusion-limited growth kinetics and parametric fits to experimental data, a simple design formulation relation is proposed and validated. This work highlights the potential of ionomer-based stabilizers for controllably producing charged NP dispersions in a scalable manner.

Effect of Charge Distribution on the Dynamics of Polyampholytic Disordered Proteins

The stability and physiological function of many biomolecular coacervates depend on the structure and dynamics of intrinsically disordered proteins (IDPs) that typically contain a significant fraction of charged residues. Although the effect of relative arrangement of charged residues on IDP conformation is a well-studied problem, the associated changes in dynamics are far less understood. In this work, we systematically interrogate the effects of charge distribution on the chain-level and segmental dynamics of polyampholytic IDPs in dilute solutions. We study a coarse-grained model polyampholyte consisting of an equal fraction of two oppositely charged residues (glutamic acid and lysine) that undergoes a transition from an ideal chain-like conformation for uniformly charge-patterned sequences to a semi-compact conformation for highly charge-segregated sequences. Changes in the chain-level dynamics with increasing charge segregation correlate with changes in conformation. The chain-level and segmental dynamics conform to simple homopolymer models for uniformly charge-patterned sequences but deviate with increasing charge segregation, both in the presence and absence of hydrodynamic interactions. We discuss the significance of these findings, obtained for a model polyampholyte, in the context of a charge-rich intrinsically disordered region of the naturally occurring protein LAF-1. Our findings have important implications for understanding the effects of charge patterning on the dynamics of polyampholytic IDPs in dilute conditions using polymer scaling theories.

Crystal Prediction via Genetic Algorithms in a Model Chiral System

Chiral crystals and their constituent molecules play a prominent role in theories about the origin of biological homochirality and in drug discovery, design, and stability. Although the prediction and identification of stable chiral crystal structures is crucial for numerous technologies, including separation processes and polymorph selection and control, predictive ability is often complicated by a combination of many-body interactions and molecular complexity and handedness. In this work, we address these challenges by applying genetic algorithms to predict the ground-state crystal lattices formed by a chiral tetramer molecular model, which we have previously shown to exhibit complex fluid-phase behavior. Using this approach, we explore the relative stability and structures of the model's conglomerate and racemic crystals, and present a structural phase diagram for the stable Bravais crystal types in the zero-temperature limit.

Single-chain and condensed-state behavior of hnRNPA1 from molecular simulations

Intrinsically disordered proteins (IDPs) are essential components for the formation of membraneless organelles, which play key functional and regulatory roles within biological systems. These complex assemblies form and dissolve spontaneously over time via liquid-liquid phase separation of IDPs. Mutations in their amino acid sequence can alter their phase behavior, which has been linked to the emergence of severe diseases. We study the conformation and phase behavior of a low-complexity domain of heterogeneous nuclear ribonucleoprotein A1 (hnRNPA1), using coarse-grained implicit solvent molecular simulations. We systematically analyze how these properties are affected by the number of aromatic residues within the examined sequences. We find a significant compaction of the chains and an increase in the critical temperature with increasing number of aromatic residues. The local persistence length is determined in single-chain simulations, revealing strong sequence-specific variations along the chain contour. Comparing single-chain and condensed state simulations, we find much more collapsed polymer conformations in the dilute systems, even at temperatures near the estimated θ-temperature of the solution. These observations strongly support the hypothesis that aromatic residues play a dominant role for condensation, which is further corroborated by a detailed analysis of the intermolecular contacts, and conversely that important properties of condensates are captured in coarse-grained simulations. Interestingly, we observe density inhomogeneities within the condensates near criticality, which are driven by electrostatic interactions. Finally, we find that the relatively small fraction of hydrophobic residues in the IDPs results in interfacial tensions which are significantly lower compared to typical combinations of immiscible simple liquids.

Self-Assembly of Amphiphilic Cubes in Suspension

We study the self-assembly of amphiphilic cubic colloids using molecular dynamics as well as rejection-free kinetic Monte Carlo simulations. We vary both the number and location of the solvophobic faces (patches) on the cubes at several colloid volume fractions and determine the resulting size and shape distributions of the self-assembled aggregates. When the binding energy is comparable to the thermal energy of the system, aggregates typically consist of only few spontaneously associating/dissociating colloids. Increasing the binding energy (or lowering the temperature) leads to the emergence of highly stable aggregates, e.g., small dimers in pure suspensions of one-patch cubes or large (system-spanning) aggregates in suspensions of multipatch colloids. In mixtures of one- and multipatch cubes, the average aggregation number increases with increasing number of solvophobic faces on the multipatch cubes as well with increasing fraction of multipatch cubes. The resulting aggregate shapes range from elongated rods over fractal objects to compact spheres, depending on the number and arrangement of solvophobic patches on the cubic colloids. Our findings establish the complex self-assembly pathways for a class of building blocks that combine both interaction and shape anisotropy, with the potential of forming hierarchically ordered superstructures.

Understanding Gas Transport in Polymer-Grafted Nanoparticle Assemblies

We rationalize the unusual gas transport behavior of polymer-grafted nanoparticle (GNP) membranes. While gas permeabilities depend specifically on the chemistry of the polymers considered, we focus here on permeabilities relative to the corresponding pure polymer which show interesting, "universal" behavior. For a given NP radius, R_c, and for large enough areal grafting densities, σ, to be in the dense brush regime we find that gas permeability enhancements display a maximum as a function of the graft chain molecular weight, M_n. Based on a recently proposed theory for the structure of a spherical brush in a melt of GNPs, we conjecture that this peak permeability occurs when the densely grafted polymer brush has the highest, packing-induced extension free energy per chain. The corresponding brush thickness is predicted to be h max = 3 R c , independent of chain chemistry and σ, i.e., at an apparently universal value of the NP volume fraction (or loading), ϕ_NP, ϕ_NP,max = [R_c/(R_c + h_max)]³ ≈ 0.049. Motivated by this conclusion, we measured CO_-2 and CH₄ permeability enhancements across a variety of R_c, M_n and σ, and find that they behave in a similar manner when considered as a function of ϕ_NP, with a peak in the near vicinity of the predicted ϕ_NP,max. Thus, the chain length dependent extension free energy appears to be the critical variable in determining the gas permeability for these hybrid materials. The emerging picture is that these curved polymer brushes, at high enough σ behave akin to a two-layer transport medium - the region in the near vicinity of the NP surface is comprised of extended polymer chains which speed-up gas transport relative to the unperturbed melt. The chain extension free energy increases with increasing chain length, up to a maximum, and apparently leads to an increasing gas permeability. For long enough grafts, there is an outer region of chain segments that is akin to an unperturbed melt with slow gas transport. The permeability maximum and decreasing permeability with increasing chain length then follow naturally.

Diffusion and sedimentation in colloidal suspensions using multiparticle collision dynamics with a discrete particle model

We study self-diffusion and sedimentation in colloidal suspensions of nearly hard spheres using the multiparticle collision dynamics simulation method for the solvent with a discrete mesh model for the colloidal particles (MD+MPCD). We cover colloid volume fractions from 0.01 to 0.40 and compare the MD+MPCD simulations to experimental data and Brownian dynamics simulations with free-draining hydrodynamics (BD) as well as pairwise far-field hydrodynamics described using the Rotne-Prager-Yamakawa mobility tensor (BD+RPY). The dynamics in MD+MPCD suggest that the colloidal particles are only partially coupled to the solvent at short times. However, the long-time self-diffusion coefficient in MD+MPCD is comparable to that in experiments, and the sedimentation coefficient in MD+MPCD is in good agreement with that in experiments and BD+RPY, suggesting that MD+MPCD gives a reasonable description of hydrodynamic interactions in colloidal suspensions. The discrete-particle MD+MPCD approach is convenient and readily extended to more complex shapes, and we determine the long-time self-diffusion coefficient in suspensions of nearly hard cubes to demonstrate its generality.

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 particles, which represent the gas of interest, segregate preferentially to the grafting surface, with this result being relatively insensitive to the nature of polymer-surface interactions. These brush layers therefore correspond to heterogeneous transport media: the gas molecules near the grafting surface have accelerated dynamics (presumably parallel to the wall) relative to the corresponding polymer melt, but they have slower dynamics in the central region of the brush. We therefore find that gas molecules perform hop-like motions - they spend a significant part of their time in the regions of fast transport, separated by motions where they "hop" from one surface to the other. These phenomena in combination lead to an overall speedup in gas dynamics in these brush layers relative to a polymer melt, in good agreement with the experimental data.

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 occur. It is found that the phase diagrams may exhibit a triple point or a nematic-nematic critical point, and that coexisting phases differ appreciably in their monomer densities. The properties of the two types of chains (nematic order parameters, chain radii, etc.) in the various phases are studied in detail, and predictions on the (anisotropic) critical behavior near the critical point of nematic-nematic unmixing are made.

Ultra-coarse-graining of homopolymers in inhomogeneous systems

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 entropy minimization. While the CG models with three-body interactions fail at reproducing stable polymer films and droplets, CG simulations with an LDP are able to do so. Minor quantitative differences between the reference and the CG simulations, namely a slight broadening of interfaces accompanied by a smaller surface tension in the CG simulations, can be attributed to the deformation of polymers near the interfaces, which cannot be resolved in the CG representation, where the polymers are mapped to spherical beads.

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-graining. Providing the NNs with additional information on the forces greatly improves the accuracy of the predictions since more correlations are taken into account; the predicted potentials become smoother, are significantly closer to the target potentials, and are more transferable as a result.

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 the difference in bending stiffness is increased, the two chain types separate into two coexisting nematic phases. The phase behavior is rationalized by considering the chemical potentials of the two components and the Gibbs excess free energy. The geometric properties of the chain conformations under the various conditions are also discussed.

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 our finding that GNPs are spatially heterogeneous transport media, with this free volume distribution being easily manipulated by the addition of free polymer. In particular, adding a small amount of appropriately chosen free polymer can increase the membrane gas selectivity by up to two orders of magnitude while only moderately reducing small gas permeability. Added short free chains, which are homogeneously distributed in the polymer layer of the GNP, reduce the permeability of all gases but yield no dramatic increases in selectivity. In contrast, free chains with length comparable to the grafts, which populate the interstitial pockets between GNPs, preferentially hinder the transport of the larger gas and thus result in large selectivity increases. This work thus establishes that we can favorably manipulate the selective gas transport properties of GNP membranes through the entropic effects associated with the addition of free chains.

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.

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 aggregates. At low shear rates, all mixtures exhibit a Newtonian-like response with intrinsic shear viscosities that are up to 2 times as large as of pure suspensions of nonadsorbing colloids at the same concentration. Furthermore, the mean aggregation number increases slightly due to the flow-enhanced collision of aggregates. At larger shear rates, however, the aggregates break up, the polymers align in the flow direction, and the mixtures exhibit shear-thinning. This shear-induced breakup occurs at stronger shear compared to pure JC suspensions, indicating that the adsorbed polymers reinforce the hybrid aggregates.

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 NP by spherical dry and interpenetration layers. This model quantitatively predicts that the thicknesses of the two layers depend on one universal parameter, x, the degree of overcrowding of grafted chains relative to chains in the melt. Both simulations and theory show that the chain extension free energy is nonmonotonic with increasing chain length at a fixed grafting density, with a well-defined maximum. This maximum is indicative of the crossover from the dry layer-dominated to interpenetration layer-dominated regime, and it could have profound consequences on our understanding of a variety of anomalous transport properties of these GNPs. Our theoretical approach therefore provides a facile means for understanding and designing solvent-free GNP-based materials.

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 the treatment of the solvent in the two models. We rationalize the presence (or absence) of stratification by measuring phenomenological multicomponent diffusion coefficients for the polymer mixtures. The diffusion coefficients show the importance of not only solvent backflow but also HI between polymers in controlling the dried supraparticle morphology.

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 hydrophobic homopolymers and the addition of an amphiphilic block copolymer. The final structure of the colloids depends on the mass fraction of each homopolymer as well as the concentration and composition of the block copolymer additive. To confirm the surface activity of the colloids, we demonstrate their potential to stabilize Pickering emulsions. This hybrid approach of simulations and experiments provides a pathway to designing and manufacturing complex polymeric colloids on an industrial scale.

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 exhibit a similar self-assembly behavior. At small and intermediate shear, we observe deformation and rearrangement of the micelles, accompanied by a Newtonian-like rheology with slightly higher shear viscosity compared to homoparticle dispersions at the same concentration. As the shear rate is increased further, the micelles eventually break up into small dimers and free particles, causing a distinct shear-thinning of the dispersions. The network-like structures exhibit a similar flow behavior at high shear rates, but for weak shear we find an almost threefold increase of the shear viscosity and a distinct shear-thinning behavior due to the fracturing of the intertwined networks. In general, we identify a strong correlation between the size of the aggregates and the rheology of the dispersions, allowing for the determination of dynamic properties solely based on structural information.

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-liquid interface, the solvophobic particles deform into oblate ellipsoids to maximize the occluded area at the interface. In contrast, amphiphilic Janus colloids orient their solvophobic/solvophilic parts toward the preferred liquids and form a prolate particle shape. The amphiphilic core-shell particles undergo a morphological transition to a prolate Janus-like structure as they anchor at the interface. We rationalize the deformation of the polymer colloids by considering a simple model system of spheroidal particles pinned at the liquid-liquid interface. We systematically compute the interfacial free energy for the various colloids as a function of their asphericity and find excellent qualitative agreement with the simulation findings. Our results show that solvophobic homoparticles can be almost as efficient surface-active agents as amphiphilic Janus colloids.

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

The reverse perturbation method [Phys. Rev. E 59, 4894 (1999)1063-651X10.1103/PhysRevE.59.4894] for shearing simple liquids and measuring their viscosity is extended to the Vicsek model (VM) of active particles [Phys. Rev. Lett. 75, 1226 (1995)PRLTAO0031-900710.1103/PhysRevLett.75.1226] 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 ν and the momentum amplification coefficient λ, can be obtained by fitting this decay with an analytical solution of the hydrodynamic equations for the VM. The viscosity of the VM consists of two parts, a kinetic and a collisional contribution. While analytical predictions already exist for the former, a novel expression for the collisional part is derived by an Enskog-like kinetic theory. To verify the predictions for the transport coefficients, Green-Kubo relations were evaluated and transverse current correlations were measured in independent simulations. Not too far to the transition to collective motion, we find excellent agreement between the different measurements of the transport coefficients. However, the measured values of ν and 1-λ are always slightly higher than the mean-field predictions, even at large mean free paths and at state points quite far from the threshold to collective motion, that is, far in the disordered phase. These findings seem to indicate that the mean-field assumption of molecular chaos is much less reliable in systems with velocity-alignment rules such as the VM, compared to models obeying detailed balance such as multiparticle collision dynamics.

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 "sandwiched" between A and C lamellae. Despite thorough synthetic variation of the B fraction to control the geometry of the sandwiched microphase, so far only Janus spheres, cylinders, and sheets have been obtained. In this combined theoretical and experimental work, we show that the blending of polybutadiene homopolymer (hPB) into lamella morphologies of polystyrene-block-polybutadiene-block-polymethylmethacrylate (SBM) triblock terpolymers allows the continuous tuning of the polybutadiene (PB) microphase. We systematically vary the volume fraction of hPB in the system, and we find in both experiments and simulations morphological transitions from PB-cylinders to perforated PB-lamellae and further to continuous PB-lamellae. Our simulations show that the hPB is distributed homogeneously in the PB microdomains. Through crosslinking of the PB domain and redispersion in a common solvent for all blocks, we separate the bulk morphologies into Janus cylinders, perforated Janus sheets, and Janus sheets. These studies suggest that more complex Janus nanostructures could be generated from ABC triblock terpolymers than previously expected.

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. Toward the center, the number of large colloids increases and they are packed amorphously. The extent of this stratification decreases with decreasing the evaporation rate. Complementary simulations indicate that evaporation leads to a local increase in density, which, in turn, exerts stronger inward forces on the larger colloids. A comparison between experiments and simulations suggest that hydrodynamic interactions between the suspended colloids reduce the extent of stratification. Our findings are relevant for the fabrication of supraparticles for applications in the fields of chromatography, catalysis, drug delivery, photonics, and a better understanding of spray-drying.

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 primarily depended on the balance between inertial lift forces pushing the droplet outward and elastic forces from the polymers driving it inward. The droplet shape was controlled by the local shear rate, and so its average shape depended on the droplet distribution.

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

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 enhance the physical properties of composite materials by the strategic design of the "molecular architecture" of constituents to benefit from synergistic effects relating to the organization of the polymer component.

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 controlled by the polymer feed concentration, the NP feed concentration, and their ratio in the feed streams. Selective localization of NPs by incorporating electrostatic interactions between polymers and NPs has also been investigated, and was shown to be an effective way to enhance NP loading and surface density as compared to the case with only van der Waals attractions. This work demonstrates that the FNP process is promising for the production of structured and hybrid nanocolloids in a continuous and scalable way, with independent control over particle properties such as size, NP location, loading, and surface density. Our results provide useful guidelines for experimental fabrication of such hybrid nanoparticles.

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 lp increases) indicate that the NP motions become decoupled from the bulk viscosity of the polymer solution. Likewise, polymer stiffness leads to deviations from PCT, which was developed for fully flexible chains. Independent of lp, however, the long-time diffusion behavior is well-described by MCT, particularly at high polymer concentration. We also observed that the short-time subdiffusive dynamics are strongly dependent on polymer flexibility. As lp is increased, the NP dynamics become more subdiffusive and decouple from the dynamics of the polymer chain center-of-mass. We posit that these effects are due to differences in the segmental mobility of the semiflexible chains.

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 we replaced water with a nonpolar poor solvent. The polymers aggregated into stable nanoparticles for a range of processing parameters. We investigated theoretically three possible explanations for this stability, i.e., electrostatic stabilization, conditional thermodynamic equilibrium, and steric stabilization. Our experiments and considerations suggest that steric stabilization is the most likely candidate.

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 director is well described by a Gaussian, irrespective of whether the chains are rod-like or rather flexible. However, the related concept of "deflection length" is shown to make sense only in the latter case for rather dilute solutions since otherwise the deflection length is of the order of about two bond lengths only. When the solution is semi-dilute, a substantial renormalization of the persistence length occurs, while this effect is absent in the isotropic phase even at rather high monomer densities. The effective radii of the "tubes" confining the chains in the related description of orientational ordering are significantly larger than the distances between neighboring chains, providing evidence for a pronounced collective character of orientational fluctuations. Hairpins can be identified close to the isotropic-nematic transition, and their probability of occurrence agrees qualitatively with the Vroege-Odijk theory. The corresponding theoretical predictions for the elastic constants, however, are not in good agreement with the simulations. We attribute the shortcomings of the theories to their neglect of the coupling between local density and orientational fluctuations. Finally, we detected for this model a transition to a smectic phase for reduced monomer densities near 0.7.

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 elucidate the nature of orientational ordering in different parameter regimes. Previous simulations have revealed the emergence of bipolar and quadrupolar topological defects on the surface when the entire cavity was filled with a concentrated polymer solution [Phys. Rev. Lett., 2017, 118, 217803]. In contrast, spherical shell confinement restricts the appearance of a bipolar order. Instead, only the extent of the quadrupolar order changes with chain stiffness, as evidenced by the relative motion of topological defects. In the case of monolayers, we observe a nematic to smectic transition accompanied by a change in the nematic grain-size distribution as the contour length was decreased.