Unifying Aspects of Generalized Calculus

Abstract

Non-Newtonian calculus naturally unifies various ideas that have occurred over the years in the field of generalized thermostatistics, or in the borderland between classical and quantum information theory. The formalism, being very general, is as simple as the calculus we know from undergraduate courses of mathematics. Its theoretical potential is huge, and yet it remains unknown or unappreciated.

1. Introduction

Studies of a calculus based on generalized forms of arithmetic were initiated in the late 1960s by Grossman and Katz, resulting in their little book Non-Newtonian Calculus [1,2,3]. Some twenty years later, the main construction was independently discovered in a different context and pushed in a different direction by Pap [4,5,6]. After another two decades the same idea, but in its currently most general form, was rediscovered by myself [7,8,9,10,11,12,13,14,15]. In a wider perspective, non-Newtonian calculus is conceptually related to the works of Rashevsky [16] and Burgin [17,18,19,20] on non-Diophantine arithmetics of natural numbers, and to Benioff’s attempts [21,22,23,24,25] of basing physics and mathematics on a common fundamental ground. Traces of non-Newtonian and non-Diophantine thinking can be found in the works of Kaniadakis on generalized statistics [26,27,28,29,30,31,32,33,34]. A relatively complete account of the formalism can be found in the forthcoming monograph [35].

In the paper, we will discuss links between generalized arithmetics; non-Newtonian calculus; generalized entropies; and classical, quantum, and escort probabilities. As we will see, certain constructions such as Rényi entropies or exponential families of probabilities have direct relations to generalized arthmetics and calculi. Some of the constructions one finds in the literature are literally non-Newtonian. Some others only look non-Newtonian, but closer scrutiny reveals formal inconsistencies, at least from a strict non-Newtonian perspective.

Our goal is to introduce non-Newtonian calculus as a sort of unifying principle, simultaneously sketching new theoretical directions and open questions.

2. Non-Diophantine Arithmetic and Non-Newtonian Calculus

The most general form of non-Newtonian calculus deals with functions A defined by the commutative diagram ($f_{\mathbb{X}}$ and $f_{\mathbb{Y}}$ are arbitrary bijections)

$$\begin{array}{c}\hfill \begin{array}{ccc}\hfill \mathbb{X}& \stackrel{A}{⟶}& \mathbb{Y}\hfill \\ \hfill f_{\mathbb{X}}\downarrow & & \downarrow f_{\mathbb{Y}}\hfill \\ \hfill \mathbb{R}& \stackrel{\tilde{A}}{⟶}& \mathbb{R}\hfill \end{array}\end{array}$$

(1)

The only assumption about the domain $\mathbb{X}$ and the codomain $\mathbb{Y}$ is that they have the same cardinality as the continuum $\mathbb{R}$. The latter guarantees that bijections $f_{\mathbb{X}}$ and $f_{\mathbb{Y}}$ exist. The bijections are automatically continuous in the topologies they induce from the open-interval topology of $\mathbb{R}$, even if they are discontinuous in metric topologies of $\mathbb{X}$ and $\mathbb{Y}$ (a typical situation in fractal applications, or in cases where $\mathbb{X}$ or $\mathbb{Y}$ are not subsets of $\mathbb{R}$). In general, one does not assume anything else about $f_{\mathbb{X}}$ and $f_{\mathbb{Y}}$. In particular, their differentiability in the usual (Newtonian) sense is not assumed. No topological assumptions are made about $\mathbb{X}$ and $\mathbb{Y}$. Of course, the structure of the diagram implies that $\mathbb{X}$ and $\mathbb{Y}$ may be regarded as Banach manifolds with global charts $f_{\mathbb{X}}$ and $f_{\mathbb{Y}}$, but one does not make the usual assumptions about changes of charts.

Non-Newtonian calculus begins with (generalized, non-Diophantine) arithmetics in $\mathbb{X}$ and $\mathbb{Y}$, induced from $\mathbb{R}$,

$$\begin{array}{ccc}\hfill x_1{\oplus }_{\mathbb{X}}{x}_2& =& f_{\mathbb{X}}^{-1}\left(f_{\mathbb{X}}\left(x_1\right)+f_{\mathbb{X}}\left(x_2\right)\right),\hfill \end{array}$$

(2)

$$\begin{array}{ccc}\hfill x_1{\ominus }_{\mathbb{X}}{x}_2& =& f_{\mathbb{X}}^{-1}\left(f_{\mathbb{X}}\left(x_1\right)-f_{\mathbb{X}}\left(x_2\right)\right),\hfill \end{array}$$

(3)

$$\begin{array}{ccc}\hfill x_1{\odot }_{\mathbb{X}}{x}_2& =& f_{\mathbb{X}}^{-1}\left(f_{\mathbb{X}}\left(x_1\right)·f_{\mathbb{X}}\left(x_2\right)\right),\hfill \end{array}$$

(4)

$$\begin{array}{ccc}\hfill x_1{\oslash }_{\mathbb{X}}{x}_2& =& f_{\mathbb{X}}^{-1}\left(f_{\mathbb{X}}\left(x_1\right)/f_{\mathbb{X}}\left(x_2\right)\right)\hfill \end{array}$$

(5)

(and analogously in $\mathbb{Y}$).

Example 1.

According to one of the axioms of standard quantum mechanics, states of a quantum system belong to a separable Hilbert space. All separable Hilbert spaces are isomorphic, so state spaces of any two quantum systems are isomorphic. Does it mean that all quantum systems are equivalent? No, it only shows that mathematically isomorphic structures can play physically different roles. Similarly, the arithmetic given by (2)–(5) is isomorphic to the standard arithmetic of $\mathbb{R}$, but it does not imply that the two arithmetics are physically equivalent.

Example 2.

The origin of Einstein’s special theory of relativity goes back to the observation that the velocity of a source of light does not influence the velocity of light itself, contradicting our everyday experiences with velocities in trains or football. Relativistic addition of velocities is based on a fundamental unit c and the dimensionless parameter β, related to velocity by $v=\beta c$. $\beta \in \mathbb{X}=(-1,1)$ while the bijection reads $f_{\mathbb{X}}\left(\beta \right)=arctanh\beta$. The velocities are added or subtracted by means of (2) and (3),

$$\begin{array}{c}\hfill {\beta }_1{\oplus }_{\mathbb{X}}{\beta }_2=tanh(arctanh{\beta }_1+arctanh{\beta }_2).\end{array}$$

(6)

Interestingly, (4) and (5) are not directly employed in special relativity. The presence of the fundamental unit c is a signature of a general non-Diophantine arithmetic (which typically works with dimensionless numbers). Numbers $\pm 1\in \mathbb{R}$ play the roles of infinities, $\pm 1_{\mathbb{R}}=\pm {\infty }_{\mathbb{X}}$. The velocity of light is therefore literally infinite in the non-Diophantine sense. The neutral element of multiplication, $1_{\mathbb{X}}=f_{\mathbb{X}}^{-1}\left(1\right)=tanh1=0.76$ (i.e., $v=0.76c$), does not seem to play in relativistic physics any privileged role.

Sometimes, for example in the context of Bell’s theorem, one works with mixed arithmetics of the form [13]

$$\begin{array}{c}\hfill x_1{\odot }_{\mathbb{Z}}^{\mathbb{XY}}{y}_2=f_{\mathbb{Z}}^{-1}\left(f_{\mathbb{X}}\left(x_1\right)·f_{\mathbb{Y}}\left(y_2\right)\right),{\odot }_{\mathbb{Z}}^{\mathbb{XY}}:\mathbb{X}\times \mathbb{Y}\to \mathbb{Z},\mathrm{etc}.\end{array}$$

(7)

Mixed arithmetics naturally occur in Taylor expansions of functions whose domains and codomains involve different arithmetics, and in the chain rule for derivatives (see Example 6).

In order to define calculus one needs limits “to zero”, and thus the notion of zero itself. In the arithmetic context a zero is a neutral element of addition, for example, $x{\oplus }_{\mathbb{X}}{0}_{\mathbb{X}}=x$ for any $x\in \mathbb{X}$. Obviously, such a zero is arithmetic-dependent. The same concerns a “one”, a neutral element of multiplication, fulfilling $x{\odot }_{\mathbb{X}}{1}_{\mathbb{X}}=x$ for any $x\in \mathbb{X}$. Once the arithmetic in $\mathbb{X}$ is specified, both neutral elements are uniquely given by the general formula: $r_{\mathbb{X}}=f_{\mathbb{X}}^{-1}\left(r\right)$ for any $r\in \mathbb{R}$. Therefore, in particular, $0_{\mathbb{X}}=f_{\mathbb{X}}^{-1}\left(0\right)$, $1_{\mathbb{X}}=f_{\mathbb{X}}^{-1}\left(1\right)$. One easily verifies that

$$\begin{array}{ccc}\hfill r_{\mathbb{X}}\oplus s_{\mathbb{X}}& =& {(r+s)}_{\mathbb{X}},\hfill \end{array}$$

(8)

$$\begin{array}{ccc}\hfill r_{\mathbb{X}}\odot s_{\mathbb{X}}& =& {\left(rs\right)}_{\mathbb{X}},\hfill \end{array}$$

(9)

for all $r,s\in \mathbb{R}$, which extends also to mixed arithmetics,

$$\begin{array}{ccc}\hfill r_{\mathbb{X}}{\oplus }_{\mathbb{X}}^{\mathbb{XY}}{s}_{\mathbb{Y}}& =& {(r+s)}_{\mathbb{X}},\hfill \end{array}$$

(10)

