From Molecular Constraints to Macroscopic Dynamics in Associative Networks Formed by Ionizable Polymers: A Neutron Spin Echo and Molecular Dynamics Simulations Study.

The association of ionizable polymers strongly affects their motion in solutions, where the constraints arising from clustering of the ionizable groups alter the macroscopic dynamics. The interrelation between the motion on multiple length and time scales is fundamental to a broad range of complex fluids including physical networks, gels, and polymer-nanoparticle complexes where long-lived associations control their structure and dynamics. Using neutron spin echo and fully atomistic, multimillion atom molecular dynamics (MD) simulations carried out to times comparable to that of chain segmental motion, the current study resolves the dynamics of networks formed by suflonated polystryene solutions for sulfonation fractions 0 ≤ f ≤ 0.09 across time and length scales. The experimental dynamic structure factors were measured and compared with computational ones, calculated from MD simulations, and analyzed in terms of a sum of two exponential functions, providing two distinctive time scales. These time constants capture confined motion of the network and fast dynamics of the highly solvated segments. A unique relationship between the polymer dynamics and the size and distribution of the ionic clusters was established and correlated with the number of polymer chains that participate in each cluster. The correlation of dynamics in associative complex fluids across time and length scales, enabled by combining the understanding attained from reciprocal space through neutron spin echo and real space, through large scale MD studies, addresses a fundamental long-standing challenge that underline the behavior of soft materials and affect their potential uses.

From ionic clusters dynamics to network constraints in ionic polymer solutions.

Physical networks formed by ionizable polymers with ionic clusters as crosslinks are controlled by coupled dynamics that transcend from ionic clusters through chain motion to macroscopic response. Here, the coupled dynamics, across length scales, from the ionic clusters to the networks in toluene swollen polystyrene sulfonate networks, were directly correlated, as the electrostatic environment of the physical crosslinks was altered. The multiscale insight is attained by coupling neutron spin echo measurements with molecular dynamics simulations, carried out to times typical of relaxation of polymers in solutions. The experimental dynamic structure factor is in outstanding agreement with the one calculated from computer simulations, as the networks are perturbed by elevating the temperature and changing the electrostatic environment. In toluene, the long-lived clusters remain stable over hundreds of ns across a broad temperature range, while the polymer network remains dynamic. Though the size of the clusters changes as the dielectric constant of the solvent is modified through the addition of ethanol, they remain stable but morph, enhancing the polymer chain dynamics.

Inducing stratification of colloidal mixtures with a mixed binary solvent.

Molecular dynamics simulations are used to demonstrate that a binary solvent can be used to stratify colloidal mixtures when the suspension is rapidly dried. The solvent consists of two components, one more volatile than the other. When evaporated at high rates, the more volatile component becomes depleted near the evaporation front and develops a negative concentration gradient from the bulk of the mixture to the liquid-vapor interface while the less volatile solvent is enriched in the same region and exhibit a positive concentration gradient. Such gradients can be used to drive a binary mixture of colloidal particles to stratify if one is preferentially attracted to the more volatile solvent and the other to the less volatile solvent. During solvent evaporation, the fraction of colloidal particles preferentially attracted to the less volatile solvent is enhanced at the evaporation front, whereas the colloidal particles having stronger attractions with the more volatile solvent are driven away from the interfacial region. As a result, the colloidal particles show a stratified distribution after drying, even if the two colloids have the same size.

Response of ionizable block copolymer assemblies to solvent dielectrics: A molecular dynamics study.

Ionizable copolymers assembly in solutions is driven by the formation of ionic clusters. Fast clustering of the ionizable blocks often leads to the formation of far-from equilibrium assemblies that ultimately impact the structure of polymer membranes and affect their many applications. Using large-scale atomistic molecular dynamics simulations, we probe the effects of electrostatics on the formation of ionizable copolymer micelles that dominate their solution structure, with the overarching goal of defining the factors that control the assembly of structured ionizable copolymers. A symmetric pentablock ionizable copolymer, with a randomly sulfonated polystyrene center tethered to polyethylene-r-propylene block, terminated by poly(t-butyl styrene), in solvents of varying dielectric constants from 2 to 20, serves as the model system. We find that independent of the solvents, this polymer forms a core-shell micelle with the ionizable segment segregating to the center of the assembly. The specific block conformation, however, strongly depends on the sulfonation levels and the dielectric constant and the polarity of the solvents.

