Nanorheology of Entangled Polymer Melts

Diffusion of Nanosheets in Unentangled Polymer Melts

Understanding the dynamics of nanosheets in polymer matrices is crucial for the processing of polymer nanocomposites and their applications in drug delivery. In this work, we investigate the diffusion of thin nanosheets in unentangled polymer melts using molecular dynamics simulations. We show that for nanosheets smaller than a characteristic size lc, which is a few times the polymer chain size, the continuum hydrodynamic theory based on macroscopic viscosity breaks down and significantly underestimates the diffusion coefficients. For nanosheets with sizes l < lc, we derive scaling relationships for both translational and rotational diffusion coefficients as functions of l and further reveal the dynamical coupling between nanosheet motion and the modes of the polymer melt. For l > lc, the continuum theory is recovered. Our findings reconcile the continuum and scaling theories for the diffusion of nanoparticles in polymer melts.

Dynamics of Nanoparticles in Solutions of Semiflexible Ring Polymers

We use hybrid molecular dynamics-multiparticle collision dynamics (MD-MPCD) simulations to investigate the influence of chain stiffness on the transport of nanoparticles (NPs) through solutions of semiflexible ring polymers. The NPs exhibit subdiffusive dynamics on short time scales before transitioning to normal diffusion at longer times. The terminal NP diffusivity decreases with increasing ring stiffness, similar to the behavior observed in solutions of semiflexible linear chains. The NP subdiffusive exponent is found to be strongly correlated with that of the polymer center of mass (COM) for the range of chain stiffnesses examined, which is at odds with the pronounced decoupling of the NP and polymer COM motions previously observed upon increasing the stiffness of linear chains. Our analysis indicates that these marked differences in the intermediate dynamics are rooted in distinct structural changes that emerge with increasing bending stiffness: Stiffer ring polymers adopt increasingly circular conformations and stack into transient tubes. The void space created near the ring centers is occupied by NPs and other polymers, resulting in strong dynamic coupling on short time scales.

Brownian yet non-Gaussian diffusion of a light particle in heavy gas: Lorentz-gas-based analysis

Non-Gaussian diffusion was recently observed in a gas mixture with mass and fraction contrast [F. Nakai et al., Phys. Rev. E 107, 014605 (2023)2470-004510.1103/PhysRevE.107.014605]. The mean-square displacement of a minor gas particle with a small mass is linear in time, while the displacement distribution deviates from the Gaussian distribution, which is called the Brownian yet non-Gaussian diffusion. In this work, we theoretically analyze this case where the mass contrast is sufficiently large. Major heavy particles can be interpreted as immobile obstacles, and a minor light particle behaves like a Lorentz gas particle within an intermediate timescale. Despite the similarity between the gas mixture and the conventional Lorentz gas system, the Lorentz gas description cannot fully describe the Brownian yet non-Gaussian diffusion. A successful description can be achieved through a canonical ensemble average of the statistical quantities of the Lorentz gas over the initial speed. Furhter, we show that the van Hove correlation function has a nonexponential tail, which is contrary to the exponential tail observed in various systems.

Force-driven active dynamics of thin nanorods in unentangled polymer melts

Recent advances in the functional material and biomedical applications of nanorods call for a fundamental understanding of the active motion of nanorods in a viscoelastic medium. Molecular dynamics simulations are performed to investigate a model system consisting of force-driven active thin nanorods in a melt of unentangled polymers. The activeness of a thin nanorod arises from a constant external force applied uniformly along the rod. The simulations demonstrate that the active force overcomes the randomness of the diffusive motion and results in a ballistic motion along the direction of the applied force at long timescales. The constant speed of the force-driven ballistic motion is determined by the balance of the active force and the friction from the coupling of the nanorod with the polymer viscosity. The friction coefficient, which is computed as the ratio of the active force and the speed, decreases as the active force increases. The origin of the reduction in the friction coefficient is the high speed that allows the nanorod to renew its local environment faster than the relaxation time of melt chains. A scaling theory is developed to quantify the dependence of the friction coefficient on the strength of the active force. The simulations also demonstrate that the force-driven ballistic motion suppresses the rotational diffusion of the rod and cuts off the de-correlation of the rod axis with time. On the scaling level, the long-time trajectory of a force-driven active nanorod piercing through unentangled polymers may be described as a stretched array of "active blobs", where the short-time random-walk trajectory within an active blob is unperturbed by the active force.

Diffusion of polymer-grafted nanoparticles with dynamical fluctuations in unentangled polymer melts

