ABC of multi-fractal spacetimes and fractional sea turtles

We clarify what it means to have a spacetime fractal geometry in quantum gravity and show that its properties differ from those of usual fractals. A weak and a strong definition of multi-scale and multi-fractal spacetimes are given together with a sketch of the landscape of multi-scale theories of gravitation. Then, in the context of the fractional theory with $q$-derivatives, we explore the consequences of living in a multi-fractal spacetime. To illustrate the behavior of a non-relativistic body, we take the entertaining example of a sea turtle. We show that, when only the time direction is fractal, sea turtles swim at a faster speed than in an ordinary world, while they swim at a slower speed if only the spatial directions are fractal. The latter type of geometry is the one most commonly found in quantum gravity. For time-like fractals, relativistic objects can exceed the speed of light, but strongly so only if their size is smaller than the range of particle-physics interactions. We also find new results about log-oscillating measures, the measure presentation and their role in physical observations and in future extensions to nowhere-differentiable stochastic spacetimes.

Usually, the discussion of dimensional flow is maintained at a rather technical level but, when trying to translate mathematical properties into physics, it is customary to describe these geometries as "fractal." Then the usual questions posed when talking about fractal spacetimes are: What is a fractal? (1) What is a fractal spacetime?
How would the world look like on a multi-fractal? (3) Sometimes, at this point one feels a slight embarrassment. Establishing a set of mathematical properties for a geometry is insufficient to acquire a complete understanding of what an observer would experience in such a geometry, especially when the "observer" is an abstract Planck-size probe. Moreover, the concept of fractal has been widely used in quantum gravity, but it never has received a proper definition. Part of the reason is that there is no intrinsic definition even for the popular fractals we come across in computer graphics [31]. As the mathematician Robert Strichartz said when asked (1), "I know one when I see one" [32]. At most, we can make a list of properties we would expect the archetypical fractal should obey, marking with an asterisk optional ones: 1. A fine structure: the set has details at every scale.
2. An irregular structure: ordinary continuous differential calculus cannot be applied on the set. *3. Self-similarity. *4. A non-integer dimension (Hausdorff dimension d h , spectral dimension d s or walk dimension d w ). 5. The relation d w = 2d h /d s holds with d s ≤ d h .
Furthermore, a multi-fractal (a set whose fractal properties change with the scale) should have an additional feature: 6. Properties 1, 2 and 5 hold at any given scale in the dimensional flow.
Notice that 6 implies that dimensional flow occurs for at least two of the dimensions d h , d s , and d w (otherwise, 5 would not hold at all scales).
In the first part of this paper, we revise questions (1) and (2) and the points of the above list in the generic context of quantum gravity and of classical spacetimes with a geometry scale hierarchy. We provide a number of arguments that address question (2) and replace properties 1-6 with the following ABC of multi-scale spacetimes (items are in order of importance): A. Dimensional flow occurs with three properties: [A1] At least two of the dimensions d h , d s , and d w vary.
[A2] The flow is continuous from the infrared (IR) down to an ultraviolet (UV) cut-off (possibly trivial, in the absence of any minimal length scale).
[A3] The flow occurs locally, i.e., curvature effects are ignored (this is to prevent a false positive). B. As a byproduct of A, a non-integer dimension (d h , d s or d w , or all of them) is observed during dimensional flow, except at a finite number of points (e.g., the UV and the IR extrema). C. If, in addition, the relations d w = 2d h /d s and d s ≤ d h hold at all scales in dimensional flow, then we call the ensuing geometry a weakly multi-fractal spacetime.
Notice that A is the definition of multi-scale spacetime, while the more specific notion of multi-fractal spacetime appears only in the case property C holds. In the last section of the paper, we will comment on a stronger definition of multi-fractal spacetimes, which includes properties A-C plus a refinement of property 2: D. A geometry is a strongly multi-fractal spacetime if, in addition of satisfying A-C, it is nowhere differentiable in the sense of integer-order derivatives, at all scales except at a finite number of points (e.g., the UV and the IR extrema).
Exploiting sea turtles as an example of nonrelativistic motion, we also give an intuitive answer to question (3) with specific reference to the multifractional theory with q-derivatives [33][34][35][36][37][38][39][40][41]: on a multifractal space, a sea turtle would travel slower at microscopic scales, while it would swim faster if time were multi-fractal. Intuitively, the structure of a fractal space would hinder the motion of the poor animal due to its irregularity (an analogy is water slowly percolating through a porous rock [42]). On the other hand, a fractal time direction would be measured by clocks that tick more slowly than ordinary ones, so that a body would have more ordinary time to cover the same distance. Similarly, it is argued that a relativistic observer on a time-like fractal can reach a velocity superior to the speed of light.
If the real world has an anomalous geometry, then there must exist some critical time and length scales t * and ℓ * below which fractal properties begin to show up. In other words, geometry must be multi-scale or multifractal, as opposed to just fractal. The theory with qderivative is of this form and upper bounds on t * and ℓ * have been derived recently [40,41]. These bounds (ℓ * < 10 −19 m, t * < 10 −27 s) are about 20-30 orders of magnitude smaller than the scales involved in the seaturtle thought experiment (ℓ turtle > 1 m, t turtle > 1 s) and there would be no way to discriminate a turtle on such a multi-fractal from one on plain earth. For the same reason, super-luminal motion would be possible only for an object of size ℓ * in an experiment with characteristic time t * , i.e., below the scales of Standard-Model interactions. Therefore, even if the multi-fractional theory with q-derivatives can avoid the side effects that super-luminal travel [43] entails in Einstein gravity (existence of exotic matter [43][44][45] and quantum instability due to a thermal flux of Hawking particles [46]) or in generic Lorentz-violating setups [47], to curb the enthusiasm of Sci-Fi aficionados we already anticipate that our model, or possibly any multi-fractal spacetime in general, cannot be used as a practical base for a hyperdrive.
Questions (1) and (2) are addressed in Sect. 2. Sections 3.1 and 3.2 are a self-contained review of the theory with q-derivatives, with a new result concerning a reduction of the parameter space of measures with logarithmic oscillations. A novel discussion of the problem of presentation is given in Sect. 3.3. Since this part of the paper is somewhat technical, we summarize its content here in a few intuitive points.
(i) In order to fully define a multi-fractional theory, one must choose a frame where geometric coordinates q(x) are written down explicitly. This is not a return to a pre-relativistic view of space and time be-cause the frame choice affects the integro-differential structure of the theory, not the metric structure [39]. (ii) There is a limited number of presentation choices; we will see four below. (iii) Different presentations of the same measure correspond to different theories in the same geometric class (i.e., they show the same scaling properties). In Sect. 6.2, this feature is restated under a new perspective based on the famous Itô-Stratonovich dilemma in stochastic processes. (iv) Although they describe the same class of geometries (iii) and they are not many (ii), different presentations may have profound consequences for physical properties such as the propagation speed of bodies or elementary particles, but only in extreme situations of high energy, high curvature or small scales. (v) Therefore, even if the q-theory is invariant under Poincaré transformations on the coordinates q meant as non-composite objects, the physics is not completely independent of the choice of coordinates x used to describe the system in the frame where physical observables are extracted.
The motion of a non-relativistic and a relativistic body are studied in, respectively, Sects. 4 and 5. In Sect. 4, we also clarify the relation between measurement units and frame choice by noting that a discrimination between a fractal and a normal spacetime is possible when we can determine dimensionless quantities such as the ratio of two observables of the same kind. Section 6 is devoted only to open threads and future developments, including the effect of log oscillations and an interesting connection between presentation choice and stochastic processes on multi-scale spacetimes, which motivates the introduction of property D.

Spacetime multi-fractals in quantum gravity
Let us now examine the list of properties 1-6 and how they apply to quantum gravity or, more generally, to anomalous spacetimes.

Dimensions
To characterize a set or a geometry, we have various operational definitions of dimension. 1 In this subsection, we first revisit these definitions in the case of a continuum space, then commenting on theories with discrete structures and finally including also time. In the process, we will get a first glimpse of the type of phenomena we would experience if we lived in a multi-fractal world, and of how to detect them. 1 Mathematical statements can be found, e.g., in [31,34].

