Connecting Complex Electronic Pattern Formation to Critical Exponents

Scanning probes reveal complex, inhomogeneous patterns on the surface of many condensed matter systems. In some cases, the patterns form self-similar, fractal geometric clusters. In this paper, we advance the theory of criticality as it pertains to those geometric clusters (defined as connected sets of nearest-neighbor aligned spins) in the context of Ising models. We show how data from surface probes can be used to distinguish whether electronic patterns observed at the surface of a material are confined to the surface, or whether the patterns originate in the bulk. Whereas thermodynamic critical exponents are derived from the behavior of Fortuin-Kasteleyn (FK) clusters, critical exponents can be similarly defined for geometric clusters. We find that these geometric critical exponents are not only distinct numerically from the thermodynamic and uncorrelated percolation exponents, but that they separately satisfy scaling relations at the critical fixed points discussed in the text. We furthermore find that the two-dimensional (2D) cross-sections of geometric clusters in the three-dimensional (3D) Ising model display critical scaling behavior at the bulk phase transition temperature. In particular, we show that when considered on a 2D slice of a 3D system, the pair connectivity function familiar from percolation theory displays more robust critical behavior than the spin-spin correlation function, and we calculate the corresponding critical exponent. We discuss the implications of these two distinct length scales in Ising models. We also calculate the pair connectivity exponent in the clean 2D case. These results extend the theory of geometric criticality in the clean Ising universality classes, and facilitate the broad application of geometric cluster analysis techniques to maximize the information that can be extracted from scanning image probe data in condensed matter systems.


I. INTRODUCTION
Since their invention in 1982, scanning probes have revolutionized our understanding of materials and their surfaces, yielding an ever increasing wealth of data available for in-depth analysis, on a wide variety of systems [1].Theory has been sprinting to catch up in order to interpret it all, and also to discover new ways of extracting information from the data, in order to fully realize this promise of new knowledge from the increasing variety of scanning probes and their ever-increasing experimental capabilities.To date, the majority of theoretical treatments have focused on microscopic physics [2], with few theoretical treatments offering guidance for how to interpret the wealth of information available in the multiscale pattern formation often observed on surfaces.We have recently pioneered a new set of techniques for analyzing scanning probe data by mapping two-component image data to random Ising models [3], based on geometric cluster methods imported from disordered statistical mechanics.A geometric cluster is defined as a set of aligned nearest neighbor sites.The key insight is that near criticality, the spatial configurations of geometric clusters are controlled by the critical fixed point, and therefore the geometric properties encode critical exponents.The method is capable of extracting information from the data about disorder, interactions, and dimension.We have already successfully applied this new technique to * liu305@purdue.eduuncover a unification of the fundamental physics governing the multiscale pattern formation observed in two disparate strongly correlated electronic materials (cuprate superconductors [3,4] and vanadium dioxide [5]).
However, because it is the Fortuin-Kasteleyn (FK) clusters [6] which encode thermodynamic criticality [7], rather than the geometric clusters which are directly accessible experimentally via scanning probes and which we employ in our method, little is known about the general theoretical structure of geometric clusters in random Ising models, and the critical exponents associated with the geometric clusters are unknown for many of the fixed points which are key to interpreting experimental data.In this paper, we advance the theory of criticality as it pertains to geometric clusters in clean Ising models to further develop geometric cluster analysis techniques [3], in order to maximize the information that can be extracted from experiments using these new methods.Although the geometric clusters do not encode thermodynamic criticality, we conjecture that when the geometric clusters percolate, whether at or below the thermodynamic critical temperature, the geometric clusters do encode geometric criticality, complete with its own set of critical exponents, which we further conjecture are distinct from the exponents of uncorrelated percolation when arising in the context of an interacting model.This paper is organized as follows.We first describe in Sec.II the model under consideration.We then ask in Sec.III whether the pair connectivity function can be a power law in any system other than uncorrelated percolation (to which we will answer "yes").In Sec.IV, we present some conjectures considering geometric criti-cality as distinct from thermodynamic criticality.Next, we present in Sec.V our results for the critical geometric cluster exponents on 2D slices of the clean 3D Ising model (denoted by C-3Dx, where x means cross-section), and in Sec.VI our simulations of the pair connectivity function of the clean 2D Ising model (denoted by C-2D), In Sec.VII we discuss the relation of the order parameter in a 3D ferromagnetic Ising model to the percolation of geometric clusters on 2D slices of the 3D system and derive that the position of the 2D slice geometric cluster percolation point coincides with bulk critical point.In Sec.VIII we discuss the implications of the pair connectivity length scale in Ising models.Finally, in Sec.IX, we discuss and summarize our findings about geometric criticality, and present conclusions in Sec.X.