Fractal dimensions of jammed packings with power-law particle size distributions in two and three dimensions.

Static structure factors are computed for large-scale, mechanically stable, jammed packings of frictionless spheres (three dimensions) and disks (two dimensions) with broad, power-law size dispersity characterized by the exponent -ß. The static structure factor exhibits diverging power-law behavior for small wave numbers, allowing us to identify a structural fractal dimension d_{f}. In three dimensions, d_{f}≈2.0 for 2.5≤ß≤3.8, such that each of the structure factors can be collapsed onto a universal curve. In two dimensions, we instead find 1.0≲d_{f}≲1.34 for 2.1≤ß≤2.9. Furthermore, we show that the fractal behavior persists when rattler particles are removed, indicating that the long-wavelength structural properties of the packings are controlled by the large particle backbone conferring mechanical rigidity to the system. A numerical scheme for computing structure factors for triclinic unit cells is presented and employed to analyze the jammed packings.

Unexpected Slow Relaxation Dynamics in Pure Ring Polymers Arise from Intermolecular Interactions.

Ring polymers have fascinated scientists for decades, but experimental progress has been challenging due to the presence of linear chain contaminants that fundamentally alter dynamics. In this work, we report the unexpected slow stress relaxation behavior of concentrated ring polymers that arises due to ring-ring interactions and ring packing structure. Topologically pure, high molecular weight ring polymers are prepared without linear chain contaminants using cyclic poly(phthalaldehyde) (cPPA), a metastable polymer chemistry that rapidly depolymerizes from free ends at ambient temperatures. Linear viscoelastic measurements of highly concentrated cPPA show slow, non-power-law stress relaxation dynamics despite the lack of linear chain contaminants. Experiments are complemented by molecular dynamics (MD) simulations of unprecedentedly high molecular weight rings, which clearly show non-power-law stress relaxation in good agreement with experiments. MD simulations reveal substantial ring-ring interpenetrations upon increasing ring molecular weight or local backbone stiffness, despite the global collapsed nature of single ring conformation. A recently proposed microscopic theory for unconcatenated rings provides a qualitative physical mechanism associated with the emergence of strong inter-ring caging which slows down center-of-mass diffusion and long wavelength intramolecular relaxation modes originating from ring-ring interpenetrations, governed by the onset variable N/ND, where the crossover degree of polymerization ND is qualitatively predicted by theory. Our work overcomes challenges in achieving ring polymer purity and by characterizing dynamics for high molecular weight ring polymers. Overall, these results provide a new understanding of ring polymer physics.

Polydots, soft nanoparticles, at membrane interfaces.

Soft nanoparticles (NPs) are emerging candidates for nano medicine, particularly for intercellular imaging and targeted drug delivery. Their soft nature, manifested in their dynamics, allows translocation into organisms without damaging their membranes. A crucial step towards incorporating soft dynamic NPs in nano medicine, is to resolve their interrelation with membranes. Here using atomistic molecular dynamics (MD) simulations we probe the interaction of soft NPs formed by conjugated polymers with a model membrane. These NPs, often termed polydots, are confined to their nano dimensions without any chemical tethers, forming dynamic long lived nano structures. Specifically, polydots formed by dialkyl para poly phenylene ethylene (PPE), with a varying number of carboxylate groups tethered to the alkyl chains to tune the interfacial charge of the surface of the NP are investigated at the interface with a model membrane that consists of di-palmitoyl phosphatidylcholine (DPPC). We find that even though polydots are controlled only by physical forces, they retain their NP configuration as they transcend the membrane. Regardless of their size, neutral polydots spontaneously penetrate the membrane whereas carboxylated polydots must be driven in, with a force that depends on the charge at their interface, all without significant disruption to the membrane. These fundamental results provide a means to control the position of the nanoparticles with respect to the membrane interfaces, which is key to their therapeutic use.

