ABSTRACT
Conformationally fluctuating, globally compact macromolecules such as polymeric rings, single-chain nanoparticles, microgels, and many-arm stars display complex dynamic behaviors due to their rich topological structure and intermolecular organization. Synthetic rings are hybrid objects with conformations that display both ideal random walk and compact globular features, which can serve as models of genomic DNA. To date, emphasis has been placed on the effect of ring molecular weight on their unusual behaviors. Here, we combine simulations and a microscopic force-level theory to build a unified understanding for how key aspects of ring dynamics depend on different tunable molecular properties including backbone rigidity, monomer concentration, degree of traditional entanglement, and molecular weight. Our large-scale molecular dynamics simulations of ring melts with very different backbone stiffnesses reveal unanticipated behaviors which agree well with our generalized theory. This includes a universal master curve for center-of-mass diffusion constants as a function of molecular weight scaled by a chemistry and thermodynamic state-dependent critical molecular weight that generalizes the concept of an entanglement cross-over for linear chains. The key physics is how backbone rigidity and monomer concentration induced changes of the entanglement length, interring packing, degree of interpenetration, and liquid compressibility slow down space-time dynamic-force correlations on macromolecular scales. A power law decay of the center-of-mass diffusion constant with inverse molecular weight squared is the first consequence, followed by an ultraslow activated hopping transport regime. Our results set the stage to address slow dynamics and kinetic arrest in different families of compact synthetic and biological polymeric systems.
ABSTRACT
Understanding the activated transport of penetrant or tracer atoms and molecules in condensed phases is a challenging problem in chemistry, materials science, physics, and biophysics. Many angstrom- and nanometer-scale features enter due to the highly variable shape, size, interaction, and conformational flexibility of the penetrant and matrix species, leading to a dramatic diversity of penetrant dynamics. Based on a minimalist model of a spherical penetrant in equilibrated dense matrices of hard spheres, a recent microscopic theory that relates hopping transport to local structure has predicted a novel correlation between penetrant diffusivity and the matrix thermodynamic dimensionless compressibility, S0(T) (which also quantifies the amplitude of long wavelength density fluctuations), as a consequence of a fundamental statistical mechanical relationship between structure and thermodynamics. Moreover, the penetrant activation barrier is predicted to have a factorized/multiplicative form, scaling as the product of an inverse power law of S0(T) and a linear/logarithmic function of the penetrant-to-matrix size ratio. This implies an enormous reduction in chemical complexity that is verified based solely on experimental data for diverse classes of chemically complex penetrants dissolved in molecular and polymeric liquids over a wide range of temperatures down to the kinetic glass transition. The predicted corollary that the penetrant diffusion constant decreases exponentially with inverse temperature raised to an exponent determined solely by how S0(T) decreases with cooling is also verified experimentally. Our findings are relevant to fundamental questions in glassy dynamics, self-averaging of angstrom-scale chemical features, and applications such as membrane separations, barrier coatings, drug delivery, and self-healing.
Subject(s)
Glass , Physics , Diffusion , Glass/chemistry , Phase Transition , ThermodynamicsABSTRACT
We extend the force-level elastically collective nonlinear Langevin equation theory to treat the spatial gradients of the alpha relaxation time and glass transition temperature, and the corresponding film-averaged quantities, to the geometrically asymmetric case of finite thickness supported films with variable fluid-substrate coupling. The latter typically nonuniversally slows down motion near the solid-liquid interface as modeled via modification of the surface dynamic free energy caging constraints that are spatially transferred into the film and which compete with the accelerated relaxation gradient induced by the vapor interface. Quantitative applications to the foundational hard sphere fluid and a polymer melt are presented. The strength of the effective fluid-substrate coupling has very large consequences for the dynamical gradients and film-averaged quantities in a film thickness and thermodynamic state dependent manner. The interference of the dynamical gradients of opposite nature emanating from the vapor and solid interfaces is determined, including the conditions for the disappearance of a bulk-like region in the film center. The relative importance of surface-induced modification of local caging vs the generic truncation of the long range collective elastic component of the activation barrier is studied. The conditions for the accuracy and failure of a simple superposition approximation for dynamical gradients in thin films are also determined. The emergence of near substrate dead layers, large gradient effects on film-averaged response functions, and a weak non-monotonic evolution of dynamic gradients in thick and cold films are briefly discussed. The connection of our theoretical results to simulations and experiments is briefly discussed, as is the extension to treat more complex glass-forming systems under nanoconfinement.
ABSTRACT
Based on integrating microscopic statistical mechanical theories for structure and ideal kinetic arrest at the naive mode coupling level, we study dynamic localization, the linear elastic shear modulus, applied stress induced modulus softening, and the absolute yielding of simple biphasic binary mixtures composed of equal diameter hard and attractive spheres. The kinetic arrest map is a rich function of total packing fraction, strength of attraction, and mixture composition. The gel to attractive ideal glass transition, the degree of glass melting re-entrancy, and the crossover boundary separating repulsive glasses from attractive glasses vary with the mixture composition. Exponential and/or apparent (high) power law dependences of the elastic shear modulus on the total packing fraction are predicted with effective exponents or exponential prefactors that are sensitive to mixture composition and location in the kinetic arrest map. An analysis of the effective mean square force on a tagged particle that induces dynamic localization reveals a compensation effect between structural correlations and degree of particle localization, resulting in the emergence of a weaker dependence of the shear modulus on mixture composition at very high attraction strengths. Based on a microrheologically inspired formulation of how external stress weakens particle localization and the shear modulus, we analyze mechanical-induced modulus softening and absolute yielding, defined as a discontinuous solid-to-fluid stress-induced transition that can occur in either one or two steps. Estimates of the corresponding yield strains predict that the binary mixture becomes more brittle with increasing sticky particle composition and/or attraction strength.
ABSTRACT
We employ polymer integral equation theory to study a simplified model of semiflexible polymerized ionic liquids (PolyILs) that interact via hard core repulsions and short range screened Coulomb interactions. The multi-scale structure in real and Fourier space of PolyILs (ions chosen to mimic Li, Na, K, Br, PF6, and TFSI) are determined as a function of melt density, Coulomb interaction strength, and ion size. Comparisons with a homopolymer melt, a neutral polymer-solvent-like athermal mixture, and an atomic ionic liquid are carried out to elucidate the distinct manner that ions mediate changes of polymer packing, the role of excluded volume effects, and the influence of chain connectivity, respectively. The effect of Coulomb strength depends in a rich manner on ion size and density, reflecting the interplay of steric packing, ion adsorption, and charge layering. Ion-mediated bridging of monomers is found, which intensifies for larger ions. Intermediate range charge layering correlations are characterized by a many-body screening length that grows with PolyIL density, cooling, and Coulomb strength, in disagreement with Debye-Hückel theory, but in accord with experiments. Qualitative differences in the collective structure, including an ion-size-dependent bifurcation of the polymer structure factor peak and pair correlation function, are predicted. The monomer cage order parameter increases significantly, but its collective ion counterpart decreases, as ions become smaller. Such behaviors allow one to categorize PolyILs into two broad classes of small and large ions. Dynamical implications of the predicted structural results are qualitatively discussed.
ABSTRACT
Understanding in a unified manner the generic and chemically specific aspects of activated dynamics in diverse glass-forming liquids over 14 or more decades in time is a grand challenge in condensed matter physics, physical chemistry, and materials science and engineering. Large families of conceptually distinct models have postulated a causal connection with qualitatively different "order parameters" including various measures of structure, free volume, thermodynamic properties, short or intermediate time dynamics, and mechanical properties. Construction of a predictive theory that covers both the noncooperative and cooperative activated relaxation regimes remains elusive. Here, we test using solely experimental data a recent microscopic dynamical theory prediction that although activated relaxation is a spatially coupled local-nonlocal event with barriers quantified by local pair structure, it can also be understood based on the dimensionless compressibility via an equilibrium statistical mechanics connection between thermodynamics and structure. This prediction is found to be consistent with observations on diverse fragile molecular liquids under isobaric and isochoric conditions and provides a different conceptual view of the global relaxation map. As a corollary, a theoretical basis is established for the structural relaxation time scale growing exponentially with inverse temperature to a high power, consistent with experiments in the deeply supercooled regime. A criterion for the irrelevance of collective elasticity effects is deduced and shown to be consistent with viscous flow in low-fragility inorganic network-forming melts. Finally, implications for relaxation in the equilibrated deep glass state are briefly considered.
ABSTRACT
Molecular, polymeric, colloidal, and other classes of liquids can exhibit very large, spatially heterogeneous alterations of their dynamics and glass transition temperature when confined to nanoscale domains. Considerable progress has been made in understanding the related problem of near-interface relaxation and diffusion in thick films. However, the origin of "nanoconfinement effects" on the glassy dynamics of thin films, where gradients from different interfaces interact and genuine collective finite size effects may emerge, remains a longstanding open question. Here, we combine molecular dynamics simulations, probing 5 decades of relaxation, and the Elastically Cooperative Nonlinear Langevin Equation (ECNLE) theory, addressing 14 decades in timescale, to establish a microscopic and mechanistic understanding of the key features of altered dynamics in freestanding films spanning the full range from ultrathin to thick films. Simulations and theory are in qualitative and near-quantitative agreement without use of any adjustable parameters. For films of intermediate thickness, the dynamical behavior is well predicted to leading order using a simple linear superposition of thick-film exponential barrier gradients, including a remarkable suppression and flattening of various dynamical gradients in thin films. However, in sufficiently thin films the superposition approximation breaks down due to the emergence of genuine finite size confinement effects. ECNLE theory extended to treat thin films captures the phenomenology found in simulation, without invocation of any critical-like phenomena, on the basis of interface-nucleated gradients of local caging constraints, combined with interfacial and finite size-induced alterations of the collective elastic component of the structural relaxation process.
ABSTRACT
We generalize and apply the microscopic self-consistent cooperative hopping theory for activated penetrant dynamics in polymer melts and crosslinked networks to address the role of highly variable non-spherical molecular shape. The focus is on vastly different shaped penetrants that have identical space filling volume in order to isolate how non-spherical shape explicitly modifies dynamics over a wide range of temperature down to the kinetic glass transition temperature. The theory relates intramolecular and intermolecular structure and kinetic constraints, and reveals how different solvation packing of polymer monomers around variable shaped penetrants impact penetrant hopping. A highly shape-dependent penetrant activated relaxation, including alpha time decoupling and trajectory level cooperativity of the hopping process, is predicted in the deeply supercooled regime for relatively larger penetrants which is sensitive to whether the polymer matrix is a melt or heavily crosslinked network. In contrast, for smaller size penetrants or at high/medium temperatures the effect of isochoric penetrant shape is relatively weak. We propose an aspect ratio variable that organizes how penetrant shape influences the activated relaxation times, leading to a (near) collapse or master curve. The relative absolute values of the penetrant relaxation time (inverse hopping rate) in polymer melts versus in crosslinked networks are found to be opposite when compared at a common absolute temperature versus when they are compared at a fixed value of distance from the glass transition based on the variable Tg/T with Tg the glass transition temperature. Quantitative comparison with recent diffusion experiments on chemically complex molecular penetrants of variable shape but fixed volume in crosslinked networks reveals good agreement, and testable new predictions are made. Extension of the theoretical approach to more complex systems of high experimental interest are discussed, including applications to realize selective transport in membrane separation applications.
ABSTRACT
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.
ABSTRACT
Biomolecular condensates can form through the liquid-liquid phase separation (LLPS) of proteins and RNAs in cells. However, other states of organization, including mesostructured network microstructures and physical gels, have been observed, the physical mechanism of which are not well understood. We use the Polymer Reference Interaction Site Model liquid state integral equation theory to study the equilibrium behavior of (generally aperiodic in sequence) biomolecular condensates based on a minimal sticker-spacer associating polymer model. The role of polymer packing fraction, sequence, and the strength and range of intermolecular interactions on macromolecular scale spatial organization and phase behavior is studied for typical sticker-spacer sequences. In addition to the prediction of conventional LLPS, a sequence-dependent strongly fluctuating polymeric microemulsion homogeneous state is predicted at high enough concentrations beyond the so-called Lifshitz-like point, which we suggest can be relevant to the dense phase of microstructured biomolecular condensates. New connections between local clustering and the formation of mesoscopic microdomains, the influence of attraction range, compressibility, and the role of spatial correlations across scales, are established. Our results are also germane to understanding the polymer physics of dense solutions of nonperiodic and unique sequence synthetic copolymers and provide a foundation to create new theories for how polymer diffusion and viscosity are modified in globally isotropic and homogeneous dense polymeric microemulsions.
Subject(s)
Biomolecular Condensates , Polymers , Cluster Analysis , Diffusion , RNAABSTRACT
The first paper of this series [J. Chem. Phys. 158, 034103 (2023)] demonstrated that excess entropy scaling holds for both fine-grained and corresponding coarse-grained (CG) systems. Despite its universality, a more exact determination of the scaling relationship was not possible due to the semi-empirical nature. In this second paper, an analytical excess entropy scaling relation is derived for bottom-up CG systems. At the single-site CG resolution, effective hard sphere systems are constructed that yield near-identical dynamical properties as the target CG systems by taking advantage of how hard sphere dynamics and excess entropy can be analytically expressed in terms of the liquid packing fraction. Inspired by classical equilibrium perturbation theories and recent advances in constructing hard sphere models for predicting activated dynamics of supercooled liquids, we propose a new approach for understanding the diffusion of molecular liquids in the normal regime using hard sphere reference fluids. The proposed "fluctuation matching" is designed to have the same amplitude of long wavelength density fluctuations (dimensionless compressibility) as the CG system. Utilizing the Enskog theory to derive an expression for hard sphere diffusion coefficients, a bridge between the CG dynamics and excess entropy is then established. The CG diffusion coefficient can be roughly estimated using various equations of the state, and an accurate prediction of accelerated CG dynamics at different temperatures is also possible in advance of running any CG simulation. By introducing another layer of coarsening, these findings provide a more rigorous method to assess excess entropy scaling and understand the accelerated CG dynamics of molecular fluids.
Subject(s)
Computer Simulation , Entropy , DiffusionABSTRACT
Coarse-grained (CG) models facilitate an efficient exploration of complex systems by reducing the unnecessary degrees of freedom of the fine-grained (FG) system while recapitulating major structural correlations. Unlike structural properties, assessing dynamic properties in CG modeling is often unfeasible due to the accelerated dynamics of the CG models, which allows for more efficient structural sampling. Therefore, the ultimate goal of the present series of articles is to establish a better correspondence between the FG and CG dynamics. To assess and compare dynamical properties in the FG and the corresponding CG models, we utilize the excess entropy scaling relationship. For Paper I of this series, we provide evidence that the FG and the corresponding CG counterpart follow the same universal scaling relationship. By carefully reviewing and examining the literature, we develop a new theory to calculate excess entropies for the FG and CG systems while accounting for entropy representability. We demonstrate that the excess entropy scaling idea can be readily applied to liquid water and methanol systems at both the FG and CG resolutions. For both liquids, we reveal that the scaling exponents remain unchanged from the coarse-graining process, indicating that the scaling behavior is universal for the same underlying molecular systems. Combining this finding with the concept of mapping entropy in CG models, we show that the missing entropy plays an important role in accelerating the CG dynamics.
Subject(s)
Molecular Dynamics Simulation , Water , Entropy , Water/chemistry , Methanol/chemistryABSTRACT
The diffusion of small molecular penetrants through polymeric materials represents an important fundamental problem, relevant to the design of materials for applications such as coatings and membranes. Polymer networks hold promise in these applications because dramatic differences in molecular diffusion can result from subtle changes in the network structure. In this paper, we use molecular simulation to understand the role that cross-linked network polymers have in governing the molecular motion of penetrants. By considering the local, activated alpha relaxation time of the penetrant and its long-time diffusive dynamics, we can determine the relative importance of activated glassy dynamics on penetrants at the segmental scale vs entropic mesh confinement on penetrant diffusion. We vary several parameters, such as the cross-linking density, temperature, and penetrant size, to show that cross-links primarily affect molecular diffusion through the modification of the matrix glass transition, with local penetrant hopping at least partially coupled to the segmental relaxation of the polymer network. This coupling is very sensitive to the local activated segmental dynamics of the surrounding matrix, and we also show that penetrant transport is affected by dynamic heterogeneity at low temperatures. To contrast, only at high temperatures and for large penetrants or when the dynamic heterogeneity effect is weak, does the effect of mesh confinement become significant, even though penetrant diffusion more broadly empirically follows similar trends as established models of mesh confinement-based transport.
ABSTRACT
We generalize a microscopic statistical mechanical theory of the activated dynamics of dilute spherical penetrants in glass-forming liquids to study the influence of crosslinking in polymer networks on the penetrant relaxation time and diffusivity over a wide range of temperature and crosslink fraction (fn). Our calculations are relevant to recent experimental studies of a nm-sized molecule diffusing in poly-(n-butyl methacrylate) networks. The theory predicts the penetrant relaxation time increases exponentially with the glass transition temperature, Tg(fn), which grows roughly linearly with the square root of fn due to the coupling of local hopping to longer-range collective elasticity. Moreover, Tg is also found to be proportional to a geometric confinement parameter defined as the ratio of the penetrant diameter to the mean network mesh size. The decoupling ratio of the penetrant and Kuhn segment alpha times displays a complex non-monotonic dependence on fn and temperature that is well collapsed based on the variable Tg(fn)/T. A model for the penetrant diffusion constant that combines activated relaxation and entropic mesh confinement is proposed, which results in a significantly stronger suppression of mass transport with degree of effective supercooling than predicted for the penetrant alpha time. This behavior corresponds to a new network-based type of "decoupling" of diffusion and relaxation. In contrast to the diffusion of larger nanoparticles in high temperature rubbery networks, our analysis in the supercooled regime suggests that for the penetrants studied the mesh confinement effects are of secondary importance relative to the consequences of crosslink-induced slowing down of activated hopping of glassy physics origin.
ABSTRACT
We combine simulation and Elastically Collective Nonlinear Langevin Equation (ECNLE) theory to study the activated relaxation in monodisperse atomic and polymeric Weeks-Chandler-Andersen (WCA) liquids over a wide range of temperatures and densities in the supercooled regime under isochoric conditions. By employing novel crystal-avoiding simulations, metastable equilibrium dynamics is probed in the absence of complications associated with size polydispersity. Based on a highly accurate structural input from integral equation theory, ECNLE theory is found to describe well the simulated density and temperature dependences of the alpha relaxation time of atomic fluids using a single system-specific parameter, ac, that reflects the nonuniversal relative importance of local cage and collective elastic barriers. For polymer fluids, the explicit dynamical effect of local chain connectivity is modeled at the fundamental dynamic free energy trajectory level based on a different parameter, Nc, that quantifies the degree of intramolecular correlation of bonded segment activated barrier hopping. For the flexible chain model studied, a physically intuitive value of Nc ≈ 2 results in good agreement between simulation and theory. A direct comparison between atomic and polymeric systems reveals that chain connectivity can speed up activated segmental relaxation due to weakening of equilibrium packing correlations but can slow down relaxation due to local bonding constraints. The empirical thermodynamic scaling idea for the alpha time is found to work well at high densities or temperatures but fails when both density and temperature are low. The rich and subtle behaviors revealed from simulation for atomic and polymeric WCA fluids are all well captured by ECNLE theory.
ABSTRACT
We employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant-matrix pair correlation function and the penetrant-matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.
ABSTRACT
We generalize the self-consistent cooperative hopping theory for a dilute spherical penetrant or tracer activated dynamics in dense metastable hard sphere fluids and glasses to address the effect of external stress, the consequences of which are systematically established as a function of matrix packing fraction and penetrant-to-matrix size ratio. All relaxation processes speed up under stress, but the difference between the penetrant and matrix hopping (alpha relaxation) times decreases significantly with stress corresponding to less time scale decoupling. A dynamic crossover occurs at a critical "slaving onset" stress beyond which the matrix activated hopping relaxation time controls the penetrant hopping time. This characteristic stress increases (decreases) exponentially with packing fraction (size ratio) and can be well below the absolute yield stress of the matrix. Below the slaving onset, the penetrant hopping time is predicted to vary exponentially with stress, differing from the power law dependence of the pure matrix alpha time due to system-specificity of the stress-induced changes in the penetrant local cage and elastic barriers. An exponential growth of the penetrant alpha relaxation time with size ratio under stress is predicted, and at a fixed matrix packing fraction, the exponential relation between penetrant hopping time and stress for different size ratios can be collapsed onto a master curve. Direct connections between the short- and long-time activated penetrant dynamics and between the penetrant (or matrix) alpha relaxation time and matrix thermodynamic dimensionless compressibility are also predicted. The presented results should be testable in future experiments and simulations.
ABSTRACT
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, Tg changes, and collective NP dynamics in PNCs based on a self-consistent treatment of the cooperative activated motions of segments and NPs.
ABSTRACT
We theoretically study the effect of external deformation on activated structural relaxation and aspects of the nonlinear mechanical response of glassy hard sphere fluids in the context of elastically collective nonlinear Langevin equation theory. This microscopic force-based approach describes activated relaxation as a coupled local-nonlocal event involving caging and longer range collective elasticity, with the latter becoming more important and ultimately dominant with increasing packing fraction under equilibrium conditions. The central new question we address is how this physical picture of activated relaxation, and the relative importance of local caging vs collective elasticity physics, depends on external deformation. Theoretical predictions are presented for deformation-induced enhancement of mobility, the onset of relaxation speed up at remarkably low values of stress, strain, or shear rate, apparent power law thinning of the steady state structural relaxation time and viscosity, a non-vanishing activation barrier in the shear thinning regime, an apparent Herschel-Bulkley form of the rate dependence of the steady state shear stress, exponential growth of different measures of a dynamic yield or flow stress with the packing fraction, and reduced fragility and dynamic heterogeneity under deformation. The results are contrasted with experiments and simulations, and qualitative or better agreement is found. An overarching conclusion is that deformation strongly reduces the importance of longer range collective elastic effects relative to the local caging aspect for most, but not all, physical questions, with deformation-dependent fragility and dynamic heterogeneity phenomena being qualitatively sensitive to collective elasticity. Overall, nonlinear rheology is predicted to be a more local problem than quiescent structural relaxation, albeit with deformation-modified activated processes still important.
ABSTRACT
We propose a microscopic theory for the decoupling of self-diffusion and structural relaxation in glass-forming liquids within the Elastically Collective Nonlinear Langevin Equation (ECNLE) activated dynamics framework. Our central hypothesis is that the heterogeneity relevant to this problem is static fluctuations of local density on the scale of 3-4 particle diameters and how this changes local packing correlations. These fluctuations modify the degree of dynamical cage expansion that mechanistically couples intracage large amplitude hopping and longer range collective elasticity in ECNLE theory. Decoupling only emerges in the deeply supercooled regime where the strongly temperature dependent elastic barrier becomes non-negligible relative to its noncooperative local analog. The theory makes predictions for various aspects of the decoupling phenomenon, including apparent fractional power law Stokes-Einstein behavior, that appear to be consistent with experiments and simulations on hard sphere fluids and molecular liquids. Of central importance is a microscopic connection between the barrier fluctuation variance and most probable barrier height. Sensible results are also obtained for the nonexponential stretching of a generic relaxation time correlation function and its temperature evolution. Nonuniversality can arise from the relative importance of the local and collective barriers (related to fragility) and the precise magnitude of the length scale that defines the transition from local cage to elastic physics. Comparison is made with a traplike model based on a Gaussian distribution of barriers.