Equivalence of the void percolation problem for overlapping spheres and a network problem

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 pore radius and transport properties of disordered hard- and overlapping-sphere models

Transport properties of porous media are intimately linked to their pore-space microstructures. We quantify geometrical and topological descriptors of the pore space of certain disordered and ordered distributions of spheres, including pore-size functions and the critical pore radius δ_{c}. We focus on models of porous media derived from maximally random jammed sphere packings, overlapping spheres, equilibrium hard spheres, quantizer sphere packings, and crystalline sphere packings. For precise estimates of the percolation thresholds, we use a strict relation of the void percolation around sphere configurations to weighted bond percolation on the corresponding Voronoi networks. We use the Newman-Ziff algorithm to determine the percolation threshold using universal properties of the cluster size distribution. The critical pore radius δ_{c} is often used as the key characteristic length scale that determines the fluid permeability k. A recent study [Torquato, Adv. Wat. Resour. 140, 103565 (2020)10.1016/j.advwatres.2020.103565] suggested for porous media with a well-connected pore space an alternative estimate of k based on the second moment of the pore size 〈δ^{2}〉, which is easier to determine than δ_{c}. Here, we compare δ_{c} to the second moment of the pore size 〈δ^{2}〉, and indeed confirm that, for all porosities and all models considered, δ_{c}^{2} is to a good approximation proportional to 〈δ^{2}〉. However, unlike 〈δ^{2}〉, the permeability estimate based on δ_{c}^{2} does not predict the correct ranking of k for our models. Thus, we confirm 〈δ^{2}〉 to be a promising candidate for convenient and reliable estimates of the fluid permeability for porous media with a well-connected pore space. Moreover, we compare the fluid permeability of our models with varying degrees of order, as measured by the τ order metric. We find that (effectively) hyperuniform models tend to have lower values of k than their nonhyperuniform counterparts. Our findings could facilitate the design of porous media with desirable transport properties via targeted pore statistics.

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.

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.

CAVD, towards better characterization of void space for ionic transport analysis

Geometric crystal structure analysis using three-dimensional Voronoi tessellation provides intuitive insights into the ionic transport behavior of metal-ion electrode materials or solid electrolytes by mapping the void space in a framework onto a network. The existing tools typically consider only the local voids by mapping them with Voronoi polyhedra vertices and then define the mobile ions pathways using the Voronoi edges connecting these vertices. We show that in some structures mobile ions are located on Voronoi polyhedra faces and thus cannot be located by a standard approach. To address this deficiency, we extend the method to include Voronoi faces in the constructed network. This method has been implemented in the CAVD python package. Its effectiveness is demonstrated by 99% recovery rate for the lattice sites of mobile ions in 6,955 Li-, Na-, Mg- and Al-containing ionic compounds extracted from the Inorganic Crystal Structure Database. In addition, various quantitative descriptors of the network can be used to identify and rank the materials and further used in materials databases for machine learning.

Percolation phase transition by removal of k^{2}-mers from fully occupied lattices

Numerical simulations and finite-size scaling analysis have been carried out to study the problem of inverse site percolation by the removal of k×k square tiles (k^{2}-mers) from square lattices. The process starts with an initial configuration, where all lattice sites are occupied and, obviously, the opposite sides of the lattice are connected by occupied sites. Then the system is diluted by removing k^{2}-mers of occupied sites from the lattice following a random sequential adsorption mechanism. The process finishes when the jamming state is reached and no more objects can be removed due to the absence of occupied sites clusters of appropriate size and shape. The central idea of this paper is based on finding the maximum concentration of occupied sites, p_{c,k}, for which the connectivity disappears. This particular value of the concentration is called the inverse percolation threshold and determines a well-defined geometrical phase transition in the system. The results obtained for p_{c,k} show that the inverse percolation threshold is a decreasing function of k in the range 1≤k≤4. For k≥5, all jammed configurations are percolating states, and consequently, there is no nonpercolating phase. In other words, the lattice remains connected even when the highest allowed concentration of removed sites is reached. The jamming exponent ν_{j} was measured, being ν_{j}=1 regardless of the size k considered. In addition, the accurate determination of the critical exponents ν, β, and γ reveals that the percolation phase transition involved in the system, which occurs for k varying between one and four, has the same universality class as the standard percolation problem.

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.