II. THE MODEL
In strongly correlated electronic systems, the combination of disorder and strong correlations can drive complex pattern formation [8], but disentangling correlations from disorder in the experimental system is an open problem.Recently, we have developed new cluster analysis techniques for interpreting scanning image probe data in cases where the spatial data can be abstracted to two components (and thus may be mapped to an Ising variable), and where the resulting cluster patterns display structure on multiple length scales [3].
Scanning probe experiments often reveal complex pattern formation at the surface of strongly correlated electronic systems [8,9].For example, charge stripe orientations display complex geometric patterns at the surface of some cuprate superconductors, as revealed by scanning tunneling microscopy [3,10,11].Complex patterns have also been observed in thin films of VO 2 as it transitions from metal to insulator, via scanning near-field infrared microscopy [5,12].Using the cluster analysis techniques we developed, we showed that in both cases the dominant type of disorder driving the pattern formation is in the random field universality class, and we argued that this is the origin of nonequilibrium behavior in both systems [3][4][5].
In this paper, we expand the diagnostic toolset of the geometric cluster analyses we have developed to include the functional form and qualitative behavior of the connectivity function, derived directly from geometric clusters.We find that the pair connectivity function can be a power law at more than just uncorrelated percolation points, and we compute the corresponding critical exponent via Monte Carlo simulations.
We consider a general short-range Ising model: where the sum runs over the sites of a cubic lattice, chosen with spacing at least as small as the resolution of the images to be studied.The tendency for neighboring regions in the data image to be of like character is modeled as a nearest neighbor ferromagnetic interaction.J sets the overall strength of the in-plane ferromagnetic coupling between nearest-neighbor Ising variables, and J ⊥ represents the overall coupling strength between Ising variables in neighboring planes.The field h represents a generalized external field which couples with the local Ising variables.Microscopic models of non-frustrated disorder flow to two broad classes of disorder under the renormalization group transformation: random bond disorder in the coupling strength (through the term δJ ij ), also known as random T c disorder, and random field disorder couplng with the Ising variable (through the term h i ) [13].
The Ising variable σ i = ±1 on each coarse-grained site represents one of the two possible states, such as the two possible electron nematic orientations in the cuprates, or the two states of conductivity (metallic or insulating) in VO 2 .For the two examples above, our cluster analysis [3][4][5] of the image data shows that the geometric clusters (which are defined as the connected set of the nearest neighbor sites with Ising variables being the same value) extracted from the multi-scale pattern formation display universal scaling behavior over multiple decades, suggesting criticality and universality as the origin of the spatial complexity revealed by scanning probe microscopy in strongly correlated electronic systems.By comparing the data-extracted critical exponents derived from the self-similarity of the geometric clusters with the theoretical values for the fixed points contained in Eqn. 1, we have shown [3] that it is possible to identify the universality class governing the scaling behavior of the geometric clusters observed on surfaces of novel materials.
Once the universality class driving the pattern formation has been identified, this in turn yields information about the relative importance of disorder and interactions, the dominant type of disorder in the system, and the dimension of the phenomenon studied [3][4][5]11].For example, by extracting the critical exponents associated with the geometric clusters appearing at the surface of a material, it is possible to understand whether those clusters are forming merely due to surface effects, or whether the clusters form throughout the bulk of the material and then intersect the surface, because the critical exponents are sensitive to dimension.

III. WHERE IS THE CONNECTIVITY FUNCTION POWER LAW?
Although geometric clusters do not encode thermodynamic criticality, we claim that whenever the geometric clusters percolate, the resulting power law behavior encodes a type of geometric criticality [14,15].We further claim that such geometric critical points in interacting models are new, distinct fixed points, which are different from uncorrelated percolation, as supported by our results to be presented in this paper.The pair connectivity function, derived from geometric clusters, is known to be power law at uncorrelated percolation points.The pair connectivity function g conn (r) is defined as the probability that two sites separated by a distance r belong to the same connected finite cluster.Can it also be power law at certain fixed points in interacting models?The most natural clusters to define in Ising models are the geometric clusters, defined by nearest-neighbor sets of like Ising "spins."The geometric clusters have the advantage that they are directly accessible in image probe data in a 2component system, and they are the clusters we employ in our method.
However, geometric clusters in Ising models are still poorly understood theoretically, presumably because they do not necessarily percolate at the thermodynamic transition temperature, T c .Rather, it is the FK clusters which percolate at the thermodynamic transition temperature T c , and which encode the thermodynamic critical behavior [7].The FK clusters also constitute the critical "droplets" of the Fisher droplet model [16].In the case of the clean 2D Ising model (C-2D), it has been shown that the geometric clusters also percolate at the thermodynamic transition, where T p = T c [14].In other cases, such as the clean and random bond Ising models in three dimensions, it has been shown that the bulk geometric clusters percolate at a temperature T p < T c , inside the ordered phase [14,15,17,18].
For the thermodynamic fixed points of Eqn. 1 at which T p = T c (as happens at C-2D), the geometric clusters have well-defined fractal dimensions [19] both fractal volume dimension d v and fractal hull dimension d h .In addition, geometric clusters defined on a 2D slice of the clean 3D Ising model also display fractal behavior [20] At these points, because of their fractal structure, the large scale geometric clusters are self-similar, i.e. they look the same on all length scales.This is a tell-tale characteristic of power law behavior, which brings us leading to our first conjecture: Conjecture #1: The connectivity function is a power law at all critical points points for which the geometric clusters have fractal dimensions, d v and d h .
Our results in Sec.VI indicate that indeed the connectivity function is a power law at the C-2D critical fixed point, as shown in Figs. 6 and 7, in support of Conjecture #1.This is of course implicit in the pioneering work of Coniglio and coworkers on the percolation of geometric clusters in this model [14,21].Our results in Sec.V show that the connectivity function is also a power law on a 2D slice as the clean 3D Ising system passes through thermodynamic criticality at T 3D c .(See Figs. 4 and 5.) It was already known from the work of Ref. [20] that geometric clusters defined on a 2D slice have well-defined fractal dimensions at T 3D c in the clean 3D Ising model, lending further support to Conjecture #1.