$$\begin{array}{ccc}\hfill r_{\mathbb{X}}{\oplus }_{\mathbb{Y}}^{\mathbb{XY}}{s}_{\mathbb{Y}}& =& {(r+s)}_{\mathbb{Y}},\hfill \end{array}$$

(11)

$$\begin{array}{ccc}\hfill r_{\mathbb{X}}{\oplus }_{\mathbb{Z}}^{\mathbb{XY}}{s}_{\mathbb{Y}}& =& {(r+s)}_{\mathbb{Z}},\mathrm{etc}.\hfill \end{array}$$

(12)

If there is no danger of ambiguity one can simplify the notation by ${\oplus }_{\mathbb{X}}^{\mathbb{XX}}={\oplus }_{\mathbb{X}}$ or ${\oplus }_{\mathbb{X}}^{\mathbb{XY}}={\oplus }_{\mathbb{X}}^{\mathbb{Y}}$. Mixed arithmetics can be given an interpretation in terms of communication channels. Mixed multiplication is in many respects analogous to a tensor product [13].

Example 3.

Consider $\mathbb{X}={\mathbb{R}}_{+}$, $\mathbb{Y}=-{\mathbb{R}}_{+}$, $f_{\mathbb{X}}\left(x\right)=lnx$, $f_{\mathbb{X}}^{-1}\left(r\right)=e^r$, $f_{\mathbb{Y}}\left(x\right)=ln(-x)$, $f_{\mathbb{Y}}^{-1}\left(r\right)=-e^r$. “Two plus two equals four” looks here as follows,

$$\begin{array}{ccc}\hfill 2_{\mathbb{X}}{\oplus }_{\mathbb{X}}{2}_{\mathbb{X}}& =& f_{\mathbb{X}}^{-1}(2+2)=4_{\mathbb{X}}=e^4,\hfill \end{array}$$

(13)

$$\begin{array}{ccc}\hfill 2_{\mathbb{X}}{\oplus }_{\mathbb{X}}^{\mathbb{Y}}{2}_{\mathbb{Y}}& =& f_{\mathbb{X}}^{-1}(2+2)=4_{\mathbb{X}}=e^4,\hfill \end{array}$$

(14)

$$\begin{array}{ccc}\hfill 2_{\mathbb{Y}}{\oplus }_{\mathbb{Y}}{2}_{\mathbb{Y}}& =& f_{\mathbb{Y}}^{-1}(2+2)=4_{\mathbb{Y}}=-e^4,\hfill \end{array}$$

(15)

$$\begin{array}{ccc}\hfill 2_{\mathbb{X}}{\oplus }_{\mathbb{Y}}^{\mathbb{X}}{2}_{\mathbb{Y}}& =& f_{\mathbb{Y}}^{-1}(2+2)=4_{\mathbb{Y}}=-e^4,\hfill \end{array}$$

(16)

where $2_{\mathbb{X}}=f_{\mathbb{X}}^{-1}\left(2\right)=e^2$, $2_{\mathbb{Y}}=f_{\mathbb{Y}}^{-1}\left(2\right)=-e^2$. From the point of view of communication channels the situation is as follows. There are two parties (“Alice” and “Bob”), each computing by means of her/his own rules. They communicate their results and agree the numbers they have found are the same, namely, “two” and “four”. However, for an external observer (an eavesdropper “Eve”), their results are opposite, say $e^4$ and $-e^4$. Mixed arithmetic plays a role of a “connection” relating different local arithmetics. This is why, in the terminology of Burgin, these types or arithmetics are non-Diophantine (from Diophantus of Alexandria who formalized the standard arithmetic). Similarly to nontrivial manifolds, non-Diophantine arithmetics do not have to admit a single global description (which we nevertheless assume in this paper).

A limit such as ${lim}_{x^{\prime }\to x}A\left(x^{\prime }\right)=A\left(x\right)$ is defined by the diagram (1) as follows,

$$\begin{array}{c}\hfill \underset{x^{\prime }\to x}{lim}A\left(x^{\prime }\right)=f_{\mathbb{Y}}^{-1}\left(\underset{r\to f_{\mathbb{X}}\left(x\right)}{lim}\tilde{A}\left(r\right)\right)\end{array}$$

(17)

i.e., in terms of an ordinary limit in $\mathbb{R}$. A non-Newtonian derivative is then defined by

$$\begin{array}{ccc}\hfill \frac{\mathrm{D}A\left(x\right)}{\mathrm{D}x}& =& \underset{\delta \to 0}{lim}\left(A\left(x{\oplus }_{\mathbb{X}}{\delta }_{\mathbb{X}}\right){\ominus }_{\mathbb{Y}}A\left(x\right)\right){\oslash }_{\mathbb{Y}}{\delta }_{\mathbb{Y}}=f_{\mathbb{Y}}^{-1}\left(\frac{\mathrm{d}\tilde{A}\left(f_{\mathbb{X}}\left(x\right)\right)}{\mathrm{d}{f}_{\mathbb{X}}\left(x\right)}\right),\hfill \end{array}$$

(18)

if the Newtonian derivative $\mathrm{d}\tilde{A}\left(r\right)/\mathrm{d}r$ exists. It is additive,

$$\begin{array}{ccc}\hfill \frac{\mathrm{D}\left[A\left(x\right){\oplus }_{\mathbb{Y}}B\left(x\right)\right]}{\mathrm{D}x}& =& \frac{\mathrm{D}A\left(x\right)}{\mathrm{D}x}{\oplus }_{\mathbb{Y}}\frac{\mathrm{D}B\left(x\right)}{\mathrm{D}x},\hfill \end{array}$$

(19)

and satisfies the Leibniz rule,

$$\begin{array}{ccc}\hfill \frac{\mathrm{D}\left[A\left(x\right){\odot }_{\mathbb{Y}}B\left(x\right)\right]}{\mathrm{D}x}& =& \left(\frac{\mathrm{D}A\left(x\right)}{\mathrm{D}x}{\odot }_{\mathbb{Y}}B\left(x\right)\right){\oplus }_{\mathbb{Y}}\left(A\left(x\right){\odot }_{\mathbb{Y}}\frac{\mathrm{D}B\left(x\right)}{\mathrm{D}x}\right).\hfill \end{array}$$

(20)

A general chain rule for compositions of functions involving arbitrary arithmetics in domains and codomains can be derived [12] (see Example 6). It implies, in particular, that the bijections defining the arithmetics are themselves always non-Newtonian differentiable (with respect to the derivatives they define). The resulting derivatives are “trivial”,

$$\begin{array}{c}\hfill \frac{\mathrm{D}{f}_{\mathbb{X}}\left(x\right)}{\mathrm{D}x}=1=\frac{\mathrm{D}{f}_{\mathbb{Y}}\left(y\right)}{\mathrm{D}y},\frac{\mathrm{D}{f}_{\mathbb{X}}^{-1}\left(r\right)}{\mathrm{D}r}=1_{\mathbb{X}},\frac{\mathrm{D}{f}_{\mathbb{Y}}^{-1}\left(r\right)}{\mathrm{D}r}=1_{\mathbb{Y}}.\end{array}$$

(21)

A non-Newtonian integral is defined by the requirement that, under typical assumptions paralleling those from the fundamental theorem of Newtonian calculus, one finds

$$\begin{array}{ccc}\hfill \frac{\mathrm{D}}{\mathrm{D}x}{\int }_y^{x}A\left(x^{\prime }\right)\mathrm{D}{x}^{\prime }& =& A\left(x\right),\hfill \end{array}$$

(22)

$$\begin{array}{ccc}\hfill {\int }_y^x\frac{\mathrm{D}A\left(x^{\prime }\right)}{\mathrm{D}{x}^{\prime }}\mathrm{D}{x}^{\prime }& =& A\left(x\right){\ominus }_{\mathbb{Y}}A\left(y\right),\hfill \end{array}$$

(23)

which uniquely implies that

$$\begin{array}{ccc}\hfill {\int }_y^{x}A\left(x^{\prime }\right)\mathrm{D}{x}^{\prime }& =& f_{\mathbb{Y}}^{-1}\left({\int }_{f_{\mathbb{X}}\left(y\right)}^{f_{\mathbb{X}}\left(x\right)}\tilde{A}\left(r\right)\mathrm{d}r\right).\hfill \end{array}$$

(24)

Here, as before, $\tilde{A}$ is defined by (1) and $\mathrm{d}r$ denotes the usual Newtonian (Riemann, Lebesgue, etc.) integration. To have a feel of the potential inherent in this simple formula, let us mention that for a Koch-type fractal (24) turns out to be equivalent to the Hausdorff integral [12,36,37]. In applications, typically the only nontrivial element is to find the explicit form of $f_{\mathbb{X}}$. It should be stressed that (24) reduces any integral to the one over a subset of $\mathbb{R}$. The fact that such a counterintuitive possibility exists was noticed already by Wiener in his 1933 lectures on Fourier analysis [38].

3. Non-Newtonian Exponential Function and Logarithm

Once we know how to differentiate and integrate, we can turn to differential equations. The so-called exponential family plays a crucial role in thermodynamics, both standard and generalized [39,40,41,42,43]. Many different deformations of the usual $e^x$ can be found in the literature. However, from the non-Newtonian perspective, the exponential function $\mathrm{Exp}:\mathbb{X}\to \mathbb{Y}$ is defined by

$$\begin{array}{c}\hfill \frac{\mathrm{D}\mathrm{Exp}\left(x\right)}{\mathrm{D}x}=\mathrm{Exp}\left(x\right),\mathrm{Exp}\left(0_{\mathbb{X}}\right)=1_{\mathbb{Y}}.\end{array}$$

(25)