Void distributions reveal structural link between jammed packings and protein cores

Dense packing of hydrophobic residues in the cores of globular proteins determines their stability. Recently, we have shown that protein cores possess packing fraction ϕ≈0.56, which is the same as dense, random packing of amino-acid-shaped particles. In this article, we compare the structural properties of protein cores and jammed packings of amino-acid-shaped particles in much greater depth by measuring their local and connected void regions. We find that the distributions of surface Voronoi cell volumes and local porosities obey similar statistics in both systems. We also measure the probability that accessible, connected void regions percolate as a function of the size of a spherical probe particle and show that both systems possess the same critical probe size. We measure the critical exponent τ that characterizes the size distribution of connected void clusters at the onset of percolation. We find that the cluster size statistics are similar for void percolation in packings of amino-acid-shaped particles and randomly placed spheres, but different from that for void percolation in jammed sphere packings. We propose that the connected void regions are a defining structural feature of proteins and can be used to differentiate experimentally observed proteins from decoy structures that are generated using computational protein design software. This work emphasizes that jammed packings of amino-acid-shaped particles can serve as structural and mechanical analogs of protein cores, and could therefore be useful in modeling the response of protein cores to cavity-expanding and -reducing mutations.

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.

Computer simulation of the effect of deformation on the morphology and flow properties of porous media

We report on the results of extensive computer simulation of the effect of deformation on the morphology of a porous medium and its fluid flow properties. The porous medium is represented by packings of spherical particles. Both random and regular as well as dense and nondense packings are used. A quasistatic model based on Hertz's contact theory is used to model the mechanical deformation of the packings, while the evolution of the permeability with the deformation is computed by the lattice-Boltzmann approach. The evolution of the pore-size and pore-length distributions, the porosity, the particles' contacts, the permeability, and the distribution of the stresses that the fluid exerts in the pore space are all studied in detail. The distribution of the pores' lengths, the porosity, and the particles' connectivity change strongly with the application of an external strain to the porous media, whereas the pore-size distribution is not affected as strongly. The permeability of the porous media strongly reduces even when the applied strain is small. When the permeabilities and porosities of the random packings are normalized with respect to their predeformation values, they all collapse onto a single curve, independent of the particle-size distribution. The porosity reduces as a power law with the external strain. The fluid stresses in the pore space follow roughly a log-normal distribution, both before and after deformation.

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.

Dimensional study of the dynamical arrest in a random Lorentz gas

The random Lorentz gas (RLG) is a minimal model for transport in heterogeneous media. Upon increasing the obstacle density, it exhibits a growing subdiffusive transport regime and then a dynamical arrest. Here, we study the dimensional dependence of the dynamical arrest, which can be mapped onto the void percolation transition for Poisson-distributed point obstacles. We numerically determine the arrest in dimensions d=2-6. Comparison of the results with standard mode-coupling theory reveals that the dynamical theory prediction grows increasingly worse with d. In an effort to clarify the origin of this discrepancy, we relate the dynamical arrest in the RLG to the dynamic glass transition of the infinite-range Mari-Kurchan-model glass former. Through a mixed static and dynamical analysis, we then extract an improved dimensional scaling form as well as a geometrical upper bound for the arrest. The results suggest that understanding the asymptotic behavior of the random Lorentz gas may be key to surmounting fundamental difficulties with the mode-coupling theory of glasses.

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.

Percolation through voids around overlapping spheres: a dynamically based finite-size scaling analysis

The percolation threshold for flow or conduction through voids surrounding randomly placed spheres is calculated. With large-scale Monte Carlo simulations, we give a rigorous continuum treatment to the geometry of the impenetrable spheres and the spaces between them. To properly exploit finite-size scaling, we examine multiple systems of differing sizes, with suitable averaging over disorder, and extrapolate to the thermodynamic limit. An order parameter based on the statistical sampling of stochastically driven dynamical excursions and amenable to finite-size scaling analysis is defined, calculated for various system sizes, and used to determine the critical volume fraction ϕc=0.0317±0.0004 and the correlation length exponent ν=0.92±0.05.

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.

Characterization of void space in polydisperse sphere packings: Applications to hard-sphere packings and to protein structure analysis

