## From Explosive to Infinite-Order Transitions on a Hyperbolic Network

Vijay Singh, C. T. Brunson and Stefan Boettcher Department of Physics, Emory University, Atlanta, GA, 30322; USA

We analyze the phase transitions that emerge from the recursive design of certain hyperbolic networks that includes, for instance, a discontinuous ("explosive") transition in ordinary percolation. To this end, we solve the q -state Potts model in the analytic continuation for non-integer q with the real-space renormalization group. We find exact expressions for this one-parameter family of models that describe the dramatic transformation of the transition. In particular, this variation in q shows that the discontinuous transition is generic in the regime q &lt; 2 that includes percolation. A continuous ferromagnetic transition is recovered in a singular manner only for the Ising model, q = 2 . For q &gt; 2 the transition immediately transforms into an infinitely smooth order parameter of the Berezinskii-Kosterlitz-Thouless (BKT) type.

PACS numbers: 64.60.ah, 64.60.ae, 64.60.aq

Real-world networks [1-3] exhibit dramatically distinct phenomenology, with a profound imprint of their geometry on the dynamics, when compared with lattices or mean-field systems. What we now call complex networks, aside from being random, possess geometries dominated by small-world bonds and scale-free degree distributions [4, 5]. These lead to novel, and often nonuniversal, scaling behaviors unknown for lattices. For example, they exhibit discontinuous transition even in ordinary percolation [6] as exactly solvable realizations of the so-called 'explosive' percolation transitions [7-13]. Unlike the Achlioptas process, it is easy to prove that this discontinuity exists [6, 14], however, its origin results from the network structure, not from a correlated process [7].

Hyperbolic networks combine a lattice geometry with a hierarchy of small-world connections. They were originally proposed as solvable models for complex networks [15-19]. These recursively defined structures provide deeper insights into small-world effects compared to random networks that otherwise require approximate or numerical methods. Work on percolation [6, 17, 2022], the Ising model [15, 23-25], and the q -state Potts model [19, 26] have shown that critical behavior, once thought to be exotic and model-specific [3], can be universally categorized near the transition point [26, 27] for a large class of hyperbolic networks, such as those discontinuous percolation transitions described in Refs. [6, 14]. They exemplify a non-trivial category of the explosive cluster-growth mechanisms leading to the discontinuous transition [28]. Yet, it has been noted that this discontinuity arises only as a non-generic case of the general theory [27], making its prevalence in percolation on hyperbolic networks somewhat surprising.

Here, we illuminate the origin of the explosive percolation transition by investigating the q -state Potts model on a simple hyperbolic network, MK1 [6]. The analytic continuation in q [29] reveals that the discontinuity that characterizes the explosive percolation transition related to the limit q → 1 , rather than being a special case, per-

FIG. 1. Depiction of hyperbolic networks MK1 of generation k = 5 . The recursive, hierarchical pattern is evident. The network features a regular geometric structures, in form of a one-dimensional backbone, and a distinct set of small-world links (shaded green arched lines). Only the backbone bonds renormalize.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0669v2-figure-cee74ab1-10c0-4b09-b655-3a679ae3d87e.jpg)

sists for all q &lt; 2 . Instead, the Ising ferromagnet ( q = 2 ) emerges as the singular case between all discontinuous ( q &lt; 2 ) and all infinite-order transitions ( q &gt; 2 ), with an infinitely smooth order parameter of the BerezinskiiKosterlitz-Thouless (BKT) type that is consistent with the theory in Ref. [27]. In particular, with real-space renormalization group (RG), we trace the existence of the discontinuous transition to the non-generic behavior of the critical exponent y h that describes the response of the system to a conjugate external field [6, 19]. This exponent, that is temperature-dependent in these hyperbolic networks, happens to reach its maximum exactly at the critical temperature for q &lt; 2 , resulting in merely quadratic corrections in its local expansion. We argue that this exceptional behavior for q &lt; 2 is intimately connected with the manner in which geometric, latticelike features are coupled with the hyperbolic structure.

As a simple and generic member of the class of hyperbolic networks, we consider here MK1, as depicted in Fig. 1. MK1 is recursively generated, starting with two sites connected by a single edge at generation k = 0 . Each new generation combines two sub-networks of the previous generation and adds single edge connecting the

