Percolation of aligned rigid rods on two-dimensional triangular lattices

Effect of shape anisotropy on percolation of aligned and overlapping objects on lattices

We investigate the percolation transition of aligned, overlapping, anisotropic shapes on lattices. Using the recently proposed lattice version of excluded volume theory, we show that shape-anisotropy leads to some intriguing consequences regarding the percolation behavior of anisotropic shapes. We consider a prototypical anisotropic shape-rectangle-on a square lattice and show that, for rectangles of width unity (sticks), the percolation threshold is a monotonically decreasing function of the stick length, whereas, for rectangles of width greater than two, it is a monotonically increasing function. Interestingly, for rectangles of width two, the percolation threshold is independent of its length. We show that this independence of threshold on the length of a side holds for d-dimensional hypercubiods as well as for specific integer values for the lengths of the remaining sides. The limiting case of the length of the rectangles going to infinity shows that the limiting threshold value is finite and depends upon the width of the rectangle. This "continuum" limit with the lattice spacing tending to zero only along a subset of the possible directions in d dimensions results in a "semicontinuum" percolation system. We show that similar results hold for other anisotropic shapes and lattices in different dimensions. The critical properties of the aligned and overlapping rectangles are evaluated using Monte Carlo simulations. We find that the threshold values given by the lattice-excluded volume theory are in good agreement with the simulation results, especially for larger rectangles. We verify the isotropy of the percolation threshold and also compare our results with models where rectangles of mixed orientation are allowed. Our simulation results show that alignment increases the percolation threshold. The calculation of critical exponents places the model in the standard percolation universality class. Our results show that shape anisotropy of the aligned, overlapping percolating units has a marked influence on the percolation properties, especially when a subset of the dimensions of the percolation units is made to diverge.

Standard and inverse site percolation of triangular tiles on triangular lattices: Isotropic and perfectly oriented deposition and removal

Numerical simulations and finite-size scaling analysis have been carried out to study standard and inverse percolation of triangular tiles of side k (k-tiles) on triangular lattices. In the case of standard percolation, the lattice is initially empty. Then, k-tiles are randomly and sequentially deposited on the lattice. In the case of inverse percolation, the process starts with an initial configuration where all lattice sites are occupied by single monomers (each monomer occupies one lattice site) and, consequently, the opposite sides of the lattice are connected by nearest-neighbor occupied sites. Then, the system is diluted by randomly removing k-tiles [composed by k(k+1)/2 monomers] from the lattice. Two schemes are used for the depositing and removing process: the isotropic scheme, where the deposition (removal) of the objects occurs with the same probability in any lattice direction; and the anisotropic (perfectly oriented or nematic) scheme, where one lattice direction is privileged for depositing (removing) the tiles. The study is conducted by following the behavior of four critical concentrations with the size k: (i) [(ii)] standard isotropic (oriented) percolation threshold θ_{c,k} (ϑ_{c,k}), which represents the minimum concentration of occupied sites at which an infinite cluster of occupied nearest-neighbor sites extends from one side of the system to the other. θ_{c,k} (ϑ_{c,k}) is reached by isotropic (oriented) deposition of k-tiles on an initially empty lattice; and (iii) [(iv)] inverse isotropic (oriented) percolation threshold θ_{c,k}^{i} (ϑ_{c,k}^{i}), which corresponds to the maximum concentration of occupied sites for which connectivity disappears. θ_{c,k}^{i} (ϑ_{c,k}^{i}) is reached after removing isotropic (completely aligned) k-tiles from an initially fully occupied lattice. The obtained results indicate that (1)θ_{c,k} (θ_{c,k}^{i}) is an increasing (decreasing) function of k in the range 1≤k≤6. For k≥7, all jammed configurations are nonpercolating (percolating) states and, consequently, the percolation phase transition disappears. (2)ϑ_{c,k} (ϑ_{c,k}^{i}) show a behavior qualitatively similar to that observed for isotropic deposition. In this case, the minimum value of k at which the phase transition disappears is k=5. (3) For both isotropic and perfectly oriented models, the curves of standard and inverse percolation thresholds are symmetric to each other with respect to the line θ(ϑ)=0.5. Thus, a complementary property is found θ_{c,k}+θ_{c,k}^{i}=1 (and ϑ_{c,k}+ϑ_{c,k}^{i}=1), which has not been observed in other regular lattices. (4) Finally, in all cases, the jamming exponent ν_{j} was measured, being ν_{j}=1 regardless of the orientation (isotropic or nematic) or 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 five (three) for isotropic (nematic) deposition scheme, has the same universality class as the standard percolation problem.

Universal Critical Behavior of Percolation in Orientationally Ordered Janus Particles and Other Anisotropic Systems