Entropic Mixing of Ring/Linear Polymer Blends.

The topological constraints of nonconcatenated ring polymers force them to form compact loopy globular conformations with much lower entropy than unconstrained ideal rings. The closed-loop structure of ring polymers also enables them to be threaded by linear polymers in ring/linear blends, resulting in less compact ring conformations with higher entropy. This conformational entropy increase promotes mixing rings with linear polymers. Here, using molecular dynamics simulations for bead-spring chains, ring/linear blends are shown to be significantly more miscible than linear/linear blends and that there is an entropic mixing, negative χ, for ring/linear blends compared to linear/linear and ring/ring blends. In analogy with small angle neutron scattering, the static structure function S(q) is measured, and the resulting data are fit to the random phase approximation model to determine χ. In the limit that the two components are the same, χ = 0 for the linear/linear and ring/ring blends as expected, while χ < 0 for the ring/linear blends. With increasing chain stiffness, χ for the ring/linear blends becomes more negative, varying reciprocally with the number of monomers between entanglements. Ring/linear blends are also shown to be more miscible than either ring/ring or linear/linear blends and stay in single phase for a wider range of increasing repulsion between the two components.

Molecular dynamics simulations of binary sphere mixtures.

Explicit simulations of fluid mixtures of highly size-dispersed particles are constrained by numerical challenges associated with identifying pair-interaction neighbors. Recent algorithmic developments have ameliorated these difficulties to an extent, permitting more efficient simulations of systems with many large and small particles of disperse sizes. We leverage these capabilities to perform molecular dynamics simulations of binary sphere mixtures with elastically stiff particles approaching the hard sphere limit and particle size ratios of up to 50, approaching the colloidal limit. The systems considered consist of 500 large particles and up to nearly 3.6×10^{6} small particles with total particle volume fractions up to 0.51. Our simulations confirm qualitative predictions for correlations between large particles previously obtained analytically and for simulations employing effective depletion interactions, but also reveal additional insights into the near-contact structure that result from the explicit treatment of the small particle solvent. No spontaneous crystal nucleation was observed during the simulations, suggesting that nucleation rates in the fluid-solid coexistence region are too small to observe crystal nucleation for feasible simulation system sizes and timescales.

Modeling Solution Drying by Moving a Liquid-Vapor Interface: Method and Applications.

A method of simulating the drying process of a soft matter solution with an implicit solvent model by moving the liquid-vapor interface is applied to various solution films and droplets. For a solution of a polymer and nanoparticles, we observe "polymer-on-top" stratification, similar to that found previously with an explicit solvent model. Furthermore, "polymer-on-top" is found even when the nanoparticle size is smaller than the radius of gyration of the polymer chains. For a suspension droplet of a bidisperse mixture of nanoparticles, we show that core-shell clusters of nanoparticles can be obtained via the "small-on-outside" stratification mechanism at fast evaporation rates. "Large-on-outside" stratification and uniform particle distribution are also observed when the evaporation rate is reduced. Polymeric particles with various morphologies, including Janus spheres, core-shell particles, and patchy particles, are produced from drying droplets of polymer solutions by combining fast evaporation with a controlled interaction between the polymers and the liquid-vapor interface. Our results validate the applicability of the moving interface method to a wide range of drying systems. The limitations of the method are pointed out and cautions are provided to potential practitioners on cases where the method might fail.

Large-scale frictionless jamming with power-law particle size distributions.

Due to significant computational expense, discrete element method simulations of jammed packings of size-dispersed spheres with size ratios greater than 1:10 have remained elusive, limiting the correspondence between simulations and real-world granular materials with large size dispersity. Invoking a recently developed neighbor binning algorithm, we generate mechanically stable jammed packings of frictionless spheres with power-law size distributions containing up to nearly 4 000 000 particles with size ratios up to 1:100. By systematically varying the width and exponent of the underlying power laws, we analyze the role of particle size distributions on the structure of jammed packings. The densest packings are obtained for size distributions that balance the relative abundance of large-large and small-small particle contacts. Although the proportion of rattler particles and mean coordination number strongly depend on the size distribution, the mean coordination of nonrattler particles attains the frictionless isostatic value of six in all cases. The size distribution of nonrattler particles that participate in the load-bearing network exhibits no dependence on the width of the total particle size distribution beyond a critical particle size for low-magnitude exponent power laws. This signifies that only particles with sizes greater than the critical particle size contribute to the mechanical stability. However, for high-magnitude exponent power laws, all particle sizes participate in the mechanical stability of the packing.

