BUDAPEST UNIVERSITYOF TECHNOLOGYAND ECONOMICS
DOCTORAL THESIS
Fractals in dimension theory and complex networks
Author: Supervisor: István KOLOSSVÁRY Dr. Károly SIMON
Doctoral School of Mathematics and Computer Science Faculty of Natural Sciences
2019
iii Acknowledgements
It is a great pleasure to thank the many people who have had a direct impact on my academic life in the past years. Foremost, I thank my supervisor Károly Simon. Beyond his deep knowledge of the field, his enthusiasm for mathematics and the unique way he gives it on has truly had a great influence on how I think about math and the world around us;
Secondly, I thank my co-authors Balázs Bárány, Gergely Kiss, Júlia Komjáthy and Lajos Vágó. I gained valuable experience from each and every joint project;
I thank the referees of my home defense for their thorough work. In particular, the comments of István Fazekas, which greatly improved the clarity of one of the chapters.
Moreover, I thank the support and generous hospitality of all the colleagues at the Department of Stochastics. I feel very fortunate to have landed here after my bachelor years, thanks to Doma Szász. The atmosphere is inspiring for research and at the same time very friendly;
and last but not least all the support I get from all my family and friends, who are in some way part of my life, even though many don’t know much about what I actually do on a daily basis.
I acknowledge the financial support of different grants and scholarships without which all my travels, meeting many new colleagues, presenting results at conferences and learning many interesting topics would not have been possible.
The template used for the thesis can be found at https://www.latextemplates. com/template/masters-doctoral-thesis. I am very satisfied, it saved me from lots of unnecessary headaches.
v
Contents
Acknowledgements iii
List of Figures vii
List of Symbols ix
List of Acronyms xi
1 Introduction1 1.1 Fractal geometry...... 1 1.2 Self-affine sets...... 2 1.2.1 Planar carpets...... 4 1.2.2 Fractal curves...... 5 1.3 Apollonian networks...... 6 1.4 Informal explanation of contribution...... 8 1.4.1 Dimension of triangular Gatzouras–Lalley-type planar carpets8 1.4.2 Pointwise regularity of zipper fractal curves...... 13 1.4.3 Distances in Random Apollonian Networks...... 15
2 Triangular Gatzouras–Lalley-type planar carpets with overlaps 19 2.1 Triangular Gatzouras–Lalley-type carpets...... 19 2.1.1 Results of Gatzouras and Lalley...... 22 2.1.2 Separation conditions...... 23 Separation of the cylinder parallelograms...... 23 Separation of the columns...... 25 2.2 Main results...... 26 2.2.1 Hausdorff dimension...... 26 2.2.2 Box dimension...... 27 2.2.3 Diagonally homogeneous carpets...... 29 2.3 Preliminaries...... 31 2.3.1 Symbolic notation...... 31 2.3.2 Atypical parallelograms...... 33 2.3.3 Ledrappier–Young formula...... 34 2.4 Upper bound for dimH Λ ...... 35 2.5 Proof of Theorem 2.2.2...... 38 2.5.1 The proof of Theorem 2.2.2 assuming Claim 2.5.1 and Proposi- tion 2.5.2...... 40 2.5.2 The proof of Claim 2.5.1...... 40 2.5.3 Proof of Proposition 2.5.2...... 41 2.6 Proof of results for box dimension...... 45 2.6.1 Diagonally homogeneous subsystems...... 46 2.6.2 Counting intersections...... 49 2.6.3 Proof of Theorem 2.2.7...... 51 vi
2.6.4 Proof of Theorem 2.2.8...... 52 2.7 Examples...... 53 2.7.1 The self-affine smiley: a non diagonally homogeneous example 53 2.7.2 Example for dimH Λ = dimB Λ ...... 53 2.7.3 Overlapping example...... 54 2.7.4 Example "X X"...... 55 ≡ 2.7.5 Negative entries in the main diagonal...... 55 2.7.6 A family of self-affine continuous curves...... 56 2.8 Three-dimensional applications...... 58
3 Pointwise regularity of parameterized affine zipper fractal curves 63 3.1 Self-affine zippers satisfying dominated splitting...... 63 3.2 Main results...... 66 3.3 Pressure for matrices with dominated splitting of index-1...... 68 3.4 Pointwise Hölder exponent for non-degenerate curves...... 73 3.5 Zippers with Assumption A...... 80 3.6 An example, de Rham’s curve...... 85
4 Distances in random and evolving Apollonian networks 89 4.1 Definitions and notations...... 89 4.2 Main results...... 91 4.2.1 Related literature...... 92 4.3 Structure of RANs and EANs...... 93 4.3.1 Tree-like structure of RANs and EANs...... 93 4.3.2 Distances in RANs and EANs: the main idea...... 97 4.3.3 Combinatorial analysis of shortcut edges...... 98 4.4 Distances in RANs and EANs...... 100 4.4.1 A continuous time branching process...... 100 4.4.2 Proof of Theorem 4.2.1 and 4.2.5...... 103 4.4.3 Proof of Theorem 4.2.3...... 106
A Basic dimension theoretic definitions 113
B No Dimension Drop is equivalent to Weak Almost Unique Coding 115
Bibliography 119 vii
List of Figures
1.1 A Bedford–McMullen carpet...... 4 1.2 A de Rham curve...... 5 1.3 Linearly parametrized de Rham curve with parameter ω = 1/10....6 1.4 An Apollonian gasket and network...... 7 1.5 Random Apollonian Networks...... 8 1.6 A GL and TGL carpet...... 9 1.7 The "self-affine smiley"...... 9 1.8 Example of Falconer and Miao together with overlapping version... 10 1.9 TGL carpets with different overlaps...... 10 1.10 Example "X X"...... 11 ≡ 1.11 A family of self-affine continuous fractal curves...... 11 1.12 Three-dimensional application...... 12 1.13 Affine zippers in the plane...... 13 1.14 A RAN after a few steps...... 15 1.15 Impact of shortcut edges on diameter of RAN...... 17
2 2.1 The skewness of Ri1...in := fi1...in ([0, 1] ) ...... 21 2.2 The IFS , where z and (1, z) are identified...... 24 T 2.3 Intersecting parallelograms Rı and R in the proof of Lemma 2.5.4.... 44 2.4 Intersecting parallelograms Rı and R in the proof of Lemma 2.6.8.... 51 2.5 Orientation reversing maps generally destroy the column structure.. 56
3.1 An affine zipper...... 64 3.2 Local neighbourhood of points in Bn,l,m ...... 74 3.3 Well ordered property...... 82
4.1 Coding RANs, initial stpes...... 94 4.2 Coding RANs, induction step...... 94 4.3 Tree like structure of RANs...... 96 4.4 Shortest path from u to v ...... 97
ix
List of Symbols
α (A) ... α (A) singular values of a d d matrix A 1 ≥ ≥ d × α(x), αr(x) pointwise and regular Hölder exponent (3.1.3) , , , Iterated Function Systems F G H T Γ, Λ, Ω attractor of an Iterated Function System Σ, Σ symbolic spaces H Π, Π , Π , π natural projections from a symbolic space s Hs H δ, (δ-approximate) s-dimensional Hausdorff measure (A.0.1) H H2 µd, σd e expectation and variance of Yd (4.1.3) p, q, λ probability vectors p∗ optimal vector for Hausdorff dimension p optimal vector for box dimension µp Bernoulli measure on a symbolic space νep push forward of µp φs singular value function (A.0.4)
D(p) formula for dimH µp (2.2.1) E(β), Er(β) β-level set of α(x) and αr(x) (3.1.5) Id(x) large deviation rate function of Yd (4.1.4) P(t) matrix pressure function (3.2.1) M(x) cone centered at x (3.1.6) Yd full coupon collector block with d + 1 coupons (4.1.2) Diam, Flood, Hop diameter, flooding time, hopcount in a graph (4.1.1) dim an unspecified dimension dimA affinity dimension (A.0.5) dimB, dimB lower and upper box dimension (A.0.3) dimH Hausdorff dimension (A.0.2) dimP packing dimension projx orthogonal projection to x-axis
xi
List of Acronyms
i.i.d. independent and identically distributed ...... 93 u.a.r. uniformly at random ...... 89 w.h.p. with high probability...... 91
AN Apollonian network ...... 7
AUC Almost Uniqe Coding ...... 25
CLT Central Limit Theorem ...... 91
CTBP continuous time branching process ...... 100
EAN evolving Apollonian network ...... 89
GL Gatzouras–Lalley ...... 20
HESC Hochman’s Exponential Separation Condition...... 25
IFS Iterated Function System...... 2
NDD No Dimension Drop ...... 26
OSC Open Set Condition ...... 3
RAN random Apollonian network ...... 89
ROSC Rectangular Open Set Condition ...... 3
SOSC Strong Open Set Condition ...... 3
TGL triangular Gatzouras–Lalley-type ...... 19
WAUC Weak Almost Uniqe Coding...... 25
1
Chapter 1
Introduction
1.1 Fractal geometry
The word "fractal" comes from the Latin fractus¯ meaning "broken" or "fractured". Benoit B. Mandelbrot coined this term when he wrote in his book Fractals: Form, Chance and Dimension in 1977 that
Many important spatial patterns of Nature are either irregular or frag- mented to such an extreme degree that ... classical geometry ... is hardly of any help in describing their form. ... I hope to show that it is possible in many cases to remedy this absence of geometric representation by using a family of shapes I propose to call fractals – or fractal sets. []
Roughly speaking, any object detailed on arbitrarily small scales that resembles itself in some way on different magnification scales can be called a fractal. Since the 1970s, 80s, fractal geometry has become an important area of mathematics with many connections to theory and practice alike. Fractals have found application in geomet- ric measure theory, dimension theory, dynamical systems, number theory, analysis, differential equations or probability theory to name a few. Also there is increasing interest in more applied areas of mathematics and natural sciences such as network theory, wavelets, percolation problems all the way to computer graphics, image com- pression, financial markets, fluid turbulence or fractal antenna, etc. The list is bound to expand further in the coming years.
Aim and structure of thesis The main aim of the Thesis is to demonstrate this diverse applicability of fractals in different areas of mathematics. Namely,
1. widen the class of planar self-affine carpets for which we can calculate the dif- ferent dimensions especially in the presence of overlapping cylinders,
2. perform multifractal analysis for the pointwise Hölder exponent of a family of continuous parameterized fractal curves in Rd including deRham’s curve,
3. show how hierarchical structure can be used to determine the asymptotic growth of the distance between two vertices and the diameter of a random graph model, which can be derived from the Apollonian circle packing problem.
Since these topics are not confined to a narrow area of mathematics, a great deal of effort has been put into making the presentation accessible to a wider mathematical audience. 2 Chapter 1. Introduction
This Chapter contains a brief (far from exhaustive) introduction to the topics • and is concluded with a section containing informal explanations of the contri- butions made to the different topics, together with the various methods used in the proofs.
For those not familiar with the terminologies used, AppendixA contains a brief • background material with the basic definitions and results.
Besides the theoretical results, several examples with illustrative pictures are • provided for explanation.
Plenty of figures assist the Reader through proofs. • Chapters2,3 and4 contain the precise definitions and rigorous formulations of our re- sults, together with the proofs. They are based on the papers [KS18; BKK18; KKV16], respectively.
1.2 Self-affine sets
Let (X, d) be a compact metric space. Usually, we will work on a compact subset of Euclidean space Rd. A very natural technique to construct fractals is via Iterated Function Systems (IFSs). An IFS consists of a finite collection of contracting maps F f : X X for i = 1, . . . , N. Huthinson proved in his seminal paper [Hut81], that for i → every IFS there exists a unique non-empty compact set Λ, called the attractor, which satisfies Λ = fi(Λ), i [N] ∈[ where [N] = 1, 2, . . . , N . An extensively studied class of IFSs are the self-affine sets, { } in which case each map of the IFS is an affine transformation, i.e.
fi(x) = Aix + ti, where A Rd d is a contracting, invertible matrix and t Rd is a translation vector. i ∈ × i ∈ An important further subclass consists of the self-similar sets, when the matrices can be written in the form Ai = riOi, where Oi is an orthogonal matrix and 0 < ri < 1 is the contracting ratio in every direction. A natural way to depict an IFS is to provide the images fi(R), where R is the small- est rectangle which contains Λ. Without loss of generality we may assume throughout this thesis that R = [0, 1]2. The correspondence between the IFS and a figure showing the collection of images of R will be unique in our study, since the maps do not contain any rotations or reflections (except in Subsection 2.7.5). See Figure 1.1 for example. Perhaps the most fundamental question in fractal geometry is to determine the dimension of a set. Roughly speaking, dimension indicates how much space a set occupies near to each of its points. Several different types of fractal dimension are used. For basic dimension theoretic definitions such as the Hausdorff, packing and (lower and upper) box dimension of a set and the Hausdorff and local dimension of measures, see AppendixA. Throughout, the Hausdorff, packing, lower and upper box dimension will be denoted by dimH, dimP, dimB, dimB and dimB, respectively. The relative position of cylinder sets fi(R) greatly influence the degree of difficulty to calculate the dimension. The simplest case is the Strong Separation Property (SSP) when all sets fi(Λ) are pairwise disjoint. Somewhat weaker is the (Strong) Open 1.2. Self-affine sets 3
Set Condition. The most difficult is when there is heavy overlapping between the cylinders. This will be central in Chapter2. Definition 1.2.1. An IFS with attractor Λ satisfies the Strong Open Set Condition (SOSC), F if there exists a non-empty open set U, with Λ U = ∅ and such that ∩ 6 f (U) U with f (U) f (U) = ∅ for i = j. (1.2.1) i ⊆ i ∩ j 6 i [N] ∈[ In particular, if U can be chosen to be R, then we say that the Rectangular Open Set Condition (ROSC) holds. This will usually be the case. If U above and Λ can be disjoint, then the Open Set Condition (OSC) holds. In case of self-similar sets satisfying the OSC, all mentioned dimensions are equal to s, often called the similarity dimension, which is the solution of the Hutchinson– Moran formula s ∑ ri = 1, i [N] ∈ see [Fal90, Section 9.2]. Regardless of overlaps, in the self-similar case, the similarity dimension is always an upper bound for the dimensions considered in this thesis. The analog upper bound for self-affine sets is the affinity dimension dimAff, introduced by Falconer [Fal88a], which comes from the "most natural" cover of the set, see (A.0.5) for the definition. All self-affine sets satisfy
dim Λ dim Λ dim Λ min dim Λ, d . H ≤ P ≤ B ≤ { Aff } A central question for the past 30 years has been to determine what can cause the drop of dimension (from the natural upper bounds). An obvious cause is the presence of an exact overlap, i.e. there are two distinct sequences i ,..., i = j ,... j such that 1 n 6 1 k f ... f (Λ) = f ... f (Λ). Another cause in higher dimensions can be the i1 ◦ ◦ in j1 ◦ ◦ jk highly regular alignment of cylinder sets, planar carpets are great examples. The full picture is not completely understood even in the simplest self-similar case on the real line. The transversality method of Pollicott and Simon [PS95], further developed in [Sol95; PS96; PS00] etc., has proven to be a useful tool in determining the dimension of a parametrized family if IFSs for almost all parameter values. Recently, Hochman [Hoc14; Hoc15] made a big breakthrough for the Hausdorff dimension of self-similar measures, which in particular implies that if an IFS on the real line is defined by algebraic parameters, then the drop of dimension is equivalent to having exact overlaps. In a generic sense, equality of dimensions is typical for self-affine sets. Falconer proved in his seminal paper [Fal88a] that for fixed linear parts A ,..., A if A < { 1 N} k ik 1/3 and the translations are chosen randomly according to N d dimensional Lebesgue × measure then all the aforementioned dimensions of the self-affine set are equal. The 1/3 bound was later relaxed by Solomyak [Sol98a] to 1/2, which is sharp due to an example of Przytycki and Urba´nski[PU89]. Building on the mentioned result of Hochman and results about the Ledrappier–Young formula for self-affine measures [BK17] (see Subsection 2.3.3), very recently Bárány, Hochman and Rapaport [BHR17] greatly improved these results in two dimensions by giving specific, but mild condi- tions on A ,..., A under which the dimensions are equal. { 1 N} However, in specific cases, which do not fall under these conditions, strict in- equality is possible. Planar carpets form a large class of examples in R2 for which this exceptional behavior is typical. The highly regular column and/or row structure causes the drop of the Hausdorff dimension. 4 Chapter 1. Introduction
1.2.1 Planar carpets Independently of each other, Bedford [Bed84] and McMullen [McM84] were the first to study planar carpets. They split the unit square R into m columns of equal width and n rows of equal height for some integers n > m 2 and considered iterated ≥ function systems of the form
1/m 0 x i/m f (x) := + (i,j) 0 1/n y j/n for (i, j) A 0, . . . , m 1 0, . . . , n 1 , see Figure 1.1. They gave explicit ∈ ⊆ { − } × { − } formula for the Hausdorff and box-counting dimension of the corresponding attractor Λ. It turns out that dimH Λ = dimB Λ is atypical, namely equality holds if and only if all non-empty columns have the same number of elements.
FIGURE 1.1: A Bedford–McMullen carpet: IFS on left, attractor on right. In this case dimH Λ < dimB Λ.
Later Gatzouras and Lalley [GL92] generalized the results to IFSs, which kept the (non-overlapping) column structure, however the width of the columns could vary and no row structure is required. They still assumed ROSC (1.2.1) and that the width of a rectangle is larger than its height, we say that direction-x dominates. The precise definition is given in Definition 2.1.3, see also Figure 1.6 for an example. More recently, Bara´nsky[Bar07] kept the row and column structure, but relaxed the direction-x dominates assumption by allowing an arbitrary subdivision of the hor- izontal and vertical axis. After appropriately choosing which direction is "dominant", the results resemble that of [GL92]. Diagonal systems assuming only ROSC (1.2.1) and no further restrictions on the translations were studied by Feng–Wang [FW05] q and Fraser [Fra12]. Former determined the L spectrum of self-affine measures νp and in particular the box dimension of the attractor. In [Fra12] linear isometries which map [ 1, 1]2 to itself are allowed and the box dimension is determined. Fraser called − these box-like sets. Observe that in all the mentioned papers the ROSC (1.2.1) was assumed. Carpets with overlaps were not studied until the last few years. Fraser and Shmerkin [FS16] shift the columns of Bedford–McMullen carpets to get overlaps, while Pardo-Simón [PS] allows shifts in both directions of Bara´nskicarpets. Relying on a recent breakthrough by Hochman [Hoc14] on the dimension of self-similar measures on the line, both pa- pers show that apart from a small exceptional set of parameters the results in [Bed84; McM84] and [Bar07] remain valid in the overlapping case. This is the type of shifted columns that can be seen in Figure 1.9. 1.2. Self-affine sets 5
In Chapter2 triangular Gatzouras–Lalley carpets are introduced: the column struc- ture of Gatzouras–Lalley carpets are kept, however instead of diagonal matrices, lower triangular ones are used to define the linear parts of the maps. Hence, the usual rectangular cylinder sets become parallelograms with two vertical sides, see Figure 1.6. Moreover, we also allow different types of overlaps in the construction, see Figure 1.9. The results of [GL92] are generalized to this setting. The parallelograms and overlaps give a much more flexible framework in which existing and many new examples can be treated in a unified way.
1.2.2 Fractal curves Fractal curves have already appeared in the 19th century, although at the time the un- usual constructions were considered by many to be mathematical "monsters". Peano created the first space-filling curve, Weierstrass presented the first example for a con- tinuous but nowhere differentiable function and von Koch gave a more geometric definition of a fractal curve simply referred to today as the von Koch snowflake. Since then, fractal curves have found abundant applications in areas such as wavelets, in- terpolation functions or signal processing. The starting point of our work was the curve introduced and studied by de Rham [Rha47; Rha56; Rha59]. The geometric construction of the curve goes as follows. Starting from the square [0, 1] [ 1, 0], it can be obtained by trisecting each side with ratios ω : (1 2ω) : ω × − − and "cutting the corners" by connecting each adjacent partitioning point to get an octagon. Again, each side is divided into three parts with the same ratio and adjacent partitioning points are connected, and so on. The de Rham curve is the limit curve of this procedure. With a more analytic approach, in the language of iterated function systems we can say that de Rham’s curve is the attractor Γ of the IFS
ω 0 0 1 2ω ω 2ω f0(x) = x and f1(x) = − x + , (1.2.2) ω 1 2ω − 2ω 0 ω 0 − where ω (0, 1/2) is the parameter. More precisely, the self-affine curve Γ defined ∈ by (1.2.2) gives the segment between two midpoints of the original square. Observe that the fixed point of f is z = [0, 1]T, while for f it is z = [1, 0]T. Furthermore, 0 0 − 1 2 f0(z2) = z1 = f1(z0), we say that the cross-condition holds. This ensures that Γ is continuous. Figure 1.2 shows both approaches: first four steps of the "corner cutting" on left and the images of f0 and f1 after one (red), two (green) and three (black) levels of iteration.
FIGURE 1.2: A de Rham curve: geometric construction (left), IFS con- struction (right). 6 Chapter 1. Introduction
The main goal of Chapter3 is to analyze the pointwise regularity of such self- affine curves given by a linear parametrization in the much more general framework of zipper IFSs as defined in [ATK03], see Definition 3.1.1. In case of de Rham’s curve, a linear parametrization v : [0, 1] R2 is of the form 7→ i i + 1 v(x) = f (v(2x i)), for x , , i = 0, 1. (1.2.3) i − ∈ 2 2 h Figure 1.3 shows a visualisation of a linearly parametrized de Rham curve. Such linear parameterizations occur in the study of wavelet functions in a natural way, see for example Protasov [Pro06], Protasov and Guglielmi [PG15], and Seuret [Seu16].
FIGURE 1.3: Linearly parametrized de Rham curve with parameter ω = 1/10. Left: The image of the unit cube w.r.t the IFS generating the graph of the de Rham curve. Middle: The second iteration. Right: The curve itself.
The analysis is done by performing multifractal analysis, i.e. to study the possible values, which occur as (regular) pointwise Hölder exponents, and determine the mag- nitude of the sets (in terms of Hausdorff dimension), where it appears. This property was studied for several types of singular functions, for example for wavelets by Barral and Seuret [BS05], Seuret [Seu16], for Weierstrass-type functions by Otani [Ota17], for complex analogues of the Takagi function by Jaerisch and Sumi [JS17] or for different functional equations by Coiffard, Melot and Willer [CMT14], by Okamura [Oka16] and by Slimane [BS03] etc. Our results applied to de Rham’s curve give finer results than existing ones in the literature.
1.3 Apollonian networks
Since the seminal work of Erd˝osand Rényi [ER60], the theory of random graphs has received immense interest in many areas of science. Social networks, the World Wide Web, traffic/shipping routes etc. have made it essential to create models which cap- ture the main features and driving forces of these real-world networks, often also called complex networks. These main features include large clustering coefficient, power law decay of degree distribution and small distances between nodes (logarithmic with size) meaning that networks tend to contain smaller, dense communities; there are few really large "hubs" and many-many much smaller ones; finally "it’s a small world after all". A particularly successful class of models are the Preferential attachment models [BA99; BR04; Bol+01; FP16]. In these dynamically evolving graphs, the probability that a new node attaches to an existing one is proportional to the degree of the existing one. 1.3. Apollonian networks 7
Thus, nodes with already large degree tend to get new edges more frequently, hence the term preferential attachment. This simple mechanism gives rise to the power law decay. Several variants of the model have been defined, here we focus on the shortest paths between vertices of Apollonian networks. The construction of deterministic and random Apollonian networks originates from the problem of Apollonian circle packing: starting with three mutually tangent circles, we inscribe in the interstice formed by the three initial circles the unique circle that is tangent to all of them: this fourth circle is known as the inner Soddy-circle. It- eratively, for each new interstice its inner Soddy-circle is drawn, see Figure 1.4. After infinite steps the result is a fractal set, an Apollonian gasket [Boy82; Gra+03]. An Apollonian network (AN) is the resulting graph if we place a vertex in the center of each circle and connect two vertices if and only if the corresponding circles are tangent, see Figure 1.4. This model was introduced independently by Andrade et al. [And+05] and Doye and Massen [DM05] as a model for real-world networks. These networks have the important properties mentioned above, in addition, due to the construction, Apollonian networks also show hierarchical structure: another property very commonly observed.
FIGURE 1.4: An Apollonian gasket (left), the deterministic network generated by it (right).
It is straightforward to generalize Apollonian packings to Rd for d 2 with ≥ mutually tangent d dimensional hyperspheres. Analogously, if each d-hypersphere corresponds to a vertex and vertices are connected by an edge if the correspond- ing d-hyperspheres are tangent, then we obtain a d-dimensional AN (see [Zha+08; Zha+06]). The network arising by this construction is deterministic. Zhou et al. [ZYW05] proposed to randomize the dynamics of the model such that in one step only one interstice is picked uniformly at random and filled with a new circle. This construc- tion yields a d dimensional random Apollonian network (RAN) [ZRC06], see Figure 1.5. Using heuristic and rigorous arguments the results in [AM08; CFU14; DS07; Ebr+13; FT12; ZRC06; ZYW05] show that RANs have the above mentioned main features of real-world networks. A different random version of the original Apollonian network was introduced by Zhang et al. [ZRZ06], called Evolutionary Apollonian networks (EAN) where in every step every interstice is picked and filled independently of each other with probability q. If an interstice is not filled in a given step, it can be filled in the next step again. We call q the occupation parameter. For q = 1 we get back the deterministic AN model. It is 8 Chapter 1. Introduction
FIGURE 1.5: An array of Random Apollonian Networks in two dimen- sions with n = 50, 250, 750 vertices. conjectured in [ZRZ06] that an EAN with parameter q, as q 0, should show similar → topological behaviour to RANs. To make this statement rigorous, instead of looking at a sequence of evolving EAN-s with decreasing parameters, we slightly modify the model and investigate the asymptotic behaviour of a single EAN when q might differ in each step of the dynamics. That is, we consider a series q ∞ of occupation { n}n=1 parameters so that qn applies for step n of the dynamics, and assume that qn tends to 0. In this setting, the interesting question is to determine the correct rate for qn that achieves the observation that EAN shows similar behaviour as RAN when the parameter tends to zero. The main goal of Chapter4 is to determine the precise asymptotic growth of the shortest path between two vertices and the diameter of the graph.
1.4 Informal explanation of contribution
The main body of the Thesis is based on three articles [KS18; BKK18; KKV16]. We now informally explain the contribution and methods used in each of these papers.
1.4.1 Dimension of triangular Gatzouras–Lalley-type planar carpets Gatzouras–Lalley (GL) carpets [GL92] are the attractors of self-affine IFSs on the plane whose first level cylinders are aligned into columns using orientation preserving maps with linear parts given by diagonal matrices. In Chapter2, we consider a nat- ural generalization of such carpets by replacing the diagonal matrices with lower triangular ones so that the column structure is preserved, see Definition 2.1.1. We call them Triangular Gatzouras–Lalley-type (TGL) planar carpets, indicating that the linear part of the maps defining the IFS are triangular matrices and it is a natural generalization of the Gatzouras–Lalley construction. The shaded rectangles and parallelograms in Figure 1.6 show the images of R un- der the maps defining a GL carpet on the left and a TGL carpet on the right. These are typical examples which satisfy the ROSC (1.2.1). Furthermore, there is a corre- spondence between the rectangles and parallelograms so that the height and width of corresponding ones coincide. We call the Gatzouras–Lalley carpet the GL-brother of the TGL carpet, see Definition 2.1.4. Even though the ROSC holds, it is not immedi- ate that the dimension of the two attractors should be the same. The parallelograms can be placed in a way that there is no bi-Lipschitz map between the two attrac- tors. Nevertheless, Bara´nskiessentially shows in [Bar08] that assuming the ROSC the 1.4. Informal explanation of contribution 9
Hausdorff and box dimension of a TGL carpet is equal to the respective dimension of its GL brother.
FIGURE 1.6: The IFS defining a Gatzouras–Lalley carpet on left and triangular Gatzouras–Lalley-type carpet on right, which are brothers.
The IFS in Figure 1.6 is an example for which dimH Λ < dimB Λ < dimAff Λ. If the orthogonal projection of Λ to the x-axis is the whole [0, 1] interval, then the box- and affinity dimensions are equal. Figure 1.7 shows such an example, where the outlines of fi(R) are shown together with the attractor, which we call the "self-affine smiley". A special class of examples consists of affine IFSs in which all matrices have the same main diagonal. We call them diagonally homogeneous carpets. Such a particu- lar TGL carpet (see the left hand side of Figure 1.8) was introduced by Falconer and Miao [FM07, Figure 1 (a)], where the box dimension of the attractor was calculated. Later, Bárány [B15´ , Subsection 4.3] showed that for this example the box and the Hausdorff dimensions are the same. This is due to the fact that all columns have the same number of maps. The well-known Bedford–McMullen carpets [Bed84; McM84] form a proper subclass of these TGL carpets. In all these examples only the boundary of the cylinder sets fi(R) could intersect. However, the main contribution of Chapter2 is to continue the works [Bar08; GL92] to allow different types of overlaps in the construction. Figure 1.9 illustrates the overlaps we consider. All three examples are brothers of the GL carpet in Fig- ure 1.6. On the left, the columns are shifted in a way that the IFS on the x-axis gen- erated by the columns, denoted by , satisfies Hochman’s Exponential Separation H Condition (HESC), see Definition 2.1.9. This type of shifted columns were considered
FIGURE 1.7: The "self-affine smiley", whose dimH = 1.20665 . . . < 1.21340 . . . = dimB = dimA, see Subsection 2.7.1. 10 Chapter 1. Introduction
FIGURE 1.8: The attractor Λa from Subsection 2.7.2 with parameter a = 3/10 (left) and Subsection 2.7.3 with parameter a = 3/20 (right), shown together with the outlines of the images of fi(R).
by Fraser and Shmerkin [FS16] and Pardo-Simón [PS] on different carpets. In the cen- ter, columns do not overlap, however, parallelograms within a column may do so if a certain transversality like condition holds. The one on the right on Figure 1.9 has both types of overlaps.
FIGURE 1.9: Triangular Gatzouras–Lalley-type carpets with different overlaps. Left: shifted columns satisfying Hochman’s exponential sep- aration condition. Center: non-overlapping columns, transversality condition. Right: mixture of both.
By modifying the translation vectors in the example on the left hand side of Fig- ure 1.8, we got a brother with overlaps seen on the right hand side, for which we show in Subsection 2.7.3 that transversality holds. Another concrete overlapping ex- ample satisfying transversality is "X X" in Figure 1.10, for which there is strict ≡ inequality between the Hausdorff, box and affinity dimensions. If instead the con- struction would be "X = X", then the Hausdorff and box dimensions would be equal. Moreover, if there were no empty columns in this example, then the box and affinity dimensions would coincide.
Results Section 2.2 contains the formal statements of all the main results. Roughly speaking, we show that for any TGL carpet Λ
dim Λ dim Λ, ≤ e 1.4. Informal explanation of contribution 11
FIGURE 1.10: Example "X X" from Subsection 2.7.4 for which ≡ dimH Λ = 1.13259 . . . < dimB Λ = 1.13626 . . . < dimA Λ = 1.2170 . . . .
where Λ is the GL brother of Λ and dim means either box or Hausdorff dimension, see Theorems 2.2.1 and 2.2.4. When ROSC holds and the IFS generated by the H columnse satisfies Hochman’s condition, then equality can be deduced from recent works [BK17; Fra12]. Also, equivalent conditions are given for the equality of the different dimensions. Our main contribution is that in the presence of overlaps described above, we give sufficient conditions under which dim Λ does not drop below dim Λ, see Theo- rems 2.2.2 and 2.2.7. In particular, for the Hausdorff dimension we allow both types of overlaps simultaneously (like the third figure in Figure 1.9), howevere for the box dimension it is either one or the other type (like the first two in Figure 1.9). For a discussion on generalizing towards orientation reversing maps, see Subsec- tion 2.7.5. In particular, we calculate the dimension of a family of self-affine continu- ous curves Λ , which is generated by an IFS containing a map that reflects on the a Fa y-axis, see Figure 1.11. They are examples for the type of fractal curves we study from a different perspective in Chapter3. The formal treatment of all the mentioned examples is done in Section 2.7.
FIGURE 1.11: Left: first (red), second (green) and third (black) level cylinders of a. Right: Λa rotated 90 degrees from Subsection 2.7.6 with parameterF a = 0.2 (black), 0.12 (red) and 0.08 (green).
One motivation to study self-affine fractals of overlapping construction is that 12 Chapter 1. Introduction
sometimes the dimension of a higher dimensional fractal of non-overlapping con- struction coincides with its lower dimensional orthogonal projection which can be a self-affine fractal of overlapping construction. We obtain such a set in 3D by starting from a TGL carpet with overlaps on the xy-plane and then "lift" it to 3D so that the in- teriors of the first level cylinders are disjoint. Figure 1.12 shows such an example with the first level cylinders (left), the attractor (center) and the projection of the cylinders and attractor to the xy-plane (right). Section 2.8 contains the formal treatment of this type of construction.
FIGURE 1.12: A three-dimensional fractal whose dimension is equal to the dimension of its orthogonal projection to the xy-plane.
Methods used to handle overlaps For each type of overlap and dimension we used different methods:
The upper bounds on the Hausdorff and box dimensions (after some simple ob- • servations) follow from proper adaptations of the results of Gatzouras–Lalley [GL92] and Fraser [Fra12], respectively.
To estimate the Hausdorff dimension from below we use the Ledrappier-Young • formula of Bárány and Käenmäki [BK17] (cited in Theorem 2.3.4) for self-affine measures. We show that this lower bound equals the upper bound
– in case of overlapping like on the second figure of Figure 1.9 by an argu- ment inspired by the transversality method introduced in [BRS]; – in case of overlapping like on the third figure of Figure 1.9, we introduce a new separation condition for the self-similar IFS obtained as the projection of the TGL carpet under consideration to the horizontal line. This sepa- ration condition is a non-trivial consequence of Hochman’s Exponential Separation Condition [Hoc14]. We prove this in AppendixB since it could be of separate interest.
To estimate the box dimension from below we could not simply use the Haus- • dorff dimension of the attractor, because, in our case, it is (typically) strictly smaller than the box dimension. Therefore,
– in case of overlapping like on the first figure of Figure 1.9, we used the method of Fraser and Shmerkin [FS16]: the main idea is to pass to a spe- cial subsystem of a higher iterate of the IFS which has non-overlapping columns; 1.4. Informal explanation of contribution 13
– in case of overlapping like on the second figure of Figure 1.9, we intro- duced a new argument to count overlapping boxes. It uses transversality and a result of Lalley [Lal88] based on renewal theory, which gives the pre- cise asymptotics of the number of boxes needed to cover the projection of the attractor to the horizontal line.
1.4.2 Pointwise regularity of zipper fractal curves We postpone the general definition of zipper fractal curves, as given in [ATK03], un- til Chapter3. For now it is enough to consider a self-affine IFS = f0,..., fN 1 F { − } in which all matrices have strictly positive entries and map [0, 1]2 into itself. Fur- T T thermore, the fixed points of f0 and fN 1 are [0, 0] and [1, 1] , respectively, and T T − fi([1, 1] ) = fi+1([0, 0] ). Figure 1.13 shows such examples in the plane with the first (red), second (green) and third (black) level cylinders visible.
FIGURE 1.13: Affine zippers in the plane defined by matrices with strictly positive entries.
Let λ = (λ0,..., λN 1) be a probability vector, which defines a subdivision of − [0, 1]. Let gi be the affine function mapping the unit interval [0, 1] to the ith subinterval of the division (indexing from 0 to N 1). Then a linear parametrization of the attractor − Γ of is a function v : [0, 1] Rd defined by the functional equation F → 1 v(x) = f v(g− (x)) if x g ([0, 1]). i i ∈ i In comparison, we call F : Rd R a self-similar function if there exists a bounded 7→ open set U Rd, and contracting similarities g ,..., g of Rd such that g (U) ⊂ 1 k i ∩ g (U) = ∅ and g (U) U for every i = j, and a smooth function g : Rd R, j i ⊂ 6 7→ and real numbers ρ < 1 for i = 1, . . . , k such that | i| k 1 F(x) = ∑ ρi F(gi− (x)) + g(x), (1.4.1) i=1 see [Jaf97b, Definition 2.1]. The graph of F over the attractor of g ,..., g can be { 1 k} written as the unique, non-empty, compact invariant set of the family of functions d+1 S1,..., Sk in R , where
Si(x, y) = (gi(x), ρiy + g(gi(x))). 14 Chapter 1. Introduction
The multifractal formalism of the pointwise Hölder exponent of self-similar func- tions was studied in several aspects, see for example Aouidi and Slimane [AS04], Slimane [BS01; BS03; BS12] and Saka [Sak05]. The main difference between the self-similar function F defined in (1.4.1) and v defined in (3.1.2) is the contraction part. Namely, while F is a real valued function rescaled by only a real number, the function v is Rd valued and a strict affine trans- formation is acting on it. This makes the study of such functions more difficult. When studying the pointwise Hölder exponent
log v(x) v(y) α(x) = lim inf k − k, y x log x y → | − | we need a good control over the distance v(x) v(y) as y x. This is much k − k → simpler for the self-similar function F, since F was scaled only by a constant. Roughly speaking
F(g (x)) F(g (y)) ρ ρ F(x) F(y) . k i1,...,in − i1,...,in k ≈ | i1 ··· in |k − k Whereas, in the case of self-affine systems, this is not true anymore. That is,
v(g (x)) v(g (y)) A A (v(x) v(y)) . k i1,...,in − i1,...,in k ≈ k i1 ··· in − k However, in general A A (v(x) v(y)) A A v(x) v(y) . In order k i1 ··· in − k 6≈ k i1 ··· in kk − k to be able to compare the distance v(g (x)) v(g (y)) with the norm of the k i1,...,in − i1,...,in k product of matrices, we need an extra assumption on the family of matrices. This notion is called dominated splitting, see Definition 3.1.2. The example of matrices with strictly positive entries satisfies this condition.
Results The key technical tool is the matrix pressure function P(t) defined in (3.2.1). Under a mild non-degeneracy condition (3.1.7) and dominated splitting we show that α(x) is equal to a constant α = P (0) for Lebesgue almost all x [0, 1]. Furthermore, the 0 ∈ multifractal formalism holds for the Hausdorff dimension of the level sets E(β) = x [0, 1] : α(x) = β of a linear parametrization of an affine zipper in a small { ∈ b } interval of possible β values. In other words dimH E(β) is equal to the Legendre- transform of P(t), see Theorem 3.2.1. Moreover, an important contribution was to find an additional rather natural as- sumption, see Assumption A (3.2), which is equivalent for the lim inf in the definition of α(x) to exist as a limit for Lebesgue almost every point, see Theorem 3.2.3. We call the limit the regular Hölder exponent and denote it αr(x). Under these conditions Theorem 3.2.2 states that the multifractal formalism for αr(x) can be extended to the full spectrum of possible β values. For example, de Rham’s curve and matrices with strictly positive entries satisfy Assumption A, hence these results can be applied to get stronger results than the ones in the literature. In general, our results do not imply differentiability of the linear parametrization v, because the lower bound we obtain for α can be strictly smaller than 1. However, in the particular case of de Rham’s curve we show in Proposition 3.6.2 that v is differ- entiable for Lebesgue-almost every x withb derivative vector equal to zero and the set where it is not differentiable has strictly positive Hausdorff dimension. 1.4. Informal explanation of contribution 15
Unfortunately, even in the simplest case of positive matrices, there is not much hope for the computation of the precise values of P(t). Only fast approximation algo- rithms exist, see works of Pollicott and Vytnova [PV15] and Morris [Mor18].
Methods The most important technical contribution was to generalize the results of Feng [Fen03], and Feng and Lau [FL02] for the pressure function of infinite products of positive ma- trices to the more general dominated splitting setting. This is based on results by Bochi and Gourmelon [BG09] characterizing systems with dominated splitting and general ergodic theoretic machinery of Bowen [Bow08]. The proof of the main results use the properties of the pressure function and rely on handling the points which are far away symbolically but close on the self-affine curve, which is done by a fine study of the underlying dynamical system associated with the matrices Ai.
1.4.3 Distances in Random Apollonian Networks We postpone the precise definition of RANs and EANs in arbitrary dimensions until Chapter4, instead just focus on RANs in two dimensions. There is a natural representation of RANs as evolving triangulations in two di- mensions: take a planar embedding of the complete graph on four vertices as in Fig- ure 1.14 and in each step pick a face of the graph uniformly at random, insert a vertex and connect it with the vertices of the chosen triangle (face). The result is a maxi- mal planar graph. Hence, a (RAN2(n))n N is equivalent to an increasing family of ∈ triangulations by successive addition of faces in the plane, called stack-triangulations. Stack-triangulations were investigated in [AM08] where the authors also considered typical properties under different weighted measures, (e.g. ones picked u.a.r. having n faces). Under a certain measure stack-triangulations with n faces are an equivalent formulation of RAN2(n), see [AM08] and references therein.
FIGURE 1.14: A RAN2(n) after n = 0, 2, 8 steps
Let u and v be two vertices of a RAN with n vertices. To study distances we denote by Hop(n, u, v) the hopcount between the vertices, i.e., the number of edges on (one of) the shortest paths between u and v. The flooding time Flood(n, u) is the maximal hopcount from u, while the diameter Diam(n) is the maximal flooding time (4.1.1). The terminology comes from the intuitive picture that if we were to open a source of fluid at u that spreads with unit speed along the edges, then the flooding time is exactly the time when the fluid reaches all other vertices of the graph. 16 Chapter 1. Introduction
Results We determine the exact asymptotic growth of Hop, Flood and Diam as the number of vertices tend to infinity. The results are valid in all dimensions d and all constants are explicit, only depending on d. All three quantities scale as c log n for different constants, which confirms the small-world property. More precisely, we show that for both RANs and EANs the hopcount obeys a cen- tral limit theorem (CLT), when the two vertices are chosen randomly according to a well-defined probability distribution. The centralizing constant is of order log n while the normalizing one is of order log n, see Theorems 4.2.1 and 4.2.5. Furthermore, for RANs the flooding time and diameter divided by log n tend in probability to well p defined constants, see Theorem 4.2.3. The constants that arise have very intuitive meaning once the structure of these graphs is understood. Understanding the hierarchical structure using the fractal viewpoint is perhaps the most important contribution of this work. It gave the opportunity to analyze all the quantities simultaneously, independently of the dimension in a unified way. The CLT for the hopcount is novel, it is especially interesting, because to the best of our knowledge in all previous cases, when a CLT was proved, the underlying graph model had random edgeweights, while RANs and EANs do not. Determining the diameter is always a greater challenge than the hopcount. The flooding time is very rarely studied, in particular for RANs it was not known before. The exact relation of our results compared to existing ones in the literature is surveyed in Subsection 4.2.1. We remark that the article [KKV16] where these results appeared also contain re- sults about the degree distribution and clustering coefficient, however they are not presented in this thesis.
Methods The proofs rely on a few key observations. The first is to realize that RANs have a nice tree-like structure: each vertex can be assigned a unique finite code and all edges can be grouped into two categories, either tree edges or shortcut edges. The tree edges go "down" in the hierarchy of the graph, while a shortcut edge goes back "upward" toward the initial graph along a branch of the tree. See Figure 4.3 for an illustration of the graph in Figure 1.14. All neighbors of a vertex higher in the hierarchy can be determined simply by looking at its unique code. This allows a combinatorial analysis and the well-known coupon collector problem comes up naturally. The idea of using codes can be attributed to my co-author Lajos Vágó, which was then further developed. Secondly, there is a natural embedding of the evolution of a RAN into the evolu- tion of a continuous time branching process (CTBP). This allows for the application of results for CTBPs [AN04; BD06; Büh71] to determine the average and maximal depth of the tree created by the tree edges. The CLT for the hopcount can be derived from here. However, for the flooding time and diameter an extra important observation is needed. Due to the shortcut edges, the deepest branch of the tree created by the tree edges is not necessarily the farthest from the initial graph in graph distance. There can be many branches slightly shorter, however farther from the initial graph, because the shortcut edges climb to the top slower. Figure 1.15 illustrates this phenomena. To determine the exact constants, we achieve the maximal distance by an entropy vs energy argument. 1.4. Informal explanation of contribution 17
branch farthest in graph distance
branch with most edges
other branches
shortcut edge
FIGURE 1.15: Impact of shortcut edges on diameter of RAN: branch with most edges is not the one farthest in graph distance.
The precise proofs use techniques from renewal theory, CTBP theory, large devia- tions, second moment method, symbolic dynamics and combinatorics.
19
Chapter 2
Triangular Gatzouras–Lalley-type planar carpets with overlaps
This chapter is based on [KS18] written jointly with Károly Simon. Recall Figure 1.6.
2.1 Triangular Gatzouras–Lalley-type carpets
Denote the closed unit square by R = [0, 1] [0, 1]. Let = A ,..., A be a family × A { 1 N} of 2 2 invertible, strictly contractive, real-valued lower triangular matrices. The × corresponding self-affine IFS is the collection of affine maps
N bi 0 ti,1 = fi(x) := Aix + ti i=1, where Ai = and ti = , (2.1.1) F { } di ai ti ,2 for translation vectors t , with t , t 0. We assume that a , b (0, 1). i i,1 i,2 ≥ i i ∈ Orthogonal projection of to the horizontal x-axis, denoted proj , generates an F x important self-similar IFS on the line
= h (x) := b x + t N . (2.1.2) H { i i i,1}i=1 We denote the attractor of eand eby Λ = Λ and Λ respectively. F H F H Definition 2.1.1. We say that an IFSe of the form (2.1.1) isetriangular Gatzouras–Lalley- type (TGL) and we call its attractor Λ a TGL planar carpet if the following conditions hold: (a) direction-x dominates, i.e.
0 < a < b < 1 for all i [N] := 1, 2, . . . , N , (2.1.3) i i ∈ { } (b) column structure: there exists a partition of [N] into M > 1 sets ,..., with cardi- I1 IM nality = N > 0 so that |Iıˆ| ıˆ
1 = 1, . . . , N1 and ıˆ = N1 + ... + Nıˆ 1 + 1, . . . , N1 + ... + Nıˆ (2.1.4) I { } I { − } for ıˆ = 2, . . . , M. Assume that for two distinct indices k and ` 1, . . . , N ∈ { }
bk = b` =: rıˆ, if there exists ıˆ 1, . . . , M such that k, ` ıˆ, then (2.1.5) ∈ { } ∈ I (tk,1 = t`,1 =: uıˆ.
We also introduce = h (x) := r x + u M , (2.1.6) H { ıˆ ıˆ ıˆ}ıˆ=1 and we observe that the attractor Λ of is identical with Λ . H H H e 20 Chapter 2. Triangular Gatzouras–Lalley-type planar carpets with overlaps
(c) we assume that ∑j aj 1 holds for every ıˆ 1, . . . , M and the non-overlapping ∈Iıˆ ≤ ∈ { } column structure
u + r u for ıˆ = 1, . . . , M 1 and u + r 1. (2.1.7) ıˆ ıˆ ≤ ıˆ+1 − M M ≤ (d) Without loss of generality we always assume in this paper that
(A1) f (R) R for all i [N] and i ⊂ ∈ (A2) The smallest and the largest fixed points of the functions of are 0 and 1 respec- H tively.
Observe that the definition allows overlaps within columns (like the second figure in Fig- ure 1.9), but columns do not overlap. We say that Λ is a shifted TGL carpet if we drop the assumption (2.1.7), that is non- M overlapping column structure is NOT assumed, we require only that ∑ rıˆ 1 (like the ıˆ=1 ≤ first figure in Figure 1.9).
We often consider the following special cases:
Definition 2.1.2. We say that a shifted TGL carpet Λ has uniform vertical fibres if
s s ∑ aj− H = 1 for every ıˆ [M], (2.1.8) j ∈ ∈Iıˆ where s = dimB Λ and s = dimB Λ . H H Furthermore, we call Λ a diagonally homogeneous shifted TGL carpet if
b b and a a for every i [N]. i ≡ i ≡ ∈ In particular, a diagonally homogeneous carpet has uniform vertical fibres if N/M N and ∈ N = N/M for every ıˆ 1, . . . , M . ıˆ ∈ { } The special case when Nıˆ = 1 for all ıˆ = 1, . . . , M is treated in the paper of Bárány, Rams and Simon [BRS, Lemma 3.1]. GL carpets are just special cases of TGL carpets.
Definition 2.1.3. A self-affine IFS is a Gatzouras–Lalley (GL) IFS and its attractor Λ is F a GL carpet if is a TGL IFS as in Definition 2.1.1 with the additional assumptions that F all off-diagonal elements di = 0 ande the rectangular open set condition (ROSC) holds, recalle Definition 1.2.1.e
Definition 2.1.4. Let Λ be a shifted TGL carpet generated from the IFS of the form (2.1.1). F We say that the Gatzouras–Lalley IFS
˜ ˜ ˜ ˜ ˜ N ˜ bi 0 ˜ ti,1 = fi(x) := Aix + ti i=1, where Ai = and ti = , F { } 0 a˜i t˜i ,2 e and its attractor Λ is the GL brother of and Λ, respectively, if a˜ = a and b˜ = b for F i i i i every i [N], furthermore, has the same column structure (2.1.10) as . If the shifted ∈ F F TGL carpet Λ is actuallye a TGL (that is Λ has non-overlapping column structure) then we also require that ti,1 = ti,1 holdse for all i [N]. ∈ M There always exists such a brother since we assume Definition 2.1.1 (c) and ∑ rıˆ 1. ıˆ=1 ≤ Throughout, thee GL brother of Λ will always be denoted with the extra tilde Λ.
e 2.1. Triangular Gatzouras–Lalley-type carpets 21
Some notation The map f is indexed by i [N]. To indicate which column i belongs to in the i ∈ partition (2.1.4) we use the function
φ : 1, 2, . . . , N 1, 2, . . . , M , φ(i) := ıˆ if i . (2.1.9) { } → { } ∈ Iıˆ With this notation we can formulate the column structure (2.1.5) as
if φ(k) = φ(`) = ıˆ, then bk = b` =: rıˆ and tk,1 = t`,1 =: uıˆ. (2.1.10)
Throughout, i is an index from [N], while ıˆ with the hat is an index corresponding to a column from 1, . . . , M . We use analogous notation for infinite sequences i = i i ... { } 1 2 and ıˆ = ıˆ1ıˆ2 . . ., see Subsection 2.3.1 for details. For compositions of maps we use the standard notation f := f f ... f , i1...in i1 ◦ i2 ◦ ◦ in where i 1, . . . , N . Similarly, for products of matrices we write j ∈ { }
bi1...in 0 Ai1...in := Ai1 ... Ain := . · · di i ai i 1... n 1... n Immediate calculations give b = b ... b , a = a ... a and i1...in i1 · · in i1...in i1 · · in n n d = d a b i1...in ∑ i` ∏ ik ∏ ir , (2.1.11) `=1 · k<` · r=`+1
n where by definition ∏ aik := 1 and ∏ bir := 1. The image Ri1...in := fi1...in (R) is a k<1 r=n+1 parallelogram with two vertical sides, see Figure 2.1. We refer to bi1...in as the width,
ai1...in as the height and γi1...in as the angle of the longer side of the parallelogram Ri1...in , in other words di1...in tan γi1...in := . (2.1.12) bi1...in n i ... 1 i a
in |
... n R i 1 i ... 1 i d | γi1...in
bi1...in
2 FIGURE 2.1: The skewness of Ri1...in := fi1...in ([0, 1] )
Since direction-x dominates, Ri1...in is extremely long and thin for large n. A simple argument gives that tan γ remains uniformly bounded away from +∞. | i1...in | 22 Chapter 2. Triangular Gatzouras–Lalley-type planar carpets with overlaps
Lemma 2.1.5. There exists a non-negative constant K0 < ∞ such that for every n and every finite length word i1 ... in d i1...in K . b ≤ 0 i1...in
Proof. Since direction-x dominates, max a /b < 1, hence using (2.1.11) i{ i i} n k 1 d d di − aij max d /b i1...in | i1 | + | k | i{| i| i} < ∞. b ≤ b ∑ b ∏ b ≤ 1 max a /b i1...in i1 k=2 ik j=1 ij i i i − { }
2.1.1 Results of Gatzouras and Lalley A standard technique to give a lower bound for the Hausdorff dimension of the at- ˜ tractor Λ = i [N] fi(Λ) is to study self-affine measures νp, i.e. compactly supported ∈ measures withS support Λ satisfying e e N e ˜ 1 νp = ∑ piνp fi− , i=1 ◦ for some probability vector p = (p ,..., p ). Let be the set of all probability distri- 1 N P butions on the set [N] and be the subset when all p > 0. By definition P0 i
sup dimH νp dimH Λ. p ≤ ∈P e Gatzouras and Lalley proved that there always exists a p∗ for which the supremum is attained, furthermore p . Let ∗ ∈ P0
α∗ := dimH νp∗ = sup dimH νp. p ∈P They explicitly calculated
N N M ∑i=1 pi log pi ∑i=1 pi log bi ∑ıˆ=1 qıˆ log qıˆ dimH νp = N + 1 N N , (2.1.13) ∑i=1 pi log ai − ∑i=1 pi log ai ! ∑i=1 pi log bi
where qıˆ = ∑j ıˆ pj. This formula is a special case of the Ledrappier–Young formula, see Subsection∈I 2.3.3 for details and references. For Bedford–McMullen carpets the optimal p∗ can be given by routine use of the Lagrange multipliers method. The main result of [GL92] is that for a GL carpet the α∗ bound is sharp, i.e.
α∗ = dimH Λ.
In [GL92] Gatzouras and Lalley also gave ane implicit formula to calculate the box M sx dimension of their carpet. Let sx be the unique real such that ∑ıˆ=1 rıˆ = 1 (rıˆ was defined in (2.1.5)). Then dimB Λ = s is the unique real such that
N e sx s sx ∑ bi ai − = 1. i=1 2.1. Triangular Gatzouras–Lalley-type carpets 23
Again, equality of dimH Λ and dimB Λ is highly atypical. It holds if and only if the α∗- dimensional Hausdorff measure of Λ, denoted α∗ (Λ), is positive and finite, which H is equivalent to the conditione e e e α∗ sx ∑ aj − = 1, for every ıˆ = 1, . . . , M. j ∈Iıˆ
For Bedford–McMullen carpets Peres showed in [Per94] that α∗ (Λ) = ∞ when H dimH Λ < dimB Λ. e 2.1.2e Separatione conditions In our main results, we assume different extents of separation for the parallelograms fi(R), recall Figures 1.6 and 1.9. This will be considered in Subsection 2.1.2. In Sub- section 2.1.2 we consider separation conditions for which are actually conditions H about the extent of separation of the column structure.
Separation of the cylinder parallelograms Definition 2.1.6 (Separation conditions for a shifted TGL IFS ). We say that F satisfies the rectangular open set condition (ROSC): recall Definition 1.2.1. • F each column independently satisfies the ROSC if for every ıˆ [M] and • ∈ k, ` we have f (U) f (U) = ∅. In other words, if the interior of two first level ∈ Iıˆ k ∩ ` cylinders intersects, then they are from different columns.
satisfies the transversality condition if there exists a K1 > 0 such that for • F n every n and words (i ... i ), (j ... j ) 1, . . . , N with φ(i ) = φ(j ) for k = 1 n 1 n ∈ { } k k 1, . . . , n and i = j (φ was defined in (2.1.9)), we have 1 6 1 proj (int(R ) int(R )) < K max a ... a , a ... a . (2.1.14) | x i1...in ∩ j1...jn | 1 · { i1 · · in j1 · · jn } Given two finite words i ... i and j ... j , i = j , the angle of the two corre- 1 n 1 n 1 6 1 sponding parallelograms Ri1...in and Rj1...jn can be defined as the angle between their non-vertical sides. The transversality condition ensures that any such pair of paral- lelograms in the same column have either disjoint interior or have an angle uniformly separated from zero. Observe that this definition of transversality coincides in the diagonally homoge- neous case with the one in [BRS]. In [BRS, Section 1.5] a sufficient condition for the transversality condition was given. Namely, the authors introduced a self-affine IFS in R3 which is (in our setup) S N b 2 := Si(x, z) := ( fi(x), Ti(z)) , (x, z) [0, 1] R, S i=1 ∈ × n o where f N bwas definedb in (2.1.1) and T : R R is given by { i}i=1 i → a b N := T (z) := i z + i . T i b · d i i i=1 The relevance of the IFS is that T 24 Chapter 2. Triangular Gatzouras–Lalley-type planar carpets with overlaps
ξ ) (` A i
Ti(z) T
i
` z
(1, 0)
ai di Ti(z) := z + bi · bi
FIGURE 2.2: The IFS , where z and (1, z) are identified T
tan γi1...in = Ti1...in (0). (2.1.15)
Indeed, from the definition (2.1.12) of tan γi1...in and formula (2.1.11) it immediately follows that n ` 1 di1 di` − aik tan γi1...in = + ∑ ∏ = Ti1...in (0). bi1 `=2 bi` · k=1 bik Using the same argument as in the proof of [BRS, Lemma 1.2] we obtain that
Lemma 2.1.7. If satisfies the Strong Separation Property (that is S (Λ) S (Λ) = ∅ if S i ∩ j i = j and Λ is the attractor of the IFS ) then the transversality condition holds. 6 b S b b b b The nextb lemma gives a different,b easy-to-check sufficient condition for transver- sality.
Lemma 2.1.8. Let := (k, `) : k, ` , k = `, R R = ∅ , where A denotes Pˆ ∈ Iˆ 6 k◦ ∩ `◦ 6 ◦ the interior of a set A. Moreover, we introduce
dk ak sk := , rk := , r∗ := max rk, bmin := min bk and s := min min sk s` . 1 ˆ M bk bk 1 k N 1 k N ∗ (k,`) ˆ | − | ≤ ≤ ≤ ≤ ≤ =≤∅ Pˆ6 ∈P Assume that 1 r∗ s bmin s > 2 or equivalently ∗ > r∗. ∗ bmin · 1 r∗ 2 + s bmin − ∗ Then the transversality condition holds. In particular, in the diagonally homogeneous case transversality holds if
a d < ∗ , (2.1.16) b 2 + d ∗ where d := min min dk d` . 1 ˆ M ∗ (k,`) ˆ | − | ≤ =≤∅ Pˆ6 ∈P Proof. Using that R [0, 1]2 we obtain that d < 1. Hence k ⊂ | k| 1 sk . (2.1.17) | | ≤ bmin 2.1. Triangular Gatzouras–Lalley-type carpets 25
For an m 1, . . . , M let Σ := j Σ, j . The transversality condition holds ∈ { } m { ∈ 1 ∈ Im} if there exists c > 0 such that for every n, for all m M with = ∅ and ≤ Pm 6
for all i, j Σ with (i1, j1) m, we have: γi n γj n > c. (2.1.18) ∈ ∈ P | − |
It follows from (2.1.15) that (2.1.18) holds whenever for all such pair of i, j and for all n n ` 1 ` 1 s s s − r s − r i1 j1 ∑ i` ∏ ik j` ∏ jk | − | − `=2 · k=1 − · k=1 ! is greater than the same positive constant uniformly. However by (2.1.17) this holds if 1 r s > 2 ∗ . ∗ b · 1 r min − ∗
Separation of the columns We will also need some separation conditions for the column structure which are represented by separation properties of , recall (2.1.6). H The symbolic spaces for and are F H Σ := 1, . . . , N N and Σ := 1, . . . , M N . { } H { } The natural projections form Σ Λ and Σ Λ are Π and Π respectively, → H → H H see Subsection 2.3.1 for details. Whenever we are given a probability vector p on 1, . . . , N , we always associate to it another probability vector q on 1, . . . , M such { } { } that qıˆ := ∑ pj. (2.1.19) j ∈Iıˆ Slightly abusing the notation we write for both the set of the probability vectors P0 of positive components on 1, . . . , N and 1, . . . , M . The Bernoulli measure pN on { } { } 1 Σ is denoted µp and its push forward is νp = Π µp = µp Π− . Analogously for µq ∗ ◦ and νq.
Definition 2.1.9 (Separation conditions for ). We say that satisfies H H Hochman’s Exponential Separation Condition (HESC) (see [Hoc14, p. 775]) if • there exist an ε > 0 and n ∞ such that for k ↑
hı(0) h(0) , if hı0(0) = h0(0); ∆n := min | − | ı, 1...M n ∞, otherwise. ∈{ı= } 6
ε nk we have ∆nk > e− · . Here h0 denotes the derivative of the function h. Weak Almost Uniqe Coding (WAUC) if for all Bernoulli measures µ there exists • q Σ (may depend on q) for which BH ⊂ H 1 µq( ) = 0 and for every ıˆ Σ :#(Π− Π (ıˆ) ) = 1. BH ∈ H \BH H H \BH
Almost Uniqe Coding (AUC) holds if for every Bernoulli measure µq and for µq-a.e. 1 ıˆ Σ :#Π− Π (ıˆ) = 1. ∈ H H H 26 Chapter 2. Triangular Gatzouras–Lalley-type planar carpets with overlaps
No Dimension Drop (NDD) if for all push forward measures νq = (Π ) µq • H ∗ M ∑ıˆ=1 qıˆ log qıˆ dimH νq = − M . ∑ qıˆ log rıˆ − ıˆ=1 The following implications hold between these conditions
HESC = NDD WAUC. (2.1.20) ⇒ ⇐⇒ HESC = NDD follows from Hochman’s work [Hoc14, Theorem 1.1]. AUC im- ⇒ plies NDD from Feng–Hu [FH09, Theorem 2.8 and Corollary 4.16], but we do not know if the reverse direction NDD = AUC holds or not. Feng informed us [Fen19, ⇒ Corollary 4.7] that he can prove the equivalence NDD WAUC for ergodic mea- ⇐⇒ sures. This result just appeared on the arXiv. However, we use it only for Bernoulli measures. For completeness, we give our own complete (much simpler) proof of NDD WAUC for Bernoulli measures in AppendixB. ⇐⇒ The set of translations (u ,..., u ) defining for which HESC does not hold U 1 M H is small. It is stated in [PS, Proposition 2.7] that it essentially follows from [Hoc15, Theorem 1.10] that the Hausdorff and packing dimension of is M 1, in particular U − has 0 M-dimensional Lebesgue measure. Moreover, [Hoc14, Theorem 1.5] states U that if the parameters (r ,..., r , u ,..., u ) defining are all algebraic, then HESC 1 M 1 M H does not hold if and only if there is an exact overlap, i.e. ∆n = 0 for some n.
2.2 Main results
We now state our main results for the Hausdorff dimension of shifted TGL carpets in Subsection 2.2.1, the box dimension in Subsection 2.2.2 and discuss diagonally homo- geneous carpets in Subsection 2.2.3. For a discussion on generalizing towards nega- tive entries in the main diagonal, see Subsection 2.7.5.
2.2.1 Hausdorff dimension
For any vector c = (c1,..., cK) with strictly positive entries and a probability vector p = (p1,..., pK) we write K pi c p := ∏ ci . h i i=1 When no confusion is made, we suppress p and write c = c . Throughout, we h i h ip use this notation for the vectors a = (a1,..., aN), b = (b1,..., bN), p = (p1,..., pN), N = (N1,..., NM) and q = (q1,..., qM), where q is derived from p via (2.1.19). Using this notation let us denote the function on the right-hand side of (2.1.13) by
log p p log b p log q q log q log p log q D(p) := h i + 1 h i h i = h i + h i − h i. (2.2.1) log a − log a log b log b log a h ip h ip h ip h i h i Theorem 2.2.1 (Upper bound). Regardless of overlaps, for any shifted triangular Gatzouras– Lalley-type planar carpet Λ
dimH Λ sup D(p) =: α∗. ≤ p ∈P Furthermore, there always exists a p for which D(p ) = α . ∗ ∈ P0 ∗ ∗ 2.2. Main results 27
The proof is given in Section 2.4. Throughout, let q∗ denote the vector qı∗ˆ =
∑j ıˆ p∗j . The next theorem states sufficient conditions under which the Hausdorff ∈I dimension of a self-affine measure νp on Λ is equal to D(p).
N Theorem 2.2.2. Let p 0, µp := p and νp := Π µp. For a shifted triangular ∈ P ∗ Gatzouras–Lalley-type planar carpet Λ we have
dimH νp = D(p) if the horizontal IFS satisfies Hochman’s Exponential Separation Condition (in particular, H always holds for non-overlapping columns) and
(i) either each column independently satisfies the ROSC or
(ii) Λ satisfies transversality (see Definition 2.1.6) and the following inequality holds:
log a p log N q h i > 1 + h i . (2.2.2) log b log q h ip − h iq We remark that Proposition 2.2.10 provides a simple way to check condition (2.2.2) in the diagonally homogeneous case. Section 2.5 is devoted to the proof of this theo- rem. As an immediate corollary, we get
Corollary 2.2.3 (Sufficient conditions). Whenever a shifted TGL carpet Λ satisfies the con- ditions of Theorem 2.2.2 with replacing p and q in (2.2.2) by p∗ and q∗, then
dimH Λ = α∗.
2.2.2 Box dimension
Recall the IFSs (2.1.2) and (2.1.6) obtained by projecting to the x-axis. Recall sx H M sx H F N s˜x was defined so that ∑ıˆ=1 rıˆ = 1 and let s˜x be the unique real such that ∑i=1 bi = 1. Furthermore, introducee s := dimB Λ = dimB Λ . H H H Since Λ is a self-similar set, s is well defined. If Λ is a TGL carpet then s = sx, H H e H otherwise s sx. The affinity dimension dimAff of Λ can be deduced from the H ≤ result of Falconer–Miao [FM07, Corollary 2.6] together with the description in [BRS, Subsection 1.3] and the fact that direction-x dominates: dimAff Λ = sA is the unique real such that N min s˜x,1 sA min s˜x,1 ∑ bi { }ai − { } = 1. (2.2.3) i=1
N sA 1 In particular, if s˜x < 1 then sA = s˜x, otherwise sA solves ∑i=1 biai − = 1. So sA only depends on the main diagonals (bi, ai), but not on the off-diagonal elements di. So, the affinity dimension of a shifted TGL carpet Λ is and its GL brother coincide. The following theorem gives an upper bound for dimB Λ, which can be strictly smaller than sA. It was proved for diagonal iterated function systems by Feng–Wang in [FW05, Corollary 1] and also follows from Fraser’s work [Fra12, Theorem 2.4, Corollary 2.7]. Here we extend its scope to triangular IFSs. In a different context, Hu [Hu98] studied a related problem, where a version of Bowen’s formula deter- mines the box dimension. 28 Chapter 2. Triangular Gatzouras–Lalley-type planar carpets with overlaps
Theorem 2.2.4 (Upper bound). Regardless of overlaps, for any shifted triangular Gatzouras– Lalley-type planar carpet Λ
dim Λ = dim Λ s s , P B ≤ ≤ A where s is the unique solution of the equation
N s s s ∑ bi H ai − H = 1. (2.2.4) i=1
In particular, if Λ satisfies the ROSC, then dimP Λ = dimB Λ = s. Corollary 2.2.5 (Equality of box- and affinity dimension). For any shifted TGL carpet
s = sA s = min s˜x, 1 . ⇐⇒ H { }
Proof. Follows immediately from comparing equations (2.2.3) and (2.2.4) defining sA and s, respectively, together with the fact that ai < 1 and bi/ai > 1 for every i = 1. . . . , N. Remark 2.2.6. a) The proof of Fraser [Fra12] does not make use of any column structure (2.1.10). Hence, Theorem 2.2.4 immediately extends to an IFS of the form (2.1.1) as long as direction-x F dominates (0 < ai < bi < 1) and the ROSC holds. b) Since Λ is compact and every open set intersecting Λ contains a bi-Lipschitz image of Λ, we get that dimP Λ = dimBΛ, see [Fal90, Corollary 3.9]. Handling overlaps to calculate the box dimension is a greater challenge, since typically dimH Λ < dimB Λ and thus the usual technique of giving a lower bound by bounding the Hausdorff dimension from below does not suffice. Hence, a new counting argument was necessary. Theorem 2.2.7 (Box dimension with overlaps). For a shifted TGL carpet Λ we have dim Λ s, hence dim Λ = s, if either of the following hold: B ≥ B (i) satisfies HESC and each column independently satisfies the ROSC or H (ii) Λ is a TGL carpet, satisfies transversality and the following inequality:
log p p + log q q < s (log b p log a p), (2.2.5) − h i h i H h i − h i where p := (p1,..., pN)eande q := (eq1,...,e qM) are definede by equatione (2.2.4): s s s s s s e e e pi = bei H ai − He and eqıˆ = ∑ bj H aj− H . (2.2.6) j ∈Iıˆ e e The analogue of the following sufficient and necessary condition for the equality of the box- and Hausdorff dimensions was proved in [GL92, Theorem 4.6]. Theorem 2.2.8 (Equality of box- and Hausdorff dimension). Assume the shifted TGL carpet Λ satisfies ROSC and satisfies No Dimension Drop. Then the following three con- H ditions are equivalent,
s s dimH Λ = dimB Λ s = dimH νq ∑ aj− H = 1 for every ıˆ [M]. (2.2.7) ⇐⇒ H ⇐⇒ j ∈ ∈Iıˆ e All results for box dimension are proved in Section 2.6. 2.2. Main results 29
2.2.3 Diagonally homogeneous carpets We show how the conditions and formulas of our main results simplify in the diag- onally homogeneous case. Recall the easy-to-check sufficient condition (2.1.16) for transversality in Lemma 2.1.8. Moreover, observe that the vector p becomes the uni- form vector pi = 1/N and thus qıˆ = Nıˆ/N. A routine use of the Lagrange multipliers method gives the optimal p∗ e e e log b 1 M log b 1 log a log a − p∗ = N − N if k . (2.2.8) k ıˆ · ∑ ˆ ∈ Iıˆ ˆ=1 Thus, conditions (2.2.2) and (2.2.5) become
log p∗ p log a log N log q q log a h i ∗ < and + 1 + h i < , (2.2.9) log q∗ q log b log M log M log b h i ∗ e e respectively. If in addition, the system has uniform vertical fibres, then pi = pi∗ = 1/N also qıˆ = qı∗ˆ = 1/M. Hence, both conditions (2.2.2) and (2.2.5) become e log N log a e < . (2.2.10) log M log b
Next, we give an equivalent explicit formulation of condition (2.2.2). Let ϕ(y) := y log y and for x (0, 1) define ∈ M x 1 ϕ ∑ıˆ=1 Nıˆ R(x) := x + (r(x) 1)− , where r(x) = . M x − ∑ ˆ=1 ϕ(Nˆ )
Lemma 2.2.9. R(x) is a continuous, strictly monotone increasing function.
Proof. Continuity is obvious. It is enough to show that r(x) is strictly monotone de- creasing. Let r0 denote the derivative. Then
=:A =:B
M 2 M 2 M M 2 x x x x x ϕ(N ) r0(x) = ϕ(N ) + log N ϕ(N ) · ∑ ˆ · z ∑ }| ˆ { z ∑ ıˆ }|· ∑ ˆ { ˆ=1 ˆ=1 ıˆ=1 ˆ=1 >0 M M M M | {z } ϕ Nx ϕ(Nx) ϕ Nx ϕ(Nx) log Nx . − ∑ ıˆ · ∑ ˆ − ∑ ıˆ ∑ ˆ · ˆ ıˆ=1 ˆ=1 ıˆ=1 ˆ=1 =:C =:D
We claim that C > A and D| B, which{z will conclude} | the proof{z of the lemma.} For ≥ brevity, write y := Nx. y = 1 N = 1, otherwise y > 1. ıˆ ıˆ ıˆ ⇔ ıˆ ıˆ To show that C > A, it is enough to prove that for 1 u v ≤ ≤ ϕ(u) + ϕ(v) < ϕ(u + v). (2.2.11) 30 Chapter 2. Triangular Gatzouras–Lalley-type planar carpets with overlaps
Then a simple induction implies that ∑ ϕ(yıˆ) < ϕ(∑ yıˆ). Recall ϕ(1) = 0. The mean value theorem implies that
ϕ(u + v) ϕ(v) = u ϕ0(ξ), for some ξ (v, u + v) − · ∈ ϕ(u) ϕ(1) = (u 1) ϕ0(ζ), for some ζ (1, u). − − · ∈
Since the derivative ϕ0(y) = 1 + log y is strictly increasing and ζ < ξ, we have ϕ0(ζ) < ϕ0(ξ). This implies (2.2.11). To prove the other inequality
M M M M M 2 D B = y log y ϕ(y ) log y log y ϕ(y ) . − ∑ ıˆ · ∑ ıˆ ∑ ˆ · ˆ − ∑ ıˆ · ∑ ˆ ıˆ=1 ıˆ=1 ˆ=1 ıˆ=1 ˆ=1
We can pull out log ∑ yıˆ > 0 and divide by it. This gives