JAMPJournal of Applied Mathematics and Physics2327-4352Scientific Research Publishing10.4236/jamp.2019.712214JAMP-97065ArticlesPhysics&Mathematics Generalized Fibonacci Sequence: Possible Template for the Constants of Nature SebastianoTosto1*(Retired Physicist) ENEA Casaccia, via Anguillarese 123, Roma, Italy0212201907123037306613, November 20199, December 2019 12, December 2019© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

The paper concerns the fundamental constants of physics. The fine structure constant is inferred first with the help of a generalized Fibonacci sequence. Next this generalized sequence is also implemented in order to show the interconnection between various fundamental constants. The analysis of an extended definition of Fibonacci sequence reveals that the fundamental constants take the meaning of archetypal templates that model the physical appearance of the observable Universe. Once more appears in this conceptual frame the natural and serendipitous link between quantum and relativistic theories.

Fibonacci Sequence Fundamental Constants
1. Introduction

The standard definition of Fibonacci sequence is

F = 1 , 1 , 2 , 3 , 5 , 8 , 13 , ⋯   (1.1)

quoted in the following as F = F ( 1.1 ) to emphasize the first two terms of the sequence. The first two numbers are crucial in determining all terms of the sequence, being by definition each term sum of the two preceding ones. This sequence was originally aimed to calculate the reproduction rate of rabbits under appropriate fertility hypotheses: the first two terms are required equal to one to signify the first couple of rabbits. If each couple becomes fertile at the end of the first month and gives birth to a new couple when the second month is completed, then the meaning of each term of the sequence is the number of fertile couples per month. From a mere mathematical point of view, it is well known why the Fibonacci numbers are related to the so-called “golden ratio”. In principle, however, no conceptual reason requires input values just equal to 1 if the sequence is not essentially restricted to the original purpose to which it was early aimed. In other words, nothing hinders to generalize this sequence in order to describe a wider range of natural phenomena simply changing the first two terms. Here is proposed the generalization of F ( 1.1 ) putting

F * = F * ( f 0 , f 1 ) , (1.2)

where f 0 and f 1 are two arbitrary numbers in principle not necessarily integers; this affects all successive terms. Clearly, the standard Fibonacci sequence F is uniquely defined; instead the mathematical implications of the modified sequence F * , and thus their related physical meaning as well, just depending on how are defined the initial terms. Defining indeed

F * = f 0 , f 1 , ( f 1 + f 0 ) , ( 2 f 1 + f 0 ) , ( 3 f 1 + 2 f 0 ) , ( 5 f 1 + 3 f 0 ) , ( 8 f 1 + 5 f 0 ) , ⋯ (1.3)

i.e.

F 1 * = f 0 ,   F 2 * = f 1 ,   F 4 * = 2 f 1 + f 0 ,   F n * = F n ( f 0 , f 1 ) = f 1 F n − 1 + f 0 F n − 2 , (1.4)

the last equation evidences that now the standard values F n still appear in F n * as a linear combination of F n − 1 and F n − 2 with coefficients f 0 and f 1 . As it actually means that the ( n − 1 ) -th and ( n − 2 ) -th terms of the sequence (1.2) are all multiplied by f 1 and f 0 respectively, (1.4) reads identically

F n * = F n − 1 ( 1 ) + F n − 2 ( 0 ) ,   F n ( 1 ) = f 1 F n − 1 ,   F n − 2 ( 0 ) = f 0 F n − 2 ,   n > 2 (1.5)

the inequality reminds that the standard coefficients of (1.3) appear indeed for n > 2 , where are in fact calculable all F n * as sum of two previous F ( 0 ) and F ( 1 ) now defined by the arbitrary f 0 and f 1 . The recursive rule is therefore unchanged. Since the single terms of the new sequence F n * are actually linear combinations of two terms with different n of the standard sequence F n , reasonably each addend of the former could have in principle its own physical meaning likewise as the number calculated in the latter. It is therefore reasonable the guess why the aforesaid ability of F in describing the natural phenomena should still hold and be even enhanced by the additional freedom degrees inherent F * . This makes significant the appropriate definition of selected specific terms F 1 * and F 2 * suitable to provide information of interest about specific physical effects, likewise as each F n yields the number of couples of rabbits reproduced in the n-th generations from the initial couple. To exemplify this point return back to the original Fibonacci sequence, where each term represents the result of all previous reproduction steps. There is no way to introduce in this (1.1) the chance that either partner generated from the initial couple, fertile by definition, could actually be sterile; in other words, all successive generations of rabbits are assumed healthy and fertile exactly as the former one. The common sense suggests however that in practice, for any biological reason, this is not realistic; is missing in the early Fibonacci calculation any reference to possible discontinuity caused by an external action affecting the reproduction regularity of rabbits. From a physical point of view this chance could instead take place for example during the evolution of a system formed by an increasing number of particles from its initial configuration, e.g. early number of allowed states, to its current configuration after an appropriate number of intermediate steps; thus the evolutionary implication inherent the Fibonacci sequence can be advantageously modified and tailored to specific physical systems through the first two terms, as suggested by (1.4). In principle the deterministic character of (1.1) is inconsistent with the chance of describing quantum phenomena; also, the early F necessarily concerns phenomena describable via pure numbers. Instead the addends of F n * , as proposed here, skip this requirement, i.e. the various terms are required to be neither integers nor dimensionless only; is crucial in this respect the numerical and dimensional choice of the coefficients f 0 and f 1 . In practice, the first few terms of F * read for example

F 3 * = 2 f 1 + f 0 ,   F 7 * = 13 f 1 + 8 f 0 , F 15 * = 610 f 1 + 377 f 0 ,   F 20 * = 6765 f 1 + 4181 f 0 ; (1.6)

i.e. the coefficients that multiply f 1 and f 2 are just the respective numbers F n − 1 and F n − 2 of the standard sequence. So it seems worth investigating the chances opened by the kind of generalization prospected by (1.4), in particular the choice of the terms F n * pertinent to the specific physical problem of interest. For example, it is possible to multiply F n * by an arbitrary dimensional factor k c in order that (1.6) reads

E 3 * = ( 2 f 1 + f 0 ) k c ,   E 7 * = ( 13 f 1 + 8 f 0 ) k c , E 15 * = ( 610 f 1 + 377 f 0 ) k c ,   E 20 * = ( 6765 f 1 + 4181 f 0 ) k c ; (1.7)

the notation exemplifies the particular case where

k c = unit   energy . (1.8)

So the terms of (1.6) and (1.7) are numerically identical although the latter represents a sequence of energies expressed in energy units corresponding to k c , whereas the pure numbers of the standard Fibonacci sequence turn into energies directly related to the evolution of a physical system: the steps that describe the total number of the n-th generation of rabbits turn into the allowed energies of the system at various time steps. In the following we implement (1.6) for simplicity of notation, subtending however that when necessary the numerical terms of the generalized sequence F n * represent identically the physical meaning of the sequence

F n * ≡ E n * k c . (1.9)

So (1.4) will be also implemented in a form immediately consequent putting

f 0 = ζ ″ f ,   f 1 = ζ ′ f ; (1.10)

owing to (1.5) it is possible to write (1.4) as exemplified in (1.6)

F n * = ( ζ ′ F n − 1 + ζ ″ F n − 2 ) f ,   f = f ( ζ ) , (1.11)

where ζ along with ζ ′ and ζ ″ are arbitrary parameters. Moreover, without loss of generality, it is sensible the chance of regarding the n-th terms of F n * via the definitions

f ( ζ ) = σ 1 α ∘ + σ 2 α ∘ 2 + ⋯ , ζ ′ = 1 + σ ′ 1 α ∘ + σ ′ 2 α ∘ 2 + ⋯ , ζ ″ = 1 + σ ″ 1 α ∘ + σ ″ 2 α ∘ 2 + ⋯ , (1.12)

i.e. replacing the first two numbers of the standard sequence (1.1) with the series expansions (1.12) defining the parameters ζ ′ and ζ ″ , which just for this reason are both defined with 1 as zero order approximation value. The new information introduced via (1.12) is that now appears a unique parameter α ∘ as a function of which is in principle calculable each F n * via suitable values of the various coefficients σ ′ j and σ ″ j , by definition uniquely fixed once for all and focused on the physical problem of interest. Note that if σ ′ j and σ ″ j are defined as

∑ j σ ′ j a ∘ j ≥ 0 ,   ∑ j σ ″ j a ∘ j ≥ 0 ,   ∑ j σ j a ∘ j = k c ,

then, whatever the specific value of a ∘ might be, (1.11) reduces to either (1.6) or (1.7) depending on whether k c = 1 or for example k c = u n i t   e n e r g y ; clearly the latter chance implies regarding σ j as dimensional factors with physical dimensions of energy according to (1.9).

So the positions (1.12) seem reasonable to define (1.11) consistently with particular cases of standard sequence (1.1) or extended sequence (1.6) or even dimensional extended sequence (1.7).

Also, σ j = 0 imply f ( ζ ) = 0 and F n * = 0 , which means that for α ∘ = 0 the terms F n > 2 * = 0 are identically null. This makes sense because (1.11) and (1.12) hold for n > 2 ; so, with reference to the standard (1.1), F 1 * = 1 and F 2 * = 1 are enough to return back to the case of an initial couple of infertile rabbits.

Nevertheless, depending on the specific physical problem, it is possible that the general (1.11) admits useful simplifications; putting for example ζ ′ = ζ ″ , one finds

F n * * = ( F n − 1 + F n − 2 ) ζ ′ f . (1.13)