IV. GEOMETRIC CRITICALITY AND THERMODYNAMIC CRITICALITY: TWO TYPES OF CRITICAL EXPONENTS
It is known that at the clean 2D Ising (C-2D) fixed point, there are two distinct sets of critical exponents: one set associated with the thermodynamic criticality and encoded by the FK clusters, and a different set associated with "geometric criticality", encoded by the geometric clusters [7,14,19].Because both types of clusters percolate right at the transition temperature, T p = T c , both types of clusters exhibit power law behavior, yielding two distinct sets of critical exponents, one set derived from each type of cluster.While the scaling laws governing the exponent relations at this fixed point are the same for the two types of clusters, the values of the geometric exponents are not equal to the values of the thermodynamic exponents, and neither set are equal to the values of uncorrelated percolation exponents.In this sense, the C-2D fixed point is therefore both a thermodynamic critical point and a type of geometric critical point.This also implies that there are two order parameters for the 2D clean Ising model (the magnetization and also the infinite cluster); as we will show below, corresponding to the two order parameters are two distinct length scales.
The set of critical exponents derived from the thermodynamic order parameter (in this case, the magnetization), is associated with percolation of the FK clusters, and encodes thermodynamic criticality -we will denote these exponents by the subscript "c" (to match the "c" in the thermodynamic critical temperature T c .)This type of exponent has been widely discussed in the literature.The other set of exponents is derived from the percolation order parameter, which is the infinite network strength, defined as the ratio of the number of sites in the infinite connected geometric cluster to the total number of sites in the system.This second set of critical exponents is associated with percolation of geometric clusters, and encodes geometric criticality -we will denote these exponents by the subscript "p" (for percolation).The fact that there can be two distinct sets of exponents has been largely ignored in the literature, except at the C-2D [14,19] fixed point where T p = T c , and in the clean 3D Ising model [17,22] where T p < T c .
There are two cases within Eqn. 1 for which T p < T c : the clean and random bond 3D Ising models [14,15,17,18].In these two cases, we claim that as the geometric clusters percolate at T p , they display a type of criticality, even though the system is not at the thermodynamic critical point.To understand why, it is first helpful to have an intuition about how T p in an interacting model is related to uncorrelated percolation.The relation of the percolation temperature T p to T c in random Ising models is constrained by the site percolation thresholds in the uncorrelated case: on a square lattice, the percolation threshold is p c = 0.59, while on a cubic lattice, it is p c = 0.31, where p is the fraction of sites that are occupied.(In the Ising case, p becomes the fraction of sites with aligned spins [14].) The high temperature limit of Eqn. 1 maps to uncorrelated percolation with p = 0.5.So, in 2D, neither up nor down spins percolate at T → ∞, since p(T → ∞) = 0.5 < p c = 0.59 on a square lattice.Since a majority geometric cluster spans the system in the ordered phase, it must pass through a percolation point at T p = T c if the transition is continuous.Indeed, it has been proven rigorously that T p = T c in clean 2D Ising models [14].Consistent with the argument above, it is known that ∆ p = ∆ c in the 2D random field Ising model [23] where ∆ is the disorder strength.We expect similar physics to obtain in the random bond Ising model in two dimensions.
In three dimensions, the high temperature limit of Eqn. 1 still maps to uncorrelated site percolation with p = 0.5, but now both up and down spins span the system, since p(T → ∞) = 0.5 > p c = 0.31 on a cubic lattice, where p c is the site percolation threshold.Upon decreasing temperature, the majority clusters grow below T c , but since they were already spanning the system, the majority clusters do not pass through a percolation point.However, the minority geometric clusters must pass through a percolation point as they begin to shrink below T c , although there is no obvious reason why they should do so at T c .In fact, in the clean 3D case, it has been rigorously shown that T p < T c [17].This inequality also holds in the random bond case as well [18].
At C-2D, the two sets of critical exponents separately satisfy the same scaling relations.For example, as discussed in Refs.[14] and [19] Although we have not found an explicit, direct calculation in the literature of the exponent η p at C-2D, which in this case should be derived from the pair connectivity function familiar from percolation theory, our results based on numerical simulations of the connectivity function g conn (r) at C-2D indicate that d − 2 + η p = 0.104 ± 0.002 as shown in Figs. 6, consistent within error bars with the scaling relation d − 2 + η p = 2β p /ν p (0.104 ± 0.002 ≈ 2 * (5/96)).Thus we see that at C-2D, there are two distinct sets of exponents, and that the geometric exponents also satisfy the scaling relations of criticality.This brings us to our second conjecture: Conjecture #2:: There are two distinct sets of critical exponents which both satisfy scaling relations at all critical points for which T p = T c : one derived from the FK clusters, and the other derived from the geometric clusters.

