This work is licensed under the Creative Commons Attribution 4.0 International License.
Introduction
The multifractal formalism introduced by Halsey et al [1] can be understood in a simple way by applying a similar reasoning that was used by Boltzmann for obtaining the thermodynamics of an ideal gas using statistical arguments instead of the microscopic description of a system conformed by 1023 particles. The central idea was to introduce the relation between the entropy and the probability associated with a macrostate [2, 3]. In the second section recalls briefly the Boltzmann fundamental ideas for obtaining the statistical description of an ideal gas in thermodynamic equilibrium. In the third section, we obtain the Eggleston’s theorem, which relates the Hausdorff dimension with the Shannon entropy; this theorem plays a similar role in fractals like that the relation between entropy and probability in the Boltzmann treatment. The theorem is showed using a multiplicative process to decompose the unitary interval in fractals M(φ⃗), conformed by points with the same frequency of digits φ⃗; evaluate the Hausdorff dimension of M(φ⃗) and obtain that it is related to the Shannon entropy. In the fourth section, we introduce a Bernoulli measure with a probability vector p⃗; make the multifractal decomposition in terms of sets Jα, conformed by points with the same pointwise dimension α, and show that they are conformed by an infinite number of sets M(φ⃗). Therefore, each Jα has a multifractal structure, to determine the Hausdorff dimension D(α), we use a basic property of the Hausdorff Dimension which plays the rule of the Boltzmann principle of maximum entropy.
In the fifth section, the Boltzmann procedure is extended to determine the distribution P⃗(q), which maximize the Eggleston’s relation, given a value of α, D(α) is determined using P⃗(q), and D(α) singularity spectrum is obtained.
In the sixth section, we introduce a family of Bernoulli q-measures with a probability vector P⃗(q). We evaluate the q-measure of the sets Jα, finding that for a particular value of q the measure is concentrated in the set Jα*, as a consequence of the singular behaviour of the q-measure, we obtain that τ(q) and -D(α) are Legendre transform of each other. The explicit values α* = α(q) and the Hausdorff dimension D(α(q)) are the functions obtained in the Boltzmann procedure to determine the D(α) singularity spectrum.
Fundamental Boltzmann ideas
Boltzmann analyzed a system conformed by N particles which interact only through elastic collision. For establishing a relation between the mechanics and thermodynamics, Boltzmann introduced a probabilistic description of this system dividing in K ⋘ N cells the phase space of a particle, and described the state of the system by the occupation number of particles of these cells:
The time evolution of the probability distribution is governed by the Boltzmann’s equation. When it is introduced in (4), it can be proved that the entropy increases until it obtains its maximum value for a stationary distribution, which is the Maxwell-Boltzmann distribution. However, it is possible to obtain this result without invocating the Boltzmann equation, using that in the equilibrium state the entropy of the system obtains its maximum value under the restrictions:
Then maximizing (4) with the lateral conditions (5), it is found that the stationary distribution is given by:
$$\begin{array}{}
\displaystyle \tilde{p}_i = \frac{1}{Z} e ^{-\beta\varepsilon_i}; \quad with \quad Z=\sum\limits_{i=1}^{K} e ^{-\beta\varepsilon_i}
\end{array} $$
The value of parameter $\begin{array}{}
\beta=\frac{1}{kT}
\end{array} $ is determined by the thermodynamic information that the internal energy of the ideal gas is given by $\begin{array}{}
E=\frac{3}{2}NkT,
\end{array} $ and (6) reduces to the Maxwell-Boltzmann distribution.
The geometric multifractal decomposition of the unitary interval 𝕀
In this section we discuss the multifractal decomposition of the unitary interval 𝕀. Any real number x ∈ 𝕀, is expressed in s-base as:
Let ni(x, K) denote the number of times the digit i ∈ (0, 1, …, s − 1) occurs among the first K digits of x. Then, the frequency in which this digit appears in x is given by:
$$\begin{array}{}
\displaystyle \lim_{ K \to \infty} \frac{n_i(x,K)}{K} = \lim_{ K \to \infty} f_i (x,K)={\boldsymbol\varphi}_i(x), \quad i=0,1,\ldots,s-1
\end{array} $$
where $\begin{array}{}
0 \leq {\boldsymbol\varphi}_i \leq 1, \sum\limits_{i=0}^{s-1} {\boldsymbol\varphi}_i =1;
\end{array} $ thus, x has associated a frequency vector:
Let M(φ⃗) be the set of points in 𝕀 with the same frequency vector. For obtaining the multifractal decomposition of 𝕀, we separate the unitary interval in the different sets M(φ⃗), and using the Eggleston’s theorem [4], we evaluate their Hausdorff Dimension. This can be done using a multiplicative process that consists of dividing 𝕀 in s-cylinders of first order:
where fr is the frequency that shows that digit r = (0, 1, …, s − 1) occurs in σK. The number of K-cylinders with the same frequency vector f⃗ is given by
due to the fact that when K → ∞, each K-cylinder goes to a point x of the unitary interval with a frequency vector given by (8). The Hausdorff dimension of M(φ⃗) is
This relation is the Eggleston’s theorem [4]. In the Appendix A, we show how (15), can be found using the definition of Hausdorff measure applied for dyadic intervals [6]; we find the value D for which this measure is non-singular, and it corresponds with the Hausdorff dimension of the set.
On the other hand, the relation (15) establishes a non-trivial connection between Hausdorff dimension and the Shannon entropy, which is discussed in the Appendix B.
Statistical multifractal decomposition of the unitary interval 𝕀
When a statistical measure is assigned to each point of 𝕀, it is decomposed in subsets Jα conformed by points with the same pointwise dimension. A simple case is when we introduce the Bernoulli measureμ on the unitary interval with probability vector p⃗ = (p0, p1, …, ps−1); assigning to each digit j belongs to x a probability pj, it introduces a singular measure that can be characterized by the pointwise dimension of μ at x [5]:
We note that all the points that belong to M(φ⃗) have the same pointwise dimension. However, there are an infinite number of setsM(φ⃗) with the same value of the pointwise dimension, because given a particular value of dμ (x) = α and the normalization condition $\begin{array}{}
\sum\limits_{j=0}^{s-1} {\boldsymbol\varphi}_j(x)=1,
\end{array} $ we cannot determine the s components of the frequency vector.
As each M(φ⃗) is a fractal with the Hausdorff dimension given by Eggleston’s theorem, then Jα, the set of points with dμ (x) = α, is a multifractal:
where was defined: $\begin{array}{}
\vec{\boldsymbol{a}} \cdot \ln \vec{\boldsymbol{b}} = \sum\limits_{j=0}^{s-1} a_j \ln b_j
\end{array} $ as the Hausdorff dimension satisfies that [6]
The statistical multifractal decomposition of 𝕀 consists of grouping the points x in subsets with the same value of the pointwise dimension, and each subset Jα is characterized by its Hausdorff dimension D(α).
Boltzmann scheme for Multifractals
We determine D(α) using the Eggleston’s theorem and the extremal principle given by (22); this procedure is similar to the maximum entropy principle used by Boltzmann for obtaining the stationary distribution characterizing the equilibrium state. The Hausdorff dimension D(α) is determined by the frequency distribution φ⃗* that maximizes:
The dimension spectra for the statistical multifractal decomposition of the unitary interval is found when the q parameter is eliminated for (25) and (26). In thermodynamics, the entropy is one of the relevant functions, but there are several functions that contain the same thermodynamic information, they are the thermodynamic potentials, we show that in multifractals, a similar situation occurs. From (24) we have:
This result infers that dD = qdα, therefore d(D − qα) = −dτ = −αdq, which implies (30) and hence τ (q) and -D(α) are the Legendre transform of each other.
Statistical q-measures of Jα
In the statistical multifractal decomposition of the unitary interval, we focused our attention to the Hausdorff dimension of the sets Jα, which is a geometrical property of these sets. However, they also have statistical properties, because they are the support of the q-measures, in such way that given a q value, this measure is supported by only one of the Jα sets. Given a probability vector p⃗, a set of probability vectors P⃗(q) can be constructed, these vectors have the ability to scan the structure of the multifractal decomposition [8]. We obtained the escort probability vector P⃗(q) in the Boltzmann scheme for multifractals, they are equivalent to the Maxwell-Boltzmann distribution in statistical physics, where the q value plays the rule of the inverse of the temperature. The statistical q-measures are Bernoulli measures in the unit interval generated by P⃗(q), which is defined by (24), i.e.
The set of all Kα-cylinders contains all the points of 𝕀 with dμ (x) = α, therefore this set conforms the cover CK(Jα) of the set Jα. For large values of K we have that:
where the infimum is taken with respect to α; thus, (42) defines τ(q) as the Legendre transform of −D(α). We note that this result is obtained from (40), which is the generalization of the random weighted curdling proposed by Mandelbrot [9] and [10].
The relations (41) and (42) imply that α*, satisfies the following conditions:
Therefore, the q-measure is concentrated in the set Jα(q) with α(q), which is given by (45), and it can be rewritten as the following average on P⃗(q):
When a Bernoulli statistical measure characterized by a probability vector p⃗ is introduced in a fractal, it is decomposed into sets Jα, which are multifractal. The determination of their Hausdorff dimension D(α) requires to use an extremal property of the Hausdoff dimension, similar to the maximum entropy principle. D(α) is determined in terms of a probability distribution P⃗(q), we find that each set Jα is an statistical attractor set where the q-measure defined in terms of P⃗(q) is concentrated.
As a consequence of the singular behaviour of the q-measure on the sets Jα, given by (40), we obtain the following:
τ(q) and −D(α) are Legendre transform of each other and,
the information on the set Jα where the q-measure is supported, and its Hausdorff dimension, are given by (25) and (26), respectively.
An important case of (40) is when q = 1. The 1-measure is generated by the probability vector p⃗, the statistical attractor or the curdling set is conformed by the points with their pointwise dimension is identical with the Hausdorff dimension of the set, and it is given by Shannon entropy of p⃗, i.e.
Then, D(q = 1) is the Hausdorff dimension of the measure theoretical support of p⃗, which was found by Billingsley [6] in his work about the Hausdorff dimension in probability theory. Chhabra and Jensen [11] applied this result to P⃗(q) and found (47), after using heuristic arguments introduces (46), and with these expressions they found an alternative method for obtaining the singularity spectrum. On the other hand, Mandelbrot [10] introduced the curdling set for explaining the energy dissipation in fully developed turbulence using a multiplicative cascade process and identified this set with the Besicovitch fractal [9] extended the Mandelbrot suggestions, Feder [7] obtained and showed that (48) characterizes the set where the 1-measure is concentrated. The result (40) can be showed for a singular measure, and obtain an unified description of the multifractal decomposition can be obtained, which relates the Halsey et al [2] and Chhabra and Jensen [11] methods for obtaining the spectral singularity (see sections 5 to 7 of reference [12]).