This particular case is interesting because F n * * * ≈ F n − 1 + F n − 2 is compliant just with the original Fibonacci intuition; indeed, recalling (1.9), ζ ′ f is the general form of the factor k c , which in turn can become itself mere dimensional factor. Furthermore, it is also possible that the series expansions (1.12) can be truncated to the respective first order approximation terms; if so, then (1.11)

F n * * * ≈ ( F n − 1 + F n − 2 ) ( 1 + σ ′ 1 α ∘ + σ ′ 2 α ∘ 2 + ⋯ ) ( σ 1 α ∘ + σ 2 α ∘ 2 + ⋯ ) (1.14)

reduces to

F n * * * ≈ ( F n − 1 + F n − 2 ) ( 1 + σ ′ 1 α ∘ ) σ 1 α ∘ . (1.15)

Finally it is clear that just for this reason, despite all manipulations so far introduced, holds anyway

F n * * * > 0 (1.16)

this condition is necessary in order that F n * * * and the initial F n is conceptually consistent. It is not enough that both are defined by the same sum rule of successive terms, it is also necessary that all terms are in both cases still positive although numerically different. If so, then (1.16) can be legitimately regarded as generalized Fibonacci sequence even for α ∘ < 0 . Note that (1.13) is simply a particular case of (1.11) but in principle it does not contain any approximation, as instead (1.15) does. Thus these last equations must be motivated by and justified for the specific physical problem of interest case by case. Nevertheless both contain the essential feature that substantiates the aforesaid premise of the present paper, i.e.: 1) the physical dimensions of the parameter α ∘ , and thus of σ ′ too, to fit the specific physical problem and 2) the actual chance of generalizing (1.1) through appropriate definitions of both parameters. Remains still evident the link between (1.15) and (1.1), which is identically reproduced putting in particular α ∘ = 1 and σ ′ 1 = σ 1 = 0 . It is clear however that (1.15) does not require all higher order approximation coefficients σ ′ j > 1 = 0 , as it concerns the first order of approximation only. So, despite the mere numerical character of the approximation (1.15), any deviation from these boundary conditions is itself a sensible generalization of (1.1). If F n * * * is known for any n, (1.15) can be trivially solved with respect to α ∘ ; i.e.

α ± o = − σ 1 ± σ 1 2 + 4 σ ′ 1 σ 1 R 2 σ ′ 1 σ 1 ,   R = F n * * * F n − 1 + F n − 2 , (1.17)

where σ 1 > 0 and R > 0 by definition whereas σ ′ 1 can be in principle positive or negative.

The aim of the present paper is to show what has to do F * with the fine structure constant α , whose acknowledged numerical value is

α = 0.007297352664. (1.18)

As concerns this number, recall that implementing

e = 4.8025 × 10 − 10 e .s .u . ,   G = 6.67430   10 − 8 cm 3 ⋅ g − 1 ⋅ s − 2

one finds that 

e = k g α G ,   k g = 0.986   g 3 / 2 ⋅ s ⋅ cm − 3 / 2   i . e .   e G ℏ c   ≈ 1.01   cm 3 / 2 ⋅ g − 3 / 2 ⋅ s − 1 , (1.19)

being k g dimensional proportionality constant. The fact that this constant is numerically very close to 1 cannot be accidental, rather it shows that reasonably α provides a direct physical link between e and G. The explanation and related implications have been concerned in  ; however here it is worth emphasizing that even this relationship must be someway inferable in a general conceptual frame aimed to understand profoundly the physical meaning of α and to calculate its numerical value.

The standard Fibonacci sequence is, per se, a mere recursive list of pure numbers; nevertheless it actually appears in several events and objects occurring in Nature    . For example it is known that it describes the shape of nautilus shell and spirals during hurricane formation, undergoing cancer cell division and filament profile during spiral galaxy formation, chicken egg contour and curls of numerous insects   , to name just a few. In principle nothing requires the correspondence between sequence of pure numbers and natural occurrences quantifiable through lengths, energies and momenta; as a matter of fact, however, the spirals consisting of arcs of circles inscribed into contiguous squares with size ratios progressively increasing in agreement with the Fibonacci numbers effectively overlap the observed shapes of galaxy arms and insect curls. The chance of examining these dimensionless ratios is just an example of how even pure numbers can take the worth of observable physical evidence. A huge literature is known about how and why this sequence is surprisingly suitable in describing natural phenomena . The chance of acknowledging the existence of well-defined fingerprints while observing largest and smallest natural phenomena cannot be accidental; rather it suggests that, although surprisingly, a mere sequence of pure numbers becomes the key criterion underlying occurrences usually measured through dimensional quantities like lengths or times or velocities and explained through familiar concepts of physics like minimum energy or maximum probability.

In other words, here is the first leading idea of the present paper: it seems sensible to expect that the numbers (1.1) are actually a sort of core boundary condition that controls even the fundamental constants of Nature.

In this way, regardless of how these constants or their combinations govern fields or interactions, these numbers can be identified as templates to which conform from time to time any physical observable; the measure units and their physical dimensions become therefore simply factors to turn the experimental observation into quantifiable data and predicting ability of standard scientific research.

Here is thus the second leading idea of this paper: to demonstrate the existence of hidden templates to which Nature conforms through its fundamental constants.

In this respect the indirect correspondence between pure numbers and natural events seems a limiting restriction to the actual worth of (1.1) and its full implementability; nevertheless it suggests the realistic chance of describing observable phenomena directly through measurable parameters someway related to (1.4), which requires in turn to convert the pure numbers of the sequence (1.1) into dimensional quantities like (1.9) with specific physical meaning. The basic step in this respect is the idea of extending the early Fibonacci sequence in order to calculate both dimensionless and dimensional constants of the Nature, like α or proton to electron mass ratio r p e or Avogadro number N A on the one hand and electron charge e or Boltzmann constant K B or G on the other hand.

In fact, this will be done implementing appropriately the coefficients f 0 and f 1 .

To understand this point in agreement with (1.7) and (1.9), note that according to (1.5) F * is sum of two sequences obtained multiplying (1.1) a first time by f 0 and a second time by f 1 and next collecting the particular terms f 1 F n − 1 and f 0 F n − 2 as in (1.4). Analogously multiply all terms of F a first time by an arbitrary dimensional factor k a f 1 , so that all F n turn into k a f 1 F n , and next by another arbitrary dimensional factor k b f 0 , so that all F n turn now into k b f 0 F n ; the two sequences obtained in this way consist of terms with well defined physical dimensions, in general different when k a ≠ k b . If for example the dimensional factor k a is defined as cm−1∙g−1∙s−1, then the f 1 k a means f 1 per unit length, unit mass and unit time. If more realistically k a is unit energy and k b unit momentum, then the sequence (1.1) early introduced to describe uniquely the population increment of rabbits now generates two separate sequences concerning possible time evolutions of energy and momentum of a physical system; note however that nothing changes from a mere numerical point of view, i.e. energy and momentum sequences are still defined by (1.1). Eventually, multiply again these sequences by f ′ 0 and f ′ 1 whose physical dimensions are defined by k c / k a and k c / k b as follows

F ′ n ∗ = f ′ 0 F ′ n − 2 + f ′ 1 F ′ n − 1 ,   f ′ 0 = f 0 k c k a ,   f ′ 1 = f 1 k c k b , F ′ n − 2 = F n − 2 k a ,   F ′ n − 1 = F n − 1 k b (1.20)

or identically

F ″ n ∗ = f ″ 0 F ″ n − 2 + f ″ 1 F ″ n − 1 ,   f ″ 0 = f 0 k c k a ,   f ″ 1 = f 1 k c k b , F ″ n − 2 = F n − 2 k a ,   F ″ n − 1 = F n − 1 k b , (1.21)

being k c a third dimensional factor. So the addends f ′ 0 F ′ n − 2 and f ′ 1 F ′ n − 1 have both physical dimension k c and therefore can be summed up likewise in (1.5), which yields indeed F ′ n ∗ = k c F n ∗ in this example. Since all factors just introduced are in principle arbitrary, but actually defined here for convenience unitary with appropriate physical dimensions, in fact (1.20) and (1.21) keep identically the numerical features of (1.7) although describing through k a and k b , for example, unit energy and unit momentum. The primed notations emphasize that F ′ n ∗ at the left hand side and right hand side of (1.20) are numerically identical to F n ∗ of (1.4), being in turn f ′ 0 and f ′ 1 numerically identical to f 0 and to f 1 as well; so the numerical values of the early addends f 0 F n − 2 and f 1 F n − 1 of F n ∗ are now directly referable to different physical quantities concerning the physical features of any system. The same considerations hold for (1.21).

In fact the current literature about the physical and biological implications of the standard sequence (1.1) concerns essentially contour profiling, shape matching and considerations on the importance of the golden ratio inherent its recursive rule only. All of this seems astonishing on the one hand, but also reductive on the other hand; the sequence of numbers (1.1) represents a preferential principle of Nature, quantitative formulas should be inferable likewise as from any physical law. The present paper proposes a mathematical model to overcome the gap between quantitative physical data, i.e. the acknowledged values of some physical constants, and results predictable via (1.3) through direct and detailed calculations assigning specific physical meaning to f0 and f1.

At this point however the true challenge is to demonstrate how far an appropriate choice of dimensional units really allows describing in fact a physical system; the goal is to implement explicitly F ′ n ∗ and F ″ n ∗ without changing the numerical value of F n ∗ , anyway governed by the Fibonacci recursive rule. The chance of obtaining physical information must still hold while preserving the numerical values of two generalized Fibonacci sequences merged together.

