WJCMP  Vol.10 No.3 , August 2020
Phase Transitions Governed by the Fifth Power of the Golden Mean and Beyond
Abstract: In this contribution results from different disciplines of science were compared to show their intimate interweaving with each other having in common the golden ratio φ respectively its fifth power φ5. The research fields cover model calculations of statistical physics associated with phase transitions, the quantum probability of two particles, new physics of everything suggested by the information relativity theory (IRT) including explanations of cosmological relevance, the ε-infinity theory, superconductivity, and the Tammes problem of the largest diameter of N non-overlapping circles on the surface of a sphere with its connection to viral morphology and crystallography. Finally, Fibonacci anyons proposed for topological quantum computation (TQC) were briefly described in comparison to the recently formulated reverse Fibonacci approach using the Janičko number sequence. An architecture applicable for a quantum computer is proposed consisting of 13-step twisted microtubules similar to tubulin microtubules of living matter. Most topics point to the omnipresence of the golden mean as the numerical dominator of our world.

1. Introduction

The golden mean or golden ratio is an omnipresent number in nature, found in the architecture of living creatures as well as human buildings, music, finance, medicine, philosophy, and of course in physics and mathematics including quantum computation [1] [2] [3]. It is the most irrational number known with the simplest continued fraction representation at all [4] and a number-theoretical chameleon with a self-similarity property. All these properties render it to be suitable for quantum computer application. According to Olson et al. [2], the whole universe functions as a golden supercomputer. The seminal ε-infinity approach developed by El Naschie years ago and applied to cosmological questions based on this simple principle of nature [5]. El Naschie and his scientific fellows condensed these ideas to a grand unification of the sciences, arts and consciousness [2]. Before also the great women scientists Mae-Wan Ho [6] and Leila Marek-Crnjac [2] [7] made valuable contributions to this topic as well as Ji-Huan He [2]. Another contribution of Klee Irwin about unification of physics and number theory takes a similar approach using statements from quantum physics [8].

However, recently the unification of physics was advanced by Suleiman’s information relativity theory (IRT) correcting an overlooked flaw in Newton’s theory with respect to time displacements between observer and moving bodies [9]. Many formal physical constructs are overcome by this new insight.

Whereas the unification of physics is intimately connected with the golden number system, in this contribution the main focus is simply on phase transformation governed by the fifth power of the golden mean and insights beyond this topic.

Herein the golden mean is denoted by φ. Its inverse number is just φ 1 = 1 + φ . One yields numerically

φ = 5 1 2 = 0.6180339887 (1)

φ 1 = 5 + 1 2 = 1 + φ = 1.6180339887 (2)

It is important to notice that in numerous papers, for instance such about

topological quantum computation (TQC), the golden mean is taken as 5 + 1 2 .

However, in emerging matrix representations for Fibonacci anyon models actually the inverse of this number is always found in the matrix elements. So the choice of the golden ratio according to Equation (1) is well founded and should be maintained.

While in nature this number emerges as result of chaotic-statistical evolution, on the other hand, the fifth power of φ is frequently linked to phase transitions from microscopic to cosmic scale. The secrete behind both number may be revealed by considering their infinitely continued fraction representations [4] indicating how simple nature’s most effective evolutionary forecast works

φ = 5 1 2 = 1 1 + 1 1 + 1 1 + (3)

φ 5 = 1 11 + 1 11 + 1 11 + (4)

Rogers-Ramanujan presented a continued fraction of the form [4] [10]

r ( τ ) = q 1 / 5 1 + q 1 + q 2 1 + q 3 1 + (5)

where q = e 2 π i τ . For τ = 0 one gets the golden mean φ, and for τ = i = 1 the interesting result is [4]

r ( i ) = 5 + 5 2 1 + 5 2 = 2 + φ 1 φ 1 π φ 5 + ( π φ 5 10 ) 2 (6)

which can be approximated by terms containing φ 5 . About a nested infinitely continued fraction calculating φ ˜ = 0.619071 see [4]. Due to the property of the golden number its fifth power can be traced back to the simple relation φ 5 = 5 φ 3 (see also Appendix).

However, if one deals especially with the behavior of matter suffering from phase transitions, then the ideas of Baxter caused the deepest impact on physics [11] [12] [13] and will be reviewed in some detail in the following chapter. Then Hardy’s quantum probability function of two particles was investigated in more detail and connected with results of mathematical statistics and statistical mechanics [14] [15]. Surprisingly, Hardy’s function is found in the matter energy density relation resulting from Suleiman’s information relativity theory (IRT), too [9] [16]. The maximum of this energy density relation at a recession velocity of φ yielded exactly φ5 and was attributed to criticality of a transition at cosmic scale [9]. Another phase transition of economic relevance, the superconducting transition, also indicates a φ5 relation [17] [18] [19]. After this a further topic is included, namely Tammes problem of the largest diameter of N non-overlapping circles on the surface of a sphere with its connection to viral morphology and crystallography [20]. Finally, Fibonacci anyons proposed for topological quantum computation (TQC) were briefly described [21] [22] in comparison to the recently formulated reverse Fibonacci approach using the Janičko number sequence [23]. All these ideas are pioneering for applications including, besides mathematics and geometry, cosmology, crystallography, electrostatics, music, botany, viral morphology as well as quantum computation in connection with human mind and consciousness, and point to the omnipresence of the golden mean as the numerical dominator of our world.

It remains to be noted that it was a long way before Baxters ideas were fully accepted by down-to-earth physicists, and it is hoped that the new ideas formulated by Suleiman respectively El Naschie and scientific fellows will be soon taken over by the scientific community.

2. Phase Transition in the Hard Hexagon Lattice Gas Model

Simplistic models of statistical mechanics covering a large system of particles are important in linking microscopic and macroscopic average properties and behavior of matter undergoing a transition from disordered states to ordered ones by density increase or temperature reduction. The Ising model of a ferromagnet was the first exactly solved two-dimensional model defined on a square lattice [24].

