Dualities and Asymptotic Mixtures Using Functional-Order Differentiation

: New deﬁnitions for fractional integro-differential operators are presented and referred to as delayed fractional operators. It is shown that delayed fractional derivatives give rise to the notion of functional order differentiation. Functional differentiation can be used to establish dualities and asymptotic mixtures between unrelated theories, something that conventional fractional or integer operators cannot do. In this paper, dualities and asymptotic mixtures are established between arbitrary functions, probability densities, the Gibbs–Shannon entropy and Hellinger distance, as well as higher-dimensional particle geometries in quantum mechanics.


Introduction
One of the goals of mathematics, science and engineering is to try and find dualities between two or more different theories or approaches that provide the same solutions even though they appear to be independent mathematical entities. In this paper, it will be shown that a modification of the differential and integral operators of fractional calculus can be used to establish dualities and asymptotic solutions between two (or more) unrelated mathematical functions or theories. What is referred to as functional-order differentiation is used to show that two probability densities, such as the one-parameter Exponential and the two-parameter Pareto, can form mathematical dualities, since either one can provide exact solutions of the other. Asymptotic solutions, which arise from dualities between these densities, are, in fact, mixtures of the two. Functional differentiation is then used to establish a Gibbs-Shannon entropy duality with the proper distance metric due to Hellinger. The duality means that both provide the same results for arbitrary probability densities, although one measures entropy and the other measures metric distance, which is symmetric and obeys the triangle inequality. Once again, asymptotic solutions between the two approaches are obtained, which show that an optimal minimum separation between densities is provided by a mixture model consisting of the Gibbs-Shannon entropy and Hellinger's distance equations. Finally, functional differentiation is used to study the relationship between hyper-geometries and the duality between functional order and fractal dimension as applied to quantum mechanics. The solution to such problems requires a modification of the operators that appear in fractional calculus, since the latter and the conventional integer operators cannot achieve these results.
Fractional calculus goes back as far as 1695, when l'Hôspital pondered about the meaning of a derivative to half order in a letter to Leibniz. The word 'fractional' is a historical misnomer, given that it refers to 'general' or 'generalised' calculus. Hence, fractional calculus generalises conventional integer calculus such that the latter is a special case of the former. Throughout the centuries, fractional calculus has been developed by mathematicians such as Lacroix, Fourier, Liouville, Riemann, Grunwald, Letnikov, Abel, Riesz, Weyl, Caputo, and many others. As a result, various formulations appear in the literature at present [1]. Fractional calculus has become an intense area of research, starting from the second half of the last century, as many physical phenomena were better described by this. The continuous development of fractional calculus has led to a wider range of applications, not only in mathematics but also in the study and analysis of physical and engineering problems . The parameter at the core of fractional differentiation or fractional integration is the fractional order α. The fractional order is a constant that, unlike its integer counterpart, can also take real or complex values, i.e., α ∈ R and α ∈ C.
The natural outcome of having a constant value for the fractional order is to consider extending it to cases where it is a function of one or more variables. While this is somewhat obvious, it has not gained much traction, even after it was officially proposed in 1993 by Samko and Ross [24]. The primary issue has been the difficulty in obtaining closed-form solutions when the fractional order is also a function of the integration variable appearing in the fractional derivative and integral formulations, respectively. In other words, if the fractional order is a function of the integration variable x, i.e., α = α(x), then the integration can become intractable, which is why choosing the correct variable dependence for the fractional order is very important. If it is a requirement that the fractional orders contain the same variable as integration, then it can be series-expanded and a sufficient number of terms can be retained for accuracy and to aid integration. In most cases, the fractional order is a function of other parameters that appear in the integrands of the fractional derivative and integral operators. The case where these parameters strictly represent space or time has been extensively studied in [25]. Due to the difficulty in obtaining closed-form solutions, another approach consists of using a modified version of the Grunwald-Letnikov series method with a parameter-dependent fractional order [26][27][28]. The other alternative is to numerically calculate such fractional derivatives and integrals [29,30].
Differentiation or integration operators that contain a fractional order that is a function of arbitrary variables has come to be known as variable-order (VO) fractional differentiation and integration in the literature. Variable-order fractional calculus is used to solve many physical problems, with impressive results, even better than those obtained via the standard case, where the fractional order is constant. Some of these applications include linear-to-nonlinear responses [31] and a transition from sub-diffusive to super-diffusive flows [32]. Other problems that are solved via the use of VO fractional calculus appear in mechanics [33], viscoelasticity [34], the modelling of anomalous transport processes [35], control theory [36], signal filters [37] and wave propagation [38]. A recent review of the latest developments in the area of variable-order fractional calculus can be found in [39] and the references therein. In addition, the reader can also refer to [40].
The ability to establish dualities and mixed solutions between independent functions or even physical theories cannot be performed via conventional integer-order or fractionalorder operators. In this paper, a slight variation in the variable-order fractional operators will be considered, which facilitates these mathematical dualities and/or mixed states. Instead of including the variable order α(t, x, ..) as part of the integration, the fractional operators will contain a delayed variable order. What this means is that the fractional operators will be used, as in the standard approach, where the fractional order is a constant and the latter will be mapped to a variable order α → α(t, x, ..). Not only does this avoid integration complexities but it will be shown that it can be used to solve a number of interesting problems, some of which will be presented in the paper. Given that the fractional operator is 'delayed' and not part of integration, as is the case with the variableorder definition, it will be referred to as functional differentiation (or integration), since the fractional order is a function of parameters that are not directly dependent on the initial integrations.
The paper is organised as follows. Section 2 derives the fractional derivatives of some functions that will be required later in the paper when discussing functional differentiation. In Section 3, functional differentiation is presented. This is then used to find dualities and asymptotic solutions between arbitrary functions, two probability densities, the Gibbs-Shannon entropy and Hellinger distance and hypergeometries, as applied to quantum mechanical particles. Conclusions are presented in Section 4.

