Effect of dimensionality on the continuum percolation of overlapping hyperspheres and hypercubes. II. Simulation results and analyses

• Published in: The Journal of chemical physics Vol. 137; no. 7; p. 074106
• Main Authors: Torquato, S, Jiao, Y
• Format: Journal Article
• Language: English
• Published: United States 21.08.2012
• ISSN: 1089-7690
• DOI: 10.1063/1.4742750

In the first paper of this series [S. Torquato, J. Chem. Phys. 136, 054106 (2012)], analytical results concerning the continuum percolation of overlapping hyperparticles in d-dimensional Euclidean space R(d) were obtained, including lower bounds on the percolation threshold. In the present investigation, we provide additional analytical results for certain cluster statistics, such as the concentration of k-mers and related quantities, and obtain an upper bound on the percolation threshold η(c). We utilize the tightest lower bound obtained in the first paper to formulate an efficient simulation method, called the rescaled-particle algorithm, to estimate continuum percolation properties across many space dimensions with heretofore unattained accuracy. This simulation procedure is applied to compute the threshold η(c) and associated mean number of overlaps per particle N(c) for both overlapping hyperspheres and oriented hypercubes for 3 ≤ d ≤ 11. These simulations results are compared to corresponding upper and lower bounds on these percolation properties. We find that the bounds converge to one another as the space dimension increases, but the lower bound provides an excellent estimate of η(c) and N(c), even for relatively low dimensions. We confirm a prediction of the first paper in this series that low-dimensional percolation properties encode high-dimensional information. We also show that the concentration of monomers dominate over concentration values for higher order clusters (dimers, trimers, etc.) as the space dimension becomes large. Finally, we provide accurate analytical estimates of the pair connectedness function and blocking function at their contact values for any d as a function of density.