Coexistence of Crumpling and Flat Sheet Conformations in Two-Dimensional Polymer Networks: An Understanding of Aggrecan Self-Assembly.

We investigate the conformational properties of self-avoiding two-dimensional (2D) ideal polymer networks with tunable mesh sizes as a model of self-assembled structures formed by aggrecan. Polymer networks having few branching points and large enough mesh tend to crumple, resulting in a fractal dimension of d_{f}≈2.7. The flat sheet behavior (d_{f}=2) emerges in 2D polymer networks having more branching points at large length scales; however, it coexists with crumpling conformations at intermediate length scales, a feature found in scattering profiles of aggrecan solutions. Our findings bridge the long-standing gap between theories and simulations of polymer sheets.

Bottlebrush polymers in the melt and polyelectrolytes in solution share common structural features.

Uncharged bottlebrush polymer melts and highly charged polyelectrolytes in solution exhibit correlation peaks in scattering measurements and simulations. Given the striking superficial similarities of these scattering features, there may be a deeper structural interrelationship in these chemically different classes of materials. Correspondingly, we constructed a library of isotopically labeled bottlebrush molecules and measured the bottlebrush correlation peak position [Formula: see text] by neutron scattering and in simulations. We find that the correlation length scales with the backbone concentration, [Formula: see text], in striking accord with the scaling of ξ with polymer concentration cP in semidilute polyelectrolyte solutions [Formula: see text] The bottlebrush correlation peak broadens with decreasing grafting density, similar to increasing salt concentration in polyelectrolyte solutions. ξ also scales with sidechain length to a power in the range of 0.35-0.44, suggesting that the sidechains are relatively collapsed in comparison to the bristlelike configurations often imagined for bottlebrush polymers.

Molecular dynamics study of the swelling and osmotic properties of compact nanogel particles.

Owing to their great importance in materials science and other fields, we investigate the solution and osmotic properties of uncharged compact nanogel particles over a wide range of solvent quality and particle concentration by molecular dynamics (MD) simulations. We characterize the osmotic pressure by estimating the second and third virial coefficients, and by extension, we identify the θ-point where the second virial coefficient vanishes. Calculations of the structure factor indicate that these particles are similar to macrogels in that the particle-like scattering profile disappears at moderate concentrations. We also find that improving the solvent quality enhances the spatial segmental uniformity, while significant heterogeneous structure arises near the θ-point. Well below the θ-point where the second osmotic virial coefficient vanishes, these heterogeneous structures become less prevalent as the particles tend to collapse. We also investigate the degree of swelling and structure of compact nanogel particles with a variable excluded volume interaction and gel particle concentration. The osmotic modulus and the scaling exponents in good and θ-point conditions of these gels are characteristic of interacting randomly branched polymers, i.e., "lattice animals".

Influence of network defects on the conformational structure of nanogel particles: From "closed compact" to "open fractal" nanogel particles.

We propose an approach to generate a wide range of randomly branched polymeric structures to gain general insights into how polymer topology encodes a configurational structure in solution. Nanogel particles can take forms ranging from relatively symmetric sponge-like compact structures to relatively anisotropic open fractal structures observed in some nanogel clusters and in some self-associating polymers in solutions, such as aggrecan solutions under physiologically relevant conditions. We hypothesize that this broad "spectrum" of branched polymer structures derives from the degree of regularity of bonding in the network defining these structures. Accordingly, we systematically introduce bonding defects in an initially perfect network having a lattice structure in three and two topological dimensions corresponding to "sponge" and "sheet" structures, respectively. The introduction of bonding defects causes these "closed" and relatively compact nanogel particles to transform near a well-defined bond percolation threshold into "open" fractal objects with the inherent anisotropy of randomly branched polymers. Moreover, with increasing network decimation, the network structure of these polymers acquires other configurational properties similar to those of randomly branched polymers. In particular, the mass scaling of the radius of gyration and its eigenvalues, as well as hydrodynamic radius, intrinsic viscosity, and form factor for scattering, all undergo abrupt changes that accompany these topological transitions. Our findings support the idea that randomly branched polymers can be considered to be equivalent to perforated sheets from a "universality class" standpoint. We utilize our model to gain insight into scattering measurements made on aggrecan solutions.

Structure and conformational properties of ideal nanogel particles in athermal solutions.