Derivative of Functions to Fractional-Order
Conventional differentiation and integration operators have been generalised using their corresponding fractional forms. The order of either differentiation or integration is not discrete but continuous, so that the order α can be either real or complex, i.e., α ∈ R or α ∈ C. There are a number of representations in the literature for these generalised operators, with those due to Riemann-Liouville, Caputo and Grunwald-Letnikov being some of the most familiar. In what follows, the Riemann-Liouville definitions will be used, but other formulations can be used instead. The fractional derivative of functions due to Riemann and Liouville is defined by the operator a D α t , acting on some function f (t). The lower terminal a = −∞ is due to Liouville, and when it is a = 0, this is due to Riemann. Thus, the fractional derivative operator acting on a function takes the following form: where n only takes on integer values: n ∈ Z + . The fractional order α is bounded by n < α < n , where n is the floor function and n is the ceiling function. The fractional integration of functions can be obtained via the Riemann-Liouville derivative operator and is defined as follows: The fractional derivative (1) and fractional integral (2) can be re-written so that, instead of the fractional order only being a constant, it can also be a continuous function α = α(t, x): It is typical in the literature to refer to these two formulations (3) and (4) as the variable-order fractional derivative and fractional integral, respectively. The fractional order becomes a function of the variables t and x but, as can be surmised, performing the integrations with respect to these is very complicated, which is why the form for the fractional order is usually only a function of t, i.e., α(t, x) = α(t). What will be referred to as differentiation or integration to functional order is essentially the same as variational order differentiation and integration, but the functional order, which can be a function of arbitrary variables α(α 0 , x, y, t), will be a 'delayed' variational order. This means that normal fractional derivatives and integrals are obtained, as given by (1) and (2), and then the fractional order α is mapped to the functional order α(α 0 , x, y, t); hence, this is a 'delayed' order rather than immediately being part of the differentiation or integration, as given by the traditional variable-order expressions (3) and (4). As a result, there is less restriction, which allows for a considerable amount of flexibility when solving problems, as will be seen later in the paper.
Since the fractional derivative of monomials and exponential functions will be required in Section 3 as a first step in the process of differentiating functions to functional order; next, their fractional derivatives will be obtained, starting with the case of a monomial t β for some power β > −1.