V. NUMERICAL RESULTS FOR THE CRITICAL EXPONENTS DEFINED ON 2D SLICES OF A 3D SYSTEM
We are ultimately interested in understanding the origin of the complex, 2-component pattern formation observed at the surfaces of some strongly correlated electronic systems [3,5,8] via scanning image probes.In cases where the data show multiscale clusters of two distinct types, one of the drivers of such behavior can be proximity to a critical point of an Ising model, Eqn 1.We focus in this paper on the clean Ising case; future work will include the effects of disorder.The problem at hand, then, is to map the properties of the observed geometric clusters to critical exponents of the model.We focus here on five exponents that that can be directly extracted from the microscopy imaging data, and can be used to compare with different models to identify the universality class (explained more fully below): The Fisher exponent τ which describes the cluster size distribution; the volume fractal dimension d v and the hull fractal dimension d h of geometric clusters; the anomalous dimension η || c which controls the spin-spin correlation function; and as we will see below, η || p which controls the pair connectivity function.
In this paper we keep using Potts model with q = 2, since a Potts model with 2 states is an Ising model with half the critical temperature [26].We are interested in this section in the critical behavior displayed by the 2D slice geometric clusters near the 3D clean Ising phase transition (C-3Dx) which occurs at the bulk 3D transition temperature T 3D c .(In Sec VII we will demonstrate that this temperature indeed coincides with the 2D slice percolation point T slice p .)The simulations of the clean Ising system are done near T 3D c for the study of geometric criticality on 2D slices, and the Wolff single cluster algorithm [27] is known to do a good job in mitigating critical slowing down effects near criticality.Therefore, we use this well-established cluster updating algorithm in our simulations.We adopt open boundary condition, because this representation is consistent with our context of microscopy experiments on real materials with finite scale.Throughout the paper, we use the logarithmic binning method which is a standard technique for analyzing the scaling behavior [28].For the 3D cubic lattice, we use T 3D c = 2.25576393 [29,30], and all simulations are performed using system sizes from L 3 = 80 3 to L 3 = 192 3 averaged over a large number of configurations in the Monte Carlo simulation in the order of magnitude 10 4 .We take samples of configurations every 1000 Wolff steps after the system reaches equilibrium, in order to mitigate explicit spatial correlations between different images.Because one whole cluster is updated for each Wolff iteration, 1000 steps are empirically a large enough interval to make two neighboring sample images look totally different.
In the simulations, we have found that the percolation critical exponents extracted from (+)-clusters and (−)clusters respectively are consistent with each other.This is expected, since at T slice p = T 3D c there is (+)/(−) symmetry (see Section.VII).Therefore we present the results with both clusters considered together for better statistics.
The cluster size distribution D(s) measures the number of geometric clusters with a given size s, and scales with s −τ as geometric clusters become critical.It is important to remember that when considered in the bulk context, geometric clusters do not display criticality at the bulk transition temperature T 3D c , but rather the minority clusters experience a percolation transition (and therefore display criticality) at a lower temperature T 3D p < T 3D c [31].However, we find that when a 2D slice of the system is considered, the geometric clusters as defined on the slice do display critical behavior in the form of power law scaling at the bulk transition temperature T 3D c .For a 2D slice embedded in the 3D bulk system, we find that D(s) ∼ s −τ at T 3D c , with the Fisher exponent τ specific for geometric clusters defined on a 2D cross section of a 3D system.With a finite field of view (FOV), D(s) in the bulk is known to have a scaling bump, which skews τ to a lower value [32].In order to mitigate this effect, we analyze D(s) using only internal clusters, i.e. those which do not touch a boundary.D(s) of the internal clusters shows unskewed power law behavior within a cutoff [33].
Figure 1 shows the cluster size distribution D(s) of internal 2D cross-sectional clusters for different system sizes at T slice p = T 3D c .For each system size, we exclude the last two points, which deviate from power law scaling, and fit all the other points within this cutoff [33].Based on the observed power law behaviors in Fig. 1, we use the discrete logarithmic derivative method (DLD) [34] to extract τ .Since we are only interested in the slope −τ , we normalize all the D(s) with D(s = 1) = 1 for better comparison between different system sizes.As can be seen in the figure, depending on the system size, power law scaling persists for about 2.5 decades of scaling.We extrapolate τ extracted at T 3D c to L → ∞, and we find that a linear fit of τ (L) yields τ = 2.001 ± 0.013 as L → ∞.Here, the extrapolation follows the linear re- gression method with error-in-variables [35].In general, the critical exponents cannot be directly obtained from a rigorous parametric fitting because the exact analytical form of the finite scaling functions are usually unknown.Under this case, linear extrapolation is a standard firstorder approximation method for exponents extractions [36], and we use this for critical exponents extrapolations throughout the whole paper.The natural finitesize scaling hypothesis for the cluster size distribution at criticality reads D(s, L) = (L 2 ) −τ D(s/L 2 ), where D is a universal scaling function [36].The inset (b) in Fig. 1 shows that using the exponent τ derived as above, the results from the main panel exhibit scaling collapse to a universal scaling function D as expected.
The cluster volumes s and hulls h become fractal near a percolation point.The fractal nature of the 2D crosssectional cluster volumes at T slice p = T 3D c can be described by s ∼ R dv s with s ≫ 1.Here, d v is the volume fractal dimension, and R s is the radius of gyration of the cluster, defined by 2R 2 s = i,j |r i − r j | 2 /s 2 , where the sum is over sites in the cluster, and r i and r j are the positions of the ith and jth sites [37].Figure 2 shows a DLD power law fit of s ∼ R dv s , using only the internal clusters to mitigate boundary effects.Since the scaling relation is for s ≫ 1, we exclude the first 3 points which correspond to the non-universal short length regime.All the other points belong to the power law scaling regime.The power law fits at T slice p = T 3D c persist for about two decades of scaling, depending on the system size.A straightforward fit extrapolating d v to the thermodynamic limit yields d v = 1.856 ± 0.018.The inset Fig. 2(b) illustrates the scaling collapse for different system sizes to the universal function ŝ using the extrapolated values of d v , based on the finite-size scaling hypothesis s = L dv ŝ(R s /L) [36].
The cluster hulls h scale with h ∼ R d h h near a perco- lation point.In this case, we are only tracking the outer (externally accessible) surfaces of each cluster.Thus one surface might contain not only the cluster itself, but also the subclusters inside.Therefore, R h here refers to the radius of gyration of all the sites enclosed by the hull, including any subclusters.To mitigate the boundary effect, we still use only the internal clusters in the DLD fit.We exclude the first 2 point from the fit (h is a smaller measure than s for a cluster), which deviate from the power law scaling regime and correspond to non-universal short distance physics.As shown in Fig. 3, the power law behavior extends over about 2 decades, depending on the system size.A straightforward fit extrapolating to the thermodynamic limit at T slice p = T 3D c yields d h = 1.714 ± 0.022.The inset Fig. 3(b) shows the scaling collapse for different system sizes to the universal function ĥ using the extrapolated values of d h , based on the finite-size scaling hypothesis h = L d h ĥ(R h /L) [36].
As introduced in the previous section, the pair connectivity function g conn (r) is defined as the probability that two sites separated by a distance r belong to the same finite cluster.This correlation function scales with g conn (r) ∼ r −(d−2+ηp) close to the percolation point.Fig. 4 shows g conn (r) at T slice p = T 3D c with different system sizes.As shown by the figure, in addition to the power law behavior, associated with g conn (r) there appears to be an exponential decay characterized by the correlation length scale of the finite system.So we infer the scaling form g conn (r) ∼ r −(d−2+ηp) e −r/ξp with a finite correlation length, and then fit the curves in the main panel using this scaling form.From the fits, we find that the curves of g conn (r) is consistent with the shape of scaling form with some large ξ p comparable with the finite system size.The inset (a) of Fig. 4 shows the extracted anomalous dimension d − 2 + η p from the scaling  In Fig. 5, we show the C-3Dx pair connectivity function g conn (r) simulated with system size L = 192 for different temperatures near T 3D c .It is evident that for T T 3D c , g conn (r) is consistent with the shape of the scaling form g conn (r) ∼ r −(d−2+ηp) e −r/ξp ; and for T < T 3D c , g conn (r) turns up at large r.For g conn (r) calculated using all the possible pairs in the system, this up-turn is a characteristic behavior of the connectivity function when the system transitions from the "disordered side" (zero infinite network strength P = 0) to the "ordered side" (non zero infinite network strength P > 0) through a percolation point, indicating the presence of an infinite cluster on a 2D slice.(Here P is the percolation order parameter defined on the 2D slice.)Therefore, this behavior can be used as a diagnostic tool to estimate the percolation point in an experimental situation.In the present case, this up-turn happens in the vicinity of T 3D c , corroborating the idea that T slice p coincides with the bulk thermodynamic transition temperature, T slice p = T 3D c .We have studied the scaling behavior of cross-sectional 2D slice geometric clusters embedded in the 3D bulk clean Ising system.We have found robust power law scaling (about two decades) for the measures related to the self-similarity of geometric clusters at T 3D c , as well as the existence of universal scaling functions.All of these phenomena corroborate the conjecture that geometric criticality on a 2D slice occurs at the bulk critical temperature T slice p = T 3D c , and we have extracted the corresponding critical cluster exponents, as summarized in Table I.

