Unconventionally fast transport through sliding dynamics of rodlike particles in macromolecular networks

Transport of rodlike particles in confinement environments of macromolecular networks plays crucial roles in many important biological processes and technological applications. The relevant understanding has been limited to thin rods with diameter much smaller than network mesh size, although the opposite case, of which the dynamical behaviors and underlying physical mechanisms remain unclear, is ubiquitous. Here, we solve this issue by combining experiments, simulations and theory. We find a nonmonotonic dependence of translational diffusion on rod length, characterized by length commensuration-governed unconventionally fast dynamics which is in striking contrast to the monotonic dependence for thin rods. Our results clarify that such a fast diffusion of thick rods with length of integral multiple of mesh size follows sliding dynamics and demonstrate it to be anomalous yet Brownian. Moreover, good agreement between theoretical analysis and simulations corroborates that the sliding dynamics is an intermediate regime between hopping and Brownian dynamics, and provides a mechanistic interpretation based on the rod-length dependent entropic free energy barrier. The findings yield a principle, that is, length commensuration, for optimal design of rodlike particles with highly efficient transport in confined environments of macromolecular networks, and might enrich the physics of the diffusion dynamics in heterogeneous media.


Abstract
Transport of rodlike particles in confinement environments of macromolecular networks plays crucial roles in many important biological processes and technological applications.The relevant understanding has been limited to thin rods with diameter much smaller than network mesh size, although the opposite case, of which the dynamical behaviors and underlying physical mechanisms remain unclear, is ubiquitous.Here, we solve this issue by combining experiments, simulations and theory.We find a nonmonotonic dependence of translational diffusion on rod length, characterized by length commensuration-governed unconventionally fast dynamics which is in striking contrast to the monotonic dependence for thin rods.Our results clarify that such a fast diffusion of thick rods with length of integral multiple of mesh size follows sliding dynamics and demonstrate it to be "anomalous yet Brownian".Moreover, good agreement between theoretical analysis and simulations corroborates that the sliding dynamics is an intermediate regime between hopping and Brownian dynamics, and provides a mechanistic interpretation based on the rod-length dependent entropic free energy barrier.The findings yield a principle, that is, length commensuration, for optimal design of rodlike particles with highly efficient transport in confined environments of macromolecular networks, and might enrich the physics of the diffusion dynamics in heterogeneous media.
The transport of guest particles in confinement environments of macromolecular networks holds significant importance in various fields, including biological science, materials science, and soft matter physics [1][2][3][4][5][6][7][8] .It serves as a fundamental problem that underlies the behavior of diverse systems.For instance, in biological systems, the selective permeation of bacteria or antibiotics across mucus is influenced by the transport dynamics occurring in confined spaces [1][2][3] .Besides, it determines the organization of nanoparticles within polymer-network-based nanocomposites, enabling the development of advanced materials with characteristic properties 4,5 .Developing efficient drug delivery systems that target tissues covered with extracellular matrix also relies on understanding and manipulating guest particle transport in such environments [6][7][8] .Thus, understanding the dynamics of particle transport is essential for advancements in nanocomposite design, medical treatments, and drug delivery strategies with improved efficiency and specificity.
Physically, the dynamics is governed by the corresponding effective free energy landscape contributed by the interaction between particle and networks strands as well as the elastic deformation of network strands 5,[9][10][11][12][13] .For spherical particles, the interplay of these two factors strongly depends on the size ratio between the network mesh and particles, resulting in various diffusion regimes including Brownian, hopping and trapped dynamics [10][11][12][13] .However, this simple picture breaks down for rodlike particles because anisotropic shape gives rise to additional competing length scales [14][15][16][17] .As a matter of fact, numerous particles in biological systems, such as bacteria, viruses, nucleic acids, proteins, and polypeptides, hold rodlike geometry 18 .As various length scales associated with the structure significantly increase the complexity of the dynamics, the fundamental dynamics of such systems is far from established.
In most studies on the transport of rodlike particles in macromolecular networks, rod diameters are much smaller than typical length scales of network mesh [19][20][21][22][23][24] ; the elastic energy due to strand deformation thereby plays a minor role in the effective free energy landscape.However, the opposite limit, with rod diameter comparable with or even larger than the latter, is ubiquitous.For example, to demonstrate this point, we perform a statistical analysis of the ratio of diameter d of rodlike bacteria to averaged mesh size ax, d/ax, for many experiments of typical bacteria in different types of mucus consisting of biomacromolecular networks [24][25][26][27][28][29][30][31][32][33][34] .As demonstrated in Fig. S1 and listed in Table S1, surprisingly, d/ax ranges from about 0.8 to almost 2.0.In this case, the entropic free energy barrier due to the conformational penalty of strands deformed by particles may overwhelm the interactions, leading to new dynamical regimes 35 .Unfortunately, little is known about the diffusion mechanisms regarding the transport of such thick rodike particles in macromolecular networks, leaving an urgent and critical issue to be addressed.
Here, by combining experiment, simulation and theory, we provide, for the first time, the systematical exploration of the dynamics of rodlike particles of diameters comparable to mesh size in a cross-linked macromolecular network, in order to unravel the fundamental mechanisms underlying the transport in the essential systems.The merits of the results are listed as follow: (1) It reports the nonmonotonic dependence of translational diffusion on rod length, where exceptional fast dynamics occurs once rod length reaches around an integral multiple of network mesh size, in stark contrast to the monotonic dependence reported previously for thin rods 20 .(2) It not only reveals that the fast transport follows sliding dynamics but also gives the analytical expressions of the time-displacement distribution of sliding dynamics, clarifying its physical relationship with hopping and Brownian dynamics.(3) This work provides new principles for the optimal design of particle transport in various networks, biological or synthetic.