Dimensions of continuous spaces
In an ambient space with D − 1 topological dimensions, the Hausdorff dimension d h is the scaling of the volume V D−1 (ℓ) of a (D − 1)-ball with respect to the radius ℓ, d h := d ln V D−1 (ℓ)/d ln ℓ. For a set with constant dimension, while for a multi-scale set h + . . . (the actual coefficients are a bit more involved; see [42,48]). Depending on the relative size of the ball with respect to the lengths ℓ 1 , ℓ 2 ,. . . , the Hausdorff dimension will be ≃ d (1) h , d (2) h , . . . . In other words, if one tries to measure the volume V D−1 of a ball, its scaling with the radius is different depending on whether ℓ is larger or smaller than the characteristic scale at which "fractal" effects become apparent. An observer in a space with d h = D − 1 at large scales ℓ ≫ ℓ * and 0 < d h < D − 1 at small scales ℓ ≪ ℓ * can make several balls of radius R 1 + δR close to some average value R 1 ≫ ℓ * (where δR ≪ ℓ * ), submerge each ball in a container of water and measure the volume of displaced liquid, noting a distribution of volumes with average R D−1 1 and width ∼ R D−2 1 δR. Making another set of balls of average radius δR < R 2 ≪ ℓ * with the same fluctuation δR, they find an average volume R dh 2 and (for D ≥ 3 and d h ≥ 1) a narrower distribution, since 1 ≪ (R 1 /ℓ * ) D−2 > (R 2 /ℓ * ) dh−1 ≪ 1. The inequality may change direction for d h < 1 but, in any case, by comparing these dimensionless observables the experimenter realizes that they are living in a space with dimensional flow.
The spectral dimension d s is the scaling of the return probability in a diffusion process. Let K(∇) be the Laplacian of a theory on a continuum Riemannian manifold; in the standard case, K(∇) = ∇ 2 . Placing a point-wise test particle at point x ′ on a spatial geometry and letting it diffuse, its motion will obey the non- where κ is a diffusion coefficient, σ is an abstract diffusion time parametrizing the process and g is the determinant of the metric. If σ = t is the proper time of the particle or a viable global time variable, then κ is measured in m 2 s −1 . Integrating the heat kernel P for coincident points over all points of the geometry, one obtains a function P(σ) (the volume factor makes the normalization finite). Then d s := −2d ln P(σ)/d ln σ. For a set with constant spectral dimension, while for a multi-scale set the structures of K and of a generalization of the operator ∂ σ determine two or more asymptotic regimes [42]. A particle in a space with d s = D − 1 in the IR and 0 < d s < D − 1 in the UV diffuses slower in the ultraviolet. 2 Finally, the walk dimension is the scaling of the mean-square displacement of a random walker X(σ), d w := 2(d ln . For a set with constant walk dimension, In a space with walk dimension d w = 2 in the IR and d w > 2 in the UV, the ratio (d X 2 /dσ)/( X 2 /σ) ∼ 1/d w between the differential velocity of the particle (measuring its local random motion) and the total finite-difference velocity decreases in the UV. This means that the trajectory of the probe becomes less ragged than usual in the UV; a pictorial demonstration of such behavior can be found in [37].

Dimensions of discrete and combinatorial structures
In several quantum-gravity approaches, there is no fundamental continuous spacetime. Nevertheless, it is possible to generalize the above operational definitions to a discrete set and to extract, only in certain regimes, sensible multi-scale profiles for the Hausdorff, spectral and walk dimension. A proof of concept is given in [4,10,11,[25][26][27][28]. For instance, a generalization of discrete exterior calculus [52][53][54] allows one to construct Laplacians on combinatorial structures and hence define diffusion processes thereon. The spectral dimension as well as d w and d h can be computed for a class of quantum-gravity states, in particular those appearing in loop quantum gravity, spin foams, and group field theory [10,11].
The main requirement we ask in order to have dimensional flow in discrete (pre-)geometries is that there exist regimes where all three geometric indicators d h , d s , and d w are real-valued and positive. These regimes should extend from the infrared down to some effective UV scale below which discreteness or quantum effects destroy some or all of the indicators, for instance if the expectation values on the chosen states become complex. In general, the effective UV scale is determined by the choice of states. The existence of a regime where quantum geometry has well-defined dimensions translates into regularity assumptions on the quantum states [11]. Beyond these assumptions, one can plunge into a wild jungle of quantum-geometry configurations with properties completely different from classical or semi-classical spacetimes.

Dimensions of spacetimes
In a continuous spacetime with D topological dimensions, the definition of Hausdorff dimension is unchanged but for the addition of the Euclideanized time direction. Then the Hausdorff dimension of spacetime is the scaling of the volume of the D-ball. For the spectral dimension, one includes (imaginary) time in the operator K(∇), while in the walk dimension time is included in Euclideanized distances X 2 (σ) = T 2 (σ) + X 2 1 (σ) + · · · . Similar considerations apply to discrete geometric or pre-geometric structures.
In the extension of all the above definitions, one takes the time direction in Euclidean signature. This step is fairly standard when one wants to define the dimension of a geometry with Lorentzian signature. If, for any reason, one cannot or does not want to Euclideanize time, then it is necessary to consider the dimensionality of spatial slices and the time line separately, instead of the whole spacetime. The reader may adopt whichever point of view they might prefer; this does not affect the following.
The spectral dimension is sometimes regarded as a theoretical parameter useful to classify spacetimes but that does not correspond to a physical observable. Elsewhere [8,20], the author and collaborators had already the occasion to advance a different view: the spectral dimension should be a meaningful observable just like the topological and Hausdorff dimension are. In that case, however, its definition must be well posed at all stages to make sense physically: if we want d s to be a physical observable (our working hypothesis here), its definition must provide also an operational way to measure it. For instance, how can we interpret the parameter σ if time is in the operator K? Also, in certain cases the form of K is such that P (x, x ′ , σ) is not a probability and there is no well-defined underlying diffusion process at all (this is a well-known problem in transport theory with higher-order or non-local operators [42] and in quantum gravity [8,20]). Mathematically there is no issue whatsoever. If σ = t, one can enact a fictitious diffusion process with some Monte-Carlo time on the geometric or pre-geometric structure one wants to explore, let it be a continuous manifold or the graph ensembles of discretized gravity. Even when P is not positive semidefinite and the picture of a diffusing probe fails, to determine d s one only needs to consider closed paths and integrate over them with a certain measure. However, this is insufficient to characterize an operational way to physically measure d s . In two different interpretations of the spectral dimension, valid in any regime where an effective field theory can be formulated, the diffusion equation is a renormalization-group running equations depending on the IR cut-off scale k = f (σ) [8] or, alternatively, it stems from the Schwinger representation of the particle propagator and √ κσ = ℓ is a length scale determining the resolution 1/ℓ at which the geometry is probed [20]. The interpretation of the parameter σ is unimportant in the following but it is worth to mention these caveats anyway.

Fine structure and dimensional flow
The first property of the list 1-6 is that a (multi-)fractal set F should have "a fine structure." By this, one means that it is possible to find points of F at all scales of observation, no matter how deeply one zooms into the set. For continuous spacetimes (among others: asymptotic safety, Hořava-Lifshitz gravity, non-local gravity and multi-fractional spacetimes), this requirement seems trivially satisfied and not very useful. However, a careful inspection shows that it is neither trivial nor satisfied in general.
One of greatest Einstein's intuitions was that spacetime points do not have a physical meaning per se unless one attaches an event to them. A spacetime devoid of particle interactions, test particles, light rays or whatever event announced by matter is an empty mathematical construct. To make sense of the idea of "finding spacetime points at all scales," one should be able to concoct an experiment where the physical probe can be utilized at all scales. Of course, sometimes the same device can give us information on the physics at very different scales, as is the case with the Planck satellite or similar observatories of the cosmic microwave background. But, in general, we do not have a universal instrument and we need to resort to different set-ups (a telescope, a particle accelerator, . . . ) to probe the physics at different scales.
Once having defined our ideal probe as a patchwork of instruments and experiments covering all scales of interest, the problem remains to test the spacetime structure at arbitrarily small scales. Apart from obvious technical limitations we have now and probably forever after (we cannot probe the Planck scale directly, nor energies near grand unification), it may even be theoretically impossible to reach an infinite resolution, mainly because of quantum uncertainty. This is the case of asymptotic safety, where, despite the absence of any fundamental length in the theory, a minimal length appears below which one cannot separate two points by a dynamical probe [55]. 3 The plethora of theories based on discrete structures is also unaffected by property 1 because there are no details below the discreteness scale. For instance, causal dynamical triangulations are a discretization of a continuum but, for any practical purpose, one cannot trust any probing of the geometry at scales comparable with the size of the triangulation cell. Loop quantum gravity and spin foams are defined on complexes that induce a minimum physical Plancksize length in the spectra of volume operators [56]. Also, both the underlying discreteness and the combinatorial structure impose an effective UV cut-off limiting the range of scales where one can make sense of the concept of spacetime dimension, while at scales larger than the cut-off they render such dimension anomalous [10,11].
A much more important property than the fine structure is that the effective geometry must have some quantum-to-classical regime where dimensional flow takes place, otherwise one could not reach a semi-classical continuum limit where the dimension of spacetime is 4. For instance, suppose to find a dimensional flow from d s ∼ 2 in the UV to d s ∼ 4 in the IR (examples of this abound in the literature [1-4, 7, 11, 12, 30, 33]), while below the UV scale one finds a non-geometric phase where one cannot define the spectral dimension, possibly for discreteness or combinatorial effects (as in [10,11]). Then below the UV scale the geometry certainly does not show a fine structure (zooming in too much, we enter "inside" the building blocks of the theory, let them be lattice cells, tetrahedra or something else). During dimensional flow, the fundamental degrees of freedom (e.g., quanta of geometry, labeled complexes, and so on) group together into collective modes such that the notion of spacetime dimension makes sense and is measurable. When coarse graining the fundamental degrees of freedom, the resulting effective structure is most likely to be "fine," which can be tested by finding effective dynamical equations on an effective continuum. However, this test is nontrivial and few are the cases where it can be carried on [4]. Usually, the only datum we know, corroborated by a numerical or analytical study of dimensional flow through all scales from the effective UV cut-off to the IR, is that discreteness effects are present but not dominant in that interval.
From this discussion, we see that theories which have dimensional flow may or may not have a fine structure at all scales. Also, theories which do not have a fine structure do not necessarily have dimensional flow (example: a canonical second-order scalar field theory on a cubic lattice), while theories which have a fine structure can describe most boring geometries (example: any canonical second-order field theory on Minkowski spacetime). We conclude that property 1 is not adequate in the context of anomalous (quantum or classical) spacetimes, many of which are not fractal in the standard sense because they do not have a fine structure.

Irregular structure
An ordinary (multi-)fractal set F has "an irregular structure" in the sense that it cannot be described by Euclidean geometry. A Euclidean ruler would fail to measure the total length of the Western coast of Britain [57]. Clearly, in a physical context the geometry in the infrared must be "regular," so that we should consider property 2 only at the microscopic scales of a multi-fractal spacetime. In gravitational theories, geometric probes are local and curvature effects are usually ignored when one determines the dimension of space (which would be modified by curvature even in a purely classical setting, 4 according to the Seeley-DeWitt formula [58]). However, even locally there are other effects that make the geometry non-Euclidean, for instance if gravity is quantized or in the presence of a non-trivial integro-differential structure. In the first case (quantum gravity), the collective effect of quanta of geometry is to push around the probe in an anomalous way, not experienced in a classical space. Often this induces effective operators in the dynamics, which leads to the second case (multi-fractional spacetimes). A third case consists in frameworks with an underlying discrete nonregular structure, such as the complexes found in loop quantum gravity, spin foams, and group field theory.
All three cases can be realized in so many different ways that establishing the "irregularity" of a geometry is a moot point. If a spacetime shows dimensional flow locally (i.e., ignoring curvature corrections), then there must be some mechanism making it irregular. Conversely, an irregular spacetime does not have to be a (multi-)fractal unless it also has a fine structure, just like the rugged surface of a rock may not be a fractal (if we zoom in, we may discover that locally it is smooth).
In Sect. 6, we will consider a more precise characterization of irregularity as one of the requirements to reproduce certain microscopic properties of stochastic processes on fractals.

Self-similarity
Self-similarity and self-affinity are what defines all deterministic fractals. A deterministic fractal F = i S i (F ) is the union of the image of some maps S i which take the set F and produce smaller copies of it (possibly deformed, if the S i are affinities). Not all fractals are deterministic, yet they are fractals indeed; sets with similarity ratios randomized at each iteration are of this sort and they are called random fractals.
Since self-similarity and self-affinity are shown by a huge but non-exhaustive class of fractals, it is clear that we cannot use them to characterize spacetimes in an efficient way. The standard Poincaré transformations x ′ µ = Λ µ ν x ν + a µ are affinity maps and the dynamics of a covariant field theory on Minkowski spacetime is selfaffine. Yet, it is not a fractal because it has no irregular structure. On the other hand, multi-fractional spacetimes with ordinary derivatives have dimensional flow but they are neither self-similar nor self-affine [3,34] and for this reason they can be used only as effective models [39].
Any theory of particle physics and quantum gravity worth of this name is both under analytic control and potentially predictive provided symmetries are enforced. There are no known exceptions to this rule. Whatever these symmetries are (diffeomorphism invariance, conformal invariance, supersymmetry, modular invariance, and more), they constitute a guiding principle and the backbone of the theory; they may or may not give rise to dimensional flow, which is an accidental property of geometry. This point of view is not very different from what happens when a mathematician wants to construct a fractal: first some maps are defined and then the geometry of the set is studied. However, the connection between symmetry and fractality is much more tenuous in physics and symmetry takes precedence over virtually anything else.

Non-integer dimension
One of the most popular features of fractals is that they have non-integer dimensions. For instance, the Hausdorff dimension d h of the middle-third Cantor set is equal to the capacity d c := − ln N/ ln λ = ln 2/ ln 3 ≈ 0.63. Each iteration is made of N = 2 copies rescaled by λ = 1/3. However, there are many fractals with integer dimension, e.g., the Mandelbrot set and its boundary (both with d h = 2). We refer the reader to [31,34] for definitions and more examples and counter-examples. Conversely, a set with integer dimension is a fractal only if it has an irregular structure. For instance, we can tell apart the string world-sheet from the boundary of the Mandelbrot set because there is no Virasoro algebra of operators acting on the latter [30].
On the other hand, if we have a continuous dimensional flow we expect to sample over all values of the dimension between the UV and IR terminal points, which implies that the dimension is integer only at a finite number of scales, from a minimum of one (in the in- A back-of-the-envelope argument shows the existence of a relation between the dimensions d h , d s , and d w of a fractal set. Let us denote by ℓ a length scale, be it the average displacement X 2 ∼ ℓ or the linear size of a volume V ∼ ℓ dh . We saw that d s is defined as the scaling of the return probability, Eq. (5). The latter is a probability per unit volume, so that it scales as an inverse volume. Then which implies This formula holds independently of the interpretation of the parameter σ and the only assumption one makes is about the volume scaling of the return probability. This assumption is what characterizes fractals or non-multi-scale sets such as ordinary manifolds. Furthermore, on fractals the spectral dimension is always smaller than the Hausdorff dimension and d w > 2 (subdiffusion). The specific type of sub-diffusion on fractals is called labyrinthine [59] because the probe is hindered by "obstacles" and "dead ends" in the geometry. Diffusion on fractals can be approximated by a diffusion equation with fractional differential operator ∂ β σ , an xdependent diffusion coefficient and a friction term [60].
We will not use the form of this diffusion equation as a criterion to define a fractal as it seems too restrictive and not as robust as the above model-independent arguments. However, we will come back to this point in Sect. 6.
There are various examples of multi-scale processes or geometries similar to fractals but which do not obey Eq. (8) or for which d w < 2. Lévy processes are an example well known to mathematicians (see [42] and references therein). In the case of spacetime geometries, certain non-commutative and non-local spacetimes have a spectral dimension which grows in the UV and d s > d h [18,23]. Also the spacetimes of the multi-fractional theories with weighted and standard derivatives (Hermitian dual to each other) are not multi-fractals, since d w = 2D/d s and d h = D [37]. On the other hand, the theory with q-derivatives (discussed below) respects Eq. (8) [37] and so does asymptotically safe quantum gravity (where d h = D) [61] (see also [8]). We do not have data about the walk dimension in the other approaches mentioned in the introduction.

ABC of spacetime fractals
From what seen in this section, the most general and powerful characterization of anomalous spacetimes is dimensional flow. In theoretical physics, having a fine or irregular structure is not so much important as having a set of fundamental symmetries in action, but the connection between symmetries and dimensional flow is rarely immediate (exceptions are multi-fractional theories). Also, relation (8) is not obeyed in several cases.
Since all the main properties of fractals are either violated or modified in quantum gravity and spacetime theories with anomalous geometry, it is highly recommended to shift the attention to the more practical notion of multi-scale spacetimes, defined with property A (and the ancillary feature B) in the introduction. If Eq. (8) and d s ≤ d h hold, then a multi-scale spacetime is also multi-fractal (property C).
In Fig. 1, we depict a first snapshot of the landscape of multi-scale theories with dimensional flow. Apart from the sub-class of multi-fractal spacetimes, some quantum-gravity frameworks are also indicated: asymptotically safe quantum gravity (AS) and loop quantum gravity (LQG). While AS realizes multi-fractal spacetimes in all known cases (see, e.g., [8]; but the situation may change as the framework evolves), the type of geometry produced by LQG depends on the states chosen in the expectation values of the operators defining the dimensions. In [11], one can see several examples of states corresponding to a multi-fractal geometry (region A), to multi-scale but not multi-fractal quantum geometries, and to highly quantum geometries which cannot be classified by conventional geometric indicators. The latter case is the corner "?" lying outside the multiscale landscape. A third class of scenarios is the one of multi-fractional spacetimes, which are not necessarily of quantum gravity. 5 Of the four multi-fractional models proposed (with ordinary, weighted, q-and fractional derivatives), two do not realize multi-fractal spacetimes (theories with ordinary and weighted derivatives), one has not been analyzed in full detail yet (theory with fractional derivatives) and the fourth lives on multifractal spacetimes. Region B includes the theory with q-derivatives and probably also the one with fractional derivatives.
The theory with q-derivatives is useful to describe the renormalization-group flow of asymptotic safety in an alternative way [36]; this is represented by an overlap between the AS set and the multi-fractional one. Used as effective descriptions of geometry, multifractional models can reproduce the dimensional flow of other theories. This connection has not been shown for LQG and is indicated here with the intersections "?" inside the landscape. Many other well-studied multiscale theories are not shown either (including noncommutative spacetimes and Hořava-Lifshitz gravity, both of which do have an overlap with multi-fractional models [17,36], and dynamical triangulations), because the walk dimension has not been calculated yet and we are presently unable to verify that d w = 2d h /d s . However, the overwhelming majority of quantum-gravity cases have d s ≤ d h , which would put them inside the multi-fractal region if property C were confirmed.
In the following, we will concentrate on the multifractional theory with q-derivatives to illustrate what we could expect to see in a multi-fractal spacetime.
3 Multi-fractional theory with q-derivatives

Sketch
Multi-fractional theories are realizations of anomalous geometries rather than frameworks for quantum gravity, but they can be used also as models describing effective regimes of other proposals and, as in this paper (where gravity does not play any role), to clarify what we mean by fractal spacetimes. Figure 1 The landscape of multi-scale theories with anomalous geometry.
We begin with a brief review of the multi-fractional theory with q-derivatives; more information can be found in [39,41]. In D topological dimensions (D = 1+3 in our case), the dynamics is defined with respect to a geometry endowed with characteristic scales. In practice, one takes their favorite action here the first entry in the left-hand side is the measure weight in the action, φ i are some generic degrees of freedom and we are ignoring gravity) and makes the replacement x µ → q µ (x µ ) everywhere (including the derivatives, ∂ x → ∂ q(x) , hence the name of the theory). The profiles q µ (x µ ) are called geometric coordinates and the theory is invariant under the non-linear q-Poincaré transformations where Λ µ ν are the usual Lorentz matrices and a µ is a constant vector. The system can be written in two different ways, where v = det |dq µ /dx µ |. The left-hand side is the frame described by the geometric coordinates q, called the integer picture, and it is the starting point to formulate the theory. The right-hand side is the frame, called the fractional picture, where the x-dependence of the geometric coordinates q(x) is manifest.
To complete the definition of the theory, we still need two data: the choice of profiles q(x) and the choice of frame. If we want our continuous spacetime to change dimension with the scale, we must be able to tell the difference between "large" and "small" distances, or between "early" and "late" times. For this purpose, we can introduce at least one characteristic length ℓ * and one characteristic time t * in the choice of q(x) (a more general hierarchy is discussed in [33,35]). It turns out that, in D = 1 dimension, dimensional flow can be achieved by exactly the same type of measure q(x) ∼ |x| α (called fractional measure, where 0 < α < 1) one would obtain if one approximated a fractal dust on a continuum line [35]. To get a multi-fractal, it is sufficient to add several power-law contributions |x| α l with different α l , each multiplied by a characteristic scale ℓ (l) * .
For a binomial measure (just two terms, one scale ℓ * ) and requiring to have a standard geometry in the infrared (α 1 = 1, α 2 = α), we have Generalizing to D − 1 spatial directions, one takes the Cartesian product of D − 1 anomalous lines and the multi-fractional spatial measure dq * (x 1 ) dq * (x 2 ) . . . dq * (x D−1 ), possibly with different exponents α i [35]. This Ansatz for a factorizable measure is not only sufficient for our aims but it may also be necessary for technical reasons explained elsewhere [62].
To complete the definition of the measure, we take a copy of (11) also in time, In the following, we will consider an isotropic configuration where all spatial α i = α and ℓ i * = ℓ * are the same, so that the parameters of the theory in position space are ℓ * , t * and the two fractional exponents α 0 and α.
Having chosen the profiles q µ (x µ ), let us consider the choice of frame. We must select in which picture physical observables are computed. On an ordinary manifold, the properties of clocks, rods and detectors are the same independently of the scale at which measurements are taken. In contrast, multi-fractional spacetimes are a framework where physical measurements are performed with instruments which do not adapt with the observation scale even if the geometry does [36]. This adaptation is encoded in the structure of the fractional coordinates (i.e., of the integration measure and of differential operators), where characteristic time, length and energy scales appear.
Specifying units for the coordinates clarifies the point. In c = 1 units, time and spatial coordinates scale as [x µ ] = −1 and so do the characteristic scales ([ℓ * ] = −1 = [t * ]) and the geometric coordinates ([q µ (x µ )] = −1). However, in the ultraviolet the variable dependence of q µ has an anomalous scaling ≃ [|x µ | αµ ] = −α µ , which implies that q-clocks and q-rods adapt with the scale of the experiment. Since our actual clocks and rods are non-adaptive rigid apparatus, observables should be compared with experiments in the fractional picture. A more detailed discussion can be found in [41] but in Sect. 4 we will make an important observation so far overlooked in the literature. One may be confused by the above argument relying on the anomalous scaling of the variable part of q. However, even if x and q have the same measurement units [x] = −1 = [q] exactly, it is possible to recognize a standard spacetime from an anomalous one by measuring dimensionless quantities such as ratios of observables.

Multi-fractal properties
We now show, by recalling prior results and finding new ones, how the theory knows the ABC of spacetime fractals. We will do this by a heuristic computation of the dimensions for a Euclideanized no-scale geometry with isotropic profiles q µ ∝ |x µ | α , where α 0 = α. (The dimension of space is obtained by ignoring the time contribution.) Exact results and the multi-scale case can be found in [35,37,42]. The Hausdorff dimension of a D-ball with radius ℓ centered at the origin is V D (ℓ) ∝ L D (ℓ), where L 2 (ℓ) = µ [q µ (ℓ)] 2 . From that, one finds d h = Dα. The spectral dimension descends from the diffusion equation where q β ∝ σ β , and reads d s = −2d ln P/d ln σ = −2βd ln P/d ln q β (σ) = Dβ. The walk dimension d w = α/β is read from X 2 ∝ q 1/α β . Then Eq. (8) holds for any β. If β ≤ α, then d s ≤ d h and the geometry is fractal; one has equality if σ is a genuine diffusion time or if one fixes the ambiguity parameter β to be the average fractional charge α [37]. 6 In Sect. 2.6, we mentioned that the effective diffusion equation on fractals is modified not only in the spatial part (Laplacian, diffusion coefficient, friction terms) 6 The anisotropic case where some or all αµ are different is tricky because the components of the random walk X 2 = X 2 0 + X 2 1 + · · · would have inhomogeneous scaling and each direction should be considered separately. This complication is responsible of the fact that the dimension of the Cartesian product of fractals may not coincide with the sum of their dimensions [31]. but also in the time part, via a fractional diffusion operator. Although we have not included these modifications as part of the definition of fractals (to the best of our knowledge, this type of diffusion equation is only an empirical tool to describe transport on fractal media), they teach us that sub-diffusion on a space-like fractal may come from a diffusive process parametrized by an anomalous clock σ. From what we know about diffusion in multi-scale spacetimes with q-derivatives [37], we recognize a diffusion operator ∂ q β (σ) in Eq. (14) which is not fractional but it is anomalous nevertheless. Moreover, expanding κ∇ 2 q(x) in x coordinates we find both an effective space-dependent diffusion coefficient ∼ κ/(∂ x q) 2 and a first-order friction term. All the ingredients of fractal diffusion are here, albeit modified with respect to the phenomenological models of [60].
The geometry of the q-theory is a random fractal, namely, a fractal endowed with symmetries whose parameters are randomized each time they are applied over the set [35,63]. To get a deterministic fractal where the symmetry parameters are fixed (the Cantor sets and the Koch curve are examples), it is sufficient to include logarithmic oscillations of the coordinates in Eqs. (11) and (12) [35,64]. Then, for each direction and in dimensionless units, one replaces q(x) = l q α l (x) with q log (x) = l q α l (x)F ω l (ln |x|), where ω l is a frequency parameter and F ω l = 1 + A cos(ω l ln |x|) + B sin(ω l ln |x|). The fractal F = ⊗ µ F µ represented by a measure with only one frequency ω > 0 is given, for each direction, by the union of N copies of itself rescaled by a factor λ ω = exp(−2π/ω) at each iteration. Since the capacity of F µ is equal to the Hausdorff dimension and reads d c = − ln N/ ln λ ω = d h = α, the number of copies is This formula is implicit in the results of [35] but here we recognize a new element that shrinks the parameter space of the theory considerably: since N is a positive integer, then ω can only take the irrational values For α = 1/2 and N = 2, 3, . . ., we have λ ω = 1/N 2 and N = 2 , ω 2 ≈ 4.53 , λ ω = 1 4 , N = 3 , ω 3 ≈ 2.86 , λ ω = 1 9 , . . . These are sets that differ from one another only by the values of the parameters λ 1,2 and a in the similarity maps S 1 (x) = λ 1 x and S 2 (x) = λ 2 x + a; the middlethird (or ternary) Cantor set is only one member of the class, with λ 1 = λ 2 = 1/3 and a = 2/3.

The problem of presentation
In this subsection, we analyze the effect of changes in the presentation of the geometric coordinates q µ (x µ ).
By construction, the symmetries of the system are the q-Poincaré transformations (9), not the ordinary ones, and the laws of physics are invariant accordingly. However, physical observables are determined in the fractional frame and, therefore, they are not invariant under these transformations. Suppose one wishes to measure the distance ∆x of two points A and B in a sheet of paper. If the paper is charted by a Cartesian system, then the distance is given by the Euclidean norm Then we make a coordinate transformation is a function of the new coordinates. For instance, going to polar coordinates {x 1 , x 2 } → {̺, θ} conveniently centered at x A , one has ∆x = r. The observed value of the distance is insensitive to the coordinates we choose to represent ∆x with. In the theory with q-derivatives, we repeat exactly the same discussion in the fractional picture, which is one of the coordinate frames {x} where the distance ∆x is calculated. However, to each of these fractional frames we must associate an integer frame described by geometric coordinates. Thus, the Cartesian fractional frame {x 1 , x 2 } is mapped into the integer frame {q 1 (x 1 ), q 2 (x 2 )} and, after inverting to x i = x i (q i ) (assuming it possible, which is not always the case) the Euclidean norm (17) is mapped into some complicated expression ∆x(q i A , q i B ) which differs from the geometric Euclidean norm ∆q : Below we will calculate the difference and encode it in functions X i . But now we redo the mapping to geometric coordinates starting from polar fractional coordinates. The new integer frame is {q r (r), q θ (θ)}, where the relations between q r and the q i are q 1 = q r cos q θ and q 2 = q r sin q θ .
Having recalled the rather self-evident fact that arbitrary changes of chart {x µ } → {x ′ µ } modify q(x), the question is: On which chart are Eqs. (11) and (12) represented? In the example of the paper sheet, is Eq. (11) the form of q in the integer frame {q 1 (x 1 ), q 2 (x 2 )} based on Cartesian coordinates {x 1 , x 2 } or the form of q in the integer frame {q 1 (r), q 2 (θ)} based on polar coordinates {r, θ} (so that q 1 (r) = r + (ℓ * /α)(r/ℓ * ) α ), or something else? Ordinary Poincaré invariance is violated by factorizable multi-scale measures. A change of presentation such as a translation, a rotation of the coordinates or an ordinary Lorentz transformation modify the size of the multi-scale corrections X and T defined below and one realizes that different choices of the fractional frame lead to a different theory in the integer frame. Clearly, q 1 (r) = [q 1 (x 1 )] 2 + [q 2 (x 2 )] 2 due to the nonlinear terms in expressions such as Eq. (11).
Starting from [35], the tacit assumption has been that Eqs. (11) and (12) are based upon the Minkowski frame where all coordinate axes are orthogonal. So far this assumption has not been discussed in detail. We fill this gap here.
First and foremost, a change of presentation changes the theory (i.e., the magnitude of the corrections X and T ) but not its qualitative features. It is well known that inequivalent presentations leave the anomalous scaling of the measure and the dimension of spacetime untouched [34,35]. Therefore, multi-fractional scenarios are robust across different presentations. Picking a presentation allows us to make predictions which will change in another presentation, but not by much.
Second, the choice of the Cartesian or Minkowski fractional frame is not so restrictive as it might seem. In the physical examples studied in the literature, the observations studied in the theory with q-derivatives involved: (a) the decay rate of the muon [40,41]; (b) the Lamb shift in the spectrum of hydrogenic atoms [40,41]; (c) the cosmic-microwave-background black-body spectrum [65]; (d) the cosmic-microwavebackground temperature spectrum [65]. In the theory with weighted derivatives, we studied (b ′ )=(b), (c ′ )=(c) and (e ′ ) the fine-structure constant determined from the light of quasars [66]. In (a), the multi-scale correction is only time dependent and t is the muon lifetime. In (b), (c), and (c ′ ) the multi-scale correction is energy-or temperature-dependent and, as we will argue below, this poses no problem of presentation. In (d), the spectrum is written as a function of the absolute value |k| of comoving spatial momentum; since we use Cartesian momenta {p 1 (k 1 ), p 2 (k 2 ), p 3 (k 3 )}, the expression in the fractional frame is in terms ofk := [p 1 (k 1 )] 2 + [p 2 (k 2 )] 2 + [p 3 (k 3 )] 2 = |k| + · · · but, to leading order in the multi-scale correction, it is not different from what one would have obtained using a profile p(|k|). In (b ′ ), the multi-scale correction depends on the characteristic time t of the electromagnetic processes involved in the Lamb shift. In (e ′ ), the correction depends on the cosmic time t of emission of light of distant objects since the big bang. All these settings are characterized by an effective one-dimensional multiscale correction, either of rest-frame energies or of a well-defined time variable. There is not much arbitrariness here and the Cartesian or Minkowski chart fits the purpose.
However, there is one last bit of ambiguity which deserves our attention: a translation Consider one dimension, the spatial interval ∆x = |x B −x A | between two points A and B and its geometric analog for a binomial measure: Extending this result to time intervals ∆t := |t B − t A |, one has a similar expression: If we do not fix the presentation (i.e.,x andt), we place ourselves in a quandary. Take for definiteness x A = 100 m and x B = 200 m. Then, for α = 1/2 (realistically found in the theory) and ℓ * = 1 mm (unrealistically large), the constraint |X | ≪ 1 holds for any presentation (consistently with ℓ * /∆x ≪ 1; this is the universal behavior of different presentations announced above) but the actual value and sign of X change. Conversely, fixing x B − x A andx to some random numbers while varying with respect to x A leads to a similar spectrum of values of both signs. A multi-fractional theory without a prescription on the measure presentation is not defined completely and, therefore, would not be predictive even if we knew the scales ℓ * and t * . We will call null presentation the one withx = 0 andt = 0 and denote by other possible presentations, there are three such that ) and, therefore, do not depend on translations of the coordinate frame. In [40,41], we made the natural identification ofx with the starting point of the experiment, in this case x A . Setting instead x = x B would lead to exactly the same result but with X → −X . Let us call these prescriptions initial-point presentation and final-point presentation, respectively. The most uninteresting case is the symmetrized presentationx = (x B + x A )/2, which yields X ≡ 0 and a trivial theory. Calling q ± the geometric coordinates in the initialand final-point presentation, we obtain the expressions ∆q ± (x) = ∆x|1 + X ± | and ∆q ± (t) = ∆t|1 + T ± |, where The sign depends on the choice between initial-point presentation (+) and final-point presentation (−). In D dimensions, for each spatial direction x i one has a copy of X i = X (x i ).
Let us now discuss the problem of presentation in Fourier space and consider the momentum (13) dual to the binomial geometric coordinateq * (t). Including also an arbitrary constantĒ, one has the energy measurē The energy in this equation is simply E = k 0 , the time component of the D-momentum, and can take either positive or negative values. The parameterĒ can be identified with the energy of the ground state of the system. The natural choice in quantum field theory is E = 0 [40,41], but in certain phenomenological situations one might want to measure the energy with respect to the ground state. For instance, a constraint on the fundamental energy scale E * was found in [40] by asking that effects of the multi-fractal geometry be smaller than the experimental error δE on the 2S − 2P Lamb shift ∆E in the hydrogen atom: WhenĒ = 0, this lower limit is E * > 35 MeV for generic α 0 and E * > 450 GeV for α 0 = 1/2. If we choose insteadĒ = E 1S , the last factor in Eq. (23) changes from E 2S ≈ −3.4 eV to E 2S → E 2S − E 1S ≈ 10.4 eV and the bound increases about three times. From this example, we can imagine that, in general, only a positive detection of multi-scale effects would be able to rule out one presentation instead of another. This expectation is confirmed by the following analysis.

Fractional sea turtles
To illustrate the effects of both the presentation choice and the picture selection (i.e., the problem of using non-adaptive instruments in a scale-dependent environment), let us consider a (1 + 1)-dimensional nonrelativistic experiment where a wildlife ranger wants to check whether they live in a smooth Minkowski spacetime or, more interestingly, in a multi-scale spacetime with binomial measure q * (x −x). To this purpose, they observe an adult sea turtle of size L ∼ 1 m that, after laying her eggs on a beach, enters the waters at point x A at time t A (event A) and reaches a buoy at point x B at a later time t B (event B). The ranger knows that, according to the theory, measurements are performed in the fractional picture, which is the right-hand side of Eq. (10). Here, coordinates x are non-composite (i.e., their scaling is one and the same at all scales) and the action S = d 2 x v(x) L of the system has a non-trivial measure weight v(x) = det |dq µ /dx µ | breaking Lorentz invariance and deforming kinetic terms. The observer will use clocks and rods which do not adapt with the scale of the observed object, a standard analog wrist watch measuring time intervals ∆t and a rigid rod measuring distances ∆x in metric units. Recall that coordinates are "adaptive" or "non-adaptive" depending on whether they are composite objects or not. We call the ranger consulting the "fractional twatch" a fractional observer O x , to distinguish them from an integer observer O q which would use an adaptive "integer q-watch" measuring intervals ∆q * (t). For O q , spacetime is ordinary Minkowski when expressed in terms of the composite (adaptive) coordinates q and the action is S = d 2 q L. In other words, an integer observer is an ordinary observer in an ordinary, nonmulti-fractal world.
Below, we omit the subscript * and use the symbols v x and v q to indicate the velocities not to be confused with the measure weight v = ∂ x q of the rest of the literature.

Preliminaries
The ranger notices that their t-watch ticks 20 times between the events A and B, so that ∆t = |t B − t A | = 20 s. On the other hand, the geometric time interval ∆q(t) passed from A to B for an integer observer O q (an ideal ranger using a q-clock in a plain world) would not be of 20 seconds: Similarly, the spatial distance of the buoy from the beach was previously measured to be ∆x = |x B − x A | = 100 m, while the geometric distance would be The speed of the sea turtle derived from the t-clock and the x-rod is where in the last step we assumed that T , X ≪ 1 (∆t ≪ t * , ∆x ≪ ℓ * ). Thus, the x-speed of the animal is 5 m s −1 but the speed v q in the integer picture is different. Clearly, if all experiments took place at the same spacetime and energy scales, the difference between the fractional and the integer picture would only be in the convention of the observer's measurement units. Living in a fractional world where a sea turtle takes 20 s to reach a buoy 100 m away would not be physically different from an integer world where the same event takes place in, say, 22 or 18 s. Integer turtles would be slower or faster (albeit not tremendously so) than fractional turtles, but that would just be the normality for integer observers.
Suppose now that, at dawn, the same ranger observes a hatchling (of size L ′ ∼ 5 cm = O(10 −1 ) L much greater than the characteristic length ℓ * ) getting out from a nest in the sand, reaching the sea line unhindered by waves or curious tourists, entering the surf at the same point A ′ = A of the adult and reaching the same buoy B with a speed v x ′ = ∆x ′ /∆t ′ . Again, an integer observer O q would disagree with the measurements of the fractional observer O x , but by a different relative amount because T = T ′ and X = X ′ . Then, using Eq. (27), one can determine the ratio of the speed of the adult over the speed of the hatchling in both frames: where Before interpreting these results, we can indulge a bit more in our thought experiment and see what would happen if the fundamental scales were as large as ℓ * 1 m and t * 10 s. Then the corrections T 1 and X 1 would be comparable with, or even dominate over, the standard term in the measure. It is convenient to separate between two cases, a time-like multi-fractal where X = 0 and a space-like multi-fractal where T = 0. In the time-like case, one has in general, except when −2 ≤ T ≤ 0. In the space-like case otherwise, except when −2 ≤ X ≤ 0. Multi-fractals along both spatial and time directions can display a more complex behavior.

Initial-point presentation and dimensionless observables
In the initial-point presentation, ∆q . According to Eqs. (22), (25), and (26), the geometric time interval ∆q + (t) and the spatial distance ∆q + (x) between A to B for an integer observer O q would be longer than for O x because α 0 , α > 0 and T + , X + > 0: Depending on which between the spatial and the temporal correction dominates, we will have v q < v x (timelike multi-fractal) or v q > v x (space-like multi-fractal): just as in the general case (30), (31). As x A ′ ,B ′ = x A,B , one has X = X ′ and the spatial correction in Eq. (28) cancels exactly. For the sake of the argument, let ∆t ′ = 10∆t. Then, since T ′ + = T + /10 1−α0 < T + , the temporal correction in Eq. (28) is positive and one has r q < r x = 10.
Here is the crucial point at last. The ranger is acquainted with the fact that the ratio r x = (speed) adult /(speed) hatchling of the swimming speed in the fractional picture depends on details such as the physiology of these animals, the temperature and time at the moment of hatching, the temperature of water, and so on, but that the average ratio reaches the universal, empirical value r x ≃ 10: in general, adults are about ten times faster than hatchlings. 7 Repeating their measurements every night and morning of the hatching season and finding a distribution of results definitely peaked at values close to r x = 10, the ranger agrees qualitatively with the turtle-speed law both as a fractional and as an integer observer, r x ≃ r q . However, imagine now that the relative experimental uncertainty reached by the ranger is better than in previous experiments, to the point where it is smaller than the correction T ′ − T in (28). Then they also find a systematic discrepancy of the data points in the integer picture and a deviation of r q to values smaller than 10. Then the observer must conclude that they are living in a fractional world.
Thus, even if geometric coordinates q µ (x µ ) have the same units as coordinates x µ , the existence of measurable dimensionless quantities such as the ratio r in our fictitious example allows an observer to recognize whether the underlying geometry is standard or multiscale. This is an aspect of the relation between measurements and frame choice not covered in previous discussions [36,39,41] and constitutes one of the main messages of the paper.
In practice, the outcome of experiments is less exciting, especially in macroscopic physics such as that describing turtle sea-faring. The magnitude of the corrections X and T are unknown since ℓ * and t * are free parameters of the theory. Therefore, the error bar can at most place an upper bound on such correction and the experiment cannot distinguish between a fractional and an integer world. Unhappy with this situation, the wildlife ranger may decide to change clothes and go to some laboratory or particle accelerator to perform altogether different experiments, this time involving atoms and quantum particles. For instance, they might want to check out the spectral lines of light atoms or the relativistic quantum particles generated in scattering events at high energies. This is precisely the type of observations considered in [40,65], where upper bounds on t * and ℓ * have been derived recently. These bounds are about 20-30 orders of magnitude smaller than the scales involved in the sea-turtle observation (t * < 10 −27 s, ℓ * < 10 −19 m) and there would be no way to discriminate a turtle on such a multi-fractal from one on ordinary earth.

Final-point presentation
Since T − = −T + < 0 and X − = −X + < 0, the situation for time-like and space-like fractals is reversed in the final-point presentation. Intervals (25) and (26) are shorter with respect to the fractional picture, while the velocities obey Also, r x < r q in the turtle example and there is a systematic excess in the geometric ratio r q with respect to the observed value.

Space-and time-like fractals
We now pause for a moment and discuss a subtle point connecting some independent data we collected so far. One is that particles diffuse at a slower speed in a "traditional" space-like fractal, d s space < d h space . Another is that non-relativistic velocities are slower on a space-like multi-fractional geometry in the initial-point presentation and faster in the final-point presentation. Then, if we reasonably assume that the local non-relativistic velocities of a coarse-grained random motion 8 follow the same trend discovered for a macroscopic body, we can infer that spacetimes with q-derivatives are space-like fractals when the measure is in the initial-point presentation [Eq. (34)]. Moreover, if T = 0 then one could represent a genuine diffusion process on space with Eq. (14) where q β (σ) = t and ∇ 2 q is only spatial. However, this would not match with the discussion below Eq. (14) about the possibility of having an anomalous clock σ when diffusing on traditional fractals. This observation does not place any strong constraint on the theory because we have seen that spacetime fractals do not have the same properties as traditional ones, but it leads us to consider as fractal also spacetimes with a timelike multi-scale geometry in the final-point configuration [Eq. (36)]. However, we should also take into account that the general behavior (30), (31) is reproduced without pathologies only in the initial-point presentation, while in the final-point presentation it is possible to hit the intervals and to reach the opposite regime (36), (37).
In reality, the three pairs of cases (30)-(31), (33)- (34), and (36)- (37) are an idealization of more complicated configurations. The main and most obvious reason is that space and time are entangled and, if the time direction is multi-scale, then the corrections T and X in Eq. (27) compete with opposite signs and they can produce velocities v x = v q |1 − |T − | + |X − || < v q even in the final-point presentation, provided |T − | > |X − |. The sign of the overall multi-scale correction in Eq. (27) can be generic also in the null presentation (21), depending on the details of the problem.

Relativistic motion
In the theory with q-derivatives without gravity, the dynamics is invariant under the q-Poincaré transformations (9), which are linear in the geometric coordinates but non-linear in the fractional coordinates. The components of the Lorentz matrices Λ µ ν are standard and so is all the apparatus of special relativity when written in geometric coordinates. In particular, the Lorentz transformations of time and space are where the Lorentz factor γ is [38] γ = 1 and c is the geometric speed of light (i.e., the speed of light in the integer picture), which is constant. The line element and dynamics of a relativistic particle are discussed in [38], while the speed of light c x in the fractional picture is computed in a parallel work [69]. From the dispersion relation of photons in electrodynamics, it can be shown that the magnitude of the difference ∆c = c x − c depends non-trivially on the energy of the photon and that ∆c = −O(1) (E/E * ) 1−α < 0 for a space-like fractal spacetime (the actual speed of light is smaller than c) and ∆c = O(1) (E/E * ) 1−α0 > 0 for a time-like fractal spacetime (c x > c) [69]. Also, for a specific presentation of the measure in momentum space ∆c ≡ 0 provided α 0 = α. In any case, ∆c is experimentally constrained to be so small that it can be neglected here, c x ≈ c. The Lorentz factor (40) in the integer picture yields a classic result of special relativity: the Galilean velocity v q of a body cannot exceed the geometric speed of light. However, from Eq. (27) we have seen that v q = v x /f , where f = |(1+T )/(1+X )| is the multi-scale correction computed above. Then the upper bound on v x is not c but For a space-like fractal spacetime, f = |1 + X | −1 < 1 in the initial-point presentation for any X > 0 or in the final-point presentation for X < −2. At microscopic scales, the speed limit is suppressed by a tiny factor f . This is consistent with both the sub-diffusion effect and the non-relativistic result. For a time-like fractal spacetime, f = |1 + T | > 1 in the initial-point presentation for any T > 0 or in the final-point presentation for T < −2. The Galilean speed v x in the fractional picture can exceed the geometric speed of light by a factor f . For a macroscopic object this factor is mild and very close to 1, while for an object at scales L ∼ ℓ * and t ∼ t * one can break the c-limit in a more spectacular way. Since ℓ * and t * are at least as small as particle-physics scales, these results imply that one cannot use this multi-fractal theory to construct useful faster-than-light spacecraft.

Future developments
Our main results have already been outlined in Sect. 1 and we will not repeat them here. We rather comment upon two open subjects to be tackled in the future: the effect of logarithmic oscillations and the role of the presentation in the microscopic structure of multi-scale spacetimes.

Log oscillations
When log oscillations are taken into account, multifractional corrections become modulated and can change sign within the same time-like or space-like configuration. To maximize these effects, we assume that they all come from the oscillatory part of the leading term in the measure and that fractional power-law corrections are negligible. Assuming for simplicity one frequency and a vanishing amplitude for the sine contribution, for each spatial direction we have plus a similar expression in time. Here A > 0 and ℓ ∞ is a fundamental scale at the bottom of the whole hierarchy ℓ ∞ < ℓ 1 < ℓ 2 < · · · of the most general measure. Given an interval ∆x = |x B − x A | and assuming x B > x A , the correction to the geometric distance can be approximated as Naively, the easiest way to simplify this expression is to consider distances small enough to be sensitive to the fractal properties of the background. In this case, ∆x = ǫℓ ∞ , where ǫ ≪ 1. Then F ω (x B ) ≃ F ω (x A ) and the last term in Eq. (43) cancels out. However, since we are on a fractal we can liberally apply the discrete scale invariance x → λ n ω x of the measure and note that We can satisfy this condition by preparing the experiment. If x A and x B lie in the same copy of spacetime, then (44) is automatic. If, by analogy with Cantor dusts, "to lie in the same copy" means to be in the same connected component at any given iteration (i.e., to belong to either S 1 (F ) or S 2 (F )), then by definition all causal experiments take place on the same copy and (44) is guaranteed. If this tentative interpretation were not correct, one could fine tune the initial and final point by hand to obtain (44). Such a fine tuning is not strong; for instance, when α = 1/2 one has λ ω = 1/N 2 and the minimum tuning is of 1/4. In both cases, one gets (approximately or exactly) Regardless of the magnitude of the amplitude A, the correction F ω can take either sign. 9 In general, due to the modulation of the oscillations we would not be able to connect relations of the sort v x < v q for non-relativistic velocities with the local velocities in a coarse-grained sub-diffusive stochastic process, even in a purely space-like or time-like fractal. The discussion of Sect. 4.4 would then need a revision. All these features will deserve to be explored in greater detail.
6.2 Nowhere differentiability: toward stochastic spacetimes Implicitly, in this paper we have begun to collect some evidence that there is a connection between the presentation choice and the stochastic properties of diffusion in these spacetimes. This interesting point went unnoticed in extant studies [37,42] and it is worth looking into it in detail.
We have seen in Sect. 4.4 that, depending on the presentation of the measure, in certain time-or spacelike systems the velocity of a non-relativistic body is slower than in an ordinary spacetime. Next, we have argued that this property is plausibly compatible with the microscopic sub-diffusion on fractals. The relation v x < v q is valid outside the "box" (38) (i.e., the range of values for T and X ), which is not a region in the parameter space of the theory: both T and X depend on the measurements taken in the given system. If we require property C, then the only system-independent configurations are (33), (34) and (36), (37). In that case, we set either the time or the spatial directions to be ordinary (this is part of the definition of the model) and the sign of the correction is unique for any system under consideration and for any regime. The box (38) should then be abandoned as a robust criterion for subdiffusion.
However, we saw in Sect. 3.2 that property C holds independently of the presentation of the measure. The presentation, in fact, does not affect the scaling of q, nor any of the scaling relations which define the dimensions d h , d s , and d w . The only element of ambiguity in these relations, fixed by an educated guess but never dispelled completely so far [37], was on the scaling q β ∼ σ β of the diffusion operator ∂/∂q β (σ) (and its multi-scale generalization) but not on its presentation. The universality of the spacetime dimension with respect to presentation choices is almost in accordance with the results (30), (31), valid everywhere in the parameter space except in the finite box (38). Therefore, at least in the theory with q-derivatives, we must agree that: -Spacetime space-like fractals have similar properties than ordinary spatial fractals, for instance the microscopic origin of sub-diffusion (i.e., the v x < v q relation), but not in all presentations of the measure. Having already analyzed the first three points, we can be more quantitative about the last. The choice of presentation of the integration measure discussed in Sect. 3.3 is intriguingly similar to the dilemma between the Itô [70] and Stratonovich [71] interpretation of stochastic integrals (see [72,73] for extensive accounts). Given a stochastic process X(t) (i.e., a sequence of random variables) and a Wiener process B(t) (also known as standard Brownian motion), the solution of the stochastic differential equationẊ(t) = a[X(t), t] + f [X(t), t]Γ (t) with initial condition X(t i ) is where dB = dt Γ . The first term is the initial point of the process. The second term drives the deterministic evolution of X and poses no conceptual problem, contrary to the third term. The latter introduces a stochastic noise in such evolution and a source of ambiguity not apparent in the formal expression (46). In fact, B(t) is nowhere differentiable, 10 meaning that ordinary integer-order derivativesḂ,B, . . . are ill defined. Splitting the time interval ∆t = t − t i as t i = t 0 < t 1 < · · · < t n−1 = t, without the pretense of being rigorous we can approximate the stochastic component in Eq. For the Riemann-Stieltjes integral of an ordinary differentiable function B(t) = b(t), the choice of the point t j ∈ [t j , t j+1 ] inside the interval ∆t j = t j+1 − t j is immaterial for the evaluation of the sum in the right-hand side: the result is unique and the Riemann sum is well defined. However, the Wiener process B(t) fluctuates so much in ∆t j that different choices oft j lead to inequivalent outcomes.
In the Itô interpretation,t j = t j is the initial point in ∆t j and the function f only depends on the behavior of B(t) up to the time t j . In this case, the stochastic process X(t) is a martingale: the expectation value of an event at some future timet j ∈ (t j , t j+1 ] is equal to the value observed at the present time t j . In other words, the knowledge of all previously observed values does not help to predict future outcomes. 11 In the Stratonovich interpretation,t j = (t j+1 + t j )/2 is taken in the middle of the interval and one symmetrizes between past and future (see the end of Sect. 4.3.6 of [73] for caveats). In this case, X(t) is not a martingale: knowledge of prior outcomes may help to determine future events.
Both interpretations are valid but not in absolute terms: they describe systems with different stochastic properties (e.g., [74]). Similarly, both the initialand the final-point presentations of the multi-fractional measure correspond to the same class of spacetimes but with different prescriptions on the volume of unit balls [34]. Multi-fractional integrals are of a form similar to Eq. (47). Considering a class of measures q(x−x) all with the same scaling, we have where x 0 = x i −x. On a discontinuous genuine fractal, we would have a nowhere differentiable measure q(x) (see, e.g., [75] and the discussion in [34]) and the same interpretation dilemma as in Eq. (47). In the theory with q-derivatives, on the other hand, q(x) is differentiable and there is no ambiguity inx j ∈ [x j , x j+1 ]. What we have, instead, is a presentation ambiguity which reflects in the choice of boundary conditions, i.e., the integration interval [x 0 , x n−1 ]. In Table 1 we list the choices for the four different presentations discussed in the text.

Presentationx
x 0 x n−1 Null 0 x i x Initial-point  The change in the integration domain is more trivial than the ambiguity in the stochastic term (47). We can easily understand why. When we consider, for instance, volcano-like measure weights ∂ x q(x) ∼ |x −x| α−1 , we find an integrable singularity at x =x. However, in a fractal we would expect to find a volcano at each and every point in the set, not just at a specific locationx. On a curved background, this issue is solved by general relativity: a local inertial frame centered on the observer is locally isomorphic to multi-scale Minkowski spacetime and each and every local inertial frame has its own volcano [39]. Still, it may be desirable to have also a global notion of irregularity, for instance when we work out particle-physics models on a flat Minkowski background. Such a global notion is not available in the theory with q-derivatives. A naive attempt to integrate over all possiblex leads to ill-defined expressions of the form )]L, unless the ambient spacetime is compact. But even in a a compact space, this modification of the measure would not correspond to "volcanoes everywhere." This is clear in the discretized version of the integrals: The last line is exactly what we would need: picking different presentations would amount to a non-trivial selection of the argumentx j +x j . We should bear in mind that the theory with qderivatives is a simplified version of a genuine nowhere differentiable fractal. Barring explicit and most challenging constructions of field theories on fractals [76][77][78] or field theories constructed with stochastic measures, the closest thing mimicking nowhere differentiability is fractional calculus. The multi-scale theory with fractional derivatives [34,35] incarnates precisely this possibility and it may be the only multi-fractional theory obeying property D in the introduction. In fact, the differential operators in the theory with ordinary derivatives are the usual partial derivatives ∂ x . In the theory with weighted derivatives, they are weighted version of the same operators, ∂ x → (∂ x q) −1/2 ∂ x [(∂ x q) 1/2 · ]. In the theory with q-derivatives, they are again of integer order but with a different weight distribution, ∂ x → (∂ x q) −1 ∂ x . These modifications of ∂ x and the non-trivial integration measure give all these models an "irregular" geometry but in a rather simple-minded way. On the other hand, the derivatives of the fourth multi-fractional theory are non-local integro-differential operators, which are known to capture the properties of sets not differentiable in the ordinary sense (see [34,75] and references therein).
Therefore, the theory with fractional derivatives may well be the only one to describe a multi-fractal geometry in the strong sense. Due to its higher technical challenges, this framework has not been explored as extensively as the other multi-fractional theories but the preliminary analysis in [34,35] and work in progress show enticing properties that include an exotic particle content and an improved perturbative renormalizability (absent in the other cases [79]). The arguments presented in this section add fuel to our curiosity and strongly suggest that spacetimes with a fractional integro-differential structure would be intrinsically stochastic. We hope to report on that soon.