We combine percolation theory and Monte Carlo simulation to study in two dimensions the connectivity of an equilibrium lattice model of interacting Janus disks which self-assemble into an orientationally ordered stripe phase at low temperature. As the patch size is increased or the temperature is lowered, clusters of patch-connected disks grow, and a percolating cluster emerges at a threshold. In the stripe phase, the critical clusters extend longer in the direction parallel to the stripes than in the perpendicular direction, and percolation is thus anisotropic. It is found that the critical behavior of percolation in the Janus system is consistent with that of standard isotropic percolation, when an appropriate spatial rescaling is made. The rescaling procedure can be applied to understand other anisotropic systems, such as the percolation of aligned rigid rods and of the q-state Potts model with anisotropic interactions.

Random adsorption process of linear k-mers on square lattices under the Achlioptas process

We study the explosive percolation with k-mer random sequential adsorption (RSA) process. We consider both the Achlioptas process (AP) and the inverse Achlioptas process (IAP), in which giant cluster formation is prohibited and accelerated, respectively. By employing finite-size scaling analysis, we confirm that the percolation transitions are continuous, and thus we calculate the percolation threshold and critical exponents. This allows us to determine the universality class of the k-mer explosive percolation transition. Interestingly, the numerical simulation suggests that the universality class of the explosive percolation transition with the AP alters when the k-mer size changes. In contrast, the universality class of the transition with the IAP is independent of k, but it differs from that of the RSA without the IAP.

Standard and inverse site percolation of straight rigid rods on triangular lattices: Isotropic and perfectly oriented deposition and removal

Numerical simulations and finite-size scaling analysis have been carried out to study standard and inverse percolation of straight rigid rods on triangular lattices. In the case of standard percolation, the lattice is initially empty. Then, linear k-mers (particles occupying k consecutive sites along one of the lattice directions) are randomly and sequentially deposited on the lattice. In the case of inverse percolation, the process starts with an initial configuration where all lattice sites are occupied by single monomers (each monomer occupies one lattice site) and, consequently, the opposite sides of the lattice are connected by nearest-neighbor occupied sites. Then the system is diluted by randomly removing sets of k consecutive monomers (linear k-mers) from the lattice. Two schemes are used for the depositing/removing process: an isotropic scheme, where the deposition (removal) of the linear objects occurs with the same probability in any lattice direction, and an anisotropic (perfectly oriented) scheme, where one lattice direction is privileged for depositing (removing) the particles. The study is conducted by following the behavior of four critical concentrations with size k: (i) [(ii)] standard isotropic[oriented] percolation threshold θ_{c,k}[ϑ_{c,k}], which represents the minimum concentration of occupied sites at which an infinite cluster of occupied nearest-neighbor sites extends from one side of the system to the other. θ_{c,k}[ϑ_{c,k}] is reached by isotropic[oriented] deposition of straight rigid k-mers on an initially empty lattice; and (iii) [(iv)] inverse isotropic[oriented] percolation threshold θ_{c,k}^{i}[ϑ_{c,k}^{i}], which corresponds to the maximum concentration of occupied sites for which connectivity disappears. θ_{c,k}^{i}[ϑ_{c,k}^{i}] is reached after removing isotropic [completely aligned] straight rigid k-mers from an initially fully occupied lattice. θ_{c,k}, ϑ_{c,k}, θ_{c,k}^{i}, and ϑ_{c,k}^{i} are determined for a wide range of k (2≤k≤512). The obtained results indicate that (1)θ_{c,k}[θ_{c,k}^{i}] exhibits a nonmonotonous dependence on the size k. It decreases[increases] for small particle sizes, goes through a minimum[maximum] at around k=11, and finally increases and asymptotically converges towards a definite value for large segments θ_{c,k→∞}=0.500(2) [θ_{c,k→∞}^{i}=0.500(1)]; (2)ϑ_{c,k}[ϑ_{c,k}^{i}] depicts a monotonous behavior in terms of k. It rapidly increases[decreases] for small particle sizes and asymptotically converges towards a definite value for infinitely long k-mers ϑ_{c,k→∞}=0.5334(6) [ϑ_{c,k→∞}^{i}=0.4666(6)]; (3) for both isotropic and perfectly oriented models, the curves of standard and inverse percolation thresholds are symmetric to each other with respect to the line θ(ϑ)=0.5. Thus a complementary property is found θ_{c,k}+θ_{c,k}^{i}=1 (and ϑ_{c,k}+ϑ_{c,k}^{i}=1) which has not been observed in other regular lattices. This condition is analytically validated by using exact enumeration of configurations for small systems, and (4) in all cases, the critical concentration curves divide the θ space in a percolating region and a nonpercolating region. These phases extend to infinity in the space of the parameter k so that the model presents percolation transition for the whole range of k.