Fractional Derivative of a Monomial
Using (1), let the integer order be n = 1 and the lower integration limit be a = 0. Then, the fractional derivative of a monomial f (t) = t β becomes Subtitute the transformation y = t − x in (5) to obtain the fractional derivative of f (t) = t β as follows: The fractional derivative of a monomial is obtained as (6) using the Riemann-Liouville method, which is equivalent to Lacroix's approach of generalising multiple integer derivatives of monomials to fractional form.

Fractional Derivative of Exponentials
Consider the fractional derivative of f (t) = exp(κt) for some κ > 0. Proceeding similarly to the monomial case in the previous section: Using transformations and simplifying gives the following result: From this, one can observe that From the fact that Γ(−α)/Γ(1 − α) = −1/α, we can see that the fractional derivative of exp(κt) becomes where (1) α = 1 ≡ e 2πi . Similarly, the fractional derivative of exp(−κt) is determined, which takes the form: Once again, via transformations and simplifications, the final result becomes: Equation (12) has a dominant term: This means that the fractional derivative has the final form, where Observe that when α = n ∈ Z + , corresponding to conventional integer differentiation, all terms in the above expressions collapse to the conventional integer derivatives. This is because, where the denominators contain the gamma-function, this becomes 1/Γ(1 − α) = 0 whenever α = n = 1, 2, 3, 4. . . Following the same mathematical process as above, it can be shown that the fractional derivative of f (t) = exp(±iκt) is given as: and Using combinations of (15) and (16), it can be shown that the fractional derivatives of f (t) = sin(t) and f (t) = cos(t) are given by and 0 D α t cos(t) = cos t + α π 2 (18) respectively. Thus, the fractional derivatives are generalised forms that contain the conventional integer derivatives as integer limits. At this point, it is worth addressing the approximations in (10), (14)- (16). The results contain terms that are very small and negligible, so that the exponential terms dominate in each case. The other terms are so small that they can be disregarded and only appear if the Riemann form for the operator is used, where the lower terminal in (7) is taken as zero. The correct approach is to use the Liouville operator, i.e., when the lower terminal is −∞. Then, the negligible terms cancel out to exactly zero and the exponentials are the terms that survive. In the strict sense, the latter approach has to be used, but the Riemann version was used to ensure consistency with all other presented cases. In the end, however, both approaches are correct.

