P-adic Properties of Time in the Bernoulli Map

Similar documents
Do we need Number Theory? Václav Snášel

A Geometry in which all Triangles are Isosceles

Valuations. 6.1 Definitions. Chapter 6

p-adic Analysis Compared to Real Lecture 1

Euclidean Models of the p-adic Integers

Nonconservative Abelian sandpile model with the Bak-Tang-Wiesenfeld toppling rule

Zair Ibragimov CSUF. Talk at Fullerton College July 14, Geometry of p-adic numbers. Zair Ibragimov CSUF. Valuations on Rational Numbers

(x x 0 ) 2 + (y y 0 ) 2 = ε 2, (2.11)

Controlling chaos in random Boolean networks

arxiv:chao-dyn/ v1 5 Mar 1996

Some topics in analysis related to Banach algebras, 2

Absolute Values and Completions

Mechanisms of Chaos: Stable Instability

LE JOURNAL DE PHYSIQUE - LETTRES

Chapter 8. P-adic numbers. 8.1 Absolute values

AN INTRODUCTION TO FRACTALS AND COMPLEXITY

arxiv:cond-mat/ v2 [cond-mat.stat-mech] 8 Sep 1999

Quantitative Description of Robot-Environment Interaction Using Chaos Theory 1

Algorithms for Picture Analysis. Lecture 07: Metrics. Axioms of a Metric

The p-adic numbers. Given a prime p, we define a valuation on the rationals by

Are chaotic systems dynamically random?

Introduction to Arithmetic Geometry Fall 2013 Lecture #5 09/19/2013

a = a i 2 i a = All such series are automatically convergent with respect to the standard norm, but note that this representation is not unique: i<0

arxiv:cond-mat/ v1 8 Jan 2004

1. Continuous Functions between Euclidean spaces

Takens embedding theorem for infinite-dimensional dynamical systems

ON THE ARROW OF TIME. Y. Charles Li. Hong Yang

Exercises for Unit V (Introduction to non Euclidean geometry)

MA257: INTRODUCTION TO NUMBER THEORY LECTURE NOTES

k-protected VERTICES IN BINARY SEARCH TREES

Classification of plane germs: metric and valorative properties

Workshop on Heterogeneous Computing, 16-20, July No Monte Carlo is safe Monte Carlo - more so parallel Monte Carlo

arxiv: v1 [math.co] 5 Apr 2019

NOTES ON DIOPHANTINE APPROXIMATION

Integer Factorization using lattices

QUANTUM CHAOS IN NUCLEAR PHYSICS

PHY411 Lecture notes Part 5

The Mathematica Journal p-adic Arithmetic

MATH 614 Dynamical Systems and Chaos Lecture 6: Symbolic dynamics.

ITERATED FUNCTION SYSTEMS WITH CONTINUOUS PLACE DEPENDENT PROBABILITIES

Chapter 2 Metric Spaces

Small worlds, mazes and random walks

Fractal Geometry Time Escape Algorithms and Fractal Dimension

AN INTRODUCTION TO FRACTALS AND COMPLEXITY

Móstoles, Spain. Keywords: complex networks, dual graph, line graph, line digraph.

NONLINEAR TIME SERIES ANALYSIS, WITH APPLICATIONS TO MEDICINE

UNM - PNM STATEWIDE MATHEMATICS CONTEST XLVI - SOLUTIONS. February 1, 2014 Second Round Three Hours

P-adic Functions - Part 1

No. 6 Determining the input dimension of a To model a nonlinear time series with the widely used feed-forward neural network means to fit the a

Global Attractors in PDE

The Behaviour of a Mobile Robot Is Chaotic

Smooth Ergodic Theory and Nonuniformly Hyperbolic Dynamics

1 Functions of Several Variables 2019 v2

2. Two binary operations (addition, denoted + and multiplication, denoted

1 Convexity, Convex Relaxations, and Global Optimization

1 Adeles over Q. 1.1 Absolute values

arxiv:cond-mat/ v1 [cond-mat.stat-mech] 13 May 2005

arxiv: v1 [math.mg] 5 Nov 2007

Fractals: A Mathematical Framework

= 1 2x. x 2 a ) 0 (mod p n ), (x 2 + 2a + a2. x a ) 2

