Differences between lattice and continuum percolation transport exponents

Relation between critical exponent of the conductivity and the morphological exponents of percolation theory

A central unsolved problem in percolation theory over the past five decades has been whether there is a direct relationship between the critical exponents that characterize the power-law behavior of the transport properties near the percolation threshold, particularly the effective electrical conductivity σ_{e}, and the exponents that describe the morphology of percolation clusters. The problem is also relevant to the relation between the static exponents of percolation clusters and the critical dynamics of spin waves in dilute ferromagnets, the elasticity of gels and composite solids, hopping conductivity in semiconductors, solute transport in porous media, and many others. We propose an approach to address the problem by showing that the contributions to σ_{e} can be decomposed into several groups representing the structure of percolation networks, including their mass and tortuosity, as well as constrictivity that describes the fluctuations in the driving potential gradient along the transport paths. The decomposition leads to a relationship between the critical exponent t of σ_{e} and other percolation exponents in d dimensions, t/ν=(d-D_{bb})+2(D_{op}-1)+d_{C}, where ν, D_{bb}, D_{op}, and d_{C} are, respectively, the correlation length exponent, the fractal dimensions of the backbones and the optimal paths, and the exponent that characterizes the constrictivity. Numerical simulations in two and three dimensions, as well as analytical results in d=1 and d=6, the upper critical dimension of percolation, validate the relationship. We, therefore, believe that the solution to the 50-year-old problem has been derived.

Structural, mechanical, and vibrational properties of particulate physical gels

Our lives are surrounded by a rich assortment of disordered materials. In particular, glasses are well known as dense, amorphous materials, whereas gels exist in low-density, disordered states. Recent progress has provided a significant step forward in understanding the material properties of glasses, such as mechanical, vibrational, and transport properties. In contrast, our understanding of particulate physical gels is still highly limited. Here, using molecular dynamics simulations, we study a simple model of particulate physical gels, the Lennard-Jones (LJ) gels, and provide a comprehensive understanding of their structural, mechanical, and vibrational properties, all of which are markedly different from those of LJ glasses. First, the LJ gels show sparse, heterogeneous structures, and the length scale ξ_s of the structures grows as the density is lowered. Second, the LJ gels are extremely soft, with both shear G and bulk K moduli being orders of magnitude smaller than those of LJ glasses. Third, many low-frequency vibrational modes are excited, which form a characteristic plateau with the onset frequency ω_* in the vibrational density of states. Structural, mechanical, and vibrational properties, characterized by ξ_s, G, K, and ω_*, respectively, show power-law scaling behaviors with the density, which establishes a close relationship between them. Throughout this work, we also reveal that LJ gels are multiscale, solid-state materials: (i) homogeneous elastic bodies at long lengths, (ii) heterogeneous elastic bodies with fractal structures at intermediate lengths, and (iii) amorphous structural bodies at short lengths.

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.

Permeability of packs of polydisperse hard spheres

The permeability of packs of spheres is important in a wide range of physical scenarios. Here, we create numerically generated random periodic domains of spheres that are polydisperse in size and use lattice-Boltzmann simulations of fluid flow to determine the permeability of the pore phase interstitial to the spheres. We control the polydispersivity of the sphere size distribution and the porosity across the full range from high porosity to a close packing of spheres. We find that all results scale with a Stokes permeability adapted for polydisperse sphere sizes. We show that our determination of the permeability of random distributions of spheres is well approximated by models for cubic arrays of spheres at porosities greater than ∼0.38, without any fitting parameters. Below this value, the Kozeny-Carman relationship provides a good approximation for dense, closely packed sphere packs across all polydispersivity.

Metal-Insulator Transition of Ultrathin Sputtered Metals on Phenolic Resin Thin Films: Growth Morphology and Relations to Surface Free Energy and Reactivity

