Localization transition of the three-dimensional lorentz model and continuum percolation

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.

Percolation through voids around toroidal inclusions

In the case of media comprised of impermeable particles, fluid flows through voids around impenetrable grains. For sufficiently low concentrations of the latter, spaces around grains join to allow transport on macroscopic scales, whereas greater impenetrable inclusion densities disrupt void networks and block macroscopic fluid flow. A critical grain concentration ρ_{c} marks the percolation transition or phase boundary separating these two regimes. With a dynamical infiltration technique in which virtual tracer particles explore void spaces, we calculate critical grain concentrations for randomly placed interpenetrating impermeable toroidal inclusions; the latter consist of surfaces of revolution with circular and square cross sections. In this manner, we study continuum percolation transitions involving nonconvex grains. As the radius of revolution increases relative to the length scale of the torus cross section, the tori develop a central hole, a topological transition accompanied by a cusp in the critical porosity fraction for percolation. With a further increase in the radius of revolution, as constituent grains become more ringlike in appearance, we find that the critical porosity fraction converges to that of high-aspect-ratio cylindrical counterparts only for randomly oriented grains.

First-passage times in conical varying-width channels biased by a transverse gravitational force: Comparison of analytical and numerical results

We study the crossing time statistic of diffusing point particles between the two ends of expanding and narrowing two-dimensional conical channels under a transverse external gravitational field. The theoretical expression for the mean first-passage time for such a system is derived under the assumption that the axial diffusion in a two-dimensional channel of smoothly varying geometry can be approximately described as a one-dimensional diffusion in an entropic potential with position-dependent effective diffusivity in terms of the modified Fick-Jacobs equation. We analyze the channel crossing dynamics in terms of the mean first-passage time, combining our analytical results with extensive two-dimensional Brownian dynamics simulations, allowing us to find the range of applicability of the one-dimensional approximation. We find that the effective particle diffusivity decreases with increasing amplitude of the external potential. Remarkably, the mean first-passage time for crossing the channel is shown to assume a minimum at finite values of the potential amplitude.

Heterogeneous Nanostructures Cause Anomalous Diffusion in Lipid Monolayers

The diffusion and mobility in biomembranes are crucial for various cell functions; however, the mechanisms involved in such processes remain ambiguous due to the complex membrane structures. Herein, we investigate how the heterogeneous nanostructures cause anomalous diffusion in dipalmitoylphosphatidylcholine (DPPC) monolayers. By identifying the existence of condensed nanodomains and clarifying their impact, our findings renew the understanding of the hydrodynamic description and the statistical feature of the diffusion in the monolayers. We find a universal characteristic of the multistage mean square displacement (MSD) with an intermediate crossover, signifying two membrane viscosities at different scales: the short-time scale describes the local fluidity and is independent of the nominal DPPC density, and the long-time scale represents the global continuous phase taking into account nanodomains and increases with DPPC density. The constant short-time viscosity reflects a dynamic equilibrium between the continuous fluid phase and the condensed nanodomains in the molecular scale. Notably, we observe an "anomalous yet Brownian" phenomenon exhibiting an unusual double-peaked displacement probability distribution (DPD), which is attributed to the net dipolar repulsive force from the heterogeneous nanodomains around the microdomains. The findings provide physical insights into the transport of membrane inclusions that underpin various biological functions and drug deliveries.

Normal and anomalous diffusion in the disordered wind-tree model

Ehrenfests' wind-tree model (EWTM) refers to a two-dimensional system where noninteracting point tracer particles move through a random arrangement of overlapping or nonoverlapping square-shaped scatterers. Here, extensive event-driven molecular dynamics simulations of the EWTM at different reduced scatterer densities ρ are presented. For nonoverlapping scatterers, the asymptotic motion of the tracer particles is diffusive. We compare their diffusion coefficient D, as obtained from the simulation, with that predicted by kinetic theory where D^{-1} is expanded up to the second order in the scatterer density. While at low density quantitative agreement between theory and simulation is found, we show that beyond the low-density regime deviations to the theory are associated with the emergence of a maximum in the non-Gaussian parameter at intermediate times. For the case of overlapping scatterers, in agreement with a theoretical prediction, the asymptotic motion of the tracer particles is subdiffusive, i.e., the mean-squared displacement at long times t grows like t^{1-2ρ/3}. We propose a model of the van Hove correlation function that describes the density dependence of the tracer particles' asymptotic subdiffusive transport on a quantitative level.

Local Dynamical Heterogeneity in Simple Glass Formers

We study the local dynamical fluctuations in glass-forming models of particles embedded in d-dimensional space, in the mean-field limit of d→∞. Our analytical calculation reveals that single-particle observables, such as squared particle displacements, display divergent fluctuations around the dynamical (or mode-coupling) transition, due to the emergence of nontrivial correlations between displacements along different directions. This effect notably gives rise to a divergent non-Gaussian parameter, α_{2}. The d→∞ local dynamics therefore becomes quite rich upon approaching the glass transition. The finite-d remnant of this phenomenon further provides a long sought-after, first-principle explanation for the growth of α_{2} around the glass transition that is not based on multiparticle correlations.

A simulation study on the subdiffusion of polymer chains in crowded environments containing nanoparticles