Single-particle tracking of rods in synthetic networks
We begin by examining the diffusion behavior of rods with different lengths in a synthetic macromolecular network, where the diameters of rods d are comparable with the mesh size ax of the network.The experimental network is polyethylene glycol diacrylate (PEGDA) network, which possesses excellent biocompatibility 36 .The PEGDA network is synthesized under ultraviolet (UV) irradiation (see Methods and Supplementary Information I for more details) 37 .ax of this network is estimated to be 21.0 ± 1.8 nm (Table S2) 38,39 .The experimental rods are PEG-capped Au nanorods (Au-NRs), with lengths L ranging from 30.6 ± 3.4 nm to 61.4 ± 7.7 nm and diameters, including grafted-layer thickness, of around 19.1 nm.Rod sizes are measured from TEM images (Fig. S2), and the averaged values are listed in Table S3.The images of trajectories are taken by using dark field techniques, as schemed in Fig. 1a [40][41][42][43] .Particularly, in the confinement environment due to the scales of such large diameter and length of the rods considered in the current work, the waiting time for a rotating event possesses a period much larger than the observation time, underscoring that the rotational dynamics plays a trivial role in the transport of rods in the experiments (see Supplementary Information V for more details).Thus, in the following discussion, we focus on the transitional dynamics of the rods.
To detect the particles in the image stack, the Crocker-Grier algorithm is applied in the experiments 44,45 .Representative snapshots of rods are illustrated in Fig. 1b, where the initial and spontaneous positions of the rod centers are marked by the dashed cyan and magenta lines, respectively.Through measuring the displacement between magenta and cyan dashed lines, it can be found that the rod with L/ax = 2.63±0.31manifests the discontinuous "hopping" motion punctuated by waiting periods; that is, the rod remains almost stationary for ~ 300s and then hops to a new position.However, the rod with L/ax = 2.97 ± 0.35 exhibits the continuous diffusion, yielding an exceptional situation where different dynamic behaviors appear.Strikingly, such a change of the dependence of diffusivity on rod length for the thick rod is distinct from the monotonic situation for the thin rod, as reported in previous works 20 .This can be confirmed through examining the trajectories obtain by Nearest Neighbor Search tracker 45 (Fig. S3).Fig. 1c presents the representative trajectories for rods with different lengths corresponding to the hop and continuous motions.Specifically, one can identify that the trajectories of the rods with L/ax = 2.02±0.18and L/ax = 2.97±0.35undergo remarkably faster diffusion than that with 2.51 ±0.27 where the "hop" can be definitely observed, indicating unconventional dynamics of the thick rods with lengths of integral multiple mesh size.
For a more quantitative analysis, we first examine the mean square displacement , where z(t) is the displacement of the center of mass along the direction of the rod contour that is determined by the eigendecomposition algorithm 46 , S6, respectively.For the noncommensurate rod, Gs(z, t) lines exhibit regular peaks, indicating that the rod transport undergoes hopping-like dynamics.In contrast, Gs(z, t) lines of the commensurate rod become shallow and irregular peaks.Such an exceptional dependence of the dynamics clarifies that the length commensuration is critical for the dynamical behaviors of thick rods, yielding new principles for the optimal design of rodlike particles with highly efficient transport in macromolecular networks.For instance, as the networks and surface chemistry of rods used in the experiments possess excellent biocompatibility, it can be anticipated that the findings provide useful guidelines for designing the shape of drug delivery systems 36,48 .

