Theory of kinetic arrest, elasticity, and yielding in dense binary mixtures of rods and spheres

Theoretical study of the impact of dilute nanoparticle additives on the shear elasticity of dense colloidal suspensions

Motivated by basic issues in soft matter physics and new experimental work on granule-nanoparticle mixtures, we systematically apply naive mode coupling theory with accurate microstructural input to investigate the elastic shear modulus of highly size asymmetric, dense, chemically complex, colloid-nanoparticle mixtures. Our analysis spans four equilibrium microstructural regimes: (i) entropic depletion induced colloid clustering, (ii) discrete adsorbed nanoparticle layers that induce colloid spatial dispersion, (iii) nanoparticle-mediated tight bridging network formation, and (iv) colloidal contact aggregation via direct attractions. Each regime typically displays a distinctive mechanical response to changing colloid-nanoparticle size ratio, packing fractions, and the strength and spatial range of interparticle attractive and repulsive interactions. Small concentrations of nanoparticles can induce orders of magnitude elastic reinforcements typically involving single or double exponential growth with increasing colloid and/or nanoparticle packing fraction. Depending on the system, the elementary stress scale can be controlled by the colloid volume, the nanoparticle volume, or a combination of both. Connections between local microstructural organization and the mixture elastic shear modulus are established. The collective structure factor of the relatively dilute nanoparticle subsystem exhibits strong spatial ordering and large osmotic concentration fluctuations imprinted by the highly correlated dense colloidal subsystem. The relevance of the theoretical results for experimental mixtures with large size asymmetry, particularly in the context of 3D ink printing and additive manufacturing, are discussed.

Gum Arabic induced assembly of cellulose nanocrystals in aqueous media

Entropy-driven assembly of nematic liquid-crystal phases of cellulose nanocrystals (SCNCs) in aqueous suspensions results in the emergence of a cholesteric liquid crystalline phase (N* phase). We report that a solvated, non-adsorbing, highly branched natural polysaccharide, Gum Arabic (GA), strongly affects the assembly of the SCNCs and modifies the phase diagram: GA leads to significant crowding of the SCNC rods and induces a new liquid-liquid phase transition, where SCNC-rich and GA-rich droplets coexist. The solvated GA does not induce coagulation or gelation of the suspended SCNCs (at low concentrations of 1-3 wt% of GA). In the SCNC-rich droplets, finite-sized nematic nano-islands assemble and further evolve into cholesteric tactoids and nucleate the formation of the N* phase at significantly lower concentration (about 1.5 wt%) than in GA-free suspensions. We observe that the inter-particle distance and the chiral pitch of the N* phase are determined by the concentration of GA (for a given SCNC concentration). The resulting mesophases are characterized via transmission electron microscopy at cryogenic temperatures (cryo-TEM), small-angle X-ray scattering (SAXS), and polarized optical microscopy (POM). Our findings indicate that GA can be used to tune the phase diagram and optical properties of SCNC suspensions, and overcome kinetic barriers that lead to gelation or kinetic arrest.

Depletion forces in dense mixtures of spheres and rods

We evaluate depletion forces in molecular dynamics simulation of a binary mixture of spheres (depleted particles) and rods (depletant particles) for a wide range of densities for both species. This evaluation was carried out using a recently proposed least squares fitting algorithm. We found that the restriction of the rods' rotational degrees of freedom, when the distance between two spheres is less than the rods length, creates a shallow, and apparently linear, attractive force ramp. For intersphere distances smaller than the rods' diameter, a much stronger attractive force is found, and a large repulsive barrier appears between these aforementioned regimes, roughly at the distance of the rods' thickness. The evaluated forces are validated via a comparison of the pairwise correlation functions obtained from molecular dynamics simulation of a mono-disperse sphere fluid, using the evaluated effective forces, against the original (full system) pairwise correlation functions. Agreement is excellent. We also record the angular pairwise correlation function, using the P2(x) Legendre polynomial, and find that for high densities of both species, a local nematic ordering starts to appear. This nematic order may be a factor in the small differences found between original and effective pairwise correlation functions at high densities of rods.

Theory of the center-of-mass diffusion and viscosity of microstructured and variable sequence copolymer liquids

Biomolecular condensates formed through the phase separation of proteins and nucleic acids are widely observed, offering a fundamental means of organizing intracellular materials in a membrane-less fashion. Traditionally, these condensates have been regarded as homogeneous isotropic liquids. However, in analogy with some synthetic copolymer systems, our recent theoretical research has demonstrated that model biomolecular condensates can exhibit a microemulsion-like internal structure, contingent upon the specific sequence, inter-chain site-site interactions, and concentrated phase polymer density. Motivated by these considerations, here we present a microscopic dynamical theory for the self-diffusion constant and viscosity of a simpler class of model systems - concentrated unentangled A/B regular multiblock copolymer solutions. Our approach integrates static equilibrium local and microdomain scale structural information obtained from PRISM integral equation theory and the time evolution of the autocorrelation function of monomer scale forces at the center-of-mass level to determine the polymer diffusion constant and viscosity in a weak caging regime far from a glass or gel transition. We focus on regular multi-block systems both for simplicity and for its relevance to synthetic macromolecular science. The impact of sequence and inter-chain attraction strength on the slowing down of copolymer mass transport and flow due to local clustering enhanced collisional friction and emergent microdomain scale ordering are established. Analytic analysis and metrics employed in the study of biomolecular condensates are employed to identify key order parameters that quantity how attractive forces, packing structure, multiblock sequence, and copolymer density determine dynamical slowing down above and below the crossover to a fluctuating polymeric microemulsion state.