The Section 2 introduces some preliminary hints about the fine structure α; this allows to calculate a first approximation numerical value of this constant in the Section 3. The Section 4 shows that with the help of the extended Fibonacci sequence the numerical value of α can be calculated with improved approximation, one order of magnitude better. The Section 5 shows how to reveal via (1.3) with the help of (1.20) and (1.21) the subtle link between various fundamental constants, thus clarifying the concept of template itself. The importance of (1.20) and (1.21) will appear in the Section 5: owing to the numerical coincidence of f ′ 0 and f ′ 1 with f 0 and f 1 , and thus that of F ′ n ∗ with F n ∗ , the calculations will be carried out implementing (1.4) although f ′ 0 F ′ n − 2 and f ′ 1 F ′ n − 1 will actually concern dimensional physical quantities. Accordingly, the paper consists in fact of two parts: the former one describes a “standard” approach based on (1.17) to infer α , the latter one clarifies how actually the fundamental constants of Nature conform themselves to hidden templates compliant with a unique recursive rule, still that guessed by Fibonacci.

The text is organized with the main intention of making the exposition as self-contained as possible. For this reason also a few basic concepts of classical electromagnetism are shortly quoted below.

2. The Fine Structure Constant: Preliminary Considerations

Before obtaining an approximate numerical estimate of the dimensionless constant α that characterizes the electromagnetic interaction, introduce first a few well-known concepts indicating how to start from first principles. Consider preliminarily the possibility of writing

α = e 2 ℏ c = ξ ε δ r ℏ c , (2.1)

where ξ / ℏ c is the proportionality constant linking ε δ r to α . Owing to this physical meaning of α , try to identify an energy ε and length δ r considering a bound system of charges. In particular it is sensible to think ε δ r at the right hand side as product of hydrogen-like atom properties, whose analytical expressions are known and simple. As are useful short reminds about outcomes of elementary wave mechanics for a charged reduced mass m in a central field of nuclear charge Z e , introduce

ε = − Z 2 e 4 m 2 n 2 ℏ 2 ,   δ r = n 2 ℏ 2 Z m e 2 ,   m = m e A m e + A , (2.2)

where m e and A are for example electron reduced and nucleus masses. In several papers, e.g.  , these formulas have been found without solving any wave function, but as straightforward consequence of the statistical formulation of the quantum uncertainty only

δ p δ x = n ℏ = δ ε δ t , (2.3)

which in turn is a corollary itself of an operative definition ℏ G / c 2 of space time  in the frame of an evolutionary quantum model of Universe . Anyway, multiplying side by side the first two (2.2) one finds

ε δ r = − ( Z e ) e 2 , (2.4)

whence | ξ | = 2 / Z so that

α = 2 Z ε δ r ℏ c . (2.5)

The remarkable simplicity of (2.2) is the reason of having introduced just the non-relativistic hydrogen-like atom, and not for example a complex many electron atoms: the latter implies complex many electron correlations and thus the lack of analytical formulas able to include various forms of interaction, the former is implemented as electromagnetic interaction between two charges only compliant with | ( Z e ) e | defining α . This kind of interaction appears explicitly because it also follows that

| ε | = 1 2 ( α Z n ) 2 m c 2 = ( Z e ) e 2 δ r , δ r = n ℏ m c ( n α Z ) = n λ C ( n α Z ) ,   ℏ m c = λ C ; (2.6)

the systematic presence of α Z / n suggests that just α Z is the key quantity to understand α . Despite (2.2) include neither spin/orbit nor spin-spin electron/nucleus interaction or the Lamb shift, conceptually crucial even for Z = 1 although in fact numerically negligible, valuable information is in fact hidden in these equations; this suggests implementing (2.2) as an acceptable basis in defining both ε and δ r of (2.1). The factor 1/2 is the fingerprint of the quantum uncertainty: although m and nucleus are δ r apart, the Coulomb-like form requires the total uncertainty range of the charge − e with respect to Z e that indeed is 2 δ r i.e. the total radial delocalization range of the electron around the nucleus. Of course δ r is itself arbitrary and unknowable because n is arbitrary  ; thus the physical content of (2.6) is well beyond that of the mere Coulomb law. Moreover it is easy to realize with the help of (2.3) that

δ r δ t = n λ C δ t n α Z = λ C δ ε ℏ n α Z = δ ε m c n α Z = v r ,   v r = δ r δ t , (2.7)

being v r the average radial velocity of m. So, owing to (2.6),

α Z n = δ ε / m v r c = δ v r c (2.8)

having defined δ v r as δ ε / m v r by dimensional reasons; indeed

δ ε = m v r δ v r = m δ ( v r 2 ) / 2   i . e .   ε = m v r 2 / 2 + c o n s t (2.9)

is trivially the classical kinetic energy of a free particle of mass equal to the reduced mass m moving at radial velocity v r . However the non-trivial consequence of this short analysis is that the early Bohr energy can be expressed not only via Coulomb energy, as it is obvious, but also via an expression not explicitly referred to charge interaction; rather it is simply related to the rest mass m c 2 via α Z / n , because now e 4 of (2.2) becomes hidden in α 2 of (2.6). However, although ε is even expressed as classical kinetic energy of reduced mass m moving at velocity v r , the quantization is still evident through (2.7).

At this point, before implementing these results to evaluate α , note that just (2.6) suggest further chances to manipulate the early Bohr formulas to highlight explicitly the electromagnetic character of the Coulomb interaction via the electromagnetic field associated to a radial electromagnetic wave: if it is true that magnetic field H is generated by moving charges, then this H and the electric field E itself should be someway both hidden in (2.2). Also, the force Φ acting on m due to these fields should also be contextually inferable. To show these points note first of all that trivial manipulations of (2.2) yield

| ε | = Z 2 e 5 m 2 ℏ 2 2 n 2 ℏ 4 e m = E ℏ 2 c 2 e m c = E ℏ λ C c 2 e = E μ B α , E = Z 2 e 5 m 2 n 2 ℏ 4 = ( α Z n ) 2 e λ C 2 ,   μ B = e λ C 2 (2.10)

and also

e λ C μ B = 2 ,   | ε | = μ B e λ C ( α Z n ) 2 m c 2 = ( Z e ) e 2 δ r = H μ B , H = m c 2 e λ C ( α Z n ) 2 = Z e λ C δ r , (2.11)

which define the moduli E and H of electric and magnetic fields inside the atom; contextually is even defined the Bohr magneton μ B . In fact (2.11) suggest what is well known, i.e. that charges randomly moving around the nucleus induce H, which however at the atom scale is quantized itself. Multiplying side by side (2.10) and (2.11), implement thus

ε 2 = E H μ B 2 α . (2.12)

Moreover, since

μ B 2 = e 2 λ C 2 4 = e 2 λ C 3 4 λ C = 3 16 π e 2 λ C V C = 3 16 π η V C ,   V C = 4 π 3 λ C 3 ,   η = e 2 λ C , (2.13)

it follows then

ε 2 = 3 16 π E H V C η α ; (2.14)

hence, according to the first (2.6),

ε 2 η = Z 2 e 4 4 δ r 2 λ C e 2 = Z 2 λ C δ r Z e 2 2 δ r = Z 2 λ C δ r | ε | . (2.15)

So, merging (2.14) and (2.15), from

λ C δ r | ε | = 3 8 π α Z E H V C (2.16)

and owing to (2.6)

α Z n 2 = λ C δ r (2.17)

one finds

| ε | = 3 8 π ( n α Z ) 2 E H V C ,   E H | ε | = ( α Z n ) 2 2 λ C 3 . (2.18)

Note now that if the electric and magnetic fields of electromagnetic waves are orthogonal, then in fact numerically E H = | E × H | while remaining true of course that the modulus at right hand side is related to cross vector product that points towards the propagation direction of the wave. In this sense a running plane wave is also described by and equivalent to the quantized energy in the volume V C of hydrogenlike atom; i.e. photons generated by electron decay between discrete energy levels propagate as steady e.m. waves. In other words the energy density field HE of e.m. wave in V C is quantized, whereas photons are just the energy quanta of field corresponding to

δ | ε | = δ | H E V C | (2.19)

between n-th allowed levels. Note in this respect that the second couple of (2.3) reads

δ | ε | = n ℏ / δ t = n ℏ ω , (2.20)

in agreement with Planck, whereas the first couple reads δ p = n h / 2 π δ x . Then

2 π δ x = n λ   if     p = p 0 + h λ (2.21)

according to De Broglie wave definition of momentum with p 0 arbitrary constant.

So the uncertainty requires on the one hand the corpuscular nature of the charge e, see Equations (2.9) and (2.2) itself, but also wavelike energy and momentum on the other hand, see oscillating fields E and H that describe an e.m. wave propagating within V C . So the wave confined in this volume cannot be radiated without changing n, whereas photons can be however emitted/absorbed compatibly with the Pauli principle. Obviously analogous conclusion holds for any quantum system of charges bound with interaction strength constant α : given | ε | , the greater the numerical value of α , the greater the energy density EH. In conclusion, at the left hand side of (2.19) still appears the well known quantized Bohr energy, at the right hand side appears the modulus EH of the Poynting vector E × H ; the proportionality constant between the energies involves α . As the last (2.11) shows that eH has physical dimensions of force, the question at this point is whether this force can be guessed and inferred itself straightforwardly from the basic Equations (2.2). While (2.2) are well known, they have marked the birth of the old quantum mechanics, nothing in principle excludes the chance of defining also a further equation alternative to (2.4); note indeed that appears at left hand side of (2.16) the ratio | ε | / δ r . So, owing to the definitions of λ C and μ B , write with the help of second (2.2)