Polymer chains in crowded environments often show subdiffusive behavior. We adopt molecular dynamics simulations to study the conditions for the subdiffusion of polymer chains in crowded environments containing randomly distributed, immobile, attractive nanoparticles (NPs). The attraction is strong enough to adsorb polymer chains on NPs. The results show that subdiffusion occurs at a low concentration of polymer chains (c_p). A transition from subdiffusion to normal diffusion is observed when c_p exceeds the transition concentration , which increases with increasing concentration of NPs while decreases with increasing size of NPs. The high concentration and small size of NPs exert a big effect on the subdiffusion of polymer chains. The subdiffusive behavior of polymer chains can be attributed to the strong adsorption of polymer chains on the attractive NPs. For the subdiffusion case, polymer chains are adsorbed strongly on multiple NPs, and they diffuse via the NP-exchange diffusion mechanism. However for the normal diffusion case, polymer chains are either free or weakly adsorbed on one or a few NPs, and they diffuse mainly via the adsorption-and-desorption diffusion mechanism.

High-dimensional percolation criticality and hints of mean-field-like caging of the random Lorentz gas

The random Lorentz gas (RLG) is a minimal model for transport in disordered media. Despite the broad relevance of the model, theoretical grasp over its properties remains weak. For instance, the scaling with dimension d of its localization transition at the void percolation threshold is not well controlled analytically nor computationally. A recent study [Biroli et al., Phys. Rev. E 103, L030104 (2021)2470-004510.1103/PhysRevE.103.L030104] of the caging behavior of the RLG motivated by the mean-field theory of glasses has uncovered physical inconsistencies in that scaling that heighten the need for guidance. Here we first extend analytical expectations for asymptotic high-d bounds on the void percolation threshold and then computationally evaluate both the threshold and its criticality in various d. In high-d systems, we observe that the standard percolation physics is complemented by a dynamical slowdown of the tracer dynamics reminiscent of mean-field caging. A simple modification of the RLG is found to bring the interplay between percolation and mean-field-like caging down to d=3.

Critical Indices and Self-Similar Power Transform

“Odd” factor approximants of the special form suggested by Gluzman and Yukalov (J. Math. Chem. 2006, 39, 47) are amenable to optimization by power transformation and can be successfully applied to critical phenomena. The approach is based on the idea that the critical index by itself should be optimized through the parameters of power transform to be calculated from the minimal sensitivity (derivative) optimization condition. The critical index is a product of the algebraic self-similar renormalization which contributes to the expressions the set of control parameters typical to the algebraic self-similar renormalization, and of the power transform which corrects them even further. The parameter of power transformation is, in a nutshell, the multiplier connecting the critical exponent and the correction-to-scaling exponent. We mostly study the minimal model of critical phenomena based on expansions with only two coefficients and critical points. The optimization appears to bring quite accurate, uniquely defined results given by simple formulas. Many important cases of critical phenomena are covered by the simple formula. For the longer series, the optimization condition possesses multiple solutions, and additional constraints should be applied. In particular, we constrain the sought solution by requiring it to be the best in prediction of the coefficients not employed in its construction. In principle, the error/measure of such prediction can be optimized by itself, with respect to the parameter of power transform. Methods of calculation based on optimized power-transformed factors are applied and results presented for critical indices of several key models of conductivity and viscosity of random media, swelling of polymers, permeability in two-dimensional channels. Several quantum mechanical problems are discussed as well.

Mean-Field Caging in a Random Lorentz Gas

The random Lorentz gas (RLG) is a minimal model of both percolation and glassiness, which leads to a paradox in the infinite-dimensional, d → ∞ limit: the localization transition is then expected to be continuous for the former and discontinuous for the latter. As a putative resolution, we have recently suggested that, as d increases, the behavior of the RLG converges to the glassy description and that percolation physics is recovered thanks to finite-d perturbative and nonperturbative (instantonic) corrections [Biroli et al. Phys. Rev. E 2021, 103, L030104]. Here, we expand on the d → ∞ physics by considering a simpler static solution as well as the dynamical solution of the RLG. Comparing the 1/d correction of this solution with numerical results reveals that even perturbative corrections fall out of reach of existing theoretical descriptions. Comparing the dynamical solution with the mode-coupling theory (MCT) results further reveals that, although key quantitative features of MCT are far off the mark, it does properly capture the discontinuous nature of the d → ∞ RLG. These insights help chart a path toward a complete description of finite-dimensional glasses.

Optimized Factor Approximants and Critical Index

Based on expansions with only two coefficients and known critical points, we consider a minimal model of critical phenomena. The method of analysis is both based on and inspired with the symmetry properties of functional self-similarity relation between the consecutive functional approximations. Factor approximants are applied together with various natural optimization conditions of non-perturbative nature. The role of control parameter is played by the critical index by itself. The minimal derivative condition imposed on critical amplitude appears to bring the most reasonable, uniquely defined results. The minimal difference condition also imposed on amplitudes produces upper and lower bound on the critical index. While one of the bounds is close to the result from the minimal difference condition, the second bound is determined by the non-optimized factor approximant. One would expect that for the minimal derivative condition to work well, the bounds determined by the minimal difference condition should be not too wide. In this sense the technique of optimization presented above is self-consistent, since it automatically supplies the solution and the bounds. In the case of effective viscosity of passive suspensions the bounds could be found that are too wide to make any sense from either of the solutions. Other optimization conditions imposed on the factor approximants, lead to better estimates for the critical index for the effective viscosity. The optimization is based on equating two explicit expressions following from two different definitions of the critical index, while optimization parameter is introduced as the trial third-order coefficient in the expansion.