Cooperative activated hopping dynamics in binary glass-forming liquids: effects of the size ratio, composition, and interparticle interactions

Slow dynamics in supercooled and glassy liquids is an important research topic in soft matter physics. Compared to the traditionally focused one-component systems, glassy dynamics in mixture systems adds in a rich set of new complexities, which are fundamentally interesting and also relevant for many technological applications. In this paper, we apply the recently developed self-consistent cooperative hopping theory (SCCHT) to systematically investigate the effects of the size ratio, composition and interparticle interactions on the cooperative activated hopping dynamics of matrix (in larger size) and penetrant (in smaller size) particles in varied binary sphere mixture model systems, with a specific focus on ultrahigh mixture packing fractions that mimic the deeply supercooled glass transition conditions for molecular/polymeric mixture materials. Analysis shows that in these high activation barrier cases, the long-range elastic distortion associated with a matrix particle hopping over its cage confinement always generates an elastic barrier of a nonnegligible magnitude, although the ratio between the elastic barrier and local barrier contribution is sensitively dependent on all three mixture-specific system factors considered in this work. SCCHT predicts two general scenarios of penetrant-matrix cooperative activated hopping dynamics: matrix/penetrant co-hopping (regime 1) or the penetrant mean barrier hopping time shorter than that of the matrix (regime 2). Increasing the penetrant-to-matrix size ratio or the penetrant-matrix cross-attraction strength is found to universally enlarge the composition window of regime 1. Diverse dynamical properties characterising different aspects of the cooperative activated hopping process, including the penetrant and matrix transient localization lengths, penetrant and matrix hopping jump distances, different types of local and elastic activated barriers, and matrix long-time diffusivity, relaxation time and dynamic fragility are quantitatively studied against a wide range of variations over the three system factors. Of particular interest is the universal "anti-plasticization" phenomenon achievable for sufficiently strong cross-attractive interactions. The prospects this work opens for the exploration of a wide variety of polymer-based mixture materials are briefly discussed at the end.

Theory of polymer diffusion in polymer-nanoparticle mixtures: effect of nanoparticle concentration and polymer length

The dynamics of polymer-nanoparticle (NP) mixtures, which involves multiple scales and system-specific variables, has posed a long-standing challenge on its theoretical description. In this paper, we construct a microscopic theory for polymer diffusion in mixtures based on a combination of the generalized Langevin equation, mode-coupling approach, and polymer physics ideas. The parameter-free theory has an explicit expression and remains tractable on a pair correlation level with system-specific equilibrium structures as input. Taking a minimal polymer-NP mixture as an example, our theory correctly captures the dependence of polymer diffusion on NP concentration and average interparticle distance. Importantly, the polymer diffusion exhibits a power law decay as the polymer length increases at dense NPs and/or a long chain, which marks the emergence of entanglement-like motion. The work provides a first-principles theoretical foundation to investigate dynamic problems in diverse polymer nanocomposites.

Theory of microstructure-dependent glassy shear elasticity and dynamic localization in melt polymer nanocomposites

We present an integrated theoretical study of the structure, thermodynamic properties, dynamic localization, and glassy shear modulus of melt polymer nanocomposites (PNCs) that spans the three microstructural regimes of entropic depletion induced nanoparticle (NP) clustering, discrete adsorbed layer driven NP dispersion, and polymer-mediated bridging network. The evolution of equilibrium and dynamic properties with NP loading, total packing fraction, and strength of interfacial attraction is systematically studied based on a minimalist model. Structural predictions of polymer reference interaction site model integral equation theory are employed to establish the rich behavior of the interfacial cohesive force density, surface excess, and a measure of free volume as a function of PNC variables. The glassy dynamic shear modulus is predicted to be softened, reinforced, or hardly changed relative to the pure polymer melt depending on system parameters, as a result of the competing and qualitatively different influences of interfacial cohesion (physical bonding), free volume, and entropic depletion on dynamic localization and shear elasticity. The localization of polymer segments is the dominant factor in determining bulk PNC softening and reinforcement effects for moderate to strong interfacial attractions, respectively. While in the athermal entropy-dominated regime, the primary origin of mechanical reinforcement is the stress stored in the aggregated NP subsystem. The PNC shear modulus is often qualitatively correlated with the segment localization length but with notable exceptions. The present work provides the foundation for developing a theory of segmental relaxation, T_g changes, and collective NP dynamics in PNCs based on a self-consistent treatment of the cooperative activated motions of segments and NPs.