Derivative of Functions to Functional-Order
The Riemann-Liouville formulation will be defined slightly differently to (1): where [α → α(x)] is the delayed-fractional order, which implies that, after the integral has been evaluated, the fractional order α is mapped to α = α(x). Notation-wise, the fractional order derivative operator (19) will be written as A similar scheme applies to the fractional integral operator (2), which can be written as d −α(x) /dx −α(x) . For the rest of the paper, only the derivative as defined by (20) will be examined, for reasons of brevity, but the integral version of the integration of functions follows a similar argument. The proposition is that the fractional order can take a number of forms, such as: In the conventional use of the fractional order, it is always taken as a constant α = α 0 , with values that are real or complex, α 0 ∈ R or α 0 ∈ C. The second definition is a function that is independent of the integration variable x. The third definition is a function of integration variable x; however, when using this definition in (21), it should be understood that this is not part of the integration. Only after the fractional derivative is obtained can α = α(x) be used, as indicated by (19) or (20), i.e., this is a delayed order. The first case where α = α 0 has always been dealt with in the literature. The second case α = α(y) has recently been discussed in [6]. The case of α = α(x) will be examined here. On this basis, the definition for differentiation to functional order is straightforward: the fractional derivative of a function is obtained and the fractional order is replaced with an arbitrary continuous function α = α(x). Hence, the following notation is used: (22) This will represent the differentiation of a function f (x) to functional order. This simple definition has some very powerful mathematical properties, which allow solutions to be obtained to a number of problems. In addition, differentiation to functional order allows for dualities to be found between distinct functions. In other words, two independent functions f (x) and g(x) can be made equivalent, in the sense that Aside from functional orders α(x) that satisfy the relation (23), it should be noted that (23) can also be used to obtain asymptotic solutions which can be thought of as 'transitional' functions or limits that are mixtures of the two functions f (x) and g(x). Put another way, if the variable in the functional order is strictly time, α(t), then this has memory properties, encompassing past present and future states. If the functional order contains spatial dimensions, e.g., α(x), then it has spatial correlation properties, a kind of spacial dimensional fading effect. In both cases, these are equivalent to or represented by the asymptotic solutions.
If the constant parameter α 0 is added to the mathematical form of the functional order α(x), it can be used as the mechanism that picks out these asymptotic solutions. When a 0 = 1 is set in the functional order, then f (x) ≡ g(x), and when α 0 = 0, then α(x) = 0, and the solutions provide the exact form of f (x). This requires solving the expression (23) for α(x), which then makes it possible, through variations to α 0 , to obtain f (x) and g(x) or any combination of them. Amongst other things, this establishes a duality between the two independent functions. In addition, the functional order α(x) is an arbitrary continuous function with validity throughout the x domain. If the entire domain is not required, then a more localised derivative order (approaching a constant or integer value) can be obtained by taking the series of α(x) to any desired number of terms. In the following sections, differentiation to functional order will be demonstrated and, via the use of (23), used to solve some interesting problems.

Applications of Functional Order Differentiation
In this section, how functional differentiation can form dualities and asymptotic mixtures between various independent functions can be shown. Relating such functions is only possible via functional differentiation (integration). This is not possible with conventional fractional or integral differentiation. Consider the case of a monomial f (x) = x β , as discussed in Section 2.1. Its derivative to functional order is given as follows: for some functional orders, α(x). Equation (24) can posses a number of interesting properties depending on the mathematical form of the functional order α(x). Suppose that the monomial is quadratic f (x) = x 2 with β = 2, and let a second function be g(x) = exp(x). It is clear from Figure 1 that f (x) and g(x) do not intersect for x ≥ 0. As they do not have any common solutions, they are independent functions. If the functional order is taken to be equal to the function g(x), for example, i.e., α(x) = g(x) = exp(x), then it is possible to relate the two independent functions so that they intersect. There are four solutions that are related by the derivative of f (x) to functional order, as can be seen in Figure 1. This is achieved via (24) so that the derivative to functional order becomes: Note that when x = 0, (25) reduces to the conventional integer derivative while the original function f (x) = x 2 corresponds to the case α(x) = 0 in (24). Hence, the function to the right of (25) now encompasses solutions that belong to the two independent functions f (x) = x 2 and g(x) = exp(x).
Next, allow two alternative functions f (x) = sin(x) and g(x) = sin(20x)e −x to be considered. The solution that makes both of these functions equal and provides their asymptotic 'transitions' will be derived. Firstly, from (23), The left-hand side of (26) has been calculated before; see (17). Hence, (26) becomes: Equation (27) holds if the functional order is given by: with α 0 = 1. The constant c ∈ Z, but here this can be set to zero: c = 0. Thus, the derivative of f (x) = sin(x) to functional order α(x), as given by (28) with α 0 = 1, provides the exact same result as function g(x) = sin(20x)e −x . This is shown in Figure 2, which also shows the function f (x) = sin(x), as well as asymptotic solutions corresponding to α 0 = 0.1, α 0 = 0.5 and α 0 = 0.9. When α 0 = 0, α(x) = 0, so original function f (x) = sin(x) is recovered. As the value of α 0 increases from zero to α 0 = 1, the asymptotic solutions first resemble the function f (x) = sin(x), and then they approach and resemble g(x) = sin(20x)e −x until the two are exact when α 0 = 1. Note that the asymptotic solutions approach the functions f (x) and g(x) from 'above', since the values for α 0 are all less than α 0 = 1. However, they can also approach these functions asymptotically from 'below' if their values were α 0 > 1. Another example of the use of differentiation to functional order is considered. From Section 2, use can be made of the functional derivative of an exponential function: Expanding the right-hand side of (29) to the first order and setting the result to the original exponential function gives the following: The functional order α(x) that allows for expression (30) to be exact is obtained as: Substituting (31) into (30) and transposing gives the result: Multiplying both sides by a and setting b = abx gives: Taking the n-th power of both sides and setting x = 1 gives, (a + b) n = a n e nb a b a n log(1+b/a)−nb/a log(b/a) Define the parameter z as and (34) can be written as: Equation (36) is an equivalent way of writing the binomial expansion, i.e., it is an identity or there is a duality between the binomial expansion and the expression on the right of (34). This duality may also allow for certain integrals to be performed more easily if, for example, integrals take the form: (a + b) n e −nb a dx ≡ a n(1−z) b nz dx (37) or some variation of (36). Here, a and b are functions or a combination of functions and constants, such that a = b. As an example, consider the relationship (37) and its use to solve: by setting a = x and b = 1 and where x = 0 and E i (·) is the exponential integral function. Thus, the complicated integral on the left is equivalent to a simple integral on the right, via the manipulation of (37), which, in turn, can be determined in closed form. As another example, it can be shown that the integral on the left: is equal to the closed form on the right, where erf(·) is the error-function, n ∈ Z + and a = x and b = 1. Similar expressions can be derived for other mathematical relations using the same approach as shown above. This highlights the power of using functional-order differentiation (integration) to form dualities between arbitrary and uniquely different functions. In what is to follow, functional order differentiation will be used to derive dualities and general relations in some physics and mathematics problems via the functional order.