end sites. As a result, the k th generation contains N k = 2 k + 1 vertices, 2 k backbone bonds, and 2 k -1 small-world bonds. It is an effectively one-dimensional version of the small-world Migdal-Kadanoff hierarchical diamond lattice [15, 17], which has been used previously to prove the existence of the discontinuous transition in ordinary percolation [6]. Recently, this discontinuity has been studied for a variety of other hyperbolic networks for bond and site percolation in Ref. [14]. Refs. [19, 26] have used MK1 to study the q -state Potts model for certain integer values of q ≥ 2 . In principle, to relate the Potts model to percolation at q → 1 first requires a q -derivative of the partition function [29]. However, we find that this derivative hardly affects the analytic properties, as we will show below by comparing each quantity obtained for q → 1 with the exactly known percolation results.

In close correspondence to the Ising model on hyperbolic networks discussed in Ref. [24] (see also Ref. [26]), we introduce the following couplings. All variables x i interact with their nearest-neighbors along the backbone with a coupling K 0 (solid links in Fig. 1), while smallworld neighbors interact with a coupling K 1 (shaded links in Fig. 1). Every x i also experiences a uniform external field B . Then, in preparation for applying the renormalization group, the Potts-Hamiltonian can be written as

$$- \beta \mathcal { H } = \sum _ { n = 1 } ^ { N _ { k } / 2 } \left ( - \beta \mathcal { H } _ { n } \right ) + \mathcal { R } \left ( K _ { 1 } \right ), \quad \ \ ( 1 ) \quad \text{For} \quad \text{and}$$

where R contains all remaining coupling terms of higher level in the hierarchy while the x i on the backbone can be sectioned into a sequence of three-variable graphlets. These consist of two adjacent lattice-backbone bonds bridged by an arched small-world bond in Fig. 1, each with their own 'sectional' Hamiltonian,

$$- \beta \mathcal { H } _ { n } & = 2 I + K _ { 0 } \left ( \delta _ { x _ { n - 1 }, x _ { n } } + \delta _ { x _ { n }, x _ { n + 1 } } \right ) + K _ { 1 } \delta _ { x _ { n - 1 }, x _ { n + 1 } } \\ & + B \left [ \left ( \delta _ { 1, x _ { n - 1 } } + \delta _ { 1, x _ { n } } \right ) + \left ( \delta _ { 1, x _ { n } } + \delta _ { 1, x _ { n + 1 } } \right ) \right ]. \ \ ( 2 ) \quad \text{ where}$$

Here, I is a constant that fixes the overall energy scale. We conveniently choose new variables, similar to (inverse) 'activities" [30], C = e -2 I , κ = e -K 0 , and θ = e -B . Similarly defined are the control parameters of temperature µ = e -βJ and field η = e -βh , with β = 1 /kT and h as the external field. The 'raw' (unrenormalized) couplings are considered as uniform: K 0 = K 1 = βJ , where we fix the energy scale via J = 1 , while B = βh . However, we find that K 1 does not change under renormalization and retains its raw value, e -K 1 = µ , a key distinguishing feature of these hyperbolic networks [24, 27]. Thus, the initial values of the renormalizing activities are C 0 = 1 , κ 0 = µ , and θ 0 = η . Then, we can rewrite the sectional Hamiltonian as

$$e ^ { - \beta \mathcal { H } _ { n } } = C ^ { - 1 } \kappa ^ { - \left ( \delta _ { x _ { n - 1 }, x _ { n } } + \delta _ { x _ { n }, x _ { n + 1 } } \right ) } \mu ^ { - \delta _ { x _ { n - 1 }, x _ { n + 1 } } } \\ \theta ^ { - \left [ \left ( \delta _ { 1, x _ { n - 1 } } + \delta _ { 1, x _ { n } } \right ) + \left ( \delta _ { 1, x _ { n } } + \delta _ { 1, x _ { n + 1 } } \right ) \right ] }.$$

The RG consists of successively tracing out the central variable x n and expressing the renormalized activities ( C , κ , θ ′ ′ ′ ) in terms of their priors ( C, κ, θ, µ ) :