The implementation of a method for the exact evaluation of the volume and surface area of cavities and free volumes in polydisperse sphere packings is described. The generalization of an algorithm for Voronoi tessellation by Tanemura et al. is presented, employing the radical plane construction, as a part of the method. We employ this method to calculate the equation of state for monodisperse and polydisperse hard-sphere fluids, crystals, and for the metastable amorphous branch up to random close packing or jamming densities. We compute the distribution of free volumes, and compare with previous results employing a heuristic definition of free volume. We show the efficacy of the method for analyzing protein structure, by computing various quantities such as the distribution of sizes of buried cavities and pockets, the scaling of solvent accessible area to the corresponding occupied volume, the composition of residues lining cavities, etc.

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.

Structure of void space in polymer solutions

The structure of void space in two- and three-dimensional (3D) polymer solutions is studied using Voronoi tessellation and percolation theory. The polymer molecules are modeled as freely jointed chains of N tangent hard disks (two dimensions) or spheres (three dimensions). Polymer chains are equilibrated via Monte Carlo simulations and the pore space in configurations of equilibrated chains is mapped using Voronoi tessellation. In d dimensions a Voronoi vertex is the center of the sphere tangent to the d+1 nearest monomers. An edge of the Voronoi diagram is the shortest route between two neighboring vertices. The edge is considered connected if a monomer can pass through and disconnected otherwise. The Voronoi construction is used to calculate the percolation threshold of the void space. The most interesting result is that the polymer area fraction at the percolation threshold is a nonmonotonic function of N in two dimensions but monotonically reaches a constant value in three dimensions. The crossover behavior of the percolation threshold is also observed in pseudo-3D. The pore size distribution decreases monotonically with increasing pore size. This is markedly different from that in configurations of hard disks (monomeric fluid) where the pore size distribution is peaked at finite size.

Critical dynamics of ballistic and Brownian particles in a heterogeneous environment

The dynamic properties of a classical tracer particle in a random, disordered medium are investigated close to the localization transition. For Lorentz models obeying Newtonian and diffusive motion at the microscale, we have performed large-scale computer simulations, demonstrating that universality holds at long times in the immediate vicinity of the transition. The scaling function describing the crossover from anomalous transport to diffusive motion is found to vary extremely slowly and spans at least five decades in time. To extract the scaling function, one has to allow for the leading universal corrections to scaling. Our findings suggest that apparent power laws with varying exponents generically occur and dominate experimentally accessible time windows as soon as the heterogeneities cover a decade in length scale. We extract the divergent length scales, quantify the spatial heterogeneities in terms of the non-Gaussian parameter, and corroborate our results by a thorough finite-size analysis.

The effect of matrix structure on the diffusion of fluids in porous media

The effect of matrix structure on the transport properties of adsorbed fluids is studied using computer simulations and percolation theory. The model system consists of a fluid of hard spheres diffusing in a matrix of hard spheres fixed in space. Three different arrangements of the fixed spheres, random, templated, and polymeric, are investigated. For a given matrix volume fraction the diffusion coefficient of the fluid, D, is sensitive to the manner in which the matrix is constructed, with large differences between the three types of matrices. The matrix is mapped onto an effective lattice composed of vertices and bonds using a Voronoi tessellation method where the connectivity of bonds is determined using a geometric criterion, i.e., a bond is connected if a fluid particle can pass directly between the two pores the bond connects, and disconnected otherwise. The percolation threshold is then determined from the connectivity of the bonds. D displays universal scaling behavior in the reduced volume fraction, i.e., D approximately (1-phi(m)phi(c))(gamma), where phi(m) is the matrix volume fraction and phi(c) is the matrix volume fraction at the percolation threshold. We find that gamma approximately 2.2, independent of matrix type, which is different from the result gamma approximately 1.53 for diffusion in lattice models, but similar to that for conduction in Swiss cheese models. Lattice simulations with biased hopping probabilities are consistent with the continuous-space simulations, and this shows that the universal behavior of diffusion is sensitive to details of local dynamics.

Void percolation and conduction of overlapping ellipsoids