Duality and Asymptotic Mixtures of Probability Densities
It is well-known that many physical, mathematical and engineering systems can be modelled via probability distributions. Using differentiation to functional order, a direct connection between different probability densities can be achieved when this is otherwise not possible. Consider the one-parameter Exponential-density with parameter λ: where the inverse of the parameter λ is equal to the mean. In addition, consider the two-parameter Pareto-density, which has been extensively used in many research areas to model various processes. The ability to model some very important phenomena using the Pareto distribution has prompted its generalisation to a fractional form. Among other things, the fractional Pareto distribution has been successfully used to model the statistical properties of radar clutter returns in [3], for example. The conventional Pareto density is given as follows: The parameter β regulates the shape of the distribution and x 0 is the support parameter. It is obvious that the two densities will have different distributions when compared to each other and cannot be equal, except for some cases where there is trivial overlap, typically at one or two intersections. No other solutions that make them both equal are possible. However, by using the differentiation to functional order approach, it will be shown that not only can they be made equal, but one can asymptotically approach the other. Let, The left-hand side of (42) for the Exponential-density involves the fractional differentiation of a negative exponential term, which was performed previously; see (7). Then, replacing κ with λ provides the following form for (42), Observe that the conventional Exponential-density is recovered in (43) when α(x) = α 0 = 0. Furthermore, (43) holds if, and only if, there is a functional order α(x) that allows for this equivalence. In fact, this is true when The substitution of (44) into the left of (43) collapses to the right side, i.e., the Pareto density. The functional-order (44) can be rewritten as follows: If α 0 = 1 in (45), this gives (44) and, as a result, (43) holds. Varying α 0 above and below unity provides the asymptotic solutions to the Exponential density that, in the limit, α 0 → 1, tend towards and exactly match the Pareto density. Figure 3 shows plots of the Pareto density for parameters x 0 = 1/2, β = 2 and the Exponential density parameter λ = 3 for varying values of α 0 , providing the exact curve that is produced by the Pareto density when α 0 = 1, and asymptotic solutions, as shown for α 0 = 0.8 and α 0 = 1.2. Hence, via differentiation to functional order, a one-parameter density (Exponential model) produces the same mathematical properties as a two-parameter density (Pareto model). In this case, the Exponential density was not only generalised, but can also produce specific solutions corresponding to a different density, namely the Pareto, while also producing the asymptotic solutions that approach the latter in the limit α 0 → 1. Any α 0 value that is greater than α 0 = 1, such as the one considered in Figure 3, namely, α 0 = 1.2, will provide solutions that asymptotically approach the Pareto density from the right. As can be seen, solutions with α 0 < 1 values approach the Pareto density from the left and, as α 0 → 0, solutions asymptotically approach the Exponential density, and equal it when α 0 = 0 (not shown in Figure 3 for reasons of brevity). The asymptotic solutions, as given by different values of α 0 , except for α 0 = 0 and α 0 = 1, are, in fact, mixed-probability densities consisting of a mixture of both the Exponential and Pareto densities, respectively. Figure 3. Plots of the Pareto density and the Exponential density. The latter is exact with the Pareto when α 0 = 1, while the asymptotic solutions of the Exponential density that approach the Pareto density are also shown to correspond to α 0 = 0.8 and α 0 = 1.2, respectively. Note that if α 0 = 0, then the curve would be equivalent to the traditional Exponential density (not shown).