Φ r = | ε | δ r = 1 2 ( α Z n ) 3 m 2 c 3 n ℏ = 1 2 ( α Z n ) 3 m c 2 n λ C = m a r ,   a r = 1 2 ( α Z n ) 3 c 2 n λ C , (2.22)

where Φ r has in fact physical dimensions of force; in effect this force appears reasonably expressed as m times a new amount a r having physical dimensions of acceleration, quantized itself too. If so, then the force here acknowledged acting on m in the field of nucleus is the radial component of a more general quantized electromagnetic interaction force of modulus Φ = | Φ | between nucleus and electron at Bohr radial distance δ r = | δ r | around it, thus delocalized in an ideal sphere of diameter 2 δ r centered on the nucleus. This explains why a r results defined via c 2 over n times the length 2 λ C compliant with the quantum uncertainty. Replacing in (2.16), one finds therefore owing to (2.17)

Φ r δ r = 3 8 π α 2 E H V C . (2.23)

The fact that at the left hand side appears Φ r δ r suggests putting

Φ ⋅ δ r = ( Φ δ r ) cos φ ,   Φ r = | Φ | cos φ , (2.24)

being φ an appropriate angle. Before concerning in some more detail (2.23) and φ of (2.24), note that the definitions (2.10) and (2.11) allow inferring immediately from the radial force component Φ r further well-known concepts that confirm the validity of the present considerations. Trivial manipulations of (2.22) yield with the help of (2.10)

Φ r = 1 2 ( α Z n ) 3 m c 2 e 2 n μ B = 1 2 ( α Z n ) 3 m 2 c 2 e 2 m n μ B = 1 2 ( α Z n ) 3 ( m c ) 2 c γ n μ B           = 1 2 ( α Z n ) 3 ℏ 2 c γ n λ C 2 μ B , γ = e 2 m c , (2.25)

whence

Φ r = 1 2 ( α Z n ) 3 ( m c ) 2 γ c n μ B = 1 2 ( α Z n ) 3 M 2 c n λ C 2 j ( j + 1 ) γ μ B , M 2 = ℏ 2 j ( j + 1 ) ,   e λ C μ B = 2 (2.26)

once having introduced the square angular momentum M 2 of m around the nucleus as a function of j = l ± s ; then it is possible to write

Φ r = ( α Z n ) 3 | M | 3 c γ 2 n j ( j + 1 ) λ C 2 μ B 2 ,   1 2 = γ | M | μ B (2.27)

and thus, recalling the last (2.11) and that μ B / γ = ℏ ,

Φ r = ( α Z n ) 3 ℏ c j ( j + 1 ) n λ C 2 . (2.28)

As concerns the second (2.27), especially interesting is the particular case where | M | = | L + S | reduces to the angular momentum of spin | S | of m only, in which case the second (2.27) yields the so-called Land? factor of m; this equation reads indeed

2 = 1 γ μ B | S | ,   1 γ = 2 m c e , (2.29)

Note however that usually γ and μ B are expressed through the rest mass m e corresponding to m rather than through the reduced m itself. For example, in the case of a hydrogen-like atom with one electron, (2.10) and (2.25) yield

1 γ = 2 m c e = 2 ( m / m e ) m e c e = m m e 1 γ e , μ B = e ℏ 2 m c = e ℏ m e 2 m e m c = μ e B m e m ,   m e m = m e + A A ;

hence replacing γ and μ B so far defined with γ e and μ e B , (2.29) turns into

2 = 1 γ μ B | S | = m m e 1 γ e m e m μ e B | S |

so that

2 = 1 γ e μ e B | S | ,   1 γ e = 2 m e c e ,

as expected in non-relativistic approach. Moreover, since from (2.6)

1 λ C ( α Z n ) 2 = H e m c 2 ,

it also follows

m c 2 λ C H ( α Z n ) 2 = e ,   γ = 1 2 c λ C H ( α Z n ) 2 = 2 π ν H ,   ν = c 2 λ C ( α Z n ) 2 , (2.30)

being clearly ν a frequency. Thus, results defined

ω = H γ = e H 2 m c ,

being ν = H γ / 2 π the Larmor frequency.

Furthermore, comparing last (2.30) and first (2.6) one finds | ε | = ν λ C m c = ( ν h ) ( λ C / h ) m c , where ν λ C has physical dimensions of velocity; so being λ C = h / m c it follows that by definition | ε | = h ν whereas λ C / h = p by dimensional reasons, which are the old Planck and De Broglie relationships already identified while commenting (2.18). Note that all relations are quantized, being inferred from early (2.2) only, i.e. without any consideration about the classical definitions of angular and magnetic momenta. Moreover, dividing side by side (2.22) and (2.11), one finds

Φ r H = 1 2 ( α Z n ) e n ; (2.31)

so according to (2.30)

n α Z v λ c = 1 2 ( α Z n ) ,   v λ = ν λ C

and thus replacing in (2.31)

Φ r = H v λ e c ,   v λ = λ C ν α Z ,

being v λ the modulus of velocity v λ calculated via Compton length of m times the frequency defined in (2.30). The vector equation corresponding to this scalar relationship can be nothing else but

Φ = e v λ c × H

as stated before about (2.18), the magnetic field H of the wave along with the electric field E , must be orthogonal to its direction propagation velocity v λ . Hence, v λ × H = v λ H n , with the unit vector n that defines the propagation direction of a transversal e.m. plane wave made by orthogonal magnetic and electric oscillating fields, whose quantization implies the concept of photon.

Nevertheless this elementary approach has a further interesting implication. Define the function

G Φ r = 1 2 G ( α Z n ) 3 m 2 c 3 n ℏ = X c 4 , (2.32)

where G is the gravity constant and thus X a dimensionless function to be found; indeed the physical dimensions of the product G times force are velocity4. Since

X = G 1 2 ( α Z n ) 3 m 2 n ℏ c ,

then

X = G m 2 ϵ 0 λ ,   ϵ 0 λ = n ℏ c 1 2 ( α Z n ) 3

being λ an arbitrary length and ϵ 0 an energy defined as

ϵ 0 = n ℏ ω 0 1 2 ( α Z n ) 3 ,   ω 0 = c λ .

So

ϵ 0 = G m 1 m 2 λ ,   m 1 m 2 = m 2 X . (2.33)

Note that being λ arbitrary, are also arbitrary m 1 m 2 and thus m 1 and m 2 themselves. The interesting fact is that even the gravitational law appears to be nested in a natural way in the present elementary approach.

This result supports the validity of (1.19). Eventually calculate also

G Φ r = G m a r = ( l m ) 2 ,

where l is another arbitrary length; the last equality in inferred by dimensional reasons. Thus, owing to the second (2.22),

l = ( m G a r ) 1 / 2 ,   a r = m G l 2 = 1 2 ( α Z n ) 3 c 2 n λ C

that defines l . Consider now the particular value l b h of the length l defined as follows

l b h = 2 m G c 2 , (2.34)

which implies

2 m G c 2 = ( m G a r ) 1 / 2 ,   4 ( m G ) 2 c 4 = m G a r ,   a r ( b h ) = c 4 4 m G = 1 2 c 2 l b h ;

the form of the particular acceleration a r ( b h ) is the same as that reported in (2.22) for a r . Note that

a r ( b h ) l c 2 = l c 2 4 m G = 1 δ ϕ (2.35)

the last equality is obvious, thinking that in general v 2 / r is the classical centripetal force acting on a point mass traveling along a circular path at tangential speed v. So l δ ϕ is the length of an arc of circumference defined by the angle δ ϕ traveled by the point mass. The fact that here appears in particular v ≡ c means that this elementary interpretation holds for a photon. Thus

δ ϕ = 4 m G l c 2 (2.36)

In effect, this result is nothing else that the deflection angle δ ϕ of a photon moving in the gravity field of m at distance l , which corresponds to the famous light deflection of light beam traveling in a curved space time; the first order approximation of Einstein approach corresponds here to having approximated via (2.35) the actual photon path just along a circular arc.

This short and elementary discussion has shown that basic concepts of classical electromagnetism are hidden in the simple Bohr results (2.2), in turn inferable as straightforward and serendipitous corollaries of (2.3). Moreover this result shows once more what has been emphasized several times in previous papers    , i.e. the intimate and natural merging of quantum physics and relativity.

It is useful therefore to show shortly here how all results hitherto inferred are related to classical physics on the one hand and relativistic physics on the other hand; of course here this last aspect of the problem is shortened as much as possible, it has been more thoroughly examined in  . Nevertheless, simple and short final remarks allow extending these non-relativistic concepts simply implied by the early Bohr approach. Define indeed more in general Φ r acting on the electron in the field of nucleus as Φ is no longer required to be radial component as prospected by the second (2.24), whereas the electron delocalization around the nucleus is described by δ r cos φ . Then assuming φ = φ ( α ) , as it is natural owing to (2.24) and (2.22), and expanding in series cos φ around φ = 0 one finds recalling (2.22)

E = ε B ( 1 + ζ 1 α + ζ 2 α 2 + ⋯ ) = ε B + ⋯ ,   ε B = Φ δ r ,   ζ j = ζ j ( Z , A , m , N ) , (2.37)

being clearly ζ j the coefficients of series expansion; in this way the early Bohr energy ε B introduced in (2.2) and so far implemented results to be in fact the first order term only of a more complex energy E in which are hidden further forms of interaction. This latter statememt can be explicitly evidenced writing

E = ε B + ( ζ 1 α + ζ 2 α 2 + ζ ′ 3 α 3 + ⋯ ) ε B + ζ ″ 3 α 3 ε B + ⋯ ,   ζ 3 = ζ ′ 3 + ζ ″ 3 ; (2.38)