Detailed microscopic dynamics revealed by molecular simulations
To elucidate the physical origin of the unconventionally fast dynamics, we turn to the detailed microscopic dynamics of thick rods in macromolecular networks.We simulate the transport of a rodlike particle in a cross-linked network using dissipative particle dynamics (DPD) 49 .Full technical details on the simulation model are described in the Methods and Supplementary Information II and briefly introduced here.The configuration of the network is taken to be a hexa-functional network 50 , as illustrated in Fig. 2a.A general particlebuilding model 51 is adopted to build a set of rod particles fabricated by numbers of beads, which move as rigid bodies 52 .To bring out the entropic nature of the interplay between rods and network strands, the rod-strand interaction is set to be the same as that between like beads, capturing the physical nature of the interaction in this system where the free energy change is predominately contributed by the entropy 53 .For all systems, the average network mesh size is at around ax = 3.35rc, while various values of d and L are set to consider their roles in the rod dynamics.rc is the cutoff radius of the soft potentials and is used as the length unit of the system.Particularly, d is set to be comparable to the mesh size, ranging from 1.3ax to 1.9ax, and L ranges from 1.5ax to 4.4ax to avoid the effect of rod rotation.The normalized sizes d/ax and L/ax are used, representing the size matching between a rod and a network mesh.Typical MSDs for a set of L/ax at d/ax = 1.4 are shown in Fig. 2b.At short time scales, all rods move ballistically and 2 ( ) z t Δ scales as t 2 .At longer time scales, for rods with a length noncommensurate with L/ax = 1.5, 2.6 and 3.5, a plateau apparently emerges between the short-and long-time diffusion regimes, indicating that the rod transport undergoes slow hopping-like dynamics 10 .In striking contrast, at L/ax = 2.1, 3.1 and 4.1, they evolve directly from the ballistic regime to the Fickian regime with 2 ( ) z t t Δ demonstrating that fast transport dynamics occurs upon that L reaches around an integral multiple of ax.The diffusion coefficients corresponding to the MSDs are illustrated in Fig. 2c, which also demonstrates the nonmonotonic dependence of the diffusivity on L. This yields optimal values for the fast diffusion of such thick rods in macromolecular networks.
In order to provide a detailed insight into the fast transport mentioned above, we compare dynamical behaviors in the regimes of hopping and fast dynamics through calculating different parameters, i.e., typical trajectories in the direction of z-axis (Fig. 2d) and the corresponding Gs(z, t) (Fig. 2e and f).The trajectory of the rod in the regime of hopping dynamics shows that the rod undergoes constrained motion punctuated by significant large-scale full hops (the red line in Fig. 2d).This can be confirmed from the Gs(z, t) lines where regular peaks emerge and the distance between adjacent peaks is around 1.0ax, indicating that the rod hops from a network cell to its neighboring one (Fig. 2e).By contrast, in the regime of fast dynamics, the rod that experiences the fast transport dynamics does not exhibit hopping-type events as sharp as those observed for the noncommensurate rods.Instead, only random and local waiting intervals of small scales appear in the trajectory (the green line in Fig. 2d), resulting in shallow and irregular peaks in lines of Gs(z, t) (Fig. 2f).By comparing Fig. 2e, f to Fig. 1f, g respectively, a good agreement between simulation and experimental results is identified, corroborating the hopping and fast dynamics for the noncommensurate and commensurate rods.Clearly, the Gs(z, t) for the commensurate rod is non-Gaussian, but is different from those of hopping dynamics with regular peaks or Brownian dynamics with Gaussian distribution (Fig. S7).For the purpose of evaluating the non-Gaussianity, we further calculate the corresponding non-Gaussian parameter α1(t) = (1/3)⟨Δz 4 (t)⟩/⟨Δz 2 (t)⟩ 2 − 1 in the direction of z-axis as displayed in Fig. 2g.Indeed, the dynamics in the fast regime presents the weak non-Gaussianity intermediating between the strong dynamical heterogeneity of hopping diffusion and the Gaussian behavior of Brownian dynamics.
To evaluate the generality of such dynamic dependence on the length scales of rod and network, we systematically explore the diffusion behaviors of thick rods as a function of L/ax and d/ax, allowing us to construct a diagram of diffusion dynamics in the two-parameter space, as depicted in Fig. 2h.It can be found that for each size ratio interval L/ax∈[n, n + 1] with n = 0, 1, 2, . . ., the diffusion of rod possesses similar dynamic types.That is, three characteristic regimes can be discriminated as denoted by the colored circles in Fig. 2h: when L is noncommensurate with ax, the rod diffusion is featured by hopping between neighboring network cells for small d/ax while it turns to the trapped dynamics 10 characterized by slight fluctuation around its equilibrium position upon about d/ax>1.6;interestingly, between two neighboring hopping regimes there indeed exists the regime of faster longitudinal dynamics, where L reaches around an integral multiple of ax, corresponding to Fig. 2b.The result in the diagram clearly corroborates the nonmonotonic dependence of the diffusivity on L for thick rods as well as the regime of the unconventionally fast dynamics.The more random fluctuation in the displacement of the faster dynamics implies that the rod undergoes a lower free energy barrier, and the thermal noise as well as the strong effect of the local environment may thereby play a nontrivial role.