Interplay between percolation and glassiness in the random Lorentz gas

The random Lorentz gas (RLG) is a minimal model of transport in heterogeneous media that exhibits a continuous localization transition controlled by void space percolation. The RLG also provides a toy model of particle caging, which is known to be relevant for describing the discontinuous dynamical transition of glasses. In order to clarify the interplay between the seemingly incompatible percolation and caging descriptions of the RLG, we consider its exact mean-field solution in the infinite-dimensional d→∞ limit and perform numerics in d=2...20. We find that for sufficiently high d the mean-field caging transition precedes and prevents the percolation transition, which only happens on timescales diverging with d. We further show that activated processes related to rare cage escapes destroy the glass transition in finite dimensions, leading to a rich interplay between glassiness and percolation physics. This advance suggests that the RLG can be used as a toy model to develop a first-principle description of particle hopping in structural glasses.

Softness, anomalous dynamics, and fractal-like energy landscape in model cell tissues

Epithelial cell tissues have a slow relaxation dynamics resembling that of supercooled liquids. Yet, they also have distinguishing features. These include an extended short-time subdiffusive transient, as observed in some experiments and recent studies of model systems, and a sub-Arrhenius dependence of the relaxation time on temperature, as reported in numerical studies. Here we demonstrate that the anomalous glassy dynamics of epithelial tissues originates from the emergence of a fractal-like energy landscape, particles becoming virtually free to diffuse in specific phase space directions up to a small distance. Furthermore, we clarify that the stiffness of the cells tunes this anomalous behavior, tissues of stiff cells having conventional glassy relaxation dynamics.

Nanoparticle Tracking to Probe Transport in Porous Media

From the granular and fractured subsurface environment to highly engineered polymer membranes used in pharmaceutical purification, porous materials are ubiquitous in nature and industrial applications. In particular, porous media are used extensively in processes including water treatment, pharmaceutical sterilization, food/beverage processing, and heterogeneous catalysis, where hindered mass transport is either essential to the process or a necessary but undesirable limitation. Unfortunately, there are currently no universal models capable of predicting mass transport based on a description of the porous material because real porous materials are complex and because many coupled dynamic mechanisms (e.g., adsorption, steric effects, hydrodynamic effects, electrostatic interactions, etc.) give rise to the observed macroscopic transport phenomena.While classical techniques, like nuclear magnetic resonance and dynamic light scattering, provide useful information about mass transport in porous media at the ensemble level, they provide limited insight into the microscopic mechanisms that give rise to complex phenomena such as anomalous diffusion, hindered pore-space accessibility, and unexpected retention under flow, among many others. To address this issue, we have developed refractive index matching imaging systems, combined with single-particle tracking methods, allowing the direct visualization of single-particle motion within a variety of porous materials.In this Account, we summarize our recent efforts to advance the understanding of nanoparticle transport in porous media using single-particle tracking methods in both fundamental and applied scenarios. First, we describe the basic principles for two-dimensional and three-dimensional single-particle tracking in porous materials. Then, we provide concrete examples of nanoparticle transport in porous materials from two perspectives: (1) understanding fundamental elementary particle transport processes in porous media, including pore accessibility and cavity escape, which limit transport in porous media, and (2) facilitating applications in industrial processes, e.g., by understanding the mechanisms of particle fouling and remobilization in filtration membranes. Finally, we provide an outlook of opportunities associated with investigating other types of mass transport in confined environments using single-particle tracking methods, including electrophoretic and self-propelled motion.

Diffusion of Anisotropic Colloids in Periodic Arrays of Obstacles

Colloidal suspensions in confined geometries exhibit rich diffusion dynamics governed by particle shapes and particle-confinement interactions. Here, we propose a colloidal system, consisting of ellipsoids in periodic array of obstacles, to investigate the confined diffusion of anisotropic colloids. From the obstacle density-dependent diffusion, we discover a decoupling of translational and rotational diffusion in which only rotational motion is localized while translational motion remains diffusive. Moreover, by evaluating the probability distributions of displacements, we found Brownian but non-Gaussian diffusion behaviors with increasing the obstacle densities, which originates from the shape anisotropy of the colloid and the multiplicity of the local configurations of the ellipsoids with respect to the obstacle. Our results suggest that the shape anisotropy and spatial confinements play a vital role in the diffusion dynamics. It is important for understanding the transportations of anisotropic objects in complex environments.

Molecular transport in systems containing binding obstacles

We studied the movement of particles in crowded environments by means of extensive Monte Carlo simulations. The dynamic lattice liquid model was employed for this purpose. It is based on the cooperative movement concept and allows the study of systems at high densities. The cooperative model of molecular transport is assumed: the motion of all moving particles is highly correlated. The model is supposed to mimic lateral motion in a membrane and therefore the system is two-dimensional with moving objects and traps placed on a triangular lattice. In our study the interaction (binding) of traps with moving particles was assumed. The conditions in which subdiffusive motion appeared in the system were analysed. The influence of the strength of binding on the dynamic percolation threshold was also shown. The differences in the dynamics compared to systems with impenetrable obstacles and with systems without correlation in motion were presented and discussed. It was shown that in the case of correlated motion the influence of deep traps is similar to that of impenetrable obstacles. If the traps are shallow a recovery to normal diffusion was observed for longer time periods.

Understanding biochemical processes in the presence of sub-diffusive behavior of biomolecules in solution and living cells

