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 φ⁵ and was attributed to criticality of a transition at cosmic scale [9]. Another phase transition of economic relevance, the superconducting transition, also indicates a φ⁵ 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.

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 = z_c

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

Quite surprisingly, the result for the critical z_c 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)

where

ρ 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 = ς⁵, the golden ratio φ and the fifth root of unity ω 1 (see Appendix).

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

and applying the h₅ 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)

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 z_c

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 CoNb₂O₈ 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].

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 e_M, dark matter e_DM, entire dark constituents e_ED, and pure dark energy e_PD 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 e_M and dark matter one e_DM 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 e_PD 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].

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φ⁴ at x = φ, underscoring the magic importance of number 12, found for instance as the number of vertices of the Platonic solid icosahedron.

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-T_c 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-T_c 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 H₃S, 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.

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.

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

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].

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 j_1,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 d_i. 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).

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 number¹. 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 link² 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.

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 a_sub 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.

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.

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.

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

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

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).

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 φ⁵: 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 .