Superstretchable Elastomer from Cross-linked Ring Polymers.

The stretchability of polymeric materials is critical to many applications such as flexible electronics and soft robotics, yet the stretchability of conventional cross-linked linear polymers is limited by the entanglements between polymer chains. We show using molecular dynamics simulations that cross-linked ring polymers are significantly more stretchable than cross-linked linear polymers. Compared to linear polymers, the entanglements between ring polymers do not act as effective cross-links. As a result, the stretchability of cross-linked ring polymers is determined by the maximum extension of polymer strands between cross-links, rather than between trapped entanglements as in cross-linked linear polymers. The more compact conformation of ring polymers before deformation also contributes to the increase in stretchability.

Overlap Concentration of Sodium Polystyrene Sulfonate in Solution.

The overlap concentration c* of sodium polystyrene sulfonate in water is calculated using multichain atomistic and coarse grained (CG) simulations for a range of chain lengths. Fully atomistic molecular dynamics simulations are carried out for N = 32-192 monomers. The CG model was parameterized to match the end-to-end distance from the atomistic simulations at small N and allows us to simulate a much larger N. Treating the hydrophobic backbone by inclusion of attraction between monomers is an essential addition to the CG model. The simulation c* are in agreement with experimental data, yet at c*, the chains are not fully stretched, even for N as large as 1200. This implies that none of the experimental systems are in the scaling regime and to reach the scaling regime for NaPSS chains much longer than N = 1200 are required.

Flow and arrest in stressed granular materials.

Flowing granular materials often abruptly arrest if not driven by sufficient applied stresses. Such abrupt cessation of motion can be economically expensive in industrial materials handling and processing, and is significantly consequential in intermittent geophysical phenomena such as landslides and earthquakes. Using discrete element simulations, we calculate states of steady flow and arrest for granular materials under the conditions of constant applied pressure and shear stress, which are also most relevant in practice. Here the material can dilate or compact, and flow or arrest, in response to the applied stress. Our simulations highlight that under external stress, the intrinsic response of granular materials is characterized by uniquely-defined steady states of flow or arrest, which are highly sensitive to interparticle friction. While the flowing states can be equivalently characterized by volume fraction, coordination number or internal stress ratio, to characterize the states of shear arrest, one needs to also consider the structural anisotropy in the contact network. We highlight the role of dilation in the flow-arrest transition, and discuss our findings in the context of rheological transitions in granular materials.

Effects of interaction strength of associating groups on linear and star polymer dynamics.

A small number of associating groups incorporated onto a polymer backbone have dramatic effects on the mobility and viscoelastic response of the macromolecules in melts. These associating groups assemble, driving the formation of clusters, whose lifetime affects the properties of the polymers. Here, we probe the effects of the interaction strength on the structure and dynamics of two topologies, linear and star polymer melts, and further investigate blends of associative and non-associating polymers using molecular dynamics simulations. Polymer chains of approximately one entanglement length are described by a bead-spring model, and the associating groups are incorporated in the form of interacting beads with an interaction strength between them that is varied from 1 to 20 kBT. We find that, for all melts and blends, interaction of a few kBT between the associating groups drives cluster formation, where the size of the clusters increases with increasing interaction strength. These clusters act as physical crosslinkers, which slow the chain mobility. Blends of chains with and without associating groups macroscopically phase separate for interaction strength between the associating groups of a few kBT and above. For weakly interacting associating groups, the static structure function S(q) is well fit by functional form predicted by the random phase approximation where a clear deviation occurs as phase segregation takes place, providing a quantitative assessment of phase segregation.

Diffusion of Thin Nanorods in Polymer Melts.