Nanostructured metal assemblies on thin and ultrathin polymeric films enable state of the art technologies and have further potential in diverse fields. Rational design of the structure–function relationship is of critical importance but aggravated by the scarcity of systematic studies. Here, we studied the influence of the interplay between metal and polymer surface free energy and reactivity on the evolution of electric conductivity and the resulting morphologies. In situ resistance measurements during sputter deposition of Ag, Au, Cu and Ni films on ultrathin reticulated polymer films collectively reveal metal–insulator transitions characteristic for Volmer–Weber growth. The different onsets of percolation correlate with interfacial energy and energy of adhesion weakly but as expected from ordinary wetting theory. A more pronounced trend of lower percolation thickness for more reactive metals falls in line with reported correlations. Ex situ grazing incidence small angle X-ray scattering experiments were performed at various thicknesses to gain an insight into cluster and film morphology evolution. A novel approach to interpret the scattering data is used where simulated pair distance distributions of arbitrary shapes and arrangements can be fitted to experiments. Detailed approximations of cluster structures could be inferred and are discussed in view of the established parameters describing film growth behavior.

Ant collective cognition allows for efficient navigation through disordered environments

The cognitive abilities of biological organisms only make sense in the context of their environment. Here, we study longhorn crazy ant collective navigation skills within the context of a semi-natural, randomized environment. Mapping this biological setting into the ‘Ant-in-a-Labyrinth’ framework which studies physical transport through disordered media allows us to formulate precise links between the statistics of environmental challenges and the ants’ collective navigation abilities. We show that, in this environment, the ants use their numbers to collectively extend their sensing range. Although this extension is moderate, it nevertheless allows for extremely fast traversal times that overshadow known physical solutions to the ‘Ant-in-a-Labyrinth’ problem. To explain this large payoff, we use percolation theory and prove that whenever the labyrinth is solvable, a logarithmically small sensing range suffices for extreme speedup. Overall, our work demonstrates the potential advantages of group living and collective cognition in increasing a species’ habitable range.

Electromagnetic Properties of Carbon Gels

The electromagnetic properties of various carbon gels, produced with different bulk densities, were investigated in a wide frequency range (20 Hz-36 GHz). The values of dielectric permittivity and electrical conductivity at 129 Hz were found to be very high, i.e., more than 10⁵ and close to 100 S/m, respectively. Both strongly decreased with frequency but remained high in the microwave frequency range (close to 10 and about 0.1 S/m, respectively, at 30 GHz). Moreover, the dielectric permittivity and the electrical conductivity strongly increased with the bulk density of the materials, according to power laws at low frequency. However, the maximum of microwave absorption was observed at lower densities. The DC conductivity slightly decreased on cooling, according to the Arrhenius law. The lower activation energies are typical of carbon gels presenting lower DC electrical conductivities, due to a higher number of defects. High and thermally stable electromagnetic properties of carbon gels, together with other unique properties of these materials, such as lightness and chemical inertness, open possibilities for producing new electromagnetic coatings.

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.

Fluid flow through packings of elastic shells

Fluid transport in porous materials is commonly studied in geological samples (soil, sediments, etc.) or idealized systems, but the fluid flow through compacted granular materials, consisting of substantially strained granules, remains relatively unexplored. As a step toward filling this gap, we study a model of liquid transport in packings of deformable elastic shells using finite-element and lattice-Boltzmann methods. We find that the fluid flow abruptly vanishes as the porosity of the material falls below a critical value, and the flow obstruction exhibits features of a percolation transition. We further show that the fluid flow can be captured by a simplified permeability model in which the complex porous material is replaced by a collection of disordered capillaries, which are distributed and shaped by the percolation transition. To that end, we numerically explore the divergence of hydraulic tortuosity τ_{H} and the decrease of a hydraulic radius R_{h} as the percolation threshold is approached. We interpret our results in terms of scaling predictions derived from the percolation theory applied to random packings of spheres.

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.

Glassy relaxation slows down by increasing mobility