VI. NUMERICAL RESULTS FOR THE CONNECTIVITY FUNCTION IN THE CLEAN 2D ISING MODEL
In the clean 2D Ising model, the percolation point coincides with the critical point, and both up and down geometric clusters percolate symmetrically at T 2D p = T 2D c under zero external field [14,38].The values of the exponents τ , d v and d h are already given in the work by Janke et.al. [19], as summarized in Table I.In the present work, we simulate the connectivity function in order to obtain the value of the anomalous dimension 2) for the square lattice [26].The inset (a) shows d − 2 + η p extracted from the scaling form fit of the connectivity function g conn (r) ∼ r −(d−2+ηp) e −r/ξp at T 2D c for different system sizes, and a straightforward fit extrapolating to the thermodynamic limit gives d − 2 + η p = 0.104 ± 0.002.With value, all of the curves in the main panel collapse on top of each other, as shown by the inset (b). Figure 7 shows the behavior of g conn (r) around T 2D c with system size fixed at L = 256.As for the case of C-3Dx, we find that for T > T 2D c , the pair connectivity function g conn (r) is consistent with the scaling form shape and for T < T 2D c g conn (r) turns up.This change of behavior happens in the close vicinity of T 2D c , consistent with the idea that

VII. RELATION OF INFINITE CLUSTER TO BULK MAGNETIZATION ON A 2D SLICE
The bulk minority geometric clusters for the clean 3D Ising model are known to percolate inside the ordered phase, with the percolation temperature [17,22].Therefore the bulk geometric clusters do not display fractal properties at the critical temperature T c , although they do at T 3D p .At first glance, one might think that the geometric clusters on a 2D slice, which are cross-sections of the 3D bulk clusters, should also display critical scaling behavior at the bulk percolation temperature T 3D p .However, as shown by Fig. 8(a), the geometric clusters defined on a 2D cross-section at T 3D p do not display self-similar scaling behavior, since the minority clusters are small, and a single majority cluster spans the slice.This shows that T 3D p is not the percolation point for 2D slice geometric clusters.On the other hand, inspired by the fact that the percolation point coincides with the critical point for 2D Ising model [14], we find that at T 3D c , as shown by Fig. 8(b), the 2D crosssectional geometric clusters indeed show self-similarity and fractal behavior over all lengthscales in the field of view.This is consistent with the findings of Ref. 20, which argues that T slice p = T 3D c .Also, our numerical simulation in Section V corroborates T slice p = T 3D c from various perspectives, including robust power law behaviors of 2D slice geometric clusters at T 3D c and change of behavior of pair-connectivity function through T 3D c .In the clean 2D Ising model, Coniglio and co-workers [14] showed that the percolation temperature coincides with the thermodynamic transition temperature, T p = T c , by showing that M < R, where R is the infinite network strength, and M is the net magnetization.This comes from the very physical idea that in a 2D system with thermal fluctuations, a net magnetization requires that at least one geometric cluster spans the system.From Fig. 8, it appears that something similar must be going on in the 3D Ising model, when clusters are defined an a 2D slice.As will be discussed later in this section, our numerical results indeed show that M < R slice , where R slice is the infinite network strength defined on a plane of a 3D system [39].
To rigorously prove that T slice p = T 3D c would require a proof the following lemma: For a 2D slice of the 3D ferromagnetic Ising model with nearest-neighbor interactions at zero external field (H = 0) and below the critical temperature(T < T c ), we have The symbol 0 + (0 − ) means H = 0 with (+)-boundary ((−)-boundary) condition.For our discussion, we limit ourselves to the (+)-boundary by convention, and the results of (−)-boundary follow those of (+)-boundary under symmetry.M (0 + , T ) is the reduced spontaneous magnetization on the 2D slice (which is identical to the reduced spontaneous magnetization of the bulk system at thermodynamic limit), and is given by M (0 + , T ) = π0 + − π 0 + .Here π i = 1 2 (1 − σ i ) and πi = 1 2 (1 + σ i ) are the lattice gas variables relative to the ith site, and σ i = ±1 is the usual Ising spin variable.πi (π i ) is equivalently the characteristic function of the event that spin in i is up (down).The characterstic function of an event is that when the event happens, the function takes the value 1, otherwise it takes the value 0. The subscript 0 refers to the origin, and + means the thermal average for the Ising model with (+)-boundary.Here we only take one slice for each Ising configuration to be counted in the thermal average, and we do not average over slices belonging to the same configuration in order to reduce spatial dependence.Therefore there is a one-to-one correspondence between the slice configuration and bulk system configuration, and by convention we fix the slice to be the z = 0 plane in the Z 3 cubic lattice so that the slice contains the origin.R slice is the 2D slice infinite network strength for up(down) spins, defined as the probability that the origin belongs to the infinite 2D slice (+)-cluster ((−)-cluster), or equivalently it is the weight of the infinite 2D slice (+)-cluster((−)-cluster).Mathematically, it is given by R slice is the characteristic function of the event that spin in i belongs to the infinite 2D slice (+)-cluster((−)-cluster).
For the 2D slice of a cubic lattice system, which is a simple planar graph admitting an elementary cell and two axes of symmetry, we have the statement that (i) an infinite 2D slice (+)-cluster cannot coexist with an infinite 2D slice (−)-cluster, and (ii) the number of each type of infinite 2D slice cluster must be either 0 or 1 [14,38,[40][41][42].This means that and thus there can only be at most one infinite (+)cluster on the 2D slice with the (+) boundary condition.Therefore, Eqn. ( 2) equivalently becomes as intuitively discussed in the previous paragraph.The rigorous demonstration of (4) requires that both the Markov property [43] and the Fortuin-Kasteleyn-Ginibre (FKG) inequality [44] both hold [14].The FKG inequality [45] may be directly extended to the 2D slice since it is applicable to any finite set in the bulk system.However the Markov property cannot be restricted to a 2D slice, since the Markov chain may extend off of the slice in a 3D system.Rather than delving into a mathematical physics approach, for the purpose of this paper, we numerically demonstrate the fact that the inequality (4) holds for the 2D slice embedded in the 3D Ising model, as revealed by Fig. 9.
By symmetry we have R slice 3) and ( 4) yield R slice ↓ (0 + , T ) = 0 and R slice ↑ (0 + , T ) ≥ M (0 + , T ) > 0. If we suppose that the finite cluster weight is continuous and has a limit at T 3D c , then due to symmetry we have lim we have lim T →T 3D c γ∞ 0 + = lim T →T 3D c γ ∞ 0 + , which then leads to lim T →T 3D c R slice ↑ (0 + , T ) = 0, so that no jump occurs at T 3D c .This is intuitively consistent with a continuous phase transition.In other words, under the (+) boundary condition with zero external field, the 2D slice (+)-cluster percolates at T 3D c , as characterized by the percolation order parameter R slice ↑ (0 + , T ) which continuously changes from 0 to non-zero.
All of the above statements apply to the (−)-cluster with (−) boundary condition under symmetry.Thus, for zero applied field H = 0, the (+)-cluster and (−)-cluster defined on a 2D slice both approach the percolation point at the limit T → T 3D c .The spontaneous infinite cluster type in the ordered phase T < T 3D c is decided by the boundary condition (whether (+) or (−)), and is symmetric under H = 0 + ↔ H = 0 − .Therefore with H = 0 at the limit T → T 3D c , both types of clusters percolate on the slice, and due to symmetry they present the same universal scaling behavior with the same critical exponents.
As supported by our numerical simulation as summarized in Table I, this critical percolation point of geometric clusters defined on 2D slices of the bulk system in c , averaged over a number of configurations in the order of magnitude ∼ 10 5 .M and R both approach 1 when T decreases, with M < R.
the 3D clean Ising model is a new universality class of geometric criticality, distinct from that of the clean 2D model and uncorrelated percolation in 2D, as well as that of the bulk percolation which occurs inside the ordered phase of the 3D system, since T 3D p = T slice p .

