Relativity of arithmetic as a fundamental symmetry of physics

Arithmetic operations can be defined in various ways, even if one assumes commutativity and associativity of addition and multiplication, and distributivity of multiplication with respect to addition. In consequence, whenever one encounters ‘plus’ or ‘times’ one has certain freedom of interpreting this operation. This leads to some freedom in definitions of derivatives, integrals and, thus, practically all equations occurring in natural sciences. A change of realization of arithmetic, without altering the remaining structures of a given equation, plays the same role as a symmetry transformation. An appropriate construction of arithmetic turns out to be particularly important for dynamical systems in fractal space-times. Simple examples from classical and quantum, relativistic and nonrelativistic physics are discussed, including the eigenvalue problem for a quantum harmonic oscillator. It is explained why the change of arithmetic is not equivalent to the usual change of variables.


Introduction
Symmetries of physical systems can be rather obvious or very abstract. Lorentz transformations, discovered as a formal symmetry of Maxwell's equations, seemed abstract until their physical meaning was understood by Einstein. Theory of group representations, the cornerstone of quantum mechanics and field theory, had its roots in Lie's studies of abstract symmetries of differential equations. It has taught us that differences in mathematical realizations of a symmetry may directly reflect physical differences.
Einstein's relativity, gauge invariance, Noether's theorems, Darboux-Bäcklund transformations, or supersymmetry are prominent examples of symmetry principles in physics. In the paper, we discuss a new type of principle, occurring in any physical theory: The symmetry of mathematical equations under modifications of arithmetic operations, the induced modifications of derivatives and integrals included. Similarly to other physical symmetries, the symmetry maintains the form of relevant equations, but may possess different mathematical realizations.
The approach developed in this paper is not in any sense close to those generalizations of physics that involve padic mathematics [1,2], although the context of fractals and the Cantor set might create such an impression. Another formalism that can be easily confused with what we do is based on 'arithmetic dynamics' [3] which involves the so-called linear conjugates F f = f −1 • F • f (which indeed will play a role here), but where the calculus is defined in the standard way.
Our formalism is a by-product of attempts of understanding the dynamics of fractals, but turned out to have more general consequences. In particular, it opens some new perspectives on quantum mechanics of processes occurring in sets of Lebesgue measure zero, invisible from the point of view of standard quantum mechanics, leading to phenomena analogous to dark energy, 'coming out of nowhere'.
In a wider context, what we do is closest to the works of Burgin [4] on non-Diophantine arithmetic, and Benioff on physical consequences of rescaled multiplication [5,6]. However, neither Burgin nor Benioff made the essential step of formulating a non-Diophantine calculus.

From arithmetic to calculus
To begin with, let us consider a bijection f : X → Y ⊂ R, where X is some set. The map f allows us to define addition, multiplication, subtraction, and division in X , One easily verifies the standard properties. Detailed proofs can be found in the preprint version of the paper, see arXiv:1412.8583v2 [math-ph]: (1) In fact, this is precisely an example of a non-Diophantine arithmetic in the sense of [4]. The rescaled-multiplication formalism of Benioff [5] corresponds to f (x) = px, p = 0, and in its generalized versions [6] to bundles of functions f X (x) where X is a space-time point. Our f (x) may be also regarded as an example of Benioff's value function, associating a number x with its value f (x).
Elements 0, 1 ∈ X are defined by 0 ⊕ x = x, 1 x = x, which implies f (0) = 0, f (1) = 1. One further finds x x = 0, x x = 1, as expected. Keeping the same symbols for 0, 1 ∈ X and 0, 1 ∈ R will not lead to ambiguities. However, in the approach of Benioff [5,6], where f (x) = px, p = 0, one finds f −1 (1) = 1/ p and thus 1/ p = 1 p is the unit element of multiplication. So here one has to carefully distinguish between 1 p ∈ X and 1 ∈ Y . In general, it is better to define subtraction independently of addition since it may happen that f (−x) is undefined. An important nontrivial example of f is provided by the Cantor function [7], or more precisely the Cantor-line function defined below, where X is the Cantor subset of [0, 1] and Y = [0, 1]. If 0 x exists, one can denote it by x.
Practically the only difference between ⊕, , , and +, −, ·, and / is that in general multiplication is not just a repeated addition: typically x ⊕ x = 2 x. Multiplication and addition are now truly independent.
Having all these arithmetic operations one can define a derivative of a function A : Employing (1) and the fact that f (0) = 0 one finds where F (y) = dF/dy is the usual derivative in Y , defined in terms of +, −, ·, and /. It is extremely important to note that (3) has been derived with no need of differentiability of f . f (0) = 0 is enough to obtain a well defined derivative. (3) is not the standard formula known for composite functions since no derivatives of f occur. To appreciate the difference between d/dx and . Then .