The dynamics of polymer-grafted nanoparticles (PGNPs) in melts of unentangled linear chains were investigated by means of coarse-grained molecular dynamics simulations. The results demonstrated that the graft monomers closer to the particle surface relax more slowly than those farther away due to the constraint of the grafted surface and the confinement of the neighboring chains. Such heterogeneous relaxations of the surrounding environment would perturb the particle motion, making them fluctuating around their centers before they can diffuse through the melt. During such intermediate-time stage, the dynamics is subdiffusive while the distribution of particle displacements is Gaussian, which can be described by the popular fractional Brownian motion model. For the long-time Fickian diffusion, we found that the diffusivity D decreases with increasing grafting density Σ_g, grafted chain length N_g, and matrix chain length N_m. This is due to the fact that the diffusivity is controlled by the viscous drag of an effective core, consisting of the NP and the non-draining layer of graft segments, and that of the free-draining graft layer outside the "core". With increasing Σ_g, the PGNPs become harder with greater effective size and thinner free draining layer, resulting in a reduction in D. At extremely high Σ_g, the diffusivity can even be estimated by the diameter-renormalized Stokes-Einstein (SE) relation. With increasing N_g, both the effective core size and the thickness of the free-draining layer increase, leading to a reduction in diffusivity by D ∼ N-γg with 0.5 < γ < 1. Increasing N_m would lead to the enlargement of the effective core size but meanwhile result in the reduction of the free-draining layer thickness due to autophobic dewetting. The counteraction between these two opposite effects leads to only a slight reduction in the diffusivity, significantly different from the typical SE behavior where D ∼ N_m^-1. These findings bear significance in unraveling the fundamental physics of the anomalous dynamics of PGNPs in various polymers, including biological and synthetic.

Molecular Characterization of Polymer Networks

Polymer networks are complex systems consisting of molecular components. Whereas the properties of the individual components are typically well understood by most chemists, translating that chemical insight into polymer networks themselves is limited by the statistical and poorly defined nature of network structures. As a result, it is challenging, if not currently impossible, to extrapolate from the molecular behavior of components to the full range of performance and properties of the entire polymer network. Polymer networks therefore present an unrealized, important, and interdisciplinary opportunity to exert molecular-level, chemical control on material macroscopic properties. A barrier to sophisticated molecular approaches to polymer networks is that the techniques for characterizing the molecular structure of networks are often unfamiliar to many scientists. Here, we present a critical overview of the current characterization techniques available to understand the relation between the molecular properties and the resulting performance and behavior of polymer networks, in the absence of added fillers. We highlight the methods available to characterize the chemistry and molecular-level properties of individual polymer strands and junctions, the gelation process by which strands form networks, the structure of the resulting network, and the dynamics and mechanics of the final material. The purpose is not to serve as a detailed manual for conducting these measurements but rather to unify the underlying principles, point out remaining challenges, and provide a concise overview by which chemists can plan characterization strategies that suit their research objectives. Because polymer networks cannot often be sufficiently characterized with a single method, strategic combinations of multiple techniques are typically required for their molecular characterization.

Interface and Interphase in Polymer Nanocomposites with Bare and Core-Shell Gold Nanoparticles

Metal nanoparticles are used to modify/enhance the properties of a polymer matrix for a broad range of applications in bio-nanotechnology. Here, we study the properties of polymer/gold nanoparticle (NP) nanocomposites through atomistic molecular dynamics, MD, simulations. We probe the structural, conformational and dynamical properties of polymer chains at the vicinity of a gold (Au) NP and a functionalized (core/shell) Au NP, and compare them against the behavior of bulk polyethylene (PE). The bare Au NPs were constructed via a systematic methodology starting from ab-initio calculations and an atomistic Wulff construction algorithm resulting in the crystal shape with the minimum surface energy. For the functionalized NPs the interactions between gold atoms and chemically adsorbed functional groups change their shape. As a model polymer matrix we consider polyethylene of different molecular lengths, from the oligomer to unentangled Rouse like systems. The PE/Au interaction is parametrized via DFT calculations. By computing the different properties the concept of the interface, and the interphase as well, in polymer nanocomposites with metal NPs are critically examined. Results concerning polymer density profiles, bond order parameter, segmental and terminal dynamics show clearly that the size of the interface/interphase, depends on the actual property under study. In addition, the anchored polymeric chains change the behavior/properties, and especially the chain density profile and the dynamics, of the polymer chain at the vicinity of the Au NP.

The Longitudinal Superdiffusive Motion of Block Copolymer in a Tight Nanopore