Duality and Asymptotic Mixtures between the Gibbs-Shannon Entropy and the Hellinger Distance
As in the previous examples, consider another important application of differentiation to functional order that is of interest to physics, mathematics and engineering, including information geometry and information theory. Using the previously discussed differentiation method to functional-order, it will be shown that there is a duality between the Gibbs-Shannon entropy for a continuous probability density p(x; ξ m 1 ) and a metric distance between the same probability density and another density q(x; ξ m 2 ) with parameter space ξ m 1 = (x 1 1 , x 2 1 , . . . ) and ξ m 2 = (x 1 2 , x 2 2 , . . . ), respectively, with m = 1, 2, 3, . . . , M and M providing the total number of parameters for each density. It has been shown that the fractional entropy for probability densities is given by [6]: where Ω is the domain of integration of the density and the modulus is used because the fractional-order α is not only real but can also be imaginary, α ∈ R and α ∈ C. It is worth pointing out that the analogue of (46) for the discrete case of probability masses was derived by Ubriaco [41]. The dissimilarity or divergence of the two probability densities can be determined using a metric distance, referred to as the Hellinger distance: where (47) is the square of the Hellinger distance. It should be obvious that both (46) and (47) will not provide the same result, as one is a measure of entropy and the other is a measure of the distance between the densities. Using differentiation to functional-order, it will be shown that there is a duality between the two formulations and both can give the same results. This duality will be established even though they each measure different properties between density/s: one being entropy and the other distance or divergence. In practice, this duality does not exist because the Gibbs-Shannon entropy is not a true distance metric as it is not symmetric and only measures entropy. On the other hand, the Hellinger formulation is symmetric and does obey the triangle inequality but measures distance and not entropy. Let the two densities be given as the Exponentials: where their respective parameter spaces consist of only one parameter, namely, (ξ 1 1 = u) and (ξ 1 2 = v). Equation (46) is already in fractional form, so there is no need to obtain the fractional derivative. If α = α(x) is substituted, this gives: for the fractional entropy. A duality between the entropy and the Hellinger distance can be claimed for some functional-order α(x), such that: In fact this is true and (50) holds if when α 0 = 1. It can be shown that when (51) is substituted into the entropy on the left of (50), it collapses to the Hellinger expression on the right, hence establishing a duality between them. This duality was achieved via the application of functional-order differentiation to conventional entropy. Numerical solutions can be obtained for the parameter values u = 2.5 and v = 4.0, for example. The conventional entropy then gives: On the other hand, the square of the Hellinger distance between the densities p(x; u) and q(x; v) cannot be the same as the entropy for the density p(x; u), as expected. The square of the Hellinger distance becomes: Using the fractional entropy and the functional-order α(x), as given by (51), with α 0 = 1: gives the value: in agreement with the Hellinger distance (53). Hence, a duality is established between them. Notice that if the square root of the right side is taken in (53), then the same can be taken for (55); they are equal either way. Asymptotic solutions between the Gibbs-Shannon and Hellinger formulations can also be obtained using (51) for varying α 0 , which will be referred to as asymptotic divergences. When α 0 = 1, the fractional entropy is equivalent to the Hellinger distance or divergence; that is, when α 0 = 1, as established above. Setting α 0 = 3.099 in (51) provides the conventional entropy result, i.e., Figure 4 shows numerical values for the asymptotic divergences between the conventional entropy and the square of the Hellinger distance as a function of α 0 . The numerical asymptotic divergence values that are bounded by the conventional entropy and Hellinger distance were plotted for α 0 ∈ [1, 3.099] in steps of δα 0 = 0.1. The first value on the left, corresponding to α 0 = 1, gives the exact solution as obtained from the Hellinger distance. As α 0 varies and becomes equal to the last value on the right, i.e., when it is α 0 = 3.099, the solution is exactly the same as that given by the conventional entropy. For α 0 values that lie between them, asymptotic divergence solutions are obtained, which are due to a 'mixture' of the two, depending on how close α 0 is to one or the other. As a result, a minimum asymptotic divergence value exists between them when α 0 = 1.4. This indicates, among other things, where the conventional entropy and Hellinger distance have similar mathematical characteristics, while also minimising their respective probability densities. Finally, returning to the fact that (49) and (51) establish a duality between the conventional entropy and the Hellinger distance when α 0 = 1, it stands to reason that there can be other valid solutions, depending on the mathematical form of α(x). For example, if the functional order is a constant value instead of a function, i.e., if α(x) = 6.9286 + 3i, for example, then it is possible to obtain the same result for the fractional entropy as when the functional order α(x) was a continuous function (51). Substituting this into (55) gives: as before. This means that the entropy can also obtain the same result as the Hellinger formulation when the functional order is equal to a complex constant and not just a function. By differentiating the entropy to functional order, it was possible to find a duality between the conventional entropy and the Hellinger distance, as well as their asymptotic divergence solutions (mixture), not withstanding the fact that the two are completely different mathematical formulations.