Integrating (25) (in a non-Newtonian way) one finds the unique solution

$$\begin{array}{c}\hfill \mathrm{Exp}\left(x\right)=f_{\mathbb{Y}}^{-1}\left(e^{f_{\mathbb{X}}\left(x\right)}\right),\mathrm{Exp}\left(x_1{\oplus }_{\mathbb{X}}{x}_2\right)=\mathrm{Exp}\left(x_1\right){\odot }_{\mathbb{Y}}\mathrm{Exp}\left(x_2\right).\end{array}$$

(26)

In thermodynamic applications, one often encounters exponents of negative arguments, $e^{-x}$. In a non-Newtonian context the correct form of a minus is ${\ominus }_{\mathbb{X}}x=0_{\mathbb{X}}{\ominus }_{\mathbb{X}}x=f_{\mathbb{X}}^{-1}\left(-f_{\mathbb{X}}\left(x\right)\right)$. The example discussed in the next section will involve $\mathbb{X}=\mathbb{R}$ and $f_{\mathbb{X}}^{-1}(-r)=-f_{\mathbb{X}}^{-1}$(r). In consequence, it will be correct to write ${\ominus }_{\mathbb{X}}x=-x$, but in general such a simple rule may be meaningless (because “−”, as opposed to ${\ominus }_{\mathbb{X}}$, may be undefined in $\mathbb{X}$).

Example 4.

Let $\mathbb{X}=({\mathbb{R}}_{+},\oplus ,\odot )$, with the arithmetic defined by $f_{\mathbb{X}}:{\mathbb{R}}_{+}\to \mathbb{R}$, $f_{\mathbb{X}}\left(x\right)=lnx$, $f_{\mathbb{X}}^{-1}\left(r\right)=e^r$. Then

$$\begin{array}{c}\hfill {\ominus }_{\mathbb{X}}x=f_{\mathbb{X}}^{-1}\left(-f_{\mathbb{X}}\left(x\right)\right)=e^{-lnx}=1/x\in {\mathbb{R}}_{+}.\end{array}$$

(27)

The same number can be both positive and negative, depending on the arithmetic.

A (natural) logarithm is the inverse of $\mathrm{Exp}$, namely, $\mathrm{Ln}:\mathbb{Y}\to \mathbb{X}$,

$$\begin{array}{c}\hfill \mathrm{Ln}\left(y\right)=f_{\mathbb{X}}^{-1}\left(ln{f}_{\mathbb{Y}}\left(x\right)\right),\mathrm{Ln}\left(y_1{\odot }_{\mathbb{Y}}{y}_2\right)=\mathrm{Ln}\left(y_1\right){\oplus }_{\mathbb{X}}\mathrm{Ln}\left(y_2\right).\end{array}$$

(28)

Expressions such as $\mathrm{Exp}x+\mathrm{Ln}y$ are in general meaningless even if $\mathbb{X}\subset {\mathbb{R}}_{+}$ and $\mathbb{Y}\subset {\mathbb{R}}_{+}$. However, formulas such as

$$\begin{array}{c}\hfill \left(\mathrm{Exp}x\right){\oplus }_{\mathbb{Z}}^{\mathbb{Y}\mathbb{X}}\left(\mathrm{Ln}y\right)=f_{\mathbb{Z}}^{-1}\left(e^{f_{\mathbb{X}}\left(x\right)}+ln{f}_{\mathbb{Y}}\left(y\right)\right)\end{array}$$

(29)

make perfect sense. For example, if $p_k\in \mathbb{X}$, then an entropy can be defined as

$$\begin{array}{ccc}\hfill S& =& \underset{k}{⨁}{}_{\mathbb{Z}}{p}_k{\odot }_{\mathbb{Z}}^{\mathbb{X}\mathbb{Y}}\mathrm{Ln}\left(1_{\mathbb{X}}{\oslash }_{\mathbb{X}}{p}_k\right)\hfill \end{array}$$

(30)

$$\begin{array}{ccc}& =& f_{\mathbb{Z}}^{-1}\left[\sum _k{f}_{\mathbb{Z}}\left(p_k{\odot }_{\mathbb{Z}}^{\mathbb{X}\mathbb{Y}}\mathrm{Ln}\left(1_{\mathbb{X}}{\oslash }_{\mathbb{X}}{p}_k\right)\right)\right]\hfill \end{array}$$

(31)

$$\begin{array}{ccc}& =& f_{\mathbb{Z}}^{-1}\left[\sum _k{f}_{\mathbb{X}}\left(p_k\right)ln\left(1/f_{\mathbb{X}}\left(p_k\right)\right)\right].\hfill \end{array}$$

(32)

Many intriguing questions occur if one asks about normalization of probabilities. We will come to it later.

Non-Newtonian constructions of Exp and Ln are systematic, general, and flexible. There seems to exist a relation between the arithmetic formalism and the method of monotone embedding discussed in information geometry [44], but the problem requires further studies.

Example 5.

In order to appreciate the difference between Newtonian and non-Newtonian differentiation let us differentiate the function $A\left(x\right)=x$, $A:\mathbb{X}\to \mathbb{Y}$, but in two cases. The first one is trivial, $\mathbb{X}=\mathbb{Y}=(\mathbb{R},+,·)$, with the arithmetic defined by the identity $f_{\mathbb{X}}=f_{\mathbb{Y}}={\mathrm{id}}_{\mathbb{R}}$. Then, the non-Newtonian and Newtonian derivatives coincide, so

$$\begin{array}{c}\hfill \frac{\mathrm{D}A\left(x\right)}{\mathrm{D}x}=\frac{\mathrm{d}A\left(x\right)}{\mathrm{d}x}=1.\end{array}$$

(33)

The second case involves, as before, the codomain $\mathbb{Y}=(\mathbb{R},+,·)$, with the arithmetic defined by the identity $f_{\mathbb{Y}}={\mathrm{id}}_{\mathbb{R}}$. However, as the domain we choose $\mathbb{X}=({\mathbb{R}}_{+},\oplus ,\odot )$, with the arithmetic defined by $f_{\mathbb{X}}:{\mathbb{R}}_{+}\to \mathbb{R}$, $f_{\mathbb{X}}\left(x\right)=lnx$, $f_{\mathbb{X}}^{-1}\left(r\right)=e^r$. Now,

$$\begin{array}{ccc}\hfill \frac{\mathrm{D}A\left(x\right)}{\mathrm{D}x}& =& \underset{\delta \to 0}{lim}\left(A\left(x{\oplus }_{\mathbb{X}}{\delta }_{\mathbb{X}}\right){\ominus }_{\mathbb{Y}}A\left(x\right)\right){\oslash }_{\mathbb{Y}}{\delta }_{\mathbb{Y}}=\underset{\delta \to 0}{lim}\frac{\left(x{\oplus }_{\mathbb{X}}{f}_{\mathbb{X}}^{-1}\left(\delta \right)\right)-x}{\delta }\hfill \\ & =& \underset{\delta \to 0}{lim}\frac{e^{lnx+\delta }-x}{\delta }=x=A\left(x\right).\hfill \end{array}$$

(34)

As, $0_{\mathbb{X}}=f_{\mathbb{X}}^{-1}\left(0\right)=e^0=1$, we find $A\left(0_{\mathbb{X}}\right)=0_{\mathbb{X}}=1=1_{\mathbb{Y}}$, and conclude that $A\left(x\right)=x$, $A:{\mathbb{R}}_{+}\to \mathbb{R}$ belongs to the exponential family. Indeed,

$$\begin{array}{c}\hfill A\left(x_1{\oplus }_{\mathbb{X}}{x}_2\right)=x_1{\oplus }_{\mathbb{X}}{x}_2=e^{ln{x}_1+ln{x}_2}=x_1·x_2=A\left(x_1\right){\odot }_{\mathbb{Y}}A\left(x_2\right).\end{array}$$

(35)