We find a striking trend reversal in the relaxation dynamics of mixtures with strong dynamical asymmetry. Simulations by both Brownian and Newtonian dynamics reveal that in mixtures of fast and slow hard spheres, above a critical density, the dynamics becomes slower upon increasing the mobility of the fast particles. Below that density, the same increase in mobility speeds up the dynamics. The critical density itself can be identified with the glass transition of the mode-coupling theory that does not depend on the dynamical asymmetry. The asymptotic dynamics close to the critical density is universal, but strong pre-asymptotic effects prevail in particular when the dynamical asymmetry also involves size asymmetry. Our observations reconcile earlier findings, where a strong dependence on kinetic parameters was found for the glassy dynamics, with the paradigm that the glass transition is determined by the properties of configuration space alone.

Thermoplastic Elastomer Systems Containing Carbon Nanofibers as Soft Piezoresistive Sensors

Soft, wearable or printable strain sensors derived from conductive polymer nanocomposites (CPNs) are becoming increasingly ubiquitous in personal-care applications. Common elastomers employed in the fabrication of such piezoresistive CPNs frequently rely on chemically cross-linked polydiene or polysiloxane chemistry, thereby generating relatively inexpensive and reliable sensors that become solid waste upon application termination. Moreover, the shape anisotropy of the incorporated conductive nanoparticles can produce interesting electrical effects due to strain-induced spatial rearrangement. In this study, we investigate the morphological, mechanical, electrical, and electromechanical properties of CPNs generated from thermoplastic elastomer (TPE) triblock copolymer systems containing vapor-grown carbon nanofiber (CNF). Modulus-tunable TPE gels imbibed with a midblock-selective aliphatic oil exhibit well-behaved properties with increasing CNF content, but generally display nonlinear negative piezoresistance at different strain amplitudes and stretch rates due to nanofiber mobility upon CPN strain-cycling. In contrast, a neat TPE possessing low hard-block content yields a distinctive strain-reversible piezoresistive response, as well as low electrical hysteresis, upon cyclic deformation. Unlike their chemically cross-linked analogs, these physically cross-linked and thus environmentally benign CPNs are fully reprocessable by thermal and/or solvent means.

Percolation probability in a system of cylindrical particles

A broad variety of materials, ranging from composites and heat transfer nano-fluids to electrochemical energy storage electrodes, widely employ cylindrical particles of various aspect ratios, such as carbon nanotubes. These particles are generally excellent conductors of heat and electricity and when dispersed in a continuous medium influence dramatically the transport properties of the heterogeneous material by forming a percolating network. Numerous theories exist to predict key parameters such as particle concentration at the percolation threshold and transport properties at concentrations beyond the threshold. The microstructure formed by connecting particles in the material is an important determinant toward such parameters but often requires complex numerical models to resolve. In this paper, we present an analytical, probabilistic model capturing the microstructure of a system of randomly positioned, soft-core, cylindrical particles with a finite aspect ratio, valid at arbitrary particle concentration. Our analytical framework allows for the calculation of the particle contact number distribution and percolation probability of the particle system. We show that our analytical model is more accurate than excluded volume theory for predicting the percolation threshold for spherocylinders of finite aspect ratios, and agrees well with the corresponding numerical results. Our theory describes the percolating network topology above the percolation threshold and can serve as the foundation for analytical composition-structure-property relationships for heterogeneous materials with conducting cylindrical particles.

Diffusion of small particles in polymer films

The motion of small probe molecules in a two-dimensional system containing frozen polymer chains was studied by means of Monte Carlo simulations. The model macromolecules were coarse-grained and restricted to vertices of a triangular lattice. The cooperative motion algorithm was used to generate representative configurations of macromolecular systems of different polymer concentrations. The remaining unoccupied lattice sites of the system were filled with small molecules. The structure of the polymer film, especially near the percolation threshold, was determined. The dynamic lattice liquid algorithm was then employed for studies of the dynamics of small objects in the polymer matrix. The influence of chain length and polymer concentration on the mobility and the character of motion of small molecules were studied. Short- and long-time dynamic behaviors of solvent molecules were also described. Conditions of anomalous diffusions' appearance in such systems are discussed. The influence of the structure of the matrix of obstacles on the molecular transport was discussed.

