Feliz N M De La Cruz, Pimentel F M L, Félix N De La Cruz, Ramirez-Pastor A J
Instituto de Física (IFIS), Facultad de Ciencias, Universidad Autónoma de Santo Domingo-FONDOCYT, Av. Alma Mater, Santo Domingo 10105, Dominican Republic.
Departamento de Física, Instituto de Física Aplicada (INFAP), Universidad Nacional de San Luis - CONICET, Ejército de Los Andes 950, D5700HHW, San Luis, Argentina.
Phys Rev E. 2024 Mar;109(3-1):034107. doi: 10.1103/PhysRevE.109.034107.
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.
已经进行了数值模拟和有限尺寸标度分析,以研究三角形晶格上边长为k的三角形瓷砖(k - 瓷砖)的标准渗流和逆渗流。在标准渗流的情况下,晶格最初是空的。然后,k - 瓷砖被随机且顺序地沉积在晶格上。在逆渗流的情况下,过程从初始构型开始,其中所有晶格位点都被单个单体占据(每个单体占据一个晶格位点),因此,晶格的相对两侧通过最近邻的占据位点相连。然后,通过从晶格中随机移除k - 瓷砖(由k(k + 1)/2个单体组成)来使系统稀释。沉积和移除过程使用两种方案:各向同性方案,其中物体在任何晶格方向上以相同概率进行沉积(移除);以及各向异性(完全取向或向列)方案,其中一个晶格方向在沉积(移除)瓷砖时具有优先权。该研究通过跟踪四个临界浓度随尺寸k的行为来进行:(i) [ (ii) ] 标准各向同性(取向)渗流阈值θ_{c,k} (ϑ_{c,k}),它表示被占据的最近邻位点的无限簇从系统的一侧延伸到另一侧时被占据位点的最小浓度。θ_{c,k} (ϑ_{c,k}) 通过在最初为空的晶格上各向同性(取向)沉积k - 瓷砖来达到;以及(iii) [ (iv) ] 逆各向同性(取向)渗流阈值θ_{c,k}^{i} (ϑ_{c,k}^{i}),它对应于连通性消失时被占据位点的最大浓度。θ_{c,k}^{i} (ϑ_{c,k}^{i}) 通过从最初完全被占据的晶格中移除各向同性(完全对齐)的k - 瓷砖来达到。获得的结果表明:(1) θ_{c,k} (θ_{c,k}^{i}) 在1≤k≤6范围内是k的递增(递减)函数。对于k≥7,所有堵塞构型都是非渗流(渗流)状态,因此,渗流相变消失。(2) ϑ_{c,k} (ϑ_{c,k}^{i}) 表现出与各向同性沉积观察到的行为在定性上相似。在这种情况下,相变消失时k的最小值为k = 5。(3) 对于各向同性和完全取向模型,标准和逆渗流阈值曲线相对于θ(ϑ)=0.5线彼此对称。因此,发现了一个互补性质θ_{c,k}+θ_{c,k}^{i}=1(以及ϑ_{c,k}+ϑ_{c,k}^{i}=1),这在其他规则晶格中未被观察到。(4) 最后,在所有情况下,测量了堵塞指数ν_{j},无论考虑的取向(各向同性或向列)或尺寸k如何,ν_{j}=1。此外,对临界指数ν、β和γ的精确确定表明,对于各向同性(向列)沉积方案,当k在一到五(三)之间变化时系统中发生的渗流相变与标准渗流问题具有相同的普适类。