In order to maintain cellular function, biomolecules like protein, DNA, and RNAs have to diffuse to the target spaces within the cell. Changes in the cytosolic microenvironment or in the nucleus during the fulfillment of these cellular processes affect their mobility, folding, and stability thereby impacting the transient or stable interactions with their adjacent neighbors in the organized and dynamic cellular interior. Using classical Brownian motion to elucidate the diffusion behavior of these biomolecules is hard considering their complex nature. The understanding of biomolecular diffusion inside cells still remains elusive due to the lack of a proper model that can be extrapolated to these cases. In this review, we have comprehensively addressed the progresses in this field, laying emphasis on the different aspects of anomalous diffusion in the different biochemical reactions in cell interior. These experiment-based models help to explain the diffusion behavior of biomolecules in the cytosolic and nuclear microenvironment. Moreover, since understanding of biochemical reactions within living cellular system is our main focus, we coupled the experimental observations with the concept of sub-diffusion from in vitro to in vivo condition. We believe that the pairing between the understanding of complex behavior and structure-function paradigm of biological molecules would take us forward by one step in order to solve the puzzle around diseases caused by cellular dysfunction.

Dynamically asymmetric binary glass formers studied by dielectric and NMR spectroscopy

We investigate the component dynamics in asymmetric binary glass formers. Focusing on the dielectric spectra of the high-T_g components m-tricresyl phosphate and quinaldine mixed with toluene as low-T_g component, the broadend spectra cannot be described by Kohlrausch or Cole-Davidson (CD) functions. Instead, we apply a generalized CD function which allows to control the width of the susceptibility independently of its high-frequency flank. The spectra show a common broadening and failure of the frequency-temperature superposition with increasing toluene concentration. This is confirmed by stimulated echo experiments showing an increased stretching of the probed orientational correlation function. In analogy to the definition of T_g, we consider "isodynamic points". For each component, a different but linear concentration dependence of 1/T_iso is revealed, indicating different time scales. Qualitativly, we do not find significant differences for the present mixtures with T_g-contrasts of 63-89K compared to those with larger T_g-contrast ( [Formula: see text] K): Whereas the high-T_g component shows relaxation features similar to those of neat glass formers, yet, with "atypical" weak relaxation broadening, the faster low-T_g component displays pronounced dynamic heterogeneities. This is supported by scrutinizing NMR relaxation data of several mixtures investigated previously as a function of concentration. A universal evolution of the dynamics of the high-T_g as well as the low-T_g component is suggested for mixtures with high [Formula: see text]T_g .

Anomalous transport in the soft-sphere Lorentz model

The sensitivity of anomalous transport in crowded media to the form of the inter-particle interactions is investigated through computer simulations. We extend the highly simplified Lorentz model towards realistic natural systems by modeling the interactions between the tracer and the obstacles with a smooth potential. We find that the anomalous transport at the critical point happens to be governed by the same universal exponent as for hard exclusion interactions, although the mechanism of how narrow channels are probed is rather different. The scaling behavior of simulations close to the critical point confirm this exponent. Our result indicates that the simple Lorentz model may be applicable to describing the fundamental properties of long-range transport in real crowded environments.

Anomalous Diffusion in Random Dynamical Systems

Consider a chaotic dynamical system generating diffusionlike Brownian motion. Consider a second, nonchaotic system in which all particles localize. Let a particle experience a random combination of both systems by sampling between them in time. What type of diffusion is exhibited by this random dynamical system? We show that the resulting dynamics can generate anomalous diffusion, where in contrast to Brownian normal diffusion the mean square displacement of an ensemble of particles increases nonlinearly in time. Randomly mixing simple deterministic walks on the line, we find anomalous dynamics characterized by aging, weak ergodicity breaking, breaking of self-averaging, and infinite invariant densities. This result holds for general types of noise and for perturbing nonlinear dynamics in bifurcation scenarios.

Binary colloidal glasses: linear viscoelasticity and its link to the microscopic structure and dynamics

We study the relation between the microscopic structure and dynamics and the macroscopic rheological response of glass-forming colloidal suspensions, namely binary colloidal hard-sphere mixtures with large size asymmetry (1 : 5) that span a large range of mixture compositions close to the glass transition. The dynamical shear moduli are measured by oscillatory rheology and the structure and dynamics on the single-particle level by confocal microscopy. The data are compared with Brownian Dynamics simulations and predictions from mode-coupling theory based on the Percus-Yevick approximation. Experiments, simulations and theory consistently observe a strong decrease of the intermediate-frequency mechanical moduli combined with faster dynamics at intermediate mixing ratios and hence a non-monotonic dependence of these parameters but a localization of the large particles which decreases monotonically as the fraction of small particles is increased. We find that the Generalized-Stokes Einstein relation applied to the mean square displacements of the two components leads to a reasonable estimate of the shear moduli of the mixtures and hence links the rheological response to the particle dynamics which in turn reflects the microscopic structure.

Non-universality of the dynamic exponent in two-dimensional random media