The structure and dynamic properties of polymer chains in a confined environment were studied by means of the Monte Carlo method. The studied chains were represented by coarse-grained models and embedded into a simple 3D cubic lattice. The chains stood for two-block linear copolymers of different energy of bead-bead interactions. Their behavior was studied in a nanotube formed by four impenetrable surfaces. The long-time unidirectional motion of the chain in the tight nanopore was found to be correlated with the orientation of both parts of the copolymer along the length of the nanopore. A possible mechanism of the anomalous diffusion was proposed on the basis of thermodynamics of the system, more precisely on the free energy barrier of the swapping of positions of both parts of the chain and the impulse of temporary forces induced by variation of the chain conformation. The mean bead and the mass center autocorrelation functions were examined. While the former function behaves classically, the latter indicates the period of time of superdiffusive motion similar to the ballistic motion with the autocorrelation function scaling with the exponent t^5/3. A distribution of periods of time of chain diffusion between swapping events was found and discussed. The influence of the nanotube width and the chain length on the polymer diffusivity was studied.

Short-time dynamics of a tracer in an ideal gas

A small tagged particle immersed in a fluid exhibits Brownian motion and diffuses on a long timescale. Meanwhile, on a short timescale, the dynamics of the tagged particle cannot be simply described by the usual generalized Langevin equation with Gaussian noise, since the number of collisions between the tagged particle and fluid particles is rather small. On such a timescale, we should explicitly consider individual collision events between the tagged particle and the surrounding fluid particles. In this study we analyze the short-time dynamics of a tagged particle in an ideal gas, where we do not have static or hydrodynamic correlations between fluid particles. We perform event-driven hard-sphere simulations and show that the short-time dynamics of the tagged particle is correlated even under such an idealized situation. Namely, the velocity autocorrelation function becomes negative when the tagged particle is relatively light and the fluid density is relatively high. This result can be attributed to the dynamical correlation between collision events. To investigate the physical mechanism which causes the dynamical correlation, we analyze the correlation between successive collision events. We find that the tagged particle can collide with the same ideal-gas particle several times and such collisions cause a strong dynamical correlation for the velocity.

Effects of Tethered Polymers on Dynamics of Nanoparticles in Unentangled Polymer Melts

Polymer-tethered nanoparticles (NPs) are commonly added to a polymer matrix to improve material properties. Critical to the fabrication and processing of such composites is the mobility of the tethered NPs. Here we study the motion of tethered-NPs in unentangled polymer melts using molecular dynamics simulations, which offer a precise control of the grafted chain length N _g and the number z of grafted chains per particle. As N _g increases, there is a crossover from particle-dominated to tethered-chain-dominated terminal diffusion of NPs with the same z. The mean squared displacement of loosely tethered NPs in the case of tethered-chain dominated terminal diffusion exhibits two sub-diffusive regimes at intermediate time scales for small z. The first one at shorter time scales arises from the dynamical coupling of the particle and matrix chains, while the one at longer time scales is due to the participation of the particle in the dynamics of the tethered chains. The friction of loosely grafted chains in unentangled melts scales linearly with the total number of monomers in the chains, as the frictions of individual monomers are additive in the absence of hydrodynamic coupling. As more chains are grafted to a particle, hydrodynamic interactions between grafted chains emerge. As a result, there is a non-draining layer of hydrodynamically coupled chain segments surrounding the bare particle. Outside the non-draining layer is a free-draining layer of grafted chain segments with no hydrodynamic coupling. The boundary of the two layers is the stick surface where the shear stress due to the relative melt flow is balanced by the friction between grafted and melt chains in the interpenetration layer. The stick surface is located further away from the bare surface of the particle with higher grafting density.

Diffusion of Nanoparticles with Activated Hopping in Crowded Polymer Solutions

A long-distance hop of diffusive nanoparticles (NPs) in crowded environments was commonly considered unlikely, and its characteristics remain unclear. In this work, we experimentally identify the occurrence of the intermittent hops of large NPs in crowded entangled poly(ethylene oxide) (PEO) solutions, which are attributed to thermally induced activated hopping. We show that the diffusion of NPs in crowded solutions is considered as a superposition of the activated hopping and the reptation of the polymer solution. Such activated hopping becomes significant when either the PEO molecular weight is large enough or the NP size is relatively small. We reveal that the time-dependent non-Gaussianity of the NP diffusion is determined by the competition of the short-time relaxation of a polymer entanglement strand, the activated hopping, and the long-time reptation. We propose an exponential scaling law τ_hop/τ_e ∼ exp(d/d_t) to characterize the hopping time scale, suggesting a linear dependence of the activated hopping energy barrier on the dimensionless NP size. The activated hopping motion can only be observed between the onset time scale of the short-time relaxation of local entanglement strands and the termination time scale of the long-time relaxation. Our findings on activated hopping provide new insights into long-distance transportation of NPs in crowded biological environments, which is essential to the delivery and targeting of nanomedicines.

Cytoskeleton-inspired artificial protein design to enhance polymer network elasticity