In the following we want to pay tribute to Baxter’s important research in the field of statistical mechanics [11] [12] [13], dealing with the geometry and phase transition of a hard two-dimensional hexagonal lattice gas model, where hard means non-overlapping hexagons. Figure 1 illustrates a possible distribution of such non-overlapping hexagons in a triangular lattice cutout. A synopsis of results besides selected references for models of rigid gas molecules on a regular lattice is given by Finch [25]. Baxter was able to give an almost exact solution of this hard-hexagon model. For such a model the grand-canonical partition function Z is given by

Z ( z ) = n = 0 N 3 z n g ( n , N ) (7)

where the positive real variable z denotes activity or fugacity, and g ( n , N ) defines the number of allowed ways of placing n particles (molecules) on the lattice. One can write the partition function per site of an infinite lattice as

κ = lim n Z 1 / N (8)

Low activity is equivalent to homogeneity, whereas high activity complies with heterogeneity. In this way, a phase transition is expected. Baxter discovered two distinct regions of positive fugacity with criticality at the point z = zc

0 z < z c and z c < z < (9)

Figure 1. Statistical sites of hard (non-overlapping) hexagons in a triangular lattice cutout.

Quite surprisingly, the result for the critical zc is numerically equal to the reciprocal fifth power of the golden ratio

z c = φ 5 = 11 + φ 5 = 11.0901699 (10)

A second (non-physical) singularity was found for z n p = φ 5 = 0.090169943

The two singularities are the roots of the simple quadratic equation

z c 2 a z c 1 = 0 , (11)

where a = 11. For a = −1 one simply gets the golden mean φ = 5 1 2 as one of

the solutions. So we are once again faced with the importance of the golden mean in statistical mechanics.

The behavior of the mean density ρ(z) and the order parameter R(z) at the

critical point has been calculated, too. By using the normalized value z r = z z c one confirms

ρ ( z ) ~ ρ c 5 3 2 ( 1 z r ) 2 3 = ρ c 1 φ 5 + φ 5 ( 1 z r ) 2 3 , z z c (low density branch) (12)


ρ c = 5 5 10 = φ 5 = 0.276393202 (13)

R ( z ) = 3 5 [ 1 5 5 ( z r 1 ) ] 1 9 = 3 5 [ 1 φ 5 + φ 5 ( z r 1 ) ] 1 9 , z z c + (high density branch) (14)

It is worth to note that the critical exponents 2/3 and 1/9 are conjectured to be universal. Further one should notice that the pre-factors given in the original publication are simply the reciprocal of the sum of the fifth power of the golden mean and its inverse

5 3 2 = 1 5 5 = 1 φ 5 + φ 5 = 11 + 2 φ 5 . (15)

For a more precise study of the hard hexagon partition function in the complex fugacity plane see the valuable contribution published by Assis et al. [26].

Following Tracy et al. [27], the fugacity z(τ) can be derived from Klein’s icosahedron function ς(τ) [27] [28] as a branched cover (analytic bijection) of the Riemann sphere giving