$$\text{cmac} \text{using} \quad \sum _ { x _ { n } = 1 } ^ { q } e ^ { - \beta \mathcal { H } _ { n } } = ( C ^ { \prime } ) ^ { - \frac { 1 } { 2 } } \left ( \kappa ^ { \prime } \right ) ^ { - \delta _ { x _ { n - 1 }, x _ { n + 1 } } } ( \theta ^ { \prime } ) ^ { - \left ( \delta _ { 1, x _ { n - 1 } } + \delta _ { 1, x _ { n + 1 } } \right ) }$$

Although the remaining variables x n ± 1 present us with up to q 2 relations between new and old quantities, only three relations are independent. We can obtain the RGflow for ( θ , κ ′ ′ ) at size N k +1 in terms of ( θ, κ ) at size N k ,

$$\underset { \text{tiny} } { \text{top} } \quad \theta ^ { \prime } = \theta \sqrt { \frac { \theta ^ { 2 } + \kappa ^ { 2 } + ( q - 2 ) \, \theta ^ { 2 } \kappa ^ { 2 } } { 1 + ( q - 1 ) \, \theta ^ { 2 } \kappa ^ { 2 } } }, \quad \text{ (3a)}$$

$$\underset { \substack { \kappa \nmid \\ \text{per-} } } { \kappa ^ { \prime } } = \frac { \cdot } { \kappa \mu \left [ 1 + \theta ^ { 2 } + ( q - 2 ) \theta ^ { 2 } \kappa \right ] } \, \left ( 3 \right )$$

Note the explicit dependence on the control parameter µ (temperature), characteristic for these networks [27].

First, we analyze the RG-flow in Eq. (3) in the absence of an external field. For η = 1 = θ 0 , it remains θ k = 1 for all k &gt; 0 and the RG-flow for the bond-coupling becomes

$$\kappa ^ { \prime } = \kappa \mu \, \frac { 2 + ( q - 2 ) \kappa } { 1 + ( q - 1 ) \kappa ^ { 2 } }. \quad \quad ( 4 )$$

For the Ising model, q = 2 , Eq. (4) reduces with κ ˆ = jk and µ ˆ = j to the equivalent of Eq. (6) in Ref. [26], and for percolation, q = 1 , it reduces with κ ˆ1 = -T and µ ˆ1 = -p to Eq. (1) in Ref. [6].

In the thermodynamic limit ( N k →∞ ) , Eq. (3) provides the fixed points κ ′ = κ = κ ∞ :

$$\kappa _ { \infty } ^ { 0 } = 0, \quad \kappa _ { \infty } ^ { \pm } = \frac { \mu } { 2 } \frac { ( q - 2 ) } { ( q - 1 ) } \left [ 1 \pm \frac { \sqrt { \mathcal { D } _ { q } \left ( \mu \right ) } } { \mu \left ( q - 2 \right ) } \right ], \ \ ( 5 )$$

where κ + ∞ &gt; κ -∞ , with

$$\mathcal { D } _ { q } \left ( \mu \right ) = \mu ^ { 2 } \left ( q - 2 \right ) ^ { 2 } + 4 \left ( 2 \mu - 1 \right ) \left ( q - 1 \right ). \quad \left ( 6 \right )$$

The strong-coupling limit κ 0 ∞ = 0 always exists for any q . However, the Ising case and the limit of Eq. (5) that corresponds to percolation, both have to be considered with some care. Yet, we easily reproduce Eq. (2) in Ref. [6]. We plot the behavior of κ ± ∞ in Fig. 2 as a function of µ for a range of q that highlights the peculiarities of both cases.

To obtain the thermal exponent y t , we calculate the relevant eigenvalues along each branch of the fixed-point lines. The Jacobian of the RG-flow in Eq. (4), evaluated at κ = κ ∞ , provides the eigenvalue

$$\lambda _ { t } = \frac { \partial \kappa ^ { \prime } } { \partial \kappa } \Big | _ { \kappa _ { \infty } } = \frac { 2 \mu \left [ 1 + \left ( q - 2 \right ) \kappa _ { \infty } - \left ( q - 1 \right ) \kappa _ { \infty } ^ { 2 } \right ] } { \left [ 1 + \left ( q - 1 \right ) \kappa _ { \infty } ^ { 2 } \right ] ^ { 2 } }. \ \ ( 7 )$$

