A Note on the Connection between Non-Additive Entropy and h-Derivative

In order to study as a whole a wide part of entropy measures, we introduce a two-parameter non-extensive entropic form with respect to the h-derivative, which generalizes the conventional Newton–Leibniz calculus. This new entropy, Sh,h′, is proved to describe the non-extensive systems and recover several types of well-known non-extensive entropic expressions, such as the Tsallis entropy, the Abe entropy, the Shafee entropy, the Kaniadakis entropy and even the classical Boltzmann–Gibbs one. As a generalized entropy, its corresponding properties are also analyzed.


Introduction
Since it was proposed over one hundred years ago, the conventional Boltzmann-Gibbs (BG) statistics has been developed very delicately and successfully with wide applications in many disciplines. During the last few decades, however, people noticed that more and more systems are difficult to be described by this simple BG distribution, such as the longrange interactions [1], the gravitational systems [2], the Lévy flights and fractals [3], and so on [4]. In order to cope with this challenge, some attempts have been made to generalize the BG statistics. Among them, the most investigated formalism is the non-extensive entropy. It was inspired by the geometrical theory of multi-fractals and its systematic use of powers of probabilities by C. Tsallis [5]: where k B is the Boltzmann constant (hereafter we assume k B = 1 for simplicity) and q stands for the Tsallis non-extensive parameter. It describes the departure of non-extensive statistics from the BG one. This entropy goes back to the usual BG form when q → 1. For more than two decades of researches and developments, the Tsallis entropy has been successfully applied to various domains: physics, chemistry, economics, computer science, biosciences, linguistics, and so on [5][6][7][8][9][10]. For the average charged-hadron yields in inelastic non-singlediffractive events, V. Khachatryan et al. observe it as the Tsallis distribution [11,12] E d 3 N ch dp 3 = where E d 3 N ch dp 3 is for the function of spectra with E for the total energy of the particle and p for its momentum, η denotes the pseudorapidity with y for the rapidity, p T stands for the transverse momentum, C is for its normalization constant, T, a variational parameter representing the temperature when the system reaches equilibrium, n is the fitting parameter which connects with Tsallis' q by n = 1/(1 − q), y = 0.5 ln[(E + p z )/(E − p z )], E T = m 2 + p 2 T − m, and m is the charged pion mass. The data fitting results show that the Tsallis distribution can well-describe both the low-p T exponential and the high-p T power-law behaviors [11,12]. One application in astrophysics is the study of the distribution of asteroid rotation periods from different regions of the solar system and diameter distributions of near-Earth asteroids (NEAs) [13]. A. S. Betzler and E. P. Borges analyze two samples from different years. They discover that the distribution of diameters of NEAs obeys a Tsallis-like distribution, and the rotation periods of asteroids can be well-approximated by a Tsallis-Gaussian function. According to the first conclusion, there should be 994 ± 30 NEAs with diameters greater than 1 km [13]. In another example, Y. Wang and J. Du study the viscosity of light charged particles in weakly ionized plasma with the power-law Tsallis-distributions using the generalized Boltzmann equation of transport and the motion equation of hydrodynamics [14].
The Tsallis entropy is indeed not unique. By now, a lot of different expressions of the non-additive entropies have been proposed, for instance, the Kaniadakis entropy [15], the Shafee entropy [16], the q − q −1 symmetric modification of Tsallis entropy [17], and the two-parameter (q, q )-entropy [18]. These expressions were obtained in quite different ways and investigated by distinct motivations. Therefore, it will be of great interest to find the relationship among these formulas or to find a simple formula to study them as a whole.
In this paper, we first introduce a two-parameter non-additive entropy, S h,h , based on the h-derivative. The h-derivative is known as a mathematical generalization of the normal Newton-Leibniz calculus. We address that S h,h unifies different types of expressions of non-extensive entropies; namely, it can connect a family of non-extensive entropies. On the other hand, we also discuss its properties in order to better understand this newly established non-additive entropic function.

h-Derivative
In the conventional mathematical theory, the Newton-Leibniz derivative is defined as: Classically, most of the physical quantities are continuous, and it is natural to apply the Newton-Leibniz derivative. In quantum physics, on the other hand, all the physical quantities will be quantized; people then try to develop quantum calculus, which utilizes the discrete forms of derivatives instead and presents a generalization of this Newton-Leibniz derivative. One formalism of quantum calculus is the h-derivative [19]. For an arbitrary function f (x), its h-differential is defined as follows: It is easily verified that and Thus, can we obtain the h-derivative of f (x): When f (x) is differentiable, the following property is obviously obtained: which is nothing but the definition of the conventional Newton-Leibniz derivative. Note that we need the function f (x) to be continuous for the Newton-Leibniz derivative, but this requirement becomes unnecessary for the h-derivative. Next, some basic rules of this h-derivative are listed: 1. Sum and difference Considering the sum and difference rules of the h-derivative, we have 2.
Product and quotient rules As for the product and quotient rules, 3. h-derivative of elementary functions Some other basic calculations of it are expressed: In Figure 1, we illustrate the behavior of D h e x at different values of h as an example. We could see that it behaves as an exponential when h = 0. For any fixed values of h, D h e x is a monotonically increasing function with respect to the variable x. The values of this derivative also increase when the parameter h becomes larger.
With the definition of h-derivative, V. Kac and P. Cheung [19] developed a type of quantum calculus, known as h-calculus. As a matter of fact, an operator such as h-derivative is called the forward difference quotient operator. Analogously, it also has the backward difference quotient operator ∇ h and the central difference quotient operator δ h , defined as [20] Note that the regular vector differential operator ∇ has been generalized based on hderivative. We then explore the connection between the h-derivative entropy and its modified forms.