The diffusion of solutes in two-dimensional random media is important in diverse physical situations including the dynamics of proteins in crowded cell membranes and the adsorption on nano-structured substrates. It has generally been thought that the diffusion constant, D, should display universal behavior near the percolation threshold, i.e., D ~ (ϕ − ϕ_c)^μ, where ϕ is the area fraction of the matrix, ϕ_c is the value of ϕ at the percolation threshold, and μ is the dynamic exponent. The universality of μ is important because it implies that very different processes, such as protein diffusion in membranes and the electrical conductivity in two-dimensional networks, obey similar underlying physical principles. In this work we demonstrate, using computer simulations on a model system, that the exponent μ is not universal, but depends on the microscopic nature of the dynamics. We consider a hard disc that moves via random walk in a matrix of fixed hard discs and show that μ depends on the maximum possible displacement Δ of the mobile hard disc, ranging from 1.31 at Δ ≤ 0.1 to 2.06 for relatively large values of Δ. We also show that this behavior arises from a power-law singularity in the distribution of transition rates due to a failure of the local equilibrium approximation. The non-universal value of μ obeys the prediction of the renormalization group theory. Our simulations do not, however, exclude the possibility that the non-universal values of μ might be a crossover between two different limiting values at very large and small values of Δ. The results allow one to rationalize experiments on diffusion in two-dimensional systems.

Ideal circle microswimmers in crowded media

Microswimmers are exposed in nature to crowded environments and their transport properties depend in a subtle way on the interaction with obstacles. Here, we investigate a model for a single ideal circle swimmer exploring a two-dimensional disordered array of impenetrable obstacles. The microswimmer moves on circular orbits in the freely accessible space and follows the surface of an obstacle for a certain time upon collision. Depending on the obstacle density and the radius of the circular orbits, the microswimmer displays either long-range transport or is localized in a finite region. We show that there are transitions from two localized states to a diffusive state each driven by an underlying static percolation transition. We determine the non-equilibrium state diagram and calculate the mean-square displacements and diffusivities by computer simulations. Close to the transition lines transport becomes subdiffusive which is rationalized as a dynamic critical phenomenon.

A general phenomenological relation for the subdiffusive exponent of anomalous diffusion in disordered media

This work numerically investigates the diffusion of finite inert tracer particles in different types of fixed gels. The mean square displacement (MSD) of the tracers reveals a transition to subdiffusive motion MSD ∼ tα as soon as the accessible volume fraction p in the gel decreases from unity. Individual tracer dynamics reveals two types of particles in the gels: mobile tracers cross the system through percolating pores following subdiffusive dynamics MSDmob ∼ tαmob, while a fraction ptrap(p) of the particles remain trapped in finite pores. Below the void percolation threshold p < pc all the particles get trapped and α → 0. By separately studying both populations we find a simple phenomenological law for the mobile tracers αmob(p) ≈ a ln p + c where c ≈ 1 and a ∼ 0.2 depends on the gel type. On the other hand, a cluster-analysis of the gel accessible volume reveals a power law for the trapping probability ptrap ∼ (p/pc)-γ, with γ ≃ 2.9. This yields a prediction for the ensemble averaged subdiffusion exponent α = αmob(1 - ptrap). Our predictions are successfully validated against the different gels studied here and against numerical and experimental results in the literature (silica gels, polyacrylamide gels, flexible F-actin networks and in different random obstacles). Notably, the parameter a ∼ 0.2 presents small differences amongst all these cases, indicating the robustness of the proposed relation.

Percolation through Voids around Randomly Oriented Polyhedra and Axially Symmetric Grains

Porous materials made up of impermeable grains constrain fluid flow to voids around the impenetrable inclusions. A percolation transition marks the boundary between densities of grains permitting bulk transport and concentrations blocking traversal on macroscopic scales. With dynamical infiltration of void spaces using virtual tracer particles, we treat inclusion geometries exactly. We calculate the critical number density per volume ρ_{c} for a variety of axially symmetric shapes and faceted solids with the former including cylinders, ellipsoids, cones, and tablet shaped grains from highly oblate (platelike) to highly prolate (needlelike) extremes. For the latter, results suggest a common asymptotic value identical to the counterpart for aligned cylindrical grains. We find percolation thresholds for each of the five platonic solids (i.e., tetrahedra, cubes, octahedra, dodecahedra, and icosahedra) as well as truncated icosahedra. For each polyhedron type, we consider aligned and randomly oriented grains, finding distinct percolation thresholds for the former versus the latter only for cubes. The anomalous diffusion exponents we find differ from those of the universality class for discrete models on 3D lattices.

Generalized Langevin equation formulation for anomalous diffusion in the Ising model at the critical temperature

We consider the two- (2D) and three-dimensional (3D) Ising models on a square lattice at the critical temperature T_{c}, under Monte Carlo spin flip dynamics. The bulk magnetization and the magnetization of a tagged line in the 2D Ising model, and the bulk magnetization and the magnetization of a tagged plane in the 3D Ising model, exhibit anomalous diffusion. Specifically, their mean-square displacements increase as power laws in time, collectively denoted as ∼t^{c}, where c is the anomalous exponent. We argue that the anomalous diffusion in all these quantities for the Ising model stems from time-dependent restoring forces, decaying as power laws in time-also with exponent c -in striking similarity to anomalous diffusion in polymeric systems. Prompted by our previous work that has established a memory-kernel based generalized Langevin equation (GLE) formulation for polymeric systems, we show that a closely analogous GLE formulation holds for the Ising model as well. We obtain the memory kernels from spin-spin correlation functions, and the formulation allows us to consistently explain anomalous diffusion as well as anomalous response of the Ising model to an externally applied magnetic field in a consistent manner.

Nonmonotonic Classical Magnetoconductivity of a Two-Dimensional Electron Gas in a Disordered Array of Obstacles