Can two chaotic systems make an order?

THREE DIMENSIONAL SYSTEMS. Lecture 6: The Lorenz Equations

(x 1, y 1 ) = (x 2, y 2 ) if and only if x 1 = x 2 and y 1 = y 2.

Introduction to Real Analysis Alternative Chapter 1

Complex Behaviors of a Simple Traffic Model

Introduction to Dynamical Systems

Point Vortex Dynamics in Two Dimensions

Nonarchimedean Cantor set and string

14.1 INTRODUCTION. Nature. Architecture. Engineering. Compose a picture-album showing symmetry. Make some symmetrical paper-cut designs.

Chaos in GDP. Abstract

Mathematical Modelling Lecture 14 Fractals

GEOMETRY OF BINOMIAL COEFFICIENTS. STEPHEN WOLFRAM The Institute jor Advanced Study, Princeton NJ 08540

Local Fields. Chapter Absolute Values and Discrete Valuations Definitions and Comments

Liapunov Stability and the ring of P-adic integers

Non-reactive strategies in decision-form games

Addendum: Lyapunov Exponent Calculation

1 Shortest Vector Problem

1 Differentiable manifolds and smooth maps

Chapter 23. Predicting Chaos The Shift Map and Symbolic Dynamics

Topological properties of Z p and Q p and Euclidean models

SPECIAL CASES OF THE CLASS NUMBER FORMULA

Chapitre 4. Transition to chaos. 4.1 One-dimensional maps

One dimensional Maps

Lecture Notes DRE 7007 Mathematics, PhD

BASINS OF ATTRACTION FOR TRIPLE LOGISTIC MAPS. 1. Introduction

Biological Networks: Comparison, Conservation, and Evolution via Relative Description Length By: Tamir Tuller & Benny Chor

LYAPUNOV EXPONENT AND DIMENSIONS OF THE ATTRACTOR FOR TWO DIMENSIONAL NEURAL MODEL

Breaking an encryption scheme based on chaotic Baker map

DOUBLE SERIES AND PRODUCTS OF SERIES

Dynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325

[STRAIGHT OBJECTIVE TYPE] Q.4 The expression cot 9 + cot 27 + cot 63 + cot 81 is equal to (A) 16 (B) 64 (C) 80 (D) none of these

Material covered: Class numbers of quadratic fields, Valuations, Completions of fields.

Existence is Difference and is expressed and transformed by Information

Characterization of phonemes by means of correlation dimension

UNCC 2001 Comprehensive, Solutions

CONSTRUCTION OF THE REAL NUMBERS.

Dynamics of partial discharges involved in electrical tree growth in insulation and its relation with the fractal dimension

BAIBEKOV S.N., DOSSAYEVA A.A. R&D center Kazakhstan engineering Astana c., Kazakhstan

CHAOS -SOME BASIC CONCEPTS

Alignment indices: a new, simple method for determining the ordered or chaotic nature of orbits

Transcription:

Apeiron, Vol. 10, No. 3, July 2003 194 P-adic Properties of Time in the Bernoulli Map Oscar Sotolongo-Costa 1, Jesus San-Martin 2 1 - Cátedra de sistemas Complejos Henri Poincaré, Fac. de Fisica, Universidad de La Habana, Habana 10400, Cuba. (oscarso@ff.oc.uh.cu) 2 - Dpto. Fisica-Matemática y Fluidos, UNED, a.p. 60141, 28080 Madrid, Spain. The Bernoulli Map is analyzed with an ultrametric approach, showing the adequacy of the non-archimedean metric to describe in a simple and direct way the chaotic properties of this map. Lyapunov exponent and Kolmogorov entropy appear to yield a better understanding. In this way, a p-adic time emerges as a natural consequence of the ultrametric properties of the map. PACS: 05.45.+b, 02.90.+p Keywords: Non-Archimedean, ultrametricity, chaos Introduction Time is, perhaps, the most primary notion for all of us. Nevertheless (and maybe because of that) its properties, and specially its metric properties, have not been exhaustively analyzed. Contrary to the metric properties of the space and even of the space-time, where a