Entropic effect examined through evaluating diverse factors
In order to explore the physical mechanism underlying the nonmonotonic dependence of the diffusivity of such thick rods, we conduct comprehensive investigations into the influence of diverse factors attributed to the large diameters of rods.To this end, first, the diameter is reduced and thereby representative thin rods, with d/ax = 0.18, are considered.Fig. 2i illustrates the dependence of MSDs and diffusion coefficients on the length of the thin rods.The simulation results demonstrate that increasing the rod length monotonically reduces the diffusion coefficient, which is consistent with some previous studies 20 but in sharp contrast to the nonmonotonically decreased diffusivity with respect to the increase of rod length for thick rods (Fig. 2c).This implies that the entropic contribution, originated from the conformational penalty of loops, which consist of a group of end-to-end connected polymer strands 9 , deformed by the thick rods with large diameter, plays a key role in the free energy barrier.Basically, the number of loops around a thick rod is determined by the rod length, leading to the rod-length dependent free energy barrier for thick rods.
Next, to further examine the aforementioned entropic nature of the length-dependent, nonmonotonic diffusion dynamics of thick rods, we evaluate the influence of the enthalpy contributions through performing simulations by fixing d/ax and varying the particle-strand interaction parameter, arp, which measures the relative strength of the interaction between rod and network strand (see Supplementary Information II for more details).A larger arp corresponds to a stronger repulsion between these both species, while between like species it is chosen to be 25 49 .As shown in Figs.3a, c and S8, the faster longitudinal dynamics continues once the rod length reaches around an integral multiple of ax for both attractive (arp = 20) and repulsive (arp = 30) particle-strand interactions, resembling the above results for the systems governed by almost purely entropic effects.It underscores that the entropic contribution still dominates the behaviors over a wide range of the energetic interactions.
Last, we quantify the distribution of mesh sizes for some typical biological networks and thereby examine the entropic effect regarding the fast dynamics of thick rods in them, in view of the fact that there is a distribution of mesh sizes in a real network which causes the polydispersity of molecular structures of the networks.Specifically, the distribution of mesh sizes is quantified through calculating the coefficient of variation of mesh sizes, defined as CV = σax/ax, where σax is the standard deviation of mesh sizes.The CV of some typical biological mucus approximately ranges from 0.1 to 0.8 23,[55][56][57] , as listed in Table S4.
Thus, we perform simulations to study the thick rod diffusion in the macromolecular network with the similar distribution of mesh sizes.Based on the analysis of CV in mucus, we set CV = 0.3, 0.7 and 0.8 by varying the distribution of strand lengths in our simulations.
As shown in Figs.3b, d and S9, for d/ax=1.4,the fast dynamics indeed occurs when the commensurate rods at CV = 0.3 and 0.7.Even at CV = 0.8, the fast dynamics emerges for the thicker rods with d/ax =1.5, indicating that the thicker the rod, the wider pore size distribution permitted.Nevertheless, the simulation results of the macromolecular networks, whose polydispersity approximates to those of typical biological networks, can basically fall to the physical principles revealed based on the regular network, underscoring that the neat model can be applied to mimic the dynamical behaviors in mucus, especially for the thick rods concerned in our work.Moreover, this indicates again that the entropic effect induced by the thick rods accounts for the unconventional dynamic behaviors, as the influence of the mesh polydispersity is significantly impaired upon the remarkable structural deformation.