Magnetotransport measurements in combination with molecular dynamics simulations on two-dimensional disordered Lorentz gases in the classical regime are reported. In quantitative agreement between experiment and simulation, the magnetoconductivity displays a pronounced peak as a function of the perpendicular magnetic field B which cannot be explained by existing kinetic theories. This peak is linked to the onset of a directed motion of the electrons along the contour of the disordered obstacle matrix when the cyclotron radius becomes smaller than the size of the obstacles. This directed motion leads to transient superdiffusive motion and strong scaling corrections in the vicinity of the insulator-to-conductor transitions of the Lorentz gas.

Diffusion, subdiffusion, and localization of active colloids in random post lattices

Combining experiments and theory, we address the dynamics of self-propelled particles in crowded environments. We first demonstrate that motile colloids cruising at constant speed through random lattices undergo a smooth transition from diffusive to subdiffusive to localized dynamics upon increasing the obstacle density. We then elucidate the nature of these transitions by performing extensive simulations constructed from a detailed analysis of the colloid-obstacle interactions. We evidence that repulsion at a distance and hard-core interactions both contribute to slowing down the long-time diffusion of the colloids. In contrast, the localization transition stems solely from excluded-volume interactions and occurs at the void-percolation threshold. Within this critical scenario, equivalent to that of the random Lorentz gas, genuine asymptotic subdiffusion is found only at the critical density where the motile particles explore a fractal maze.

Nonscaling displacement distributions as may be seen in fluorescence correlation spectroscopy

A continuous time random walk (CTRW) model with waiting times following the Lévy-stable distribution with exponential cutoff in equilibrium is a simple theoretical model giving rise to normal, yet non-Gaussian, diffusion. The distribution of the particles' displacements is explicitly time dependent and does not scale. Since fluorescent correlation spectroscopy (FCS) is often used to investigate diffusion processes, we discuss the influence of this lack of scaling on the possible outcome of the FCS measurements and calculate the FCS autocorrelation curves for such equilibrated CTRWs. The results show that although the deviations from Gaussian behavior may be detected when analyzing the short- and long-time asymptotic behavior of the corresponding curves, their bodies are still perfectly fitted by the fit forms used for normal diffusion. The diffusion coefficients obtained from the fits may however differ considerably from the true tracer diffusion coefficients as describing the time dependence of the mean squared displacement.

Dynamic heterogeneities and non-Gaussian behavior in two-dimensional randomly confined colloidal fluids

A binary mixture of superparamagnetic colloidal particles is confined between glass plates such that the large particles become fixed and provide a two-dimensional disordered matrix for the still mobile small particles, which form a fluid. By varying fluid and matrix area fractions and tuning the interactions between the superparamagnetic particles via an external magnetic field, different regions of the state diagram are explored. The mobile particles exhibit delocalized dynamics at small matrix area fractions and localized motion at high matrix area fractions, and the localization transition is rounded by the soft interactions [T. O. E. Skinner et al., Phys. Rev. Lett. 111, 128301 (2013)PRLTAO0031-900710.1103/PhysRevLett.111.128301]. Expanding on previous work, we find the dynamics of the tracers to be strongly heterogeneous and show that molecular dynamics simulations of an ideal gas confined in a fixed matrix exhibit similar behavior. The simulations show how these soft interactions make the dynamics more heterogeneous compared to the disordered Lorentz gas and lead to strong non-Gaussian fluctuations.

Nonergodic subdiffusion from transient interactions with heterogeneous partners

Spatiotemporal disorder has been recently associated to the occurrence of anomalous nonergodic diffusion of molecular components in biological systems, but the underlying microscopic mechanism is still unclear. We introduce a model in which a particle performs continuous Brownian motion with changes of diffusion coefficients induced by transient molecular interactions with diffusive binding partners. In spite of the exponential distribution of waiting times, the model shows subdiffusion and nonergodicity similar to the heavy-tailed continuous time random walk. The dependence of these properties on the density of binding partners is analyzed and discussed. Our work provides an experimentally testable microscopic model to investigate the nature of nonergodicity in disordered media.

Comparison of different models of motion in a crowded environment: a Monte Carlo study

In this paper we investigate the motion of molecules in crowded environments for two dramatically different types of molecular transport. The first type is realized by the dynamic lattice liquid model, which is based on a cooperative movement concept and thus, the motion of molecules is highly correlated. The second one corresponds to a so-called motion of a single agent where the motion of molecules is considered as a random walk without any correlation with other moving elements. The crowded environments are modeled as a two-dimensional triangular lattice with fixed impenetrable obstacles. Our simulation results indicate that the type of transport has an impact on the dynamics of the system, the percolation threshold, critical exponents, and on molecules' trajectories.

Molecular finite-size effects in stochastic models of equilibrium chemical systems

The reaction-diffusion master equation (RDME) is a standard modelling approach for understanding stochastic and spatial chemical kinetics. An inherent assumption is that molecules are point-like. Here, we introduce the excluded volume reaction-diffusion master equation (vRDME) which takes into account volume exclusion effects on stochastic kinetics due to a finite molecular radius. We obtain an exact closed form solution of the RDME and of the vRDME for a general chemical system in equilibrium conditions. The difference between the two solutions increases with the ratio of molecular diameter to the compartment length scale. We show that an increase in the fraction of excluded space can (i) lead to deviations from the classical inverse square root law for the noise-strength, (ii) flip the skewness of the probability distribution from right to left-skewed, (iii) shift the equilibrium of bimolecular reactions so that more product molecules are formed, and (iv) strongly modulate the Fano factors and coefficients of variation. These volume exclusion effects are found to be particularly pronounced for chemical species not involved in chemical conservation laws. Finally, we show that statistics obtained using the vRDME are in good agreement with those obtained from Brownian dynamics with excluded volume interactions.

