arXiv:1801.09970v1 [cond-mat.stat-mech] 30 Jan 2018 iceiainshmshv enue oextrapolate to hand, used other been the have On schemes [2–4]. in- discretization analysis tesselation an Voronoi to amenable volving are spheres placed randomly not elements. often discrete spaces into interstitial partitioned with readily implement more around to are voids methods difficult for infiltration cluster However, grains, overlapping exists. spanning system clus- network a if finding conducting determine to for grains [1] connected of algorithm ters such Hoshen-Kopelman techniques the to amenable as is situation net- in latter conducting The random as work. a into grains, sintered charge particles conducting to metallic contiguous complementary through scenario a transport is particles barrier type. any of polyhedra thresh- oriented percolation randomly calculate to around first olds knowl- the our is To study our orientations. edge polyhe- random with interpenetrating inclusions for dral calcu- thresholds we media, percolation granular late occurring naturally often disorder dodecahedra, in strong present the octahedra, incorporate To cubes, icosahedra). tetrahedra, and Platonic the of (i.e. each solids for barriers ran- polyhedral of placed cal- configurations domly by for shapes thresholds grain percolation of culating continuum percolation angular this examine along to We phenomena geometries crystallites. round polyhedral weathered faceted from shapes of hall- transition. the phase all those continuous with and a transition of percolation flow from marks a fluid transition is not bulk the do allowing marking which grains latter, of The den- assemblies threshold grains. a a of beyond on exist sity networks to void cease system scale densities, form macroscopic low spaces at void networks contiguous spanning in- ion the Whereas or surrounding voids flow clusions. of fluid networks with through grains transport impermeable of comprised sebiso ihysmerclicuin uhas such inclusions symmetrical highly of Assemblies impenetrable around voids through flow ion or Fluid spectrum a occupy media porous up making Grains are materials porous situations, relevant of variety a In eclto hog od rudRnol retdPlaton Oriented Randomly around Voids through Percolation 1 ue,othdda oeaer,adioaer) neach In fo icosahedra). inclusions concentrations exact and polyhedral an dodecahedra, give octahededra, for we cubes, densities sca particles, critical macroscopic tracer on calculate virtual traversed using net be spaces void not void spanning may system which contain configurations which grains percola of A assemblies particles. barrier constituent impenetrable the esstelte nyfrcb-hpdgan.W calculate We value grains. subdiffusive finding cube-shaped threshold, for percolation only latter the versus ASnmes 64.60.ah,61.43.Gt,64.60.F- numbers: PACS eateto hsc srnm,Yugtw tt Univer State Youngstown Astronomy, & Physics of Department oosmtrasmd po memal oyerlgan c grains polyhedral impermeable of up made materials Porous ρ c o lge n admyoine ris nigdistinct finding grains, oriented randomly and aligned for .J ror radN .McGuigan J. N. and Jr Priour, J. D. Dtd coe ,2018) 7, October (Dated: z 0 = ntecs fcbcicuin;fralohrplhdaun- polyhedra other only all the for scenarios consideration inclusions; former der cubic percola- of and distinct case latter the finding the in barriers, for aligned thresholds oriented for tion randomly as are well as polyhedra constituent the which densities critical calculating polyhedra. included of surfaces reflections planar follow- specular the inclusions, by from placed punctuated randomly trajectories of linear center ing sample the at large location a initial of an from tracer propagate implementation, particles plausible computationally of physically a context and For the efficient [10–13]. in spheres described dynamical placed been of randomly have examples 3D Previous in simulations dynamically spaces. to void particles tracer infiltrate virtual barrier of using geometry by the particles of void treatment of exact an description with explicit Toshapes an circumvent [9]. we cubes polyhe- grains, placed oriented aligned dral randomly randomly and for of 8], thresholds percolation [7, case find ellipsoids the 6], in [5, limit spheres continuum the to oeo h ita rcr rvlfreog obreach to enough that far insist travel we tracers size, virtual in the infinite of assemblies effectively For none be to inclusions. grains placed of randomly of volumes n find and is techniques error. two Carlo the Monte the of to results comparable among agreement sity; itclflcutos n osdra es 4 least at consider and fluctuations, tistical itnttm u-nevl samast eemn fone if below determine for or to above means displacements is a is RMS sub-intervals time the distinct from calculated exponents z n ad aaclassi h otx ffiiesz scal- size finite of yield context the analysis the On ing in collapses data methods. hand, complementary one two using thresholds error. statistical Carlo lentvl,cmaigeetv yaia scaling dynamical effective comparing Alternatively, . . 91 omnt l ri hpsconsidered. shapes grain all to common 19(1) eeaieteipc foinainldsre by disorder orientational of impact the examine We nodrt xrplt otetemdnmclimit thermodynamic the to extrapolate to order In si eetsuy[4 ecluaepercolation calculate we [14] study recent a in As h v ltncsld ie tetrahedra, (i.e. solids platonic five the r intasto ak h onaybetween boundary the marks transition tion ae ecluaeproainthreshold percolation calculate we case, ramn fgangoere,adwe and geometries, grain of treatment h yaia cln xoeta the at exponent scaling dynamical the ok,amtigbl rnpr,and transport, bulk admitting works, ρ c e.Wt yaia nlrto of infiltration dynamical With les. eaeaeoe iodrt upessta- suppress to disorder over average we , ntanfli o ovisaround voids to flow fluid onstrain iy onson H455 USA 44555, OH Youngstown, sity, 1 ρ ρ c c n hrb n h hehl den- threshold the find thereby and , n h yaia rtclexponent critical dynamical the and ρ ρ c c ausaeietclu oMonte to up identical are values ausfrteformer the for values ρ c o sebiso risin grains of assemblies for cSolids ic × 10 4 cubic 2 the simulation volume. To this end, we use system sizes ranging from L = 90 to L = 170 in units of the radius of spheres circumscribed about the polyhedra; smaller sizes are appropriate for the more densely packed angu- lar tetrahedral grains, while larger systems are needed for the bulkier dodecahedra and icosahedra. To allow for at least on the order of 106 collisions with faceted inclu- sions, tracer dwell times are 2 105 to 5 105 time units (with normalized inter-collision× velocities× where ~v = 1), depending on the type of barrier particle. | | z To mitigate the computational burden of locating the nearest polyhedral surface, the cubic simulation volumes are partitioned into voxels a unit on a side, with disorder realizations generated voxel by voxel for optimal memory usage. The probability for the geometric centers of n n ρv grains in a voxel, P (n) = (ρv) e− /n! (v = 1 being the voxel volume), may be efficiently and accurately sampled with a robust random number generator. x In identifying collisions of tracer particles with grains, y neighboring voxels are checked for penetration of spheres circumscribed about polyhedra; in the event of contact FIG. 1: (Color online) Configurations of randomly oriented with a sphere, each planar facet of the polyhedron within tetrahedral grains with a subset (i.e. 100) of the tetrahedra is examined for an intersection with the tracer trajectory. shown. The fact that regular convex polyhedra are the set of all points confined within the planes containing the m facets is tantamount to the constraintn ˆ (V~ V~ ) ǫ i · − 0 ≤ with V~0 being the geometric center of the solid and V~ a point within the polyhedron;n ˆi is the unit vector normal to the ith facet, and ǫ the distance of each facet plane to the polyhedron center. For a collision with the ith facet, the relation becomes an equality which may be solved for a transit time, with the polyhedron constraint imposed for each of the m 1 remaining faces to reject z spurious scattering events.− Appealing to the geometries of the corresponding polyhedra gives ǫ = 1/3 for the tetrahedron, ǫ =1/√3 for the cube and octahedron, and ǫ = [(5 + 2√5)/15]1/2 for dodecahedra and icosahedra. The implementation of an aligned or randomly ori- ented scheme for the polyhedral inclusions is determined by the choice of the local coordinate system for each grain. For aligned polyhedra, normalized vectors cor- responding to each planar facet are specified in relation x to a common set of Cartesian axes. On the other hand, y for randomly oriented grains, distinct orthogonal coor- dinate systems are selected randomly for each inclusion. FIG. 2: (Color online) Configurations of randomly oriented To prevent an orientational bias, we isotropically choose icosahedral grains with a subset (i.e. 100) of the icosahedra two normalized vectorsv ˆ1 andv ˆ2. The latter are con- shown. verted into an orthogonal pair by puttingu ˆ1 =v ˆ1, with uˆ obtained by normalizingv ˆ (ˆv vˆ )ˆv to remove 2 2 − 1 · 2 1 the component parallel tov ˆ1. A third perpendicular axis are shown in Fig. 1 for tetrahedra and in Fig. 2 for icosa- is then fixed withu ˆ3 =u ˆ1 uˆ2. The components for hedra. ′ × the vectorsn ˆi normal to the facets of a polyhedron are A trajectory sequence for a tracer particle passing ′ 3 hence nik = Pj=1 nij ujk with j and k corresponding to through void networks among randomly oriented icosahe- the randomly oriented and common coordinate systems, dral grains is shown in Fig. 3, with each blue segment rep- respectively. Assemblies of randomly oriented polyhedra resenting 400 collisions for the sake of clarity in highlight- 3

z(t t ) 4.5 (a) 0.26 16 , 4 z (t4 , t1) 0.22 4.0 0.18

Dynamical Exponent 0.14

RMS 3.5 5.3 5.4 5.5 5.6 δ ρ z −

t 3.0 t1 t7 t13 t2 t8 t14 2.5 t3 t9 t15 z t4 t10 t16 t5 t11 2.0 t6 t12 −0.4 −0.2 0.0 0.2 0.4 y(ρ − ρ ) t c

12.0 0.30 z (t t ) (b) 16 , 4 11.0 0.26 z (t4 ,t1) 10.0 0.22 9.0 0.18 Dynamical Exponent y 0.14 RMS 8.0 δ 1.15 1.17 1.19 1.21

x z ρ − 7.0 t t t t 1 7 13 FIG. 3: (Color online) Trajectories for tracer particles moving t2 t8 t14 6.0 t t t among randomly oriented icosahedral grains are show for ρ 3 9 15 6 t4 t10 t16 c × 5.0 just below ρ for 4 10 time units. Blue segments represent t5 t11 t 400 collision events, and the boxed image in the upper right 4.0 6 t12 hand corner is a magnified view of the indicated section of the void network. −0.08 −0.04 0.0 0.04 y (ρ − ρ t c ) ing void networks. In the red box is an expanded view FIG. 4: (Color online) Data Collapse plots are shown for ran- showing the narrow regions joining larger void spaces. domly oriented tetrahedra and randomly oriented dodecahe- dra in panel (a) and panel (b), respectively (red lines represent Critical indices, such as the threshold density ρc and the dynamical critical exponent z may be calculated us- analytical scaling curves). Legends in the plots show symbols corresponding to the 16 subintervals of the total dwell time. ing the framework of single parameter finite size scaling The insets in each graph show the intersection of effective theory [15]; with systems effectively infinite in size, RMS exponents z(t16,t4) and z(t4,t1). displacements are calculated for a range of finite times. In terms of elapsed times (we consider 16 time subin- tervals in this work) rather than sizes, one anticipates a ln[δRMS(t2)/δRMS(t1)]/ ln(t2/t1). Critical densities ob- scaling form δ (ρ)= tzf[ty(ρ ρ )] [15], where z is the tained from the crossings are in agreement with data RMS − c dynamical critical exponent and y a scaling exponent. collapse results up to Monte Carlo error. z y Plotting t− f with respect to t (ρ ρc) for ρ near ρc in principle yields a data collapse with− Monte Carlo data Grain Type ρc vB φc z y Aligned Tetrahedra 5.47(2) 0.5132 0.0605(6) 0.20(1) .23(1) falling on a single curve. We use the data collapse phe- Rotated Tetrahedra 5.47(2) 0.5132 0.0605(6) 0.19(1) .24(1) nomenon as a quantitative tool by considering a scaling Aligned Cubes 2.122(5) 1.5396 0.0381(3) 0.20(1) .25(1) n j Rotated Cubes 2.012(8) 1.5396 0.0452(6) 0.18(1) .23(1) function g(x)= Pj=0 Aj x , and optimizing with respect to the A coefficients as well as ρ , z, and y using non- Aligned Octahedra 2.402(6) 1.3333 0.0407(3) 0.19(1) .25(1) j c Rotated Octahedra 2.419(9) 1.3333 0.0398(5) 0.19(1) .23(1) linear least squares fitting. The polynomial form with Aligned Dodecahedra 1.198(3) 2.7852 0.0356(3) 0.18(1) .24(1) relatively few coefficients (n = 5 in this work) is war- Rotated Dodecahedra 1.194(4) 2.7852 0.0360(3) 0.18(1) .22(1) ranted by the approximate linearity of f(x) near x = 0. Aligned Icosahedra 1.326(3) 2.5362 0.0346(3) 0.19(1) .24(1) Rotated Icosahedra 1.338(8) 2.5362 0.0336(7) 0.18(1) .23(1) Sample data collapses are shown in panel (a) and panel Spherical Grains 0.837(1) 4.1888 0.0301(1) 0.17(1) .22(1) (b) of Fig. 4 for randomly oriented tetrahedra and dodec- ahedra, respectively. Data collapses yield the percolation TABLE I: Table with critical indices for randomly oriented threshold concentrations ρc to within a few tenths of a and aligned Platonic solids as well as spheres. percent, as well as the dynamical critical exponent z and the scaling exponent y. Results for randomly oriented Platonic solids and the The graph insets show crossings of effective dynam- aligned counterparts are summarized in Table I. In dis- ical exponents z(t ,t ) and z(t ,t ) where z(t ,t ) cussing void percolation phenomena, one often speci- 16 4 4 1 2 1 ≡ 4

c versal dynamical exponent z = 0.19(1), which is also in

φ rotated (collapse) rotated (crossing) accord with a z obtained in the context of a 3D discrete 0.05 aligned (collapse) system [15]. aligned (crossing) Threshold densities are represented graphically in Fig. 5 with φc and ρc data in the main graph and in-

5.0 set, respectively; diamonds or squared represent results c

0.03 ρ for polyhederal inclusions per the main graph legend, 3.0 sphere data while round symbols indicate spherical grains. Apart

Critical Density 1.0 sphere from cube-shaped grains, φc and ρc data match closely. 0.01 data 0.0 0.2 0.4 0.6 0.8 1.0 The threshold density in the inset is monotonic in increas- Fractional Polyhedron Volume ing polyhedral volume fraction, a trend slightly disrupted Critical Excluded Volume 0.0 0.2 0.4 0.6 0.8 1.0 for the excluded volume, where φc decreases until rising Fractional Polyhedron Volume slightly among icosahedra and dodecahedra. In conclusion, we have calculated percolation thresh- FIG. 5: (Color online) The main graphs shows critical ex- olds ρc for randomly oriented interpenetrating polyhedral cluded volumes φc with respect to the polyhedron volume grains for each of the Platonic solids, the first calculation fraction relative to that of the circumscribed sphere, and the of critical densities and indices characterizing critical be- inset graph displays the threshold density ρc versus the poly- hedron volume fraction. The legend corresponds to both the havior for percolation phenomena through voids around main graph and the inset plot. randomly oriented polyhedra. The diminished ρc for ran- domly oriented cubes relative to aligned counterparts is not replicated for the other solids considered. Results for ρcvB fies the excluded volume φc = e− with vB being the dynamical critical exponent z are consistent with a the inclusion volume. Standard geometric arguments z =0.19(1) candidate for a universal value. give in terms of the radius of the circumscribed sphere We acknowledge helpful discussion with Michael 8√3 8√3 4 2√3 √ 1/2 Crescimanno. Calculations in this work have benefitted vT = 27 , vC = 9 , vO = 3 , vI = 9 (5 + 5) , and 2√2 √ from use of the Ohio Supercomputer facility (OSC) [16]. vD = 3 (5+ 5) for the tetrahedron, octahedron, cube, icosahedron, and dodecahedron respectively. 3 Accordingly, the critical density ρc (in units of R− ) and the threshold excluded volumes φc are specified for aligned and randomly oriented grains. One might antic- [1] J. Hoshen and R. Kopelman, Phys. Rev. B, 14, 3438 ipate orientational disorder to be disruptive to void net- (1976). works and thus lower the percolation threshold. However, [2] W. T. Elam, A. R. Kerstein, and J. J. Rehr, Phys. Rev. the results are only partly consistent with this intuition. Lett. 52, 1516 (1984). Introducing random orientations in the case of cube- [3] S. C. van der Marck, Phys. Rev. Lett. 77, 1785 (1996). [4] M. D. Rintoul, Phys. Rev. E 62, 68 (2000). shaped grains reduces ρc by 5% from 2.122(5) to 2.012(8). [5] N. S. Martys, S. Torquato, and D. P. Bentz, Phys. Rev. The φ =0.0381(1) result for assemblies of aligned cubes c E 50, 403 (1994). is nearly within error of the value φc =0.36(1) reported [6] R. S. Maier, D. N. Kroll, H. T. Davis, and R. S. Bernard, previously [9]. Nevertheless, up to Monte Carlo error, we J. Colloid Interface Sci. 217, 341 (1999). do not find distinct percolation thresholds for aligned and [7] Y. B. Yi, Phys. Rev. E 74, 031112 (2006). randomly oriented configurations for non-cubic grains. In [8] Y. B. Yi and K. Esmail, J. Appl. Phys. 111, 124903 spite of the robustness of ρc under these circumstances, (2012). the fact that critical densities even for dodecahedra and [9] Z. Koza, G. Kondrat, and K. Suszczy´nski, J. Stat. Mech.: Th. Exp. 11, P11005 (2014). icosahedra differ significantly from ρ for spherical grains c [10] A. Kammerer, F. H¨ofling, and T. Franosch, EPL 84, indicates that all of the polyhedra we consider are dis- 66002 (2008). tinguished in our results from spherical inclusions. On [11] F. H¨ofling, T. Munk, E. Frey, and T. Franosch, J. Chem. the other hand, the prevalence of adjacent parallel non- Phys. 128, 164517 (2008). penetrating faces in the case of aligned cubes, which is [12] M. Spanner, F. H¨ofling, G. E. Schr¨oder-Turk, K. Mecke, not replicated in aligned non-cubic polyhedra, could ac- and T. Franosch, J. Phys.: Condens. Matter 23, 234120 (2011). count for the greater likelihood of intact void networks 89 than in randomly oriented configurations; alignment of [13] D. J. Priour, Jr., Phys. Rev. E, , 012148 (2014). [14] D. J. Priour and N. J. McGuigan,cond-mat/1712.10241. non-cubic inclusions may be of less benefit in mitigating [15] D. Stauffer and A. Aharony, Introduction to Percolation the truncation of connected void spaces. Theory, 2nd ed. (Taylor and Francis, Bristol, UK, 1992). The critical indices in Table I include dynamical ex- [16] Ohio Supercomputer Center. 1987. Ohio Supercomputer ponents z and scaling exponents y; up to Monte Carlo Center. Columbus, OH: Ohio Supercomputer Center. error the results are compatible with a candidate uni- http://osc.edu/ark:/19495/f5s1ph73.