Generalized Fibonacci numbers, cosmological analogies, and an invariant

Continuous generalizations of the Fibonacci sequence satisfy ODEs that are formal analogues of the Friedmann equation describing spatially homogeneous and isotropic cosmology in general relativity. These analogies are presented, together with their Lagrangian and Hamiltonian formulations and with an invariant of the Fibonacci sequence.


I. INTRODUCTION
Fibonacci, also known as Leonardo Pisano or Leonardo Bonacci, introduced Hindu-Arabic numerals to Europe with his book Liber Abaci in 1202 [1]. He posed and solved a well-known problem involving the growth of a population of rabbits in idealized situations. The solution, now known as the Fibonacci sequence 0, 1, 1, 2, 3, 5, 8, 13, ... is written as where F n is the n-th Fibonacci number and F 0 = 0, F 1 = 1. Furthermore, the ratio of two consecutive terms F n+1 /F n approaches the golden ratio as n → +∞. The Fibonacci sequence can be generalized to the continuum using Binet's formula Furthermore, the analytic function reproduces part of the Fibonacci numbers, F n = F (e) (n) for even x = n ∈ N, while reproduces the other Fibonacci numbers for odd x = n ∈ N. The function reproduces the entire Fibonacci sequence for x = n ∈ N. Here we focus on F (e,o) (x), which admit analogies with relativistic cosmology, while no such analogy exists for F (x). The functions F (e,o) satisfy the dual relations (where a prime denotes differentiation with respect to x) of which F (e) and F (o) are two linearly independent solutions. In physics, this equation describes the one-dimensional motion of a particle of position F in the inverted harmonic oscillator potential V (F ) = −kF 2 /2 (with K = 2 ln 2 ϕ), which is used as an example of an unstable mechanical system (this property corresponds to the fact that the Fibonacci numbers F n increase without bound as n → ∞).
We will also use which is a Fibonacci function according to the definition of Ref. [2], i.e., However, F (e,o) (x) are not Fibonacci functions since they contain ϕ −x , which is not a Fibonacci function (it is easy to prove [2] that, among power-law functions y(x) = b x , the only Fibonacci function is the one with base equal to the golden ratio, b = ϕ).
In the following section, we briefly recall the basics of spatially homogeneous and isotropic cosmology in general relativity, and then we present the formal analogy with Eqs. (1.7)-(1.9) in Sec. III. We then deduce Lagrangian and Hamiltonian formulations for the Fibonacci ODEs and derive from these an invariant of the (discrete) Fibonacci sequence.