The formalism from [3] employs functions of the form
, hold true. The explicit form reads where F(y)dy is the standard (say, Lebesgue) integral in R.

Correspondence principle
Let us consider the two plots of the same function f (x) shown in Fig. 1. The right one is practically indistinguishable from f (x) = x, which means that numbers that are sufficiently large will satisfy an effectively Diophantine arithmetic and calculus, with Functions f that lead to such a correspondence principle can be characterized by Anticipating results of the next section let us show, in Fig. 2, the exponential function resulting from the solution of (7) and with arithmetic defined by f (x) from Fig. 1. As one can see, the resulting exponent involves two standard exponential regimes separated by a 'non-exponential' part. Still, one has to keep in mind that this awkward looking function does nevertheless satisfy the fundamental property (9), so this is just the exponential function corresponding to the non-Diophantine arithmetic and calculus implied by this concrete f . The exponent exp f from Fig. 2 has striking similarity to the expansion factor occurring in the CDM standard model of cosmology. The two 'standard' exponential regimes thus can be termed the inflationary and dark-energy parts of exp f x.
x and thus the resulting arithmetic is effectively Diophantine

Exponential function
To understand why functions of the form (2) are so essential let us solve the differential equation by assuming that A(x) = ⊕ ∞ n=0 a n x n , where x n = x · · · x (n times). Then, comparing term by term, one finds the unique solution Its inverse is 3 , in spite of its nonlinearity, keeps (7) unchanged. Accordingly, this is a symmetry of (7), a fact explaining the title of the paper.
Taking f (x) = x n with sufficiently large n, we would find a dynamics looking like a motion along a square. Note that the oscillator is described by the standard linear equation, even if we take a highly nonlinear function such as f (x) = x n . This shows that a change of As our next example consider a classical harmonic oscillator where ω 2 = ω ω. Setting x(t) = ⊕ ∞ n=0 a n t n , one obtains and C 1 , C 2 are constants.
An instructive exercise is to plot phase-space trajectories of the harmonic oscillator corresponding to various choices of f . Figure 3 shows the trajectories for the Cantor-line function, defined below, and f (x) = x n , with n = 1, 3, 5. All these trajectories represent a classical harmonic oscillator that satisfies the usual law of 'force oppositely proportional to displacement', with conserved energy 'ẋ 2 + ω 2 x 2 ', but with different meanings of 'plus' and 'times'. The resulting trigonometric functions are essentially the chirp signals [8] known from signal analysis.

Fractals: Cantor set
If one still has impression that what we do is just standard physics in nonstandard coordinates, consider the problem of a fractal Universe of dimension 4 − , analogous to the one arising in causal dynamical triangulation theory [9]. Our physical equations have to be formulated in terms of notions that are intrinsic to the Universe, but what should be meant by a velocity, say? We have to subtract positions and divide by time, but we have to do it in a way that is intrinsic to the Universe we live in. Moreover, from our perspective positions and flow of time seem continuous even if they would appear discontinuous from an exactly 4-dimensional perspective. We should not make the usual step and turn to fractional derivatives [10], since for inhabitants of (4 − )-dimensional Universe the velocity is just the first derivative of position with respect to time, and not some derivative of order 0 < α < 1.