in this way the first addend ε B of E is the Bohr zero order term, the second addend yields the electron correlation terms as a function of α for a number of electrons N ≥ 1 , which is certainly possible in principle determining appropriately the series coefficients ζ j of (2.37). In particular the third addend ∝ α 5 is the fingerprint of the Lamb energy due to the “vacuum” polarization. Even for N = 1 appear thus expectable in principle the electron spin-orbit and spin-spin interactions missing in (QQQ). These higher order effects have been considered in   simply starting from the statistical formulation of quantum uncertainty, but will not be considered here for brevity. However it is enough noticing that rewriting (2.3) as δ ε = v δ p , with velocity modulus defined by v = δ x / δ t , and multiplying both sides by ε , one finds ε δ ε = ε v δ p . To shorten as much as possible this discussion, put now by definition p = σ ε v , where σ is a proportionality constant; this position is purposeful to obtain ε δ ε = σ ε v δ ( ε v ) that reads then

1 2 δ ( ε 2 ) = 1 2 σ δ ( ε v ) 2 ,   σ = c o n s t a n t .

Hence δ ( ε 2 ) = δ σ ( ε v ) 2 yields identically

δ ( ε 2 ) = δ ( σ ( ε v ) 2 + c o n s t ) , (2.39)

because of course δ c o n s t = 0 by definition. So the last equation reads

ε 2 = σ ( p / σ ) 2 + c o n s t .

Since by dimensional reasons σ = v e l o c i t y − 2 , this constant velocity can be nothing else but σ = c − 2 . So

ε 2 = ( p c ) 2 + c o n s t ,   p = ε v c 2 . (2.40)

Moreover implement the boundary condition for v → 0 to find const, which has physical dimensions of square energy: as

lim v → 0 p v = m = ε 0 c 2 ,

then (2.39) reads

ε 0 2 = ( m c 2 ) 2 ,   lim v → 0 ε 2 = c o n s t = ε 0 2 = ( m c 2 ) 2 . (2.41)

Clearly (2.41) and (2.40) imply the Lorentz factor 1 − v 2 / c 2 . Calculate with the help of (2.2) and (2.40) via (2.3), then

( δ ε ) 2 − ( δ p c ) 2 = n 2 ℏ 2 ( δ t ) 2 − n 2 ℏ 2 c 2 δ x 2 = ( n ℏ ) 2 ( δ x ) 2 − ( δ c t ) 2 ( δ x δ t ) 2 ;

if one implements arbitrary ranges δ ε = ε 2 − ε 1 and δ p = p 2 − p 1 at the left hand side with ε and p fulfilling (2.40), it follows that the ratio at the right hand side must be invariant. This means

( δ x ) 2 − ( c δ t ) 2 = i n v ,   δ x δ t = i n v .

Of course the same holds considering ( δ p c ) 2 − ( δ ε ) 2 that would yield ( c δ t ) 2 − ( δ x ) 2 = i n v . These invariant results are actually the signature of the quantum uncertainty; they hold not only for a free particle, but also for a particle of reduced mass m subjected to e.m. interaction. It is known that the 4D invariant interval, here symbolized by δ x 2 only instead of δ x i 2 , i = 1 , 2 , 3 for 3 space coordinates, is the conceptual basis of special relativity .

Finally implement Equation (2.8) dividing both sides by c 2 δ t and write according to (2.7)

δ ε / m c 2 v r δ t = − δ ϕ / c 2 δ r = − δ v r δ t ,   ϕ = ε m , (2.42)

so that, replacing ε = ℏ ω to describe in particular a photon in the field ϕ , one finds the relationship δ ( ℏ ω / m c 2 ) = δ ϕ / c 2 ; thus the respective energy changes yield

δ ω ω 0 = δ ϕ c 2 ,   ω 0 = m c 2 ℏ . (2.43)

Being δ ω arbitrary, define in general δ ω = ω 1 − ω 0 . Note that (2.42) has the form of one dimensional space gradient component δ ϕ / δ r of the function ϕ equal to velocity component change rate v ˙ r = δ v r / δ t , i.e. in the usual 3D form ∇ ϕ = − v ˙ , typical of the energy potential . So, regarding the left hand side of (2.43) as the energy change of a photon in the field gradient ϕ , the result

δ ω ω 0 = ϕ s − ϕ o c 2

concerns energy loss δ ω < 0 of a photon moving against the field strength gradient; the subscripts stand for source and observer respectively. This equation is in particular the red shift of a photon beam moving from the emission point of a source, ϕ s , towards the observation point, ϕ o , with decreased energy potential ϕ s < ϕ o .

These short considerations have evidenced the importance of the constant α in determining the physical properties of hydrogenlike atoms and the existence of quantized electromagnetic field required by the elementary Bohr formulas. Mostly important, this elementary analysis of the e.m. interaction shows that there is no conflict between QM and relativity; rather, trivial algebraic steps are enough to bridge both theories.

With this introductory background in mind, it is possible to tackle self-consistently in the next sections the problem of determining the numerical value of α .

3. First Approximate Calculation of α

After having implemented (2.2) to emphasize a few basic concepts of electromagnetism, is reasonable the attempt to implement also Φ r to get a preliminary estimate of α . Note that the Equations (2.22) contain explicitly m, whose value determines whether these equations describe specifically hydrogenlike atoms or mesic atoms or even a proton/antiproton system; thus it is convenient in principle to rewrite (2.22) in dimensionless form

Φ r * = Φ r n ℏ m 2 c 3 = Φ r 2 π n h m 2 c 3 = n Φ r π λ C 2 m c 2 = 1 2 ( α Z n ) 3 (3.1)

in order that Φ r * holds for any bound system governed by electromagnetic interaction with central charge Z e regardless of the particular reduced mass m of the charge e. First of all (3.1) suggests identifying at the left hand side of the last equality the amount ( Φ r / 2 ) ( h / m 2 c 3 ) n such that

n Φ r m a r = n Φ r Φ r = n = π 2 ( α Z n ) 3 ,   Φ r = m a r ,   a r = 2 c 2 λ C , (3.2)

where clearly a r is the radial acceleration of m due to the force Φ r . So comparing with (2.22)

a r = 2 c 2 λ C = 1 2 ( α Z n ) 3 c 2 n λ C = π ( α Z n ) 3 c 2 λ C

i.e.

n 4 2 π = ( α Z ) 3 ,   α = α ( Z , n ) .

In this way, α is a function of n and Z. To define uniquely the value of α Z is interesting to calculate this result in particular for the fundamental state n = 1 and Z = Z max , which actually calculates the minimum value allowed for α . Introducing thus the boundary condition n = 1 and Z = Z max that defines Φ max * = Φ r * ( n = 1 , Z = Z max ) , the last equation reads

α = ( 2 π ) 1 / 3 1 Z max . (3.3)

This result can be calculated recalling that in  it has been found that the upper limit of high atomic number heavy nuclei is Z max = 118 . Then

α = ( 2 π ) 1 / 3 1 118 = 0.0072902883 ,   d e v = − 0.0968 % . (3.4)

The deviation of α calculated in this way from the true value (1.18) is due to the non-relativistic approximation of the Equations (2.2); even so, however, the approach hitherto followed appears basically correct.

The agreement provided by this kind of preliminary approach, although carried out through classical concepts of acceleration and force, explains the stability of the quantum system of interacting charges. The definition of acceleration in (3.2) reminds formally that v 2 / r of the circular motion of m tethered around a fixed point, which however before the birth of QM has posed the problem of explaining why there is no radiation loss in the case of a bound system of a charge m around a fixed nucleus. The reason is that classically, being in general v = v ( t ) and r = r ( t ) , the emission of radiation is allowed to occur contextually to the spiral motion of the charge m approaching closer and closer to the nucleus till to the collapse of the system. However the definition (3.2) of a r prevents such behavior, because both c and λ C are constant quantities: here instead no change of c and λ C , and then no collapse. Hence the atomic system does not radiate because m cannot spiral towards the nucleus; the Bohr postulate is actually the constancy of light speed in vacuum. Nevertheless the system can still emit discrete amounts of energy, we call them photons, changing n by integer amounts.

Anyway, a much more effective approach to obtain a better value of α is next carried out implementing the Fibonacci sequence.

4. The Generalized Fibonacci Sequence

Often an approximation is useful to check immediately and easily whether or not any reasoning or algebraic steps point to the right direction in calculating something, provided that however the basic requirements of the calculation model are still valid; for this reason this section starts just from (1.16) and (1.17), where σ 1 > 0 and R > 0 by definition whereas σ ′ 1 can be in principle positive or negative.

When switching from (1.1) to (1.3), two considerations arise about f 0 and f 1 replacing 1. The first problem concerns the chances f 0 > 1 and f 1 > 1 or 0 < f 0 < 1 and 0 < f 1 < 1 , both possible in principle. The second problem is that just the most general meaning of (1.3) requires determining two initial unknowns, a sort of boundary condition without which the sequence is in fact undefined. These requirements are compliant with (1.10) regarding either σ ′ j > 0 and σ ″ j > 0 or σ ′ j < 0 and σ ″ j < 0 in the definitions (1.12). Actually both requirements are fulfilled in the simplified form (1.15) at the first order only, where σ ′ 1 = σ ″ 1 ; even neglecting the higher order coefficients, it is enough to regard either σ ′ 1 > 0 or σ ′ 1 < 0 , in the following indicated with the notation − | σ ′ 1 | . For σ ′ 1 > 0 holds the upper sign only, i.e. the first (1.17) reads