We investigate the conformational properties of "ideal" nanogel particles having a lattice network topology by molecular dynamics simulations to quantify the influence of polymer topology on the solution properties of this type of branched molecular architecture. In particular, we calculate the mass scaling of the radius of gyration (Rg), the hydrodynamic radius, as well as the intrinsic viscosity with the variation of the degree of branching, the length of the chains between the branched points, and the average mesh size within these nanogel particles under good solvent conditions. We find competing trends between the molecular characteristics, where an increase in mesh size or degree of branching results in the emergence of particle-like characteristics, while an increase in the chain length enhances linear polymer-like characteristics. This crossover between these limiting behaviors is also apparent in our calculation of the form factor, P(q), for these structures. Specifically, a primary scattering peak emerges, characterizing the overall nanogel particle size. Moreover, a distinct power-law regime emerges in P(q) at length scales larger than the chain size but smaller than Rg of the nanogel particle, and the Rg mass scaling exponent progressively approaches zero as the mesh size increases, the same scaling as for an infinite network of Gaussian chains. The "fuzzy sphere" model does not capture this feature, and we propose an extension to this popular model. These structural features become more pronounced for values of molecular parameters that enhance the localization of the branching segments within the nanogel particle.

Comparative experimental and computational study of synthetic and natural bottlebrush polyelectrolyte solutions.

We systematically investigate model synthetic and natural bottlebrush polyelectrolyte solutions through an array of experimental techniques (osmometry and neutron and dynamic light scattering) along with molecular dynamics simulations to characterize and contrast their structures over a wide range of spatial and time scales. In particular, we perform measurements on solutions of aggrecan and the synthetic bottlebrush polymer, poly(sodium acrylate), and simulations of solutions of highly coarse-grained charged bottlebrush molecules having different degrees of side-branch density and inclusion of an explicit solvent and ion hydration effects. While both systems exhibit a general tendency toward supramolecular organization in solution, bottlebrush poly(sodium acrylate) solutions exhibit a distinctive "polyelectrolyte peak" in their structure factor, but no such peak is observed in aggrecan solutions. This qualitative difference in scattering properties, and thus polyelectrolyte solution organization, is attributed to a concerted effect of the bottlebrush polymer topology and the solvation of the polymer backbone and counterions. The coupling of the polyelectrolyte topological structure with the counterion distribution about the charged polymer molecules along with direct polymer segmental hydration makes their solution organization and properties "tunable," a phenomenon that has significant ramifications for biological function and disease as well as for numerous materials applications.

Design of nearly perfect hyperuniform polymeric materials.

Disordered hyperuniform materials are exotic amorphous systems that simultaneously exhibit anomalous suppression of long-range density fluctuations, comparable in amplitude to that of crystals and quasi-crystalline materials, while lacking the translational order characteristic of simple liquids. We establish a framework to quantitatively predict the emergence of hyperuniformity in polymeric materials by considering the distribution of localized polymer subregions, instead of considering the whole material. We demonstrate that this highly tunable approach results in arbitrarily small long-range density fluctuations in the liquid state. Our simulations also indicate that long-ranged density fluctuation of the whole polymeric material is remarkably insensitive to molecular topology (linear chain, unknotted ring, star, and bottlebrush) and depends on temperature in an apparently near universal fashion. Our findings open the way for the creation of nearly perfect hyperuniform polymeric materials.

Influence of polymer topology on crystallization in thin films.

We investigate how varying molecular topology of polymers influences crystallization in thin polymer films. In particular, we simulate linear and star polymers of fixed mass having a progressively increasing number of arms (f ≤ 16) in a system where the linear polymer exhibits crystallization in a thin film geometry, but no apparent crystallization in the corresponding bulk material. The degree of crystallization of the polymer film at long times decreases progressively with increasing f, and no crystallization is observed beyond f = 8. Crystallization for smaller values of f develops as a sigmoidally shaped wavefront initiating from the supporting crystalline interface. We suggest that large shape fluctuations and the competition of length scales of star polymers with high f lead to inhibited crystallization.

Dynamic heterogeneity and collective motion in star polymer melts.