Rotations in fractal plane
Now let us switch to higher-dimensional examples. First consider the plane, i.e., the Cartesian product of two lines. One checks that sin 2 It is an appropriate place to stress that our approach does not seem to be related to exponential, trigonometric, and hyperbolic functions defined (in various ways) in the context of time-scales dynamics [17][18][19], or in nonextensive thermodynamics [20,21], but there are links to Kolmogorov-Nagumo averages employed in Rényi information theory [21,24]  x = x cos f α ⊕ y sin f α, defines a rotation. The rotation satisfies the usual group composition rule, a fact immediately implying that one can work with generalized-arithmetic matrix equations.

Fractal Minkowski space
In an analogous way one arrives at Lorentz transformations in Cantorian Minkowski space, the Cartesian product of four Cantor lines with the invariant form 3 . Such Lorentz transformations are unrelated to those occurring in exactly 4-dimensional fractal space-time of scale relativity [14]. In 1 + 1 dimensions, employing, In consequence, the fact that f (1) = 1 sets the limit |β| ≤ 1 for maximal velocity independently of the choice of f . In principle, problems such as clock synchronization, composition of velocities, or the twin paradox, may lead to direct experimental tests of f and some insights into a putative fractal structure of space-time.

Harmonic oscillator and quantum mechanics
Arithmetic of complex numbers requires some care. One should not just take f : C → C due to the typical multi-valuedness of f −1 and the resulting ill-definiteness of ⊕ and . Definition of i as a π/2 rotation also does not properly work since one cannot guarantee a correct behavior of i n for a general f . The correct solution is the simplest one: one should treat complex numbers as pairs of reals satisfying the following arithmetic supplemented by conjugation (x, y) * = (x, y). The modulus reads As stressed in [15], the resulting complex structure is just the standard one, but no mysterious 'imaginary number' is employed.
In this way we have arrived at quantum mechanics. The most obvious example is the eigenvalue problem for a 1-dimensional harmonic oscillator. Consider where α, β are parameters. Since , then acting with f on both sides of (11) one finds where y = f (x). This is the standard harmonic oscillator equation with eigenvalues f (E n, f ) = f (α) f (β)(2n +1), and thus The eigenvalues of the quantum harmonic oscillator explicitly depend on the choice of arithmetic. The normalized ground state is , The excited states can be derived in the usual way.
There are two peculiarities of the resulting quantum mechanics one should be aware of. First of all, if f is a Cantor-like function representing a fractal whose dimension is less than 1, then the real-line Lebesgue measure of the fractal is zero. Keeping in mind that states in quantum mechanics are represented by equivalence classes of wave functions that are identical up to sets of measure zero, we can remove the Cantor line from R without altering standard quantum mechanics. Having removed the Cantor line C from R we still can do ordinary quantum mechanics on R \ C, whereas C itself can become a universe for its own, Cantorian theory. Removing C from R does not mean that we impose some fractal-like boundary conditions or consider a Schrödinger equation with a delta-peaked potential of Cantor-set support [25]. We just use the freedom to modify wave functions on sets of measure zero. So we can keep the standard Gaussian f (x) = x ground state on R \ C, and employ the Cantorian ψ 0 f (x) on C. According to quantum mechanics the resulting wave function belongs to the same equivalence class as the usual Gaussian, and thus represents the same state. However, now the energy ishω/2 + α β, with α β 'appearing from nowhere'. The analogy to dark energy is evident. The additional energy is a real number so it can be added tohω/2, similarly to many other energies that occur in physics and are additive in spite of unrelated origins.
The second subtlety concerns physical dimensions of various quantities occurring in f -generalized arithmetic. Even the simple case of ω t may imply a necessity of dimensionless ω and t if f is sufficiently nontrivial.
Functions of the form f (x) = x q , q ∈ R + , are in this respect exceptional since then x y = x y, and one can work with dimensional quantities. For example, the harmonic oscillator for f (x) = x q has energy levels In the formalism of Benioff [5,6], which employs a linear f (x) = px, one obtains E n, f = phω(n + 1/2).
In general, for example in a fractal space-time, we have to work with dimensionless variables x in order to make f (x) meaningful. It is thus simplest to begin with reformulating all the 'standard' theories in dimensionless forms, similarly to c = 1 andh = 1 conventions often employed in relativity and quantum theory.