Simple Rules Govern the Patterns of Arctic Sea Ice Melt Ponds

Climate change, amplified in the far north, has led to rapid sea ice decline in recent years. In the summer, melt ponds form on the surface of Arctic sea ice, significantly lowering the ice reflectivity (albedo) and thereby accelerating ice melt. Pond geometry controls the details of this crucial feedback; however, a reliable model of pond geometry does not currently exist. Here we show that a simple model of voids surrounding randomly sized and placed overlapping circles reproduces the essential features of pond patterns. The only two model parameters, characteristic circle radius and coverage fraction, are chosen by comparing, between the model and the aerial photographs of the ponds, two correlation functions which determine the typical pond size and their connectedness. Using these parameters, the void model robustly reproduces the ponds' area-perimeter and area-abundance relationships over more than 6 orders of magnitude. By analyzing the correlation functions of ponds on several dates, we also find that the pond scale and the connectedness are surprisingly constant across different years and ice types. Moreover, we find that ponds resemble percolation clusters near the percolation threshold. These results demonstrate that the geometry and abundance of Arctic melt ponds can be simply described, which can be exploited in future models of Arctic melt ponds that would improve predictions of the response of sea ice to Arctic warming.

Fabrication of Graphene-Polyimide Nanocomposites with Superior Electrical Conductivity

We report on the fabrication of a novel class of lightweight materials, polyimide-graphene nanocomposites (0.01-5 vol %), with tunable electrical conductivity. The graphene-polyimide nanocomposites exhibit an ultra-low graphene percolation threshold of 0.03 vol % and maximum dc conductivity of 0.94 S/cm, which we attribute to excellent dispersion, extraordinary electron transport in the well-dispersed graphene, high number density of graphene nanosheets, and the π-π interactions between the aromatic moieties of the polyimide and the carbon rings in graphene. The dc conductivity data are shown to follow the power-law dependence on the graphene volume fraction near the percolation threshold. The ac conductivity of the nanocomposites is accurately represented by the extended pair-approximation model. The exponent s of the approximation is estimated to be 0.45-0.61, indicating anomalous diffusion of charge particles and a fractal structure for the conducting phase, lending support to the percolation model. Low-temperature dc conductivity of the nanocomposites is well-approximated by the thermal fluctuation-induced tunneling. Wide-angle X-ray scattering and transmission electron microscopy were utilized to correlate the morphology with the electrical conductivity. The lack of maxima in X-ray indicates the loss of structural registry and short-range ordering.

Behavior of the supercritical phase of a continuum percolation model on ℝd

A continuum percolation model onis considered. Using a renormalization technique developed by Grimmett and Marstrand, we show a continuum analogue of their results. We prove the critical value of the percolation equals the limit of the critical value of a slice as the thickness of the slice tends to infinity. We also prove that the effective conductivity in the model is bounded from below by a positive constant in the supercritical case.

Percolative switching in transition metal dichalcogenide field-effect transistors at room temperature

We have addressed the microscopic transport mechanism at the switching or 'on-off' transition in transition metal dichalcogenide (TMDC) field-effect transistors (FETs), which has been a controversial topic in TMDC electronics, especially at room temperature. With simultaneous measurement of channel conductivity and its slow time-dependent fluctuation (or noise) in ultrathin WSe2 and MoS2 FETs on insulating SiO2 substrates where noise arises from McWhorter-type carrier number fluctuations, we establish that the switching in conventional backgated TMDC FETs is a classical percolation transition in a medium of inhomogeneous carrier density distribution. From the experimentally observed exponents in the scaling of noise magnitude with conductivity, we observe unambiguous signatures of percolation in a random resistor network, particularly, in WSe2 FETs close to switching, which crosses over to continuum percolation at a higher doping level. We demonstrate a powerful experimental probe to the microscopic nature of near-threshold electrical transport in TMDC FETs, irrespective of the material detail, device geometry, or carrier mobility, which can be extended to other classes of 2D material-based devices as well.

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.