Simulation of Molecular Transport in Systems Containing Mobile Obstacles

In this paper, we investigate the movement of molecules in crowded environments with obstacles undergoing Brownian motion by means of extensive Monte Carlo simulations. Our investigations were performed using the dynamic lattice liquid model, which was based on the cooperative movement concept and allowed to mimic systems at high densities where the motion of all elements (obstacles as well as moving particles) were highly correlated. The crowded environments are modeled on a two-dimensional triangular lattice containing obstacles (particles whose mobility was significantly reduced) moving by a Brownian motion. The subdiffusive motion of both elements in the system was analyzed. It was shown that the percolation transition does not exist in such systems in spite of the cooperative character of the particles' motion. The reduction of the obstacle mobility leads to the longer caging of liquid particles by mobile obstacles.

Time- and ensemble-averages in evolving systems: the case of Brownian particles in random potentials

Anomalous diffusion is a ubiquitous phenomenon in complex systems. It is often quantified using time- and ensemble-averages to improve statistics, although time averages represent a non-local measure in time and hence can be difficult to interpret. We present a detailed analysis of the influence of time- and ensemble-averages on dynamical quantities by investigating Brownian particles in a rough potential energy landscape (PEL). Initially, the particle ensemble is randomly distributed, but the occupancy of energy values evolves towards the equilibrium distribution. This relaxation manifests itself in the time evolution of time- and ensemble-averaged dynamical measures. We use Monte Carlo simulations to study particle dynamics in a potential with a Gaussian distribution of energy values, where the long-time limit of the diffusion coefficient is known from theory. In our experiments, individual colloidal particles are exposed to a laser speckle pattern inducing a non-Gaussian roughness and are followed by optical microscopy. The relaxation depends on the kind and degree of roughness of the PEL. It can be followed and quantified by the time- and ensemble-averaged mean squared displacement. Moreover, the heterogeneity of the dynamics is characterized using single-trajectory analysis. The results of this work are relevant for the correct interpretation of single-particle tracking experiments in general.

Characterizing anomalous diffusion in crowded polymer solutions and gels over five decades in time with variable-lengthscale fluorescence correlation spectroscopy

The diffusion of macromolecules in cells and in complex fluids is often found to deviate from simple Fickian diffusion. One explanation offered for this behavior is that molecular crowding renders diffusion anomalous, where the mean-squared displacement of the particles scales as 〈r(2)〉∝t(α) with α < 1. Unfortunately, methods such as fluorescence correlation spectroscopy (FCS) or fluorescence recovery after photobleaching (FRAP) probe diffusion only over a narrow range of lengthscales and cannot directly test the dependence of the mean-squared displacement (MSD) on time. Here we show that variable-lengthscale FCS (VLS-FCS), where the volume of observation is varied over several orders of magnitude, combined with a numerical inversion procedure of the correlation data, allows retrieving the MSD for up to five decades in time, bridging the gap between diffusion experiments performed at different lengthscales. In addition, we show that VLS-FCS provides a way to assess whether the propagator associated with the diffusion is Gaussian or non-Gaussian. We used VLS-FCS to investigate two systems where anomalous diffusion had been previously reported. In the case of dense cross-linked agarose gels, the measured MSD confirmed that the diffusion of small beads was anomalous at short lengthscales, with a cross-over to simple diffusion around ≈1 μm, consistent with a caged diffusion process. On the other hand, for solutions crowded with marginally entangled dextran molecules, we uncovered an apparent discrepancy between the MSD, found to be linear, and the propagators at short lengthscales, found to be non-Gaussian. These contradicting features call to mind the "anomalous, yet Brownian" diffusion observed in several biological systems, and the recently proposed "diffusing diffusivity" model.

Anomalous dynamics of intruders in a crowded environment of mobile obstacles

Many natural and industrial processes rely on constrained transport, such as proteins moving through cells, particles confined in nanocomposite materials or gels, individuals in highly dense collectives and vehicular traffic conditions. These are examples of motion through crowded environments, in which the host matrix may retain some glass-like dynamics. Here we investigate constrained transport in a colloidal model system, in which dilute small spheres move in a slowly rearranging, glassy matrix of large spheres. Using confocal differential dynamic microscopy and simulations, here we discover a critical size asymmetry, at which anomalous collective transport of the small particles appears, manifested as a logarithmic decay of the density autocorrelation functions. We demonstrate that the matrix mobility is central for the observed anomalous behaviour. These results, crucially depending on size-induced dynamic asymmetry, are of relevance for a wide range of phenomena ranging from glassy systems to cell biology.

The classical Lorentz gas model is widely used to describe constrained transport, but its assumption of an immobile environment is not applicable to many biological and industrial processes. Here, the authors show that the mobility of the matrix induces anomalous, logarithmic dynamics of the confined particles.

Splitting of the Universality Class of Anomalous Transport in Crowded Media