While glass formation of linear chain polymer melts has widely been explored, comparatively little is known about glass formation in star polymer melts. We study the segmental dynamics of star polymer melts via molecular dynamics simulations and examine the cooperative nature of segmental motion in star melts. In particular, we quantify how the molecular architecture of star polymers, i.e., the number of arms and the length of those arms, affects the glass transition temperature Tg, the non-Gaussian nature of molecular displacements, the collective string-like motion of monomers, and the role of chain connectivity in the cooperative motion. Although varying the number of arms f and the molecular mass Ma of the star arms can significantly influence the average star molecular shape, all our relaxation data can be quantitatively described in a unified way by the string model of glass formation, an activated transport model that derives from the Adam-Gibbs model, where the degree of cooperative motion is identified with the average length L of string-like particle exchange motions observed in our simulations. Previous work has shown the consistency of the string model with simulations of linear polymers at constant volume and constant pressure, as well as for thin supported polymer films and nanocomposites with variable polymer-surface interactions, where there are likewise large mobility gradients as in the star polymer melts studied in the present paper.

Systematic investigation of synthetic polyelectrolyte bottlebrush solutions by neutron and dynamic light scattering, osmometry, and molecular dynamics simulation.

There is a great interest in the synthesis and characterization of polyelectrolytes that mimic naturally occurring bottlebrush polyelectrolytes to capitalize on the unique properties of this class of macromolecules. Charged bottlebrush polymers form the protective mucus layer in the lungs, stomach, and orifices of animals and provide osmotic stabilization and lubrication to joints. In the present work, we systematically investigate bottlebrush poly(sodium acrylates) through a combination of measurements of solution properties (osmometry, small-angle neutron scattering, and dynamic light scattering) and molecular dynamics simulations, where the bottlebrush properties are compared in each case to their linear polymer counterparts. These complementary experimental and computational methods probe vastly different length- and timescales, allowing for a comprehensive characterization of the supermolecular structure and dynamics of synthetic polyelectrolyte bottlebrush molecules in solution.

Superionic UO2: A model anharmonic crystalline material.

Crystalline materials at elevated temperatures and pressures can exhibit properties more reminiscent of simple liquids than ideal crystalline materials. Superionic crystalline materials having a liquidlike conductivity σ are particularly interesting for battery, fuel cell, and other energy applications, and we study UO2 as a prototypical superionic material since this material is widely studied given its commercial importance as reactor fuel. Using molecular dynamics, we first investigate basic thermodynamic and structural properties. We then quantify structural relaxation, dynamic heterogeneity, and average ion mobility. We find that the non-Arrhenius diffusion and structural relaxation time of this prototypical superionic material can be described in terms of a generalized activated transport model ("string model") in which the activation energy varies with the average string length. Our transport data can also be described equally well by an Adam-Gibbs model in which the excess entropy density of the crystalline material is estimated from specific heat and thermal expansion data, consistent with the average scale of stringlike collective motion scaling inversely with the excess entropy of the crystal. Strong differences in the temperature dependence of the interfacial mobility from nonionic materials are observed, and we suggest that this difference is due to the relatively high cohesive interaction of ionic materials. In summary, the study of superionic UO2 provides insight into the role of cooperative motion in enhancing ion mobility in ionic materials and offers design principles for the development of new superionic materials for use in diverse energy applications.

Hidden Hyperuniformity in Soft Polymeric Materials.

We investigate the nature of long-range density fluctuations in melts of model "soft" polymers, specifically stars and bottlebrushes, over a wide temperature range by molecular dynamics simulation. The cores of the stars and the backbones of bottlebrush polymers are found to have a hyperuniform distribution; i.e., they exhibit anomalously small density fluctuations over a wide temperature range above the glass transition temperature. The hyperuniformity of these substituent polymer subregions is hidden since the fluid as a whole does not exhibit this property. These findings offer a strategy for the practical design of hyperuniform polymeric materials.

Polyelectrolyte association and solvation.