Rounding of the localization transition in model porous media

The generic mechanisms of anomalous transport in porous media are investigated by computer simulations of two-dimensional model systems. In order to bridge the gap between the strongly idealized Lorentz model and realistic models of porous media, two models of increasing complexity are considered: a cherry-pit model with hard-core correlations as well as a soft-potential model. An ideal gas of tracer particles inserted into these structures is found to exhibit anomalous transport which extends up to several decades in time. Also, the self-diffusion of the tracers becomes suppressed upon increasing the density of the systems. These phenomena are attributed to an underlying percolation transition. In the soft potential model the transition is rounded, since each tracer encounters its own critical density according to its energy. Therefore, the rounding of the transition is a generic occurrence in realistic, soft systems.

A scalable and facile synthesis of alumina/exfoliated graphite composites by attrition milling

We present a facile one-pot synthesis of alumina/exfoliated graphite composite having excellent electrical conductivity (>1,000 S m⁻¹), fracture toughness (5.6 MPa m^0.5), and wear resistance, which is enhanced by 7.7 times compared to pure alumina.

Simulation of diffusion in a crowded environment

We performed extensive and systematic simulation studies of two-dimensional fluid motion in a complex crowded environment. In contrast to other studies we focused on cooperative phenomena that occurred if the motion of particles takes place in a dense crowded system, which can be considered as a crude model of a cellular membrane. Our main goal was to answer the following question: how do the fluid molecules move in an environment with a complex structure, taking into account the fact that motions of fluid molecules are highly correlated. The dynamic lattice liquid (DLL) model, which can work at the highest fluid density, was employed. Within the frame of the DLL model we considered cooperative motion of fluid particles in an environment that contained static obstacles. The dynamic properties of the system as a function of the concentration of obstacles were studied. The subdiffusive motion of particles was found in the crowded system. The influence of hydrodynamics on the motion was investigated via analysis of the displacement in closed cooperative loops. The simulation and the analysis emphasize the influence of the movement correlation between moving particles and obstacles.

Electric conductivity percolation in naturally dehydrating, lightly wetted, hydrophilic fumed silica powder

In studying the dehydration of surface-moistened fumed silica Aerosil powders, we found a conductivity percolation transition at low hydration levels. Both the percolation exponent and the threshold are typical for correlated site-bond transitions in complex two-dimensional (2D) systems. The exponent values, 0.94-1.10, are indicative of severe heterogeneity in the conducting medium. The surface moisture at the percolation threshold takes on a universal value of 0.65 mg([H2O])/m(2)([silica]), independent of the silica grain size, and equivalent to twice the first hydration monolayer. This level is just sufficient to sustain a quasi-2D, hydrogen-bonded water network spanning the silica surface.

Anomalous transport in the crowded world of biological cells

A ubiquitous observation in cell biology is that the diffusive motion of macromolecules and organelles is anomalous, and a description simply based on the conventional diffusion equation with diffusion constants measured in dilute solution fails. This is commonly attributed to macromolecular crowding in the interior of cells and in cellular membranes, summarizing their densely packed and heterogeneous structures. The most familiar phenomenon is a sublinear, power-law increase of the mean-square displacement (MSD) as a function of the lag time, but there are other manifestations like strongly reduced and time-dependent diffusion coefficients, persistent correlations in time, non-Gaussian distributions of spatial displacements, heterogeneous diffusion and a fraction of immobile particles. After a general introduction to the statistical description of slow, anomalous transport, we summarize some widely used theoretical models: Gaussian models like fractional Brownian motion and Langevin equations for visco-elastic media, the continuous-time random walk model, and the Lorentz model describing obstructed transport in a heterogeneous environment. Particular emphasis is put on the spatio-temporal properties of the transport in terms of two-point correlation functions, dynamic scaling behaviour, and how the models are distinguished by their propagators even if the MSDs are identical. Then, we review the theory underlying commonly applied experimental techniques in the presence of anomalous transport like single-particle tracking, fluorescence correlation spectroscopy (FCS) and fluorescence recovery after photobleaching (FRAP). We report on the large body of recent experimental evidence for anomalous transport in crowded biological media: in cyto- and nucleoplasm as well as in cellular membranes, complemented by in vitro experiments where a variety of model systems mimic physiological crowding conditions. Finally, computer simulations are discussed which play an important role in testing the theoretical models and corroborating the experimental findings. The review is completed by a synthesis of the theoretical and experimental progress identifying open questions for future investigation.