We investigate the emergence of subdiffusive transport by obstruction in continuum models for molecular crowding. While the underlying percolation transition for the accessible space displays universal behavior, the dynamic properties depend in a subtle nonuniversal way on the transport through narrow channels. At the same time, the different universality classes are robust with respect to introducing correlations in the obstacle matrix as we demonstrate for quenched hard-sphere liquids as underlying structures. Our results confirm that the microscopic dynamics can dominate the relaxational behavior even at long times, in striking contrast to glassy dynamics.

Anomalous Magnetotransport in Disordered Structures: Classical Edge-State Percolation

By event-driven molecular dynamics simulations we investigate magnetotransport in a two-dimensional model with randomly distributed scatterers close to the field-induced localization transition. This transition is generated by percolating skipping orbits along the edges of obstacle clusters. The dynamic exponents differ significantly from those of the conventional transport problem on percolating systems, thus establishing a new dynamic universality class. This difference is tentatively attributed to a weak-link scenario, which emerges naturally due to barely overlapping edge trajectories. We make predictions for the frequency-dependent conductivity and discuss implications for active colloidal circle swimmers in a hetegogeneous environment.

Horizons and free-path distributions in quasiperiodic Lorentz gases

We study the structure of quasiperiodic Lorentz gases, i.e., particles bouncing elastically off fixed obstacles arranged in quasiperiodic lattices. By employing a construction to embed such structures into a higher-dimensional periodic hyperlattice, we give a simple and efficient algorithm for numerical simulation of the dynamics of these systems. This same construction shows that quasiperiodic Lorentz gases generically exhibit a regime with infinite horizon, that is, empty channels through which the particles move without colliding, when the obstacles are small enough; in this case, the distribution of free paths is asymptotically a power law with exponent -3, as expected from infinite-horizon periodic Lorentz gases. For the critical radius at which these channels disappear, however, a new regime with locally finite horizon arises, where this distribution has an unexpected exponent of -5, previously observed only in a Lorentz gas formed by superposing three incommensurable periodic lattices in the Boltzmann-Grad limit where the radius of the obstacles tends to zero.

Geometry controlled anomalous diffusion in random fractal geometries: looking beyond the infinite cluster

We investigate the ergodic properties of a random walker performing (anomalous) diffusion on a random fractal geometry. Extensive Monte Carlo simulations of the motion of tracer particles on an ensemble of realisations of percolation clusters are performed for a wide range of percolation densities. Single trajectories of the tracer motion are analysed to quantify the time averaged mean squared displacement (MSD) and to compare this with the ensemble averaged MSD of the particle motion. Other complementary physical observables associated with ergodicity are studied, as well. It turns out that the time averaged MSD of individual realisations exhibits non-vanishing fluctuations even in the limit of very long observation times as the percolation density approaches the critical value. This apparent non-ergodic behaviour concurs with the ergodic behaviour on the ensemble averaged level. We demonstrate how the non-vanishing fluctuations in single particle trajectories are analytically expressed in terms of the fractal dimension and the cluster size distribution of the random geometry, thus being of purely geometrical origin. Moreover, we reveal that the convergence scaling law to ergodicity, which is known to be inversely proportional to the observation time T for ergodic diffusion processes, follows a power-law ∼T(-h) with h < 1 due to the fractal structure of the accessible space. These results provide useful measures for differentiating the subdiffusion on random fractals from an otherwise closely related process, namely, fractional Brownian motion. Implications of our results on the analysis of single particle tracking experiments are provided.

Diffusion in Jammed Particle Packs

Using random walk simulations we explore diffusive transport through monodisperse sphere packings over a range of packing fractions ϕ in the vicinity of the jamming transition at ϕ(c). Various diffusion properties are computed over several orders of magnitude in both time and packing pressure. Two well-separated regimes of normal "Fickian" diffusion, where the mean squared displacement is linear in time, are observed. The first corresponds to diffusion inside individual spheres, while the latter is the long-time bulk diffusion. The intermediate anomalous diffusion regime and the long-time value of the diffusion coefficient are both shown to be controlled by particle contacts, which in turn depend on proximity to ϕ(c). The time required to recover normal diffusion t* scales as (ϕ-ϕ(c))(-0.5) and the long-time diffusivity D(∞)∼(ϕ-ϕ(c))0.5, or D(∞)∼1/t*. It is shown that the distribution of mean first passage times associated with the escape of random walkers between neighboring particles controls both t* and D(∞) in the limit ϕ→ϕ(c).

Superdiffusion dominates intracellular particle motion in the supercrowded cytoplasm of pathogenic Acanthamoeba castellanii

Acanthamoebae are free-living protists and human pathogens, whose cellular functions and pathogenicity strongly depend on the transport of intracellular vesicles and granules through the cytosol. Using high-speed live cell imaging in combination with single-particle tracking analysis, we show here that the motion of endogenous intracellular particles in the size range from a few hundred nanometers to several micrometers in Acanthamoeba castellanii is strongly superdiffusive and influenced by cell locomotion, cytoskeletal elements, and myosin II. We demonstrate that cell locomotion significantly contributes to intracellular particle motion, but is clearly not the only origin of superdiffusivity. By analyzing the contribution of microtubules, actin, and myosin II motors we show that myosin II is a major driving force of intracellular motion in A. castellanii. The cytoplasm of A. castellanii is supercrowded with intracellular vesicles and granules, such that significant intracellular motion can only be achieved by actively driven motion, while purely thermally driven diffusion is negligible.