z ( τ ) = { ς 5 ( τ ) disordered regime ς 5 ( τ ) ordered regime (16)

which is compatible with Baxters results for the hard hexagon approach [11] [12] [13]. The close connection between an icosahedron as a regular solid and the theory of groups should be noted here. Re-expressing the approach of Assis et al. [26] by likewise using the complex variable z = ς5, the golden ratio φ and the fifth root of unity ω 1 (see Appendix).

ω 1 = φ 2 + i 2 1 + φ 2 , | ω 1 | = 1 , (17)

and applying the h5 order-5 transformation

ς h 5 ( ς ) = ω 1 φ ς ω 1 φ + ς (18)

then due to this transformation the Hauptmodul H(ς) is invariant by the involution ς ς 1 and as well by the order-5 transformation ς ω 1 ς

H ( ς ) = H ( h 5 ( ς ) ) = H ( ς 1 ) = H ( ω 1 ς ) (19)

3. Hard Square Lattice Model

It was also Baxter [12] [13] who introduced the two-dimensional hard-square lattice model. Figure 2 depicts a statistical site occupation cutout based on an illustration given by Finch [25].

Baxter showed that this model is integrable resulting in two critical points at fugacities zc

z c = ± φ 5 (20)

Again the fifth power of the golden mean is involved in the result. Baxter’s famous work was seemingly overlooked and therefore not quoted in the Science article about quantum criticality, which was experimentally verified in the golden ratio bearing spin dynamics next a phase critical point of the quasi-one-dimensional Ising ferromagnet CoNb2O8 published by Coldea et al. [29].

Beyond it should be mentioned that the quantum Hamiltonian of a one-dimensional hard-boson model was found to map the transfer matrix of Baxter’s hard-square model with a quantum critical point at [30]

V w = φ 5 (21)

where V represents the nearest-neighbor (two sites apart) interaction, U the chemical potential of the bosons, and w 2 = V 2 + U V . The Hilbert space of the chosen one-dimensional quantum approach is shown to be identical of the space of states along a line of the hard-square model [30].

Figure 2. Hard squares sites on the square lattice [25].

4. Results of ε-Infinity Theory and Beyond

El Naschie’s E-infinity (ε) theory [5] originates from a fractal Cantorian set theory [31] as a number-theoretical route of physics for explaining the dualism between particles and waves that can serve solving cosmological mysteries such as dark matter and dark energy [5]. The approach delivers effective quantum gravity formulas for the cosmological mass (energy) constituents of baryonic matter eM, dark matter eDM, entire dark constituents eED, and pure dark energy ePD as follows

e M = φ 5 2 = 0.04508497 (22)

e E D = 1 e M = 5 2 φ 2 = 0.9549150 (23)

e D M = 3 2 φ 4 = 0.218847 (24)

e P D = 2 φ 1 2 = 0.736068 (25)

e M + e D M + e P D = 1 (26)

The readers are faced in the resulting formulas with the golden mean, especially with its fifth power in the beautiful relation

e M + e E D = φ 5 2 + 5 2 φ 2 = 1 (27)

where matter is completed by its dark wavy surrounding.

In a previously published contribution the present author recast both the matter amount eM and dark matter one eDM into a suitable form by combining φ 5 and the inverse number to suggest a reciprocity relation between both matter entities according to [17] [19]

e M = 1 10 5 φ 5 = 0.04508 (28)

e D M = 1 10 ( 5 φ 5 ) 1 = 0.22180 (29)

Then a persuasive equation for the pure dark energy ePD can be written down [17] [19]

e P D = 1 1 10 ( 5 φ 5 + ( 5 φ 5 ) 1 ) = 0.7331 ( 73.31 % ) (30)

Such quantum entanglement based coincidence means that the constituents of the cosmos should not be considered independent of each other, which was impressively confirmed later by the IRT theory of Suleiman [9].

5. Hardy’s Quantum Probability and Suleiman’s Matter Energy Density Relation

Hardy’s maximum quantum probability of two quantum particles [14] [15] exactly equals the fifth power of φ (Figure 1). This asymmetric probability distribution function with pτ as entanglement variable, running from not entangled states to completely entangled ones, is given by

P = p τ 2 1 p τ 1 + p τ (31)

The maximum of P yielded

P max = 1 φ 1 + φ φ 2 = φ 5 = 0.090169943 (32)

The Hardy function, displayed in Figure 3, turns out to be a central topic of the scale-free Information Relativity theory (IRT) of Suleiman [9] by mapping the transformation of his relative matter energy density. Suleiman characterized the behavior at the critical recession velocity β c r = φ as phase criticality at cosmic scale [9].

Now we want to study this function in more detail, replacing the variable pτ by x and write the function P as

f ( x ) = x 2 1 x 1 + x (33)

Integration of f ( x ) leads to

1 2 0 1 f ( x ) d x = ln ( 2 ) 2 3 = 0.026480514 1 12 π (34)

The first derivative of f ( x ) resulted in

f ( x ) = d f ( x ) d x = 2 x ( 1 x x 2 ) ( x + 1 ) 2 (35)

By comparing the nominator term ( x 2 + x 1 ) with Equation (11) the reference to the golden ratio becomes evident. The remaining multiplicative term caused a discontinuity at x = 1 .

Higher derivatives of f ( x ) lead to the following results

f ( x ) = 2 ( 1 x 3 3 x 2 3 x ) ( x + 1 ) 3 (36)

f ( x ) = 12 ( x + 1 ) 4 (37)

Special values for f ( x ) , f ( x ) and higher derivatives were summarized in the following Table 1.

If one connects the original function with energy, then its third derivative should be somewhat like a pressure. More surprisingly, the third derivative comes along as a very simple function with an origin value of −12 at x = 0 and a value of −12φ4 at x = φ, underscoring the magic importance of number 12, found for instance as the number of vertices of the Platonic solid icosahedron.

Figure 3. Hardy’s quantum probability P = ƒ( x ) for two particles [14] respectively Suleiman’s matter energy density [9] versus the entanglement variable x (or recession velocity by Suleiman).

Table 1. Special values for the functions f ( x ) and derivatives.

6. Superconductivity

However, nature presents much more relationships to keep in mind, where the golden mean or its fifth power is involved, and superconductivity is no exception. So we must reassess the theory considering the dark matter surrounding the moving electrons, which dive into the dark after marriage, or in other words, become superconducting under special conditions. Recently, the present author suggested linking the optimum hole doping σ 0 of high-Tc superconductors with the golden mean again in the form of its fifth power [17] [18] [19]

σ 0 8 π φ 5 = 0.2293 (38)

It was suggested recently that the same is true for conventional superconductors [18]. Obviously, this optimum is confirmed again near a quantum critical point in the superconductor phase diagram. In addition, the relation of the Fermi speed to the Klitzing speed comes out as [17]

v F v K 2 π φ 5 = 0.0571 (39)

Both relations document the fractal nature of the electronic response in superconductors. Also Prester had reported before about evidence of a fractal dissipative regime in high-Tc superconductors [32].

In Figure 4 earlier results [17] that illustrate the dependence between the critical temperature and the mean cationic charge q c for unconventional superconductors were completed with data points for the lanthanum and yttrium super-hydrides respectively H3S, showing the branch of n-type superconductors besides the p-type branch.

The prospect of developing a resilient theory of superconductivity is promised. Quantum entanglement of two moving electrons is influenced by local interaction of their interwoven dark matter surroundings, quoting the cogwheel picture of Suleiman [9]. In a recent contribution the present author designed the picture of two stretched electrons that locally interact to become superconducting. Such particle stretching into strands may in the end lead to a double-helically wounded wavy entity, which escapes in the dark.

Figure 4. Mean cationic charge versus critical temperature Tco(K) of superconducting compounds. p: p-type branch, n: n-type branches having σ < σo and σ > σo (σ equivalent to q c ). For more details see [17].

7. Tammes Problem of Decorating a Unit Sphere with Hard Circles

In 1930 the Dutch biologist Tammes [20] posed the question about the number and the arrangement of exit points in pollen grains. Reformulating the problem leads to the question what is the largest diameter of N non-overlapping equal circles placed on the surface of a unit sphere, referred to as hard-spheres problem [33] [34] [35]. A sphere can’t be tiled by using only hexagons. However, a perturbation of some hexagons in the hexagonal lattice into pentagons enables such a fit on the sphere, exemplified by the soccer ball structure [35]. In the following the energy of interaction is considered resulting from the arrangements of N points on the surface of a unit sphere S 2 : = { x R 3 : | x | = 1 } . The interaction is given by the power low potential V = | x y | α , 2 < α < 2 of the Euclidian

distance between two points x , y S 2 . For α = 0 the potential V = ln ( 1 | x y | )

is used. Given a fixed configuration of N points ω N = ( x 1 , x 2 , , x N ) , the α-energy sums up as [33] [36]

E ( α , ω N ) : = { 1 i < j N ln 1 | x i x j | , if α = 0 1 i < j N | x i x j | α , if α 0 (40)

Of importance is the determination of extremal energies for N points on the sphere [34] [36]

ε ( α , N ) : = { inf ω N S 2 E ( α , ω N ) , if α 0 sup ω N S 2 E ( α , ω N ) , if α > 0 (41)

Our concern should now be the special case ε ( 0 , N ) for which a conjecture already exists [34] [36]

ε ( 0 , N ) 1 4 ln ( 4 e ) N 2 1 4 N ln N + O ( N ) (42)

lim N E ( 0 , N ) ( 1 4 ln ( 4 e ) N 2 1 4 N ln N ) N = μ (43)

with an estimate of μ = 0.026422 [36]. The bound μ may be compared with the result given for the integral over the matter energy density (or even the Hardy function) according to Equation (35). We assume that this last value should be the correct bound

μ = 2 3 ln ( 2 ) = 0.0264805 (44)

Quoting Saff et al. [37] respectively Rakhmanov et al. [36], numerical calculations indicated that logarithmic equilibrium points (α = 0) with a generalized spiral point set ω ^ N have the tendency to distribute themselves over a nearly regular spherical hexagonal net giving again a connection with μ (our interpretation)

E ( 0 , ω ^ N ) ε ( 0 , N ) 114 ln N 3 | μ | ln N (45)

The interesting connection between minimal discrete energy on the sphere and matter energy density of the cosmos according to the information relativity theory [9] may suggest far-reaching consequences and should therefore be investigated more detailed in future.

8. Viral Morphology and Beyond

Nowadays we are faced with the worldwide pandemic spread of a very mobile corona virus. Therefore, the viral self-assembly and reproduction as a thermodynamic process comes ones again into focus. Viral self-assembly is intimately related to the before mentioned Tammes problem [20] of covering the surface of a sphere by optimally packed N circular disks with the alteration that in case of viruses the coverage may be conformational-optimized by disks-like protein capsomers (subunits) having two different radii. Some time ago, Bruinsma et al. [38] showed that such a two-radius model as a form of structural free energy minimization prefers the icosahedral symmetry as well as the transition of sphere-like viruses to rod-like ones such as found for bucky-tube aggregates.

Models for the viral spread are quite complex, but may be likewise considered as noise dynamical phase transition with a critical point [39]. Here the golden number may come into play, which had dominated the development of first self-assembling precursors of life as well as later life itself and the ability to process and compress information with extreme speed.

Heretically put, the deadly virus acts as a corrective to human stupidity by clipping human dominance and restlessness for the good of the environment. Optimistically, the mankind, if it survives at all, could be able in the distant future to control respectively destroy any small pathogenic bodies via teleported information, influencing the dark matter halo of the particles.

In anticipation of the explanations in the next chapter, we may ask whether quantum computation with (stable) viral structures could be possible in future or by use of artificial microtubules build from helically rolled up hexagon lattice entities similar to tubulin proteins that then may operate under ambient conditions. The author recommends reading the seminal contributions of Penrose and Hameroff [40] [41] and the recently published remarkable contribution of Olson et al. [2] explaining in more detail parts of this topic. Tubulin protein molecules form linear chains which are able to self-assemble into 2D sheets that can roll up to microtubules (Figure 5) [42].

The concept of a special arrangement of microtubules into a “heavy” hexagonal lattice with numbers that follow the Fibonacci number sequence and with vertices that minimize frequency collisions to protect the system from de-coherence, has been effectively refined again and again in nature by evolution. Not surprisingly, it resembles somewhat the qubit architecture of IBM’s last quantum computer (QC) [43]. As in buckytubes, ballistic electron transport through the microtubules should be an essential property. Consequently, particles moving with

Figure 5. Sketch of a microtubule of helically rolled up tubulin protein subunit chains.

ballistic speed are entangled through their pronounced dark matter halo with the environment and can by that pass on and store information. This is an important result of the IRT of Suleiman [9].

9. Quantum Computation with the Reverse Fibonacci Sequence and Fibonacci Anyons

The golden mean φ as the most irrational number found in nature with the simplest continued fraction representation at all is approximated by the division of two consecutive numbers of the Fibonacci number sequence [44] [45]. However, Janičko [23] recently introduced a new number sequence that he named reverse Fibonacci or Janičko sequence. The reverse sequence may be represented by the recursive formula beginning with the first two numbers as 0 and 1

J n + 2 = 8 ( J n + 1 J n ) (46)

Comparable to the quadratic Equation (7), where solutions lead to the golden mean, we can formulate a quadratic equation for the reverse sequence, too. It yields

J 2 8 J + 8 = 0 (47)

giving fundamental number j1,2 as its two solutions [23]

j 1 = 4 + 8 = 2 ( 2 + 2 ) = 6.82842712475 (48)

j 2 = 4 8 = 2 ( 2 2 ) = 1.171572875 (49)

As in the case of the golden mean also the reverse number j 1 1 seems to be important

j 1 1 = 1 j 2 1 = j 2 8 = 0.1464466094067 (50)

Again, both numbers can be approximated by the ratio of consecutive Janičko numbers [23].

This number sequence may be important in calculations with quantum bits. The particle-wave duality is essential in the quantum information theory, where the unit of information is given by the quantum bit (qubit) coined by Schumacher [46], which exhibits the aspects of particle localization (counting) and wave interference to represent a signal with high fidelity [46]. Such two-state quantum system can be represented by the superposition principle: S U ( 2 ) H 1 Q = α | 0 + β | 1 , where α and β are complex numbers, where α 2 + β 2 = 1 [46]. A single qubit state can also be written as

| ψ = p r | 0 + e i ϕ 1 p r | 1 (51)

where 0 p r 1 is the probability of the bit being in the 0 state, and 0 ϕ 2 π is the phase of the quantum. All possible values of pr and the phase ϕ are reachable by application of quantum gates, realized by reliable and repeatable Josephson junctions. For instance, operating the quantum logical T-gate as quantum circuit with a phase of π/4 that is mapped by the unitary matrix

T = ( 1 0 0 e i π / 4 ) (52)

one gets for p r ( 0 ) = j 2 1 = 0.8535533 respectively for p r ( 1 ) = j 1 1 = 0.1464466 . In this way the connection of a single qubit with Janičko numbers is demonstrated as before suggested by Janičko itself [23]. The number j 1 1 is near the fourth power of φ : φ 4 = 0.14589803 and j 2 2 φ = 1.17557050 .

In contrast, topological quantum computation using Fibonacci anyons as qubits [47] [48] may be the ultimate approach because it naturally simulates processes which determine the speed and storage capacity of human mind. Fibonacci anyons have been proposed as two-dimensional quasi-particles exhibiting fractional quantum Hall states at 5/2 respectively 12/5 filling fractions. As an example, a qubit can be build up from four Fibonacci anyons, where the fourth “inert” anyon may serve for error correction. The world lines of such anyons pass around one another to form braids. Braids can be thought of as quantum circuits forming the logic gates, which building up the computer [49]. Their non-Abelian braiding statistics with only two particle types are shown to be suitable for quantum computation. The two particle types or topological charges are the trivial type 1 and the non-trivial type τ, which are self-dual and being their own antiparticles (1 = 1*, τ = τ*). The combination or fusion of anyons governed by fusion rules may be thought of as a measurement, graphically displayed by a fusion tree [50]. The probability of a fusion is denoted as quantum dimension di. The quantum dimension of τ is φ 1 , the inverse of the golden mean. When two τ anyons are combined, the probability pr to yield 1 is p r 0 = φ 2 , and to yield τ is p r 1 = φ [50].

For a multi-anyon system the number of fusion paths resulted exactly in the Fibonacci number, and the move accompanied by the consecutive fuse of the anyons from left to right can be represented by a F matrix. For more than three anyons consistency equations for the fusion trees in form of pentagon relations deliver the entries of the F-matrices [50]. Up to arbitrary phases the 2 × 2 unitary F matrix reads

F τ τ τ τ = ( F 11 F 1 τ F τ 1 F τ τ ) = ( φ φ φ φ ) (53)

where again φ is the golden number (Figure 6).

Figure 6. Removing crossings by F and R moves after Slingerland [51].

The braiding operator B can be derived by matrix multiplication of F-matrices with R-matrices. When the R-matrix is given by diagonal spin factors

R = d i a g ( e 4 π i 5 , e 3 π i 5 , e 3 π i 5 , e 4 π i 5 , e 3 π i 5 ) (54)

and the extended 5 × 5 F-matrix by

F = ( 1 1 1 φ φ φ φ ) (55)

then one obtains finally the B matrix as [50]

B = F 1 R F = ( e 4 π i / 5 e 3 π i / 5 e 3 π i / 5 φ e 4 π i / 5 φ e 2 π i / 5 φ e 2 π i / 5 φ ) (56)

“Pentagon” moves respectively “hexagon” ones are third order polynomial equations in many variables with much more equations than variables, but not always existing solutions.

Freedman, Larsen and Wang [49] confirmed the universal applicability of the following approximated quantity pr (probability) for quantum computation

p r = | V ( L , t ) ( t 1 2 + t 1 2 ) g 1 | 2 (57)

where V ( L , t ) is the Jones polynomial of a link L [52], t = exp ( 2 π i r ) is a principal root of unity with r = 5 or r 7 , and g is the bridge number1. Lattice roots

of unity, expressed in the parlance of a crystallographer, with r { 1 , 2 , 3 , 4 , 6 } are excluded. Field and Simula [53] reported recently about the connection between

Fibonacci anyons and the Jones polynomial at the fifth root of unity t = exp ( 2 π i 5 ) .

The magnitude of the Jones polynomial [47] [54] of the Hopf link2 is related to the probability p r ( | 0 5 ) simply by [53]

V H o p f ( e 2 π i 5 ) = p r ( | 0 5 ) φ 1 = φ (58)

with only two elementary physical braiding operations used before measuring the annihilation probability of the anyons upon fusion [53].

This short survey may serve to demonstrate the importance of the golden number for the development of quantum computation that should stimulate research in a direction along which way nature has always evolutionary and subtly worked. In this sense the seminal ideas of Penrose and Hameroff [40] [41] should be further elaborated (see Section 7). In the next chapter new golden-mean-based quantum computer architecture was proposed copying nature’s most effective principles.

10. Proposal of Quantum Computer Architecture Based on Fibonacci Net Microtubules

Besides the compounds currently applied such as GaAs hetero-structures or graphene, which show the fractional quantum Hall effect (FQHE), the investigation of other inorganic (semiconducting) compounds is recommended, for instance ferroelastic compounds or compounds that show Fibonacci-like unit-cell decoration [55] [56]. Interfacial twin walls act as sink for charge carriers or defects generating locally superconductivity or ferroelectricity as striking properties [57]. Fibonacci unit-cell tailoring was displayed in Figure 7. Geometric relations (Fibonacci relations) between these cells can be followed by studying [55] [56] respectively in the Appendix. The small black hexagonal cell (sub-cell) is rotated

by an angle of α = atan ( 5 3 3 ) 30 = 13.898 to generate the blue outlined

unit-cell with lattice parameter a that has 13 times the volume of the small cell. One may repeat this procedure to obtain a super-cell and so on. Furthermore, it should be possible to roll up this 2D net into a microtubule (Figure 8) and plug tubules with different radii or orientation into each other to produce a multi-shell structure. The starting small cell is then helically wound at the twist angle α. It needs 13-times the small cell lattice parameter asub to reach again identity with a large cell lattice point. In this way the diameter of the smallest geometrically possible microtubule is calculated to be

d 1 13 a s u b cos a s u b cos ( α ) / π = 3.5 a / π . Interesting is d 6 = 21 a π showing a further Fibonacci number. By a full turn one gains a height in filament direction of h 1 = 39 2 a s u b π a s u b . If you look through the microtubule perpendicular to the filament direction, you will see the mirror image on the back somewhat offset.

Figure 7. Fibonacci arrangement of a hexagonal net and its mirror image [55] [56]. The light-blue outlined unit-cell contains 13 sub-cells, offset by an angle of α = 13.9˚.

Figure 8. Helically twisted microtubule projected down the filament direction with 13 light-blue atoms or atom groups on sub-lattice positions.

In practice it may be a sophisticated task to produce such a tubule, but nowadays it should be possible. The future will show whether such a composite is applicable for QC, after the right compound with optimal physical properties is found. One could start with doped graphene rolled up in the explained Fibonacci twisting. The Kitaev honeycomb lattice approach may be applied in an extended form to this special lattice variant [58].

Interestingly, the helically twisted tubulin microtubules show strongest reinforcement of the ordered pattern when the protofilament number n equals 13 [59] [60]. This supports the proposal for a twisted Fibonnaci net as displayed in Figure 7.

A very stimulating question is what happens if the microtubule is compressed or dilated due to vibration in filament direction finally going through an incommensurability state represented by 13 + φ 5 = ( 2 + φ 1 ) 2 = 13.0901699 instead of 13 twist steps (see Appendix). One may tentatively characterize such a behavior also by a phase transition. A highly effective quantum computer could possibly be realized based on such a principle. Besides the explained longitudinal adaption of the twisted net, its sense of rotation, optimal diameter, and the freedom given by the decoration of the net itself, there are many more degrees of freedom, if one constructs a hexagonal lattice of such 13-step microtubules with non-hexagonal internal symmetry. Neighboring filament can be rotated one again the other violating the hexagonal symmetry. In addition, the border of the assembled filaments may be important. This is what nature has already optimized to equip its creatures with perception, memory and intelligence.

11. Conclusions

This contribution reviewed and extended results of different scientific disciplines that involve the golden mean φ as very important number of nature. It covered especially phase transitions governed by the fifth power of φ, found in Baxter’s hard lattice approaches of statistical mechanics as well as in energy density relations of the cosmos, comparing results of the ε-infinity theory of El Naschie with the conclusive findings of Suleiman’s scale-free information relativity theory (IRT) besides Hardy’s quantum entanglement probability. Also superconductivity is described as an inherently fractal property connected with the interaction of two moving electrons entangled and finally paired by way of their dark wavy surroundings. Since Klein’s research on the icosahedron many others solved and combined fundamental scientific questions beyond icosahedral structures. From Tammes problem of decorating a unit sphere with hard circles one may be guided to viral morphology as well as phase transition, self-assembly and reproduction of viruses. Are there connections between the energy density of the cosmos with Tammes problem?

The golden mean is involved in the development of modern topological quantum computation, too. The performance of topological quantum computation with Fibonacci anyons can give a vague idea of how the human brain actually processes and stores information. It may be highly interesting to develop artificial self-assembling structures based for instance on the tubulin protein or viruses that can form helically wounded microtubules with the ability to store and to process information. A proposal was made for a twisted microtubule tailors from a rolled up Fibonacci net. Quantum computation is easier understandable, if one accepts the suggestions of Suleiman that quantum entanglement is a local realistic phenomenon “where each moving particle will be permeated by the dark matter halo of the other one making the two physically entangled” [15]. Essential suggestions of Irvin and scientific fellows about the golden ratio as a fundamental constant of nature could be ennobled by considering the IRT theory [8] [61]. The author wants to encourage “down-to-earth” physicists to follow and complete these suggestions and the perspective views of the present author, including the proposed change of some important physical constants [19] [62].


The author appreciates valuable suggestions given by a competent and cooperative reviewer.


Roots of unity: By solving the quartic polynomial ϕ 5 ( x ) = n = 0 4 x n one yields the four primitive roots of unity ω i (i = 1 to 4) here notably expressed by terms of the golden mean

ω 1 , 2 = φ 2 ± i 2 1 + φ 2 = φ 2 ± i 2 2 + φ 1 (59)

ω 3 , 4 = φ 1 2 ± i 2 1 + φ 2 = φ 1 2 ± i 2 2 φ (60)

About the Fibonacci number 13: This number is very special because its square root, coming along as the unit-cell—sub-cell parameter ratio of the Fibonacci net approach, can be approximated by 2 + φ 1

13 = 3.60555 2 + φ 1 = 3.6180339887 (61)

Furthermore, the exact adjustment connects the number 13 with the fifth power of φ, which number was the main concern of this contribution

2 + φ 1 = 5 φ 1 = 13 + φ 5 = 13.0901699 = 3.6180339887 (62)

When adding to the Hardy function [14] respectively Suleiman’s matter energy density relation [9] an amount of 3 as a solely mathematical procedure, for the present without physical background, giving then the maximum of 5 φ instead of φ 5 at x = φ , the function is changed to

f ( x ) = 3 + ( 1 x ) x 2 ( 1 + x ) (63)

An exceptional result is obtained for the x value at zero of this function

x = 5 3 ( 10 2 3 + 10 1 3 + 10 0 3 ) = 5 × 2.59867450788 = 12.993372 13 (64)

Crystallography of the Fibonacci 2D-Lattice: The reciprocal sub-cell is put up according to [55] by the strongest hexagonal basis reflections 31 4 ¯ 0 and 52 7 ¯ 0 . Their Miller-Bravais index relations represent Fibonacci numbers respectively products of them

( 31 4 ¯ 0 ) : h 2 + k 2 + h k = 13 (65a)

( 52 7 ¯ 0 ) : h 2 + k 2 + h k = 39 = 3 × 13 (65b)

Also remarkable, the sine of the twist angle α = 13.898 can be approximated well by the ratio of π and number 13 + φ 5

sin ( 13.898 ) = 0.240192 π 13 + φ 5 = 0.239996 = sin ( 13.886 ) (66)

The Fibonacci net can be generated by application of symmetry operations of the two-dimensional symmetry group p6. One needs three positions to get the 13 lattice points. Two general six-fold positions with coordinates in fractions of again

13, besides the origin, yield the lattice points:;; x 3 = 5 13 , y 3 = 7 13 (Figure 9).

Figure 9. Illustration of the three differently colored atomic positions of the 2D Fibonacci net.

Cabinet of curiosities: With this title Fang et al. [63] published a result concerning the offset angle of chained aggregates of tetrahedra with coincident faces given by the relation

β = a r c c o s ( 3 φ 1 1 4 ) = 15.5224 (67)

However, one can approximate this angle very well by a relation that uses the fifth power of φ

β 3 a r c s i n ( φ 5 ) = 15.52015 (68)

Geometric interpretation of φ5: The angle of β / 3 = 5.173 is also near the difference of the half inner angles between pentagon and hexagon being 36 30 = 6 . The difference in the triangle areas between both regular polygons yielded a value near φ 5 = 0.090169943

A 6 6 A 5 5 = a 2 4 ( 3 1 + 1 φ + 0.5 ) = 0.088917221 a 2 (69)

where A 6 and A 5 are the full areas of the hexagon respectively pentagon with edge length a . When the pentagon is only slightly rotated about the edge against the hexagon plane, and then the projected area is subtracted, one can well adapt the exact value of φ 5 .


1The bridge number g = 2n is the number of -turns at each end of a capped braid (plat diagram).

2One speaks of a Hopf link, representing a (2, 2)-torus link, when two circles are linked together with a crossing number of two.

Cite this paper: Otto, H. (2020) Phase Transitions Governed by the Fifth Power of the Golden Mean and Beyond. World Journal of Condensed Matter Physics, 10, 135-158. doi: 10.4236/wjcmp.2020.103009.

[1]   Olson, S. (2006) The Golden Section: Nature’s Greatest Secret. Bloomsbury, London, 64 p.

[2]   Olsen, S., Marek-Crnjak, L., He, J.H. and El Naschie, M.S. (2020) A Grand Unification of the Sciences, Art & Consciousness: Rediscovering the Pythagorean Plato’s Golden Mean Number System. Journal of Progressive Research in Mathematics, 16, 2880-2931.

[3]   Sherbon, M.A. (2014) Fundamental Nature of the Fine-Structure Constant. International Journal of Physical Research, 2, 1-9.

[4]   Otto, H.H. (2017) Continued Fraction Representations of Universal Numbers and Approximations., 1-4.

[5]   El Naschie, M.S. (2013) Quantum Entanglement: Where Dark Energy and Negative Gravity plus Accelerated Expansion of the Universe Comes from. Journal of Quantum Information Science, 3, 57-77.

[6]   Ho, M.-W. (2014) Golden Geometry and E-Infinity Fractal Spacetime. Story of Phi Part 5. Researchgate.

[7]   Marek-Crnjac, L. (2013) Cantorian Space-Time Theory. Lambert Academic Publishing, Saarbrücken, 1-50.

[8]   Irwin, K. (2019) Towards the Unification of Physics and Number Theory. Reports in Advances of Physical Sciences 3, Article ID: 1950003.

[9]   Suleiman, R. (2019) Relativizing Newton. Nova Scientific Publisher, New York, 1-207.

[10]   Andrews, G.E. (1981) Ramanujan’s “Lost” Notebook. III. The Rogers-Ramanujan Continued Fraction. Advances in Mathematics, 41, 186-208.

[11]   Baxter, R.J. (1980) Hard Hexagons: Exact Solution. Journal of Physics A: Mathematical and General, 13, L61-L70.

[12]   Baxter, R.J., Enting, I.G. and Tsang, S.K. (1980) Hard-Square Lattice Gas. Journal of Statistical Physics, 22, 465-489.

[13]   Baxter, R.J. (1982) Exactly Solved Models in Statistical Mechanics. Academic Press, London.

[14]   Hardy, L. (1993) Nonlocality for Two Particles without Inequalities for Almost All Entangled States. Physical Review Letters, 71, 1665-1668.

[15]   Mermin, N.D. (1994) Quantum Mysteries Refined. American Journal of Physics, 62, 880-887.

[16]   Otto, H.H. (2018) Reciprocity Relation between the Mass Constituents of the Universe and Hardy’s Quantum Entanglement Probability. World Journal of Condensed Matter Physics, 8, 30-35.

[17]   Otto, H.H. (2016) A Different Approach to High-Tc Superconductivity: Indication of Filamentary-Chaotic Conductance and Possible Routes to Superconductivity above Room Temperature. World Journal of Condensed Matter Physics, 6, 244-260.

[18]   Otto, H.H. (2019) Super-Hydrides of Lanthanum and Yttrium: On Optimal Conditions for Achieving near Room Temperature Superconductivity. World Journal of Condensed Matter Physics, 9, 22-36.

[19]   Otto, H.H. (2020) Reciprocity as an Ever-Present Dual Property of Everything. Journal of Modern Physics, 11, 98-121.

[20]   Tammes, P.M.L. (1930) On the Origin of Number and Arrangement of the Places of Exit on Pollen Grains. Dissertation, J.H. De Bussy, Groningen. (See also Recueil des Travaux Botanique Néerlandais, 27, 1-84)

[21]   Trebst, S., Troyer, M., Wang, Z., Ludwig, A.W.W. (2009) A Short Introduction to Fibonacci Anyon Models. arXiv: 0902.3275v1 [cond-mat.stat-mech], 1-24.

[22]   Field, B. and Simula, T. (2018) Introduction to Topological Quantum Computation with Non-Abelian Anyons. arXiv: 1802.06176v2 [quant-ph].

[23]   Janičko, O. (2018) New Fundamental Discovery of the Reverse Fibonacci Sequence., 1-6.

[24]   Onsager, L. (1944) Crystal Statistics. I. A Two-Dimensional Model with Order-Disorder Transition. Physical Review, 65, 117-149.

[25]   Finch, S.R. (2003) Mathematical Constants. Cambridge University Press, New York.

[26]   Assis, M., Jacobsen, J.L., Jensen, I., Maillard, J.M. and McCoy, B.M. (2013) The Hard Hexagon Partition Function for Complex Fugacity. Journal of Physics A: Mathematical and Theoretical, 46, Article ID: 445202.

[27]   Tracy, C.A., Grove, L. and Newman, M.F. (1987) Modular Properties of the Hard Hexagon Model. Journal of Statistical Physics, 48, 477-502.

[28]   Klein, F. (1884) Vorlesungen über das Ikosaeder und die Auflösung der Gleichungen vom fünften Grad. Verlag B. G. Teubner, Leipzig.

[29]   Coldea, R., Tennant, D.A., Wheeler, E.M., Wawrzynska, E., Prabhakaram, D., Telling, M., Habicht, K., Smeibidl, P. and Kiefer, K. (2010) Quantum Criticality in an Ising Chain: Experimental Evidence for Emergent E8 Symmetry. Science, 327, 177-180.

[30]   Fendley, P., Sengupta, K. and Sanchev, S. (2004) Competing Density-Wave Orders in a One-Dimensional Hard-Boson Model. Physical Review B, 69, Article ID: 075106.

[31]   Cantor, G. (1932) Gesammelte Abhandlungen. Mathematischen und Philosophischen Inhalts. Springer, Berlin.

[32]   Prester, M. (1999) Experimental Evidence of a Fractal Dissipative Regime in High-Tc Superconductors. Physical Review B, 60, 3100-3103.

[33]   Wagner, W. (1990) On Means of Distances on the Surface of a Sphere (Lower Bounds). Pacific Journal of Mathematics, 114, 389-398.

[34]   Zhou, Y.M. (1995) Arrangements of Points on the Sphere. Ph.D. Thesis, University of South Florida, Tampa, FL.

[35]   Aste, T. and Weaire, D. (2000) The Pursuit of Perfect Packing. Taylor and Francis Group, Abingdon-on-Thames, 108-110.

[36]   Rakhmanov, E.A., Saff, E.B. and Zhou, Y.M. (1994) Minimal Discrete Energy on the Sphere. Mathematical Research Letters, 1, 647-662.

[37]   Saff, E.B. and Kuijlaars, A.B.J. (1997) Distributing Many Points on a Sphere. The Mathematical Intelligencer, 19, 5-11.

[38]   Bruinsma, R.F., Gelhart, W.M., Reguera, D., Rudnick, J. and Zandi, R. (2003) Viral Self-Assembly as a Thermodynamic Process. Physical Review Letters, 90, Article ID: 248101.

[39]   Drake, J.M., Brett, T.S., Chen, S., Epureanu, B.I., Ferrari, M.J., Marty, E., Miller, P.B., O’Dea, E.B., O’Regan, S.M., Park, A.W. and Rohani, P. (2019) The Statistics of Epidemic Transitions. PLoS Computational Biology, 15, e1006917.

[40]   Penrose, R. and Hamedoff, S.R. (1995) Quantum Computing in Microtubules: Self-Collapse as a Possible Mechanism for Consciousness. Toward a Science of Consciousness. MIT Press, Cambridge, MA.

[41]   Penrose, R. and Hamedoff, S.R. (2011) Consciousness in the Universe: Neuroscience, Quantum Space-Time Geometry and Orch OR Theory. Journal of Cosmology, 14, 1-50.

[42]   Sahu, S., Ghosh, S., Fujita, D. and Bandyopadhyay, A. (2014) Live Visualization of Single Isolated Tubulin Protein Self-Assembly via Tunneling Current Effect of Electromagnetic Pumping during Spontaneous Growth of Microtubule. Scientific Reports, 4, Article No. 7303.

[43]   Chamberland, C., Zhu, G., Yoder, T.J., Hertzberg, J.B. and Cross, A.W. (2019) Topological and Subsystem Codes on Low-Degree Graphs with Flag Qubits. ArXiv: 1907.09528v2 [quant-ph].

[44]   Pisano, L. (1202) Fibonacci’s Liber Abaci (Book of Calculation). Biblioteca Nazionale Centrale di Firenze, Firenze.

[45]   Lucas, E. (1891) Theorie des nombres. Gauthier-Villars, Paris.

[46]   Schumacher, B. (1995) Quantum Coding. Physical Review A, 51, 2738-2747.

[47]   Aharonov, D., Jones, V. and Landau, Z. (2006) A Polynomial Quantum Algorithm for Approximating the Jones Polynomial. Proceedings of the 38th Annual ACM Symposium on Theory of Computing, Seattle WA, May 2006, 427-436.

[48]   Nayak, C., Simon, S.H., Stern, A., Freedman, M. and Das Sarma, S. (2008) Non-Abelian Anyons and Topological Quantum Computation. Reviews of Modern Physics, 80, 1083-1159.

[49]   Freedman, M.H., Larsen, M.J. and Wang, Z. (2002) The Two-Eigenvalue Problem and Density of Jones Representation of Braid Groups. Communications in Mathematical Physics, 228, 177-199.

[50]   Trebst, S., Troyer, M., Zhenghan, W. and Ludwig, A.W.W. (2009) A Short Introduction to Fibonacci Anyon Models. arXiv: 0902.3275v1 [cond-mat.stat-mech].

[51]   Slingerland, J. (2007) Anyon Models—Theory, Interferometry, Bose Condensation. Presentation, UC Riverside, Caltech, California, USA.

[52]   Jones, V.F.R. (1985) A Polynomial Invariant for Knots via Von Neumann Algebras. Bulletin (New Series) of the American Mathematical Society, 12, 103-111.

[53]   Field, B. and Simula, T. (2018) Introduction to Topological Quantum Computation with Non-Abelian Anyons. Quantum Science and Technology, 3, Article ID: 045004.

[54]   Kuperberg, G. (2015) How Hard It Is to Approximate the Jones Polynomial? Theory of Computing, 11, 183-219.

[55]   Otto, H.H. (2015) Pyroelectric Bi5−x(Bi2S3)39I12S: Fibonacci Superstructure, Synthesis Options and Solar Cell Potential. World Journal of Condensed Matter Physics, 5, 66-77.

[56]   Otto, H.H. (1965) Zur Kristallchemie synthetischer Blei-Wismut-Spießglanze. Thesis, TU, Berlin.

[57]   Salje, E. (2012) Ferroelastic Materials. Annual Review of Materials Research, 42, 265-283.

[58]   Kitaev, A. (2008) Anyons in an Exactly Solved Model and Beyond. arXiv: 0506438v3 [cond-mat.mes-hall], 1-113.

[59]   Wade, R.H., Chrétien, D. and Job, D. (1990) Characterization of Microtubule Protofilament Numbers. How Does the Surface Lattice Accommodate? Journal of Molecular Biology, 212, 775-786.

[60]   Ti, S.C., Alushin, G.M. and Kapoor, T.M. (2018) Human β-Tubulin Isotypes Can Regulate Microtubule Protofilament Number and Stability. Developmental Cell, 47, 175-190.

[61]   Irvin, K., Amaral, M.M., Aschheim, R. and Fang, F. (2017) Quantum Walk on a Spin Network and the Golden Ratio as a Fundamental Constant of Nature.

[62]   Otto, H.H. (2020) A Primer of Important Natural Numbers and Revisited Fundamental Physical Constants.

[63]   Fang, F., Irvin, K., Kovacs, J. and Sadler, G. (2013) Cabinet of Curiosities: The Interesting Geometry of the Angle β.