VIII. TWO CORRELATION LENGTHS
In this section we elaborate on the two different correlation lengths associated with the pair connectivity function and the spin-spin correlation function.In Fig. 10 we illustrate the pair connectivity function g conn and spinspin correlation function g spin for both the 2D slice of the 3D clean Ising model and 2D clean Ising model at their respective transition temperatures with specific system sizes (L = 160 for C-3Dx and L = 256 for C-2D).Since g spin approaches 0 for large r, it can fluctuate into negative values under our simulation with finite system sizes and finite number of configurations for averaging.Once negative values appear, in Fig. 10 we cut off the remaining points with larger lengthscales on the log-log plots.
By doing scaling form fits of the correlation functions, we can extract the correlation length for each function.We refer to the correlation length derived from geometric criticality of the percolation order parameter ξ p , extracted from the pair connectivity function.The standard correlation length used for thermodynamic criticality of the spin-spin correlation function we refer to as ξ c .In Fig. 10, which is calculated at the thermodynamic critical temperature T 3D c (L → ∞) of the infinite size system for C-3Dx, these two correlation lengths are ξ p = 167.81and ξ c = 19.50.At the thermodynamic critical temperature T 2D c (L → ∞) of the infinite size system for C-2D, these correlation lengths are ξ p = 327.38 and ξ c = 94.65.The fact that the correlation lengths are not infinite at this point is due to the fact that the simulated systems are finite and that the corresponding critical temperatures are nicely shifted from its infinite value because of finite size effects.
Notice that in both cases, ξ p is much larger than ξ c .From the scaling form fits, we have d − 2 + η c = 1.078 ± 0.068 for C-3Dx, and d − 2 + η c = 0.272 ± 0.008.(The results for d − 2 + η p from the fits are already reflected in the insets (a) of Fig. 4 and Fig. 6) They are very close to the established values d − 2 + η c = 1.0336 (3D clean Ising model) [46] and d − 2 + η c = 0.25 (2D clean Ising model) [24,47], corroborating the validity of our results for the correlation lengths.As implied in our previous discussion, these two different correlation lengths correspond to two categories of criticality, geometric criticality and thermodynamic criticality, with two distinct order parameters (the infinite network strength R and the magnetization M respectively).
For C-2D, it is well established that the correlation length exponents ν p (based on geometric clusters) and ν c (based on FK clusters) are different and are related by ν p ≥ ν c [14].Therefore, near the critical point , the connectivity correlation length of geometric clusters is expected to be large compared to the spin-spin correlation length (ξ p > ξ c ).This is consistent with our simulation Fig. 10(b).Our simulation Fig. 10(a) reveals that ξ p > ξ c also applies for C-3Dx.To further validate this inequality, in Fig. 11, we compute the two correlation lengths for C-3Dx geometric criticality as a function of temperature near the critical point T slice p = T 3D c on the disordered side, where the correlation function is consistent with the scaling form shape.As shown by Fig. 11, ξ p is always much larger than ξ c , and they increase as T → T c .This suggests that the exponent inequality ν p ≥ ν c also applies for C-3Dx.