Reducing topological network defects to enhance elasticity in polymeric materials remains a grand challenge. Efforts to control network topology, primarily focused on crosslinking junctions, continue to underperform compared to theoretical estimations from idealized networks using affine and phantom network theories. Here, artificial protein technology was adapted for the design of polymer-network hydrogels with precisely defined coil-like and rod-like strands to observe the impact of strand rigidity on the mechanical properties of polymeric materials. Cytoskeleton-inspired polymer-network hydrogels incorporated with rod-like protein strands nearly tripled the gel shear elastic modulus and relaxation time compared to coil-like protein strands, indicating an enhanced effective crosslinking density. Furthermore, asymmetric rod-coil protein designs in network strands with an optimal rod:coil ratio improved the hydrogel relaxation time, enhancing the stability of physical macromolecular associations by modulating crosslinker mobility. The careful design of strand rigidity presents a new direction to reduce topological defects for optimizing polymeric materials.

A practical method to account for random phase approximation effects on the dynamic scattering of multi-component polymer systems

Investigations of polymer systems that rely on the interpretation of dynamical scattering results as, e.g., the structure factor S(Q, t) of single chains or chain sections may require the inclusion of effects, as described within the framework of the random phase approximation (RPA) for polymers. To do this in practice for the dynamic part of S(Q, t) beyond the initial slope is a challenge. Here, we present a method (and software) that allows a straightforward assessment of dynamical RPA effects and inclusion of these in the process/procedures of model fitting. Examples of applications to the interpretation of neutron spin-echo data multi-component polymer melts are shown.

Investigation of the properties of nanographene in polymer nanocomposites through molecular simulations: dynamics and anisotropic Brownian motion

The dynamical behavior of nanographene sheets dispersed in polymer matrices is investigated through united-atom molecular dynamics simulations. The Brownian motion of the sheet and the anisotropy in its translational and orientational diffusion are the topics of the current study. Different polymer matrices and pristine and functionalized graphene constitute various nanocomposite systems. Interactions between the nanographene flake and the matrix determine the dynamics of the systems. The dynamics is reduced in polyethylene oxide compared to polyethylene matrix, whereas carboxylated sheets move considerably slower than the pristine nanographene in any matrix. Diffusion is anisotropic for short times, while it becomes isotropic in the long time limit. The in-plane motion of the nanographene sheet is faster than the out-of-plane component, in agreement with the diffusion of perfectly oblate ellipsoids. In functionalized graphene, the anisotropy is suppressed. By exploring the temperature effect on both the nanographene sheet and polymer close to the surface, indications for coupling in the motion of the two components are revealed. The strong effect of edge functional groups on the dynamics can be used as a way to control the Brownian motion of nanographene sheets in polymer nanocomposites and consequently tailor the properties of the materials.

Optical Tweezers Microrheology Maps the Dynamics of Strain-Induced Local Inhomogeneities in Entangled Polymers

Optical tweezers microrheology (OTM) offers a powerful approach to probe the nonlinear response of complex soft matter systems, such as networks of entangled polymers, over wide-ranging spatiotemporal scales. OTM can also uniquely characterize the microstructural dynamics that lead to the intriguing nonlinear rheological properties that these systems exhibit. However, the strain in OTM measurements, applied by optically forcing a microprobe through the material, induces network inhomogeneities in and around the strain path, and the resultant flow field complicates the measured response of the system. Through a robust set of custom-designed OTM protocols, coupled with modeling and analytical calculations, we characterize the time-varying inhomogeneity fields induced by OTM measurements. We show that homogenization following strain does not interfere with the intrinsic stress relaxation dynamics of the system, rather it manifests as an independent component in the stress decay, even in highly nonlinear regimes such as with the microrheological large-amplitude-oscillatory-shear (MLAOS) protocols we introduce. Our specific results show that Rouse-like elastic retraction, rather than disentanglement and disengagement, dominates the nonlinear stress relaxation of entangled polymers at micro- and mesoscales. Thus, our study opens up possibilities of performing precision nonlinear microrheological measurements, such as MLAOS, on a wide range of complex macromolecular systems.

Miscibility and Nanoparticle Diffusion in Ionic Nanocomposites

We investigate the effect of various spherical nanoparticles in a polymer matrix on dispersion, chain dimensions and entanglements for ionic nanocomposites at dilute and high nanoparticle loading by means of molecular dynamics simulations. The nanoparticle dispersion can be achieved in oligomer matrices due to the presence of electrostatic interactions. We show that the overall configuration of ionic oligomer chains, as characterized by their radii of gyration, can be perturbed at dilute nanoparticle loading by the presence of charged nanoparticles. In addition, the nanoparticle's diffusivity is reduced due to the electrostatic interactions, in comparison to conventional nanocomposites where the electrostatic interaction is absent. The charged nanoparticles are found to move by a hopping mechanism.