α + o = − 1 + 1 + 4 σ ′ 1 R / σ 1 2 σ ′ 1 , (4.1)

whereas for σ ′ 1 < 0 both signs of the solutions are admissible but with the following condition

α ± o = − σ 1 ± σ 1 2 − 4 | σ ′ 1 | σ 1 R − 2 | σ ′ 1 | σ 1 = 1 ∓ 1 − 4 R | σ ′ 1 | / σ 1 2 | σ ′ 1 | ,   4 | σ ′ 1 | σ 1 R ≤ 1. (4.2)

This last case is interesting as

α + o + α − o = 1 | σ ′ 1 | (4.3)

i.e. if | σ ′ 1 | = 1 , then α + o and α − o take probabilistic meaning whose sum yields the certainty. This means that F n * * * has probabilistic meaning, as it can take here two values depending on whether α ∘ takes either value α + ∘ or α − ∘ . Clearly (1.17) admits in general several values of α ∘ depending on the number of series expansion terms in (1.14); yet it is reasonable to think that a probabilistic equation analogous to (4.3) still holds for the sum of all real solutions α j ∘ . Anyway, since | σ ′ 1 | = 1 with σ ′ 1 < 0 means σ ′ 1 = − 1 , now (1.15) reads

F n * * * σ 1 = ( F n − 1 + F n − 2 ) ( 1 − α ∘ ) α ∘ . (4.4)

As (3.4) implies that α must depend on Z max − 1 , it is reasonable to expect also now

α ∘ = c o n s t Z max ,   F * * * σ 1 ≈ Z max (4.5)

so that the solution of (4.4) with respect to a ∘ reads

α ± o = 1 ∓ 1 − 4 Z max / ( F n − 1 + F n − 2 ) 2 ,   4 Z max ≤ F n − 1 + F n − 2 . (4.6)

Since by definition two known addends F n − 1 + F n − 2 concur to the n-th term of the standard (1.1), which is thus known, this result is calculable as a function of Z max . If the present reasoning is correct, it should be true that

α = α ± o Z max = 1 ∓ 1 − 4 Z max / ( F n − 1 + F n − 2 ) 2 Z max ; (4.7)

in other words, α ± o replace the factor ( 2 / π ) 1 / 3 of (3.4). Yet the interesting fact is that now not only α ± o depends itself upon Z max but also that it can take two values. To assess this conclusion according to the condition (4.6), it is enough to put Z max = 118 , i.e. 4 Z max = 472 , and check which n of (1.1) fulfills

1 ± 1 − 472 / ( F n − 1 + F n − 2 ) 2 ≈ ( 2 π ) 1 / 3 ,   F n − 1 + F n − 2 ≥ 472 (4.8)

if the reasoning is correct, then an appropriate value of n must exist such that either sign should provide a better approximation to the true value of α . A few results for selected values of n are reported in following Table 1 that for convenience of comparison also reminds the reference value (4.8) early introduced in (3.3) to calculate (3.4).

( 2 π ) 1 / 3 = 0.86025 ⋯

It appears that in effect for n = 16 the expression (4.8) with plus sign fulfills unambiguously the given similarity condition α + o ≈ ( 2 / π ) 1 / 3 , which hopefully should provide the value of α + o Z max of interest to calculate according to (4.7) a better value of α . Indeed

α c a l c = α + o | n = 16 Z max = 0.007298074622 ,   d e v = 0.0098 % (4.9)

now the deviation from (1.18) is one order of magnitude better than that of (3.4). Apart from the enhanced numerical agreement of this result, a crucial question arises now: is it possible to improve further the accuracy of the calculations? From the mere numerical point of view, certainly the implementation of (1.14) instead of (1.15) affects the final approximation inherent the value (4.9); so, it is in principle reasonable to expect that the calculation of α ± o including in (1.17) higher order terms prospected by (1.14) would bring to a result even better than (4.9).

The calculations proposed in this section aimed merely to provide a first answer to the challenge proposed at the end of the Section 1, i.e.: the Fibonacci sequence appears in fact adequate to calculate an acceptable value of α if appropriately implemented. Nevertheless the next section will describe in this respect a much better and far reaching approach by following a completely different strategy.

5. Calculation of Fundamental Constants via (1.4)

Consider again (1.4) with the purpose of correlating the fundamental constants of nature, which are now assumed all known: this section aims indeed to find the possible interconnection between these constants.

Calculation of a + ∘ and a − ∘ according to the Equation (4.7) for various values of n. The solutions for n < 15 are imaginary
nα + oα − o
150.73781794930.2621820507
160.86117280540.1388271946
170.91965635820.0803436418
180.95203358880.0479664112

To this purpose implement in particular Avogadro’s number N A and ratio r p e of proton to electron rest masses m p / m e , here reported for completeness along with α itself quoted in (1.18)

r p e = m p m e = 1836.15267343 ,   N A = 6.02214076 × 10 23 . (5.1)

It is possible in principle to calculate f 0 and f 1 for various test values of n by defining appropriately F n * , for example through the positions

f 1 − f 0 n = α ,   F n * N A = 3 ,   F n * = f 1 F n − 1 + f 0 F n − 2 , (5.2)

and solving this system of two equations with respect to f 0 and f 1 . The result of this calculation that specifies the respective F n * ( n ) is reported in Table 2 for selected n only, the ones of interest for the next considerations.

Regarding n as arbitrary numerical parameter, this calculation has actual physical worth only if it someway identifies among the various f 0 ( n ) and f 1 ( n ) a particular value n * uniquely related to some specific quantity already introduced, like for example α + o and α − o of Table 1; if so, then (5.2) not only define a relationship linking the values of r p e and N A to α but also reveal how this n * is in turn related itself via f 0 = f 0 ( n * ) and f 1 = f 1 ( n * ) to another fundamental constant, i.e. to Z max .

This premise outlines the strategy of the present approach: instead of improving further (3.3) and (4.9) to calculate α only, which would be conceptually restricted and futile, are now implemented two further constants of Nature, N A and r p e ,to understand how their mutual correlation defines uniquely F n * * and how even other constants are possibly hidden in this F n * * along with α too.

Regarding in fact (5.2) as boundary conditions for α and N A , it follows that the index n is upper bound in a natural way by n max = 118 , beyond which appear negative values of f 0 . But these negative values are to be excluded, as they would imply F n * = f 1 F n − 1 − | f 0 | F n − 2 , i.e. F n * would no longer correspond to the sum of two previous values f 1 F n − 1 and f 0 F n − 2 . To comply with the original definition of sequence (1.3), therefore, identify n * ≡ n max as upper limit of n fulfilling (5.2) and (3.4), and thus n max ≡ Z max . In other words, the upper boundary of Z allowed to the heaviest nucleus of hydrogen-like atom is related to the meaning of n max -th term for which the generalized sequence (1.4) is self-consistently satisfied. So, if this conclusion is true, the numbers calculable through f 0 and f 1 of Table 2 corresponding to n = Z max should have their own physical meaning; indeed, accepting the positions (5.2), Z max enters in a natural way into the model via the upper boundary allowed to n.

Consider separately the addends f 0 ( n ) / n and f 1 ( n ) / n of the first (5.2) that by definition concur to calculate α . This means considering in particular the ratios

f 1 Z max = 0.007410575579 ,   f 0 Z max = 0.0001132 ,   n = n max = Z max (5.3)

Calculation of f 0 and f 1 solutions of (5.2) for selected values of n
nf 1f 0
1152.6314959911.792300446
1161.7515782840.9050853867
1171.2088239040.3550336533
1180.87444791830.01336031542
1190.6688568636−0.1995280918

that extend further the approximate (3.3) and (4.7) successfully implemented together to calculate α , here assumed instead exactly known. These initial positions are introduced thinking that in fact nothing excludes their own physical meaning even though regarded separately. To highlight the possible links of (5.3) with other fundamental constants of nature, exploit the quantities of Table 2 recalling also both (1.20) and (1.21)

f ′ 1 = f 1 k c k b ,   f ″ 1 = f 1 k c k b ,   f ′ 0 = f 0 k c k b ,   f ″ 0 = f 0 k c k b

since k a and k b and k c are dimensional unit factors, anyway the numerical values of Table 2 remain in fact unchanged. According (1.20) these multiplicative factors simply switch f 0 and f 1 to f ′ 0 and f ′ 1 and thus assign physical dimensions to F n − 1 and F n − 2 that switch to F ′ n − 1 and F ′ n − 2 . Although F n ∗ takes by consequence physical dimensions F ′ n ∗ , nothing changes from a mere numerical point of view as concerns Table 2 once considering (1.20) instead of (1.4). Analogous reasoning holds for the corresponding double primed quantities of (1.21).

To highlight the worth of these considerations about (5.2), let us examine the implications of Table 2.

1) The first implication is that just mentioned, i.e. n max = Z max indicates the maximum atomic number Z max = 118 expectable for the reasons previously exposed in the Section 3.

2) Define the following ratios directly inferred from (5.3).

r 1 p = m p m e f ′ 1 Z max = 13.60694816 k c k b , r 1 d = m p m e 1 f ″ 1 Z max = 251619.0161 k c k b (5.4)

and analogously

r 0 p = m p m e f ′ 0 Z max = 0.2078947362 k c k a , r 0 d = m p m e 1 f ″ 0 Z max = 1.621713325 × 10 7 k c k a (5.5)

indeed there is no reason to exclude that the analytical form guessed in (5.4) for f 1 should be also allowed as in (5.5) for f 0 . The ratio r p e = m p / m e is regarded in both cases as proportionality constant, purposely introduced to include and link in the following pattern of calculations also the proton to electron mass ratio. Clearly the idea is to implement fundamental constants only, and not arbitrary proportionality constants possibly regarded as additional best fit values.