The void percolation and conduction problems for equisized overlapping ellipsoids of revolution are investigated using the discretization method. The method is validated by comparing the estimated percolation threshold of spheres with the precise result found in literature. The technique is then extended to determine the threshold of void percolation as a function of the geometric aspect ratio of ellipsoidal particles. The finite element method is also applied to evaluate the equivalent conductivity of the void phase in the system. The results confirm that there are no universalities for void percolation threshold and conductivity in particulate systems, and these properties are clearly dependent on the geometrical shape of particles. As a consequence, void percolation and conduction associated with ellipsoidal particles of large aspect ratio should be treated differently from spheres.

Lateral diffusion and percolation in membranes

An algorithm based on Voronoi tessellation and percolation theory is presented to study the diffusion of model membrane components (solutes) in the plasma membrane. The membrane is modeled as a two-dimensional space with integral membrane proteins as static obstacles. The Voronoi diagram consists of vertices, which are equidistant from three matrix obstacles, joined by edges. An edge between two vertices is said to be connected if solute particles can pass directly between the two regions. The percolation threshold, pc, determined using this passage criterion is pc approximately equal to 0.53. This is smaller than if the connectivity of edges were assigned randomly, in which case the percolation threshold pr=2/3, where p is the fraction of connected edges. Molecular dynamics simulations show that diffusion is determined by percolation of clusters of edges.

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

The localization transition and the critical properties of the Lorentz model in three dimensions are investigated by computer simulations. We give a coherent and quantitative explanation of the dynamics in terms of continuum percolation theory and obtain an excellent matching of the critical density and exponents. Within a dynamic scaling ansatz incorporating two divergent length scales we achieve data collapse for the mean-square displacements and identify the leading corrections to scaling. We provide evidence for a divergent non-Gaussian parameter close to the transition.

Model system for classical fluids out of equilibrium

A model system for classical fluids out of equilibrium, referred to as a dissipative particles dynamics (DPD) solid, is studied by analytical and simulation methods. The time evolution of a DPD particle is described by a fluctuating heat equation. This DPD solid with transport based on collisional transfer (high-density mechanism) is complementary to the Lorentz gas with only kinetic transport (low-density mechanism). Combination of both models covers the qualitative behavior of transport properties of classical fluids over the full-density range. The heat diffusivity is calculated using a mean-field theory, leading to a linear-density dependence of this transport coefficient, which is exact at high densities. Subleading density corrections are obtained as well. At lower densities the model has a conductivity threshold below which heat conduction is absent. The observed threshold is explained in terms of percolation diffusion on a random proximity network. The geometrical structure of this network is the same as in continuum percolation of completely overlapping spheres, but the dynamics on this network differs from continuum percolation diffusion. Furthermore, the kinetic theory for DPD is extended to the generalized hydrodynamic regime, where the wave-number-dependent decay rates of the Fourier modes of the energy and temperature fields are calculated.

Precise determination of the void percolation threshold for two distributions of overlapping spheres

The void percolation threshold is calculated for a distribution of overlapping spheres with equal radii, and for a binary-sized distribution of overlapping spheres, where half of the spheres have radii twice as large as the other half. Using systems much larger than previous work, we determine a much more precise value for the percolation thresholds and correlation length exponent. The value of the percolation threshold for the monodisperse case is shown to be 0. 0301+/-0.0003, whereas the value for the bidisperse case is shown to be p(c)=0.0287+/-0.0005. The fact that these are significantly different is in contrast with previous, less precise works that speculated that the threshold might be universal with respect to sphere size distribution.

Voronoi-Delaunay analysis of voids in systems of nonspherical particles

The Voronoi network is known to be a useful tool for the structural description of voids in the packings of spheres produced by computer simulations. In this article we extend the Voronoi-Delaunay analysis to packings of nonspherical convex objects. Main properties of the Voronoi network, which are known for systems of spheres, are valid for systems of any convex objects. A general numerical algorithm for calculation of the Voronoi network in three dimensions is proposed. It is based on the calculation of the trajectory of the imaginary empty sphere of variable size, moving inside a system (the Delaunay empty sphere method). Analysis of voids is presented for an ensemble of random straight lines and for a molecular dynamics model of liquid crystal. The spatial distribution of voids and a simple percolation analysis are obtained. The distributions of the bottleneck radii and the radii of spheres inscribed in the voids are calculated.