There has been significant interest in the tendency of highly charged particles having the same charge to form dynamic clusters in solution, but an accepted theoretical framework that can account for this ubiquitous phenomenon has been slow to develop. The theoretical difficulties are especially great for flexible polyelectrolytes due to the additional complex coupling between the polyelectrolyte chain configurations and the spatial distribution of the ionic species in solution. For highly charged polyelectrolytes, this leads to the formation of a diffuse "polarizable" cloud of counter-ions around these polymers, an effect having significant implications for the function of proteins and other natural occurring polyelectrolytes, as emphasized long ago by Kirkwood and co-workers. To investigate this phenomenon, we perform molecular dynamics simulations of a minimal model of polyelectrolyte solutions that includes an explicit solvent and counter-ions, where the relative affinity of the counter-ions and the polymer for the solvent is tunable through the variation of the relative strength of the dispersion interactions of the polymer and ions. In particular, we find that these dispersion interactions can greatly influence the nature of the association between the polyelectrolyte chains under salt-free conditions. We calculate static and dynamic correlation functions to quantify the equilibrium structure and dynamics of these complex liquids. Based on our coarse-grained model of polyelectrolyte solutions, we identify conditions in which three distinct types of polyelectrolyte association arise. We rationalize these types of polyelectrolyte association based on the impact of the selective solvent affinity on the charge distribution and polymer solvation in these solutions. Our findings demonstrate the essential role of the solvent in the description of the polyelectrolyte solutions, as well as providing a guideline for the development of a more predictive theory of the properties of the thermodynamic and transport properties of these complex fluids.

A comparative study of thermodynamic, conformational, and structural properties of bottlebrush with star and ring polymer melts.

Thermodynamic, conformational, and structural properties of bottlebrush polymer melts are investigated with molecular dynamics simulations and compared to linear, regular star, and unknotted ring polymer melts to gauge the influence of molecular topology on polymer melt properties. We focus on the variation of the backbone chain length, the grafting density along the backbone, and the length of the side chains at different temperatures above the melt glass transition temperature. Based on these comparisons, we find that the segmental density, isothermal compressibility, and isobaric thermal expansion of bottlebrush melts are quantitatively similar to unknotted ring polymer melts and star polymer melts having a moderate number ( f = 5 to 6) of arms. These similarities extend to the mass scaling of the chain radius of gyration. Our results together indicate that the configurational properties of bottlebrush polymers in their melt state are more similar to randomly branched polymers than linear polymer chains. We also find that the average shape of bottlebrush polymers having short backbone chains with respect to the side chain length is also rather similar to the unknotted ring and moderately branched star polymers in their melt state. As a general trend, the molecular shape of bottlebrush polymers becomes more spherically symmetric when the length of the side chains has a commensurate length as the backbone chain. Finally, we calculate the partial static structure factor of the backbone segments and we find the emergence of a peak at the length scales that characterizes the average distance between the backbone chains. This peak is absent when we calculate the full static structure factor. We characterize the scaling of this peak with parameters characterizing the bottlebrush molecular architecture to aid in the experimental characterization of these molecules by neutron scattering.

Free Surface Relaxations of Star-Shaped Polymer Films.

The surface relaxation dynamics of supported star-shaped polymer thin films are shown to be slower than the bulk, persisting up to temperatures at least 50 K above the bulk glass transition temperature T_{g}^{bulk}. This behavior, exhibited by star-shaped polystyrenes with functionality f=8 arms and molecular weights per arm M_{arm}

Influence of polymer architectures on diffusion in unentangled polymer melts.

Recent simulations have indicated that the thermodynamic properties and the glassy dynamics of polymer melts are strongly influenced by the average molecular shape, as quantified by the radius of gyration tensor of the polymer molecules, and that the average molecular shape can be tuned by varying the molecular topology (e.g., ring, star, linear chain, etc.). In the present work, we investigate if the molecular shape is similarly a predominant factor in understanding the polymer center of mass diffusion D in the melt, as already established for polymer solutions. We find that all our D data for unentangled polymer melts having a range of topologies can be reasonably described as a power law of the polymer hydrodynamic radius, Rh. In particular, this scaling is similar to the scaling of D for a tracer sphere having a radius on the order of the chain radius of gyration, Rg. We conclude that the chain topology influences the molecular dynamics in as much as the polymer topology influences the average molecular shape. Experimental evidence seems to suggest that this situation is also true for entangled polymer melts.

Solution properties of star polyelectrolytes having a moderate number of arms.