II. FLRW COSMOLOGY
In the context of Einstein's relativistic theory of gravity, the most basic assumption of cosmology is the Copernican principle stating that, on average (i.e., on scales larger than a few tens of Megaparsecs), the universe is spatially homogeneous and isotropic [3][4][5][6][7][8]. The stringent symmetry requirements of spatial homogeneity and isotropy force the geometry to have constant spatial curvature [9]. The spacetime metric is necessarily given by the Friedmann-Lemaître-Robertson-Walker (FLRW) line element, written in polar comoving coordinates (t, r, ϑ, φ) as [3,4] where dΩ 2 (2) = dϑ 2 + sin 2 ϑ dφ 2 is the line element on the unit 2-sphere and the sign of the constant curvature index K classifies the 3-dimensional spatial sections t = const. K > 0 corresponds to positively curved 3-spheres, K = 0 to Euclidean flat sections, and K < 0 to hyperbolic 3-sections. The dynamics of FLRW cosmology is contained in the scale factor a(t). In FLRW cosmology, the matter source causing the spacetime to curve is usually (but not necessarily) taken to be a perfect fluid with energy density ρ(t) and pressure P (t) = wρ(t), where w is a constant equation of state parameter. The evolution of a(t) and ρ(t) is ruled by the Einstein equations adapted to the high degree of symmetry, the Einstein-Friedmann equations. They comprise [3][4][5] the Friedmann equation (a first order constraint), the acceleration or Raychaudhuri equation (contrary to FLRW cosmology, these quantities are dimensionless in the Fibonacci side of the analogy). A priori the formal equivalence with the Friedmann equation (2.2) is not sufficient for the analogy to hold and one must check that also Eqs. (2.3) and (2.4) are satisfied: this is straightforward to do using Eq. (1.9), while the conservation equation is trivially satisfied with ρ = P = 0 (and also by the Λ-effective fluid with P = −ρ = const.). The universe with scale factor analogous to the function F (e) (x) has hyperbolic 3-D spatial sections and is empty, but expands due to the repulsive cosmological constant. By squaring Eq. (3.1), one has introduced the possibility of solutions with F ′ (e) < 0, corresponding to a contracting, instead of expanding, universe. In fact, Eq. (2.2) then giveṡ where C and t 0 are integration constants (C serves the purpose of making the argument of the logarithm dimensionless since, in the units used [3][4][5], the scale factor a carries the dimensions of a length and Λ those of an inverse length squared). By exponentiating both sides, squaring, and collecting similar terms, one is left with with a 0 constant. In practice, in FLRW cosmology the dimensionless radial coordinate r can be rescaled to normalize the curvature index K (usually to ±1); here we can rescale to obtain |K|ΛC 2 = 1/3 and where, in order to keep the scale factor non-negative, the upper sign applies for t ≥ t 0 and the lower one for t ≤ t 0 . At late times t → +∞, this metric is asymptotic to the metric of de Sitter spacetime a (dS) (t) = a 0 e √ (3.13) The analogous FLRW universe is again empty and propelled by the positive cosmological constant, but the spatial sections are now closed 3-spheres. Again, by squaring one introduces the possibility of a negative sign forȧ. Proceeding in parallel with what has been done for F (e) , one obtains (3.14) and, normalizing K suitably, This scale factor describes a universe contracting from a infinite size, bouncing at the minimum value a 0 , and then expanding forever and asymptoting to the de Sitter space, a behaviour sought for in quantum cosmology to avoid the classical Big Bang singularity a = 0. Finally, we can consider the Fibonacci function f (x) = ϕ x , which trivially satisfies and is analogous to an empty (ρ = P = 0), spatially flat (K = 0) universe expanding exponentially, a (dS) = a 0 e √ Λ 3 t due to the positive cosmological constant Λ = 3 ln 2 ϕ. This is the maximally symmetric de Sitter universe, which is an attractor in inflationary models of the early universe [6,7] and in dark energy-dominated models of the late (present-day) accelerating universe [10].

IV. LAGRANGIAN AND HAMILTONIAN
Both the mechanical analogy of Eq. (1.9) and the cosmological analogy suggest the Lagrangian for F (e,o) (x) . which gives the first integral F 2 (o,e) − F 2 (e,o) on the Fibonacci side of the analogy. In terms of the discrete Fibonacci sequence we have, therefore, the

Proposition:
The quantity does not depend on m and is an invariant of the Fibonacci sequence.
To the best of our knowledge, this invariant is not related to known invariants (for example, those of the Fibonacci convolution sequences [11]).

V. CONCLUSIONS
The functions F (e,o) associated with the continuum generalization of the Fibonacci sequence exhibit analogies with certain spatially homogeneous and isotropic universes in FLRW cosmology, and also a mechanical analogy with an inverted harmonic oscillator, which we have presented. With the help of these analogies, it is easy to derive a Lagrangian and a Hamiltonian associated with the ODEs satisfied by the functions F (e,o) . Then, the conservation of the Hamiltonian (or energy) yields a quantity I, given by Eq. (4.8) that is constant across the Fibonacci sequence, i.e., does not depend on the index n. This invariant may turn out to be useful in view of the many applications of the Fibonacci numbers in mathematics and in the natural sciences (e.g., [12][13][14][15][16][17][18]).