IX. DISCUSSION
So far, we have argued that the bulk thermodynamic critical point of the 3D clean Ising model coincides with the percolation of geometric clusters defined on a 2D slice, T slice p = T 3D c .Also we have numerically investigate this geometric criticality by extracting the critical cluster exponents and studying the behaviors of the pairconnectivity functions (together with the 2D clean Ising model).These results have several important implications: (1) There are two length scales in the clean 3D Ising model.Aside from the usual length scale associated with the spin-spin correlation length ξ c , we have defined a new length scale, ξ p which is the correlation length of the pair connectivity function of geometric clusters defined on a 2D slice.We have furthermore shown that these two length scales diverge differently as T → T 3D c , and that   , we have τ = 2.001 ± .013 and d v,p = 1.856±.018,so that (d+d v,p )/d v,p = 2.078±.010 .Recall that the scaling function associated with the exponent τ suffers from a significant scaling bump, which skews the value of τ to lower values in finite size systems.The slight mismatch in this scaling relation is likely due to this well-known finite size effect for τ [32].
We find some indication that the Coniglio inequalities for these two different sets of critical exponents also hold on a 2D slice of the 3D clean Ising model.For example, our results in Fig. 11 suggest that the exponent inequality derived by Coniglio and coworkers [14] for the correlation length exponents ν p and ν c at the C-2D critical point also holds at the C-3Dx critical point, ν p ≥ ν c .
( (4) All of the above implies that there are two order parameters on the 2D slice of the 3D clean Ising model, just as there are two order parameters in the 2D clean Ising model.One of it is associated with the thermodynamic phase transition, and the other is associated with the geometric cluster percolation.

X. CONCLUSIONS AND OUTLOOK
In summary, we have studied geometric criticality associated with the correlated percolation of interacting geometric clusters on a 2D slice of clean 3D Ising models (C-3Dx).We find that as in the clean 2D Ising model (C-2D), the geometric criticality associated with the percolation of interacting geometric clusters at the C-3Dx critical fixed point corresponds to a unique geometric universality class which is distinct from that of uncorrelated percolation.
In addition, we find that at both the C-3Dx fixed point and the C-2D fixed point, the geometric clusters become fractal, leading to power law behavior in the pair connectivity function.We find that in the vicinity of thermodynamic criticality, the pair connectivity function displays power law behavior over a wider region of the phase diagram than does the spin-spin correlation function.
One consequence of this finding is that the pair connectivity function can be a useful tool for diagnosing criticality in the context of image data from probes which can take "slice" data, and also potentially from scanning surface image probes such as scanning tunneling microscopy, atomic force microscopy, and scanning infrared microscopy.In a future publication, we will explore the relation of these concepts to surface criticality, which is relevant for the application of these ideas to scanning surface probes.Especially in light of the tremendous increase in data coming from a growing number of scanning image probes [2], it is useful to have another method of analysis in hand.We expect that future studies regarding geometric criticality in random bond and random field Ising models will further facilitate the application of geometric cluster analyses to the interpretation of the many types of 2D image probe experimental techniques, especially in cases where complex pattern formation is observed [3][4][5]11].