FIG. 2. Plot of the fixed-point couplings κ ± ∞ in Eq. (5) as a function of temperature µ for values of q = 1 2 , 1 2 7 , , (black lines, from left to right on the bottom). The blue-shaded line locates the branch-point singularity in κ ± ∞ , evaluated at D q ( µ ) = 0 in Eq. (6), which starts from µ = 0 and κ = -∞ at q = 1 and rises to cross κ = 0 exactly at µ = 1 2 for q = 2 , above which it recedes back to µ = 0 for q →∞ . The greenshaded line marks the raw coupling κ 0 = µ from which the RG-flow in Eq. (4) initiates and flows vertically towards the nearest stable fixed-point, κ + ∞ &gt; 0 or κ 0 ∞ = 0 . For all q &gt; 2 , the location of the branch point provides the critical point µ c ( &lt; 1 2 ) with an infinite-order BKT transition [27]. For all q &lt; 2 , µ c = 1 2 and the transition becomes discontinuous.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0669v2-figure-14b7c488-3486-4d5a-810f-0d4b4a075003.jpg)

Hence, with Eqs. (5), we get

$$\lambda _ { t } ^ { 0 } & = 2 \mu, \\ \lambda _ { t } ^ { \pm } & = \frac { 1 - \mu } { q ^ { 2 } \mu } \left \{ 4 \left ( q - 1 \right ) + \mu \left ( q - 2 \right ) ^ { 2 } \left [ 1 \pm \frac { \sqrt { \mathcal { D } _ { q } \left ( \mu \right ) } } { \mu \left ( q - 2 \right ) } \right ] \right \}. \quad \ y _ { h } \sim$$

The behavior of the eigenvalues, depicted in Fig. 3, follows closely the description developed for general hyperbolic lattices in Ref. [27].

As Fig. 2 shows, for q &lt; 2 the critical point is located always at µ c = 1 2 , such that the relevant thermal scaling exponent [31] y t = log 2 λ + t is given by

$$y _ { t } \sim \frac { 2 \left ( \mu _ { c } - \mu \right ) } { \ln 2 } \quad \left ( q < 2, \mu \rightarrow \mu _ { c } = \frac { 1 } { 2 } \right ), \quad ( 9 )$$

which indicates that y t vanishes for all q &lt; 2 . As the generic theory suggests [27], this linear behavior is to be expected when κ + ∞ is regular at µ = µ c such that it possesses a linear slope, which is satisfied for q &lt; 2 , see Fig. 2. Eq. (7) implies that this will remain true also for percolation, q = 1 , in particular. However, for the Ising ferromagnet, and exactly only at q = 2 , λ t in Eq. (7), and hence y t , becomes a function of κ 2 ∞ . Miraculously, at q = 2 the branch-point in κ ∞ is placed exactly on the κ = 0 -axis at µ c with κ ± ∞ ∼ ± √ 2 ( µ -µ c ) . Both effects combine into a linear correction but a different coefficient,

$$y _ { t } \sim \frac { 4 \left ( \mu _ { c } - \mu \right ) } { \ln 2 } \quad \left ( q = 2, \mu \rightarrow \mu _ { c } = \frac { 1 } { 2 } \right ). \quad ( 1 0 )$$

For q &gt; 2 , κ + ∞ no longer intercepts the strong-coupling fixed-point line at κ 0 ∞ = 0 , see Fig. 2, rendering the issue moot. The critical behavior with an infinite-order BKTlike transition is now described by analyzing the RG-flow near the branch point itself, as previously described [15, 19, 24, 27], moving the critical point to µ c = µ B &lt; 1 2 for q &gt; 2 , where D q ( µ B ) = 0 in Eq. (6) gives µ B = 2 [ q √ q -1 -2 ( q -1) ] / ( q -2) 2 .

We now study the effect of an external field for η → 1 to analyze the behavior of the order parameter near the transition. In the limit η → 1 , i.e., θ → 1 , and κ → κ ∞ , Eqs. (3) provide the Jacobian matrix ∂ ( θ ,κ ′ ′ ) ∂ θ,κ ( ) , which turns out to be upper-triangular, since ∂θ ′ ∂κ = 0 at θ = 1 . Hence, the thermal exponent y t remains unaffected by the magnetic field. Accordingly, we obtain the magnetic exponent y h = log 2 λ h with