Physical origin established by theoretical analysis
To further confirm the entropic nature and thereby pinpoint the physical origin behind the unconventionally fast transport of a rod in a macromolecule network, we develop a theoretical model to analyze the free energy landscape and dynamical regimes depending on the length scales of rod and network.More details regarding this theoretical model can be found in Supplementary Information III.Considering a canonical ensemble of a rod in a macromolecule network, it is specified by (1) the set of crosslinks { } 1 , with M crosslinks between the efficiently bridged Gaussian chains; (2) the collection of linker connections marked as the tuple (i, j); (3) the continue curve path of linked strands Rij(s) . For a Gaussian chain of N bonds of Kuhn length b, the network mesh size 10 , ax = bN 1/2 = 1 is the unit length of the system.Through coupling the particle effect into the renowned theory of network elasticity 58,59 , the partition function of the rod-network system takes the form: ( ) where is the tangent direction of a network strand, rod r and rod l give the position and direction of the rod, δ is the delta function, β = 1 / kBT, is the Boltzmann constant, and T is temperature.Umr represents the hard-core monomer-rod interaction, which takes the form, / 2, / 2 ( , , ) 0 where ij ij rod ′ = − r r r .The free energy of the rod-network system has the form, ( , ) ln ( , ) To quantitatively examine the free energy experienced by the rod, the free energy change is defined as , where z  is the position of rod center as denoted by the coordinate in Fig. 4b, and ( ) min F z  represents the minimum free energy along z-axis.
We calculate some typical profiles of ΔF as a function of / x z a  for various d/ax at L/ax = 1.5 and 2.0, corresponding respectively to the noncommensurate and commensurate rods (Fig. 4a).As illustrated by the schematic in Fig. 4b, ΔF reaches the maximum at the network cell center with 0 z =  , and the minimum at the network cell edge with 0.5 x z a =  .Then, the height of the free energy barrier for the transition of the rod is  is the maximum of ( ) F z Δ  .In view of the Kramers' escape theory 60  there is a valley of 0 < Ub < kBT for commensurate rods, corresponding to the fast dynamics.
Moreover, the definite boundary of Ub = kBT is plotted, distinguishing the detailed ranges of fast dynamics around an integral multiple of L/ax.This demonstrates that the fast dynamics is attributable to the anomalous "faster-than-expected" behavior of the sliding dynamics while the length of a thick rod reaches around an integral multiple of the network mesh size, and the emergence of sliding dynamics can be ascribed to such low free energy barrier, with which thermal noise may remarkably thrust the dynamical process.Moreover, the diffusion coefficients obtained in the diffusive regime show nonmonotonic dependence on the rod length (Figs.1c and 2d), indicating the changes of the free energy landscape.
The fairly smooth landscape reveals that the fast dynamics really follows the sliding dynamics characterized by the energy barrier of order of kBT and emerging predominantly in the diffusion of proteins along DNA 61,62 .Our simulations demonstrate that in such unique dynamics, MSD is simply proportional to time (Fickian), yet the DPDF is not Gaussian as should be expected of a classical random walk, bearing great resemblance to the "anomalous yet Brownian" diffusion as found in some crowded fluids containing colloidal particles, macromolecules and filaments 63,64 .
Based on the calculations of the free energy landscape, we extend the theoretical analysis to microscopic dynamics to provide a refined picture of these dynamical regimes (see Supplementary Information IV for more details).Fundamentally, its microscopic dynamics can be described by the generalized Fokker-Planck equation 65 , where α is the anomalous coefficient.γ = 2πηL/ln(L/d) denotes the longitudinal friction coefficient of the rod, where η is the viscosity, and γκ = kBT.The distribution of the total number of network cells traversed can be obtained in Fourier-Laplace space from the Montroll-Weiss equation 66 , where S(k,s) is the Fourier-Laplace transform of , and derived by Mel'nikov 67,68 , where hop τ is the characteristic hopping time, ( , ) / 4 exp( / 4 ) ( ) where the exponential tails prevail.As the red curve in Fig. 4c, the plot of ( , ) s G z t  exhibits an identical shape as expected to the hopping dynamics shown in Figs 1f and 2e.
For low barriers of 0 < Ub < kBT, a master equation is applied to respect the nonlocal nature of irregular peaks in ( , ) where the kernel where I0 is the modified Bessel function of the first kind.Here we plot the blue curve of ( , ) in the inset of Fig. 4c, demonstrating the random distribution of peaks resembling the sliding dynamics shown in Figs.1g and 2f.
For no barrier of Ub = 0, the solution of eq. 2 takes the Gaussian form ( ) where τ0 is the characteristic waiting time in solvents.As shown by the cyan curve in Fig. 4c the theoretical result is consistent with the simulated curve at the late stage (Fig. S7), reflecting the Gaussian distribution of Brownian dynamics.