FIG. 1 .
FIG. 1. Cluster size distribution D(s) ∼ s −τ of internal C-3Dx geometric clusters at T slice p = T 3D c for system sizes from L = 80 to L = 192.The inset (a) shows the extrapolation of τ from the DLD fits in the main panel to the thermodynamic limit L → ∞, which gives τ = 2.001 ± 0.013.The inset (b) shows the scaling collapse of the curves in the main panel, using the extrapolated τ in (a).

FIG. 2 .
FIG.2.The power law fits for Rs ∼ s 1/dv , using internal C-3Dx geometric clusters at T slice p = T 3D c for system sizes from L = 80 to L = 192.The inset (a) shows the extrapolation of dv from the fits in the main panel to the thermodynamic limit L → ∞, which gives dv = 1.856 ± 0.018.The inset (b) shows the scaling collapse of curves in the main panel, using the extrapolated dv in (a).

FIG. 3 .FIG. 4 .
FIG. 3. The power law fits for Rh ∼ h 1/d h , using internal C-3Dx geometric clusters at T slice p = T 3D cfor system sizes from L = 80 to L = 192.The inset (a) shows the extrapolation of d h from the fits in the main panel to the thermodynamic limit L → ∞, which gives d h = 1.714 ± 0.022.The inset (b) shows the scaling collapse of curves in the main panel, using the extrapolated d h in (a) form fits in the main panel.A straightforward fit extrapolating to L → ∞ yields d − 2 + η p = 0.322 ± 0.002.The curves in the main panel collapse onto each other for the extrapolated d − 2 + η p , as shown by the inset (b) of Fig. 4.

FIG. 5 .
FIG. 5. Pair connectivity function gconn(r) on the 2D slice of the clean 3D Ising system with size L = 192 at different temperatures around T slice p = T 3D c .gconn(r) turns up at large r for T < T 3D c .

FIG. 6 .
FIG. 6. Pair connectivity function gconn(r) of the 2D clean Ising system with different sizes from L = 96 to L = 352 at T 2D p = T 2D c .The inset (a) shows the extrapolation of d−2+ηp from the scaling form fits of gconn(r) in the main panel to the thermodynamic limit L → ∞, which yields d − 2 + ηp = 0.104 ± 0.002.The inset (b) shows the scaling collapse of curves in the main panel, using the extrapolated d − 2 + ηp in (a)

FIG. 7 .FIG. 8 .
FIG. 7. Pair connectivity function gconn(r) of the 2D clean Ising system with size L = 256 at different temperatures around T 2D p = T 2D c .The change of behavior for gconn(r) from the scaling form shape to the turning-up shape happens in the close vicinity of T 2D p = T 2D c

FIG. 9 .
FIG.9.Numerical results of the reduced spontaneous magnetization M and the infinite network strength R of the 2D slice embedded in 3D Ising system under zero external field with L = 192 at T < T 3D c , averaged over a number of configurations in the order of magnitude ∼ 10 5 .M and R both approach 1 when T decreases, with M < R.

FIG. 10 .
FIG. 10.Pair connectivity function gconn(r) and spin-spin correlation function gspin(r) for (a) C-3Dx with L = 160 at the critical temperature T 3D c (L → ∞) and (b) C-2D with L = 256 at the critical temperature T 2D c (L → ∞).The black cross markers on the correlation functions denote the lengthscales of correlation lengths.The thick lines represent the power law behaviors of the correlation functions by setting the exponential terms to unity

( 2 )
As discussed in the context of Conjecture #2 in Section IV, for the 2D clean Ising model, the geometric cluster percolation exponents and thermodynamic critical exponents are different in definition and values, but they satisfy the same scaling relation in their own closed sets respectively.As in the purely 2D case, there are also two distinct sets of critical exponents at T 3D c for the clean 3D Ising model.One set comes from the geometric clusters defined on a 2D slice of the clean 3D Ising system, and the other from the bulk FK clusters.Using the numerical values derived in this work, the scaling relation d − 2 + η = 2(d − d v ) is satisfied within error bars.Inserting the values from Table I, d − 2 + η p = 0.322 ± .002, and d v,p = 1.856 ± .018gives 2(d − d v,p ) = 0.288 ± .036 .The scaling relation τ = (d + d v )/d v is almost satisfied: From the Table ) The C-3Dx fixed point represents a new universality class: not only are the critical exponents numerically distinct from the geometric criticality of the 2D clean Ising system and from 2D uncorrelated percolation, but they are also distinct from the values of bulk percolation in the 3D clean Ising model, which happens inside the ordered phase T 3D p < T 3D c = T slice p .

TABLE I .
Theoretical/Numerical values for the 2D geometric cluster self-similarity characterized critical percolation exponents for Ising models and standard percolation model.the generic case is that ξ p > ξ c .One consequence of this is that near this critical point, experiments which can detect clusters on a slice (for example, scanning probes like scanning tunneling microscopy) should reveal that the pair connectivity function defined on a slice is power law over a larger region of the phase diagram than the spin-spin correlation function.That is, the critical region appears to be larger when measured via a pair correlation function than via a spin-spin correlation function.The argument above also applies to 2D clean Ising model.