$$\Lambda _ { h } & = \frac { \partial \theta ^ { \prime } } { \partial \theta } = \frac { ( q - 2 ) \kappa _ { \infty } ^ { 2 } + 2 } { ( q - 1 ) \kappa _ { \infty } ^ { 2 } + 1 } \\ & = \frac { 2 } { \mu q } + \frac { q - 2 } { 2 \mu q \left ( q - 1 \right ) } \left [ \mu ( 3 q - 2 ) - \sqrt { \mathcal { D } _ { q } \left ( \mu \right ) } \right ],$$

when evaluated along κ + ∞ near µ c = 1 2 . As λ h is purely a function of κ 2 ∞ , there are no odd correction terms for the exponent near µ c ,

$$\Big \} ^ { ( 8 ) } \quad y _ { h } \sim 1 - \frac { 8 q } { \left ( 2 - q \right ) ^ { 2 } } \frac { \left ( \mu _ { c } - \mu \right ) ^ { 2 } } { \ln 2 } \quad \left ( q < 2, \mu \rightarrow \mu _ { c } \right ). \ \left ( 1 2 \right )$$

At q = 2 , we arrive at a more generic form, but merely by the fact that κ + ∞ now has its root-singularity:

$$y _ { h } \sim 1 - 2 \frac { ( \mu _ { c } - \mu ) } { \ln 2 } \quad ( q = 2, \mu \to \mu _ { c } ) \,. \quad ( 1 3 )$$

Using the standard scaling relation [31] for the behavior of the order parameter near the transition,

$$m \sim ( \mu _ { c } - \mu ) ^ { \beta } \,, \quad \beta = \frac { 1 - y _ { h } } { y _ { t } } \,, \quad \quad ( 1 4 )$$

we find β = 0 for all q &lt; 2 , while β = 1 2 for q = 2 . The later result, predicting a 2nd -order transition with a mean-field-like exponent, has been obtained previously [26]. However, the former result suggests that the discontinuous percolation transitions found on various hyperbolic networks [6, 14, 32] is not an artifact of the particular percolation limit q → 1 of the q -state Potts model but is maintained for all q &lt; 2 . As soon as q &gt; 2 (now with µ c = µ B ), 1 -y h remains finite while y t → 0 for µ → µ c such that β diverges and the transition becomes instantly infinite-order, as explained in Ref. [27].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0669v2-figure-e08eee84-0dbb-47f4-a308-1bf5105d1140.jpg)

FIG. 3. Plot of the Jacobian eigenvalues λ in Eq. (8) as a function of temperature µ for the q -state Potts model on MK1 for q = 1 1 2 2 7 , . , , . For all q at the strong-coupling fixed point κ 0 ∞ = 0 we have λ 0 = 2 µ , while the two branches λ ± belonging to non-zero κ ± ∞ change dramatically with q and vary non-trivially with µ . At q = 1 , λ -= 0 and λ + = 2 -2 µ remain disconnected but coalesce for 1 &lt; q &lt; 2 (here, q = 1 2 . ) with a branch point at unity (dashed line). At q = 2 , both branches degenerate into λ ± = (1 -µ /µ ) . For q &gt; 2 (here, q = 7 ), the branches re-open and resemble the case 1 &lt; q &lt; 2 . Stable fixed points correspond to λ &lt; 1 (black), unstable ones to λ &gt; 1 (red).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0669v2-figure-0d4ae41b-b9f6-4ca9-9ba4-ce57b6499201.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0669v2-figure-9b4c1a54-7415-4873-a2ad-17c3480ca715.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0669v2-figure-bcfe542b-5e6f-442d-8f13-9b9217ecbb91.jpg)