Eventually examine also

f 0 f 1 f ′ 1 Z max = 0.865770134 × 10 − 4 k c k b , f 0 f 1 f ′ 0 Z max = 1.322773129 × 10 − 6 k c k a , (5.6)

where f 0 f 1 = f 0 f 1 ( Z max ) is another proportionality constant compliant with the leading ideas of (5.4) and (5.5). Eventually (5.6) have been also added to the list as further examples of a different proportionality factor, f 1 ( n max ) f 0 ( n max ) , appropriately replacing m p / m e for a further check of results.

So the positions (5.4), (5.5) and (5.6) that simply add a proportionality factor to (5.3) according to (1.20) and (1.21), should be reasonably admissible taking into account the idea just exposed: i.e. any expression involving f 1 ( n max ) cannot exclude f 0 ( n max ) of Table 2 too.

3) Start from (5.4); a glance to the first numerical value suggests putting purposely k c / k b → u n i t   e n e r g y and ( k b k c ) − 1 → u n i t   e n e r g y too, as it is in fact possible. Multiplying side by side these correspondences one finds 1 / k b 2 = ( u n i t   e n e r g y ) 2 i.e.

k b − 1 = ± u n i t   e n e r g y (5.7)

whatever k c might be.

Put then k c = 1 in order that both r 1 p and r 1 d take the meaning of energies whose numerical values coincide of course with ( m p / m e ) ( f 1 / Z max ) and ( m p / m e ) / ( f 1 / Z max ) . In other words (5.4), although calculated as pure numbers solving (5.2), with this particular choice of unit dimensional factors correspond in fact to two different values of energy allowed by the generalized Fibonacci sequence.

The importance of this conclusion appear recalling the Dirac equation of hydrogen like atoms

E D m e c 2 = 1 1 + α 2 Z 2 ( n − ( j + 1 / 2 ) + ( j + 1 / 2 ) 2 − α 2 Z 2 ) 2 − 1 , m e c 2 = 0. 51099895 × 10 6 eV ; (5.8)

once calculated for n = 1 , j = l + s = 1 / 2 , this expression yields for Z = 1 and Z = 118 respectively

E D ( Z = 1 ) = 13.60601134   eV ,   E D ( Z max ) = 251178.148   eV .   (5.9)

Now it is possible to compare (5.4) and (5.9): it is immediate to acknowledge that the energy E D ( Z = 1 ) corresponds to the template value ( m p / m e ) ( f ′ 1 / Z max ) reported in the first (5.4) with numerical deviation 0.007% only. If this idea is sensible, analogous conclusion should hold also for E D ( Z = Z max ) : i.e. the upper boundary of numerical value of Dirac energy should fit its template given by the second (5.4). In effect the second (5.4) results fairly consistent with the idea that ( m p / m e ) ( f ″ 1 / Z max ) − 1 is the energy to which corresponds E D ( Z = Z max ) ; now however the Dirac value deviates from its numerical template by 1.4%. It must be reminded in this respect that the early Dirac equation quoted here is inherently approximate, as it does not account for the Lamb shift, which increases with Z; in effect this term appears in a natural way in a modified Dirac equation inferred in . So it is not surprising that the deviation of E D ( Z = 1 ) from ( m p / m e ) ( f ′ 1 / Z max ) − 1 is better than that of E D ( Z = Z max ) from ( m p / m e ) ( f ″ 1 / Z max ) . Hence, regardless of the physical background that brings up to (5.8) starting from first principles, it seems that both boundary values (5.9) are coherent with and in fact already inherent (5.4).

It is also worth noticing in this respect that the double sign allowed in (5.7) is consistent with the negative energy states of Dirac too.

As anticipated in Section 1, it appears here that that F n * ( n max ) is a sort of numerical template, such that the dimensionless values from it calculated via f 0 and f 1 represent a constrain to the actual fundamental constants governing observable physical properties of Nature like the generalized Dirac energy.

Let us verify this idea to remove the suspect that the prospective Fibonacci templates of energy fit accidentally only these favorable checks looking for further correspondences similarly as that just proposed. In other words, let us check whether or not characteristic energies or lengths or times describing natural events or properties of matter are actually allowed values relatable to and in line with distinctive pure numbers inferred from the generalized Fibonacci sequence.

To identify significant templates that govern physical constants and/or their combinations is in principle a difficult task. The fact that the values calculable via (5.4), (5.5) and (5.6) are actually abstract numbers, does not indicate “a priori” which specific physical properties they refer to; only “a posteriori” one acknowledges, by comparison like that of (5.4) and (5.9), which allowed physical constants and properties tentatively found are really correlatable. Moreover are crucial in this respect also the possible physical corollaries like (5.7) contextual to the numerical comparisons. Additional matches are shown in the following.

Consider the next template numbers (5.5) to find their possible corresponding fundamental constants and note that

3 e c = 0.2083384819   s 2 cm 5 / 2 ⋅ g ,   G h = 1.613785378 × 10 7 cm / s g 2 (5.10)

having expressed h as eVs. The dimensional unit factors k c / k a and k c k a , now appropriately regarded in (5.5) in order to match the ones in (5.10), allow comparing the respective numerical values. It is immediate to observe that these acknowledged values deviate from the corresponding two (5.5) by 0.2% and 0.5% respectively. It is worth noticing that the numerical agreement of the first (5.10) with its own numerical template (5.5) requires c and e/3.

Consider eventually, by analogy with (5.4) and (5.5), the positions (5.6) involving as further proportionality constant f 0 f 1 alternative to r p e already tested. Is crucial also now the fact that the physical dimensions of k c and k a separately allow defining new dimensions k c / k a and k c k a consistent with further that of fundamental constants of Nature. Quote then the Boltzmann constant with energy expressed in eV as well

k B = 0 .8617333262145 × 10 − 4 eV / K (5.11)

then k B identifies now the dimensional ratio k c / k a , whereas the numerical deviation of k B from the first template (5.6) is 0.5% only. Moreover calculate

3 e r p e 2 = α r p e ℏ c 2 e / 3 = 1.322912454 × 10 − 6 g 1 / 2 ⋅ cm 3 / 2 s (5.12)

this value differs from the second template (5.6) by 0.01% only, once having assumed the physical dimensions of k c k a consistent with that at the right hand side of (5.12). This agreement requires now 2 e / 3 .

Note now that the first Equation (1.19) can be better understood and enhanced comparing once more the link of G to e ratio via α through its pertinent template

G α e = 1.013971342   cm 3 / 2 g 3 / 2 ⋅ s ,   f ′ 1 α n max = f 1 α Z max k c k b = 1.015515629 k c k b ; (5.13)

the numerical deviation of the former from its own template is 0.15% only. It is significant that one calculates here just the numerical proportionality factor tentatively estimated in (1.19) with the known values of G and e / α , whereas the dimensional matching is fulfilled once more by k c / k b .

Emphasize also now the physical implications of these numerical comparisons.

On the one hand the last results (5.10), (5.12) and (5.13) provide a list of “effective charges” allowed by the respective templates

e ,   2 3 e ,   1 3 e . (5.14)

On the other hand, the link between G and e via α intuitively guessed in (1.19) appears actually to be the further result of a wider conceptual frame that according to (5.2) assigns also to f 1 / n α = f 1 / ( f 1 − f 0 ) the physical meaning of template value, thus extending the Fibonacci driven list of hidden links that in fact represent the constants of Nature.

In summary: the templates (5.4), (5.5) and (5.6) has been related respectively to (5.9), (5.5) and (5.11), (5.12).

6. Discussion

For sake of clarity the present paper has deliberately proposed a stepwise exposition about the fundamental constants of Nature; the purpose of this strategy was to show first in the Sections 3 and 4 that (1.12) and (1.6) account decently for the numerical definition of α and prospect in principle even better results simply depending on the number of higher order terms of series expansion that define the parameters of (1.11).

Next the considerations of Section 5 have explored the idea that the visible look of the Universe is controlled by a set of fundamental constants that organize its observable features analogously to the formal aspect of a typewritten sheet subjected to a pattern of rules coded in a basic template. This holds if in fact exists a subtle link underlying the numerical values of the various constants that take physical meaning via their physical dimensions; i.e. physical dimensions and numerical values cannot be regarded as separate properties. In other words, the gravity constant or the electric charge have the value G and e we know coherently with the respective physical dimensions; the link (1.19) confirmed in (5.13) must be regarded in this way, which justifies the dimensional factor k g in a natural way and a wider conceptual frame. In effect just the dimensional requirements allow the correspondences between physical constants and their own templates via the unit factors k a and k b : as shown by (1.20) and (1.21), the sequence (1.4) switches from (1.6) to k c F n * likewise as shown for example in (1.9).

The numerical agreements of Section 5 can be hardly interpreted as mere accidental coincidences.

In other words: considering “exact” by definition the known values of N A and α implemented to calculate the coefficients f 0 and f 1 of the generalized Fibonacci sequence (1.3) as a function of n, the fundamental constant m p / m e results to be itself the proportionality factor enabling not only to calculate reasonably the ground states (5.9) of the Dirac hydrogen-like atom at Z = 1 and Z = Z max but also to obtain the further correspondences (5.10) with (5.5). It is significant the fact that after having introduced (5.3) the proportionality factor defining (5.4) and (5.5) is just one among the fundamental constants of Nature and not an arbitrary numerical parameter that would appear inevitably as a best fit value to match purposely some known physical constant. The way to define (5.4) and (5.5), instead, excludes any best fit attempt; rather all correspondences found in Section 5 are self contained, in that the template model plugs physical constants into the extended Fibonacci sequence via (5.3) only. Indeed the same holds identically for the positions (5.6) that implement the factor f 0 f 1 and not some superfluous proportionality parameter to be intentionally defined.