Permeability of porous materials determined from the Euler characteristic

We study the permeability of quasi-two-dimensional porous structures of randomly placed overlapping monodisperse circular and elliptical grains. Measurements in microfluidic devices and lattice Boltzmann simulations demonstrate that the permeability is determined by the Euler characteristic of the conducting phase. We obtain an expression for the permeability that is independent of the percolation threshold and shows agreement with experimental and simulated data over a wide range of porosities. Our approach suggests that the permeability explicitly depends on the overlapping probability of grains rather than their shape.

Effect of polydispersity on diffusion in random obstacle matrices

The dynamics of tracers in disordered matrices is of interest in a number of diverse areas of physics such as the biophysics of crowding in cells and cell membranes, and the diffusion of fluids in porous media. To a good approximation the matrices can be modeled as a collection of spatially frozen particles. In this Letter, we consider the effect of polydispersity (in size) of the matrix particles on the dynamics of tracers. We study a two dimensional system of hard disks diffusing in a sea of hard disk obstacles, for different values of the polydispersity of the matrix. We find that for a given average size and area fraction, the diffusion of tracers is very sensitive to the polydispersity. We calculate the pore percolation threshold using Apollonius diagrams. The diffusion constant, D, follows a scaling relation D~(φ(c)-φ(m))(μ-β) for all values of the polydispersity, where φ(m) is the area fraction and φ(c) is the value of φ(m) at the percolation threshold.

Dispersion in porous media, continuous-time random walks, and percolation

A promising approach to the modeling of anomalous (non-Gaussian) dispersion in flow through heterogeneous porous media is the continuous-time random walk (CTRW) method. In such a formula on the waiting time distribution ψ(t) is usually assumed to be given by ψ(t)∼t-1-α, with α fitted to the experimental data. The exponent α is also related to the power-law growth of the mean-square displacement of the solute with the time t <R2(t)> ∼ tζ. Invoking percolation and using a scaling analysis, we relate α to the geometrical exponents of percolation (ν, β, and βB) as well as the exponents μ and e that characterize the power-law behavior of the effective conductivity and permeability of porous media near the percolation threshold. We then explain the cause of the nonuniversality of α in terms of the nonuniversality of μ and e in continuum systems, and in percolation models with long-range correlations, and propose bounds for it. The results are consistent with the experimental data, both at the laboratory and field scales.

Permeability, porosity, and percolation properties of two-dimensional disordered fracture networks

Using extensive Monte Carlo simulations, we study the effective permeability, porosity, and percolation properties of two-dimensional fracture networks in which the fractures are represented by rectangles of finite widths. The parameters of the study are the width of the fractures and their number density. For low and intermediate densities, the average porosity of the network follows a power-law relation with the density. The exponent of the power law itself depends on the fractures' width through a power law. For an intermediate range of the densities, the effective permeability scales with the fractures' width as a power law, with an exponent that depends on the density. For high densities the effective permeability also depends on the porosity through a power law, with an exponent that depends on the fractures' width. In agreement with the results, experimental data also indicate the existence of a power-law relationship between the effective permeability and porosity in consolidated sandstones and sedimentary rocks with a nonuniversal exponent. The percolation threshold or critical number density of the fractures depends on their width and is maximum if they are represented by squares, rather than rectangles.

Anomalous transport of a tracer on percolating clusters