FIG. 4. Plot of the exponent Ψ as a function of temperature µ for q = 1 , 3 2 , 2 , and 7. The physical branch is given by the thickened, dark lines. The red-shaded or dashed lines mark the unphysical branches of Ψ . Note that Ψ( ) µ is continuous at µ c = 1 2 for q &lt; 2 with a parabolic continuation but degenerates at q = 2 (i.e., for the Ising model) to a linear continuation. For q &gt; 2 , Ψ( ) µ is discontinuous at µ c = µ B &lt; 1 2 and continues after the drop-off on the lower branch of the root-singularity.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0669v2-figure-cd8323ec-4229-4c83-9292-2933f42143d3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0669v2-figure-0fd51139-d861-4edd-9e3f-740eb166a382.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0669v2-figure-e1d6e865-6286-446d-bead-26d5ccf4410c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0669v2-figure-257124b4-1477-4dc4-bac2-bfde8fef82d0.jpg)

Another perspective on these peculiar transitions is provided by the exponent [26]

$$\Psi = 2 y _ { h } - 1 \quad \quad \ \ ( 1 5 ) \quad \text{re}$$

for µ &gt; µ c , and Ψ ≡ 1 for µ &lt; µ c , that describes the divergence of the susceptibility with the system size,

$$\chi \sim N ^ { \Psi ( \mu ) }, \quad \ \ \ ( 1 6 ) \quad \ \ \ t i,$$

glyph[negationslash]

not only near the transition but for all temperatures 0 ≤ µ ≤ 1 . For percolation, χ is related to the average size of the largest cluster, 〈 S p ( ) 〉 , after a derivative in q . In Fig. 4, we plot Ψ( ) µ in Eq. (15) for various q . The non-generic behavior of y h also manifests itself in Ψ : no branch of Ψ ever exceeds unity, as could be generically expected, but rather Ψ maintains a global extremum exactly of unit height for all q = 2 . Also Ref. [32] has found the same 'delicate' behavior in Ψ while studying the transformation of the transition from discontinuous to infinite-order for a one-parameter family of percolation models. The generic case can be observed, for instance, for the families of Ising models studied in Ref. [24] by applying an external field [33], where y h and hence Ψ exhibit linear corrections near µ c because the branch point singularity is outside the physical regime and Ψ develops a maximum at an unphysical value Ψ &gt; 1 , rendering the magnetization order parameter continuous at µ c .

While we found that the RG results for the couplings and y t of the q -state Potts are continuously connected to those of percolation at q → 1 , this is not true for Ψ (i.e., y h or possibly other quantities that are affected by a q -derivate of the partition function). However, a direct comparison reveals that the Potts model at q = 1 in Eq. (12) and the percolation result from Eq. (14) in Ref. [6] merely differ by a factor of 2 in their 2 nd -order correction near µ c . Their difference is more pronounced for larger µ , of course. For example, Ψ remains positive for all µ &lt; 1 in percolation and only vanishes for zero bond-probability, p = 1 -µ = 0 , suggesting that diverging clusters are always possible. In turn, Ψ vanishes for non-trivial values of µ in the q -state Potts model, see Fig. 4, as is easily obtained from Eq. (11). While insignificant for q &lt; 2 , this remains true, remarkably, at q ≥ 2 , which implies another potentially interesting transition in the behavior of χ at temperatures above criticality where fluctuations become independent of system size N , as shown in Fig. 5. There, we demonstrate the collapse of the data for χ when properly rescaled with system size N according to Eq. (16). We observe a local minimum at µ c which reaches a finite value on this scale for q &lt; 2 but appears to rapidly go to zero for all larger q with increasing N . For q → 1 , a further factor of q -1 is required to converge to a thermodynamic limit. (Remember that χ at q = 1 is not directly related to the percolation cluster size, which would require an additional derivative of q in the partition function.) Most remarkable is the spike in χ at the point µ where Ψ( µ ) becomes zero, see Fig. 4. Only above such temperature do fluctuations behave accord-

FIG. 5. Plot of the susceptibility χ as a function of temperature µ , rescaled with the appropriate power of system size N according to Eq. (16) using Ψ( µ ) (and Ψ ≡ 1 for µ &lt; µ c ). From left to right, panels refer to q ↘ 1 , q = 1 2 . , q = 2 , and q = 7 . For the larger values of q , the rescaled χ are show for a sequence of increasing system sizes to demonstrate convergence in the thermodynamic limit, while for the case q → 1 we show that χ requires a further rescaling by 1 ( / q -1) to converge to a finite value throughout (at some fixed, already large system size, here N = 2 512 ). Note that in this scaling, the susceptibility develops a minimum at the corresponding value of µ c , however, it is sharply peaked at the respective value at which Ψ reaches (and then remains at) zero, see Fig. 4.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0669v2-figure-026950c1-d892-4b5f-a4b6-f068ab65615d.jpg)