Dualities between Higher-Dimensional Geometries and the Functional Order in Quantum Mechanics
Higher-dimensional spheres or hyperspheres have been of tremendous interest in the study of the effective properties of materials, since they exhibit interesting electromagnetic behaviour when they contain these hyperspheres as inclusions. The scattering effects of hyperspheres have been studied, and so has their use in metamaterial structures for example [42]. Hyperspheres, and other hypergeometries such as hypercubes, are used in the study of error rate and data correction algorithms or in packing geometries. Functionalorder differentiation will be used to examine the dualities between the functional order and the dimension of hyperspheres and hypercubes. These hypergeometries and, in particular, the geometry of hyperspheres, can represent quantum mechanical particles and their classical and relativistic properties. On this basis, consider the functional derivative V αd (R) of the hypervolume of a hypersphere, which can be calculated from: where V d (R) is the hypervolume of a hypersphere. The modulus is included because, in principle, α ∈ R and α ∈ C. From this, the functional variation in the hypervolume of a d-dimensional sphere (hypersphere) is given by: where d is the dimension, α is the functional order and R is the radius of the hypersphere. The functional change in the hypervolume of a hypercube is obtained as follows: where L is the side length of the hypercube. The radius of a hypersphere R is related to the length L of a hypercube by the following: where d = α. The functional order, which provides equal hypervolumes between hyperspheres and hypercubes is given by: Hence, for the case represented in Figure 5, the functional order can be determined from (63) and is α(3.1827) = 1.06088. The hypervolume of the hypersphere starts from above 80 on the left in magenta and the hypervolume of the hypercube starts below 40 on the left in blue. The volume of a regular cube is constant as represented by the horizontal orange line. For the given parameters, the functional order α = 1.06088 makes the hypervolumes/volume equal. The functional order α(d) given by (63) not only provides the duality or connection between the hyperspherical and hypercubic geometries, but also does the same for the relativistic energies of particles with such geometries. The relativistic energy of a hyperspherical particle is given by the following: where ρ is the density of the mass, c is the speed of light and u is the speed of the hypersphere, respectively. Similarly, the energy of a hypercube particle is given by the following: The rest energy E 0 of a hypersphere or hypercube particle is given by (64) and (65) when u = 0, respectively. For non-relativistic speeds u << c, the energies reduce to E HS (HyperSphere) and E HC (HyperCube), respectively: and where E HS 0 is the rest energy (64) and E HC 0 is the rest energy (65), i.e., when u = 0. If a hypersphere and a hypercube particle have the same density and the same velocity, the ratio of their energies E = E HS /E HC is given by the following: (68) Figure 6 shows the relationship between the rest energy, functional order and fractal dimension of hypersherical and hypercube particles as a comparison. Referring to Figures 7 and 8, consider the case where an electron is assumed to be a hypersphere (functional order α = 0). When the electron is assumed to have a spherical geometry, d = 3, the corresponding energy agrees with the theory. When the same electron is assumed to have the energy of a typical open superstring of Planck length and energy, the dimension is fractal at d = 1.33. Superstrings are considered to be approximately d = 1 dimensional from the macroscopic perspective but, in reality, they possess an internal geometry that arises from their vibration in Calabi-Yau manifolds. The figures show that the fractal dimension of the superstring (electron), namely, d = 1.33, 'drops' out from extrapolating the dimension as the energies increase and match the expected values from the theory; see Figure 7 for the particle at rest and Figure 8 when it is moving at 95% c. Revisiting the hypersphere volume (60), a hyper-particle of radius R = 1 and even dimensions provides a hypervolume variation, where d = 0, 1, 2, . . . and α ∈ R + . Summing (69) to infinity, i.e., summing all even dimensional hypervolumes to infinite dimensions, gives the closed form: where 2F2 (·) is the regularised generalised hypergeometric function. This means that the sum of all relativistic energies for even dimensional hyperparticles of radius one is given by the following: If the velocity of the hyperparticles is u << c, then: An interesting result is obtained for the sum of the rest energies of all the evendimensional hyperparticles. If the hypervolume is considered, i.e., α = 0 when u = 0 in (71) or (72), the total energy becomes: This can easily be seen if one considers the sum in (70) ∞ ∑ d=0 π d Γ(1 + d) = 1 + π + π 2 2! + π 3 3! + . . . = e π (74) In a similar way, for α = 1 and α = 2, the energies become: and E ∞ = 2πρc 2 (2π + 1)e π (76) respectively, etc. These results are very interesting and raise a number of further questions that need to be answered, such as determining their physical interpretation.
(a) (b) Figure 6. The rest energies of hypersphere (a) and hypercube (b) geometries as a function of functional order and fractal dimension. Note that the energy is given in units of EJ or exa-joules (10 18 ) Joules.

Conclusions
Variable-order fractional calculus is an extension of the standard constant order to a function containing arbitrary parameters. Unfortunately, these parameters can make the calculations intractable because they convolute the integrals that are present in the operators. The idea of a delayed variable order or differentiation/integration to functional order is an alternative that has been used to solve some problems, which avoid these issues. This allows for the easy derivation of asymptotic solutions between different functions and can establish dualities between them that are not otherwise possible. Such dualities and their asymptotic solutions beg the question as to whether there are further mathematical properties or 'connections' that exist between seemingly independent mathematical expressions or theories. If the variable in the functional order is time α(t), then the solutions have memory properties: past, present and future. If the functional order contains a space variable α(x), then it has spatial correlation properties, a kind of spacial dimensional fading effect. These memory and spatial correlation properties are equivalent to the asymptotic solutions between two or more functions. Aside from the direct dualities, asymptotic solutions can be obtained between two or more functions that are essentially nothing more than mixtures of these functions. It is hoped that the use of functional-order calculus to solve problems, as presented in this paper, will encourage more researchers to use this approach in their work.