Apeiron, Vol. 10, No. 3, July 2003 195 profusion of research has appeared, there is a void in the metric description concerning time. In our opinion, ultrametricity concepts open a path to this description, and we try to illustrate it by studying the ultrametric properties of time in the case of the Bernoulli map. In ultrametric spaces, concepts such as exponential separation of neighboring trajectories, and characteristic parameters (Lyapunov exponents and Kolmogorov entropy) seem to yield a simpler understanding than the Euclidean metric. In the last years, ultrametricity has triggered interest in a wide range of physical phenomena, due to its applications in different fields: spin glasses, mean field theory, turbulence, and nuclear physics. Also optimization theory, evolution, taxonomy, and protein folding benefit from it (for an excellent review see [1]). Wherever a hierarchical concept appears, non-archimedean analysis is an adequate tool to study the problem. As an example where Euclidean metric is not very adequate, let us consider Baker's map [4]. The interval [ 0,1] [0,1] is mapped into [ 0,1] [0,1]. Therefore, the distance between two points cannot be larger than the distance between two opposite corners in [ 0,1] [0,1]. Nonetheless, Baker's map has a Lyapunov exponent greater than one. Thus, the distance between neighboring points grows exponentially in a finite region of the phase space. In the Euclidean space we would have to define the distance as the Euclidean length of the shortest path lying entirely within the region. Ultrametricity is a promising tool in the theory of branching processes, which, at the same time, has shown promise in the study of self-organized critical processes where a branching representation can be introduced [2,3]. It seems possible to find simpler tools to describe the geometry of these processes.

