ORIGINAL RESEARCH article
Towards Multifractional Calculus
- Instituto de Estructura de la Materia, Consejo Superior de Investigaciones Científicas, Madrid, Spain
After motivating the need of a multiscale version of fractional calculus in quantum gravity, we review current proposals and the program to be carried out in order to reach a viable definition of scale-dependent fractional operators. We present different types of multifractional Laplacians and comment on their known or expected properties.
A branch of theoretical physics which has been attracting considerable attention in the last years is quantum gravity. Several independent theories, models and hypotheses are gathered under this broad name, from string theory to asymptotic safety, from non-local to loop quantum gravity, from causal dynamical triangulations to causal sets, and so on [1–4]. Most of these proposals aim to conciliate classical general relativity with the laws of quantum mechanics, in order to unify all forces of Nature under the same framework and to solve some problems left open in the traditional paradigms (for instance, the big-bang and cosmological constant problems ).
A surprising feature emerging from this variegated landscape is that the properties of spacetime geometry, such as the spectral or Hausdorff dimension and the way particles diffuse, change with the probed scale in all quantum gravities [5–7]. This so-called dimensional flow seems to be a manifestation of the impossibility to perform infinitely precise time and distance measurements in geometries with intrinsic uncertainties of quantum or stochastic origin [8, 9]. Some of these findings were made possible by assuming dimensional flow by default and treating spacetime geometry as fundamentally scale dependent. This general method can be embodied in a class of theories, called multifractional, where classical and quantum fields live on a spacetime characterized by a scale hierarchy, anomalous transport and correlation properties, and a multifractal geometry . Surprisingly, all these features emerge automatically by assuming a slow dimensional flow at large scales (dimension in the infrared almost constant) [10, 11].
One can encode a multiscale geometry in the dynamics of particles and fields in several ways. The one followed by multifractional theories is a change in the integro-differential structure . Integrals (such as dynamical actions) and derivatives (in kinetic terms) acquire a non-trivial scale dependence that can be illustrated in the prototype example of the scalar field theory
where ϱ(x) is the spacetime measure, is a kinetic operator, and V is the scalar potential. In the standard case and in the absence of gravity (which will be ignored here), ϱ(x) = dDx is the usual Lebesgue measure in D topological dimensions and is the second-order Laplace–Beltrami operator. In the presence of dimensional flow, if the measure is factorizable in the coordinates (an assumption to make the problem tractable) then it takes the unique form 
all indices μ are inert (there is no Einstein summation convention), the first factor 1 in Equation (1) is optional  (it can be set to zero in the stochastic version of the theory ), and are 2D length scales for each n, and αμ, n, Aμ, n, Bμ, n, and ωμ, n are 4D real constants for each n.
Since the measure is factorized, in the following we can focus the discussion on the one-dimensional model
where ℓ* = ℓ1, all μ indices are omitted and we can also ignore terms subleading both in the infrared and in the ultraviolet; this corresponds to consider only the n = 1 term in Equations (2) and (3).
In this paper, we will study the properties of three versions of the kinetic operator , expanding on the proposals sketched in Calcagni . Since, in the context of quantum gravity, the integration measure is uniquely defined independently of the type of derivatives in the Lagrangian as in Equation (4) , here we are not interested in the formal properties of “multiscale integrals,” the inverse of multiscale derivatives. Some of these operators are known, as is the case of Equations (5) and (35) below , while in the case of Equations (28) and (29) they are unknown and will require further work. On the other hand, there is no inverse operator for a linear combination of operators with different inverse, such as Equation (24). In all these cases, for our purposes it is sufficient to study the properties of multiscale derivatives with respect to the ordinary Lebesgue measure dx while, at the same time, taking into account the measure weight by inserting weight factors in the definitions of such derivatives to make them self-adjoint with respect to the measure.
While there is a unique parametric form of the measure q(x), there is more freedom in the choice of kinetic operator . It turns out that there are three viable possibilities. One is a theory with so-called weighted derivatives, but this can be reduced to a system with ordinary derivatives and the spectral dimension of spacetime is constant in that case . Another possibility is the second-order operator 
where and q(x) is given by Equation (2). This “q-derivative” has a number of highly desirable properties:
1. It is multiscale, since the scale hierarchy is already encoded in the measure weight v(x).
2. Its composition law is very simple:
3. It is linear. For any f and g in a suitably defined functional space,
4. Its kernel is trivial and given by a constant:
5. The Leibniz rule is extremely simple. For any f and g in a suitably defined functional space,
6. Integration by parts is straightforward. For any f and g in a suitably defined functional space,
where we threw away boundary terms. Consequently, is self-adjoint:
Notice that, in principle, these rules hold for an arbitrary q(x), although in our case this profile is fixed as in Equation (4).
3. Fractional Derivatives
The third extant multifractional theory is the least explored, but also the most interesting because it employs fractional calculus. This is by far the most obvious tool to implement an anomalous scaling in the geometry. The application of fractional derivatives to multiscale theories is not an easy task. Before seeing why, let us recall some basic aspects of fractional calculus.
There are different versions of fractional derivatives [14–16]1 and one must make a choice suitable for quantum gravity . In particular, we believe that one cannot renounce to have a trivial kernel (Equation 8). Two fractional derivatives with this property are the Liouville derivative
and the Weyl derivative
where θ is Heaviside's step function. Obviously, these operators act linearly on f and . One can also check that and (these fractional derivatives commute) and that the Leibniz rule is
and the same expression for the Weyl derivative, where ∂α−j = Ij−α are integrations for j ⩾ 1. Also, integration by parts with the Liouville derivative generates the Weyl derivative, and vice versa:
3.1. Complicated Leibniz Rule
The importance to have the standard Leibniz rule (Equation 9) can be appreciated when trying to do physics with fractional calculus. In the theory with q-derivatives, integration by parts does not produce extra contributions and the kinetic terms or −∫dq∂qϕ∂qϕ are completely equivalent. Therefore, the equation of motion can be determined easily by applying the variational principle on Equation (4). On the other hand, suppose we choose another type of multiscale derivative such that and its Leibniz rule is more complicated:
where X = X(f, g; x) is a function of f, g, their ordinary derivatives and the coordinate x. For consistency, if the kernel of is trivial (), then X(f, 1;x) = X(1, f; x) = 0 for any f. In particular, if , then
Therefore, when varying the action (Equation 4) with respect to ϕ one gets
Assuming that one could repeatedly integrate Y by parts to write it as Y = 2δϕZ(ϕ, x) up to some boundary term, we would end up with a dynamical equation
characterized by a term Z that can considerably hinder the study of solutions.
This is the main obstacle that prevented so far to consider multiscale theories with derivatives different from Equation (5) (barring the mathematically trivializable case of weighted derivatives). In fact, the only derivative with anomalous scaling such that X = 0 in the Leibniz rule (Equation 16) is the q-derivative . Genuine fractional derivatives always have X ≠ 0.
3.2. Self-Adjoint Laplacian
Although X≠0, one could still obtain a clean integration by parts if, thanks to miraculous cancellations, Y were a total derivative or Z were zero on shell. This possibility is suggested by Equation (15), which implies that, for any combination
For instance, if 2,
In the limit α → 1, and , so that . Therefore, we can define an operator self-adjoint with respect to any measure weight v(x):
Note that other, complex-valued choices of may be more convenient when studying the spectrum of eigenvalues of these operators .
4. Multifractional Derivatives: Three Proposals
At this point, we can try to extend fractional calculus to a multiscale setting. We have found three ways to do that.
4.1. Explicit Multiscaling
The most direct mean to induce a hierarchy of scales and a variable anomalous scaling is to consider a superposition of fractional derivatives of different order α . In the mathematical literature, several authors [25–32] did propose a continuous superposition, the distributed-order fractional derivatives , where m(α) is a distribution on the interval [0, 1]. However, from previous experience in quantum gravity it may be more convenient, or just sufficient, to take a sum instead of an integral:
where gn = gn(ℓn) are some constant coefficients and is defined in Equations (22) and (21). A non-trivial dimensional flow is generated by just one scale, i.e., a sum of two terms: . The equation of motion from the action (Equation 4) with kinetic operator (Equation 24) is
This formulation of a multiscale theory with fractional derivatives is not exempt from problems. The operator consists of many terms, even in the simplest case of only one scale where is made of seven pieces (ignoring weight factors), . Therefore, the dynamics (Equation 25) is deceptively clean and hides a rather messy multiorder fractional differential structure which may be very difficult to solve analytically. This eminently practical issue could be very important, or even fatal, at the time of studying the dynamics. To bypass it, one could consider another version of the kinetic operator :
where m = 2 in Equations (12) and (13):
At the classical level, the great advantage of Equation (26) is that, in the single-scale case, it consists of just three terms (again, weight factors are ignored) instead of seven. However, this is not quadratic, since , which can lead to problems when quantizing the theory in Hamiltonian formalism: the kinetic term is not the square of a momentum operator.
At present, it is not clear which definition between Equations (24) and (26) will be more viable in the long run. They differ only in transient terms that can be dropped both at large and small scales, so that classically they give rise to the same physics. However, both have the added inconvenience of leading to a virtually symmetryless dynamics , a further point of concern if we want to do field theory and gravity with this formalism.
4.2. Implicit Multiscaling
The multiscaling characterizing Equation (24) is of a twofold nature, an explicit one in the sum over αn and an implicit one in the measure weight v(x). These two structures have been combined independently and we imposed that the sum over αn in the combination of fractional derivatives is the same sum over αn inside v(x). There is nothing wrong with this construction, but there may be a more elegant formulation where the scale hierarchy is all included within the measure q(x) . Noting that the denominator (x − x′)α in Equations (12) and (13) for m = 1 (0 < α < 1) is the ultraviolet part of the profile in Equation (4), we can generalize those definitions as a left and right multifractional q-derivative:
where, again, the profile q(x) is uniquely given by Equation (2). These expressions are similar to the so-called variable-order fractional derivatives proposed by Lorenzo and Hartley , although in our case q(x − x′) is completely fixed.
The kinetic operator in Equation (4) is then
To understand the dynamics, we first need to spell out the properties of these derivatives. At short scales, q(x−x′)~|x−x′|α and Equations (28) and (29) reduce to the Liouville and Weyl derivatives, respectively:
while at large scales q(x−x′)~x−x′ and Equations (28) and (29) give
We have not made a formal proof of these statements, but it should not be difficult. The Leibniz and integration-by-parts rules are also unknown but they should coincide with those of Weyl and Liouville fractional derivatives in the limit of small scales or in any plateau region of dimensional flow ().
Therefore, we expect a complicated Leibniz and integration-by-parts rules everywhere at all scales of dimensional flow, except in plateau regions where a clean integration by parts of the type (Equation 15) emerges. For this reason, a variational principle valid at all scales may be ill defined in this case and an exact form of the equations of motion may be out of reach, although their asymptotic form at plateaux is obviously given by the limit of Equation (25) for one exponent α.
These considerations could eventually select the multifractional derivatives with explicit multiscaling as a simpler tool in quantum gravity, since they yield exact equations of motion.
4.3. Multiscale Differentials
A third alternative is to introduce a multiscale differential based on the geometric coordinate (Equation 2) or its simplified version in Equation (4) :
which is a linear combination of the usual and fractional  differentials, 𝕕q ~ 𝕕x + 𝕕|x|α + … = dx + |dx|α + … . [In D dimensions, this differential generates the line element , where gμν is the metric.] The operator 𝔻 is a superposition of ordinary and fractional derivatives of the form (to be taken as indicative; coefficients are ignored)
The following multiscale derivative and Laplacian are then defined implicitly:
These operators are invariant under translations, since , while is invariant also by “q-boosts” . Therefore, this theory has more symmetries than the theories with multifractional derivatives with explicit or implicit multiscaling.
In any plateau of dimensional flow, and . Notice that 𝔻 ≃ ∂q in the near-infrared limit 𝕕 → d where the non-linear part of q is subdominant. Therefore, the theory with q-derivatives can be regarded as an approximation of the theory with multiscale derivatives (and, presumably, also of the other two theories with fractional derivatives) when the anomalous scaling effects are weak. The experimental constraints on the scales and parameters of the theory with q-derivatives might thus miss some important effects present in the fractional versions of the multiscale paradigm.
To determine the Leibniz and integration-by-parts rules, one should first define the operator 𝔻 appearing in the differential 𝕕 = 𝕕q𝔻. Again, we expect these rules to reduce to the usual ones in the infrared and to those of fractional derivatives in the ultraviolet. Since the operator (Equation 24) with explicit multiscaling is already a well-defined linear combination of fractional derivatives, we reach the same conclusion of the previous section, namely, that the operator (Equation 24) may be the best candidate for the concrete realization of multifractional theories with fractional derivatives. However, the main problem of the definitions (33) and (35) is that they are too abstract, which is the reason why we used qualitative expressions marked by “~.” Understanding their actual properties will require more work.
In this paper, we have further analyzed the proposals of Calcagni  for a multifractional calculus with viable applications to field theory and gravity. Without the pretense of being rigorous, we have considered some properties of scale-dependent derivative operators which, in physical applications to quantum gravity, are interpreted to encode the multiscaling of the underlying anomalous geometry. The conclusion is that the most promising multifractional theory possibly is the one with explicit multiscaling in fractional derivatives. However, only a full systematic study of the properties of all these operators will be able to confirm which is the most viable Laplacian from a theoretical and practical point of view. The value of the complex coefficients in Equation (20) will be especially important to determine a well-defined calculus and spectral theory . We will analyze the associated dynamics in detail in a future publication.
The author confirms being the sole contributor of this work and approved it for publication.
Conflict of Interest Statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author is under a Ramón y Cajal contract and is supported by the I+D grants FIS2014-54800-C2-2-P and FIS2017-86497-C2-2-P.
1. ^For bibliographic references, see [12, 10]. More recent applications and solving methods can be found in Baleanu et al. , Yang et al. , Sun et al. , Baleanu et al. [20, 21] and Inc et al. .
18. Yang XJ, Tenreiro Machado JA, Baleanu D. Anomalous diffusion models with general fractional derivatives within the kernels of the extended Mittag-Leffler type functions. Rom Rep Phys. (2017) 69:115.
20. Baleanu D, Inc M, Yusuf A, Aliyu AI. Lie symmetry analysis, exact solutions and conservation laws for the time fractional modified Zakharov–Kuznetsov equation. Nonlin Anal Mod Control (2017) 22:861. doi: 10.15388/NA.2017.6.9
21. Baleanu D, Inc M, Yusuf A, Aliyu AI. Lie symmetry analysis, exact solutions and conservation laws for the time fractional Caudrey–Dodd–Gibbon–Sawada–Kotera equation. Commun Nonlin Sci Num Sim. (2018) 59:222. doi: 10.1016/j.cnsns.2017.11.015
22. Inc M, Yusuf A, Aliyu AI, Baleanu D. Dark and singular optical solitons for the conformable space-time nonlinear Schrödinger equation with Kerr and power law nonlinearity. Optik (2018) 162:65. doi: 10.1016/j.ijleo.2018.02.085
29. Chechkin AV, Gorenflo R, Sokolov IM. Retarding subdiffusion and accelerating superdiffusion governed by distributed-order fractional diffusion equations. Phys Rev E (2002) 66:046129. doi: 10.1103/PhysRevE.66.046129
Keywords: quantum gravity, fractional calculus, fractional derivatives, multiscale geometry, multifractional spacetimes
Citation: Calcagni G (2018) Towards Multifractional Calculus. Front. Phys. 6:58. doi: 10.3389/fphy.2018.00058
Received: 09 February 2018; Accepted: 23 May 2018;
Published: 21 June 2018.
Edited by:Dumitru Baleanu, University of Craiova, Romania
Reviewed by:Aliyu Isa Aliyu, Federal University, Dutse, Nigeria
Yasser Dehnavi, Islamic Azad University of Bojnourd, Iran
Copyright © 2018 Calcagni. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gianluca Calcagni, email@example.com