To understand the result, write $A\left(x\right)=f_{\mathbb{Y}}^{-1}\left(\tilde{A}(f_{\mathbb{X}}\left(x\right)\right)=\tilde{A}(lnx)=x$, so that $\tilde{A}\left(r\right)=e^r$. Then, by the second form of derivative in (18),

$$\begin{array}{ccc}\hfill \frac{\mathrm{D}A\left(x\right)}{\mathrm{D}x}& =& f_{\mathbb{Y}}^{-1}\left(\frac{\mathrm{d}\tilde{A}\left(f_{\mathbb{X}}\left(x\right)\right)}{\mathrm{d}{f}_{\mathbb{X}}\left(x\right)}\right)=\frac{\mathrm{d}{e}^{f_{\mathbb{X}}\left(x\right)}}{\mathrm{d}{f}_{\mathbb{X}}\left(x\right)}=e^{f_{\mathbb{X}}\left(x\right)}=e^{lnx}=x.\hfill \end{array}$$

(36)

The map A does not affect the value of x, but changes its arithmetic properties. It behaves as if it assigned a different meaning to the same word. The example becomes even more intriguing if one realizes that logarithm is known to approximately relate stimulus with sensation in real-life sensory systems (hence the logarithmic scale of decibels and star magnitudes) [35].

Example 6.

Many calculations in thermodynamics reduce to formulas of the form

$$\begin{array}{c}\hfill \mathrm{d}U(S,V)={\left(\frac{\partial U}{\partial S}\right)}_V\mathrm{d}S+{\left(\frac{\partial U}{\partial V}\right)}_S\mathrm{d}V,\end{array}$$

(37)

being equivalent to the derivative $\mathrm{d}U\left(S\right(t),V(t\left)\right)/\mathrm{d}t$ of a composite function of several variables. The latter has a unique formulation in non-Newtonian calculus: One only needs to specify the arithmetics. For example, let U be a map $U:\mathbb{S}\times \mathbb{V}\to \mathbb{U}$, and let $S:\mathbb{T}\to \mathbb{S}$, $V:\mathbb{T}\to \mathbb{V}$. Then,

$$\begin{array}{ccc}\hfill \frac{\mathrm{D}U\left(S\right(t),V(t\left)\right)}{\mathrm{D}t}& =& \underset{\delta \to 0}{lim}\left(U\left((S\left(t{\oplus }_{\mathbb{T}}{\delta }_{\mathbb{T}}\right),V\left(t{\oplus }_{\mathbb{T}}{\delta }_{\mathbb{T}}\right)\right){\ominus }_{\mathbb{U}}U\left(S\left(t\right),V\left(t\right)\right)\right){\oslash }_{\mathbb{U}}{\delta }_{\mathbb{U}}.\hfill \end{array}$$

(38)

As

$$\begin{array}{c}\hfill \underset{x^{\prime }\to x}{lim}A\left(x\right){\oplus }_{\mathbb{Y}}B\left(x\right)=\left(\underset{x^{\prime }\to x}{lim}A\left(x\right)\right){\oplus }_{\mathbb{Y}}\left(\underset{x^{\prime }\to x}{lim}B\left(x\right)\right)\end{array}$$

(39)

(see Appendix A), we rewrite (38) as

$$\begin{array}{ccc}\hfill \frac{\mathrm{D}U\left(S\right(t),V(t\left)\right)}{\mathrm{D}t}& =& \underset{\delta \to 0}{lim}\left(U\left((S\left(t{\oplus }_{\mathbb{T}}{\delta }_{\mathbb{T}}\right),V\left(t{\oplus }_{\mathbb{T}}{\delta }_{\mathbb{T}}\right)\right){\ominus }_{\mathbb{U}}U\left((S\left(t{\oplus }_{\mathbb{T}}{\delta }_{\mathbb{T}}\right),V\left(t\right)\right)\right){\oslash }_{\mathbb{U}}{\delta }_{\mathbb{U}}\hfill \\ & & {\oplus }_{\mathbb{U}}\underset{\delta \to 0}{lim}\left(U\left((S\left(t{\oplus }_{\mathbb{T}}{\delta }_{\mathbb{T}}\right),V\left(t\right)\right){\ominus }_{\mathbb{U}}U\left(S\left(t\right),V\left(t\right)\right)\right){\oslash }_{\mathbb{U}}{\delta }_{\mathbb{U}}.\hfill \end{array}$$

(40)

Under the usual assumptions about continuity of $\tilde{U}:\mathbb{R}\times \mathbb{R}\to \mathbb{R}$ in

$$\begin{array}{c}\hfill \begin{array}{ccc}\hfill \mathbb{S}\times \mathbb{V}& \stackrel{U}{⟶}& \mathbb{U}\hfill \\ \hfill f_{\mathbb{S}}\downarrow f_{\mathbb{V}}\downarrow & & \downarrow f_{\mathbb{U}}\hfill \\ \hfill \mathbb{R}\times \mathbb{R}& \stackrel{\tilde{U}}{⟶}& \mathbb{R}\hfill \end{array},\end{array}$$

(41)

we reduce (40) to

$$\begin{array}{ccc}\hfill \frac{\mathrm{D}U\left(S\right(t),V(t\left)\right)}{\mathrm{D}t}& =& \underset{\delta \to 0}{lim}\left(U\left((S\left(t\right),V\left(t{\oplus }_{\mathbb{T}}{\delta }_{\mathbb{T}}\right)\right){\ominus }_{\mathbb{U}}U\left(\left(S\right(t),V(t)\right)\right){\oslash }_{\mathbb{U}}{\delta }_{\mathbb{U}}\hfill \\ & & {\oplus }_{\mathbb{U}}\underset{\delta \to 0}{lim}\left(U\left((S\left(t{\oplus }_{\mathbb{T}}{\delta }_{\mathbb{T}}\right),V\left(t\right)\right){\ominus }_{\mathbb{U}}U\left(S\left(t\right),V\left(t\right)\right)\right){\oslash }_{\mathbb{U}}{\delta }_{\mathbb{U}},\hfill \end{array}$$

(42)

and then to two instances of the non-Newtonian chain rule,

$$\begin{array}{c}\hfill \frac{\mathrm{D}(B\circ A)\left(x\right)}{\mathrm{D}x}=f_{\mathbb{Z}}^{-1}\left[f_{\mathbb{Z}}\left(\frac{\mathrm{D}B\left(A\left(x\right)\right)}{DA\left(x\right)}\right)f_{\mathbb{Y}}\left(\frac{\mathrm{D}A\left(x\right)}{\mathrm{D}x}\right)\right]=\frac{\mathrm{D}B\left(A\left(x\right)\right)}{DA\left(x\right)}{\odot }_{\mathbb{Z}}^{\mathbb{ZY}}\frac{\mathrm{D}A\left(x\right)}{\mathrm{D}x},\end{array}$$

(43)

valid for the composition

$$\begin{array}{c}\hfill \begin{array}{ccccc}\hfill \mathbb{X}& \hfill \stackrel{A}{⟶}& \hfill \mathbb{Y}& \hfill \stackrel{B}{⟶}& \hfill \mathbb{Z}\\ \hfill f_{\mathbb{X}}\downarrow & & \hfill f_{\mathbb{Y}}\downarrow & & \hfill f_{\mathbb{Z}}\downarrow \\ \hfill \mathbb{R}& \hfill \stackrel{\tilde{A}}{⟶}& \hfill \mathbb{R}& \hfill \stackrel{\tilde{B}}{⟶}& \hfill \mathbb{R}\end{array}\end{array}$$

(44)

of maps. Finally,

$$\begin{array}{c}\hfill \frac{\mathrm{D}U\left(S\right(t),V(t\left)\right)}{\mathrm{D}t}=\frac{\mathrm{D}U\left(S\right(t),V(t\left)\right)}{\mathrm{D}S\left(t\right)}{\odot }_{\mathbb{U}}^{\mathbb{US}}\frac{\mathrm{D}S\left(t\right)}{\mathrm{D}t}{\oplus }_{\mathbb{U}}\frac{\mathrm{D}U\left(S\right(t),V(t\left)\right)}{\mathrm{D}V\left(t\right)}{\odot }_{\mathbb{U}}^{\mathbb{UV}}\frac{\mathrm{D}V\left(t\right)}{\mathrm{D}t}.\end{array}$$

(45)

Effectively,

$$\begin{array}{c}\hfill \mathrm{D}U(S,V)={\left(\frac{\mathrm{D}U}{\mathrm{D}S}\right)}_V{\odot }_{\mathbb{U}}^{\mathbb{US}}\mathrm{D}S{\oplus }_{\mathbb{U}}{\left(\frac{\mathrm{D}U}{\mathrm{D}V}\right)}_S{\odot }_{\mathbb{U}}^{\mathbb{UV}}\mathrm{D}V,\end{array}$$

(46)

is the non-Newtonian formula for a differential.

The next section shows that the above mentioned subtleties with arithmetics of domains and codomains have straightforward implications for generalized thermostatistics.

4. Kaniadakis $\mathbf{\kappa }$-Calculus Versus Non-Newtonian Calculus

Kaniadakis, in a series of papers [26,27,28,29,30,31,32,33,34], developed a generalized form of arithmetic and calculus, with numerous applications to statistical physics, and beyond. In the present section, we will clarify links between his formalism and non-Newtonian calculus. As we will see, some of the results have a straightforward non-Newtonian interpretation, but not all.

Assume $\mathbb{X}=\mathbb{R}$, with the bijection $f_{\mathbb{X}}\equiv f_{\kappa }:\mathbb{R}\to \mathbb{R}$ given explicitly by

$$\begin{array}{ccc}\hfill f_{\kappa }\left(x\right)& =& \frac{1}{\kappa }\mathrm{arcsinh}\kappa x,\hfill \end{array}$$

(47)

$$\begin{array}{ccc}\hfill f_{\kappa }^{-1}\left(x\right)& =& \frac{1}{\kappa }sinh\kappa x.\hfill \end{array}$$

(48)

Kaniadakis’ $\kappa$-calculus begins with the arithmetic,

$$\begin{array}{ccc}\hfill x\stackrel{\kappa }{\oplus }y& =& f_{\kappa }^{-1}\left(f_{\kappa }\left(x\right)+f_{\kappa }\left(y\right)\right),\hfill \end{array}$$

(49)

$$\begin{array}{ccc}\hfill x\stackrel{\kappa }{\ominus }y& =& f_{\kappa }^{-1}\left(f_{\kappa }\left(x\right)-f_{\kappa }\left(y\right)\right),\hfill \end{array}$$

(50)

$$\begin{array}{ccc}\hfill x\stackrel{\kappa }{\odot }y& =& f_{\kappa }^{-1}\left(f_{\kappa }\left(x\right)·f_{\kappa }\left(y\right)\right),\hfill \end{array}$$

(51)

$$\begin{array}{ccc}\hfill x\stackrel{\kappa }{\oslash }y& =& f_{\kappa }^{-1}\left(f_{\kappa }\left(x\right)/f_{\kappa }\left(y\right)\right).\hfill \end{array}$$

(52)

As $f_0\left(x\right)=x$, the case $\kappa =0$ corresponds to the usual field ${\mathbb{R}}_0=(\mathbb{R},+,·)$, which we will shortly denote by $\mathbb{R}$. The neutral element of addition, $0_{\kappa }=f_{\kappa }^{-1}\left(0\right)=0$, is the same for all $\kappa$s. The neutral element of $\kappa$-multiplication is nontrivial, $1_{\kappa }=f_{\kappa }^{-1}\left(1\right)\ne 1$. The fields ${\mathbb{R}}_{\kappa }=(\mathbb{R},\stackrel{\kappa }{\oplus },\stackrel{\kappa }{\odot })$ are isomorphic to one another due to their isomorphism with ${\mathbb{R}}_0$,

$$\begin{array}{ccc}\hfill f_{\kappa }\left(x\stackrel{\kappa }{\oplus }y\right)& =& f_{\kappa }\left(x\right)+f_{\kappa }\left(y\right),\hfill \end{array}$$

(53)

$$\begin{array}{ccc}\hfill f_{\kappa }\left(x\stackrel{\kappa }{\odot }y\right)& =& f_{\kappa }\left(x\right)·f_{\kappa }\left(y\right).\hfill \end{array}$$

(54)

Kaniadakis defines his $\kappa$-derivative of a real function $A\left(x\right)$ as

$$\begin{array}{ccc}\hfill \frac{\mathrm{d}A\left(x\right)}{{\mathrm{d}}_{\kappa }x}& =& \underset{\delta \to 0}{lim}\frac{A(x+\delta )-A\left(x\right)}{(x+\delta )\stackrel{\kappa }{\ominus }x}=\frac{\mathrm{d}A\left(x\right)}{\mathrm{d}x}/\frac{\mathrm{d}{f}_{\kappa }\left(x\right)}{\mathrm{d}x}=\frac{\mathrm{d}A\left(x\right)}{\mathrm{d}x}\sqrt{1+{\kappa }^2{x}^2}.\hfill \end{array}$$

(55)

We will now specify in which sense the $\kappa$-derivative is non-Newtonian. First consider a function A,

$$\begin{array}{c}\hfill \begin{array}{ccc}\hfill {\mathbb{R}}_{{\kappa }_1}& \stackrel{A}{⟶}& {\mathbb{R}}_{{\kappa }_2}\hfill \\ \hfill f_{{\kappa }_1}\downarrow & & \downarrow f_{{\kappa }_2}\hfill \\ \hfill \mathbb{R}& \stackrel{\tilde{A}}{⟶}& \mathbb{R}\hfill \end{array}\end{array}$$

(56)

Its non-Newtonian derivative

$$\begin{array}{c}\hfill \frac{\mathrm{D}A\left(x\right)}{\mathrm{D}x}=\underset{\delta \to 0}{lim}\left(A(x\stackrel{{\kappa }_1}{\oplus }{\delta }_{{\kappa }_1})\stackrel{{\kappa }_2}{\ominus }A\left(x\right)\right)\stackrel{{\kappa }_2}{\oslash }{\delta }_{{\kappa }_2},\end{array}$$

(57)

if compared with (55), suggests ${\kappa }_2=0$. Setting ${\kappa }_1=\kappa$, ${\kappa }_2=0$, we find

$$\begin{array}{c}\hfill \frac{\mathrm{D}A\left(x\right)}{\mathrm{D}x}=\underset{\delta \to 0}{lim}\frac{A(x\stackrel{\kappa }{\oplus }{\delta }_{\kappa })-A\left(x\right)}{\delta }=\underset{\delta \to 0}{lim}\frac{A[x\stackrel{\kappa }{\oplus }{f}_{\kappa }^{-1}\left(\delta \right)]-A\left(x\right)}{\delta }=\underset{\delta \to 0}{lim}\frac{A(x\stackrel{\kappa }{\oplus }\delta )-A\left(x\right)}{\delta },\end{array}$$

(58)

as $f_{\kappa }^{-1}\left(\delta \right)\approx \delta$ for $\delta \approx 0$. Denoting $x\stackrel{\kappa }{\oplus }\delta =x+{\delta }^{\prime }$ we find $\delta =(x+{\delta }^{\prime })\stackrel{\kappa }{\ominus }x$, and

$$\begin{array}{c}\hfill \frac{\mathrm{D}A\left(x\right)}{\mathrm{D}x}=\underset{{\delta }^{\prime }\to 0}{lim}\frac{A(x+{\delta }^{\prime })-A\left(x\right)}{(x+{\delta }^{\prime })\stackrel{\kappa }{\ominus }x},\end{array}$$

(59)

in agreement with the Kaniadakis formula. However, as a by-product of the calculation we have proved that $\kappa$-calculus is applicable only to functions mapping ${\mathbb{R}}_{\kappa }$ into $\mathbb{R}$. Kaniadakis exponential function satisfies

$$\begin{array}{c}\hfill \frac{\mathrm{D}\mathrm{Exp}\left(x\right)}{\mathrm{D}x}=\mathrm{Exp}\left(x\right),\mathrm{Exp}\left(0\right)=1,\end{array}$$

(60)

with $0=0_{\kappa }$, $1=1_0$. Accordingly,

$$\begin{array}{c}\hfill \mathrm{Exp}\left(x\right)=f_{\mathbb{Y}}^{-1}\left(e^{f_{\mathbb{X}}\left(x\right)}\right)=e^{f_{\kappa }\left(x\right)}=e^{\frac{1}{\kappa }\mathrm{arcsinh}\kappa x},\end{array}$$

(61)

which is indeed the Kaniadakis result. Recalling that $f_{\mathbb{Y}}\left(x\right)=x$, we find the explicit form of the logarithm, $\mathrm{Ln}:\mathbb{R}\to {\mathbb{R}}_{\kappa }$,

$$\begin{array}{c}\hfill \mathrm{Ln}\left(y\right)=f_{\mathbb{X}}^{-1}\left(ln{f}_{\mathbb{Y}}\left(y\right)\right)=\frac{1}{\kappa }sinh(\kappa lny),\end{array}$$

(62)

which again agrees with the Kaniadakis definition.

Yet, the readers must be hereby warned that it is not allowed to apply the Kaniadakis definition of derivative to $\mathrm{Ln}x$. The correct non-Newtonian form is

$$\begin{array}{c}\hfill \frac{\mathrm{D}\mathrm{Ln}\left(y\right)}{\mathrm{D}y}=\underset{\delta \to 0}{lim}\left(\mathrm{Ln}(y+\delta )\stackrel{\kappa }{\ominus }\mathrm{Ln}\left(y\right)\right)\stackrel{\kappa }{\oslash }{\delta }_{\kappa }=f_{\mathbb{X}}^{-1}\left(1/f_{\mathbb{Y}}\left(y\right)\right)=\frac{1}{\kappa }sinh(\kappa /y),\end{array}$$

(63)

because $\mathrm{Ln}$ maps $\mathbb{R}$ into ${\mathbb{R}}_{\kappa }$. Kaniadakis is aware of the subtlety and thus introduces also another derivative, meant for differentiation of inverse functions,

$$\begin{array}{c}\hfill \frac{{\mathrm{d}}_{\kappa }A\left(y\right)}{\mathrm{d}y}=\underset{u\to y}{lim}\frac{A\left(y\right)\stackrel{\kappa }{\ominus }A\left(u\right)}{y-u}=\underset{\delta \to 0}{lim}\frac{A(y+\delta )\stackrel{\kappa }{\ominus }A\left(y\right)}{\delta },\end{array}$$

(64)

a definition which, from the non-Newtonian standpoint, must be nevertheless regarded as incorrect (‘/’ should be replaced by $\stackrel{\kappa }{\oslash }$ typical of the codomain ${\mathbb{R}}_{\kappa }$). As a result,

$$\begin{array}{c}\hfill \frac{{\mathrm{d}}_{\kappa }\mathrm{Ln}\left(y\right)}{\mathrm{d}y}=\frac{1}{y}\ne \frac{\mathrm{D}\mathrm{Ln}\left(y\right)}{\mathrm{D}y}=\frac{1}{\kappa }sinh\frac{\kappa }{y}.\end{array}$$

(65)

This is probably why (64), as opposed to (55), has not found too many applications.

Let us finally check what would have happened if instead of (61) one considered the exponential function mapping ${\mathbb{R}}_{\kappa }$ into itself, $f_{\mathbb{Y}}=f_{\mathbb{X}}=f_{\kappa }$,

$$\begin{array}{c}\hfill \mathrm{Exp}\left(x\right)=f_{\mathbb{Y}}^{-1}\left(e^{f_{\mathbb{X}}\left(x\right)}\right)=f_{\kappa }^{-1}\left(e^{f_{\kappa }\left(x\right)}\right)=\frac{1}{\kappa }sinh\left(\kappa e^{\frac{1}{\kappa }\mathrm{arcsinh}\kappa x}\right).\end{array}$$

(66)

As in thermodynamic applications one typically encounters $\mathrm{Exp}$ of a negative argument, one expects that physical differences between $\mathrm{Exp}:{\mathbb{R}}_{\kappa }\to {\mathbb{R}}_{\kappa }$ and $\mathrm{Exp}:{\mathbb{R}}_{\kappa }\to \mathbb{R}$ should not be essential. Moreover, indeed, Figure 1 shows that both exponents lead to identical asymptotic tails.

5. A Cosmological Aspect of the Kaniadakis Arithmetic

Kaniadakis explored possible relativistic implications of his formalism. In particular, he noted that fluxes of cosmic rays depend on energy in a way that seems to indicate $\kappa >0$. It is therefore intriguing that essentially the same arithmetic was recently shown [14] to have links with the problem of accelerated expansion of the Universe, one of the greatest puzzles of contemporary physics.

Cosmological expansion is well described by the Friedman equation,

$$\begin{array}{c}\hfill \frac{\mathrm{d}a\left(t\right)}{\mathrm{d}t}=\sqrt{{\mathsf{\Omega }}_{\mathsf{\Lambda }}a{\left(t\right)}^2+\frac{{\mathsf{\Omega }}_M}{a\left(t\right)}},a\left(t\right)>0,\end{array}$$

(67)

for a dimensionless scale factor $a\left(t\right)$ evolving in a dimensionless time t (in units of the Hubble time $t_H\approx 13.58\times {10}^9$ yr). The observable parameters are ${\mathsf{\Omega }}_M=0.3$, ${\mathsf{\Omega }}_{\mathsf{\Lambda }}=0.7$ [45,46]. ${\mathsf{\Omega }}_{\mathsf{\Lambda }}\ne 0$ is typically interpreted as an indication of dark energy. Equation (67) is solved by

$$\begin{array}{c}\hfill a\left(t\right)={\left(\sqrt{\frac{{\mathsf{\Omega }}_M}{{\mathsf{\Omega }}_{\mathsf{\Lambda }}}}sinh\frac{3\sqrt{{\mathsf{\Omega }}_{\mathsf{\Lambda }}}t}{2}\right)}^{2/3},t>0.\end{array}$$

(68)

Now assume that

$$\begin{array}{c}\hfill \begin{array}{ccc}\hfill \mathbb{X}& \stackrel{a}{⟶}& \mathbb{R}\hfill \\ \hfill f_{\mathbb{X}}\downarrow & & \downarrow f_{\mathbb{R}}={\mathrm{id}}_{\mathbb{R}}\hfill \\ \hfill \mathbb{R}& \stackrel{\tilde{a}}{⟶}& \mathbb{R}\hfill \end{array},\end{array}$$

(69)

whereas the Friedman equation involves no ${\mathsf{\Omega }}_{\mathsf{\Lambda }}$,

$$\begin{array}{c}\hfill \frac{\mathrm{D}a\left(t\right)}{\mathrm{D}t}=\sqrt{\frac{\mathsf{\Omega }}{a\left(t\right)}},a\left(t\right)>0,\end{array}$$

(70)

for some $\mathsf{\Omega }$. Its solution by non-Newtonian techniques reads

$$\begin{array}{ccc}\hfill a\left(t\right)& =& {\left(\frac{3}{2}\sqrt{\mathsf{\Omega }}{f}_{\mathbb{X}}\left(t\right)\right)}^{2/3},\hfill \end{array}$$

(71)

so, comparing (71) with (68), we find

$$\begin{array}{c}\hfill f_{\mathbb{X}}\left(t\right)=\frac{2}{3\sqrt{0.7}}\sqrt{\frac{{\mathsf{\Omega }}_M}{\mathsf{\Omega }}}sinh\frac{3\sqrt{0.7}}{2}t=\sqrt{\frac{{\mathsf{\Omega }}_M}{\mathsf{\Omega }}}{f}_{\kappa }^{-1}\left(t\right),\mathrm{for}\kappa =1.255.\end{array}$$

(72)

Accelerated expansion of the Universe looks like a combined effect of non-Euclidean geometry and non-Diophantine arithmetic. The resulting dynamics is non-Newtonian in both meanings of this term.

The presence of the inverse bijection $f_{\kappa }^{-1}$ and $\kappa >1$ raises a number of interesting questions. It is related to the fundamental duality between Diophantine and non-Diophantine arithmetics. Namely, any equation of the form, say

$$\begin{array}{ccc}\hfill x_1\oplus x_2& =& f^{-1}\left(f\left(x_1\right)+f\left(x_2\right)\right),\hfill \end{array}$$

(73)

can be inverted by $f\left(x\right)=y$ into

$$\begin{array}{ccc}\hfill y_1+y_2& =& f\left(f^{-1}\left(y_1\right)\oplus f^{-1}\left(y_2\right)\right),\hfill \end{array}$$

(74)

suggesting that it is ⊕ and not + which is the Diophantine arithmetic operation. Having two isomorphic arithmetics we, in general, do not have any criterion telling us which of the two is “normal”, and which is “generalized”.

6. Kolmogorov–Nagumo Averages and Non-Diophantine/Non-Newtonian Probability

Another non-Diophantine/non-Newtonian aspect that can be identified in the context of information theory and thermodynamics is implicitly present in the works of Kolmogorow, Nagumo, and Rényi. Let us recall that a Kolmogorov–Nagumo average is defined as [47,48,49,50,51,52,53,54]

$$\begin{array}{ccc}\hfill {\langle a\rangle }_f& =& f^{-1}\left(\sum _k{p}_{k}f\left(a_k\right)\right).\hfill \end{array}$$

(75)

Rewriting (75) as

$$\begin{array}{c}\hfill {\langle a\rangle }_f=f^{-1}\left(\sum _{k}f\left(p_k^{\prime }\right)f\left(a_k\right)\right)=\underset{k}{⨁}{p}_k^{\prime }\odot a_k,\end{array}$$

(76)

where $p_k^{\prime }=f^{-1}\left(p_k\right)$, one interprets the average as the one typical of a non-Diophantine-arithmetic-valued probability. Apparently, neither Kolmogorov nor Nagumo nor Rényi had interpreted their results from this arithmetic point of view [7].

The lack of arithmetic perspective is especially visible in the works of Rényi [49] who, while deriving his $\alpha$-entropies, began with a general Kolmogorov–Nagumo average. Trying to derive a meaningful class of fs he demanded that

$$\begin{array}{c}\hfill {\langle a+c\rangle }_f={\langle a\rangle }_f+c\end{array}$$

(77)

be valid for any constant random variable c, and this led him to the exponential family $f_{\alpha }\left(x\right)=2^{(1-\alpha )x}$ (up to a general affine transformation $f\mapsto Af+B$, which does not affect Kolmogorov–Nagumo averages). In physical applications, it is more convenient to work with natural logarithms, so let us replace $f_{\alpha }$ by $f_q\left(x\right)=e^{(1-q)x}$, $f_q^{-1}\left(x\right)=\frac{1}{1-q}lnx$, $q\in R$. With this particular choice of f one finds

$$\begin{array}{ccc}\hfill {\langle a\rangle }_{f_q}& =& \frac{1}{1-q}ln\left(\sum _k{p}_k{e}^{(1-q)a_k}\right).\hfill \end{array}$$

(78)

As is well known, the standard linear average is the limiting case ${lim}_{q\to 1}{\langle a\rangle }_{f_q}={\sum }_k{p}_k{a}_k,$ that includes the entropy of Shannon, $S={\sum }_k{p}_{k}ln(1/p_k)=S_1$, as the limit $q\to 1$ of the Rényi entropy

$$\begin{array}{c}\hfill S_q=\frac{1}{1-q}ln\left(\sum _k{p}_k{e}^{(1-q)ln(1/p_k)}\right)=\frac{1}{1-q}ln\sum _k{p}_k^q.\end{array}$$

(79)

Still, notice that ${\langle a\oplus b\rangle }_f={\langle a\rangle }_f\oplus {\langle b\rangle }_f$ for any f, so had Rényi been thinking in arithmetic categories, he would not have arrived at his $f_{\alpha }$. Yet, $f_{\alpha }$ is an interesting special case. For example,

$$\begin{array}{c}\hfill p_k^{\prime }=f_q^{-1}\left(p_k\right)=\frac{1}{q-1}ln(1/p_k).\end{array}$$

(80)

The random variable $a_k={log}_b(1/p_k)$ is, according to Shannon [49,55], the amount of information obtained by observing an event whose probability is $p_k$. The choice of b defines units of information. Therefore, Rényi’s non-Diophantine probability $p_k^{\prime }$ is the amount of information encoded in $p_k$.

7. Escort Probabilities and Quantum Mechanical Hidden Variables

Non-Diophantine arithmetics have several properties that make them analogous to sets of values of incompatible random variables in quantum mechanics. Generalized arithmetics and non-Newtonian calculi have nontrivial consequences for the problem of hidden variables and completeness of quantum mechanics.

Example 7.

Pauli matrices ${\sigma }_1$ and ${\sigma }_2$ represent random variables whose values are $s_1=\pm 1$ and $s_2=\pm 1$, respectively. However, it is not allowed to assume that ${\sigma }_1+{\sigma }_2$ represents a random variable whose possible values are $s_1+s_2=0,\pm 2$, even though an average of ${\sigma }_1+{\sigma }_2$ ia a sum of independent averages of ${\sigma }_1$ and ${\sigma }_2$. In non-Diophantine arithmetic one encounters a similar problem. In general it makes no sense to perform additions of the form $x_{\mathbb{X}}+y_{\mathbb{Y}}$ even if $x_{\mathbb{X}}\in \mathbb{R}$ and $y_{\mathbb{Y}}\in \mathbb{R}$. One should not be surprised if non-Diophantine probabilities turn out to be analogous to quantum probabilities, at least in some respects.

Normalization of probability implies

$$\begin{array}{c}\hfill 1_{\mathbb{X}}=f^{-1}\left(1\right)=f^{-1}\left(\sum _k{p}_k\right)=f^{-1}\left(\sum _{k}f\left(p_k^{\prime }\right)\right)=\underset{k}{⨁}{p}_k^{\prime }.\end{array}$$

(81)

In principle, $1_{\mathbb{X}}\ne 1$. An interesting and highly nontrivial case occurs if both $p_k$ and $p_k^{\prime }=f^{-1}\left(p_k\right)$ are probabilities in the ordinary sense, i.e., in addition to (81) one finds $1_{\mathbb{X}}=1$, $0\le p_k^{\prime }\le 1$, and ${\sum }_k{p}_k^{\prime }=1$. What can be then said about f? We can formalize the question as follows.

Problem 1.

Find a characterization of those functions $g:[0,1]\to [0,1]$ that satisfy

$$\begin{array}{c}\hfill \sum _{k}g\left(p_k\right)=1,for any choice of probabilities p_k.\end{array}$$

(82)

In analogy to the generalized thermostatistics literature we can term $p_k^{\prime }=g\left(p_k\right)$ the escort probabilities [56,57,58]. Notice that we are not in interested in the trivial solution, often employed in the context of Tsallis and Rényi entropies, where $p_k$ is replaced by $p_k^q$ and then renormalized,

$$\begin{array}{c}\hfill P_k=\frac{p_k^q}{{\sum }_j{p}_j^q}=g_k(p_1,\dots ,p_n,\dots )\end{array}$$

(83)

as $g_k(p_1,\dots ,p_n,\dots )\ne g\left(p_k\right)$ for a single function g of one variable. As we will shortly see, the solution of (82) turns out to have straightforward implications for the quantum mechanical problem of hidden variables, and relations between classical and quantum probabilities.

The most nontrivial result is found for binary probabilities, $p_1+p_2=1$.

Lemma 1.

$g\left(p_1\right)+g\left(p_2\right)=1$ for all $p_1+p_2=1$ if and only if

$$\begin{array}{c}\hfill g\left(p\right)=\frac{1}{2}+h\left(p-\frac{1}{2}\right)\end{array}$$

(84)

where $h(-x)=-h\left(x\right)$.

Proof. 

See Appendix B. □

The lemma has profound consequences for foundations of quantum mechanics, as it allows to circumvent Bell’s theorem by non-Newtonian hidden variables. For more details the readers are referred to [13,15], but here just a few examples.

Example 8.

The trivial case $g\left(p\right)=p$ implies $h\left(x\right)=x$, where $0\le p\le 1$ and $-1/2\le x\le 1/2$.

Example 9.

Consider $g\left(p\right)={sin}^2\frac{\pi }{2}p$. Then,

$$\begin{array}{c}\hfill h\left(x\right)=g\left(x+\frac{1}{2}\right)-\frac{1}{2}=\frac{1}{2}sin\pi x.\end{array}$$

(85)

Let us cross-check,

$$\begin{array}{c}\hfill g\left(p\right)+g(1-p)={sin}^2\frac{\pi }{2}p+{sin}^2\frac{\pi }{2}(1-p)={sin}^2\frac{\pi }{2}p+{cos}^2\frac{\pi }{2}p=1.\end{array}$$

(86)

Now let $p=(\pi -\theta )/\pi$ be the probability of finding a point belonging to the overlap of two half-circles rotated by θ. Then,

$$\begin{array}{c}\hfill g\left(p\right)={sin}^2\frac{\pi }{2}\frac{\pi -\theta }{\pi }={cos}^2\frac{\theta }{2}\end{array}$$

(87)

is the quantum-mechanical law describing the conditional probability for two successive measurements of spin-1/2 in two Stern–Gerlach devices placed one after another, with relative angle θ. Escort probability has become a quantum probability.

Example 10.

Let us continue the analysis of Example 9. Function $g:[0,1]\to [0,1]$, $g\left(p\right)={sin}^2\frac{\pi }{2}p$, is one-to-one. It can be continued to the bijection $g:\mathbb{R}\to \mathbb{R}$ by the periodic repetition,

$$\begin{array}{c}\hfill g\left(x\right)=n+{sin}^2\frac{\pi }{2}(x-n),n\le x\le n+1,n\in \mathbb{Z}.\end{array}$$

(88)

Now let $f=g^{-1}$. (88) leads to a non-Diophantine arithmetic and non-Newtonian calculus. Let $\theta =\alpha -\beta$, $0\le \theta \le \pi$, be an angle between two vectors representing directions of Stern-Gerlach devices. Quantum conditional probability (87) can be represented in a non-Newtonian hidden-variable form,

$$\begin{array}{ccc}\hfill {cos}^2\frac{\alpha -\beta }{2}& =& {sin}^2\frac{\pi }{2}\frac{\pi -(\alpha -\beta )}{\pi }=f^{-1}\left(\frac{1}{\pi }{\int }_{\alpha }^{\pi +\beta }\mathrm{d}r\right)=f^{-1}\left({\int }_{f\left({\alpha }^{\prime }\right)}^{f({\pi }^{\prime }\oplus {\beta }^{\prime })}\tilde{\rho }\left(r\right)\mathrm{d}r\right)\hfill \\ & =& {\int }_{{\alpha }^{\prime }}^{{\pi }^{\prime }\oplus {\beta }^{\prime }}\rho \left(\lambda \right)\mathrm{D}\lambda ,\hfill \end{array}$$

(89)

where $x^{\prime }=f^{-1}\left(x\right)$. Here, ρ is a conditional probability density of non-Newtonian hidden-variables (the half-circle is a result of conditioning by the first measurement).

Non-Newtonian calculus shifts the discussion on relations between classical and quantum probability, or classical and quantum information, into unexplored areas.

Example 11.

In typical Bell-type experiments one deals with four probabilities, corresponding to four combinations $(\pm ,\pm )$, $(\pm ,\mp )$ of pairs of binary results. The corresponding non-Newtonian model is obtained by rescaling $g\left(p_k\right)\mapsto pg(p_k/p)$, with $p=1/2$. The rescaled bijection satisfies $g\left(p_1\right)+g\left(p_2\right)=p$ for any $p_1+p_2=p$. Explicitly,

$$\begin{array}{c}\hfill g\left(p_{++}\right)+g\left(p_{+-}\right)+g\left(p_{-+}\right)+g\left(p_{--}\right)=1=p_{++}+p_{+-}+p_{-+}+p_{--}.\end{array}$$

(90)

The resulting hidden-variable model is local, but standard Bell’s inequality cannot be proved [15]. Why? Mainly because the non-Newtonian integral is not a linear map with respect to the ordinary Diophantine addition and multiplication (unless f is linear), whereas the latter is always assumed in proofs of Bell-type inequalities.

A generalization to arbitrary probabilities, $p_1+\cdots +p_n=1$, leads to an affine deformation of arithmetic, an analogue of Benioff number scaling [21,22,23,24,25]. Affine transformations do not affect Kolmogorov–Nagumo averages.

Lemma 2.

Consider probabilities $p_1,\dots ,p_n$, $n\ge 3$. $g\left(p_k\right)$ are probabilities for any choice of $p_k$ if and only if $g\left(p_k\right)=\frac{1-a+2a{p}_k}{n+(2-n)a}$, $-1\le a\le 1$.

Proof. 

See Appendix C. □

The bijection g implied by Lemma 2 depends on n. In infinitely dimensional systems, that is when n can be arbitrary, the only option is $a=1$ and thus $g\left(p\right)=p$ is the only acceptable solution. However, in spin systems there exits an alternative interpretation of this property: The dimension n grows with spin in such a way that $g_n\left(p\right)\to p$ with $n\to \infty$ is a correspondence principle meaning that very large spins are practically classical. The transition non-Diophantine → Diophantine, non-Newtonian → Newtonian becomes an analogue of non-classical → classical.

Example 12.

Limitations imposed by Lemma 2 can be nevertheless circumvented in various ways. For example, let $g\left(1\right)=1$ for a solution g from Lemma 1, so that $1_{\mathbb{X}}=1$. Obviously,

$$\begin{array}{c}\hfill 1=1_{\mathbb{X}}\odot \cdots \odot 1_{\mathbb{X}}=1\odot \cdots \odot 1=1·\dots ·1.\end{array}$$

(91)

Replacing each of the 1s by an appropriate sum of binary conditional probabilities

$$\begin{array}{c}\hfill 1=g\left(p_{k_1\dots k_{n}1}\right)+g\left(p_{k_1\dots k_{n}2}\right)=g\left(p_{k_1\dots k_{n}1}\right)\oplus g\left(p_{k_1\dots k_{n}2}\right)\end{array}$$

(92)

we can generate various conditional classical or quantum probabilities typical of a generalized Bernoulli-type process, representing several classical or quantum filters placed one after another.

8. Non-Newtonian Maximum Entropy Principle

Let us finally discuss the implications of our non-Newtonian form (32) of entropy for maximum entropy principles. Assume probabilities belong to $\mathbb{X}$. Define the Massieu function [43] by

$$\begin{array}{ccc}\hfill \mathrm{\Phi }& =& S{\ominus }_{\mathbb{Z}}{\alpha }_{\mathbb{Z}}{\odot }_{\mathbb{Z}}N{\ominus }_{\mathbb{Z}}{\beta }_{\mathbb{Z}}{\odot }_{\mathbb{Z}}H,\hfill \end{array}$$

(93)

$$\begin{array}{ccc}\hfill N& =& \underset{k}{⨁}{}_{\mathbb{Z}}^{\mathbb{X}}{p}_k=f_{\mathbb{Z}}^{-1}\left(\sum _k{f}_{\mathbb{X}}\left(p_k\right)\right),\hfill \end{array}$$

(94)

$$\begin{array}{ccc}\hfill H& =& \underset{k}{⨁}{}_{\mathbb{Z}}{p}_k{\odot }_{\mathbb{Z}}^{\mathbb{X}\mathbb{E}}{E}_k=f_{\mathbb{Z}}^{-1}\left(\sum _k{f}_{\mathbb{X}}\left(p_k\right)f_{\mathbb{E}}\left(E_k\right)\right),\hfill \end{array}$$

(95)

where $E_k\in \mathbb{E}$, and ${\alpha }_{\mathbb{Z}}=f_{\mathbb{Z}}^{-1}\left(\alpha \right)$, ${\beta }_{\mathbb{Z}}=f_{\mathbb{Z}}^{-1}\left(\beta \right)$ are Lagrange multipliers. Explicitly,

$$\begin{array}{c}\hfill \mathrm{\Phi }=f_{\mathbb{Z}}^{-1}\left[\sum _k{f}_{\mathbb{X}}\left(p_k\right)ln\left(1/f_{\mathbb{X}}\left(p_k\right)\right)-\alpha \sum _k{f}_{\mathbb{X}}\left(p_k\right)-\beta \sum _k{f}_{\mathbb{X}}\left(p_k\right)f_{\mathbb{E}}\left(E_k\right)\right].\end{array}$$

(96)

Vanishing of the derivative of $\mathrm{\Phi }$,

$$\begin{array}{c}\hfill \frac{\mathrm{D}\mathrm{\Phi }}{\mathrm{D}{p}_l}=0_{\mathbb{Z}},\end{array}$$

(97)

is equivalent to the standard formula for probabilities $f_{\mathbb{X}}\left(p_k\right)$ (see the second form of non-Newtonian derivative in (18)),

$$\begin{array}{c}\hfill \frac{\mathrm{d}}{\mathrm{d}{f}_{\mathbb{X}}\left(p_l\right)}\left(\sum _k{f}_{\mathbb{X}}\left(p_k\right)ln\left(1/f_{\mathbb{X}}\left(p_k\right)\right)-\alpha \sum _k{f}_{\mathbb{X}}\left(p_k\right)-\beta \sum _k{f}_{\mathbb{X}}\left(p_k\right)f_{\mathbb{E}}\left(E_k\right)\right)=0.\end{array}$$

(98)

Accordingly, the solution reads

$$\begin{array}{ccc}\hfill p_k& =& f_{\mathbb{X}}^{-1}\left(e^{-\beta f_{\mathbb{E}}\left(E_k\right)}/\tilde{Z}\left(\beta \right)\right)=\mathrm{Exp}\left({\ominus }_{\mathbb{E}}{\beta }_{\mathbb{E}}{\odot }_{\mathbb{E}}{E}_k\right){\oslash }_{\mathbb{X}}{Z}_{\mathbb{X}}\left(\beta \right),\hfill \end{array}$$

(99)

$$\begin{array}{ccc}\hfill Z_{\mathbb{X}}\left(\beta \right)& =& f_{\mathbb{X}}^{-1}\left(\tilde{Z}\left(\beta \right)\right)=f_{\mathbb{X}}^{-1}\left(\tilde{Z}\left(f_{\mathbb{E}}\left({\beta }_{\mathbb{E}}\right)\right)\right)=Z\left({\beta }_{\mathbb{E}}\right),\hfill \end{array}$$

(100)

and involves the exponential function $\mathrm{Exp}:\mathbb{E}\to \mathbb{X}$ we have encountered before. The normalization,

$$\begin{array}{c}\hfill 1_{\mathbb{X}}=\underset{k}{⨁}{}_{\mathbb{X}}{p}_k=f_{\mathbb{X}}^{-1}\left(\sum _k{e}^{-\beta f_{\mathbb{E}}\left(E_k\right)}/\tilde{Z}\left(\beta \right)\right)=f_{\mathbb{X}}^{-1}\left(1\right),\end{array}$$

(101)

implies the usual relation $\tilde{Z}\left(\beta \right)={\sum }_k{e}^{-\beta f_{\mathbb{E}}\left(E_k\right)}$.

Equivalently, directly at the level of $\mathbb{X}$,

$$\begin{array}{ccc}\hfill Z_{\mathbb{X}}\left(\beta \right)& =& f_{\mathbb{X}}^{-1}\left(\sum _k{e}^{-\beta f_{\mathbb{E}}\left(E_k\right)}\right)=f_{\mathbb{X}}^{-1}\left(\sum _k{f}_{\mathbb{X}}\circ f_{\mathbb{X}}^{-1}\left(e^{f_{\mathbb{E}}\left({\ominus }_{\mathbb{E}}{\beta }_{\mathbb{E}}{\odot }_{\mathbb{E}}{E}_k\right)}\right)\right)\hfill \\ & =& f_{\mathbb{X}}^{-1}\left(\sum _k{f}_{\mathbb{X}}\left(\mathrm{Exp}\left({\ominus }_{\mathbb{E}}{\beta }_{\mathbb{E}}{\odot }_{\mathbb{E}}{E}_k\right)\right)\right)=\underset{k}{⨁}{}_{\mathbb{X}}\mathrm{Exp}\left({\ominus }_{\mathbb{E}}{\beta }_{\mathbb{E}}{\odot }_{\mathbb{E}}{E}_k\right)=Z\left({\beta }_{\mathbb{E}}\right).\hfill \end{array}$$

(102)

All the standard tricks one finds in thermodynamics textbooks will work here. For example,

$$\begin{array}{ccc}\hfill H& =& f_{\mathbb{Z}}^{-1}\left(\sum _k{f}_{\mathbb{X}}\left(p_k\right)f_{\mathbb{E}}\left(E_k\right)\right)=f_{\mathbb{Z}}^{-1}\left(\sum _k{e}^{-\beta f_{\mathbb{E}}\left(E_k\right)}{f}_{\mathbb{E}}\left(E_k\right)/\tilde{Z}\left(\beta \right)\right)=f_{\mathbb{Z}}^{-1}\left(-\frac{\mathrm{d}ln\tilde{Z}\left(\beta \right)}{\mathrm{d}\beta }\right)\hfill \\ & =& {\ominus }_{\mathbb{Z}}{f}_{\mathbb{Z}}^{-1}\left(\frac{\mathrm{d}ln\tilde{Z}\left(\beta \right)}{\mathrm{d}\beta }\right)={\ominus }_{\mathbb{Z}}{f}_{\mathbb{Z}}^{-1}\left(\frac{\mathrm{d}ln\tilde{Z}\left(\beta \right)}{\mathrm{d}{f}_{\mathbb{E}}\left({\beta }_{\mathbb{E}}\right)}\right)={\ominus }_{\mathbb{Z}}{f}_{\mathbb{Z}}^{-1}\left(\frac{\mathrm{d}\tilde{A}\left(f_{\mathbb{E}}\left({\beta }_{\mathbb{E}}\right)\right)}{\mathrm{d}{f}_{\mathbb{E}}\left({\beta }_{\mathbb{E}}\right)}\right),\hfill \end{array}$$

(103)

for some function

$$\begin{array}{c}\hfill \begin{array}{ccc}\hfill \mathbb{E}& \stackrel{A}{⟶}& \mathbb{Z}\hfill \\ \hfill f_{\mathbb{E}}\downarrow & & \downarrow f_{\mathbb{Z}}\hfill \\ \hfill \mathbb{R}& \stackrel{\tilde{A}}{⟶}& \mathbb{R}\hfill \end{array}\end{array}$$

(104)

we yet have to determine. Clearly,

$$\begin{array}{ccc}\hfill ln\tilde{Z}\left(\beta \right)& =& \tilde{A}\left(f_{\mathbb{E}}\left({\beta }_{\mathbb{E}}\right)\right)=\tilde{A}\left(\beta \right),\hfill \\ \hfill A\left(x\right)& =& f_{\mathbb{Z}}^{-1}\left(\tilde{A}\left(f_{\mathbb{E}}\left(x\right)\right)\right)=f_{\mathbb{Z}}^{-1}\left(ln\left(\tilde{Z}\left(f_{\mathbb{E}}\left(x\right)\right)\right)\right)=f_{\mathbb{Z}}^{-1}\left(ln\left(f_{\mathbb{X}}\circ f_{\mathbb{X}}^{-1}\left[\tilde{Z}\left(f_{\mathbb{E}}\left(x\right)\right)\right]\right)\right)\hfill \end{array}$$

(105)

$$\begin{array}{ccc}& =& f_{\mathbb{Z}}^{-1}\left(ln{f}_{\mathbb{X}}\left(Z\left(x\right)\right)\right)=\mathrm{Ln}Z\left(x\right),\hfill \end{array}$$

(106)

where $Z:\mathbb{E}\to \mathbb{X}$, $\mathrm{Ln}:\mathbb{X}\to \mathbb{Z}$. Ultimately,

$$\begin{array}{ccc}\hfill H& =& {\ominus }_{\mathbb{Z}}\frac{\mathrm{D}(\mathrm{Ln}\circ Z)\left({\beta }_{\mathbb{E}}\right)}{\mathrm{D}{\beta }_{\mathbb{E}}}.\hfill \end{array}$$

(107)

9. Final Remarks

Non-Newtonian calculus, and the non-Diophantine arithmetics behind it, are as simple as the undergraduate arithmetic and calculus we were taught at schools. Their conceptual potential is immense but they remain largely unexplored and unappreciated. Apparently, physicists in general do not feel any need of going beyond standard Diophantine arithmetic operations, in spite of the fact that the two greatest revolutions of the 20th century physics were, in their essence, arithmetic (i.e., relativistic addition of velocities and quantum mechanical addition of probabilities). It is thus intriguing that two of the most controversial issues of modern science—dark energy and Bell’s theorem—reveal new aspects when reformulated in generalized arithmetic terms.

One should not be surprised that those who study generalizations of Boltzmann–Gibbs statistics are naturally more inclined to accept non-aprioric rules of physical arithmetic. Anyway, the very concept of non-extensivity, the core of many studies on generalized entropies, is implicitly linked with generalized forms of addition, multiplication, and differentiation [54,59,60,61].