Apeiron, Vol. 10, No. 3, July 2003 196 Here, we illustrate the advantages of a hierarchical representation in the case of the Bernoulli shift. This will enable us, using simple geometric considerations, to determine the magnitudes governing the system, and the advantages of a p-adic metric will be stressed over the Euclidean metric. The ultrametric distance will be shown to be consistent with the characteristic behavior of this chaotic onedimensional map. In this paper we explore the application of ultrametricity in linking the Bernoulli map with a branching structure, which will reveal the possibilities of assigning an ultrametric measure to processes that, to all appearances, are not linked with a given metric (e.g., minority game and related problems) so that an adequate understanding of the ultrametric properties of a given process may lead to deeper understanding [5]. As any nontrivial norm is equivalent to the Euclidean or any of the p-adics (Ostrowski's theorem [1]), it would be convenient to measure the distance between points in Baker's map with a p-adic metric. Ultrametricity An ultrametric space is a space endowed with an ultrametric distance, defined as a distance satisfying the inequality d( A, C) Max{ d( A, B), d( B, C)} (1) (A, B and C are points of this ultrametric space), instead of the usual triangular inequality, characteristic of Euclidean geometry d ( A, C) d( A, B) + d( B, C) (2) A metric space E is a space for which a distance function d ( x, y) is defined for any pair of elements (x,y) belonging to E. A norm satisfying

Apeiron, Vol. 10, No. 3, July 2003 197 x+ y Max{ x, y } (3) is called a non-archimedean metric, because equation (3) implies that x + x x (4) holds, and equation (4) does not satisfy the Archimedes principle: x + x x (5) A metric is called non-archimedean or ultrametric, if (1) holds for any three points (x, y, z) dxz (, ) Max{ dxy (, ), d( yz, )} (6) A non-archimedean norm induces a non-archimedean metric: dxz (, ) = x z Max{ d( x, y), d( yz, )} (7) It is known that equation (7) implies a lot of surprising facts, e.g., that all triangles are isosceles or equilateral and every point inside a ball is itself at the center of the ball, while the diameter of the ball is equal to its radius. An example of ultrametric distance is given by the p-adic distance, defined as d p ( x, y) = x y (8) where the notation defines the p-adic absolute value: x p r p (9) p where p is a fixed prime number, x 0 is any integer, and r is the highest power of p dividing x. Two numbers are p-adically closer as long as r is higher, such that p r divides x y. Amazingly, for p = 5 the result is that 135 is closer to 10 than 35.

where Apeiron, Vol. 10, No. 3, July 2003 198 Any positive or negative integer can be represented by a sum = i x a i p (10) i=0 0 a i p 1 (11) If negative exponents are considered in the sum, rational numbers can also be represented. Such a representation is unique. The set of all sums Q p is the field of p-adic numbers, and contains the field of rational numbers Q but is different from it. Lyapunov exponent, Kolmogorov entropy and ultrametric time. With the above description the p-adic numbers have a hierarchical structure, whose natural representation is a tree. Let us now use this description to work with the Bernoulli map (See [4]): xn+ 1 = 2xn mod1 (12) n = 0,1,2... Here, we may note that the numbers can be represented as a set of points in a straight line or by a hierarchical structure, depending on the definition of distance (Euclidean or Archimedean) as we will see below. Let us represent the initial value (state) to be mapped into the unit interval by the sequence 0, a 1... a N... with a i = 0 or 1 to denote the initial value in binary notation. It is possible to reorder these sequences as a hierarchical tree. To do this, let us perform the following process to represent the result of the application of the Bernoulli map:

Apeiron, Vol. 10, No. 3, July 2003 199 We begin at an arbitrary point. We read, consecutively, the values of a i, from i = 1 to N, of the sequence a 1... an... When a i takes the value 0 we move to the left, and the same distance down. When a takes the value 1 we do the same, but moving to the right. The i N result is 2 branches of a hierarchical tree. Any finite path inside this branching structure represents univocally a possible finite sequence a 1... an... Thus, for instance, the sequence 0,0110 represents: left, right, right, left. The distance d ( x i, x j ) between two branches (sequences) x i, x j in this tree is given by ( m n) 2 i j d( xi, xj) = (13) 0 i = j where m is the number of levels one must move up the tree to find a common branch linking x i and x j, and N is the number of levels (the length of the sequence). This is equivalent to h 2 i j d( xi, xj) = (13a) 0 i = j where h is the position of the last block a h in which a i (i = 1,...,h) are common to the two sequences x i, x j. This means that the numbers x i and x j are near the h th binary place. This distance is ultrametric. To calculate the Lyapunov exponent it is necessary to know how neighboring points x0 + ε and x 0 evolve during the Bernoulli map. If h δ1 δ2 N ε be equal to 2 [1+ 2 + 2 +...] > 2, then the first different position between x0 = 0, aa 1 2... ah 1aN... and x0 + ε is a h.

Apeiron, Vol. 10, No. 3, July 2003 200 It is then is necessary to move up the tree N h + 1 levels from the bottom line to find the common branch in the position a h 1 (obviously, the last common figure between x 0 and x0 + ε ). Thus, and d( x, ) 2 h ( + 1) 0 + ε x0 = (14) d( f ( x + ε), f ( x )) = 2 (15) n n h ++ 1 n 0 0 because the iteration f n moves away the common branch n positions from the bottom level. To calculate the Lyapunov exponent it is necessary to express the exponential growth of the distance between two neighboring points: limlim2 ε = limlim d( f ( x + ε), f ( x )) (16) λ n n n 0 0 n ε 0 n ε 0 Since the base for measuring the p-adic distance in our space is the number 2, in the preceding equation we have expressed the exponential growth as 2 λn n instead of e λ. n n Replacing ε and d( f ( x0 + ε), f ( x0)) in the preceding equation we obtain h δ1 δ2 λn h++ 1 n limlim2 (1+ 2 + 2 +...)2 = limlim2 (17) n h n h from (17) it can be easily observed that λ = 1. As the Lyapunov exponent in the Bernoulli map is ln 2 [4], we ln2 recover this result with p-adic metric, since 2 = e. This means that each unit time interval implies a new doubling of branches in each node of the hierarchical tree. Once a unit time interval has elapsed, the number of levels one must move up the tree to find a common branch increases by one. This result will be crucial to understanding how information is lost in the course of time.

Apeiron, Vol. 10, No. 3, July 2003 201 In one-dimensional maps like the one considered here, the Kolmogorov entropy coincides with the Lyapunov exponent [4]. The expression for the Kolmogorov entropy is: 1 K = limlim pi 1... i lg2 p1... 0 n i i (18) n τ n nτ where i1... in i1... in p is the probability to reach the i n -th state of the system in the phase space following a given path i 1 i 2 i n. It can be seen that in our case this probability only depends on the final state i n, because for each state there is just one path, i.e., that given by the sequence i 1 i 2 i n. Moreover, the number of states in the n th level is 2 n, and τ is the time elapsed in passing from one state to the next. The probability of occupying one of the 2 n 1 states is pn = pii 12... i =, and this results n n 2 in 2 K = limlim (19) n τ 0 2 n τ But the distance between two successive states of the n th level is 2 1 n, because they are common up to the (n 1) th level. Since the speed v to pass from one sequence to the next is constant in the 1 n 2 Bernoulli map, i.e., v = = 1, the time τ elapsed between these τ two successive states is τ = 2 1 n. As expected, k = 1, coinciding with the Lyapunov exponent. Hence, time has ultrametric properties in this process. Notice that the existence of a p-adic proper time is essential to the coincidence of the Kolmogorov entropy and the Lyapunov exponent. The spatial p-adic structure is unavoidably joined to the p- adic structure of proper time.

Apeiron, Vol. 10, No. 3, July 2003 202 Therefore, we can say that this problem possesses a p-adic spatial and temporal geometry instead of a sole p-adic spatial geometry. To see the importance of the introduction of a p-adic time, see [6]. The Kolmogorov entropy measures the loss of information in the process. This loss of information can easily be seen from our representation, since the process of separation of trajectories is such that for any step the increase of the distance between two points duplicates the number of branches through which this increment can be reached. We are loosing information because we do not know exactly the way we are separating two states. On the other hand, we can see that in the ultrametric space the natural time of the system is also ultrametric. The time of transition between two sequences x i, x j satisfies the same expression (13) as the distance between x i,x j. Furthermore, the subsequent behavior of two states that separate at a given point in the ultrametric space depends on the point at which separation occurs, revealing that ultrametricity can be applied to decision processes (like minority games, aging effects, hierarchical processes, etc.), where ultrametric concepts have been poorly applied. The application of ultrametricity to the minority game will be treated in future work. Conclusions It has been shown that the Bernoulli map leads to a hierarchical structure in the p-adic metric. In light of the ultrametric distance, the Lyapunov exponent and the Kolmogorov entropy can be better understood and a direct geometric interpretation is supplied by the hierarchical structure. The p-adic metric seems to be the natural metric of this map. The hierarchical structure generates p-adic properties for temporal evolution.

Apeiron, Vol. 10, No. 3, July 2003 203 This work was partially supported by an Alma Mater prize, Havana University, and by the Ministerio de Educacion y Cultura, Spain. References [1] Rammal, R., G. Toulouse and M. A. Virasoro, Ultrametricity for Physicists, Rev. Mod. Phys 58, 3, (1986) 765-788. [2] Garcia-Pelayo, R., Dimension of branching processes and self-organized criticality, Phys. Rev. E 49, 6, (1994) 4903-4906. [3] Sotolongo-Costa, O., Y. Moreno-Vega, J.J. Lloveras and J. C. Antoranz, Criticality in droplet fragmentation, Phys. Rev. Lett 76, 1, (1996) 42-45. [4] H. G. Schuster Deterministic Chaos Physik-Verlag, (1984) 216pp. [5] M. Barnsley Fractals Everywhere Academic Press (1988) 394pp. [6] D.K. Arrowsmith and F. Vivaldi, Geometry of p-adic Siegel disks, Physica D 71, (1994) 222-236.