Discussions
Through conducting the experiments of single-particle tracking, performing molecular simulations and developing new theories, the transport dynamics of thick rods in macromolecular networks have been comprehensively investigated.We found that by tuning the rod length with respect to the averaged mesh size relatively fast longitudinal dynamics occurs once the rod length reaches around an integral multiple of the mesh size.
We identified that the unconventionally fast transport follows sliding dynamics, which is demonstrated to be anomalous yet, Brownian.We further gave the analytical expression of time-displacement distribution of sliding dynamics and clarified its physical relationship with hopping and Brownian dynamics.Our results revealed that the transition of these dynamical regimes is fundamentally attributed to the rod-length dependent free energy barrier originated predominantly from the entropic contribution due to the conformational penalty of strands deformed by thick rods, in contrast to the rigorous periodic potential.The findings might be of immediate interest to the optimal design of particle transport in various networks, biological or synthetic.As similar landscape of free energy is experienced by both passive and active rods in the confined environment of macromolecular networks, we speculate that the dependence of transport dynamics on size commensuration could be extended to the active case, suggesting the fundamental cornerstone for the further understanding of these nonequilibrium phenomena as well as the dissipative self-assembly of diverse building blocks 35,70 .hydrodynamic behavior emerges.The equations of motion are integrated in time with a modified velocity-Verlet algorithm.The factor kBT is taken as the characteristic energy scale.
In our simulations, kBT = 1.The characteristic time scale is then defined as τ = (mrc 2 /kBT) 1/2 = 1.The remaining simulation parameter are γ = 4.5 and Δt = 0.02τ with a total bead number density of ρ = 3.To demonstrate the dynamics of a rod-like particle in the network, we choose a cubic box with dimensions of 42.76rc × 42.76rc × 42.76rc.
Fig. 2h, allowing us to construct Ub landscape plotted by the colored contour map in the regimes are different, depending on Ub.Our following discussions are thereby based on the aforementioned regimes of Ub.For high barriers of Ub > kBT, distribution of the particle at stationary state, n is an integral number, and C is the normalization constant such that 0 can deduce that for large time scales, it takes the form,