We investigate the dynamics of a single tracer exploring a course of fixed obstacles in the vicinity of the percolation transition for particles confined to the infinite cluster. The mean-square displacement displays anomalous transport, which extends to infinite times precisely at the critical obstacle density. The slowing down of the diffusion coefficient exhibits power-law behavior for densities close to the critical point and we show that the mean-square displacement fulfills a scaling hypothesis. Furthermore, we calculate the dynamic conductivity as a response to an alternating electric field. Last, we discuss the non-gaussian parameter as an indicator for heterogeneous dynamics.

Diffusion amid random overlapping obstacles: similarities, invariants, approximations

Efficient and accurate numerical techniques are used to examine similarities of effective diffusion in a void between random overlapping obstacles: essential invariance of effective diffusion coefficients (D(eff)) with respect to obstacle shapes and applicability of a two-parameter power law over nearly entire range of excluded volume fractions (φ), except for a small vicinity of a percolation threshold. It is shown that while neither of the properties is exact, deviations from them are remarkably small. This allows for quick estimation of void percolation thresholds and approximate reconstruction of D(eff) (φ) for obstacles of any given shape. In 3D, the similarities of effective diffusion yield a simple multiplication "rule" that provides a fast means of estimating D(eff) for a mixture of overlapping obstacles of different shapes with comparable sizes.

Transport on exploding percolation clusters

We propose a simple generalization of the explosive percolation process [Achlioptas et al., Science 323, 1453 (2009)], and investigate its structural and transport properties. In this model, at each step, a set of q unoccupied bonds is randomly chosen. Each of these bonds is then associated with a weight given by the product of the cluster sizes that they would potentially connect, and only that bond among the q set which has the smallest weight becomes occupied. Our results indicate that, at criticality, all finite-size scaling exponents for the spanning cluster, the conducting backbone, the cutting bonds, and the global conductance of the system, change continuously and significantly with q. Surprisingly, we also observe that systems with intermediate values of q display the worst conductive performance. This is explained by the strong inhibition of loops in the spanning cluster, resulting in a substantially smaller associated conducting backbone.

Two-dimensional continuum percolation threshold for diffusing particles of nonzero radius

Lateral diffusion in the plasma membrane is obstructed by proteins bound to the cytoskeleton. The most important parameter describing obstructed diffusion is the percolation threshold. The thresholds are well known for point tracers, but for tracers of nonzero radius, the threshold depends on the excluded area, not just the obstacle concentration. Here thresholds are obtained for circular obstacles on the continuum. Random obstacle configurations are generated by Brownian dynamics or Monte Carlo methods, the obstacles are immobilized, and the percolation threshold is obtained by solving a bond percolation problem on the Voronoi diagram of the obstacles. The percolation threshold is expressed as the diameter of the largest tracer that can cross a set of immobile obstacles at a prescribed number density. For random overlapping obstacles, the results agree with the known analytical solution quantitatively. When the obstacles are soft disks with a 1/r(12) repulsion, the percolating diameter is approximately 20% lower than for overlapping obstacles. A percolation model predicts that the threshold is highly sensitive to the tracer radius. To our knowledge, such a strong dependence has so far not been reported for the plasma membrane, suggesting that percolation is not the factor controlling lateral diffusion. A definitive experiment is proposed.

Explosive percolation in a nanotube-based system

Using the percolation theory, we study the underlying mechanism in the formation of single-walled nanotube bundles with uniform diameter. By applying the cluster repulsion process to stick percolation, we find that the transition becomes explosive. To understand the transition nature, we first investigate the scaling behavior of transition interval Δ. By comparing the results with loopless and loop-allowed bond percolations, we find that the loops crucially affect the scaling behavior of Δ, and Δ is not universal. Moreover, the scaling behavior of Δ for the present nanostick systems is the same as that for loopless bond percolation. For more systematic studies on the transition nature, we also measure the changes in order parameter during the stick removal process and show that there exists a hysteresis. The results more clearly show that the transition of the stick system with cluster repulsion is discontinuous.