For this reason has been omitted the possible explanation about why just the factor m p / m e links reasonably f 1 / Z max to E D ( Z = 1 ) and E D ( Z = Z max ) : if so, then the boundary conditions at Z = Z max and Z = 1 of the generalized Fibonacci sequence are accordingly related even regardless of the theoretical approach leading to the Dirac equation.

Anyway, apart from the chance of justifying rationally why the energies (5.9) conform to their templates (5.4), if is correct the present idea of regarding (1.2) at n = n max as numerical template to which the properties of Nature conform likewise a boundary condition, the approach just described links and correlates in a unique frame the ratio m p / m e , the upper and lower limits of Dirac relativistic ground energy, the Avogadro number and α itself along with the value of Z max ; indeed the latter two have defined f 1 at n = Z max . It is as if the generalized Fibonacci sequence “would know” in advance that the standard Dirac equation, as such, is not fully adequate to describe the reality although being an analytical expression very close to the reality; in other words (5.4), despite its mathematical abstractness, seems suggesting that something is missing in E D to fit correctly the numerical template of Nature.

The last section has shown that selecting appropriately case by case the dimensional unit factors introduced in (1.20) and (1.21), the numerical values calculable through (1.4) take physical meaning of further fundamental constants actually governing the Nature.

The practical help of the results in Section 5 is not limited to the determination of some fundamental constants, rather it implies the knowledge of further values having the rank of fundamental constants themselves. A typical example is the second (5.9): even without knowing the result inferred in , the significant deviation of ED (Zmax) from its template value required by the second (5.4) reveals that something is missing in the Dirac equation. In this specific case, the discrepancy has a known origin, the vacuum polarization causing the Lamb shift; in a more general case, however, the values calculated via the extended sequence (1.3) set end points of calculations useful to refine and adjust the theoretical models.

Last but not the least: the few remarks in Section 2, although very shortly sketched, are enough to show once again in a straightforward and elementary way the intimate link between relativity and quantum physics.

7. Conclusion

The sequence (1.1) is defined by its own sum rule of terms only, the sequence (1.3) instead does not because are also required the initial values; it needs defining also the initial values of two arbitrary parameters f 0 and f 1 . In this sense the sequence F n * (1.4) is closer to the rational formulation of physical problems, whose solution requires indeed specifying appropriate boundary conditions. Thus it is not surprising the fact that the abstractness of pure numbers (1.1) points via (1.4) directly to the reality of natural events, regardless of how the fundamental constants enter into specific formulas calculable and verifiable by comparison with the experience. This assertion is probably the best way to summarize the concept of template hitherto concerned, because it is equivalent to say: the fundamental constants shape the observable Universe regardless of how the theoretical models infer specific formulas likewise as a template controls how a typewritten sheet appears regardless of its specific text content.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

Cite this paper

Tosto, S. (2019) Generalized Fibonacci Sequence: Possible Template for the Constants of Nature. Journal of Applied Mathematics and Physics, 7, 3037-3066. https://doi.org/10.4236/jamp.2019.712214

Appendix A

The purpose of this section is to sketch how (2.3) allow inferring (2.2) along with j of (2.26), thus making really self-contained the results of present paper.

1) Angular momentum component. Since by definition M = r × p , write its classical form and its radial component as δ M = δ r × δ p and δ M k = δ r × δ p ⋅ k : the left hand sides emphasize the range sizes resulting from the range of values of r and p within the respective range sizes, k is an arbitrary unit vector. Trivial manipulations yield δ M k = k × δ r ⋅ δ p , so that calling δ r ⊥ = k × δ r one finds δ M k = δ r ⊥ ⋅ δ p . Hence with δ p ⊥ = δ r w ⋅ δ p / | δ r w | the result is δ M k = δ r ⊥ δ p ⊥ . Therefore (2.3), which must be anyway fulfilled whatever the concerned ranges might specifically be, requires δ M k = l ℏ ; with usual notation l is the integer including 0 corresponding to the product of ranges δ r ⊥ δ p ⊥ whose sizes, whatever they might be, must be anyway regarded according to the initial (2.2). Now δ M k is the uncertainty range of values of angular momentum components M k = l ℏ corresponding to the arbitrary values of l.

2) Angular momentum. Of course it is evident that one component only of M k is definable, usually called M z : changing k does not add new physical information to δ M k . Yet it is shown in  that it is possible to calculate

M 2 = 〈 M x 2 〉 + 〈 M y 2 〉 + 〈 M z 2 〉 = ℏ 2 ∑ i = 1 3 l i 2 2 L + 1 = L ( L + 1 ) ℏ 2 ,   − L ≤ l i ≤ L

and that the allowed l i concurring to the quantum angular number L are just the values of magnetic quantum number usually denoted with m.

3) Energy. The classical energy ε of hydrogen like atom ε = ε c m + p r 2 / 2 m e + M 2 / 2 m e r 2 − Z e 2 / r , where ε c m is the center of mass energy of the atom as a whole is rewritten in order to emphasize that actually r must be intended as δ r , i.e. a range of distances of the electron from the nucleus rather than a deterministic classical distance. It is immediate to find that δ ε = ε c m + δ p r 2 / 2 m e + δ M 2 / 2 m e δ r 2 − Z e 2 / δ r yields just both (2.2) with the help of (2.3) only; the key step is that here δ p r = n ℏ / δ r , whereas n appearing explicitly itself in the final form of ε is the principal quantum number.

4) Spin. Note the identity L ( L + 1 ) ≡ ( L + s ) 2 − s 2 with s = 1 / 2 ; so M 2 = ( L + s ) 2 ℏ 2 − ℏ 2 s 2 defines a new angular momentum M 2 = M 2 + ℏ 2 s 2 + ℏ 2 ( L + s ) = ℏ 2 ( L + s ) 2 + ℏ 2 ( L + s ) i.e. M 2 = ℏ 2 j ( j + 1 ) , with j = L + s . So, the quantization implies itself the chance of the spin. Actually the existence of such half-integer angular momentum is shown starting from first principles in other papers, while being in fact j = L ± s  .

These short remarks aim merely to sketch the chance of justifying (2.2) and (5.8) via (2.3) only, in order to make the present paper as self contained as possible.

ReferencesTosto, S. (1996) An Analysis of States in the Phase Space: The Energy Levels of Quantum Systems. Il Nuovo Cimento B, 111, 193-215. https://doi.org/10.1007/BF02724645Ferrando, V., Calatayud, A., Andres, P., Torroba, R., Furlan, W.D. and Monsoriu, J.A. (2014) Imaging Properties of Kinoform Fibonacci Lenses. IEEE Photonics Journal, 6, Article ID: 6500106. https://doi.org/10.1109/JPHOT.2014.2304560Li, C.R., Zhang, X.N. and Cao, Z.X. (2005) Triangular and Fibonacci Number Patterns Driven by Stress on Core/Shell Microstructures. Science, 309, 909-911. https://doi.org/10.1126/science.1113412Linage, G., Montoya, F., Sarmiento, A., Showalter, K. and Parmananda, P. (2006) Fibonacci Order in the Period-Doubling Cascade to Chaos. Physics Letters A, 359, 638-639. https://doi.org/10.1016/j.physleta.2006.07.036Mac, N., Jagannathan, A. and Pichon, F. (2016) Fractal Dimensions of the Wave Functions and Local Spectral Measures on the Fibonacci Chain. Physical Review B, 93, Article ID: 205153. https://doi.org/10.1103/PhysRevB.93.205153Dai, H.T., Liu, Y.J. and Sun, X.W. (2012) The Focusing Property of the Spiral Fibonacci Zone Plate. Proceedings of SPIE, 8257, 27. https://doi.org/10.1117/12.908398Mitchison, G.J. (1977) Phyllotaxis and the Fibonacci Series. Science, 196, 270-275. https://doi.org/10.1126/science.196.4287.270Carbone, A., Gromov, M. and Prusinkiewicz, P. (2000) Pattern Formation in Biology, Vision and Dynamics. World Scientific, Singapore, 78. https://doi.org/10.1142/4043Pletser, V. (2017) Fibonacci Numbers and the Golden Ratio in Biology, Physics, Astrophysics, Chemistry and Technology: A Non-Exhaustive Review.Tosto, S. (1996) An Analysis of States in the Phase Space: The Diatomic Molecules. Il Nuovo Cimento D, 18, 1363-1394. https://doi.org/10.1007/BF02453780Tosto, S. (2016) Quantum and Relativistic Corollaries of an Operative Definition of Space Time. Journal of Advances in Physics, 11, 3408. https://doi.org/10.24297/jap.v11i5.362Tosto, S. (2018) Evolution, Quantization, Relativity: An Ab Initio Model. Journal of Modern Physics, 9, 2495-2593. https://doi.org/10.4236/jmp.2018.914161Tosto, S. (2011) An Analysis of States in the Phase Space: From Quantum Mechanics to General Relativity.Landau, L. and Lifchitz, E. (1966) Theorie du Champ. Ed Mir, Moscow. (In French)Tosto, S. (2009) An Analysis of States in the Phase Space: The Specific Heat of Metals. In: Pandalai, S.G., Ed., Recent Research Developments in Physics, Vol. 8, Transworld Research Network, Trivandrum, 7-92.Tosto, S. (2000) A Heuristic Approach to the Relativistic Angular Momentum.