Figure 1 |
Figure 1 | Single-particle tracking of rods in synthetic networks.a. Schematic illustration of dark-field microscopy (DFM) technique to track the nanorod diffusion.b.Representative time series of snapshots of nanorods diffused in PEGDA network, where L/ax: (I) 2.51 ± 0.27 and (II) 2.97 ± 0.35.The edge length of each snapshot is 0.5um.The cyan dashed lines mark the spontaneous position of the center of the nanorod, while the magenta dashed line denotes the initial position of the rod center (t = 0s).c.Typical trajectories of rods with different L/ax: 2.02 ± 0.18, 2.51 ± 0.27 and 2.97 ± 0.35 from top to the bottom, corresponding respectively to the fast, hopping and fast dynamics.The color and scale bars on the left top indicate the values of the temporal and spatial scales of the trajectories.d. ⟨∆z 2 (t)⟩ for different L/ax in experiments.e. D/D0 as a function of L/ax.The error bars represent the standard deviation.Experimental Gs(z, t) of rods with f.L/ax = 2.51 ± 0.27 and g.L/ax = 2.97 ± 0.35, where solid lines represent smoothed results of circles.

Figure 2 |
Figure 2 | Detailed microscopic dynamics revealed by molecular simulations.a. Schematic of a rodlike particle of diameter d and length L in a macromolecular network with mesh size ax.b. ⟨∆z 2 (t)⟩ for different L/ax at d/ax = 1.4.c.D for different L/ax at d/ax = 1.4.d.Typical trajectories of the center of mass of the corresponding rods with L/ax = 2.6 (red) and L/ax = 3.1 (green) at d/ax = 1.4 along z-axis, where thin lines represent unsmoothed original displacement while thick lines represent smoothed displacement determined by a wavelet-based method 54 .Simulated Gs(z, t) of rods with (e) L/ax = 2.6 and (f) L/ax = 3.1 at d/ax = 1.4 along z-axis, showing hopping-type diffusion and speeding-up diffusion respectively.g.  for Brownian dynamics (cyan), fast dynamics (blue) and hopping dynamics (red).h.Diagram of rod dynamics interrelating to L/ax and d/ax.Circles: simulation results of (red) hopping, (cyan) fast, and (purple) trapped diffusion beheviors.The contour map: theoretical results of Ub with various d/ax and L/ax.The green lines:

Figure 4 |
Figure 4 | Physical origin established by theoretical analysis.a. ∆F as a function of ̃ for various d/ax, where L/ax =1.5 (left) and 2.0 (right).The color bar indicates the value of d/ax.Typical positions of rod center for various L/ax are schematically illustrated in b. A. / x z a  = 0, L/ax = 1.5;B. / x z a  = 0, L/ax = 2.0; C. / x z a  = 0.5, L/ax = 2.0.c.Theoretical results of Gs(, t) for (red) hopping, (blue) sliding and (cyan) Brownian dynamics.