ing to a high-temperature regime, independent of system size N . Such a transition from finite-size to diverging but subextensive fluctuations was absent, for instance, in the Ising model studied on any of the networks considered in Ref. [24].

To understand what makes these networks behave in such a non-generic way, let us consider the case of percolation on a binary tree, the most extreme hyperbolic network. Starting from the root of the tree at level i = 0 , the average number of sites at branching level i &gt; 0 that are connected in a single cluster with the root is 〈 n i +1 〉 = 2 p 〈 n i 〉 , hence, the average size of the rooted cluster at system size N k = 2 k +1 -1 is 〈 S k ( p ) 〉 = ∑ k i =0 〈 n i 〉 = (2 ) [ p k +1 -1 ] / (2 p -1) , or 〈 S p ( ) 〉 ∼ N Ψ( ) p with Ψ( ) = log p 2 (2 p ) for 1 2 &lt; p &lt; p c = 1 . Just like on any other hyperbolic network, there are diverging, yet sub-extensive, clusters well below p c . (In fact, such diverging cluster emerge as soon as there are spanning clusters connecting the root to the perimeter of the tree, for

1 2 &lt; p &lt; p c .) However, here we find linear corrections,

$$\Psi \left ( p \right ) \sim 1 - \frac { 1 } { \ln 2 } \left ( p _ { c } - p \right ) \quad \left ( p \rightarrow p _ { c } \right ), \quad \left ( 1 7 \right )$$

remembering that µ ˆ1 = -p . We suggest that the relevant feature that distinguishes the networks with non-trivial explosive percolation transitions [28] with p c &lt; 1 from the trivial cases with p c = 1 , like this tree, is the interplay between a geometric backbone and the small-world hyperbolic links. For instance, adding lateral links turn the tree into a non-amenable graph with a non-trivial transition, studied in Refs. [18, 34]. Similarly, MK1 here, or the hyperbolic Hanoi networks in Ref. [14], become nontrivial by mixing geometric with hyperbolic structures. It is remarkable, that in Ref. [22] the tree-approximation for Ψ( ) p for percolation on a hyperbolic network, that is exact for p &lt; p l but predicts a trivial transition at p c = 1 with linear corrections in Ψ( ) p , breaks down above p l exactly where the geometric backbone begins to facilitate spanning clusters [17] and ultimately leads

to a non-trivial explosive transition at p c &lt; 1 , again with quadratic corrections in Ψ( ) p [14]. While renormalization group calculations provide these exact results, it will likely require more detailed numerical simulations or better models to elucidate the role of geometry in the dynamics of cluster formation that lead to such corrections.

Finally, we speculate on the origin of the singularity of the q -state Potts model for q → 2 . What distinguishes the q ≥ 2 case from that for q &lt; 2 is the nature of the Z q -symmetry [29]. While for q &lt; 2 (and, correspondingly, for percolation) it is sufficient to merge two clusters by merely adding a link, for q ≥ 2 (i.e., for the Ising spin model or other multi-state degrees of freedom) two cluster only merge when both share the same orientation. Apparently, the analytic continuation in q delimits these two distinct cluster-growth processes via a singularity at q → 2 .