Probability and entropy
Quantum mechanics has brought us to the issue of probability. An appropriate normalization is ⊕ k p k = 1 which, in virtue of f (1) = 1, implies k f ( p k ) = k P k = 1. We automatically obtain two coexisting but inequivalent sets of probabilities, in close analogy to probabilities P k and escort probabilities p k = P q k occurring in generalized statistics and multifractal theory [20,21]. Averages have the form of Kolmogorov-Nagumo averages [21][22][23], which implies the usual bounds a min ≤ a f ≤ a max . As a by-product of the construction, we have arrived at an entirely new interpretation of the Kolmogorov-Nagumo averaging: these are just the ordinary averages, but evaluated in a non-Diophantine arithmetic.
From the point of view of modified arithmetic the constraints one should impose on escort probabilities and Kolmogorov-Nagumo averages are, though, completely different from those employed in nonextensive statistics and Rényi's information theory [24], provided instead of additivity one has ⊕-additivity in mind. Rényi can be replaced by a much wider class of f s. In consequence, instead of Rényi's entropy one finds a general class of entropies parametrized by f , where P k = f ( p k ).

Remarks
The modified calculus is as simple as the one one knows from undergraduate education. What may be nontrivial is to find f if X is a sufficiently 'strange' object. The case of the Cantor line was relatively obvious, but the choice of f may be much less evident if X is a multifractal or a higher-dimensional fractal. Such more complicated examples are explicitly discussed in the sequel to the present paper [26].
In order to conclude, let us return to Figs. 3 and 4. All the phase-space trajectories represent the same physical system: A harmonic oscillator satisfying the Newton equation d 2 x/dt 2 = −ω 2 x, with the same physical parameters for each of the trajectories. So how come the trajectories are different? The answer is: Because the very form of Newton's equation does not tell us what should be meant by 'plus' or 'times'. This observation extends to any theory that employs arithmetic of real numbers. It would not be very surprising if some alternative arithmetic proved essential for Planck-scale physics, where fractal space-time is expected [27][28][29][30], or to biological modeling where fractal structures are ubiquitous. Links of f -generalized trigonometric functions to chirps suggest that f s even much more 'ordinary' than fractal functions may also lead to nontrivial applications. Spectrum of the harmonic oscillator is explicitly f -dependent, a fact showing that different versions of arithmetic lead to non-equivalent descriptions of physics. Fig. 4 A new type of rotational symmetry. Circles in Cantorian plane C 2 for ten different radii. The rotational symmetry is hidden in the sense that in order to reveal it one has to control the arithmetic intrinsic to the fractal set C. The circles are simultaneously phase-space trajectories of the harmonic oscillator from Fig. 3, but now with Cantor line f (x) Equation (7) and its solution (8) illustrate another important difference between the change of arithmetic and the usual change of variables, in spite of the misleading similarity of some formulas. Namely, note that a change of f from, say, f (x) = x to an arbitrary f does not influence the linearity of (7), even for a nonlinear f . This should be contrasted with nonlinear gauge transformations associated with Doebner-Goldin equations [31][32][33][34][35], that map a linear Schrödinger equation into a nonlinear one. Such a preservation of form is typical of all differential equations formulated within the generalized non-Diophantine paradigm. Another reason why one should not confuse the transformation F → f −1 • F • f = F f with a gauge transformation is the resulting rule for dF/dx → d f F f /d f x which is completely different from what one finds in gauge theories. Modification of arithmetic is a symmetry operation that works at a much more primitive level than a gauge transformation.