The diffusion of monomerically thin nanorods in polymer melts is studied by molecular dynamics simulations. We focus on the systems where chains are long enough to screen the hydrodynamic interactions such that the diffusion coefficient D ∥ for the direction parallel to the rod decreases linearly with increasing rod length l. In unentangled polymers, the diffusion coefficient for the direction normal to the rod exhibits a crossover from D ⊥ ~ l -2 to D ⊥ ~ l -1 with increasing l, corresponding to a progressive coupling of nanorod motion to the polymers. Accordingly, the rotational diffusion coefficient D R ≈ D ⊥ l -2 ~ l -4 and then D R ~ l -3 as l increases. In entangled polymers, D ⊥ and D R are suppressed for l larger than the entanglement mesh size a. D ⊥ ~ l -3 and D R ~ l -5 for l sufficiently above a in agreement with de Gennes' rod reptation model.

Shear Is Not Always Simple: Rate-Dependent Effects of Flow Type on Granular Rheology.

Despite there being an infinite variety of types of flow, most rheological studies focus on a single type such as simple shear. Using discrete element simulations, we explore bulk granular systems in a wide range of flow types at large strains and characterize invariants of the stress tensor for different inertial numbers and interparticle friction coefficients. We identify a strong dependence on the type of flow, which grows with increasing inertial number or friction. Standard models of yielding, repurposed to describe the dependence of the stress on flow type in steady-state flow and at finite rates, are compared with data.

Assembly of Polymer-Grafted Nanoparticles in Polymer Matrices.

Leibler pioneered the idea that long enough matrix polymers of length P will spontaneously dewet a chemically identical polymer layer, comprising chains of length N, densely end-grafted to a flat surface ("brush"). This entropically driven idea is routinely used to explain experiments in which 10-20 nm diameter nanoparticles (NPs) densely grafted with polymer chains are found to phase separate from chemically identical melts for P/N ≳4. At lower grafting densities, these effects are also thought to underpin the self-assembly of grafted NPs into a variety of structures. To explore the validity of this picture, we conducted large-scale molecular dynamics simulations of grafted NPs in a chemically identical polymer melt. For the NPs we consider, in the ≈5 nm diameter range, we find no phase separation even for P/N = 10 in the absence of attractions. Instead, we find behavior that more closely parallels experiments when all of the chain monomers are equally attractive to each other but repel the NPs. Our results thus imply that experimental situations investigated to date are dominated by the surfactancy of the NPs, which is driven by the chemical mismatch between the inorganic core and the organic ligands (the graft and free chains are chemically identical). Entropic effects, that is, the translational entropy of the NPs and the matrix, the entropy of mixing of the grafts and the matrix, and the conformational entropy of the chains appear to thus have a second-order effect even in the context of these model systems.

Granular packings with sliding, rolling, and twisting friction.

Intuition tells us that a rolling or spinning sphere will eventually stop due to the presence of friction and other dissipative interactions. The resistance to rolling and spinning or twisting torque that stops a sphere also changes the microstructure of a granular packing of frictional spheres by increasing the number of constraints on the degrees of freedom of motion. We perform discrete element modeling simulations to construct sphere packings implementing a range of frictional constraints under a pressure-controlled protocol. Mechanically stable packings are achievable at volume fractions and average coordination numbers as low as 0.53 and 2.5, respectively, when the particles experience high resistance to sliding, rolling, and twisting. Only when the particle model includes rolling and twisting friction were experimental volume fractions reproduced.

Topological Linking Drives Anomalous Thickening of Ring Polymers in Weak Extensional Flows.

Molecular dynamics simulations confirm recent extensional flow experiments showing ring polymer melts exhibit strong extension-rate thickening of the viscosity at Weissenberg numbers Wi≪1. Thickening coincides with the extreme elongation of a minority population of rings that grows with Wi. The large susceptibility of some rings to extend is due to a flow-driven formation of topological links that connect multiple rings into supramolecular chains. Links form spontaneously with a longer delay at lower Wi and are pulled tight and stabilized by the flow. Once linked, these composite objects experience larger drag forces than individual rings, driving their strong elongation. The fraction of linked rings depends nonmonotonically on Wi, increasing to a maximum when Wi∼1 before rapidly decreasing when the strain rate approaches 1/τ_{e}.