h-Derivative and Non-Additive Entropy
In order to generalize the non-additive entropic forms, one could utilize Equation (7) and give out the following equation: with the normalization condition ∑ W i=1 p i = 1. When h → 0, it will go back to the usual BG entropy. Note that it also recovers the Tsallis non-extensive entropy, S q , cf. Equation (1) under the transformation of h = q − 1.
Following the ways of the central difference quotient operator of Equation (20), we define a new form of two-parameter (h, h )-derivative, The corresponding (h, h )-entropy is Similarly, when h = h → 0, the entropy S h,h returns to the BG one. It is shown that, Note that L'Hospital's rule has been applied within the formula de αx /dx = αe αx for the last step in the above. It is constructive to explore the connections with the already known statistical distributions. For example, the Tsallis entropy is obtained by h = q − 1, h = 0. While taking h = q − 1, h = 1 − q −1 , we can obtain the Abe entropy Equation (A5) (see the discussion in the Appendix A) [17]. The non-extensive entropy given in Equation (A7) proposed by Borges and Roditi [18] (also see the Appendix A) can then be recovered with the relationship of h = q − 1, h = 1 − q . Although the non-extensive entropy of Borges and Roditi and our two-parameter (h, h )-entropy have similar forms, we gained them using different mathematical methods. Specifically, we used (h, h )-derivative developed by ourselves, which differs from the q-calculus used by Borges and Roditi. In addition to the difference in the form of expression between the two-type derivative, a conspicuous point is that our two-parameter (h, h )-derivative does not require the function f (x) to be continuous and differentiable at x = 0.
By assuming h = h, we could also obtain another new form of entropy S h,h , which is obviously invariant under the interchange h ↔ −h. As a matter of fact, it is nothing new but the well-known Kaniadakis non-extensive κ-entropy [15]. Last but not least, it is set that h → −h and h = −h + δ. Considering the limit of δ → 0 and h → −h, we could also cover the exact Shafee entropy [16,21] by taking the transformation of q = h + 1.
In Table 1, we summarize different entropy functions, which can be represented by this two-parameter S h,h entropy through taking different values of h and h . In addition, by choosing h = −1/h the function S h,h becomes Note that this entropic form looks much similar to Abe entropy [17], but it is totally different in fact that Abe entropy cannot be recovered only by exchanging q and h when comparing them. Hereby, we name it the modified-Abe entropy function. Except for the entropy forms listed in Table 1, there is a well-known entropy-Renyi entropy, which can be related to S h,h through the relationship between Renyi entropy and Tsallis entropy (only for q ≤ 1) [22],  [16,21,24] h → −h modified Abe h = −1/h

Properties
Now we shall address some properties of this (h, h )-entropy, S h,h . As we all know, the Boltzmann-Gibbs and the Tsallis entropy can be expressed as [6,25] where . . . ≡ ∑ W i=1 p i (. . .) is the standard mean value, and ln q is q-logarithm. Along this line, we straightforwardly obtain where ln h,h is the (h, h )-logarithm, and it can be expressed as

Non-Negativity
First of all, we consider a thermal system within any possible state. The probability distribution of each microstate i is defined as p i . If we assume p 1+h i p 1−h i , namely, 1 + h 1 − h , for 0 p i 1, thus can we obtain h + h 0 and this two-parameter entropy S h,h 0.

Extremal at Equal Probabilities
Utilizing the Tsallis entropy, , this two-parameter entropy S h,h can be expressed with it as, For the Tsallis entropies inside this formula, namely S T 1+h and S T 1−h , it is easy to know that both of them reach their extreme values when all the probabilities are equal [6]. Therefore, at the state of equal probability, our entropy S h,h will also approach to its extreme value since
This property trivially follows from the definition itself. It means when we add some events with zero probabilities, S h,h keeps invariant.

Non-Additivity
When we consider a system that can be decomposed into two independent subsystems, A and B, The values of h and h cannot be zero together in case (or h = −1 and h = 1 appear at the same time). In other words, S h,h is said to be non-additive similar to the Tsallis non-extensive entropy.

Summary and Outlook
To summarize, with the generalized h-derivative we firstly propose a two-parameter non-additive entropy, S h,h , in order to connect several non-extensive entropy functions. The h-derivative motivated non-additive entropy, S h,h , is demonstrated to recover different kinds of non-extensive entropy formalisms, such as the Tsallis entropy (h = q − 1, h = 0), the Abe entropy (h = q − 1, h = 1 − q −1 ), the Borges-Roditi entropy (h = q − 1, h = 1 − q ), the Kaniadakis κ-entropy (h = h = κ) and the Shafee (h → −h or h = −h + δ, here δ → 0) non-extensive entropy by varying values of h, h . On the other hand, the present two-parameter entropy exhibits all the relative properties as a generalized nonextensive entropy. Furthermore, the remarkable relationship between S h,h and other nonextensive entropies may cast a light on the connection of non-extensive entropy and some mathematical structures such as quantum calculus. It may lead to a deeper understanding of the mathematical and physical foundations of non-extensive statistics. We also noticed some other two-parameter distribution functions, such as the (r, q)-distribution and (α, κ)distribution [26,27], which have been well-applied to astrophysics or space plasma physics. These two-parameter distributions provide another view to investigate the non-Maxwellian systems. It will be of great interest to associate this (h, h )-entropy with them and find out the deeper connections. There are also various elegant forms of entropy, such as fractional entropy [28] and Deng entropy [29]. Our (h, h )-entropy, S h,h , is indeed unable to establish a connection with theirs. Further exploration of the inherent connections between different forms of entropy is necessary.