We like to thank T. Nogawa and T. Hasegawa for fruitful discussions. This work was supported by NSF (DMR #1207431 &amp; IOS #1208126) and James S. McDonnell Foundation (#220020321).

- [1] M. Barthelemy, EPL (Europhysics Letters) 63 , 915 (2003).
- [2] S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, and D.-U. Hwang, Phys. Rep. 424 , 175 (2006).
- [3] S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes, Rev. Mod. Phys. 80 , 1275 (2008).
- [4] D. J. Watts and S. H. Strogatz, Nature 393 , 440 (1998).
- [5] A.-L. Barabasi and R. Albert, Science 286 , 509 (1999).
- [6] S. Boettcher, V. Singh, and R. M. Ziff, Nature Communications 3 , 787 (2012).
- [7] D. Achlioptas, R. M. D'Souza, and J. Spencer, Science 323 , 1453 (2009).
- [8] R. A. da Costa, S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes, Phys. Rev. Lett. 105 , 255701 (2010).
- [9] P. Grassberger, C. Christensen, G. Bizhani, S.-W. Son, and M. Paczuski, Phys. Rev. Lett. 106 , 225701 (2011).
- [10] H. K. Lee, B. J. Kim, and H. Park, Phys. Rev. E 84 ,

020101 (2011).

- [11] O. Riordan and L. Warnke, Science 333 , 322 (2011).
- [12] N. Bastas, P. Giazitzidis, M. Maragakis, and K. Kosmidis, Physica A: Statistical Mechanics and its Applications 407 , 54 (2014).
- [13] N. A. M. Araújo, P. Grassberger, B. Kahng, K. J. Schrenk, and R. M. Ziff, (arXiv:1404.5325).
- [14] V. Singh and S. Boettcher, Physical Review E 90 , 012117 (2014).
- [15] M. Hinczewski and A. N. Berker, Phys. Rev. E 73 , 066126 (2006).
- [16] S. Boettcher, B. Gonçalves, and H. Guclu, J. Phys. A: Math. Theor. 41 , 252001 (2008).
- [17] S. Boettcher, J. L. Cook, and R. M. Ziff, Phys. Rev. E 80 , 041115 (2009).
- [18] T. Nogawa and T. Hasegawa, J. Phys. A: Math. Theor. 42 , 145001 (2009).
- [19] T. Nogawa, T. Hasegawa, and K. Nemoto, Phys. Rev. Lett. 108 , 255703 (2012).
- [20] A. N. Berker, M. Hinczewski, and R. R. Netz, Phys. Rev. E 80 , 041118 (2009).
- [21] P. Minnhagen and S. K. Baek, Phys. Rev. E 82 , 011113 (2010).
- [22] T. Hasegawa and T. Nogawa, Phys. Rev. E 87 , 032810 (2013).
- [23] M. Bauer, S. Coulomb, and S. N. Dorogovtsev, Phys. Rev. Lett. 94 , 200602 (2005).
- [24] S. Boettcher and C. T. Brunson, Phys. Rev. E 83 , 021103 (2011).
- [25] S. K. Baek, H. Mäkelä, P. Minnhagen, and B. J. Kim, Phys. Rev. E 84 , 032103 (2011).
- [26] T. Nogawa, T. Hasegawa, and K. Nemoto, Phys. Rev. E 86 , 030102 (2012).
- [27] S. Boettcher and C. T. Brunson, (arXiv:1209.3447).
- [28] Y. S. Cho and B. Kahng, (arXiv:1404.4470).
- [29] F. Y. Wu, Reviews of Modern Physics 54 , 235 (1982).
- [30] M. Plischke and B. Bergersen, Equilibrium Statistical Physics, 2nd edition (World Scientifc, Singapore, 1994).
- [31] R. K. Pathria, Statistical Mechanics, 2nd Ed. (Butterworth-Heinemann, 1996).
- [32] T. Nogawa and T. Hasegawa, Phys. Rev. E 89 , 042803 (2014).
- [33] C. T. Brunson and S. Boettcher, (in preparation).
- [34] T. Hasegawa, T. Nogawa, and K. Nemoto, EuroPhys. Lett. 104 , 16006 (2013).
Vijay Singh C. T. Brunson Stefan Boettcher
From explosive to infinite-order transitions on a hyperbolic network
We analyze the phase transitions that emerge from the recursive design of
certain hyperbolic networks that includes, for instance, a discontinuous
("explosive") transition in ordinary percolation. To this end, we solve the
$q$-state Potts model in the analytic continuation for non-integer $q$ with the
real-space renormalization group. We find exact expressions for this
one-parameter family of models that describe the dramatic transformation of the
transition. In particular, this variation in $q$ shows that the discontinuous
transition is generic in the regime $q<2$ that includes percolation. A
continuous ferromagnetic transition is recovered in a singular manner only for
the Ising model, $q=2$. For $q>2$ the transition immediately transforms into an
infinitely smooth order parameter of the Berezinskii-Kosterlitz-Thouless (BKT)
type.