We investigate polyelectrolyte stars having a moderate number of arms by molecular dynamics simulations of a coarse-grained model over a range of polyelectrolyte concentrations, where both the counter-ions and solvent are treated explicitly. This class of polymeric materials is found to exhibit rather distinct static and dynamic properties from linear and highly branched star polyelectrolyte solutions emphasized in past studies. Moderately branched polymers are particle-like in many of their properties, while at the same time they exhibit large fluctuations in size and shape as in the case of linear chain polymers. Correspondingly, these fluctuations suppress crystallization at high polymer concentrations, leading apparently to an amorphous rather than crystalline solid state at high polyelectrolyte concentrations. We quantify the onset of this transition by measuring the polymer size and shape fluctuations of our model star polyelectrolytes and the static and dynamic structure factor of these solutions over a wide range of polyelectrolyte concentration. Our findings for star polyelectrolytes are similar to those of polymer-grafted nanoparticles having a moderate grafting density, which is natural given the soft and highly deformable nature of both of these "particles."

Communication: Counter-ion solvation and anomalous low-angle scattering in salt-free polyelectrolyte solutions.

We investigate the influence of counter-ion solvation on the homogeneity of salt-free polyelectrolyte solutions based on a coarse-grained model that includes an explicit solvent. We show that the solvation of the counter-ions can cause a transformation between a nearly homogeneous to a non-uniform polymer solution, in which there is both a chain clustering and the formation of large charge-free domains, i.e., "voids." The emergence of these heterogeneous structures induced by counter-ion solvation is accompanied by the localization and formation of counter-ion rich domains that are symptomatic of emergent effective long-range attractive interchain interactions.

Counter-ion distribution around flexible polyelectrolytes having different molecular architecture.

We explore the monovalent counter-ion distribution around flexible highly-charged polyelectrolytes with different molecular architectures (linear chains, stars, and unknotted and trefoil rings) using molecular dynamics simulations that include an explicit solvent that interacts with the polyelectrolyte. In particular, we find that the molecular topology influences the fraction of counter-ions transiently associating with the polyelectrolyte on a scale of the order of the chain segments, forming a "condensed" counter-ion interfacial layer. As with the hydrogen bonding of water to proteins and other polymers, the persistence time of these interfacial "bound" counter-ions is relatively short, O(1 ps), and we characterize the fluctuations in the number of the counter-ions populating the interfacial layer. We also find that the counter-ions are distributed in a non-uniform fashion on the polyelectrolyte backbone, forming dynamical clusters whose form and average size is sensitive to molecular architecture. In addition, we find that the residual bound counter-ions, not located in either the interfacial layer or the bulk solution, form a diffuse ionic cloud around the polyelectrolyte due to the uncompensated polyelectrolyte charge along the backbone. Generally charge valence strongly influences the extent of the diffuse counter-ion cloud, but in the case of monovalent counter-ions, we find that the size of the diffuse counter-ion cloud nearly coincides with the polyelectrolyte radius of gyration, independent of molecular topology.

Self-assembly of polymer-grafted nanoparticles in solvent-free conditions.

Grafting of polymer chains onto the surface of spherical nanoparticles leads to a hybrid type of fluid that exhibits properties of both particle suspensions and melts of star polymers-these properties being controlled by the relative dimensions of the grafted polymer chains to the nanoparticle diameter, D, and the number of the number of chains grafted on the nanoparticle surface, f. While polymer-grafted nanoparticles (GNP) of this kind typically have a spherical average shape after grafting even a moderate number of chains, their instantaneous molecular shape can fluctuate significantly due to the deformation of the grafted chains. Both simulations and measurements have previously revealed that these "conformationally polarizable" particles can exhibit self-assembly into large scale polymeric structures in both solution and in polymer melts, and we simulate polymer-grafted nanoparticles with D and temperature (T) variations without a dispersing solvent to better understand the nature of this self-assembly process. We observe a reversible self-assembly into linear and branched dynamic GNP structures, where the extent of the assembly and geometry depend on D and T, and we constructed a map capturing the GNP structural behavior with D and T variations. Since the shape of the GNPs appeared to be correlated with the occurrence of the GNP self-assembly, we quantified the average shape and a measure of shape fluctuations to better understand how molecular shape influences their propensity to self-assemble into different structural forms. Based on this framework, we describe the clustering process of the GNPs as an equilibrium polymerization phenomenon and we calculate the order parameter governing the dynamic clustering behavior of GNPs, the average mass of the clusters, size distribution, and the apparent fractal dimension of the clusters.